Nucleotide substitution at major histocompatibility complex class II loci: evidence for overdominant selection

Considerable MHC diversity suggests that the functional extinction of baiji is not related to population genetic collapse

To further extend our understanding of the mechanism causing the current nearly extinct status of the baiji (Lipotes vexillifer), one of the most critically endangered species in the world, genetic diversity at the major histocompatibility complex (MHC) class II DRB locus was investigated in the baiji. Nine highly divergent DRB alleles were identified in 17 samples, with an average of 28.4 (13.2%) nucleotide difference and 16.7 (23.5%) amino acid difference between alleles. The unexpectedly high levels of DRB allelic diversity in the baiji may partly be attributable to its evolutionary adaptations to the freshwater environment which is regarded to have a higher parasite diversity compared to the marine environment. In addition, balancing selection was found to be the main mechanisms in generating sequence diversity at baiji DRB gene. Considerable sequence variation at the adaptive MHC genes despite of significant loss of neutral genetic variation in baiji genome might suggest that intense selection has overpowered random genetic drift as the main evolutionary forces, which further suggested that the critically endangered or nearly extinct status of the baiji is not an outcome of genetic collapse.

Diversity of MHC class II DAB1 in the koala (Phascolarctos cinereus)

The host immune response is an important factor determining the outcome of the host–pathogen–environment interaction. At the gateway between the innate and adaptive immune systems are MHC molecules, which facilitate antigen presentation to T lymphocytes, and initiate the adaptive immune response. Despite their integral role in adaptive immunity, the genes encoding class II MHC molecules have not been examined directly in koalas. Furthermore, indirect historical evidence suggests that this species might lack functional diversity in class II MHC genes, with potential implications for disease susceptibility. We have examined diversity in the β chain genes of the koala class II MHC DA gene family and identified 23 alleles, including several atypical alleles. The levels of diversity observed are consistent with other marsupial and eutherian species, and do not support the paucity of variation suggested by the early literature. These findings are relevant to the conservation management of koalas and provide both a benchmark for maintaining population diversity and a platform for further conservation genetic research in this species.

Impact of historical founder effects and a recent bottleneck on MHC variability in Commander Arctic foxes (Vulpes lagopus)

Populations of Arctic foxes (Vulpes lagopus) have been isolated on two of the Commander Islands (Bering and Mednyi) from the circumpolar distributed mainland population since the Pleistocene. In 1970-1980, an epizootic outbreak of mange caused a severe population decline on Mednyi Island. Genes of the major histocompatibility complex (MHC) play a primary role in infectious disease resistance. The main objectives of our study were to compare contemporary variation of MHC class II in mainland and island Arctic foxes, and to document the effects of the isolation and the recent bottleneck on MHC polymorphism by analyzing samples from historical and contemporary Arctic foxes. In 184 individuals, we found 25 unique MHC class II DRB and DQB alleles, and identified evidence of balancing selection maintaining allelic lineages over time at both loci. Twenty different MHC alleles were observed in mainland foxes and eight in Bering Island foxes. The historical Mednyi population contained five alleles and all contemporary individuals were monomorphic at both DRB and DQB. Our data indicate that despite positive and diversifying selection leading to elevated rates of amino acid replacement in functionally important antigen-binding sites, below a certain population size, balancing selection may not be strong enough to maintain genetic diversity in functionally important genes. This may have important fitness consequences and might explain the high pathogen susceptibility in some island populations. This is the first study that compares MHC diversity before and after a bottleneck in a wild canid population using DNA from museum samples.

Population genetics of duplicated alternatively spliced exons of the Dscam gene in Daphnia and Drosophila

In insects and crustaceans, the Down syndrome cell adhesion molecule (Dscam) occurs in many different isoforms. These are produced by mutually exclusive alternative splicing of dozens of tandem duplicated exons coding for parts or whole immunoglobulin (Ig) domains of the Dscam protein. This diversity plays a role in the development of the nervous system and also in the immune system. Structural analysis of the protein suggested candidate epitopes where binding to pathogens could occur. These epitopes are coded by regions of the duplicated exons and are therefore diverse within individuals. Here we apply molecular population genetics and molecular evolution analyses using Daphnia magna and several Drosophila species to investigate the potential role of natural selection in the divergence between orthologs of these duplicated exons among species, as well as between paralogous exons within species. We found no evidence for a role of positive selection in the divergence of these paralogous exons. However, the power of this test was low, and the fact that no signs of gene conversion between paralogous exons were found suggests that paralog diversity may nonetheless be maintained by selection. The analysis of orthologous exons in Drosophila and in Daphnia revealed an excess of non-synonymous polymorphisms in the epitopes putatively involved in pathogen binding. This may be a sign of balancing selection. Indeed, in Dr. melanogaster the same derived non-synonymous alleles segregate in several populations around the world. Yet other hallmarks of balancing selection were not found. Hence, we cannot rule out that the excess of non-synonymous polymorphisms is caused by segregating slightly deleterious alleles, thus potentially indicating reduced selective constraints in the putative pathogen binding epitopes of Dscam.

Hijacking of the Host SCF Ubiquitin Ligase Machinery by Plant Pathogens

The SCF (SKP1-CUL1-F-box protein) ubiquitin ligase complex mediates polyubiquitination of proteins targeted for degradation, thereby controlling a plethora of biological processes in eukaryotic cells. Although this ubiquitination machinery is found and functional only in eukaryotes, many non-eukaryotic pathogens also encode F-box proteins, the critical subunits of the SCF complex. Increasing evidence indicates that such non-eukaryotic F-box proteins play an essential role in subverting or exploiting the host ubiquitin/proteasome system for efficient pathogen infection. A recent bioinformatic analysis has identified more than 70 F-box proteins in 22 different bacterial species, suggesting that use of pathogen-encoded F-box effectors in the host cell may be a widespread infection strategy. In this review, we focus on plant pathogen-encoded F-box effectors, such as VirF of Agrobacterium tumefaciens, GALAs of Ralstonia solanacearum, and P0 of Poleroviruses, and discuss the molecular mechanism by which plant pathogens use these factors to manipulate the host cell for their own benefit.

Molecular cloning, polymorphism and tissue distribution of the MHC class IIB gene in the Chinese goose (Anser cygnoides)

1. The goose major histocompatibility complex (MHC) class IIB cDNA (Ancy-MHCII) was cloned by homology cloning and rapid amplification of cDNA ends by polymerase chain reaction (RACE-PCR), and the genomic structure and tissue expression were investigated. 2. Three different 5'-RACE sequences (Ancy-MHC II5'-1, Ancy-MHC II5'-2, Ancy-MHC II5'-3), one 3'-RACE sequence (Ancy-MHC II-3') and two different full length Ancy-MHC IIB cDNA sequences (Ancy-CD01, Ancy-CD02), which came from different alleles at one locus or different loci, were determined. 3. The genomic organisation is composed of 6 exons and 5 introns, with a longer intron region than that of the chicken. The alleles encode 259 and 260 amino acids in the mature protein. 4. The number of non-synonymous substitutions (dN) in the peptide-binding region of exon 2 from 8 alleles was higher than that of the synonymous substitutions (dS). 5. Tissue-specific expression of Ancy-MHC II mRNA was detected in an adult goose using RT-PCR. These results showed that Ancy-MHC II mRNA was expressed in the lung, spleen, liver, intestine, heart, kidney, pancreas, brain, skin and muscle. This is consistent with the expression of MHC class IIB in various tissues from the chicken. 6. Sequences from goose, snipe and duck clustered together when compared with known MHC class IIB sequences from the other species, significantly differing from mammals and aquatic species, indicating a pattern consistent with accepted evolutionary pathways.

Diversity and evolutionary patterns of immune genes in free-ranging Namibian leopards (Panthera pardus pardus)

The genes of the major histocompatibility complex (MHC) are a key component of the mammalian immune system and have become important molecular markers for fitness-related genetic variation in wildlife populations. Currently, no information about the MHC sequence variation and constitution in African leopards exists. In this study, we isolated and characterized genetic variation at the adaptively most important region of MHC class I and MHC class II-DRB genes in 25 free-ranging African leopards from Namibia and investigated the mechanisms that generate and maintain MHC polymorphism in the species. Using single-stranded conformation polymorphism analysis and direct sequencing, we detected 6 MHC class I and 6 MHC class II-DRB sequences, which likely correspond to at least 3 MHC class I and 3 MHC class II-DRB loci. Amino acid sequence variation in both MHC classes was higher or similar in comparison to other reported felids. We found signatures of positive selection shaping the diversity of MHC class I and MHC class II-DRB loci during the evolutionary history of the species. A comparison of MHC class I and MHC class II-DRB sequences of the leopard to those of other felids revealed a trans-species mode of evolution. In addition, the evolutionary relationships of MHC class II-DRB sequences between African and Asian leopard subspecies are discussed.

MHC polymorphism and disease resistance to Vibrio anguillarum in 8 families of half-smooth tongue sole (Cynoglossus semilaevis)

Background

Genes in the major histocompatibility complex (MHC) have a critical role in both the innate and adaptive immune responses because of their involvement in presenting foreign peptides to T cells. However, the nature has remained largely unknown.

Results

We examined the genetic variation in MHC class IIB in half-smooth tongue sole (Cynoglossus semilaevis) after challenge with vibrio anguillarum. Two thousand and four hundred fry from 12 half-smooth tongue sole families were challenged with Vibrio anguillarum. To determine any association between alleles and resistance or susceptibility to V. anguillarum, 160 individuals from four high-resistance (HR, < 40.55% mortality) families and four low-resistance (LR, > 73.27% mortality) families were selected for MHC IIB exon2 gene sequence analysis. The MHC IIB exon2 genes of tongue sole displayed a high level of polymorphism and were discovered at least four loci. Meanwhile, the d_N/d_S [the ratio of non-synonymous (d_N) substitutions to synonymous (d_S) substitutions] in the peptide-binding region (PBR) was higher than that in the non-peptide-binding region (non-PBR). Eighty-eight alleles were discovered among 160 individuals, and 13 out of 88 alleles were used to analyze the distribution pattern between the resistant and susceptible families. Certain alleles presented in HR and LR with a different frequency, while other alleles were discovered in only the HR or LR families, not both. Five alleles, Cyse-DBB*6501, Cyse-DBB*4002, Cyse-DBB*6102, Cyse-DBB*5601 and Cyse-DBB*2801, were found to be associated with susceptibility to V. anguillarum with a frequency of 1.25%, 1.25%, 1.25%, 1.25% and 2.5% in the HR families, and 35%, 33.75%, 27.5%, 16.25%, 15% in the LR families (p < 0.01, 0.01, 0.01, 0.01, 0.01), respectively. Four alleles, Cyse-DBB*3301, Cyse-DBB*4701, Cyse-DBB*6801 and Cyse-DBB*5901, were found to be associated with resistance to V. anguillarum, with a frequency of 13.75%, 11.25%, 11.25%, 8.75% in the HR families and 1.25%, 1.25%, 1.25%, 1.25% and 1.25% in the LR families (p < 0.01, 0.05, 0.05 and p = 0.064), respectively.

Conclusions

Elucidation of the role of MHC II B genes in half-smooth tongue sole should prove to be helpful to the in-depth development of marker-assisted selective breeding in half-smooth tongue sole.

Polymorphism and Balancing Selection of MHC Class II DAB Gene in 7 Selective Flounder (Paralichthys olivaceus) Families

In order to determine the genetic variation of the MHC class IIB exon2 allele in the offspring, 700 fry from seven families of Japanese flounder challenged with V. anguillarum were studied, and different mortality rates were found in those families. Five to ten surviving and dead fry from each of the seven families were selected to study the MHC class II B exon2 gene with PCR and a direct sequencing method. One hundred and sixteen different exon2 sequences were found and 116 different alleles were identified, while a minimum of four loci were revealed in the MHC class II B exon2 gene. The ratio (d_N/d_S) of nonsynonymous substitution (d_N) to synonymous substitutions (d_S) in the peptide-binding region (PBR) of the MHC class IIB gene was 6.234, which indicated that balancing selection is acting on the MHC class IIB genes. The MHC IIB alleles were thus being passed on to their progeny. Some alleles were significantly more frequent in surviving than dead individuals. All together our data suggested that the alleles Paol-DAB*4301, Paol-DAB*4601, Paol-DAB*4302, Paol-DAB*3803, and Paol-DAB*4101 were associated with resistance to V. anguillarum in flounder.

Characterization of MHC-I in the blue tit (Cyanistes caeruleus) reveals low levels of genetic diversity and trans-population evolution across European populations

The major histcompatibility complex (MHC) is a vital component of the adaptive immune system in all vertebrates. This study is the first to characterize MHC class I (MHC-I) in blue tits (Cyanistes caeruleus), and we use MHC-I exon 3 sequence data from individuals originating from three locations across Europe: Spain, the Netherlands to Sweden. Our phylogeny of the 17 blue tit MHC-I alleles contains one allele cluster with low nucleotide diversity compared to the remaining more diverse alleles. We found a significant evidence for balancing selection in the peptide-binding region in the diverse allele group only. No separation according to geographic location was found in the phylogeny of alleles. Although the number of MHC-I loci of the blue tit is comparable to that of other passerine species, the nucleotide diversity of MHC-I appears to be much lower than that of other passerine species, including the closely related great tit (Parus major) and the severely inbred Seychelles warbler (Acrocephalus sechellensis). We believe that this initial MHC-I characterization in blue tits provides an important step towards understanding the mechanisms shaping MHC-I diversity in natural populations.

DR haplotype diversity of the cynomolgus macaque as defined by its transcriptome

The DR region of particular primate species may display allelic polymorphism and gene copy number variation (region configuration polymorphism). The sum of these distinct types of polymorphism is defined as complexity. To date, however, the DR region of cynomolgus macaques (Macaca fascicularis) has been poorly defined. Transcriptome analysis of a pedigreed colony, comprising animals from Indonesia and Indochina, revealed a total of 15 Mafa-DRA and 57 DRB alleles, specifying 28 different region configurations. The DRA alleles can be divided into two distinct lineages. One lineage is polymorphic, but the majority of the amino acid replacements map to the leader peptide. The second lineage is at best oligomorphic, and segregates with one specific Mafa-DRB allele. The number of Mafa-DRB genes ranges from two to five per haplotype. Due to the presence of pseudogenes, however, each haplotype encodes only one to three bona fide DRB transcripts. Depending on the region configuration in which the Mafa-DRB gene is embedded, identical alleles may display differential transcription levels. Region configurations appear to have been generated by recombination-like events. When genes or gene segments are relocated, it seems plausible that they may be placed in the context of distinct transcription control elements. As such, DRB region-related transcription level differences may add an extra layer of polymorphism to this section of the adaptive immune system.

On characterizing adaptive events unique to modern humans

Ever since the first draft of the human genome was completed in 2001, there has been increased interest in identifying genetic changes that are uniquely human, which could account for our distinct morphological and cognitive capabilities with respect to other apes. Recently, draft sequences of two extinct hominin genomes, a Neanderthal and Denisovan, have been released. These two genomes provide a much greater resolution to identify human-specific genetic differences than the chimpanzee, our closest extant relative. The Neanderthal genome paper presented a list of regions putatively targeted by positive selection around the time of the human–Neanderthal split. We here seek to characterize the evolutionary history of these candidate regions—examining evidence for selective sweeps in modern human populations as well as for accelerated adaptive evolution across apes. Results indicate that 3 of the top 20 candidate regions show evidence of selection in at least one modern human population (P < 5 × 10⁵). Additionally, four genes within the top 20 regions show accelerated amino acid substitutions across multiple apes (P < 0.01), suggesting importance across deeper evolutionary time. These results highlight the importance of evaluating evolutionary processes across both recent and ancient evolutionary timescales and intriguingly suggest a list of candidate genes that may have been uniquely important around the time of the human–Neanderthal split.

Effects of natural selection and gene conversion on the evolution of human glycophorins coding for MNS blood polymorphisms in malaria-endemic African populations

Malaria has been a very strong selection pressure in recent human evolution, particularly in Africa. Of the one million deaths per year due to malaria, more than 90% are in sub-Saharan Africa, a region with high levels of genetic variation and population substructure. However, there have been few studies of nucleotide variation at genetic loci that are relevant to malaria susceptibility across geographically and genetically diverse ethnic groups in Africa. Invasion of erythrocytes by Plasmodium falciparum parasites is central to the pathology of malaria. Glycophorin A (GYPA) and B (GYPB), which determine MN and Ss blood types, are two major receptors that are expressed on erythrocyte surfaces and interact with parasite ligands. We analyzed nucleotide diversity of the glycophorin gene family in 15 African populations with different levels of malaria exposure. High levels of nucleotide diversity and gene conversion were found at these genes. We observed divergent patterns of genetic variation between these duplicated genes and between different extracellular domains of GYPA. Specifically, we identified fixed adaptive changes at exons 3-4 of GYPA. By contrast, we observed an allele frequency spectrum skewed toward a significant excess of intermediate-frequency alleles at GYPA exon 2 in many populations; the degree of spectrum distortion is correlated with malaria exposure, possibly because of the joint effects of gene conversion and balancing selection. We also identified a haplotype causing three amino acid changes in the extracellular domain of glycophorin B. This haplotype might have evolved adaptively in five populations with high exposure to malaria.

A comprehensive study of polymorphic sites along the HLA-G gene: implication for gene regulation and evolution

HLA-G molecule plays an important role on immune response regulation and has been implicated on the inhibition of T and natural killer cell cytolytic function and inhibition of allogeneic T-cell proliferation. Due to its immune-modulator properties, the HLA-G gene expression has been associated with the outcome of allograft and of autoimmune, infectious, and malignant disorders. Several lines of evidence indicate that HLA-G polymorphisms at the 5'-upstream regulatory region (5' URR) and 3'-untranslated region (3' UTR) may influence the HLA-G expression levels. Because Brazilians represent one of the most heterogeneous populations in the world with the widest HLA-G coding region variability already detected among the studied populations, a high level of variability and haplotype diversity would be expected in Brazilians. On this basis, the 5' URR, coding, and 3' UTR variability were evaluated in a Brazilian series consisting of 100 healthy bone marrow donors, as well as the linkage disequilibrium pattern along the gene and the extended haplotypes encompassing several gene segment variations. The HLA-G locus seems to present six different HLA-G lineages showing functional variations mainly in nucleotides of the regulatory regions. Differences were observed at the 5' URR in positions that either coincide with or are close to transcription factor-binding sites and at the 3' UTR mainly in positions that have already been reported to influence HLA-G mRNA availability. We report several lines of evidence for balancing selection acting on the regulatory regions, which may indicate that these HLA-G lineages may be related to the differential HLA-G expression profiles.

Immune recognition of transplacentally acquired lymphoid allografts selects for increased major histocompatibility polymorphism

The extreme polymorphism of mammalian major histocompatibility (MHC) Class I and Class II alleles has been attributed to inbreeding avoidance, heterozygote advantage and pathogen driven selection for rare MHC alleles. However, MHC alleles can be classified into a limited number of allele supertypes based on the specificity of their peptide binding grooves (about 10 supertypes in the case of human MHC Class I alleles). The paradox is that if antigen presentation can be accomplished by a limited number of binding groove motifs, why are these loci so polymorphic? An unexplored driver of this complexity may be selection pressure to enhance the antigenicity and immune recognition of transplacentally acquired lymphoid allografts during pregnancy. The exchange of lymphoid cells between mother and fetus probably occurs in all pregnancies and may lead to fetal and/or maternal lymphoid microchimerism, a known cause of autoimmune disease. Natural selection may have favoured increased polymorphism at MHC Class I and Class II loci in order to improve immune surveillance of these cells and thereby reduce the incidence of maternal and fetal autoimmune disease. At the same time, selection may have favoured the retention of a limited set of allele supertypes which optimally present immunodominant antigens.

Reduced MHC and neutral variation in the Galápagos hawk, an island endemic

BACKGROUND

Genes at the major histocompatibility complex (MHC) are known for high levels of polymorphism maintained by balancing selection. In small or bottlenecked populations, however, genetic drift may be strong enough to overwhelm the effect of balancing selection, resulting in reduced MHC variability. In this study we investigated MHC evolution in two recently diverged bird species: the endemic Galápagos hawk (Buteo galapagoensis), which occurs in small, isolated island populations, and its widespread mainland relative, the Swainson's hawk (B. swainsoni).

RESULTS

We amplified at least two MHC class II B gene copies in each species. We recovered only three different sequences from 32 Galápagos hawks, while we amplified 20 unique sequences in 20 Swainson's hawks. Most of the sequences clustered into two groups in a phylogenetic network, with one group likely representing pseudogenes or nonclassical loci. Neutral genetic diversity at 17 microsatellite loci was also reduced in the Galápagos hawk compared to the Swainson's hawk.

CONCLUSIONS

The corresponding loss in neutral diversity suggests that the reduced variability present at Galápagos hawk MHC class II B genes compared to the Swainson's hawk is primarily due to a founder event followed by ongoing genetic drift in small populations. However, purifying selection could also explain the low number of MHC alleles present. This lack of variation at genes involved in the adaptive immune response could be cause for concern should novel diseases reach the archipelago.

Adaptive molecular evolution of the Major Histocompatibility Complex genes, DRA and DQA, in the genus Equus

BACKGROUND

Major Histocompatibility Complex (MHC) genes are central to vertebrate immune response and are believed to be under balancing selection by pathogens. This hypothesis has been supported by observations of extremely high polymorphism, elevated nonsynonymous to synonymous base pair substitution rates and trans-species polymorphisms at these loci. In equids, the organization and variability of this gene family has been described, however the full extent of diversity and selection is unknown. As selection is not expected to act uniformly on a functional gene, maximum likelihood codon-based models of selection that allow heterogeneity in selection across codon positions can be valuable for examining MHC gene evolution and the molecular basis for species adaptations.

RESULTS

We investigated the evolution of two class II MHC genes of the Equine Lymphocyte Antigen (ELA), DRA and DQA, in the genus Equus with the addition of novel alleles identified in plains zebra (E. quagga, formerly E. burchelli). We found that both genes exhibited a high degree of polymorphism and inter-specific sharing of allele lineages. To our knowledge, DRA allelic diversity was discovered to be higher than has ever been observed in vertebrates. Evidence was also found to support a duplication of the DQA locus. Selection analyses, evaluated in terms of relative rates of nonsynonymous to synonymous mutations (dN/dS) averaged over the gene region, indicated that the majority of codon sites were conserved and under purifying selection (dN <dS). However, the most likely evolutionary codon models allowed for variable rates of selection across codon sites at both loci and, at the DQA, supported the hypothesis of positive selection acting on specific sites.

CONCLUSIONS

Observations of elevated genetic diversity and trans-species polymorphisms supported the conclusion that balancing selection may be acting on these loci. Furthermore, at the DQA, positive selection was occurring at antigen binding sites, suggesting that a few selected residues may play a significant role in equid immune function. Future studies in natural equid populations will be valuable for understanding the functional significance of the uniquely diverse DRA locus and for elucidating the mechanism maintaining diversity at these MHC loci.

Computational prediction of MHC II-antigen binding supports divergent allele advantage and explains trans-species polymorphism

The major histocompatibility complex (MHC), coding for antigen presenting molecules of the adaptive immune system, represents one of the most polymorphic regions in the vertebrate genome. The exceptional polymorphism, which is potentially maintained by balancing selection under host-parasite coevolution, comprises excessive sequence divergence among alleles as well as ancient allelic lineages that predate species divergence (trans-species polymorphism). Here, the mechanisms that are proposed to maintain such sequence divergence and ancient lineages are investigated. Established computational antigen-binding prediction algorithms, which are based on empirical databases, are employed to determine the overlap in bound antigens among individual MHC class IIB alleles. The results show that genetically more divergent allele pairs experience less overlap and thus present a broader range of potential antigens. These findings support the divergent allele advantage hypothesis and furthermore suggest an evolutionary advantage explaining the maintenance of divergent allelic lineages, that is, trans-species polymorphism. In addressing a quantitative rather than qualitative aspect of MHC alleles, these insights highlight a new direction for future research on MHC evolution.

Immunogenetics as a tool in anthropological studies

The genes coding for the main molecules involved in the human immune system--immunoglobulins, human leucocyte antigen (HLA) molecules and killer-cell immunoglobulin-like receptors (KIR)--exhibit a very high level of polymorphism that reveals remarkable frequency variation in human populations. 'Genetic marker' (GM) allotypes located in the constant domains of IgG antibodies have been studied for over 40 years through serological typing, leading to the identification of a variety of GM haplotypes whose frequencies vary sharply from one geographic region to another. An impressive diversity of HLA alleles, which results in amino acid substitutions located in the antigen-binding region of HLA molecules, also varies greatly among populations. The KIR differ between individuals according to both gene content and allelic variation, and also display considerable population diversity. Whereas the molecular evolution of these polymorphisms has most likely been subject to natural selection, principally driven by host-pathogen interactions, their patterns of genetic variation worldwide show significant signals of human geographic expansion, demographic history and cultural diversification. As current developments in population genetic analysis and computer simulation improve our ability to discriminate among different--either stochastic or deterministic--forces acting on the genetic evolution of human populations, the study of these systems shows great promise for investigating both the peopling history of modern humans in the time since their common origin and human adaptation to past environmental (e.g. pathogenic) changes. Therefore, in addition to mitochondrial DNA, Y-chromosome, microsatellites, single nucleotide polymorphisms and other markers, immunogenetic polymorphisms represent essential and complementary tools for anthropological studies.

Cheetah paradigm revisited: MHC diversity in the world's largest free-ranging population

For more than two decades, the cheetah (Acinonyx jubatus) has been considered a paradigm of disease vulnerability associated with low genetic diversity, particularly at the immune genes of the major histocompatibility complex (MHC). Cheetahs have been used as a classic example in numerous conservation genetics textbooks as well as in many related scientific publications. However, earlier studies used methods with low resolution to quantify MHC diversity and/or small sample sizes. Furthermore, high disease susceptibility was reported only for captive cheetahs, whereas free-ranging cheetahs show no signs of infectious diseases and a good general health status. We examined whether the diversity at MHC class I and class II-DRB loci in 149 Namibian cheetahs was higher than previously reported using single-strand conformation polymorphism analysis, cloning, and sequencing. MHC genes were examined at the genomic and transcriptomic levels. We detected ten MHC class I and four class II-DRB alleles, of which nine MHC class I and all class II-DRB alleles were expressed. Phylogenetic analyses and individual genotypes suggested that the alleles belong to four MHC class I and three class II-DRB putative loci. Evidence of positive selection was detected in both MHC loci. Our study indicated that the low number of MHC class I alleles previously observed in cheetahs was due to a smaller sample size examined. On the other hand, the low number of MHC class II-DRB alleles previously observed in cheetahs was further confirmed. Compared with other mammalian species including felids, cheetahs showed low levels of MHC diversity, but this does not seem to influence the immunocompetence of free-ranging cheetahs in Namibia and contradicts the previous conclusion that the cheetah is a paradigm species of disease vulnerability.

Polymorphism and gene organization of water buffalo MHC-DQB genes show homology to the BoLA DQB region

In cattle (Bos taurus), there is evidence of more than 50 alleles of BoLA-DQB (bovine lymphocyte antigen DQB) that are distributed across at least five DQB loci, making this region one of the most complex in the BoLA gene family. In this study, DQB alleles were analysed for the water buffalo (Bubalus bubalis), another economically important bovine species. Twelve alleles for Bubu-DQB (Bubalis bubalis DQB) were determined by nucleotide sequence analysis. A phylogenetic analysis revealed numerous trans-species polymorphisms, with alleles from water buffalo assigned to at least three different loci (BoLA-DQB1, BoLA-DQB3 and BoLA-DQB4) that are also found in cattle. These presumptive loci were analysed for patterns of synonymous (d(S)) and non-synonymous (d(N)) substitution. Like BoLA-DQB1, Bubu-DQB1 was observed to be under strong positive selection for polymorphism. We conclude that water buffalo and cattle share the current arrangement of their DQB region because of their common ancestry.

HLA DNA sequence variation among human populations: molecular signatures of demographic and selective events

Molecular differences between HLA alleles vary up to 57 nucleotides within the peptide binding coding region of human Major Histocompatibility Complex (MHC) genes, but it is still unclear whether this variation results from a stochastic process or from selective constraints related to functional differences among HLA molecules. Although HLA alleles are generally treated as equidistant molecular units in population genetic studies, DNA sequence diversity among populations is also crucial to interpret the observed HLA polymorphism. In this study, we used a large dataset of 2,062 DNA sequences defined for the different HLA alleles to analyze nucleotide diversity of seven HLA genes in 23,500 individuals of about 200 populations spread worldwide. We first analyzed the HLA molecular structure and diversity of these populations in relation to geographic variation and we further investigated possible departures from selective neutrality through Tajima's tests and mismatch distributions. All results were compared to those obtained by classical approaches applied to HLA allele frequencies.Our study shows that the global patterns of HLA nucleotide diversity among populations are significantly correlated to geography, although in some specific cases the molecular information reveals unexpected genetic relationships. At all loci except HLA-DPB1, populations have accumulated a high proportion of very divergent alleles, suggesting an advantage of heterozygotes expressing molecularly distant HLA molecules (asymmetric overdominant selection model). However, both different intensities of selection and unequal levels of gene conversion may explain the heterogeneous mismatch distributions observed among the loci. Also, distinctive patterns of sequence divergence observed at the HLA-DPB1 locus suggest current neutrality but old selective pressures on this gene. We conclude that HLA DNA sequences advantageously complement HLA allele frequencies as a source of data used to explore the genetic history of human populations, and that their analysis allows a more thorough investigation of human MHC molecular evolution.

Contrasting responses to selection in class I and class IIα major histocompatibility-linked markers in salmon

Comparison of levels and patterns of genetic variation in natural populations either across loci or against neutral expectation can yield insight into locus-specific differences in the strength and direction of evolutionary forces. We used both approaches to test the hypotheses on patterns of selection on major histocompatibility (MH)-linked markers. We performed temporal analyses of class I and class IIα MH-linked markers and eight microsatellite loci in two Atlantic salmon populations in Ireland on two temporal scales: over six decades and 9 years in the rivers Burrishoole and Delphi, respectively. We also compared contemporary Burrishoole and Delphi samples with nearby populations for the same loci. On comparing patterns of temporal and spatial differentiation among classes of loci, the class IIα MH-linked marker was consistently identified as an outlier compared with patterns at the other microsatellite loci or neutral expectation. We found higher levels of temporal and spatial heterogeneity in heterozygosity (but not in allelic richness) for the class IIα MH-linked marker compared with microsatellites. Tests on both within- and among-population differentiation are consistent with directional selection acting on the class IIα-linked marker in both temporal and spatial comparisons, but only in temporal comparisons for the class I-linked marker. Our results indicate a complex pattern of selection on MH-linked markers in natural populations of Atlantic salmon. These findings highlight the importance of considering selection on MH-linked markers when using these markers for management and conservation purposes.

MHC class II DRB diversity, selection pattern and population structure in a neotropical bat species, Noctilio albiventris

Genes of the major histocompatibility complex (MHC) have a crucial role in the immune response of vertebrates, alter the individual odour and are involved in shaping mating preferences. Pathogen-mediated selection, sexual selection and maternal-fetal interactions have been proposed as the main drivers of frequently observed high levels of polymorphism in functionally important parts of the MHC. Bats constitute the second largest mammalian order and have recently emerged as important vectors of infectious diseases. In addition, Chiroptera are interesting study subjects in evolutionary ecology in the context of olfactory communication, mate choice and associated fitness benefits. Thus, it is surprising that they belong to the least studied mammalian taxa in terms of their MHC diversity. In this study, we investigated the variability in the functionally important MHC class II gene DRB, evidence for selection and population structure in the group-living lesser bulldog bat, Noctilio albiventris, in Panama. We found a single expressed, polymorphic Noal-DRB gene. The substitution pattern of the nucleotide sequences of the 18 detected alleles provided evidence for positive selection acting above the evolutionary history of the species in shaping MHC diversity. Roosting colonies were not genetically differentiated but females showed lower levels of heterozygosity than males, which might be a sign that the sexes differ in the selection pressures acting on the MHC. This study provides the prerequisites for further investigations of the role of the individual MHC constitution in parasite resistance, olfactory communication and mate choice in N. albiventris and other bats.

Allelic polymorphism, gene duplication and balancing selection of the MHC class II DAB gene of Cynoglossus semilaevis (Cynoglossidae)

Major histocompatibility complex (MHC) genes play an important role in the immune response of vertebrates. Allelic polymorphism and evolutionary mechanism of MHC genes have been investigated in many mammals, but much less is known in teleosts. We examined the polymorphism, gene duplication and balancing selection of the MHC class II DAB gene of the half-smooth tongue sole (Cynoglossus semilaevis); 23 alleles were found in this species. Gene duplication manifested as three to six distinct sequences at each domain in the same individuals. Non-synonymous substitutions occurred at a significantly higher frequency than synonymous substitutions in the PBR domain, suggesting balancing selection for maintaining polymorphisms at the MHC II DAB locus. Many positive selection sites were found to act very intensely on antigen-binding sites of MHC class II DAB gene.

Geographic variation of the major histocompatibility complex in Eastern Atlantic grey seals (Halichoerus grypus)

Pathogen-driven balancing selection maintains high genetic diversity in many vertebrates, particularly in the major histocompatibility complex (MHC) immune system gene family, which is often associated with disease susceptibility. In large natural populations where subpopulations face different pathogen pressures, the MHC should show greater genetic differentiation within a species than neutral markers. We examined genetic diversity at the MHC-DQB locus and nine putatively neutral microsatellite markers in grey seals (Halichoerus grypus) from eight United Kingdom (UK) colonies, the Faeroe Islands and Sable Island, Canada. Five DQB alleles were identified in grey seals, which varied in prevalence across the grey seal range. Among the seal colonies, significant differences in DQB allele and haplotype frequencies and in average DQB heterozygosity were observed. Additionally, the DQB gene exhibited greater differentiation among colonies compared with neutral markers, yet a weaker pattern of isolation by distance (IBD). After correcting for the underlying IBD pattern, subpopulations breeding in similar habitats were more similar to one another in DQB allele frequencies than populations breeding in different habitats, but the same did not hold true for microsatellites, suggesting that habitat-specific pathogen pressure influences MHC evolution. Overall, the data are consistent with selection at MHC-DQB loci in grey seals with both varying selective pressures and geographic population structure appearing to influence the DQB genetic composition of breeding colonies.

Balancing selection maintains a form of ERAP2 that undergoes nonsense-mediated decay and affects antigen presentation

A remarkable characteristic of the human major histocompatibility complex (MHC) is its extreme genetic diversity, which is maintained by balancing selection. In fact, the MHC complex remains one of the best-known examples of natural selection in humans, with well-established genetic signatures and biological mechanisms for the action of selection. Here, we present genetic and functional evidence that another gene with a fundamental role in MHC class I presentation, endoplasmic reticulum aminopeptidase 2 (ERAP2), has also evolved under balancing selection and contains a variant that affects antigen presentation. Specifically, genetic analyses of six human populations revealed strong and consistent signatures of balancing selection affecting ERAP2. This selection maintains two highly differentiated haplotypes (Haplotype A and Haplotype B), with frequencies 0.44 and 0.56, respectively. We found that ERAP2 expressed from Haplotype B undergoes differential splicing and encodes a truncated protein, leading to nonsense-mediated decay of the mRNA. To investigate the consequences of ERAP2 deficiency on MHC presentation, we correlated surface MHC class I expression with ERAP2 genotypes in primary lymphocytes. Haplotype B homozygotes had lower levels of MHC class I expressed on the surface of B cells, suggesting that naturally occurring ERAP2 deficiency affects MHC presentation and immune response. Interestingly, an ERAP2 paralog, endoplasmic reticulum aminopeptidase 1 (ERAP1), also shows genetic signatures of balancing selection. Together, our findings link the genetic signatures of selection with an effect on splicing and a cellular phenotype. Although the precise selective pressure that maintains polymorphism is unknown, the demonstrated differences between the ERAP2 splice forms provide important insights into the potential mechanism for the action of selection.

It has long been known that the extremely high levels of genetic diversity present in the major histocompatibility locus (MHC) are due to balancing selection, a type of natural selection that maintains advantageous genetic diversity in populations. The MHC encodes for molecules required for a type of antigen presentation that mediates detection of infected and cancerous cells by the immune system; the genetic diversity of the MHC thus ensures an adequate response to the wide variety of pathogens that humans encounter. Here, we show that other genes involved in the same antigen-presentation pathway are also subject to balancing selection in humans. Specifically, we show that balancing selection acts to maintain two forms of the endoplasmic reticulum aminopeptidase 2 gene (ERAP2), which encodes a protein also involved in antigen presentation. Although the two ERAP2 forms are present in a similar frequency (close to 0.5), they are associated with differences with respect to the levels of MHC molecules on the cell surface of immune cells. In summary, our findings show that natural selection maintains variants of ERAP2 that affect immune surveillance; they also establish ERAP2 as one of the few examples of balancing selection in humans where the selected variant, its functional consequences, and its influence in interpersonal diversity are known.

MHC class IIalpha gene polymorphism and its association with resistance/susceptibility to Vibrio anguillarum in Japanese flounder (Paralichthys olivaceus)

Association between polymorphism of the major histocompatibility complex (MHC) genes and disease resistance has been documented for few teleosts. In this study, we first investigate the genetic variation at the MHC IIalpha gene in Japanese flounder (Paralichthys olivaceus) for survival after challenge with bacterial infection. To explore the specific allele associated with disease resistance, about 6000 individuals from 60 families challenged with Vibrio anguillarum, which causes significantly different mortality in flounder families. 15-20 individuals from each of six high-resistance (HR) and six low-resistance (LR) families were screened for their MHC class IIalpha genotypes using sequence analysis. High polymorphism of MHC IIalpha gene and at least two loci were discovered in Japanese flounder and the rate of d(N) occurred at a significantly higher frequency than that of d(S) in PBR and non-PBR, the balancing selection hypothesis could explain the high polymorphism of MHC IIalpha gene in flounder. From the 212 individuals, a total of 55 MHC class IIalpha alleles were identified, and 9 alleles were used to study association between alleles and resistance to disease. Three alleles, Paol-DAA*1301, Paol-DAA*1401 and Paol-DAA*2201 were significantly associated with resistance against V. anguillarum, and Paol-DAA*1001 and Paol-DAA*1501 alleles were significantly associated with increased susceptibility to V. anguillarum. This study confirmed the association between alleles of MHC class IIalpha gene and disease resistance or susceptibility to bacterial infection in flounder, and the disease resistance-related MHC markers could be used for molecular marker-assisted selective breeding in the flounder.

Polymorphisms in Anopheles gambiae immune genes associated with natural resistance to Plasmodium falciparum

Many genes involved in the immune response of Anopheles gambiae, the main malaria vector in Africa, have been identified, but whether naturally occurring polymorphisms in these genes underlie variation in resistance to the human malaria parasite, Plasmodium falciparum, is currently unknown. Here we carried out a candidate gene association study to identify single nucleotide polymorphisms (SNPs) associated with natural resistance to P. falciparum. A. gambiae M form mosquitoes from Cameroon were experimentally challenged with three local wild P. falciparum isolates. Statistical associations were assessed between 157 SNPs selected from a set of 67 A. gambiae immune-related genes and the level of infection. Isolate-specific associations were accounted for by including the effect of the isolate in the analysis. Five SNPs were significantly associated to the infection phenotype, located within or upstream of AgMDL1, CEC1, Sp PPO activate, Sp SNAKElike, and TOLL6. Low overall and local linkage disequilibrium indicated high specificity in the loci found. Association between infection phenotype and two SNPs was isolate-specific, providing the first evidence of vector genotype by parasite isolate interactions at the molecular level. Four SNPs were associated to either oocyst presence or load, indicating that the genetic basis of infection prevalence and intensity may differ. The validity of the approach was verified by confirming the functional role of Sp SNAKElike in gene silencing assays. These results strongly support the role of genetic variation within or near these five A. gambiae immune genes, in concert with other genes, in natural resistance to P. falciparum. They emphasize the need to distinguish between infection prevalence and intensity and to account for the genetic specificity of vector-parasite interactions in dissecting the genetic basis of Anopheles resistance to human malaria.

Anopheles gambiae is the main malaria vector in Africa, transmitting the parasite when it blood feeds on human hosts. The parasite undergoes several developmental stages in the mosquito to complete its life cycle, during which time it is confronted by the mosquito's immune system. The resistance of mosquitoes to malaria infection is highly variable in wild populations and is known to be under strong genetic control, but to date the specific genes responsible for this variation remain to be identified. The present study uncovers variations in A. gambiae immune genes that are associated with natural resistance to Plasmodium falciparum, the deadliest human malaria parasite. The association of some mosquito genetic loci with the level of infection depended on the P. falciparum isolate, suggesting that resistance is determined by interactions between the genome of the mosquito and that of the parasite. This finding highlights the need to account for the natural genetic diversity of malaria parasites in future research on vector-parasite interactions. The loci uncovered in this study are potential targets for developing novel malaria control strategies based on natural mosquito resistance mechanisms.

MHC class I and MHC class II DRB gene variability in wild and captive Bengal tigers (Panthera tigris tigris)

Bengal tigers are highly endangered and knowledge on adaptive genetic variation can be essential for efficient conservation and management. Here we present the first assessment of allelic variation in major histocompatibility complex (MHC) class I and MHC class II DRB genes for wild and captive tigers from India. We amplified, cloned, and sequenced alpha-1 and alpha-2 domain of MHC class I and beta-1 domain of MHC class II DRB genes in 16 tiger specimens of different geographic origin. We detected high variability in peptide-binding sites, presumably resulting from positive selection. Tigers exhibit a low number of MHC DRB alleles, similar to other endangered big cats. Our initial assessment-admittedly with limited geographic coverage and sample size-did not reveal significant differences between captive and wild tigers with regard to MHC variability. In addition, we successfully amplified MHC DRB alleles from scat samples. Our characterization of tiger MHC alleles forms a basis for further in-depth analyses of MHC variability in this illustrative threatened mammal.

Gene duplication and evidence for balancing selection acting on MHC class II DAA gene of the half-smooth tongue sole (Cynoglossus semilaevis)

Allelic polymorphism and evolution mechanism of major histocompatibility complex (MHC) genes has been investigated in many mammals, however, much less is known in teleost. In order to investigate the mechanisms creating and maintaining variability at the MHC class II DAA locus, we examined the polymorphism, gene duplication and balancing selection of MHC class II DAA gene of the half-smooth tongue sole (Cynoglossus semilaevis). We described 33 alleles in the C. semilaevis, recombination and gene duplication seems to play more important roles in the origin of new alleles. The rate of non-synonymous substitutions (d(N)) occurred at a significantly higher frequency than that of synonymous substitutions (d(S)) in peptide-binding region (PBR) and non-PBR, suggesting balancing selection for maintaining polymorphisms at the MHC II DAA locus. Many positive selection sites were found to act very intensively on antigen-binding sites. Our founding suggests a snapshot in an evolutionary process of MHC-DAA gene evolution of the C. semilaevis.

Towards population genomics of effector-effector target interactions

Pathogen-plant host coevolutionary interactions exert strong natural selection on both organisms, specifically on the genes coding for effectors (pathogens), as well as on those coding for effector targets and R proteins (plant hosts). Natural selection leaves behind DNA sequence signatures on such genes and on linked genomic regions. These signatures can readily be detected by studying the patterns of intraspecies polymorphisms and interspecies divergence of the DNA sequences. Recent developments in DNA sequencing technology have made whole-genome studies on patterns of DNA polymorphisms : divergence possible. This type of analysis, called 'population genomics', appears to be powerful enough to identify novel effector-effector target genes. Here, we provide an overview of the statistical tools used for population genomics and their applications. This is followed by a brief review of evolutionary studies on plant genes involved in host-pathogen interactions. Finally we provide an example from our study on Magnaporthe oryzae.

Identification of the expressed MHC class II DQB gene of the Asiatic black bear, Ursus thibetanus, in Japan

Genetic diversity estimation of the major histocompatibility complex (MHC) gene may be an important tool in the assessment of immune response ability against infectious disease. We were able to identify a near full-length expressed DQB sequence by RACE-PCR method from the Asiatic black bear, Ursus thibetanus in Japan. This is the first such full length expression in the Ursidae. The bear had at least one functional DQB locus. In phylogenetic tree analysis its DQB amino acid sequence formed a monophyletic group with DQB sequences from members of the order Carnivora and had a 90% nucleotide sequence similarity with the DQB allele of the California sea lion, Zalophus californianus. We compared the DQB amino acid composition of U. thibetanus with those of several other mammalian species including Homo sapiens. Amino acid residues known to be functionally important for human MHC genes, tended to be also conserved among other mammalian species while PBRs in the beta1 domain were heterogeneous among mammalian species. The DQB sequence obtained from the bear had not only no putative frameshifts or deletions but also no abnormal amino acid mutations such as had been observed in human DQB molecules. This suggests that the bear DQB sequence was an apparently functional DQB allele. As a preliminary study, we sequenced the exon 2 region of DQB alleles from genomic DNA, and succeeded to amplify the exon 2 of DQB loci. Our study will provide useful information for conservation genetics of the U. thibetanus as well as more generally regarding the mammalian MHC region.

Allelic diversity of the MHC class II DRB genes in brown bears (Ursus arctos) and a comparison of DRB sequences within the family Ursidae

The allelic diversity of the DRB locus in major histocompatibility complex (MHC) genes was analyzed in the brown bear (Ursus arctos) from the Hokkaido Island of Japan, Siberia, and Kodiak of Alaska. Nineteen alleles of the DRB exon 2 were identified from a total of 38 individuals of U. arctos and were highly polymorphic. Comparisons of non-synonymous and synonymous substitutions in the antigen-binding sites of deduced amino acid sequences indicated evidence for balancing selection on the bear DRB locus. The phylogenetic analysis of the DRB alleles among three genera (Ursus, Tremarctos, and Ailuropoda) in the family Ursidae revealed that DRB allelic lineages were not separated according to species. This strongly shows trans-species persistence of DRB alleles within the Ursidae.

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.

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.

Sympatric and allopatric divergence of MHC genes in threespine stickleback

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

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.

The GimA locus of extraintestinal pathogenic E. coli: does reductive evolution correlate with habitat and pathotype?

IbeA (invasion of brain endothelium), which is located on a genomic island termed GimA, is involved in the pathogenesis of several extraintestinal pathogenic E. coli (ExPEC) pathotypes, including newborn meningitic E. coli (NMEC) and avian pathogenic E. coli (APEC). To unravel the phylogeny of GimA and to investigate its island character, the putative insertion locus of GimA was determined via Long Range PCR and DNA-DNA hybridization in 410 E. coli isolates, including APEC, NMEC, uropathogenic (UPEC), septicemia-associated E. coli (SEPEC), and human and animal fecal isolates as well as in 72 strains of the E. coli reference (ECOR) collection. In addition to a complete GimA (∼20.3 kb) and a locus lacking GimA we found a third pattern containing a 342 bp remnant of GimA in this strain collection. The presence of GimA was almost exclusively detected in strains belonging to phylogenetic group B2. In addition, the complete GimA was significantly more frequent in APEC and NMEC strains while the GimA remnant showed a higher association with UPEC strains. A detailed analysis of the ibeA sequences revealed the phylogeny of this gene to be consistent with that obtained by Multi Locus Sequence Typing of the strains. Although common criteria for genomic islands are partially fulfilled, GimA rather seems to be an ancestral part of phylogenetic group B2, and it would therefore be more appropriate to term this genomic region GimA locus instead of genomic island. The existence of two other patterns reflects a genomic rearrangement in a reductive evolution-like manner.

Testing for selection on synonymous sites in plant mitochondrial DNA: the role of codon bias and RNA editing

Since plant mitochondrial genomes exhibit some of the slowest known synonymous substitution rates, it is generally believed that they experience exceptionally low mutation rates. However, the use of synonymous substitution rates to infer mutation rates depends on the implicit assumption that synonymous sites are evolving neutrally (or nearly so). To assess the validity of this assumption in plant mitochondrial genomes, we examined coding sequence for footprints of selection acting at synonymous sites. We found that synonymous sites exhibit an AT rich and pyrimidine skewed nucleotide composition compared to both non-synonymous sites and non-coding regions. We also found some evidence for selection associated with both biased codon usage and conservation of regulatory sequences involved in mRNA processing, although some of these findings are subject to alternative non-adaptive interpretations. Regardless, the inferred strength of selection appears too weak to account for the variation in substitution rates between the mitochondrial genomes of plants and other multicellular eukaryotes. Therefore, these results are consistent with the interpretation that plant mitochondrial genomes experience a substantially lower mutation rate rather than increased functional constraints acting on synonymous sites. Nevertheless, there are important nucleotide composition patterns (particularly the differences between synonymous sites and non-coding DNA) that remain largely unexplained.

Gene duplication and fragmentation in the zebra finch major histocompatibility complex

BACKGROUND

Due to its high polymorphism and importance for disease resistance, the major histocompatibility complex (MHC) has been an important focus of many vertebrate genome projects. Avian MHC organization is of particular interest because the chicken Gallus gallus, the avian species with the best characterized MHC, possesses a highly streamlined minimal essential MHC, which is linked to resistance against specific pathogens. It remains unclear the extent to which this organization describes the situation in other birds and whether it represents a derived or ancestral condition. The sequencing of the zebra finch Taeniopygia guttata genome, in combination with targeted bacterial artificial chromosome (BAC) sequencing, has allowed us to characterize an MHC from a highly divergent and diverse avian lineage, the passerines.

RESULTS

The zebra finch MHC exhibits a complex structure and history involving gene duplication and fragmentation. The zebra finch MHC includes multiple Class I and Class II genes, some of which appear to be pseudogenes, and spans a much more extensive genomic region than the chicken MHC, as evidenced by the presence of MHC genes on each of seven BACs spanning 739 kb. Cytogenetic (FISH) evidence and the genome assembly itself place core MHC genes on as many as four chromosomes with TAP and Class I genes mapping to different chromosomes. MHC Class II regions are further characterized by high endogenous retroviral content. Lastly, we find strong evidence of selection acting on sites within passerine MHC Class I and Class II genes.

CONCLUSION

The zebra finch MHC differs markedly from that of the chicken, the only other bird species with a complete genome sequence. The apparent lack of synteny between TAP and the expressed MHC Class I locus is in fact reminiscent of a pattern seen in some mammalian lineages and may represent convergent evolution. Our analyses of the zebra finch MHC suggest a complex history involving chromosomal fission, gene duplication and translocation in the history of the MHC in birds, and highlight striking differences in MHC structure and organization among avian lineages.

Olfactory receptor gene family evolution in stickleback and medaka fishes

Interaction of olfactory receptor (OR) genes with environmental odors is regarded as the first step of olfaction. In this study, OR genes of two fish, medaka (Oryzias latipes) and stickleback (Gasterosteus aculeatus), were identified and an evolutional analysis was conducted. The selection pressure of different TM regions and complete coding region were compared. Three TM regions (TM4, TM5 and TM6) were found to have higher average Ka/Ks values, which might be partly caused by positive selection as suggested by subsequent positive selection analysis. Further analysis showed that many PTSs overlap, or are adjacent to previously deduced binding sites in mammals. These results support the hypothesis that binding sites of fish OR genes may evolved under positive selection.

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.