The evolution of MHC diversity by segmental duplication and transposition of retroelements

Universal signatures of transposable element compartmentalization across eukaryotic genomes

The evolutionary mechanisms that drive the emergence of genome architecture remain poorly understood but can now be assessed with unprecedented power due to the massive accumulation of genome assemblies spanning phylogenetic diversity1,2. Transposable elements (TEs) are a rich source of large-effect mutations since they directly and indirectly drive genomic structural variation and changes in gene expression3. Here, we demonstrate universal patterns of TE compartmentalization across eukaryotic genomes spanning ~1.7 billion years of evolution, in which TEs colocalize with gene families under strong predicted selective pressure for dynamic evolution and involved in specific functions. For non-pathogenic species these genes represent families involved in defense, sensory perception and environmental interaction, whereas for pathogenic species, TE-compartmentalized genes are highly enriched for pathogenic functions. Many TE-compartmentalized gene families display signatures of positive selection at the molecular level. Furthermore, TE-compartmentalized genes exhibit an excess of high-frequency alleles for polymorphic TE insertions in fruit fly populations. We postulate that these patterns reflect selection for adaptive TE insertions as well as TE-associated structural variants. This process may drive the emergence of a shared TE-compartmentalized genome architecture across diverse eukaryotic lineages.

HLA Genetics for the Human Diseases

Highly polymorphic human leukocyte antigen (HLA) molecules (alleles) expressed by different classical HLA class I and class II genes have crucial roles in the regulation of innate and adaptive immune responses, transplant rejection and in the pathogenesis of numerous infectious and autoimmune diseases. To date, over 35,000 HLA alleles have been published from the IPD-IMGT/HLA database, and specific HLA alleles and HLA haplotypes have been reported to be associated with more than 100 different diseases and phenotypes. Next generation sequencing (NGS) technology developed in recent years has provided breakthroughs in various HLA genomic/gene studies and transplant medicine. In this chapter, we review the current information on the HLA genomic structure and polymorphisms, as well as the genetic context in which numerous disease associations have been identified in this region.

Human leukocyte antigen super-locus: nexus of genomic supergenes, SNPs, indels, transcripts, and haplotypes

The human Major Histocompatibility Complex (MHC) or Human Leukocyte Antigen (HLA) super-locus is a highly polymorphic genomic region that encodes more than 140 coding genes including the transplantation and immune regulatory molecules. It receives special attention for genetic investigation because of its important role in the regulation of innate and adaptive immune responses and its strong association with numerous infectious and/or autoimmune diseases. In recent years, MHC genotyping and haplotyping using Sanger sequencing and next-generation sequencing (NGS) methods have produced many hundreds of genomic sequences of the HLA super-locus for comparative studies of the genetic architecture and diversity between the same and different haplotypes. In this special issue on 'The Current Landscape of HLA Genomics and Genetics', we provide a short review of some of the recent analytical developments used to investigate the SNP polymorphisms, structural variants (indels), transcription and haplotypes of the HLA super-locus. This review highlights the importance of using reference cell-lines, population studies, and NGS methods to improve and update our understanding of the mechanisms, architectural structures and combinations of human MHC genomic alleles (SNPs and indels) that better define and characterise haplotypes and their association with various phenotypes and diseases.

Alternative splicing and genetic variation of mhc-e: implications for rhesus cytomegalovirus-based vaccines

Rhesus cytomegalovirus (RhCMV)-based vaccination against Simian Immunodeficiency virus (SIV) elicits MHC-E-restricted CD8+ T cells that stringently control SIV infection in ~55% of vaccinated rhesus macaques (RM). However, it is unclear how accurately the RM model reflects HLA-E immunobiology in humans. Using long-read sequencing, we identified 16 Mamu-E isoforms and all Mamu-E splicing junctions were detected among HLA-E isoforms in humans. We also obtained the complete Mamu-E genomic sequences covering the full coding regions of 59 RM from a RhCMV/SIV vaccine study. The Mamu-E gene was duplicated in 32 (54%) of 59 RM. Among four groups of Mamu-E alleles: three ~5% divergent full-length allele groups (G1, G2, G2_LTR) and a fourth monomorphic group (G3) with a deletion encompassing the canonical Mamu-E exon 6, the presence of G2_LTR alleles was significantly (p = 0.02) associated with the lack of RhCMV/SIV vaccine protection. These genomic resources will facilitate additional MHC-E targeted translational research.

Class I/Class II HLA Evolutionary Divergence Ratio Is an Independent Marker Associated With Disease-Free and Overall Survival After Allogeneic Hematopoietic Stem Cell Transplantation for Acute Myeloid Leukemia

Class I Human Leukocyte Antigen (HLA) evolutionary divergence (HED) is a metric which reflects immunopeptidome diversity and has been associated with immune checkpoint inhibitor responses in solid tumors. Its impact and interest in allogeneic hematopoietic stem cell transplantation (HCT) have not yet been thoroughly studied. This study analyzed the clinical and immune impact of class I and II HED in 492 acute myeloid leukemia (AML) recipients undergoing HCT. The overall cohort was divided into a training (n=338) and a testing (n=132) set. Univariate cox screening found a positive impact of a high class I HED and a negative impact of a high class II HED on both disease-free (DFS) and overall survival (OS). These results were combined in a unique marker, class I/class II HED ratio, and assessed in the testing cohort. The final multivariate cox model confirmed the positive impact of a high versus low class I/class II HED ratio on both DFS (Hazard Ratio (HR) 0.41 [95% CI 0.2-0.83]; p=0.01) and OS (HR 0.34 [0.19-0.59]; p<0.001), independently of HLA matching and other HCT parameters. No significant association was found between the ratio and graft-versus-host disease (GvHD) nor with neutrophil and platelet recovery. A high class I HED was associated with a tendency for an increase in NK, CD8 T-cell, and B cell recovery at 12 months. These results introduce HED as an original and independent prognosis marker reflecting immunopeptidome diversity and alloreactivity after HCT.

The New Kid on the Block: HLA-C, a Key Regulator of Natural Killer Cells in Viral Immunity

The human leukocyte antigen system (HLA) is a cluster of highly polymorphic genes essential for the proper function of the immune system, and it has been associated with a wide range of diseases. HLA class I molecules present intracellular host- and pathogen-derived peptides to effector cells of the immune system, inducing immune tolerance in healthy conditions or triggering effective immune responses in pathological situations. HLA-C is the most recently evolved HLA class I molecule, only present in humans and great apes. Differentiating from its older siblings, HLA-A and HLA-B, HLA-C exhibits distinctive features in its expression and interaction partners. HLA-C serves as a natural ligand for multiple members of the killer-cell immunoglobulin-like receptor (KIR) family, which are predominately expressed by natural killer (NK) cells. NK cells are crucial for the early control of viral infections and accumulating evidence indicates that interactions between HLA-C and its respective KIR receptors determine the outcome and progression of viral infections. In this review, we focus on the unique role of HLA-C in regulating NK cell functions and its consequences in the setting of viral infections.

Like Wings of a Bird: Functional Divergence and Complementarity between HLA-A and HLA-B Molecules

Human leukocyte antigen (HLA) genes are among the most polymorphic of our genome, as a likely consequence of balancing selection related to their central role in adaptive immunity. HLA-A and HLA-B genes were recently suggested to evolve through a model of joint divergent asymmetric selection conferring all human populations, including those with severe loss of diversity, an equivalent immune potential. However, the mechanisms by which these two genes might undergo joint evolution while displaying very distinct allelic profiles in populations are still unknown. To address this issue, we carried out extensive data analyses (among which factorial correspondence analysis and linear modeling) on 2,909 common and rare HLA-A, HLA-B, and HLA-C alleles and 200,000 simulated pathogenic peptides by taking into account sequence variation, predicted peptide-binding affinity and HLA allele frequencies in 123 populations worldwide. Our results show that HLA-A and HLA-B (but not HLA-C) molecules maintain considerable functional divergence in almost all populations, which likely plays an instrumental role in their immune defense. We also provide robust evidence of functional complementarity between HLA-A and HLA-B molecules, which display asymmetric relationships in terms of amino acid diversity at both inter- and intraprotein levels and in terms of promiscuous or fastidious peptide-binding specificities. Like two wings of a flying bird, the functional complementarity of HLA-A and HLA-B is a perfect example, in our genome, of duplicated genes sharing their capacity of assuming common vital functions while being submitted to complex and sometimes distinct environmental pressures.

Haplotype Shuffling and Dimorphic Transposable Elements in the Human Extended Major Histocompatibility Complex Class II Region

The major histocompatibility complex (MHC) on chromosome 6p21 is one of the most single-nucleotide polymorphism (SNP)-dense regions of the human genome and a prime model for the study and understanding of conserved sequence polymorphisms and structural diversity of ancestral haplotypes/conserved extended haplotypes. This study aimed to follow up on a previous analysis of the MHC class I region by using the same set of 95 MHC haplotype sequences downloaded from a publicly available BioProject database at the National Center for Biotechnology Information to identify and characterize the polymorphic human leukocyte antigen (HLA)-class II genes, the MTCO3P1 pseudogene alleles, the indels of transposable elements as haplotypic lineage markers, and SNP-density crossover (XO) loci at haplotype junctions in DNA sequence alignments of different haplotypes across the extended class II region (∼1 Mb) from the telomeric PRRT1 gene in class III to the COL11A2 gene at the centromeric end of class II. We identified 42 haplotypic indels (20 Alu, 7 SVA, 13 LTR or MERs, and 2 indels composed of a mosaic of different transposable elements) linked to particular HLA-class II alleles. Comparative sequence analyses of 136 haplotype pairs revealed 98 unique XO sites between SNP-poor and SNP-rich genomic segments with considerable haplotype shuffling located in the proximity of putative recombination hotspots. The majority of XO sites occurred across various regions including in the vicinity of MTCO3P1 between HLA-DQB1 and HLA-DQB3, between HLA-DQB2 and HLA-DOB, between DOB and TAP2, and between HLA-DOA and HLA-DPA1, where most XOs were within a HERVK22 sequence. We also determined the genomic positions of the PRDM9-recombination suppression sequence motif ATCCATG/CATGGAT and the PRDM9 recombination activation partial binding motif CCTCCCCT/AGGGGAG in the class II region of the human reference genome (NC_ 000006) relative to published meiotic recombination positions. Both the recombination and anti-recombination PRDM9 binding motifs were widely distributed throughout the class II genomic regions with 50% or more found within repeat elements; the anti-recombination motifs were found mostly in L1 fragmented repeats. This study shows substantial haplotype shuffling between different polymorphic blocks and confirms the presence of numerous putative ancestral recombination sites across the class II region between various HLA class II genes.

Major histocompatibility complex genes and locus organization in the Komodo dragon (Varanus komodoensis)

We performed a meta-analysis of the newly assembled Komodo dragon (Varanus komodoensis) genome to characterize the major histocompatibility complex (MHC) of the species. The MHC gene clusters of the Komodo dragon are gene dense, complex, and contain counterparts of many genes of the human MHC. Our analysis identified 20 contigs encompassing ~ 6.9 Mbp of sequence with 223 annotated genes of which many are predicted orthologs to the genes of the human MHC. These MHC contigs range in size from 13.2 kb to 21.5 Mbp, contain an average of one gene per 30 kb, and are thought to occur on at least two chromosomes. Eight contigs, each > 100 kb, could be aligned to the human MHC based on gene content, and these represent gene clusters found in each of the recognized mammalian MHC subregions. The MHC of the Komodo dragon shares organizational features of other non-mammalian taxa. Multiple class Iα and class IIβ genes are indicated, with linkage between classical class I and immunoproteasome genes and between framework class I genes and genes associated with the mammalian class III subregion. These findings are supported in both Komodo genome assemblies and provide new insight into the MHC organization of these unique squamate reptiles.

Viruses, cancer and non-self recognition

Virus-host interactions form an essential part of every aspect of life, and this review is aimed at looking at the balance between the host and persistent viruses with a focus on the immune system. The virus-host interaction is like a cat-and-mouse game and viruses have developed ingenious mechanisms to manipulate cellular pathways, most notably the major histocompatibility (MHC) class I pathway, to reside within infected cell while evading detection and destruction by the immune system. However, some of the signals sensing and responding to viral infection are derived from viruses and the fact that certain viruses can prevent the infection of others, highlights a more complex coexistence between the host and the viral microbiota. Viral immune evasion strategies also illustrate that processes whereby cells detect and present non-self genetic material to the immune system are interlinked with other cellular pathways. Immune evasion is a target also for cancer cells and a more detailed look at the interfaces between viral factors and components of the MHC class I peptide-loading complex indicates that these interfaces are also targets for cancer mutations. In terms of the immune checkpoint, however, viral and cancer strategies appear different.

SNP-Density Crossover Maps of Polymorphic Transposable Elements and HLA Genes Within MHC Class I Haplotype Blocks and Junction

The genomic region (~4 Mb) of the human major histocompatibility complex (MHC) on chromosome 6p21 is a prime model for the study and understanding of conserved polymorphic sequences (CPSs) and structural diversity of ancestral haplotypes (AHs)/conserved extended haplotypes (CEHs). The aim of this study was to use a set of 95 MHC genomic sequences downloaded from a publicly available BioProject database at NCBI to identify and characterise polymorphic human leukocyte antigen (HLA) class I genes and pseudogenes, MICA and MICB, and retroelement indels as haplotypic lineage markers, and single-nucleotide polymorphism (SNP) crossover loci in DNA sequence alignments of different haplotypes across the Olfactory Receptor (OR) gene region (~1.2 Mb) and the MHC class I region (~1.8 Mb) from the GPX5 to the MICB gene. Our comparative sequence analyses confirmed the identity of 12 haplotypic retroelement markers and revealed that they partitioned the HLA-A/B/C haplotypes into distinct evolutionary lineages. Crossovers between SNP-poor and SNP-rich regions defined the sequence range of haplotype blocks, and many of these crossover junctions occurred within particular transposable elements, lncRNA, OR12D2, MUC21, MUC22, PSORS1A3, HLA-C, HLA-B, and MICA. In a comparison of more than 250 paired sequence alignments, at least 38 SNP-density crossover sites were mapped across various regions from GPX5 to MICB. In a homology comparison of 16 different haplotypes, seven CEH/AH (7.1, 8.1, 18.2, 51.x, 57.1, 62.x, and 62.1) had no detectable SNP-density crossover junctions and were SNP poor across the entire ~2.8 Mb of sequence alignments. Of the analyses between different recombinant haplotypes, more than half of them had SNP crossovers within 10 kb of LTR16B/ERV3-16A3_I, MLT1, Charlie, and/or THE1 sequences and were in close vicinity to structurally polymorphic Alu and SVA insertion sites. These studies demonstrate that (1) SNP-density crossovers are associated with putative ancestral recombination sites that are widely spread across the MHC class I genomic region from at least the telomeric OR12D2 gene to the centromeric MICB gene and (2) the genomic sequences of MHC homozygous cell lines are useful for analysing haplotype blocks, ancestral haplotypic landscapes and markers, CPSs, and SNP-density crossover junctions.

A fast and reliable method for detecting SNP rs67384697 (Hsa-miR-148a binding site) by a single run of allele-specific real-time PCR

Surface expression of human leukocyte antigen (HLA)-class I molecules is critical for modulating T/natural killer lymphocytes' effector functions. Among HLA molecules, HLA-C, the most recently evolved form of class I antigens, is subjected to both transcriptional and multiple post-transcriptional regulation mechanisms affecting its cell surface expression. Among the latter a region placed in the 3' untranslated region of HLA-C transcript contains the single nucleotide polymorphism (SNP) rs67384697 "G-ins/del" that has been found to be strictly associated with surface levels of HLA-C allomorphs because of the effect on the binding site of a microRNA (Hsa-miR-148a). Higher expression of HLA-C has been proved to influence HIV-1 infection via a better control of viremia and a slower disease progression. More importantly, the analysis of SNP rs67384697 "G-ins/del" combined with the evaluation of the HLA-Bw4/-Bw6 C1/C2 supratype, as well as the killer immunoglobulin-like receptor genetic asset, has proved to be pivotal in defining the status of Elite Controllers in the Caucasian population. Here we describe a new reliable and fast method of allele-specific real-time PCR to monitor the integrity/disruption of the binding site of the microRNA Hsa-miR-148a in a high-throughput format that can be easily applied to studies involving large cohorts of individuals.

Molecular evolution of elements controlling HLA-C expression: Adaptation to a role as a killer-cell immunoglobulin-like receptor ligand regulating natural killer cell function

The regulatory elements controlling the transcription of the HLA-A, HLA-B, and HLA-C genes have been extensively studied and compared. However, few studies have considered regulatory differences in the HLA genes from the perspective of their role as ligands for the killer-cell immunoglobulin-like receptor (KIR) family of HLA receptors expressed by natural killer (NK) cells. HLA-C is the most recently evolved gene, and there is considerable evidence pointing to its emergence as a specialized KIR ligand playing a major role in the missing-self recognition system of NK cells. Here I evaluate gene-specific differences in regulatory elements of the HLA genes, showing alterations that are consistent with the adaptation of HLA-C to a role in NK cell regulation.

Probing the Association between Early Evolutionary Markers and Schizophrenia

Schizophrenia is suggested to be a by-product of the evolution in humans, a compromise for our language, creative thinking and cognitive abilities, and thus, essentially, a human disorder. The time of its origin during the course of human evolution remains unclear. Here we investigate several markers of early human evolution and their relationship to the genetic risk of schizophrenia. We tested the schizophrenia evolutionary hypothesis by analyzing genome-wide association studies of schizophrenia and other human phenotypes in a statistical framework suited for polygenic architectures. We analyzed evolutionary proxy measures: human accelerated regions, segmental duplications, and ohnologs, representing various time periods of human evolution for overlap with the human genomic loci associated with schizophrenia. Polygenic enrichment plots suggest a higher prevalence of schizophrenia associations in human accelerated regions, segmental duplications and ohnologs. However, the enrichment is mostly accounted for by linkage disequilibrium, especially with functional elements like introns and untranslated regions. Our results did not provide clear evidence that markers of early human evolution are more likely associated with schizophrenia. While SNPs associated with schizophrenia are enriched in HAR, Ohno and SD regions, the enrichment seems to be mediated by affiliation to known genomic enrichment categories. Taken together with previous results, these findings suggest that schizophrenia risk may have mainly developed more recently in human evolution.

Endogenized viral sequences in mammals

Reverse-transcribed RNA molecules compose a significant portion of the human genome. Many of these RNA molecules were retrovirus genomes either infecting germline cells or having done so in a previous generation but retaining transcriptional activity. This mechanism itself accounts for a quarter of the genomic sequence information of mammals for which there is data. We understand relatively little about the causes and consequences of retroviral endogenization. This review highlights functions ascribed to sequences of viral origin endogenized into mammalian genomes and suggests some of the most pressing questions raised by these observations.

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

Background

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

Methods

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

Results

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

Conclusions

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

Electronic supplementary material

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

Origin and evolution of yolk proteins: expansion and functional diversification of large lipid transfer protein superfamily

Vitellogenin (VTG) and apolipoprotein (APO) play a central role in animal reproduction and lipid circulation, respectively. Although previous studies have examined the structural and functional relationships of these large lipid transfer proteins (LLTPs) from an evolutionary perspective, the mechanism in generating these different families have not been addressed in invertebrates. In this study, the most comprehensive phylogenetic and genomic analysis of the LLTP superfamily genes is carried out. We propose the expansion and diversification of LLTPs in invertebrates are mediated via retrotransposon-mediated duplications, followed by either subfunctionalization or neofunctionalization in different lineages. In agreement with a previous hypothesis, our analysis suggests that all LLTPs originate from a series of duplications of a primitive yolk protein gene similar to VTG. Two early consecutive duplications of the yolk protein genes resulted in the formation of microsomal triglyceride transfer protein (MTP) and the APO gene ancestor. Gains and losses of domains and genes occurred in each of these families in different animal lineages, with MTP becoming truncated. MTP maintained only the components stabilizing the huge lipoprotein particle. Surprisingly, for the first time, two VTG-like protein families were found to independently arise in the lineages of insects. This work consolidates the reconstruction of the evolutionary roadmap of the LLTP superfamily and provides the first mechanistic explanation on the expansion of family members via retrotransposition in invertebrates.

The role of duplications in the evolution of genomes highlights the need for evolutionary-based approaches in comparative genomics

Understanding the evolutionary plasticity of the genome requires a global, comparative approach in which genetic events are considered both in a phylogenetic framework and with regard to population genetics and environmental variables. In the mechanisms that generate adaptive and non-adaptive changes in genomes, segmental duplications (duplication of individual genes or genomic regions) and polyploidization (whole genome duplications) are well-known driving forces. The probability of fixation and maintenance of duplicates depends on many variables, including population sizes and selection regimes experienced by the corresponding genes: a combination of stochastic and adaptive mechanisms has shaped all genomes. A survey of experimental work shows that the distinction made between fixation and maintenance of duplicates still needs to be conceptualized and mathematically modeled. Here we review the mechanisms that increase or decrease the probability of fixation or maintenance of duplicated genes, and examine the outcome of these events on the adaptation of the organisms.

Reviewers

This article was reviewed by Dr. Etienne Joly, Dr. Lutz Walter and Dr. W. Ford Doolittle.

Genomic plasticity of the MHC class I A region in rhesus macaques: extensive haplotype diversity at the population level as revealed by microsatellites

The Mamu-A genes of the rhesus macaque show different degrees of polymorphism, transcription level variation, and differential haplotype distribution. Per haplotype, usually one “major” transcribed gene is present, A1 (A7), in various combinations with “minor” genes, A2 to A6. In silico analysis of the physical map of a heterozygous animal revealed the presence of similar Mamu-A regions consisting of four duplication units, but with dissimilar positions of the A1 genes on both haplotypes, and in combination with different minor genes. Two microsatellites, D6S2854 and D6S2859, have been selected as potential tools to characterize this complex region. Subsequent analysis of a large breeding colony resulted in the description of highly discriminative patterns, displaying copy number variation in concert with microsatellite repeat length differences. Sequencing and segregation analyses revealed that these patterns are unique for each Mamu-A haplotype. In animals of Indian, Burmese, and Chinese origin, 19, 15, or 9 haplotypes, respectively, could be defined, illustrating the occurrence of differential block duplications and subsequent rearrangements by recombination. The haplotypes can be assigned to 12 unique combinations of genes (region configurations). Although most configurations harbor two transcribed A genes, one or three genes per haplotype are also present. Additionally, haplotypes lacking an A1 gene or with an A1 duplication appear to exist. The presence of different transcribed A genes/alleles in monkeys from various origins may have an impact on differential disease susceptibilities. The high-throughput microsatellite technique will be a valuable tool in animal selection for diverse biomedical research projects.

The major histocompatibility complex (Mhc) class IIB region has greater genomic structural flexibility and diversity in the quail than the chicken

Background

The quail and chicken major histocompatibility complex (Mhc) genomic regions have a similar overall organization but differ markedly in that the quail has an expanded number of duplicated class I, class IIB, natural killer (NK)-receptor-like, lectin-like and BG genes. Therefore, the elucidation of genetic factors that contribute to the greater Mhc diversity in the quail would help to establish it as a model experimental animal in the investigation of avian Mhc associated diseases.

Aims and approaches

The main aim here was to characterize the genetic and genomic features of the transcribed major quail MhcIIB (CojaIIB) region that is located between the Tapasin and BRD2 genes, and to compare our findings to the available information for the chicken MhcIIB (BLB). We used four approaches in the study of the quail MhcIIB region, (1) haplotype analyses with polymorphic loci, (2) cloning and sequencing of the RT-PCR CojaIIB products from individuals with different haplotypes, (3) genomic sequencing of the CojaIIB region from the individuals with the different haplotypes, and (4) phylogenetic and duplication analysis to explain the variability of the region between the quail and the chicken.

Results

Our results show that the Tapasin-BRD2 segment of the quail Mhc is highly variable in length and in gene transcription intensity and content. Haplotypic sequences were found to vary in length between 4 to 11 kb. Tapasin-BRD2 segments contain one or two major transcribed CojaIIBs that were probably generated by segmental duplications involving c-type lectin-like genes and NK receptor-like genes, gene fusions between two CojaIIBs and transpositions between the major and minor CojaIIB segments. The relative evolutionary speed for generating the MhcIIBs genomic structures from the ancestral BLB2 was estimated to be two times faster in the quail than in the chicken after their separation from a common ancestor. Four types of genomic rearrangement elements (GRE), composed of simple tandem repeats (STR), were identified in the MhcIIB genomic segment located between the Tapasin-BRD2 genes. The GREs have many more STR numbers in the quail than in the chicken that displays strong linkage disequilibrium.

Conclusion

This study suggests that the Mhc classIIB region has a flexible genomic structure generated by rearrangement elements and rapid SNP accumulation probably as a consequence of the quail adapting to environmental conditions and pathogens during its migratory history after its divergence from the chicken.

Host genetics and viral infections: immunology taught by viruses, virology taught by the immune system

Viruses are prototypic obligate intracellular pathogens, and are therefore, by necessity, highly pre-adapted to surviving the host immune response. Nevertheless, host genetic factors remain an important determinant of disease outcome, particularly in the case of viruses that have encountered humans in the more recent past (e.g. human immunodeficiency virus). Recent studies have identified an increasingly rich network of functionally relevant polymorphic immune factors, including major MHC alleles, killer immunoglobulin-like receptors and functional chemokine receptor polymorphisms. Moreover, genetic variation is increasingly appreciated beyond the single genotype level, incorporating extended haplotypes as well as regions of segmental genetic duplication. These issues can be considered within an evolutionary perspective that acknowledges the crucial role of adaptive host-viral relationships in shaping both host and pathogen genetic diversity.

Impact of transposable elements on the evolution of mammalian gene regulation

Transposable elements (TEs) are present in all organisms and nearly half of the human and mouse genome is derived from ancient transpositions. This fact alone suggests that TEs have played a major role in genome organization and evolution. Studies undertaken over the last two decades or so clearly show that TEs of various kinds have played an important role in organism evolution. Here we review the impact TEs have on the evolution of gene regulation and gene function with an emphasis on humans. Understanding the mechanisms resulting in genomic change is central to our understanding of gene regulation, genetic disease and genome evolution. Full comprehension of these biological processes is not possible without an in depth knowledge of how TEs impact upon the genome.

ERVK9, transposons and the evolution of MHC class I duplicons within the alpha-block of the human and chimpanzee

The genomic sequences within the alpha-block (approximately 288-310 kb) of the human and chimpanzee MHC class I region contains ten MHC class I genes and three MIC gene fragments grouped together within alternating duplicated genomic segments or duplicons. In this study, the chimpanzee and human genomic sequences were analyzed in order to determine whether the remnants of the ERVK9 and other retrotransposon sequences are useful genomic markers for reconstructing the evolutionary history of the duplicated MHC gene families within the alpha-block. A variety of genes, pseudogenes, autologous DNA transposons and retrotransposons such as Alu and ERVK9 were used to categorize the ten duplicons into four distinct structural groups. The phylogenetic relationship of the ten duplicons was examined by using the neighbour joining method to analyze transposon sequence topologies of selected Alu members, LTR16B and Charlie9. On the basis of these structural groups and the phylogeny of the duplicated transposon sequences, a duplication model was reconstructed involving four multipartite tandem duplication steps to explain the organization and evolution of the ten duplicons within the alpha-block of the chimpanzee and human. The phylogenetic analysis and inferred duplication history suggests that the Patr/HLA-F was the first MHC class I gene to have been fixed and not required as a precursor for further duplication within the alpha-block of the ancestral species.

Indels and imperfect duplication have driven the evolution of human Complement Receptor 1 (CR1) and CR1-like from their precursor CR1 alpha: importance of functional sets

This study examines the effects of duplication and insertions-deletions (indels) by comparing human complement receptor 1 (CR1) and human CR1-like (CR1L) with syntenic genes from four other vertebrates (chimpanzee, baboon, rat, and mouse). By phylogenetic analysis, the domains of these genes can be classified into 10 distinct subfamilies (a, b, c, d, e, f, g(-like), h, j, and k), which have been largely conserved throughout vertebrate and invertebrate evolution. In spite of many complex and diverse duplications and indels, the subfamily order of domains (a, j, e, f, b, k, d, g(-like)) has been maintained. The number of domain sets has increased progressively, thereby expanding the functional repertoire.

Genomic evolution of MHC class I region in primates

To elucidate the origins of the MHC-B-MHC-C pair and the MHC class I chain-related molecule (MIC)A-MICB pair, we sequenced an MHC class I genomic region of humans, chimpanzees, and rhesus monkeys and analyzed the regions from an evolutionary stand-point, focusing first on LINE sequences that are paralogous within each of the first two species and orthologous between them. Because all the long interspersed nuclear element (LINE) sequences were fragmented and nonfunctional, they were suitable for conducting phylogenetic study and, in particular, for estimating evolutionary time. Our study has revealed that MHC-B and MHC-C duplicated 22.3 million years (Myr) ago, and the ape MICA and MICB duplicated 14.1 Myr ago. We then estimated the divergence time of the rhesus monkey by using other orthologous LINE sequences in the class I regions of the three primate species. The result indicates that rhesus monkeys, and possibly the Old World monkeys in general, diverged from humans 27-30 Myr ago. Interestingly, rhesus monkeys were found to have not the pair of MHC-B and MHC-C but many repeated genes similar to MHC-B. These results support our inference that MHC-B and MHC-C duplicated after the divergence between apes and Old World monkeys.

Structural evolution of the BRCA1 genomic region in primates

Segmental duplications account for up to 6% of the human genome, and the resulting low-copy repeats (LCRs) are known to be associated with more than 20 genomic disorders. Many such duplication events coincided with the burgeoning of the Alu repeat family during the last 50 million years of primate evolution, and it has been suggested that the two phenomena might be causally related. In tracing the evolution of the BRCA1 17q21 region through the primate clade, we discovered the occurrence over the last 40 million years of a complex set of about eight large gene-conversion-mediated rearrangements in the approximately 4 Mb surrounding the BRCA1 gene. These have resulted in the presence of large and probably recombinogenic LCRs across the region, the creation of the NBR2 gene, the duplication of the BRCA1/NBR1 promoter, the bisection of the highly conserved ARF2 gene, and multiple copies of the KIAA0563 gene. The junctions lie within AluS repeats, members of an Alu subfamily which experienced massive expansion during the time that the rearrangements occurred. We present a detailed history of this region over a critical 40 million-year period of genomic upheaval, including circumstantial evidence for a causal link between Alu family expansion and the rearrangement-mediated destruction and creation of transcription units.

Evolutionary relationships of major histocompatibility complex class I genes in simian primates

New World monkeys (NWMs) occupy a critical phylogenetic position in elucidating the evolutionary process of major histocompatibility complex (MHC) class I genes in primates. From three subfamilies of Aotinae, Cebinae, and Atelinae, the 5'-flanking regions of 18 class I genes are obtained and phylogenetically examined in terms of Alu/LINE insertion elements as well as the nucleotide substitutions. Two pairs of genes from Aotinae and Atelinae are clearly orthologous to human leukocyte antigen (HLA) -E and -F genes. Of the remaining 14 genes, 8 belong to the distinct group B, together with HLA-B and -C, to the exclusion of all other HLA class I genes. These NWM genes are classified into four groups, designated as NWM-B1, -B2, -B3, and -B4. Of these, NWM-B2 is orthologous to HLA-B/C. Also, orthologous relationships of NWM-B1, -B2, and -B3 exist among different families of Cebidae and Atelidae, which is in sharp contrast to the genus-specific gene organization within the subfamily Callitrichinae. The other six genes belong to the distinct group G. However, a clade of these NWM genes is almost equally related to HLA-A, -J, -G, and -K, and there is no evidence for their orthologous relationships to HLA-G. It is argued that class I genes in simian primates duplicated extensively in their common ancestral lineage and that subsequent evolution in descendant species has been facilitated mainly by independent loss of genes.

Analysis of a marsupial MHC region containing two recently duplicated class I loci

A 37-kb cosmid containing two complete major histocompatibility complex (MHC) class I alpha chain loci from the opossum Monodelphis domestica was isolated, fully sequenced, and characterized. This sequence represents the largest contiguous genomic sequence reported for the MHC region of a nonplacental mammal. Based on particular conserved amino acid residues, and limited expression analyses, the two MHC-I loci, designated ModoUB and ModoUC, appear to encode functional MHC-I molecules. The two coding regions are 98% identical at the nucleotide level; however, their promoter regions differ significantly. Two CpG islands present in the cosmid sequence correspond to the two coding regions. Twelve microsatellites and six retroelements were also present in the cosmid. The retroelements share highest sequence homology to the CORE-SINE family of retroelements. Due to high sequence identity, it is very likely that ModoUB and ModoUC loci are products of recent gene duplication that occurred less than 4 million years ago.

High similarity between flanking regions of different microsatellites detected within each of two species of Lepidoptera: Parnassius apollo and Euphydryas aurinia

Microsatellite flanking regions have been compared in two butterfly species. Several microsatellite flanking regions showed high similarity to one another among different microsatellites within a same species, but very few similarities were found between species. This can be the consequence of either duplication/multiplication events involving large regions containing microsatellites or of microsatellites imbedded in minisatellite regions. The multiplication of microsatellites might also be linked to mobile elements. Furthermore, crossing over between nonhomologous microsatellites can lead to the exchange of the flanking regions between microsatellites. The same phenomenon was observed in both studied butterfly species but not in Aphis fabae (Hemiptera), which was screened at the same time using the same protocol. These findings might explain, at least partially, why microsatellite isolation in Lepidoptera has been relatively unsuccessful so far.

Evolutionary Relationships of Major Histocompatibility Complex Class I Genes in Simian Primates

New World monkeys (NWMs) occupy a critical phylogenetic position in elucidating the evolutionary process of major histocompatibility complex (MHC) class I genes in primates. From three subfamilies of Aotinae, Cebinae, and Atelinae, the 5′-flanking regions of 18 class I genes are obtained and phylogenetically examined in terms of Alu/LINE insertion elements as well as the nucleotide substitutions. Two pairs of genes from Aotinae and Atelinae are clearly orthologous to human leukocyte antigen (HLA) -E and -F genes. Of the remaining 14 genes, 8 belong to the distinct group B, together with HLA-B and -C, to the exclusion of all other HLA class I genes. These NWM genes are classified into four groups, designated as NWM-B1, -B2, -B3, and -B4. Of these, NWM-B2 is orthologous to HLA-B/C. Also, orthologous relationships of NWM-B1, -B2, and -B3 exist among different families of Cebidae and Atelidae, which is in sharp contrast to the genus-specific gene organization within the subfamily Callitrichinae. The other six genes belong to the distinct group G. However, a clade of these NWM genes is almost equally related to HLA-A, -J, -G, and -K, and there is no evidence for their orthologous relationships to HLA-G. It is argued that class I genes in simian primates duplicated extensively in their common ancestral lineage and that subsequent evolution in descendant species has been facilitated mainly by independent loss of genes.

IL-10 mediates suppression of the CD8 T cell IFN-gamma response to a novel viral epitope in a primed host

Priming to Ag can inhibit subsequent induction of an immune response to a new epitope incorporated into that Ag, a phenomenon referred to as original antigenic sin. In this study, we show that prior immunity to a virus capsid can inhibit subsequent induction of the IFN-gamma effector T cell response to a novel CD8-restricted antigenic epitope associated with the virus capsid. Inhibition does not involve Ab to the virus capsid, as it is observed in animals lacking B cells. CD8-restricted virus-specific T cell responses are not required, as priming to virus without CTL induction is associated with inhibition. However, IL-10(-/-) mice, in contrast to IL-10(+/+) mice, generate CD8 T cell and Ab responses to novel epitopes incorporated into a virus capsid, even when priming to the capsid has resulted in high titer Ab to the capsid. Furthermore, capsid-primed mice, unable to mount a response to a novel epitope in the capsid protein, are nevertheless able to respond to the same novel epitope delivered independently of the capsid. Thus, inhibition of responsiveness to a novel epitope in a virus-primed animal is a consequence of secretion of IL-10 in response to presented Ag, which inhibits local generation of new CD8 IFN-gamma-secreting effector T cells. Induction of virus- or tumor Ag-specific CD8 effector T cells in the partially Ag-primed host may thus be facilitated by local neutralization of IL-10.

Association of MHC dimorphic Alu insertions with HLA class I and MIC genes in Japanese HLA-B48 haplotypes

A large proportion of Japanese with the HLA-B48 allele have a MICA gene deletion associated with a MICB null allele within the class I region of the Major Histocompatibility Complex (MHC). Here, we report for the first time a novel positive association between the presence of a polymorphic Alu insertion, AluyMICB, within the first intron of the MICB gene and the MICAdel/MICBnull/HLA-B48 haplotype for five of six well-characterized Japanese cell-lines. The AluyMICB insertion was found to be present at a frequency of 0.242 in 86 Japanese tissue donors and in four of the five individuals with the HLA-B48 allele. The AluyMICB insertion was also associated with at least three different MICB alleles, *0102, *0107N and *0105, and three different HLA-B alleles, B13, B48 and B57, respectively, in the seven Workshop cell-lines (the 4th Asia-Oceania Histocompatibility Workshop, and the 10th International Histocompatibility Workshop) and the six Japanese cell-lines that were selected for this study. Based on the analysis of associations between different polymorphic markers within the beta block, the MICB*0102 allele was inferred to be the ancestral form of the MICB*0105 and MICB*0107N alleles. The AluyMICB polymorphism can now be used to further investigate its relationship with other MICB alleles and consequently their origins. In addition, we have examined the absence and presence of three other polymorphic Alu markers distributed within the alpha block of the class I region of the HLA-B48/AluyMICB haplotype. We conclude that the extended HLA-B haplotypes are best defined by considering multiple genomic sites including the four polymorphic Alu insertions described in this study.

Genomic organization of the mammalian MHC

The Human Genome Project transformed the quest of more than 50 years to understand the major histocompatibility complex (Mhc). The sequence of the Mhc from human and mouse, together with a large amount of sequence and mapping information from several other species, allows us to draw general conclusions about the organization and origin of this crucial part of the immune system. The Mhc is a mosaic of stretches formed by conserved and nonconserved genes. Surprisingly, of the approximately 3.6-Mb Mhc, the stretches that encode the class I and class II genes, which epitomize the Mhc, are the least conserved part, whereas the approximately 1.7-Mb stretches that encode at least 115 other genes are highly conserved. We summarize the available data to answer the questions (a) What is the Mhc? and (b) How can we define it in a general, not species-specific, way? Knowing what is essential and what is incidental helps us understand the fundamentals of the Mhc, and defining the species differences makes the model organisms more useful.

Species-specific class I gene expansions formed the telomeric 1 mb of the mouse major histocompatibility complex

We have determined the complete sequence of 951,695 bp from the class I region of H2, the mouse major histocompatibility complex (Mhc) from strain 129/Sv (haplotype bc). The sequence contains 26 genes. The sequence spans from the last 50 kb of the H2-T region, including 2 class I genes and 3 class I pseudogenes, and includes the H2-M region up to Gabbr1. A 500-kb stretch of the H2-M region contains 9 class I genes and 4 pseudogenes, which fall into two subfamilies, M1 and M10, distinct from other mouse class I genes. This M1/M10 class I gene-cluster is separated from the centromeric H2-T and the telomeric H2-M4, -5 and -6 class I genes by "nonclass I genes". Comparison with the corresponding 853-kb region of the human Mhc, which includes the HLA-A region, shows a mosaic of conserved regions of orthologous nonclass I genes separated by regions of species-specific expansion of paralogous Mhc class I genes. The analysis of this mosaic structure illuminates the dynamic evolution of the Mhc class I region among mammals and provides evidence for the framework hypothesis.

The promoter of the classical MHC class I locus in rainbow trout (Oncorhynchus mykiss)

In rainbow trout only a single classical major histocompatibility complex (MHC) class I locus has been identified. In previous studies it was shown that this locus, Onmy-UBA, is characterised by extensive polymorphism and ubiquitous expression. In the present study the sequence of the promoter region of Onmy-UBA was analysed. Within this region several fragments for transcription regulation could be detected, including interferon-sensitive enhancer elements. The promoter region is very similar to that of classical MHC class I in higher vertebrates. In agreement with the interferon sensitive elements identified in its promoter, Onmy-UBA expression in rainbow trout gonad cells (RTG-2 cell line) was up-regulated after infection with infectious haematopoietic necrosis virus (IHNV).

Comparative genomic analysis of the MHC: the evolution of class I duplication blocks, diversity and complexity from shark to man

The major histocompatibility complex (MHC) genomic region is composed of a group of linked genes involved functionally with the adaptive and innate immune systems. The class I and class II genes are intrinsic features of the MHC and have been found in all the jawed vertebrates studied so far. The MHC genomic regions of the human and the chicken (B locus) have been fully sequenced and mapped, and the mouse MHC sequence is almost finished. Information on the MHC genomic structures (size, complexity, genic and intergenic composition and organization, gene order and number) of other vertebrates is largely limited or nonexistent. Therefore, we are mapping, sequencing and analyzing the MHC genomic regions of different human haplotypes and at least eight nonhuman species. Here, we review our progress with these sequences and compare the human MHC structure with that of the nonhuman primates (chimpanzee and rhesus macaque), other mammals (pigs, mice and rats) and nonmammalian vertebrates such as birds (chicken and quail), bony fish (medaka, pufferfish and zebrafish) and cartilaginous fish (nurse shark). This comparison reveals a complex MHC structure for mammals and a relatively simpler design for nonmammalian animals with a hypothetical prototypic structure for the shark. In the mammalian MHC, there are two to five different class I duplication blocks embedded within a framework of conserved nonclass I and/or nonclass II genes. With a few exceptions, the class I framework genes are absent from the MHC of birds, bony fish and sharks. Comparative genomics of the MHC reveal a highly plastic region with major structural differences between the mammalian and nonmammalian vertebrates. Additional genomic data are needed on animals of the reptilia, crocodilia and marsupial classes to find the origins of the class I framework genes and examples of structures that may be intermediate between the simple and complex MHC organizations of birds and mammals, respectively.

Songbird genomics: analysis of 45 kb upstream of a polymorphic Mhc class II gene in red-winged blackbirds (Agelaius phoeniceus)

Here we present the sequence of a 45 kb cosmid containing a previously characterized poly-morphic Mhc class II B gene (Agph-DAB1) from the red-winged blackbird (Agelaius phoeniceus). We compared it with a previously sequenced cosmid from this species, revealing two regions of 7.5 kb and 13.0 kb that averaged greater than 97% similarity to each another, indicating a very recent shared duplication. We found 12 retroelements, including two chicken repeat 1 (CR1) elements, constituting 6.4% of the sequence and indicating a lower frequency of retroelements than that found in mammalian genomic DNA. Agph-DAB3, a new class II B gene discovered in the cosmid, showed a low rate of polymorphism and may be functional. In addition, we found a Mhc class II B gene fragment and three genes likely to be functional (encoding activin receptor type II, a zinc finger, and a putative gamma-filamin). Phylogenetic analysis of exon 2 alleles of all three known blackbird Mhc genes indicated strong clustering of alleles by locus, implying that large amounts of interlocus gene conversion have not occurred since these genes have been diverging. Despite this, interspecific comparisons indicate that all three blackbird Mhc genes diverged from one another less than 35 million years ago and are subject to concerted evolution in the long term. Comparison of blackbird and chicken Mhc promoter regions revealed songbird promoter elements for the first time. The high gene density of this cosmid confirms similar findings for the chicken Mhc, but the segment duplications and diversity of retroelements resembles mammalian sequences.

Perspective: transposable elements, parasitic DNA, and genome evolution

The nature of the role played by mobile elements in host genome evolution is reassessed considering numerous recent developments in many areas of biology. It is argued that easy popular appellations such as "selfish DNA" and "junk DNA" may be either inaccurate or misleading and that a more enlightened view of the transposable element-host relationship encompasses a continuum from extreme parasitism to mutualism. Transposable elements are potent, broad spectrum, endogenous mutators that are subject to the influence of chance as well as selection at several levels of biological organization. Of particular interest are transposable element traits that early evolve neutrally at the host level but at a later stage of evolution are co-opted for new host functions.

Identification and characterization of a beta proteasome subunit cluster in the Japanese pufferfish (Fugu rubripes)

The low molecular mass polypeptide (LMP2, LMP7, and MECL-1) genes code for beta-type subunits of the proteasome, a multimeric complex that degrades proteins into peptides as part of the MHC class I-mediated Ag-presenting pathway. These gene products are up-regulated in response to infection by IFN-gamma and replace the corresponding constitutively expressed subunits (X, Y, and Z) during the immune response. In humans, the LMP2 and LMP7 genes both reside within the class II region of the MHC (6p21.3), while MECL-1 is located at 16q22.1. In the present study, we have identified all three IFN-gamma-regulated beta-type proteasome subunits in Fugu, which are present as a cluster within the Fugu MHC class I region. We show that in this species, LMP7, LMP2, and MECL-1 are linked. Also within this cluster is an LMP2-like subunit (which seems specific to all teleosts tested to date) and a closely linked LMP7 pseudogene, indicating that within Fugu and potentially other teleosts, there has been an additional regional duplication involving these genes.

Impact of transposable elements on the human genome

Presence of transposable elements (TEs) in the human genome has profound effects on genome function, structure and evolution. TE mobility and inter-TE recombination are the origin of a large spectrum of mutations and genome reorganization leading to diseases. From the data provided by the Human Genome Project and from information on the detection and dynamics of TEs within and between species acquired during the last two decades, we now know that these elements are not only involved in mutagenesis but can also participate in many cellular functions including recombination, gene regulation, protein-coding RNA messages and, possibly, cellular stress response and centromere function. TEs also promote a general genome shuffling process that has been important for the evolution of several gene families and for the development of new regulatory pathways.

Structure and polymorphism of two stress-activated protein kinase genes centromeric of the MHC: SAPK2a and SAPK4

As MHC genes are potent determinants of susceptibility to immunopathological diseases, the mapping of SAPK2a (CSBP) and SAPK4 to chromosome 6p 21.2-21.3 suggested that these genes may mediate the effects of the MHC on disease. Here we describe the genomic structure and localisation of both genes approximately 2.3Mb centromeric of HLA-DP. Examination of the complete coding region and selected intronic regions of SAPK2a and SAPK4 from 22 human EBV-transformed B-cell lines of different MHC haplotypes and racial background revealed complete sequence conservation. There were no notable differences in levels of expression of SAPK2a and SAPK4 mRNA in cell lines of different MHC haplotypes or racial origin. Examination of the SAPK2a and SAPK4 sequences from two chimpanzees revealed 3 nucleotide differences between human and chimpanzee in each gene resulting in only one amino acid change in SAPK4, and 6 nucleotide substitutions plus 2 deletions in 600bp of intronic sequence from SAPK4. This highlights the selective pressure placed on these genes to maintain their protein sequence, but does not favour a role in genetic regulation of disease or provide evidence of linkage disequilibrium with the MHC.

Different Evolutionary Histories in Two Subgenomic Regions of the Major Histocompatibility Complex

Two subgenomic regions within the major histocompatibility complex, the alpha and beta blocks, contain members of the multicopy gene families HLA class I, human endogenous retroviral sequence (HERV-16; previously known as P5 and PERB3), hemochromatosis candidate genes (HCG) (II, IV, VIII, IX), 3.8-1, and MIC (PERB11). In this study we show that the two blocks consist of imperfect duplicated segments, which contain linked members of the different gene families. The duplication and truncation sites of the segments are associated with retroelements. The retroelement sites appear to generate the imperfect duplications, insertions/deletions, and rearrangements, most likely via homologous recombination. Although the two blocks share several characteristics, they differ in the number and orientation of the duplicated segments. On the 62.1 haplotype, the alpha block consists of at least 10 duplicated segments that predominantly contain pseudogenes and gene fragments of the HLA class I and MIC (PERB11) gene families. In contrast, the beta block has two major duplications containing the genes HLA-B and HLA-C, and MICA(PERB11.1) and MICB(PERB11.2). Given the common origin between the blocks, we reconstructed the duplication history of the segments to understand the processes involved in producing the different organization in the two blocks. We then found that the beta block contains four distinct duplications from two separate events, whereas the alpha block is characterized by multisegment duplications. We will discuss these results in relation to the genetic content of the two blocks.

Genomics of the major histocompatibility complex: haplotypes, duplication, retroviruses and disease

The genomic region encompassing the Major Histocompatibility Complex (MHC) contains polymorphic frozen blocks which have developed by local imperfect sequential duplication associated with insertion and deletion (indels). In the alpha block surrounding HLA-A, there are ten duplication units or beads on the 62.1 ancestral haplotype. Each bead contains or contained sequences representing Class I, PERB11 (MHC Class I chain related (MIC) and human endogenous retrovirus (HERV) 16. Here we consider explanations for co-occurrence of genomic polymorphism, duplication and HERVs and we ask how these features encode susceptibility to numerous and very diverse diseases. Ancestral haplotypes differ in their copy number and indels in addition to their coding regions. Disease susceptibility could be a function of all of these differences. We propose a model of the evolution of the human MHC. Population-specific integration of retroviral sequences could explain rapid diversification through duplication and differential disease susceptibility. If HERV sequences can be protective, there are exciting prospects for manipulation. In the meanwhile, it will be necessary to understand the function of MHC genes such as PERB11 (MIC) and many others discovered by genomic sequencing.

Duplication and polymorphism in the MHC: Alu generated diversity and polymorphism within the PERB11 gene family

The PERB11 gene family has at least five members within the telomeric region of the MHC. The PERB11.1 and PERB11.2 genes are approximately 40 kb and 160 kb centromeric of HLA-B, respectively. Using continuous genomic sequence encompassing PERB11.1 and PERB11.2, we have found a large (approximately 25 kb) segmental duplication extending beyond the genes themselves and other potential coding sequences. The major difference between the segments are large indels which are predominantly Alu sequences. The Alu sequences within the duplicated segments have created diversity via the internal and 3' poly A-rich region. A sequence comparison of an Alu sequence between two different human ancestral haplotypes shows a high level of polymorphism, particularly in the poly A-rich regions. This study characterises the Alu sequences within the peri-PERB11.1 and peri-PERB11.2 duplicated segments in relation to diversity and polymorphism and as evolutionary markers.