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     

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.     

Sequence and evolution of cattle MHC class I cDNAs: concerted evolution has not taken place in cattle

To explore genetic mechanisms responsible for major histocompatibility complex (MHC) class I evolution in the artiodactyls, we cloned and sequenced MHC class I cDNAs from a Bos taurus bull heterozygous for cattle MHC (BoLA) class I serological specificities w2 and w30. Four unique cDNAs were found, indicating the presence of at least two MHC class I loci. Analysis of these four cDNAs and all previously published BoLA cDNA sequences suggested that there may be three cattle MHC class I loci. Additionally, comparison of all of the BoLA class I cDNAs to MHC class I cDNAs of other artiodactyls showed that some of the BoLA class I cDNAs were more similar to certain sheep cDNAs than they were to other cattle cDNAs. These data indicate that each BoLA class I locus has evolved independently after an ancestral gene duplication event and that inter-locus segmental exchange o or concerted evolution has not occurred rapidly enough to cause extensive divergence between the orthologous MHC class I loci of sheep and cattle.The nucleotide sequence data reported in this paper have been submitted to the GenBank nucleotide sequence database and have been assigned the accession numbers L02832–L02835. Correspondence to: T. L. Garber at the present address.     

Gene organization of human NOTCH4 and (CTG)n polymorphism in this human counterpart gene of mouse proto-oncogene Int3

The cDNA and genomic clones for the human counterpart of the mouse mammary tumor gene Int3 were isolated and sequenced. We designated this human major histocompatibility complex (MHC) class III gene as NOTCH4, since very recently, by sequencing cDNA clones, the complete form of the mouse proto-oncogene Int3 has been clarified and named Notch4. The present human NOTCH4 sequence is the first example of the genomic sequence for the extracellular portion of the mammalian Notch4, and by comparing it with the mouse Notch4 cDNA sequence, the exon/intron organization was clarified. The comparison of the predicted amino acid sequence of human NOTCH4 with those of other Notch homologues of a wide range of species revealed four subfamilies for mammalian Notch. In the protein coding region of human NOTCH4, we found (CTG)_n repeats showing a variable number tandem repeat (VNTR) polymorphism for different human leukocyte antigen (HLA) haplotypes. Ten genes mapped on 6p21.3, including NOTCH4, were found to have counterparts structurally and functionally similar to those mostly mapped on 9q33-q34, indicating segmental chromosome duplication during the course of evolution. Similarity of genes on chromosomes 1, 6, 9 and 19 was also discussed.     

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     

Toward understanding MHC disease associations: partial resequencing of 46 distinct HLA haplotypes

We carried out a resequencing project that examined 552 kb of sequence from each of 46 individual HLA haplotypes representing a diversity of HLA allele types, generating nearly 27 Mb of fully phased genomic sequence. Haplotype blocks were defined extending from telomeric of HLA-F to centromeric of HLA-DP including in total 5186 MHC SNPs. To investigate basic questions about the evolutionary origin of common HLA haplotypes, and to obtain an estimate of rare variation in the MHC, we similarly examined two additional sets of samples. In 19 independent HLA-A1, B8, DR3 chromosomes, the most common HLA haplotype in Northern European Caucasians, variation was found at 11 SNP positions in the 3600-kb region from HLA-A to DR. Partial resequencing of 282 individuals in the gene-dense class III region identified significant variability beyond what could have been detected by linkage to common SNPs.     

Interchromosomal duplication of major histocompatibility complex class I regions in rainbow trout (Oncorhynchus mykiss), a species with a presumably recent tetraploid ancestry

Salmonid fishes are among the few animal taxa with a probable recent tetraploid ancestor. The present study is the first to compare large (>100 kb) duplicated genomic sequence fragments in such species. Two contiguous stretches with major histocompatibility complex (MHC) class I genes were detected in a rainbow trout BAC library, mapped and sequenced. The MHC class I duplicated regions, mapped by fluorescence in situ hybridization (FISH), were shown to be located on different metaphase chromosomes, Chr 14 and 18. Gene organization in both duplications is similar to that in other fishes, in that the class I loci are tightly linked with the PSMB8, PSMB9, PSMB10 and ABCB3 genes. Whereas one region, Onmy-IA, has a classical MHC class I locus (UBA), Onmy-IB encodes only non-classical class Ib proteins. The nucleotide diversity between the Onmy-IA and Onmy-IB noncoding regions is about 14%. This suggests that the MHC class I duplication event has occurred about 60 mya close to the time of an hypothesized ancestral tetraploid event. The present article is the first convincing report on the co-existence of two closely related MHC class I core regions on two different chromosomes. The interchromosomal duplication and the homology levels are supportive of the tetraploid model.Nucleotide sequence data reported are available in the DDBJ/EMBL/GenBank database under the accession numbers AB162342, AB162343 and from AY525774 to AY525776.     

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.     

Structure and function of human transplantation antigens

Human transplantation antigens encoded in the major histocompatibility complex (MHC) region play a key role in regulating the immune responses. Here, we will describe the summary of our analyses on the structure and function of the human MHC molecules, HLA antigens as follows. 1) The genomic organization of the HLA antigen region was examined by cosmid cloning and pulsed-field gel electrophoresis technique. The HLA antigen region spans over at least 3,000 kb, and constitutes a multigene family. 2) Genetic polymorphisms in the HLA gene region were analyzed by Southern hybridization with restriction endonuclease digested genomic DNA using the class II cDNAs as probes (RFLP) and found to be tightly associated with each allo specificity. 3) The functional expression of the HLA class II gene product were observed after transfer of their cloned genes into the mouse fibroblast and human lymphocytes. 4) Narcolepsy is completely associated with HLA-DR2 Dw2, but no difference in the sequence of the DQ beta 1 domain could be found between narcoleptic and healthy individuals. This fact suggests that narcolepsy is not caused by mutation in the DQ beta gene. Based on results, it was inferred that one or both of the two Asps within the second variable region in the first domain of the DR beta chain is directly correlated with predisposition to narcolepsy.     

Genetic analysis of completely sequenced disease-associated MHC haplotypes identifies shuffling of segments in recent human history

The major histocompatibility complex (MHC) is recognised as one of the most important genetic regions in relation to common human disease. Advancement in identification of MHC genes that confer susceptibility to disease requires greater knowledge of sequence variation across the complex. Highly duplicated and polymorphic regions of the human genome such as the MHC are, however, somewhat refractory to some whole-genome analysis methods. To address this issue, we are employing a bacterial artificial chromosome (BAC) cloning strategy to sequence entire MHC haplotypes from consanguineous cell lines as part of the MHC Haplotype Project. Here we present 4.25 Mb of the human haplotype QBL (HLA-A26-B18-Cw5-DR3-DQ2) and compare it with the MHC reference haplotype and with a second haplotype, COX (HLA-A1-B8-Cw7-DR3-DQ2), that shares the same HLA-DRB1, -DQA1, and -DQB1 alleles. We have defined the complete gene, splice variant, and sequence variation contents of all three haplotypes, comprising over 259 annotated loci and over 20,000 single nucleotide polymorphisms (SNPs). Certain coding sequences vary significantly between different haplotypes, making them candidates for functional and disease-association studies. Analysis of the two DR3 haplotypes allowed delineation of the shared sequence between two HLA class II-related haplotypes differing in disease associations and the identification of at least one of the sites that mediated the original recombination event. The levels of variation across the MHC were similar to those seen for other HLA-disparate haplotypes, except for a 158-kb segment that contained the HLA-DRB1, -DQA1, and -DQB1 genes and showed very limited polymorphism compatible with identity-by-descent and relatively recent common ancestry (<3,400 generations). These results indicate that the differential disease associations of these two DR3 haplotypes are due to sequence variation outside this central 158-kb segment, and that shuffling of ancestral blocks via recombination is a potential mechanism whereby certain DR-DQ allelic combinations, which presumably have favoured immunological functions, can spread across haplotypes and populations.     

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.     

Segmental duplications and copy-number variation in the human genome

The human genome contains numerous blocks of highly homologous duplicated sequence. This higher-order architecture provides a substrate for recombination and recurrent chromosomal rearrangement associated with genomic disease. However, an assessment of the role of segmental duplications in normal variation has not yet been made. On the basis of the duplication architecture of the human genome, we defined a set of 130 potential rearrangement hotspots and constructed a targeted bacterial artificial chromosome (BAC) microarray (with 2,194 BACs) to assess copy-number variation in these regions by array comparative genomic hybridization. Using our segmental duplication BAC microarray, we screened a panel of 47 normal individuals, who represented populations from four continents, and we identified 119 regions of copy-number polymorphism (CNP), 73 of which were previously unreported. We observed an equal frequency of duplications and deletions, as well as a 4-fold enrichment of CNPs within hotspot regions, compared with control BACs (P < .000001), which suggests that segmental duplications are a major catalyst of large-scale variation in the human genome. Importantly, segmental duplications themselves were also significantly enriched >4-fold within regions of CNP. Almost without exception, CNPs were not confined to a single population, suggesting that these either are recurrent events, having occurred independently in multiple founders, or were present in early human populations. Our study demonstrates that segmental duplications define hotspots of chromosomal rearrangement, likely acting as mediators of normal variation as well as genomic disease, and it suggests that the consideration of genomic architecture can significantly improve the ascertainment of large-scale rearrangements. Our specialized segmental duplication BAC microarray and associated database of structural polymorphisms will provide an important resource for the future characterization of human genomic disorders.     

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 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.     

Comparative analysis of the bovine MHC class IIb sequence identifies inversion breakpoints and three unexpected genes

The bovine major histocompatibility complex (MHC) or BoLA is organized differently from typical mammalian MHCs in that a large portion of the class II region, called class IIb, has been transposed to a position near the centromere on bovine chromosome 23. Gene mapping indicated that the rearrangement resulted from a single inversion, but the boundaries and gene content of the inverted segment have not been fully determined. Here, we report the genomic sequence of BoLA IIb. Comparative sequence analysis with the human MHC revealed that the proximal inversion breakpoint occurred approximately 2.5 kb from the 3' end of the glutamate-cysteine ligase, catalytic subunit (GCLC) locus and that the distal breakpoint occurred about 2 kb from the 5' end from a divergent class IIDRbeta-like sequence designated DSB. Gene content, order and orientation of BoLA IIb are consistent with the single inversion hypothesis when compared with the corresponding region of the human class II MHC (HLA class II). Differences with HLA include the presence of a single histone H2B gene located between the proteasome subunit, beta type, 9 (PSMB9) and DMB loci and a duplicated TAP2 with a variant splice site. BoLA IIb spans approximately 450 kb DNA, with 20 apparently intact genes and no obvious pseudogenes. The region contains 227 simple sequence repeats (SSRs) and approximately 167 kb of retroviral-related repetitive DNA. Nineteen of the 20 genes identified in silico are supported by bovine EST data indicating that the functional gene content of BoLA IIb has not been diminished because it has been transposed from the remainder of BoLA genes.     

Conserved haplotype blocks within the sheep MHC and low SNP heterozygosity in the Class IIa subregion

This report describes single-nucleotide polymorphisms (SNPs) in the sheep major histocompatibility complex (MHC) class II and class III regions and provides insights into the internal structure of this important genomic complex. MHC haplotypes were deduced from sheep family trios based on genotypes from 20 novel SNPs representative of the class II region and 10 previously described SNPs spanning the class III region. All 30 SNPs exhibited Hardy-Weinberg proportions in the sheep population studied. Recombination within an extended sire haplotype was observed within the class II region for 4 of 20 sheep chromosomes, thereby supporting the presence of separated IIa and IIb subregions similar to those present in cattle. SNP heterozygosity varied across the class II and III regions. One segment of the class IIa subregion manifested very low heterozygosity for several SNPs spanning approximately 120 Kbp. This feature corresponds to a subregion within the human MHC class II region previously described as a 'SNP desert' because of its paucity of SNPs. Linkage disequilibrium (LD) was reduced at the junction separating the putative class IIb and IIa subregions and also between the class IIa and the class III subregions. The latter observation is consistent with either an unmapped physical separation at this location or more likely a boundary characterized by more frequent recombination between two conserved subregions, each manifesting high within-block LD. These results identify internal blocks of loci in the sheep MHC, within which recombination is relatively rare.     

Major histocompatibility complex class I loci from the gray whale (Eschrichtius robustus)

Sequences from exons encoding the peptide binding region of MHC class I (MHC-I) molecules were isolated from California gray whale (Eschrichtius robustus) genomic DNA to initiate an investigation of variation in these genes in a cetacean. These represent the first mysticete MHC-I sequences to be reported. The analysis of gray whale MHC-I sequences suggests the presence of at least three loci, which share greatest similarity to MHC-I in the ungulates, consistent with current views on cetacean phylogenetics. The peptide binding region of MHC is the most polymorphic part of the molecule and analysis of the variation and synonymous to nonsynonymous substitution ratios in gray whale sequences found these genes to display polymorphism characteristics similar to that attributed to selection in other species.