Population biology of antigen presentation by MHC class I molecules

Soma-to-germline feedback is implied by the extreme polymorphism at IGHV relative to MHC: The manifest polymorphism of the MHC appears greatly exceeded at Immunoglobulin loci, suggesting antigen-selected somatic V mutants penetrate Weismann's Barrier

Soma-to-germline feedback is forbidden under the neo-Darwinian paradigm. Nevertheless, there is a growing realization it occurs frequently in immunoglobulin (Ig) variable (V) region genes. This is a surprising development. It arises from a most unlikely source in light of the exposure of co-author EJS to the haplotype data of RL Dawkins and others on the polymorphism of the Major Histocompatibility Complex, which is generally assumed to be the most polymorphic region in the genome (spanning ∼4 Mb). The comparison between the magnitude of MHC polymorphism with estimates for the human heavy chain immunoglobulin V locus (spanning ∼1 Mb), suggests IGHV could be many orders of magnitude more polymorphic than the MHC. This conclusion needs airing in the literature as it implies generational churn and soma-to-germline gene feedback. Pedigree-based experimental strategies to resolve the IGHV issue are outlined.

Heterozygote advantage at MHC DRB may influence response to infectious disease epizootics

The effect of MHC polymorphism on individual fitness variation in the wild remains equivocal; however, much evidence suggests that heterozygote advantage is a major determinant. To understand the contribution of MHC polymorphism to individual disease resistance or susceptibility in natural populations, we investigated two MHC class II B loci, DQB and DRB, in the New Zealand sea lion (NZSL, Phocarctos hookeri). The NZSL is a threatened species which is unusually susceptible to death by bacterial infection at an early age; it has suffered three bacterial induced epizootics resulting in high mortality levels of young pups since 1997. The MHC DQB and DRB haplotypes of dead NZSL pups with known cause of death (bacteria, enteritis or trauma) were sequenced and reconstructed, compared to pups that survived beyond 2 months of age, and distinct MHC DRB allele frequency and genotype differences were identified. Two findings were striking: (i) one DRB allele was present only in dead pups, and (ii) one heterozygous DRB genotype, common in live pups, was absent from dead pups. These results are consistent with some functional relationship with these variants and suggest heterozygote advantage is operating at DRB. We found no association between heterozygosity and fitness at 17 microsatellite loci, indicating that general heterozygosity is not responsible for the effect on fitness detected here. This result may be a consequence of recurrent selection by multiple pathogen assault over recent years and highlights the importance of heterozygote advantage at MHC as a potential mechanism for fitness differences in wild populations.

Lipid antigen presentation through CD1d pathway in mouse lung epithelial cells, macrophages and dendritic cells and its suppression by poly-dispersed single-walled carbon nanotubes

Effect of poly-dispersed acid-functionalized single-walled carbon nanotubes (AF-SWCNTs) was examined on lipid antigen presentation through CD1d pathway on three cell lines, LA4, MHS, and JAWSII used as prototype antigen presenting cells (APCs). CD1d molecule was expressed on 80-90% MHS (prototype macrophages) and JAWSII (prototype dendritic cells) cells whereas <5% LA4 cells (lung epithelial cells, non-classical APCs) expressed CD1d. Treatment with AF-SWCNTs but not with pristine SWCNTs resulted in a significant decline in the level of CD1d mRNA as well as mRNA levels of some other intracellular proteins involved in lipid antigen presentation pathway (MTP, ApoE, prosaposin, SR-BI and LDLr). Lipid antigen presentation was assessed by first incubating the cells with a prototype lipid antigen (α-Glactosylceramide or αGC) and then staining with L363 monoclonal antibody that detects αGC bound to CD1d molecule. While 100% MHS and JAWSII cells presented αGC, only 20% LA4 cells presented the CD1d antigen. Treatment with AF-SWCNTs resulted in a 30-40% decrease in αGC antigen presentation in all three cell lines. These results show that AF-SWCNT treatment down regulated the lipid antigen presentation pathway in all three cell lines and significantly lowered the ability of these cell lines to present αGC antigen.

Identification and characterization of a novel HLA-A allele, HLA-A*68:105, by genomic sequencing

A novel HLA-A allele, HLA-A*68:105, was detected by sequence-based typing (SBT) in an Italian bone marrow donor. It differs from HLA-A*68:01:02 at five nucleotides, three intronic, nt 699 T->G (intron 2), nt 705 T->C (intron 2) and nt 2770 G->A (intron 7), and two located in exon 3, at positions 726 A-G (codon 94 Ile->Val) and 733 T-G (codon 97 Arg->Met), respectively.

Variation and genetic control of gene expression in primary immunocytes across inbred mouse strains

To determine the breadth and underpinning of changes in immunocyte gene expression due to genetic variation in mice, we performed, as part of the Immunological Genome Project, gene expression profiling for CD4(+) T cells and neutrophils purified from 39 inbred strains of the Mouse Phenome Database. Considering both cell types, a large number of transcripts showed significant variation across the inbred strains, with 22% of the transcriptome varying by 2-fold or more. These included 119 loci with apparent complete loss of function, where the corresponding transcript was not expressed in some of the strains, representing a useful resource of "natural knockouts." We identified 1222 cis-expression quantitative trait loci (cis-eQTL) that control some of this variation. Most (60%) cis-eQTLs were shared between T cells and neutrophils, but a significant portion uniquely impacted one of the cell types, suggesting cell type-specific regulatory mechanisms. Using a conditional regression algorithm, we predicted regulatory interactions between transcription factors and potential targets, and we demonstrated that these predictions overlap with regulatory interactions inferred from transcriptional changes during immunocyte differentiation. Finally, comparison of these and parallel data from CD4(+) T cells of healthy humans demonstrated intriguing similarities in variability of a gene's expression: the most variable genes tended to be the same in both species, and there was an overlap in genes subject to strong cis-acting genetic variants. We speculate that this "conservation of variation" reflects a differential constraint on intraspecies variation in expression levels of different genes, either through lower pressure for some genes, or by favoring variability for others.

Genetic polymorphism of major histocompatibility complex class IIB alleles and pathogen resistance in the giant spiny frog Quasipaa spinosa

Major histocompatibility complex (MHC) genes are candidates for determining disease susceptibility due to their pivotal role in both innate and adaptive immune responses. Accordingly, the association between the genetic variation of MHC genes and the pathogen resistance has been investigated in numerous vertebrates. To date, however, little is reported in amphibians. In this study, we investigate the genetic variation at the MHC class IIB gene in the giant spiny frog Quasipaa spinosa, which has high commercial value in China. The full length of MHC class IIB cDNA was cloned from Q. spinosa by homology cloning and rapid amplification of cDNA end-polymerase chain reaction (RACE-PCR). Two MHC class IIB loci were identified in Q. spinosa. We also developed PCR primers for a portion of the second exon of the MHC class IIB gene. A total of 26 MHC class IIB alleles were identified. The dN rate was significantly higher than the dS rate in the putative peptide-binding region, thereby proving the positive selection hypothesis. In addition, individuals intraperitoneally injected with Aeromonas hydrophila were used to study the association between MHC class IIB alleles and pathogen resistance/susceptibility, to explore the specific alleles in balancing selection. Eighty frogs were used after exposure to A. hydrophila infection. Nine alleles were used to study the association between the alleles and disease resistance. Two alleles, namely, Pasa-DAB(∗)1301 and Pasa-DAB(∗)0901, were significantly associated with resistance against A. hydrophila. This study provides valuable information on the structure of the MHC class IIB gene and confirms the association between MHC class IIB gene alleles and disease resistance to bacterial infection in Q. spinosa. Moreover, pathogen resistance-related MHC markers can be used for the selective breeding of the giant spiny frog.

Advantageous diversity maintained by balancing selection in humans

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

Infectious disease, shifting climates, and opportunistic predators: cumulative factors potentially impacting wild salmon declines

Emerging diseases are impacting animals under high-density culture, yet few studies assess their importance to wild populations. Microparasites selected for enhanced virulence in culture settings should be less successful maintaining infectivity in wild populations, as once the host dies, there are limited opportunities to infect new individuals. Instead, moderately virulent microparasites persisting for long periods across multiple environments are of greatest concern. Evolved resistance to endemic microparasites may reduce susceptibilities, but as barriers to microparasite distributions are weakened, and environments become more stressful, unexposed populations may be impacted and pathogenicity enhanced. We provide an overview of the evolutionary and ecological impacts of infectious diseases in wild salmon and suggest ways in which modern technologies can elucidate the microparasites of greatest potential import. We present four case studies that resolve microparasite impacts on adult salmon migration success, impact of river warming on microparasite replication, and infection status on susceptibility to predation. Future health of wild salmon must be considered in a holistic context that includes the cumulative or synergistic impacts of multiple stressors. These approaches will identify populations at greatest risk, critically needed to manage and potentially ameliorate the shifts in current or future trajectories of wild populations.

Major histocompatibility complex class II polymorphisms in forest musk deer (Moschus berezovskii) and their probable association with purulent disease

Genes of the major histocompatibility complex (MHC) family are crucial in immune responses because they present pathogenic peptides to T cells. In this study, we analysed the genetic variation in forest musk deer (Moschus berezovskii) MHC II genes and its potential association with musk deer purulent disease. In total, 53 purulent disease-susceptible and 46 purulent disease-resistant individuals were selected for MHC II exon 2 fragment analysis. Among them, 16 DQ alleles and four additional DR alleles were identified, with DQ exon 2 fragments displaying a low level of polymorphism. The nonsynonymous substitutions exceeded the synonymous substitutions in the peptide-binding sites of DQA2, DQB1 and DQB2. Then, 28 MHC II alleles were used to analyse the distribution patterns of purulent disease between the susceptible and resistant groups. Among them, three alleles (DQA1*01, DQA1*02 and DQA2*04) were found to be resistant, and five alleles (DRB3*07, DQA1*03, DQA1*04, DQA2*05 and DQA2*06) were found to increase susceptibility. Additionally, three haplotypes were found to be putatively associated with musk deer purulent disease. However, these three haplotypes were only found in the resistant or susceptible group, and their frequencies were low. The results from our study support a contributory role of MHC II polymorphisms in the development of purulent disease in forest musk deer.

Stepwise threshold clustering: a new method for genotyping MHC loci using next-generation sequencing technology

Genes of the vertebrate major histocompatibility complex (MHC) are of great interest to biologists because of their important role in immunity and disease, and their extremely high levels of genetic diversity. Next generation sequencing (NGS) technologies are quickly becoming the method of choice for high-throughput genotyping of multi-locus templates like MHC in non-model organisms. Previous approaches to genotyping MHC genes using NGS technologies suffer from two problems:1) a "gray zone" where low frequency alleles and high frequency artifacts can be difficult to disentangle and 2) a similar sequence problem, where very similar alleles can be difficult to distinguish as two distinct alleles. Here were present a new method for genotyping MHC loci--Stepwise Threshold Clustering (STC)--that addresses these problems by taking full advantage of the increase in sequence data provided by NGS technologies. Unlike previous approaches for genotyping MHC with NGS data that attempt to classify individual sequences as alleles or artifacts, STC uses a quasi-Dirichlet clustering algorithm to cluster similar sequences at increasing levels of sequence similarity. By applying frequency and similarity based criteria to clusters rather than individual sequences, STC is able to successfully identify clusters of sequences that correspond to individual or similar alleles present in the genomes of individual samples. Furthermore, STC does not require duplicate runs of all samples, increasing the number of samples that can be genotyped in a given project. We show how the STC method works using a single sample library. We then apply STC to 295 threespine stickleback (Gasterosteus aculeatus) samples from four populations and show that neighboring populations differ significantly in MHC allele pools. We show that STC is a reliable, accurate, efficient, and flexible method for genotyping MHC that will be of use to biologists interested in a variety of downstream applications.

β2-Microglobulin-mediated signaling as a target for cancer therapy

β2-microglobulin (β2-m) has become the focus of intense scrutiny since the discovery of its undesirable roles promoting osteomimicry and cancer progression. β2-m is a well-known housekeeping protein that forms complexes with the heavy chain of major histocompatibility complex class I molecules, which are heterodimeric cell surface proteins that present antigenic peptides to cytotoxic T cells. On recognition of foreign peptide antigens on cell surfaces, T cells actively bind and lyse antigen-presenting cancer cells. In addition to its roles in tumor immunity, β2-m has two different functions in cancer cells, either tumor promoting or tumor suppressing, in cancer cell context-dependent manner. Our studies have demonstrated that β2-m is involved extensively in the functional regulation of growth, survival, apoptosis, and even metastasis of cancer cells. We found that β2-m is a soluble growth factor and a pleiotropic signaling molecule which interacts with its receptor, hemochromatosis protein, to modulate epithelial-to-mesenchymal transition (EMT) through iron-responsive pathways. Specific antibodies against β2-m have remarkable tumoricidal activity in cancer, through β2-m action on iron flux, alterations of intracellular reactive oxygen species, DNA damage and repair enzyme activities, β-catenin activation and cadherin switching, and tumor responsiveness to hypoxia. These novel functions of β2-m and β2-m signaling may be common to several solid tumors including human lung, breast, renal, and prostate cancers. Our experimental results could lead to the development of a novel class of antibody-based pharmaceutical agents for cancer growth control. In this review, we briefly summarize the recent data regarding β2-m as a promising new cancer therapeutic target and discuss antagonizing this therapeutic target with antibody therapy for the treatment of localized and disseminated cancers.

Pathogenesis of idiosyncratic drug-induced liver injury and clinical perspectives

Idiosyncratic drug-induced liver injury (DILI) is a rare disease that develops independently of drug dose, route, or duration of administration. Furthermore, idiosyncratic DILI is not a single disease entity but rather a spectrum of rare diseases with varying clinical, histological, and laboratory features. The pathogenesis of DILI is not fully understood. Standardization of the DILI nomenclature and methods to assess causality, along with the information provided by the LiverTox Web site, will harmonize and accelerate research on DILI. Studies of new serum biomarkers such as glutamate dehydrogenase, high mobility group box protein 1, and microRNA-122 could provide information for use in diagnosis and prognosis and provide important insights into the mechanisms of the pathogenesis of DILI. Single nucleotide polymorphisms in the HLA region have been associated with idiosyncratic hepatotoxicity attributed to flucloxacillin, ximelagatran, lapatinib, and amoxicillin-clavulanate. However, genome-wide association studies of pooled cases have not associated any genetic factors with idiosyncratic DILI. Whole genome and whole exome sequencing analyses are under way to study cases of DILI attributed to a single medication. Serum proteomic, transcriptome, and metabolome as well as intestinal microbiome analyses will increase our understanding of the mechanisms of this disorder. Further improvements to in vitro and in vivo test systems should advance our understanding of the causes, risk factors, and mechanisms of idiosyncratic DILI.

Peptide-binding motifs and characteristics for HLA -B*13:01 molecule

Trichloroethylene (TCE)-induced hypersensitivity dermatitis is one of the critical occupational diseases among workers in China. Our previous studies have identified a strong linkage between the disease and the HLA-B*13:01 allele. In this study, we searched for peptides bound to the HLA-B*13:01 molecule; 57 HLA-B*13:01-bound peptides in total were identified and 54 peptides were used to calculate frequency of amino acid residues to obtain binding motifs of HLA-B*13:01 molecule. The results showed P2, P3, and P9 were the primary binding anchor positions with the dominant anchor motifs of L, Q at P2, L at P9, D at P3. HLA-B*13:01-bound peptides were identified for the first time in our research, the results of which could contribute to the human leukocyte antigen (HLA)-binding peptides database.

HLA-B50 polymorphism in the Saudi population

The HLA-B50 serologic family is very frequent in people of Arabic origin. In Saudi Arabia, HLA-B50 is the most frequent HLA-B allele. The aim of this study was to investigate the distribution of HLA-B50 alleles in healthy Saudi individuals. A total of 162 healthy Saudi individuals were selected based on low-resolution HLA typing. DNA samples were typed by sequence-based typing method for exons 2, 3 and 4 of the HLA-B locus (Genome Diagnostics B.V.). The HLA-B*50 alleles were analysed using SBT engine software. HLA-B*50:01:01 was found in 161 of 162 individuals (99.4%), while HLA-B*50:09 was found in one individual (0.6%). HLA-B*50:01:01 is the most common HLA-B50 allele in Saudi Arabia.

The functional significance of cattle major histocompatibility complex class I genetic diversity

Current concerns about food security highlight the importance of maintaining productive and disease-resistant livestock populations. Major histocompatibility complex (MHC) class I genes have a central role in immunity. A high level of diversity in these genes allows populations to survive despite exposure to rapidly evolving pathogens. This review aims to describe the key features of MHC class I genetic diversity in cattle and to discuss their role in disease resistance. Discussion centers on data derived from the cattle genome sequence and studies addressing MHC class I gene expression and function. The impact of intensive selection on MHC diversity is also considered. A high level of complexity in MHC class I genes and functionally related gene families is revealed. This highlights the need for increased efforts to determine key genetic components that govern cattle immune responses to disease, which is increasingly important in the face of changing human and environmental demands.

Unique peptide-binding motif for Mamu-B*037:01: an MHC class I allele common to Indian and Chinese rhesus macaques

Indian and Chinese rhesus macaques are often used in biomedical research. Genetic analyses of the major histocompatibility class I region have revealed that these macaques display a substantial level of polymorphism at Mamu-A and Mamu-B loci, which have been subject to duplication. Only a few Mamu class I allotypes are characterised for their peptide-binding motifs, although more information of this nature would contribute to a better interpretation of T cell-mediated immune responses. Here, we present the results of the characterisation of the functional properties of Mamu-B*037:01, an allotype commonly encountered in rhesus macaques of Indian and Chinese origin. Mamu-B*037:01 is seen to have a strong preference for acidic amino acids at the third residue, and for arginine, lysine, and tyrosine at the carboxyl terminus. This peptide-binding motif is not described in the human population.

The repertoire of MHC class I genes in the common marmoset: evidence for functional plasticity

In humans, the classical antigen presentation function of major histocompatibility complex (MHC) class I molecules is controlled by the human leukocyte antigen HLA -A, HLA-B and HLA-C loci. A similar observation has been made for great apes and Old World monkey species. In contrast, a New World monkey species such as the cotton-top tamarin (Saguinus oedipus) appears to employ the G locus for its classical antigen presentation function. At present, little is known about the classical MHC class I repertoire of the common marmoset (Callithrix jacchus), another New World monkey that is widely used in biomedical research. In the present population study, no evidence has been found for abundant transcription of classical I class genes. However, in each common marmoset, four to seven different G-like alleles were detected, suggesting that the ancestral locus has been subject to expansion. Segregation studies provided evidence for at least two G-like genes present per haplotype, which are transcribed by a variety of cell types. The alleles of these Caja-G genes cluster in separate lineages, suggesting that the loci diversified considerably after duplication. Phylogenetic analyses of the introns confirm that the Caja-G loci cluster in the vicinity of HLA-G, indicating that both genes shared an ancestor. In contrast to HLA-G, Caja-G shows considerable polymorphism at the peptide-binding sites. This observation, together with the lack of detectable transcripts of A and B-like genes, indicates that Caja-G genes have taken over the function of classical class I genes. These data highlight the extreme plasticity of the MHC class I gene system.

Genotype List String: a grammar for describing HLA and KIR genotyping results in a text string

Knowledge of an individual's human leukocyte antigen (HLA) genotype is essential for modern medical genetics, and is crucial for hematopoietic stem cell and solid-organ transplantation. However, the high levels of polymorphism known for the HLA genes make it difficult to generate an HLA genotype that unambiguously identifies the alleles that are present at a given HLA locus in an individual. For the last 20 years, the histocompatibility and immunogenetics community has recorded this HLA genotyping ambiguity using allele codes developed by the National Marrow Donor Program (NMDP). While these allele codes may have been effective for recording an HLA genotyping result when initially developed, their use today results in increased ambiguity in an HLA genotype, and they are no longer suitable in the era of rapid allele discovery and ultra-high allele polymorphism. Here, we present a text string format capable of fully representing HLA genotyping results. This Genotype List (GL) String format is an extension of a proposed standard for reporting killer-cell immunoglobulin-like receptor (KIR) genotype data that can be applied to any genetic data that use a standard nomenclature for identifying variants. The GL String format uses a hierarchical set of operators to describe the relationships between alleles, lists of possible alleles, phased alleles, genotypes, lists of possible genotypes, and multilocus unphased genotypes, without losing typing information or increasing typing ambiguity. When used in concert with appropriate tools to create, exchange, and parse these strings, we anticipate that GL Strings will replace NMDP allele codes for reporting HLA genotypes.

Major histocompatibility complex genomics and human disease

Over several decades, various forms of genomic analysis of the human major histocompatibility complex (MHC) have been extremely successful in picking up many disease associations. This is to be expected, as the MHC region is one of the most gene-dense and polymorphic stretches of human DNA. It also encodes proteins critical to immunity, including several controlling antigen processing and presentation. Single-nucleotide polymorphism genotyping and human leukocyte antigen (HLA) imputation now permit the screening of large sample sets, a technique further facilitated by high-throughput sequencing. These methods promise to yield more precise contributions of MHC variants to disease. However, interpretation of MHC-disease associations in terms of the functions of variants has been problematic. Most studies confirm the paramount importance of class I and class II molecules, which are key to resistance to infection. Infection is likely driving the extreme variation of these genes across the human population, but this has been difficult to demonstrate. In contrast, many associations with autoimmune conditions have been shown to be specific to certain class I and class II alleles. Interestingly, conditions other than infections and autoimmunity are also associated with the MHC, including some cancers and neuropathies. These associations could be indirect, owing, for example, to the infectious history of a particular individual and selective pressures operating at the population level.

Extensive variation at MHC DRB in the New Zealand sea lion (Phocarctos hookeri) provides evidence for balancing selection

Marine mammals are often reported to possess reduced variation of major histocompatibility complex (MHC) genes compared with their terrestrial counterparts. We evaluated diversity at two MHC class II B genes, DQB and DRB, in the New Zealand sea lion (Phocarctos hookeri, NZSL) a species that has suffered high mortality owing to bacterial epizootics, using Sanger sequencing and haplotype reconstruction, together with next-generation sequencing. Despite this species' prolonged history of small population size and highly restricted distribution, we demonstrate extensive diversity at MHC DRB with 26 alleles, whereas MHC DQB is dimorphic. We identify four DRB codons, predicted to be involved in antigen binding, that are evolving under adaptive evolution. Our data suggest diversity at DRB may be maintained by balancing selection, consistent with the role of this locus as an antigen-binding region and the species' recent history of mass mortality during a series of bacterial epizootics. Phylogenetic analyses of DQB and DRB sequences from pinnipeds and other carnivores revealed significant allelic diversity, but little phylogenetic depth or structure among pinniped alleles; thus, we could neither confirm nor refute the possibility of trans-species polymorphism in this group. The phylogenetic pattern observed however, suggests some significant evolutionary constraint on these loci in the recent past, with the pattern consistent with that expected following an epizootic event. These data may help further elucidate some of the genetic factors underlying the unusually high susceptibility to bacterial infection of the threatened NZSL, and help us to better understand the extent and pattern of MHC diversity in pinnipeds.

Immunogenetics of HIV disease

Host genetic factors are a major contributing factor to the inter-individual variation observed in response to human immunodeficiency virus (HIV) infection and are linked to resistance to HIV infection among exposed individuals, as well as rate of disease progression and the likelihood of viral transmission. Of the genetic variants that have been shown to affect the natural history of HIV infection, the human leukocyte antigen (HLA) class I genes exhibit the strongest and most consistent association, underscoring a central role for CD8(+) T cells in resistance to the virus. HLA proteins play important roles in T-cell-mediated adaptive immunity by presenting immunodominant HIV epitopes to cytotoxic T lymphocytes (CTLs) and CD4(+) T cells. Genetic and functional data also indicate a function for HLA in natural killer cell-mediated innate immunity against HIV by interacting with killer cell immunoglobulin-like receptors (KIR). We review the HLA and KIR associations with HIV disease and discuss the mechanisms underlying these associations.

Duplication and population dynamics shape historic patterns of selection and genetic variation at the major histocompatibility complex in rodents

Genetic variation at the major histocompatibility complex (MHC) is vitally important for wildlife populations to respond to pathogen threats. As natural populations can fluctuate greatly in size, a key issue concerns how population cycles and bottlenecks that could reduce genetic diversity will influence MHC genes. Using 454 sequencing, we characterized genetic diversity at the DRB Class II locus in montane voles (Microtus montanus), a North American rodent that regularly undergoes high-amplitude fluctuations in population size. We tested for evidence of historic balancing selection, recombination, and gene duplication to identify mechanisms maintaining allelic diversity. Counter to our expectations, we found strong evidence of purifying selection acting on the DRB locus in montane voles. We speculate that the interplay between population fluctuations and gene duplication might be responsible for the weak evidence of historic balancing selection and strong evidence of purifying selection detected. To further explore this idea, we conducted a phylogenetically controlled comparative analysis across 16 rodent species with varying demographic histories and MHC duplication events (based on the maximum number of alleles detected per individual). On the basis of phylogenetic generalized linear model-averaging, we found evidence that the estimated number of duplicated loci was positively related to allelic diversity and, surprisingly, to the strength of purifying selection at the DRB locus. Our analyses also revealed that species that had undergone population bottlenecks had lower allelic richness than stable species. This study highlights the need to consider demographic history and genetic structure alongside patterns of natural selection to understand resulting patterns of genetic variation at the MHC.

Hypervariable antigen genes in malaria have ancient roots

Background

The var genes of the human malaria parasite Plasmodium falciparum are highly polymorphic loci coding for the erythrocyte membrane proteins 1 (PfEMP1), which are responsible for the cytoaherence of P. falciparum infected red blood cells to the human vasculature. Cytoadhesion, coupled with differential expression of var genes, contributes to virulence and allows the parasite to establish chronic infections by evading detection from the host’s immune system. Although studying genetic diversity is a major focus of recent work on the var genes, little is known about the gene family's origin and evolutionary history.

Results

Using a novel hidden Markov model-based approach and var sequences assembled from additional isolates and species, we are able to reveal elements of both the early evolution of the var genes as well as recent diversifying events. We compare sequences of the var gene DBLα domains from divergent isolates of P. falciparum (3D7 and HB3), and a closely-related species, Plasmodium reichenowi. We find that the gene family is equally large in P. reichenowi and P. falciparum -- with a minimum of 51 var genes in the P. reichenowi genome (compared to 61 in 3D7 and a minimum of 48 in HB3). In addition, we are able to define large, continuous blocks of homologous sequence among P. falciparum and P. reichenowi var gene DBLα domains. These results reveal that the contemporary structure of the var gene family was present before the divergence of P. falciparum and P. reichenowi, estimated to be between 2.5 to 6 million years ago. We also reveal that recombination has played an important and traceable role in both the establishment, and the maintenance, of diversity in the sequences.

Conclusions

Despite the remarkable diversity and rapid evolution found in these loci within and among P. falciparum populations, the basic structure of these domains and the gene family is surprisingly old and stable. Revealing a common structure as well as conserved sequence among two species also has implications for developing new primate-parasite models for studying the pathology and immunology of falciparum malaria, and for studying the population genetics of var genes and associated virulence phenotypes.

The HIV-1 pandemic: does the selective sweep in chimpanzees mirror humankind's future?

An HIV-1 infection progresses in most human individuals sooner or later into AIDS, a devastating disease that kills more than a million people worldwide on an annual basis. Nonetheless, certain HIV-1-infected persons appear to act as long-term non-progressors, and elite control is associated with the presence of particular MHC class I allotypes such as HLA-B*27 or -B*57. The HIV-1 pandemic in humans arose from the cross-species transmission of SIV_cpz originating from chimpanzees. Chimpanzees, however, appear to be relatively resistant to developing AIDS after HIV-1/SIV_cpz infection. Mounting evidence illustrates that, in the distant past, chimpanzees experienced a selective sweep resulting in a severe reduction of their MHC class I repertoire. This was most likely caused by an HIV-1/SIV-like retrovirus, suggesting that chimpanzees may have experienced long-lasting host-virus relationships with SIV-like viruses. Hence, if natural selection is allowed to follow its course, prospects for the human population may look grim, thus underscoring the desperate need for an effective vaccine.

HLA genetic profile of Mapuche (Araucanian) Amerindians from Chile

Amerindian Mapuche (Araucanians) are now living in Chile and Argentina at both sides of Andean Mountains. They are anthropologically and genetically different from southernmost South America Patagonian Amerindians. Most of the HLA alleles found in our Mapuche sample are frequent or very frequent in North and South America Amerindians: (1) Class I: A*02:01, A*03:01, A*68:01, B*39:09, B*51:01, (2) Class II: DRB1*03:01, DRB1*04:03, DRB1*07:01, DRB1*08:02, DRB1*14:02, DRB1*16:02. One of the nine most frequent extended haplotypes seems to be from European origin, suggesting the existence of a degree of admixture with Europeans in our Mapuche sample. It has been calculated of about 11 % admixture. Three of the extended haplotypes are also found in other Amerindians and five of them are newly found in Mapuche Amerindians: A*68:01-B*39:09-DRB1*08:02-DQB1*04:02; A*68:01-B*51:01-DRB1*04:03-DQB1*03:02; A*29:01-B*08:01-DRB1*03:01-DQB1*02:01; A*02:01-B*15:01-DRB1*04:03-DQB1*03:02; A*33:01-B*14:02-DRB1*07:01-DQB1*03:03. The medical importance of calculating HLA profile is discussed on the diagnostic (HLA and disease) and therapeutical bases of HLA pharmacogenomics and on the construction of a virtual transplantation HLA list profile. Also, anthropological conclusions are drawn.

Evolution of MHC class I genes in the endangered loggerhead sea turtle (Caretta caretta) revealed by 454 amplicon sequencing

Background

In evolutionary and conservation biology, parasitism is often highlighted as a major selective pressure. To fight against parasites and pathogens, genetic diversity of the immune genes of the major histocompatibility complex (MHC) are particularly important. However, the extensive degree of polymorphism observed in these genes makes it difficult to conduct thorough population screenings.

Methods

We utilized a genotyping protocol that uses 454 amplicon sequencing to characterize the MHC class I in the endangered loggerhead sea turtle (Caretta caretta) and to investigate their evolution at multiple relevant levels of organization.

Results

MHC class I genes revealed signatures of trans-species polymorphism across several reptile species. In the studied loggerhead turtle individuals, it results in the maintenance of two ancient allelic lineages. We also found that individuals carrying an intermediate number of MHC class I alleles are larger than those with either a low or high number of alleles.

Conclusions

Multiple modes of evolution seem to maintain MHC diversity in the loggerhead turtles, with relatively high polymorphism for an endangered species.

A glow of HLA typing in organ transplantation

The transplant of organs and tissues is one of the greatest curative achievements of this century. In organ transplantation, the adaptive immunity is considered the main response exerted to the transplanted tissue, since the main goal of the immune response is the MHC (major histocompatibility complex) molecules expressed on the surface of donor cells. Cell surface molecules that induce an antigenic stimulus cause the rejection immune response to grafted tissue or organ. A wide variety of transplantation antigens have been described, including the major histocompatibility molecules, minor histocompatibility antigens, ABO blood group antigens and endothelial cell antigens. The sensitization to MHC antigens may be caused by transfusions, pregnancy, or failed previous grafts leading to development of anti-human leukocyte antigen (HLA) antibodies that are important factor responsible for graft rejection in solid organ transplantation and play a role in post-transfusion complication Anti-HLA Abs may be present in healthy individuals. Methods for HLA typing are described, including serological methods, molecular techniques of sequence-specific priming (SSP), sequence-specific oligonucleotide probing (SSOP), Sequence based typing (SBT) and reference strand-based conformation analysis (RSCA) method. Problems with organ transplantation are reservoir of organs and immune suppressive treatments that used to decrease rate of rejection with less side effect and complications.

Clinical Pharmacogenetics Implementation Consortium guidelines for human leukocyte antigen-B genotype and allopurinol dosing

Allopurinol is the most commonly used drug for the treatment of hyperuricemia and gout. However, allopurinol is also one of the most common causes of severe cutaneous adverse reactions (SCARs), which include drug hypersensitivity syndrome, Stevens–Johnson syndrome, and toxic epidermal necrolysis. A variant allele of the human leukocyte antigen (HLA)-B, HLA-B*58:01, associates strongly with allopurinolinduced SCAR. We have summarized the evidence from the published literature and developed peer-reviewed guidelines for allopurinol use based on HLA-B genotype.

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

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

Genetic variants contribute to gene expression variability in humans

Expression quantitative trait loci (eQTL) studies have established convincing relationships between genetic variants and gene expression. Most of these studies focused on the mean of gene expression level, but not the variance of gene expression level (i.e., gene expression variability). In the present study, we systematically explore genome-wide association between genetic variants and gene expression variability in humans. We adapt the double generalized linear model (dglm) to simultaneously fit the means and the variances of gene expression among the three possible genotypes of a biallelic SNP. The genomic loci showing significant association between the variances of gene expression and the genotypes are termed expression variability QTL (evQTL). Using a data set of gene expression in lymphoblastoid cell lines (LCLs) derived from 210 HapMap individuals, we identify cis-acting evQTL involving 218 distinct genes, among which 8 genes, ADCY1, CTNNA2, DAAM2, FERMT2, IL6, PLOD2, SNX7, and TNFRSF11B, are cross-validated using an extra expression data set of the same LCLs. We also identify ∼300 trans-acting evQTL between >13,000 common SNPs and 500 randomly selected representative genes. We employ two distinct scenarios, emphasizing single-SNP and multiple-SNP effects on expression variability, to explain the formation of evQTL. We argue that detecting evQTL may represent a novel method for effectively screening for genetic interactions, especially when the multiple-SNP influence on expression variability is implied. The implication of our results for revealing genetic mechanisms of gene expression variability is discussed.

HLA-A, -B and -DRB1 allele frequencies in Cyrenaica population (Libya) and genetic relationships with other populations

The frequencies of HLA-A, HLA-B and HLA-DRB1 alleles in 118 unrelated Libyans from Benghazi (Cyrenaica) were analysed using high resolution typing and compared with other populations. Their relatedness has been tested by correspondence analyses and principal component analysis. The most frequent HLA-A alleles were A(∗)02:01:01:01 (15.7%), A(∗)01:01:01:01 (11.4%) and A(∗)03:01:01:01 (9.3%). For the HLA-B locus, the commonest allele was HLA-B(∗)50:01:01 (14.4%) followed by B(∗)51:01:01 (9.8%) and B(∗)08:01:01 (6.4%). For the HLA-DRB1 locus, the commonest was HLA-DRB1(∗)07:01:01:01 (16.9%) followed by DRB1(∗)03:01:01:01 (13.6%) and DRB1(∗)13:02:01 (9.3%). The most frequent two-locus haplotypes were HLA-A(∗)02:01:01:01-B(∗)07:02:01 (3.0%) and HLA-B(∗)50:01:01-DRB1(∗)07:01:01:01 (9.6%), and three-locus haplotypes were HLA-A(∗)02:01:01:01-B(∗)50:01:01-DRB1(∗)07:01:01:01 (4.2%) and HLA-A(∗)11:01:01-B(∗)52:01:01:01-DRB1(∗)15:02:01 (2.5%). This study is the first on the HLA status of a Libyan population. The results, when compared to similar HLA data obtained previously from African and Mediterranean populations, indicate genetic influences from several ethnic groups. Moreover, the differences in the HLA allele frequencies between the Libyan population and others reveals that significant admixture has occurred between the original Berber inhabitants and neighbouring and more distant populations, even though a strong genetic Berber substratum remains. These data will be of value to future anthropological and disease association studies involving the Libyan population.

Evidence for selection on a chordate histocompatibility locus

Allorecognition is the ability of an organism to differentiate self or close relatives from unrelated individuals. The best known applications of allorecognition are the prevention of inbreeding in hermaphroditic species (e.g., the self-incompatibility [SI] systems in plants), the vertebrate immune response to foreign antigens mediated by MHC loci, and somatic fusion, where two genetically independent individuals physically join to become a chimera. In the few model systems where the loci governing allorecognition outcomes have been identified, the corresponding proteins have exhibited exceptional polymorphism. But information about the evolution of this polymorphism outside MHC is limited. We address this subject in the ascidian Botryllus schlosseri, where allorecognition outcomes are determined by a single locus, called FuHC (Fusion/HistoCompatibility). Molecular variation in FuHC is distributed almost entirely within populations, with very little evidence for differentiation among different populations. Mutation plays a larger role than recombination in the creation of FuHC polymorphism. A selection statistic, neutrality tests, and distribution of variation within and among different populations all provide evidence for selection acting on FuHC, but are not in agreement as to whether the selection is balancing or directional.

Microfluidic devices modulate tumor cell line susceptibility to NK cell recognition

This study aims to adoptively reduce the major histocompatibility complex class I (MHC-I) molecule surface expression of cancer cells by exposure to microfluid shear stress and a monoclonal antibody. A microfluidic system is developed and tumor cells are injected at different flow rates. The bottom surface of the microfluidic system is biofunctionalized with antibodies (W6/32) specific for the MHC-I molecules with a simple method based on microfluidic protocols. The antibodies promote binding between the bottom surface and the MHC-I molecules on the tumor cell membrane. The cells are injected at an optimized flow rate, then roll on the bottom surface and are subjected to shear stress. The stress is localized and enhanced on the part of the membrane where MHC-I proteins are expressed, since they stick to the antibodies of the system. The localized stress allows a stripping effect and consequent reduction of the MHC-I expression. It is shown that it is possible to specifically treat and recover eukaryotic cells without damaging the biological samples. MHC-I molecule expression on treated and control cell surfaces is measured on tumor and healthy cells. After the cell rolling treatment a clear reduction of MHC-I levels on the tumor cell membrane is observed, whereas no changes are observed on healthy cells (monocytes). The MHC-I reduction is investigated and the possibility that the developed system could induce a loss of these molecules from the tumor cell surface is addressed. The percentage of living tumor cells (viability) that remain after the treatment is measured. The changes induced by the microfluidic system are analyzed by fluorescence-activated cell sorting and confocal microscopy. Cytotoxicity tests show a relevant increased susceptibility of natural killer (NK) cells on microchip-treated tumor cells.

Natural selection on marine carnivores elaborated a diverse family of classical MHC class I genes exhibiting haplotypic gene content variation and allelic polymorphism

Pinnipeds, marine carnivores, diverged from terrestrial carnivores ~45 million years ago, before their adaptation to marine environments. This lifestyle change exposed pinnipeds to different microbiota and pathogens, with probable impact on their MHC class I genes. Investigating this question, genomic sequences were determined for 71 MHC class I variants: 27 from harbor seal and 44 from gray seal. These variants form three MHC class I gene lineages, one comprising a pseudogene. The second, a candidate nonclassical MHC class I gene, comprises a nonpolymorphic transcribed gene related to dog DLA-79 and giant panda Aime-1906. The third is the diversity lineage, which includes 62 of the 71 seal MHC class I variants. All are transcribed, and they minimally represent six harbor and 12 gray seal MHC class I genes. Besides species-specific differences in gene number, seal MHC class I haplotypes exhibit gene content variation and allelic polymorphism. Patterns of sequence variation, and of positions for positively selected sites, indicate the diversity lineage genes are the seals' classical MHC class I genes. Evidence that expansion of diversity lineage genes began before gray and harbor seals diverged is the presence in both species of two distinctive sublineages of diversity lineage genes. Pointing to further expansion following the divergence are the presence of species-specific genes and greater MHC class I diversity in gray seals than harbor seals. The elaboration of a complex variable family of classical MHC class I genes in pinnipeds contrasts with the single, highly polymorphic classical MHC class I gene of dog and giant panda, terrestrial carnivores.

Evolution of MHC class I genes in the European badger (Meles meles)

The major histocompatibility complex (MHC) plays a central role in the adaptive immune system and provides a good model with which to understand the evolutionary processes underlying functional genes. Trans-species polymorphism and orthology are both commonly found in MHC genes; however, mammalian MHC class I genes tend to cluster by species. Concerted evolution has the potential to homogenize different loci, whereas birth-and-death evolution can lead to the loss of orthologs; both processes result in monophyletic groups within species. Studies investigating the evolution of MHC class I genes have been biased toward a few particular taxa and model species. We present the first study of MHC class I genes in a species from the superfamily Musteloidea. The European badger (Meles meles) exhibits moderate variation in MHC class I sequences when compared to other carnivores. We identified seven putatively functional sequences and nine pseudogenes from genomic (gDNA) and complementary (cDNA) DNA, signifying at least two functional class I loci. We found evidence for separate evolutionary histories of the α1 and α2/α3 domains. In the α1 domain, several sequences from different species were more closely related to each other than to sequences from the same species, resembling orthology or trans-species polymorphism. Balancing selection and probable recombination maintain genetic diversity in the α1 domain, evidenced by the detection of positive selection and a recombination event. By comparison, two recombination breakpoints indicate that the α2/α3 domains have most likely undergone concerted evolution, where recombination has homogenized the α2/α3 domains between genes, leading to species-specific clusters of sequences. Our findings highlight the importance of analyzing MHC domains separately.

Why must T cells be cross-reactive?

Clonal selection theory proposed that individual T cells are specific for a single peptide-MHC antigen. However, the repertoire of αβ T cell receptors (TCRs) is dwarfed by the vast array of potential foreign peptide-MHC complexes, and a comprehensive system requires each T cell to recognize numerous peptides and thus be cross-reactive. This compromise on specificity has profound implications because the chance of any natural peptide-MHC ligand being an optimal fit for its cognate TCR is small, as there will almost always be more-potent agonists. Furthermore, any TCR raised against a specific peptide-MHC complex in vivo can only be the best available solution from the naive T cell pool and is unlikely to be the best possible solution from the substantially greater number of TCRs that could theoretically be produced. This 'systems view' of TCR recognition provides a plausible cause for autoimmune disease and substantial scope for multiple therapeutic interventions.

The complement-mediated prozone effect in the Luminex single-antigen bead assay and its impact on HLA antibody determination in patient sera

The Luminex xMAP system has become an important tool for HLA antibody screening and identification in sera of transplant patients. Recently, the Luminex single antigen bead assay was shown to be prone to an artefact, the so called prozone phenomenon: Sera with high titer HLA antibodies gave negative results when tested neat, but reacted strongly positive after 1:10 dilution. We also observed such a phenomenon and found that it was most likely caused by the complement component 1 (C1) by competitively displacing the detection antibodies. In this article we review the complement-mediated prozone effect and other mechanisms of interference with solid phase assays, and we discuss possible consequences for HLA antibody testing with the Luminex SAB assay.

Immune self-reactivity triggered by drug-modified HLA-peptide repertoire

Human leukocyte antigens (HLAs) are highly polymorphic proteins that initiate immunity by presenting pathogen-derived peptides to T cells. HLA polymorphisms mostly map to the antigen-binding cleft, thereby diversifying the repertoire of self-derived and pathogen-derived peptide antigens selected by different HLA allotypes. A growing number of immunologically based drug reactions, including abacavir hypersensitivity syndrome (AHS) and carbamazepine-induced Stevens-Johnson syndrome (SJS), are associated with specific HLA alleles. However, little is known about the underlying mechanisms of these associations, including AHS, a prototypical HLA-associated drug reaction occurring exclusively in individuals with the common histocompatibility allele HLA-B*57:01, and with a relative risk of more than 1,000 (refs 6, 7). We show that unmodified abacavir binds non-covalently to HLA-B*57:01, lying across the bottom of the antigen-binding cleft and reaching into the F-pocket, where a carboxy-terminal tryptophan typically anchors peptides bound to HLA-B*57:01. Abacavir binds with exquisite specificity to HLA-B*57:01, changing the shape and chemistry of the antigen-binding cleft, thereby altering the repertoire of endogenous peptides that can bind HLA-B*57:01. In this way, abacavir guides the selection of new endogenous peptides, inducing a marked alteration in 'immunological self'. The resultant peptide-centric 'altered self' activates abacavir-specific T-cells, thereby driving polyclonal CD8 T-cell activation and a systemic reaction manifesting as AHS. We also show that carbamazepine, a widely used anti-epileptic drug associated with hypersensitivity reactions in HLA-B*15:02 individuals, binds to this allotype, producing alterations in the repertoire of presented self peptides. Our findings simultaneously highlight the importance of HLA polymorphism in the evolution of pharmacogenomics and provide a general mechanism for some of the growing number of HLA-linked hypersensitivities that involve small-molecule drugs.

Ancestral polymorphism at the major histocompatibility complex (MHCIIß) in the Nesospiza bunting species complex and its sister species (Rowettia goughensis)

BACKGROUND

The major histocompatibility complex (MHC) is an important component of the vertebrate immune system and is frequently used to characterise adaptive variation in wild populations due to its co-evolution with pathogens. Passerine birds have an exceptionally diverse MHC with multiple gene copies and large numbers of alleles compared to other avian taxa. The Nesospiza bunting species complex (two species on Nightingale Island; one species with three sub-species on Inaccessible Island) represents a rapid adaptive radiation at a small, isolated archipelago, and is thus an excellent model for the study of adaptation and speciation. In this first study of MHC in Nesospiza buntings, we aim to characterize MHCIIß variation, determine the strength of selection acting at this gene region and assess the level of shared polymorphism between the Nesospiza species complex and its putative sister taxon, Rowettia goughensis, from Gough Island.

RESULTS

In total, 23 unique alleles were found in 14 Nesospiza and 2 R. goughensis individuals encoding at least four presumably functional loci and two pseudogenes. There was no evidence of ongoing selection on the peptide binding region (PBR). Of the 23 alleles, 15 were found on both the islands inhabited by Nesospiza species, and seven in both Nesospiza and Rowettia; indications of shared, ancestral polymorphism. A gene tree of Nesospiza MHCIIß alleles with several other passerine birds shows three highly supported Nesospiza-specific groups. All R. goughensis alleles were shared with Nesospiza, and these alleles were found in all three Nesospiza sequence groups in the gene tree, suggesting that most of the observed variation predates their phylogenetic split.

CONCLUSIONS

Lack of evidence of selection on the PBR, together with shared polymorphism across the gene tree, suggests that population variation of MHCIIß among Nesospiza and Rowettia is due to ancestral polymorphism rather than local selective forces. Weak or no selection pressure could be attributed to low parasite load at these isolated Atlantic islands. The deep divergence between the highly supported Nesospiza-specific sequence Groups 2 and 3, and the clustering of Group 3 close to the distantly related passerines, provide strong support for preserved ancestral polymorphism, and present evidence of one of the rare cases of extensive ancestral polymorphism in birds.

Multilocus definition of MHC haplotypes in pedigreed cynomolgus macaques (Macaca fascicularis)

Cynomolgus macaques (Macaca fascicularis) are used widely in biomedical research, and the genetics of their MHC (Mhc-Mafa) has become the focus of considerable attention in recent years. The cohort of Indonesian pedigreed macaques that we present here was typed for Mafa-A, -B, and -DR, by sequencing, as described in earlier studies. Additionally, the DRB region of these animals was characterised by microsatellite analyses. In this study, full-length sequencing of Mafa-DPA/B and -DQA/B in these animals was performed. A total of 75 different alleles were observed; 22 of which have not previously been reported, plus 18 extended exon 2 alleles that were already known. Furthermore, two microsatellites, D6S2854 and D6S2859, were used to characterise the complex Mafa-A region. Sequencing and segregation analyses revealed that the length patterns of these microsatellites are unique for each Mafa-A haplotype. In this work, we present a pedigreed colony of approximately 120 cynomolgus macaques; all of which are typed for the most significant polymorphic MHC class I and class II markers. Offspring of these pedigreed animals are easily characterised for their MHC by microsatellite analyses on the Mafa-A and -DRB regions, which makes the cumbersome sequencing analyses redundant.

The impact of mutation and gene conversion on the local diversification of antigen genes in African trypanosomes

Patterns of genetic diversity in parasite antigen gene families hold important information about their potential to generate antigenic variation within and between hosts. The evolution of such gene families is typically driven by gene duplication, followed by point mutation and gene conversion. There is great interest in estimating the rates of these processes from molecular sequences for understanding the evolution of the pathogen and its significance for infection processes. In this study, a series of models are constructed to investigate hypotheses about the nucleotide diversity patterns between closely related gene sequences from the antigen gene archive of the African trypanosome, the protozoan parasite causative of human sleeping sickness in Equatorial Africa. We use a hidden Markov model approach to identify two scales of diversification: clustering of sequence mismatches, a putative indicator of gene conversion events with other lower-identity donor genes in the archive, and at a sparser scale, isolated mismatches, likely arising from independent point mutations. In addition to quantifying the respective probabilities of occurrence of these two processes, our approach yields estimates for the gene conversion tract length distribution and the average diversity contributed locally by conversion events. Model fitting is conducted using a Bayesian framework. We find that diversifying gene conversion events with lower-identity partners occur at least five times less frequently than point mutations on variant surface glycoprotein (VSG) pairs, and the average imported conversion tract is between 14 and 25 nucleotides long. However, because of the high diversity introduced by gene conversion, the two processes have almost equal impact on the per-nucleotide rate of sequence diversification between VSG subfamily members. We are able to disentangle the most likely locations of point mutations and conversions on each aligned gene pair.

Immunogenetic surveillance of HIV/AIDS

Evolutionary pressure by viruses is most likely responsible for the extraordinary allelic polymorphism of genes encoding class I human leukocyte antigens (HLA) and killer immunoglobulin-like receptors (KIR). Such genetic diversity has functional implications for the immune response to viruses and generates population-based variations in HLA class I allele frequencies and KIR gene profiles. The HIV-1 virus has relatively recently established itself as a major human pathogen, rapidly diversifying into a variety of phylogenetic subtypes or clades (A-G) and recombinants in different populations. HIV-1 clade C is the most common subtype in circulation accounting for 48% of all infections, followed by HIV-1 clades A and B which are responsible for 13% and 11% of infections in the current pandemic, respectively. Candidate gene studies of large cohorts of predominantly HIV-1 clade B but also clades C and A infected patients, have consistently shown significant associations between certain HLA class I alleles namely HLA-B*57, B*58, B*27, B*51 and relatively low viraemia. However, there is evidence that other associations between HLA-B*15, B*18 or B*53 and levels of HIV-1 viraemia are clade-specific. Recent genome-wide association studies of HIV-1 clade B exposed cohorts have confirmed that HLA-B, which is the most polymorphic locus in the human genome, is the major genetic locus contributing to immune control of viraemia. Moreover, the presence of natural killer cell receptors encoded by KIR-3DL1 and 3DS1 genes together with certain HLA class I alleles carrying the KIR target motif Bw4Ile80, provides an enhanced ability to control HIV-1 viraemia in some individuals. It is likely that rapid co-evolution of HIV-1 immune escape variants together with an adjustment of human immune response gene profiles has occurred in some exposed populations. Taken together, immunogenetic surveillance of HIV-1 exposed cohorts has revealed important correlates of natural immunity, which could provide a rational platform for the design and testing of future vaccines aimed at controlling the current AIDS pandemic.

Genetic diversity of MHC class I loci in six non-model frogs is shaped by positive selection and gene duplication

Comparative studies of major histocompatibility complex (MHC) genes across vertebrate species can reveal the evolutionary processes that shape the structure and function of immune regulatory proteins. In this study, we characterized MHC class I sequences from six frog species representing three anuran families (Hylidae, Centrolenidae and Ranidae). Using cDNA from our focal species, we amplified a total of 79 unique sequences spanning exons 2-4 that encode the extracellular domains of the functional alpha chain protein. We compared intra- and interspecific nucleotide and amino-acid divergence, tested for recombination, and identified codon sites under selection by estimating the rate of non-synonymous to synonymous substitutions with multiple codon-based maximum likelihood methods. We determined that positive (diversifying) selection was acting on specific amino-acid sites located within the domains that bind pathogen-derived peptides. We also found significant signals of recombination across the physical distance of the genes. Finally, we determined that all the six species expressed two or three putative classical class I loci, in contrast to the single locus condition of Xenopus laevis. Our results suggest that MHC evolution in anurans is a dynamic process and that variation in numbers of loci and genetic diversity can exist among taxa. Thus, the accumulation of genetic data for more species will be useful in further characterizing the relative importance of processes such as selection, recombination and gene duplication in shaping MHC loci among amphibian lineages.