Rhesus macaque class I duplicon structures, organization, and evolution within the alpha block of the major histocompatibility complex

The alpha block of the human and chimpanzee major histocompatibility complex (MHC) class I genomic region contains 10 to 11 duplicated MHC class I genes, including the HLA/Patr-A, -G, and -F genes. In comparison, the alpha block of the rhesus macaque (Macaca mulatta, Mamu) has an additional 20 MHC class I genes within this orthologous region. The present study describes the identification and analysis of the duplicated segmental genomic structures (duplicons) and genomic markers within the alpha block of the rhesus macaque and their use to reconstruct the duplication history of the genes within this region. A variety of MHC class I genes, pseudogenes, transposons, and retrotransposons, such as Alu and ERV16, were used to categorize the 28 duplicons into four distinct structural categories. The phylogenetic relationship of MHC class I genes, Alu, and LTR16B sequences within the duplicons was examined by use of the Neighbor-Joining (NJ) method. Two single-duplicon tandem duplications, two polyduplicon tandem duplications with an accompanying inversion product per duplication, eight polyduplicon tandem duplications steps, 12 deletions, and at least two recombinations were reconstructed to explain the highly complex organization and evolution of the 28 duplicons (nine inversions) within the Mamu alpha block. On the basis of the phylogenetic evidence and the reconstructed tandem duplication history of the 28 duplicons, the Mamu/Patr/HLA-F ortholog was the first MHC class I gene to have been fixed without further duplication within the alpha block of primates. Assuming that the rhesus macaque and the chimpanzee/human lineages had started with the same number of MHC class I duplicons at the time of their divergence approximately 24 to 31 MYA, then the number of genes within the alpha block have been duplicated at an approximately three times greater rate in the rhesus macaque than in either the human or chimpanzee.     

Genomic and Phylogenetic Analysis of the Human CD1 and HLA Class I Multicopy Genes

The human CD1 proteins belong to a lipid-glycolipid antigen-presenting gene family and are related in structure and function to the MHC class I molecules. Previous mapping and DNA hybridization studies have shown that five linked genes located within a cluster on human chromosome 1q22-23 encode the CD1 protein family. We have analyzed the complete genomic sequence of the human CD1 gene cluster and found that the five active genes are distributed over 175,600 nucleotides and separated by four expanded intervening genomic regions (IGRs) ranging in length between 20 and 68 kb. The IGRs are composed mostly of retroelements including five full-length L1 PA sequences and various pseudogenes. Some L1 sequences have acted as receptors for other subtypes or families of retroelements. Alu molecular clocks that have evolved during primate history are found distributed within the HLA class I duplicated segments (duplicons) but not within the duplicons of CD1. Phylogeny of the alpha3 domain of the class I-like superfamily of proteins shows that the CD1 cluster is well separated from HLA class I by a number of superfamily members including MIC (PERB11), HFE, Zn-alpha2-GP, FcRn, and MR1. Phylogenetically, the human CD1 sequences are interspersed by CD1 sequences from other mammalian species, whereas the human HLA class I sequences cluster together and are separated from the other mammalian sequences. Genomic and phylogenetic analyses support the view that the human CD1 gene copies were duplicated prior to the evolution of primates and the bulk of the HLA class I genes found in humans. In contrast to the HLA class I genomic structure, the human CD1 duplicons are smaller in size, they lack Alu clocks, and they are interrupted by IGRs at least 4 to 14 times longer than the CD1 genes themselves. The IGRs seem to have been created as "buffer zones" to protect the CD1 genes from disruption by transposable elements.     

Using Alu J Elements as Molecular Clocks to Trace the Evolutionary Relationships Between Duplicated HLA Class I Genomic Segments

The class I region of the major histocompatibility complex contains two subgenomic blocks (250–350 kb each), known as the alpha and beta blocks. These blocks contain members of multicopy gene families including HLA class I, HERV-16 (previously called P5 sequences), and PERB11 (MIC). We have previously shown that each block consists of imperfect duplicated segments (duplicons) containing linked members of different gene families, retroelements and transposons that have coevolved as part of two separate evolutionary events. Another region provisionally designated here as the kappa block is located between the alpha and the beta blocks and contains HLA-E, -30, and -92, HERV-16 (P5.3), and PERB11.3 (MICC) within about 250 kb of sequence. Using Alu elements to trace the evolutionary relationships between different class I duplicons, we have found that (a) the kappa block contains paralogous (duplicated) Alu J sequences and other retroelement patterns more in common with the beta than the alpha block; (b) the retroelement pattern associated with the HLA-E duplicon is different from all other HLA class I duplicons, indicating a more complex evolution; (c) the HLA-92 duplicon, although substantially shorter, is closely related in sequence to the HLA-B and -C duplicons; (d) two of the six paralogous Alu J elements within the HLA-B and -C duplicons are associated with the HLA-X duplicon, confirming their evolutionary relationships within the beta block; and (e) the paralogous Alu J elements within the alpha block are distinctly different from those identified within the beta and kappa blocks. The sequence conservation and location of duplicated (paralogous) Alu J elements in the MHC class I region show that the beta and kappa blocks have evolved separately from the alpha block beginning at a time before or during the evolution of Alu J elements in primates. Received: 22 September 1999 / Accepted: 24 January 2000     

Duplication and Diversification of the Apolipoprotein CI (APOCI) Genomic Segment in Association with Retroelements

We have previously shown that several multicopy gene families within the major histocompatibility complex (MHC) arose from a process of segmental duplication. It has also been observed that retroelements play a role in generating diversity within these duplicated segments. The objective of this study was to compare the genomic organization of a gene duplication within another multicopy gene family outside the MHC. Using new continuous genomic sequence encompassing the APOE-CII gene cluster, we show that APOCI and its pseudogene, APOCI′, are contained within large duplicated segments which include sequences from the hepatic control region (HCR). Flanking Alu sequences are observed at both ends of the duplicated unit, suggesting a possible role in the integration of these segments. As observed previously within the MHC, the major differences between the segments are the insertion of sequences (approximately 200–1000 bp in length), consisting predominantly of Alu sequences. Ancestral retroelements also contribute to the generation of sequence diversity between the segments, especially within the 3′ poly(A) tract of Alu sequences. The exonic and regulatory sequences of the APOCI and HCR loci show limited sequence diversity, with exon 3 being an exception. Finally, the typing of pre- and postduplication Alus from both segments indicates an estimated time of duplication of approximately 37 million years ago (mya), some time prior to the separation of Old and New World monkeys. Received: 17 July 1999 / Accepted: 6 November 1999     

Genomic and phylogenetic analysis of the S100A7 (Psoriasin) gene duplications within the region of the S100 gene cluster on human chromosome 1q21

Abstract The human S100 gene family encodes the EF-hand superfamily of calcium-binding proteins, with at least 14 family members clustered relatively closely together on chromosome 1q21. We have analyzed the most recently available genomic sequence of the human S100 gene cluster for evidence of tandem gene duplications during primate evolutionary history. The sequences obtained from both GenBank and GoldenPath were analyzed in detail using various comparative sequence analysis tools. We found that of the S100A genes clustered relatively closely together within a genomic region of 260 kb, only the S100A7 (psoriasin) gene region showed evidence of recent duplications. The S100A7 gene duplicated region is composed of three distinct genomic regions, 33, 11, and 31 kb, respectively, that together harbor at least five identifiable S100A7-like genes. Regions 1 and 3 are in opposite orientation to each other, but each region carries two S100A7-like genes separated by an 11-kb intergenic region (region 2) that has only one S100A7-like gene, providing limited sequence resemblance to regions 1 and 3. The duplicated genomic regions 1 and 3 share a number of different retroelements including five Alu subfamily members that serve as molecular clocks. The shared (paralogous) Alu S insertions suggest that regions 1 and 3 were probably duplicated during or after the phase of AluS amplification some 30–40 mya. We used PCR to amplify an indel within intron 1 of the S100A7a and S100A7c genes that gave the same two expected product sizes using 40 human DNA samples and 1 chimpanzee sample, therefore confirming the presence of the region 1 and 3 duplication in these species. Comparative genomic analysis of the other S100 gene members shows no similarity between intergenic regions, suggesting that they diverged long before the emergence of the primates. This view was supported by the phylogenetic analysis of different human S100 proteins including the human S100A7 protein members. The S100A7 protein, also known as psoriasin, has important functions as a mediator and regulator in skin differentiation and disease (psoriasis), in breast cancer, and as a chemotactic factor for inflammatory cells. This is the first report of five copies of the S100A7 gene in the human genome, which may impact on our understanding of the possible dose effects of these genes in inflammation and normal skin development and pathogenesis.     

Coevolution of PERB11 (MIC) and HLA Class I Genes with HERV-16 and Retroelements by Extended Genomic Duplication

The recent availability of genomic sequence information for the class I region of the MHC has provided an opportunity to examine the genomic organization of HLA class I (HLAcI) and PERB11/MIC genes with a view to explaining their evolution from the perspective of extended genomic duplications rather than by simple gene duplications and/or gene conversion events. Analysis of genomic sequence from two regions of the MHC (the alpha- and beta-blocks) revealed that at least 6 PERB11 and 14 HLAcI genes, pseudogenes, and gene fragments are contained within extended duplicated segments. Each segment was searched for the presence of shared (paralogous) retroelements by RepeatMasker in order to use them as markers of evolution, genetic rearrangements, and evidence of segmental duplications. Shared Alu elements and other retroelements allowed the duplicated segments to be classified into five distinct groups (A to E) that could be further distilled down to an ancient preduplication segment containing a HLA and PERB11 gene, an endogenous retrovirus (HERV-16), and distinctive retroelements. The breakpoints within and between the different HLAcI segments were found mainly within the PERB11 and HLA genes, HERV-16, and other retroelements, suggesting that the latter have played a major role in duplication and indel events leading to the present organization of PERB11 and HLAcI genes. On the basis of the features contained within the segments, a coevolutionary model premised on tandem duplication of single and multipartite genomic segments is proposed. The model is used to explain the origins and genomic organization of retroelements, HERV-16, DNA transposons, PERB11, and HLAcI genes as distinct segmental combinations within the alpha- and beta-blocks of the human MHC. Received: 5 December 1998 / Accepted: 27 January 1999     

The temporal distribution of gene duplication events in a set of highly conserved human gene families

Using a data set of protein translations associated with map positions in the human genome, we identified 1520 mapped highly conserved gene families. By comparing sharing of families between genomic windows, we identified 92 potentially duplicated blocks in the human genome containing 422 duplicated members of these families. Using branching order in the phylogenetic trees, we timed gene duplication events in these families relative to the primate-rodent divergence, the amniote-amphibian divergence, and the deuterostome-protostome divergence. The results showed similar patterns of gene duplication times within duplicated blocks and outside duplicated blocks. Both within and outside duplicated blocks, numerous duplications were timed prior to the deuterostome-protostome divergence, whereas others occurred after the amniote-amphibian divergence. Thus, neither gene duplication in general nor duplication of genomic blocks could be attributed entirely to polyploidization early in vertebrate history. The strongest signal in the data was a tendency for intrachromosomal duplications to be more recent than interchromosomal duplications, consistent with a model whereby tandem duplication-whether of single genes or of genomic blocks-may be followed by eventual separation of duplicates due to chromosomal rearrangements. The rate of separation of tandemly duplicated gene pairs onto separated chromosomes in the human lineage was estimated at 1.7 x 10(-9) per gene-pair per year.     

The Evolution of MHC Diversity by Segmental Duplication and Transposition of Retroelements

Sequence analysis of a 237 kb genomic fragment from the central region of the MHC has revealed that the HLA-B and HLA-C genes are contained within duplicated segments peri-B (53 kb) and peri-C (48 kb), respectively, and separated by an intervening sequence (IF) of 30 kb. The peri-B and peri-C segments share at least 90% sequence homology except when interrupted by insertions/deletions including Alu, L1, an endogenous retrovirus, and pseudogenes. The sequences of peri-B, IF, and peri-C were searched for the presence of Alu elements to use as markers of evolution, chromosomal rearrangements, and polymorphism. Of 29 Alu elements, 14 were identified in peri-B, 11 in peri-C, and 4 in IF. The Alu elements in peri-B and peri-C clustered phylogenetically into two clades which were classified as ``preduplication' and ``postduplication' clades. Four Alu J elements that are shared by peri-B and peri-C and are flanked by homologous sequences in their paralogous locations, respectively, clustered into a ``preduplication' clade. By contrast, the majority of Alu elements, which are unique to either peri-B or peri-C, clustered into a postduplication clade together with the Alu consensus subfamily members ranging from platyrrhine-specific (Spqxcg) to catarrhine-specific Alu sequences (Y). The insertion of platyrrhine-specific Alu elements in postduplication locations of peri-B and peri-C implies that these two segments are the products of a duplication which occurred in primates prior to the divergence of the New World primate from the human lineage (35–44 mya). Examination of the paralogous Alu integration sites revealed that 9 of 14 postduplication Alu sequences have produced microsatellites of different length and sequence within the Alu 3′-poly A tail. The present analysis supports the hypothesis that HLA-B and HLA-C genes are products of an extended segmental duplication between 44 and 81 million years ago (mya), and that subsequent diversification of both genomic segments occurred because of the mobility and mutation of retroelements such as Alu repeats. Received: 21 May 1997 / Accepted: 9 July 1997     

Evolution of major histocompatibility complex by “en bloc” duplication before mammalian radiation

Duplications are an important mechanism for the emergence of genetic novelties. Reports on duplicated genes are numerous, and mechanisms for polyploidization or local gene duplication are beginning to be understood. When a local duplication is studied, searches are usually done gene-by-gene, and the size of duplicated segments is not often investigated. Therefore, we do not know if the gene in question has duplicated alone or with other genes, implying that "en bloc" duplications are poorly studied. We propose a method for identification of "en bloc" duplication using mapping, phylogenetic and statistical analyses. We show that two segments present in the major histocompatibility complex (MHC) region of human chromosome 6 have resulted from an "en bloc" duplication that took place between divergence of amniotes and methaterian/eutherian separation. These segments contain members of the same multigenic families, namely olfactory receptors genes, genes encoding proteins containing B30.2 domain, genes encoding proteins containing immunoglobulin V domain and MHC class I genes. We will discuss the fact that olfactory receptors and MHC genes have undergone positive selection, which could have helped in fixation of the surrounding genes.     

Primate evolution of the alpha-globin gene cluster and its Alu-like repeats

The arrangement of alpha-globin genes in Old World and New World monkeys and a prosimian, galago, has been determined by restriction mapping. Recombinant DNAs containing galago and Old World monkey alpha-globin genes have been isolated and subjected to a partial sequence determination for comparison to alpha-globin genes in human, chimpanzee and non-primate mammals. The results of this extensive structural analysis are relevant to several topics concerning the evolution of primate alpha-globin genes and Alu family repeats. All orders of higher primates (i.e. Old and New World monkeys, chimpanzee and human) have the same arrangement of alpha-globin genes. In contrast, the arrangement and correction of galago alpha-globin genes differ from those of higher primates, but are similar to those of non-primate mammals. The 5' and 3'-flanking regions of the human alpha 1 gene are orthologous to the corresponding region in galago, identifying the human alpha 2 gene as the more recently duplicated gene. The human psi alpha 1 gene is found to be inactivated after divergence of the human and galago lineages but prior to the divergence of human and monkey. Orthologous Alu family members in human and monkey DNAs indicate that the dispersion of some Alu repeats occurred prior to the divergence of these lineages. However, the Alu-like repeats of prosimian and higher primates result from entirely independent events giving rise to different repeat elements inserted at distinct genomic positions.     

Polymorphism and signature of selection in the MHC class I genes of the three-spined stickleback Gasterosteus aculeatus

The role and intensity of positive selection maintaining the polymorphism of major histocompatibility complex (MHC) class I genes in the three-spined stickleback Gasterosteus aculeatus was investigated. The highly polymorphic set of MHC class I genes found was organized in a single linkage group. Between 5 and 14 sequence variants per individual were identified by single-stranded conformation polymorphism (SSCP) analysis. Segregation analysis studied in 10 three-spined stickleback families followed the expected pattern of Mendelian inheritance. The gamete fusion in three-spined stickleback thus seems to be random with respect to the MHC class I genes. The DNA sequence analyses showed that the expressed MHC class I loci are under strong selection pressure, possibly mediated by parasites. Codons that were revealed to be under positive selection are potentially important in antigen binding. MHC class I sequences did not form significant supported clusters within a phylogenetic tree. Analogous to MHC class II genes, it was not possible to assign the class I sequences to a specific locus, suggesting that the class I genes may have been generated by recent gene duplication.     

Divergent origins and concerted expansion of two segmental duplications on chromosome 16

An unexpected finding of the human genome was the large fraction of the genome organized as blocks of interspersed duplicated sequence. We provide a comparative and phylogenetic analysis of a highly duplicated region of 16p12.2, which is composed of at least four different segmental duplications spanning in excess of 160 kb. We contrast the dispersal of two different segmental duplications (LCR16a and LCR16u). LCR16a, a 20 kb low-copy repeat sequence A from chromosome 16, was shown previously to contain a rapidly evolving novel hominoid gene family (morpheus) that had expanded within the last 10 million years of great ape/human evolution. We compare the dispersal of this genomic segment with a second adjacent duplication called LCR16u. The duplication contains a second putative gene family (KIAA0220/SMG1) that is represented approximately eight times within the human genome. A high degree of sequence identity (approximately 98%) was observed among the various copies of LCR16u. Comparative analyses with Old World monkey species show that LCR16a and LCR16u originated from two distinct ancestral loci. Within the human genome, at least 70% of the LCR16u copies were duplicated in concert with the LCR16a duplication. In contrast, only 30% of the chimpanzee loci show an association between LCR16a and LCR16u duplications. The data suggest that the two copies of genomic sequence were brought together during the chimpanzee/human divergence and were subsequently duplicated as a larger cassette specifically within the human lineage. The evolutionary history of these two chromosome-specific duplications supports a model of rapid expansion and evolutionary turnover among the genomes of man and the great apes.     

Phylogenetic tests of the hypothesis of block duplication of homologous genes on human chromosomes 6, 9, and 1

There are 10 gene families that have members on both human chromosome 6 (6p21.3, the location of the human major histocompatibility complex [MHC]) and human chromosome 9 (mostly 9q33-34). Six of these families also have members on mouse chromosome 17 (the mouse MHC chromosome) and mouse chromosome 2. In addition, four of these families have members on human chromosome 1 (1q21-25 and 1p13), and two of these have members on mouse chromosome 1. One hypothesis to explain these patterns is that members of the 10 gene families of human chromosomes 6 and 9 were duplicated simultaneously as a result of polyploidization or duplication of a chromosome segment ("block duplication"). A subsequent block duplication has been proposed to account for the presence of representatives of four of these families on human chromosome 1. Phylogenetic analyses of the 9 gene families for which data were available decisively rejected the hypothesis of block duplication as an overall explanation of these patterns. Three to five of the genes on human chromosomes 6 and 9 probably duplicated simultaneously early in vertebrate history, prior to the divergence of jawed and jawless vertebrates, and shortly after that, all four of the genes on chromosomes 1 and 9 probably duplicated as a block. However, the other genes duplicated at different times scattered over at least 1.6 billion years. Since the occurrence of these clusters of related genes cannot be explained by block duplication, one alternative explanation is that they cluster together because of shared functional characteristics relating to expression patterns.      

Polymorphic Alu insertions within the Major Histocompatibility Complex class I genomic region: a brief review

Most polymorphic Alu insertions (POALINs) belong to a subgroup of the Alu multicopy retrotransposon family of short interspersed nucleotide elements (SINEs) that are categorized as AluYb8 and AluYa5. The number of AluYb8/AluYa5 members (approximately 4,492 copies) is significantly less than the approximately one million fixed Alu copies per human genome. We have studied the presence of POALINs within the Major Histocompatibility Complex (MHC) class I region on the short arm of chromosome 6 (6p21.3) because this region has a high gene density, many genes with immune system functions, large sequence variations and diversity, duplications and redundancy, and a strong association with more than 100 different diseases. Since little is known about POALINs within the MHC genomic region, we undertook to identify some of the members of the AluYb8/AluYa5 subfamily and to study their frequency of distribution and genetic characteristics in different populations. As a result of our comparative genomic analyses, we identified the insertion sites for five POALINs distributed within the MHC class I region. This brief review outlines the locations of the insertions and sequence features of the five MHC POALINs, their single site and haplotype frequencies in different geographic populations, and their association with different HLA class I genes and disease. We show that the MHC POALINs have a potential value as lineage and linkage markers for the study of human population genetics, disease associations, genomic diversity and evolution.     

Genomic Organization Around the Centromeric End of the HLA Class I Region: Large-Scale Sequence Analysis

We previously sequenced two regions around the centromeric end of HLA class I and the boundary between class I and class III. In this paper we analyze the two regions of about 385 kb and confirm, giving a new line of evidence, that the following two pairs of the genomic segments were duplicated in evolution: (i) a 43-kb genomic segment including the HLA-B gene showing the highest polymorphism among the classical HLA class I loci (class Ia) and a 40-kb segment including the HLA-C locus showing the lowest polymorphism and (ii) a 52-kb segment including the MIC (MHC class I chain related gene) B and a 35-kb segment including MICA. We also found that repetitive elements such as SINEs, LINEs, and LTRs occupy as much as 47% of nucleotides in this 385-kb region. This unusually high content of repetitive elements indicates that repeat-mediated rearrangements have frequently occurred in the evolutionary history of the HLA class Ia region. Analysis of LINE compositions within the two pairs of duplicated segments revealed that (i) LINEs in these regions had been dispersed prior to both the duplication of the HLA-B and -C loci and the duplication of the MICB and MICA loci, and (ii) the divergence of the HLA-B and -C loci occurred prior to the duplication of the MICA and MICB loci. To find novel genes responsible for HLA class I-associated or other diseases, we performed computer analysis applying GenScan and GRAIL to GenBank's dbEST. As a result, at least five as yet uncharacterized genes were newly mapped on the HLA class I centromeric region studied. These novel genes should be analyzed further to determine their relationships to diseases associated with this region. Received: 16 June 1998 / Accepted: 18 August 1998     

Expansion and contraction of major histocompatibility complex genes: a teleostean example

The MHC is a multigene family that has arisen through recurrent expansion and contraction of genes, and a continuum of the evolutionary process is observed in the teleost fishes. The number of duplicated genes observed in different phylogenetic groups of teleost fish varies from one to 42, with only a few genes observed in the primitive euteleost species, and greater numbers of genes observed in the more advanced neoteleost species. In this study, an attempt is made to isolate all of the Mhc class I genes of an early neoteleost species, Atlantic cod (Gadus morhua L.), in the superorder Paracanthopterygii. Eighty-three sequences were isolated from the cDNA of an individual G. morhua. The level of gene duplication observed within each of the lineages and sublineages was similar, and most contained an estimated two to four duplicated genes. Mhc class I gene duplication in G. morhua was independent of, and possibly more recent than, extensive duplication in the Acanthopterygian superorder. Only limited contraction of Mhc genes is observed in G. morhua. A low level of haplotype diversity is observed, with most individuals containing at least one copy of each of the lineages tested. Divergence of the conserved N- and C-terminal residues of the antigen recognition site is observed, indicative of the initial stage of degeneration from classical to non-classical genes. However, most or all of the lineages are still polymorphic, and degeneration is present both within and among lineages. Thus, the outcome (i.e., which genes will remain classical) is as yet undetermined.     

Independent evolution of functional MHC class II DRB genes in New World bat species

Genes of the major histocompatibility complex (MHC) play a pivotal role in the vertebrate immune system and are attractive markers for functional, fitness-related, genetic variation. Although bats (Chiroptera) represent the second largest mammalian order and are prone to various emerging infectious diseases, little is known about MHC evolution in bats. In the present study, we examined expressed MHC class II DRB sequences (exons 1 to 4) of New World bat species, Saccopteryx bilineata, Carollia perspicillata, Noctilio albiventris and Noctilio leporinus (only exon 2). We found a wide range of copy number variation of DRB loci with one locus detected in the genus Noctilio and up to ten functional loci observed in S. bilineata. Sequence variation between alleles of the same taxa was high with evidence for positive selection. We found statistical support for recombination or gene conversion events among sequences within the same but not between bat species. Phylogenetic relationships among DRB alleles provided strong evidence for independent evolution of the functional MHC class II DRB genes in the three investigated species, either by recent gene duplication, or homogenization of duplicated loci by frequent gene conversion events. Phylogenetic analysis of all available chiropteran DRB exon 2 sequences confirmed their monophyletic origin within families, but revealed a possible trans-species mode of evolution pattern in congeneric bat species, e.g. within the genera Noctilio and Myotis. This is the first study investigating phylogenetic relationships of MHC genes within bats and therefore contributes to a better understanding of MHC evolution in one of the most dominant mammalian order.