Gene organisation, sequence variation and isochore structure at the centromeric boundary of the human MHC

Mapping genomic loci implicates genes and synaptic biology in schizophrenia

Schizophrenia has a heritability of 60-80%1, much of which is attributable to common risk alleles. Here, in a two-stage genome-wide association study of up to 76,755 individuals with schizophrenia and 243,649 control individuals, we report common variant associations at 287 distinct genomic loci. Associations were concentrated in genes that are expressed in excitatory and inhibitory neurons of the central nervous system, but not in other tissues or cell types. Using fine-mapping and functional genomic data, we identify 120 genes (106 protein-coding) that are likely to underpin associations at some of these loci, including 16 genes with credible causal non-synonymous or untranslated region variation. We also implicate fundamental processes related to neuronal function, including synaptic organization, differentiation and transmission. Fine-mapped candidates were enriched for genes associated with rare disruptive coding variants in people with schizophrenia, including the glutamate receptor subunit GRIN2A and transcription factor SP4, and were also enriched for genes implicated by such variants in neurodevelopmental disorders. We identify biological processes relevant to schizophrenia pathophysiology; show convergence of common and rare variant associations in schizophrenia and neurodevelopmental disorders; and provide a resource of prioritized genes and variants to advance mechanistic studies.

Origin and evolution of the major histocompatibility complex class I region in eutherian mammals

Major histocompatibility complex (MHC) genes in vertebrates are vital in defending against pathogenic infections. To gain new insights into the evolution of MHC Class I (MHCI) genes and test competing hypotheses on the origin of the MHCI region in eutherian mammals, we studied available genome assemblies of nine species in Afrotheria, Xenarthra, and Laurasiatheria, and successfully characterized the MHCI region in six species. The following numbers of putatively functional genes were detected: in the elephant, four, one, and eight in the extended class I region, and κ and β duplication blocks, respectively; in the tenrec, one in the κ duplication block; and in the four bat species, one or two in the β duplication block. Our results indicate that MHCI genes in the κ and β duplication blocks may have originated in the common ancestor of eutherian mammals. In the elephant, tenrec, and all four bats, some MHCI genes occurred outside the MHCI region, suggesting that eutherians may have a more complex MHCI genomic organization than previously thought. Bat-specific three- or five-amino-acid insertions were detected in the MHCI α1 domain in all four bats studied, suggesting that pathogen defense in bats relies on MHCIs having a wider peptide-binding groove, as previously assayed by a bat MHCI gene with a three-amino-acid insertion showing a larger peptide repertoire than in other mammals. Our study adds to knowledge on the diversity of eutherian MHCI genes, which may have been shaped in a taxon-specific manner.

A genome-wide association study on chronic HBV infection and its clinical progression in male Han-Taiwanese

It is common to observe the clustering of chronic hepatitis B surface antigen (HBsAg) carriers in families. Intra-familial transmission of hepatitis B virus (HBV) could be the reason for the familial clustering of HBsAg carriers. Additionally, genetic and gender factors have been reported to be involved. We conducted a three-stage genome-wide association study to identify genetic factors associated with chronic HBV susceptibility. A total of 1,065 male controls and 1,623 male HBsAg carriers were included. The whole-genome genotyping was done on Illumina HumanHap550 beadchips in 304 healthy controls and HumanHap610 beadchips in 321 cases. We found that rs9277535 (HLA-DPB1, P = 4.87×10(-14)), rs9276370 (HLA-DQA2, P = 1.9×10(-12)), rs7756516 and rs7453920 (HLA-DQB2, P = 1.48×10(-11) and P = 6.66×10(-15) respectively) were significantly associated with persistent HBV infection. A novel SNP rs9366816 near HLA-DPA3 also showed significant association (P = 2.58×10(-10)). The "T-T-G-G-T" haplotype of the five SNPs further signified their association with the disease (P = 1.48×10(-12); OR = 1.49). The "T-T" haplotype composed of rs7756516 and rs9276370 was more prevalent in severe disease subgroups and associated with non-sustained therapeutic response (P = 0.0262). The "G-C" haplotype was associated with sustained therapeutic response (P = 0.0132; OR = 2.49). We confirmed that HLA-DPB1, HLA-DQA2 and HLA-DQB2 loci were associated with persistent HBV infection in male Taiwan Han-Chinese. In addition, the HLA-DQA2 and -DQB2 complex was associated with clinical progression and therapeutic response.

The human Major Histocompatibility Complex as a paradigm in genomics research

Since its discovery more than 50 years ago, the human Major Histocompatibility Complex (MHC) on chromosome 6p21.3 has been at the forefront of human genetic research. Here, we review from a historical perspective the major advances in our understanding of the nature and consequences of genetic variation which have involved the MHC, as well as highlighting likely future directions. As a consequence of its particular genomic structure, its remarkable polymorphism and its early implication in numerous diseases, the MHC has been considered as a model region for genomics, being the first substantial region to be sequenced and establishing fundamental concepts of linkage disequilibrium, haplotypic structure and meiotic recombination. Recently, the MHC became the first genomic region to be entirely re-sequenced for common haplotypes, while studies mapping gene expression phenotypes across the genome have strongly implicated variation in the MHC. This review shows how the MHC continues to provide new insights and remains in the vanguard of contemporary research in human genomics.

Characterization of a Novel Gene in the Extended MHC Region of Mouse, NG29/Cd320, a Homolog of the Human CD320

Background

The MHC region of the chromosome contains a lot of genes involved in immune responses. Here we have investigated the mouse NG29/Cd320 gene in the centrometrically extended MHC region of chromosome 17.

Methods

We cloned the NG29 gene by RT-PCR and confirmed the tissue distribution of its gene expression by northern blot hybridization. We generated the NG29 gene expression constructs and polyclonal antibody against the NG29 protein to perform the immunofluorescence, immunoprecipitation and flow cytometric analysis.

Results

The murine NG29 gene and its human homologue, the CD320/8D6 gene, were similar in the gene structure and tissue expression patterns. We cloned the NG29 gene and confirmed its expression in plasma membrane and intracellular compartments by transfecting its expresssion constructs into HEK 293T cells. The immunoprecipitation studies with rabbit polyclonal antibody raised against the NG29-NusA fusion protein indicated that NG29 protein was a glycoprotein of about 45 kDa size. A flow cytometric analysis also showed the NG29 expression on the surface of Raw 264.7 macrophage cell line.

Conclusion

These findings suggested that NG29 gene in mouse extended MHC class II region was the orthologue of human CD320 gene even though human CD320/8D6 gene was located in non-MHC region, chromosome 19p13.

Sequence variation within the major histocompatibility complex subregion centromeric of HLA class II in type 1 diabetes

The extended major histocompatibility complex (xMHC) has been studied intensively with regard to type 1 diabetes (T1D) predisposition. So far, little attention has been given to the subregion centromeric of MHC class II. We selected five single nucleotide polymorphisms in genes with potential immune-related functions in the genomic regions of death-domain-associated protein 6 (DAXX, apoptosis associated), TAP-binding protein (TAPBP, human leukocyte antigen class I loading) and retinoic acid receptor beta (RXRB, vitamin D receptor function) that may bear relevance to the pathogenesis of T1D. A total of 277 unrelated individuals with juvenile-onset T1D and 286 control subjects were genotyped using sequence-specific priming-polymerase chain reaction. The genotype and allelic frequencies of the markers tested were not significantly different between patients and control subjects. Subsequent haplotype analysis showed six DAXX-TAPBP-RXRB haplotypic configurations. No difference was observed between patients and control cohorts when stratified for T1D high-risk DQ2-DR17 and DQ8-DR4 haplotypes. However, the distribution of these haplotypes affected T1D susceptibility encoded by the intermediate risk haplotypes DQ5-DR1 and DQ2-DR7 by increasing and decreasing susceptibility, respectively. We propose that studying genetic variants in the xMHC may be particularly rewarding to define disease pathways in patients displaying intermediate risk DQ-DR haplotypes.

Full-length sequence analysis of the HLA-DRB1 locus suggests a recent origin of alleles

The HLA region harbors some of the most polymorphic loci in the human genome. Among them is the class II locus HLA-DRB1, with more than 400 known alleles. The age of the polymorphism and the rate at which new alleles are generated at HLA loci has caused much controversy over the years. Previous studies have mostly been restricted to the 270 base pairs that constitute the second exon and represent the most variable part of the gene. Here, we investigate the evolutionary history of the HLA-DRB1 locus on the basis of an analysis of 15 genomic full-length alleles (10-15 kb). In addition, the variation in 49 complete coding sequences and 322 exon 2 sequences were analyzed. When excluding exon 2 from the analysis, the diversity at the synonymous sites was found to be similar to the intron diversity. The overall diversity in noncoding region was also similar to the genome average. The DRB1*03 lineage has been found in human, chimpanzee, bonobo, gorilla, and orangutan. An ancestral "proto HLA-DRB1*03 lineage" appeared to have diverged in the last 5 million years into the human-specific lineages *08, *11, *13, and *14. With exception to exon 2, both the coding- and the noncoding diversity suggests a recent origin (<1 million years ago) for most of the alleles at the HLA-DRB1 locus. Sites encoding for amino acids involved in antigen binding [antigen recognizing sites (ARS)] appear to have a more ancient origin. Taken together, the recent origin of most alleles, the high diversity between allelic lineages, and the ancient origin of sequence motifs in exon 2, is consistent with a relatively rapid generation of novel alleles by gene conversion like events.

Evolutionary history of the ABCB2 genomic region in teleosts

Gene duplication, silencing and translocation have all been implicated in shaping the unique genomic architecture of the teleost MH regions. Previously, we demonstrated that trout possess five unlinked regions encoding MH genes. One of these regions harbors ABCB2 which in all other vertebrate classes is found in the MHC class II region. In this study, we sequenced a BAC contig for the trout ABCB2 region. Analysis of this region revealed the presence of genes homologous to those located in the human class II (ABCB2, BRD2, psiDAA), extended class II (RGL2, PHF1, SYGP1) and class III (PBX2, Notch-L) regions. The organization and syntenic relationships of this region were then compared to similar regions in humans, Tetraodon and zebrafish to learn more about the evolutionary history of this region. Our analysis indicates that this region was generated during the teleost-specific duplication event while also providing insight about potential MH paralogous regions in teleosts.

Comparative genomics of major histocompatibility complexes

The major histocompatibility complex (MHC) is a gene dense region found in all jawed vertebrates examined to date. The MHC contains a high percentage of immune genes, in particular genes involved in antigen presentation, which are generally highly polymorphic. The region plays an important role in disease resistance. The clustering of MHC genes could be advantageous for co-evolution or regulation, and its study in many species is desirable. Even though some linkage of MHC genes is apparent in all gnathostomes, the genomic organization can differ greatly by species, suggesting rapid evolution of MHC genes after divergence from a common ancestor. Previous reviews of comparative MHC organization have been written when relatively fragmentary sequence and mapping data were available on many species. This review compares maps of MHC gene orders in commonly studied species, where extensive sequencing has been performed.

Gene map of the extended human MHC

The major histocompatibility complex (MHC) is the most important region in the vertebrate genome with respect to infection and autoimmunity, and is crucial in adaptive and innate immunity. Decades of biomedical research have revealed many MHC genes that are duplicated, polymorphic and associated with more diseases than any other region of the human genome. The recent completion of several large-scale studies offers the opportunity to assimilate the latest data into an integrated gene map of the extended human MHC. Here, we present this map and review its content in relation to paralogy, polymorphism, immune function and disease.

Loss of heterozygosity, a frequent but a non-exclusive mechanism responsible for HLA dysregulation in non-Hodgkin's lymphomas

The frequent alteration of human leucocyte antigen (HLA) class I molecule expression observed in non-Hodgkin's lymphomas (NHL), similarly to solid tumours, has been reported to favour tumoral escape from the immune system. In order to identify the underlying mechanisms, we analysed 15 HLA defective NHL including partial (n = 10) and total class I (n = 5) loss, as well as HLA class II defects (n = 5). The HLA defect concerned HLA-A and -B antigens in 14 of 15 cases. In the cases with partial defect, the use of specific allelic monoclonal antibodies detected a defect of both alleles of A or B loci in six of seven tested cases. Allelic reverse transcription polymerase chain reaction (RT-PCR) demonstrated defects in six of nine cases, including four alterations of both A and B mRNA alleles. Real-time quantitative RT-PCR (RQ-PCR) did not detect the HLA-DR transcript in the two negative HLA-DR lymphomas, contrasting with the presence of CMH II transactivator (CIITA) transcript. Loss of heterozygosity (LOH) was detected in nine of 14 cases through variable pattern of nine microsatellites markers of the HLA locus. Taken together, these findings demonstrate the complexity and the variability of the mechanisms underlying HLA protein deficiencies with a high frequency of LOH. The diversity of these mechanisms indicates the importance of positive selection of HLA altered clones in the development of these NHL cases.

Representing GC variation along eukaryotic chromosomes

Genome sequencing now permits direct visual representation, at any scale, of GC heterogeneity along the chromosomes of several higher eukaryotes. Plots can be easily obtained from the chromosomal sequences, yet sequence releases of mammalian or plant chromosomes still tend to use small scales or window sizes that obscure important large-scale compositional features. To faithfully reveal, at one glance, the compositional variation at a given scale, we have devised a simple scheme that combines line plots with color-coded shading of the regions underneath the plots. The scheme can be applied to different eukaryotic genomes to facilitate their comparison, as illustrated here for a sample of chromosomes chosen from seven selected species. As a complement to a previously published compact view of isochores in the human genome sequence, we include here an analogous map for the recently sequenced mouse genome, and discuss the contribution of repetitive DNA to the GC variation along the plots. Supplementary information, including a database of color-coded GC profiles for all recently sequenced eukaryotes and the program draw_chromosomes_gc.pl used to obtain them, are available at.

IsoFinder: computational prediction of isochores in genome sequences

Isochores are long genome segments homogeneous in G+C. Here, we describe an algorithm (IsoFinder) running on the web (http://bioinfo2.ugr.es/IsoF/isofinder.html) able to predict isochores at the sequence level. We move a sliding pointer from left to right along the DNA sequence. At each position of the pointer, we compute the mean G+C values to the left and to the right of the pointer. We then determine the position of the pointer for which the difference between left and right mean values (as measured by the t-statistic) reaches its maximum. Next, we determine the statistical significance of this potential cutting point, after filtering out short-scale heterogeneities below 3 kb by applying a coarse-graining technique. Finally, the program checks whether this significance exceeds a probability threshold. If so, the sequence is cut at this point into two subsequences; otherwise, the sequence remains undivided. The procedure continues recursively for each of the two resulting subsequences created by each cut. This leads to the decomposition of a chromosome sequence into long homogeneous genome regions (LHGRs) with well-defined mean G+C contents, each significantly different from the G+C contents of the adjacent LHGRs. Most LHGRs can be identified with Bernardi's isochores, given their correlation with biological features such as gene density, SINE and LINE (short, long interspersed repetitive elements) densities, recombination rate or single nucleotide polymorphism variability. The resulting isochore maps are available at our web site (http://bioinfo2.ugr.es/isochores/), and also at the UCSC Genome Browser (http://genome.cse.ucsc.edu/).

Identification of isochore boundaries in the human genome using the technique of wavelet multiresolution analysis

Incorporated with the Z curve method, the technique of wavelet multiresolution (also known as multiscale) analysis has been proposed to identify the boundaries of isochores in the human genome. The human MHC sequence and the longest contigs of human chromosomes 21 and 22 are used as examples. The boundary between the isochores of Class III and Class II in the MHC sequence has been detected and found to be situated at the position 2,490,368bp. This result is in good agreement with the experimental evidence. An isochore with a length of about 7Mb in chromosome 21 has been identified and found to be gene- and Alu-poor. We have also found that the G+C content of chromosome 21 is more homogeneous than that of chromosome 22. Compared with the window-based methods, the present method has the highest resolution for identifying the boundaries of isochores, even at a scale of single base. Compared with the entropic segmentation method, the present method has the merits of more intuitiveness and less calculations. The important conclusion drawn in this study is that the segmentation points, at which the G+C content undergoes relatively dramatic changes, do exist in the human genome. These 'singularity' points may be considered to be candidates of isochore boundaries in the human genome. The method presented is a general one and can be used to analyze any other genomes.

An isochore map of the human genome based on the Z curve method

The distribution of the G+C content in the human genome has been studied by using a windowless technique derived from the Z curve method. The most important findings presented in this paper are twofold. First, abrupt variations of the G+C content along human chromosome sequences are the main variation patterns of G+C content. It is found that at some sites, the G+C content undergoes abrupt changes from a G+C-rich region to a G+C-poor region alternatively and vice versa. Second, it is shown that long domains with relatively homogeneous G+C content along each chromosome do exist. These domains are thought to be isochores, which usually have sharp boundaries. Consequently, 56 isochores longer than 3 Mb have been identified in chromosomes 1-22, X and Y. Boundaries, size and G+C content of each isochore identified are listed in detail. As an example to demonstrate the power of the method, the boundary between the Classes III and II isochores of the MHC sequence has been determined and found to be at 2,477,936, which is in good agreement with the experimental evidence. A homogeneity index is introduced to measure the homogeneity of G+C content in isochores. We emphasize that the homogeneity of G+C content is relative. The isochores in which the G+C content keeps absolutely constant do not exist. Isochore structures appear to be a basic organization of the human genome. Due to the relevance to many important biological functions, the clarification of isochore structures will provide much insight into the understanding of the human genome.

No evidence of type 1 diabetes susceptibility genes in the region centromeric of the HLA complex

There is strong evidence that DQB1, DQA1, and DRB1 alleles are not the only contributors to the human leukocyte antigen (HLA) linked type 1 diabetes (T1D) predisposition. Although the HLA complex is much studied for disease association, little is known about the neighboring centromeric region. We have previously found suggestive association on DQ2-DR3 haplotypes for marker D6S291, located 3.6-Mb centromeric of HLA-DQB1. This region on human chromosome 6 is syntenic to a part of the region adjacent to the mouse major histocompatibility complex (MHC) on chromosome 17, which has been suggested to harbor a susceptibility gene in mouse (Idd16). To evaluate a possible role of the region centromeric of HLA-DQB1 in human T1D, we have scanned the region with nine microsatellite markers in 267 T1D families from five different populations. Our results indicate that the characteristic strong linkage disequilibrium in the HLA complex does not extend into this region. Furthermore, we did not detect any consistent T1D association for the markers analyzed in the study. In conclusion, our data argue against the presence of any strong genetic susceptibility factors for T1D in the region centromeric of the HLA complex.

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.

Comparative genome organization of human, murine, and feline MHC class II region

To study comparative molecular dynamics in the genesis of the major histocompatibility complex (MHC), we determined a complete nucleotide sequence spanning 758,291 bp of the domestic cat (Felis catus) extended and classical class II region. The feline class II MHC includes 44 genes (31 predicted to be expressed) which display DNA sequence homology and ordered gene synteny with human HLA and mouse H2, in extended class II and centromere proximal regions (DM to DO) of the classical class II region. However, remarkable genomic alterations including gene gain and loss plus size differentials of 250 kb are evident in comparisons of the cat class II with those of human and mouse. The cat MHC lacks the entire DQ region and retains only relict pseudogene homologs of DP genes, compensated by expansion and reorganization of seven modern DR genes. Repetitive gene families within the feline MHC comprise 35% of the feline MHC with very different density and abundance of GC levels, SINES, LINES, STRs, and retro-elements from the same repeats in human and mouse MHC. Comparison of the feline MHC with the murine and human MHC offers a detailed view of the consequences of genome organization in three mammalian lineages.

Isochores merit the prefix 'iso'

The isochore concept in the human genome sequence was challenged in an analysis by the International Human Genome Sequencing Consortium (IHGSC). We argue here that a statement in the IHGSC's analysis concerning the existence of isochores is misleading, because the homogeneity was not examined at a large enough length scale and consequently an inappropriate statistical test was applied. A test of the existence of isochores should be equivalent to a test of homogeneity or equality of windowed GC%. The statistical test applied in the IHGSC's analysis, the binomial test, is a test of whether individual bases are independent and identically-distributed (iid). For testing the existence of isochores, or homogeneity in windowed GC%, we propose to use another statistical test: the analysis of variance (ANOVA). It can be shown that DNA sequences that are rejected by the binomial test may not be rejected by the ANOVA test.

Comparative and evolutionary analysis of the rhesus macaque extended MHC class II region

The sequence-based map of a part of the rhesus macaque major histocompatibility complex (MHC) extended class II region is presented. The sequenced region encompasses 67,401 bp and contains the SACM2L, RING1, FABGL and KE4 genes, as well as the HTATSF1-like and ZNF-like pseudogenes. Similar to human, but different from rat and mouse, no class I genes are found in the SACM2L- RING1 interval. The rhesus macaque extended MHC class II region shows a high degree of conservation of exonic as well as intronic and intergenic sequences compared with the respective human region. It is concluded that this particular genomic organization of the extended class II region-i.e., the absence of class I genes and the presence of the HTATSF1-like and ZNF-like pseudogenes-can be traced back to a common ancestor of humans and rhesus macaques about 23 million years ago.

Genetics of the immune response: identifying immune variation within the MHC and throughout the genome

With the advent of modern genomic sequencing technology the ability to obtain new sequence data and to acquire allelic polymorphism data from a broad range of samples has become routine. In this regard, our investigations have started with the most polymorphic of genetic regions fundamental to the immune response in the major histocompatibility complex (MHC). Starting with the completed human MHC genomic sequence, we have developed a resource of methods and information that provide ready access to a large portion of human and nonhuman primate MHCs. This resource consists of a set of primer pairs or amplicons that can be used to isolate about 15% of the 4.0 Mb MHC. Essentially similar studies are now being carried out on a set of immune response loci to broaden the usefulness of the data and tools developed. A panel of 100 genes involved in the immune response have been targeted for single nucleotide polymorphism (SNP) discovery efforts that will analyze 120 Mb of sequence data for the presence of immune-related SNPs. The SNP data provided from the MHC and from the immune response panel has been adapted for use in studies of evolution, MHC disease associations, and clinical transplantation.

Are isochore sequences homogeneous?

Three statistical/mathematical analyses are carried out on isochore sequences: spectral analysis, analysis of variance, and segmentation analysis. Spectral analysis shows that there are GC content fluctuations at different length scales in isochore sequences. The analysis of variance shows that the null hypothesis (the mean value of a group of GC contents remains the same along the sequence) may or may not be rejected for an isochore sequence, depending on the subwindow sizes at which GC contents are sampled, and the window size within which group members are defined. The segmentation analysis shows that there are stronger indications of GC content changes at isochore borders than within an isochore. These analyses support the notion of isochore sequences, but reject the assumption that isochore sequences are homogeneous at the base level. An isochore sequence may pass a homogeneity test when GC content fluctuations at smaller length scales are ignored or averaged out.

Isochore chromosome maps of the human genome

The human genome is a mosaic of isochores, which are long DNA segments (z.Gt;300 kbp) relatively homogeneous in G+C. Human isochores were first identified by density-gradient ultracentrifugation of bulk DNA, and differ in important features, e.g. genes are found predominantly in the GC-richest isochores. Here, we use a reliable segmentation method to partition the longest contigs in the human genome draft sequence into long homogeneous genome regions (LHGRs), thereby revealing the isochore structure of the human genome. The advantages of the isochore maps presented here are: (1) sequence heterogeneities at different scales are shown in the same plot; (2) pair-wise compositional differences between adjacent regions are all statistically significant; (3) isochore boundaries are accurately defined to single base pair resolution; and (4) both gradual and abrupt isochore boundaries are simultaneously revealed. Taking advantage of the wide sample of genome sequence analyzed, we investigate the correspondence between LHGRs and true human isochores revealed through DNA centrifugation. LHGRs show many of the typical isochore features, mainly size distribution, G+C range, and proportions of the isochore classes. The relative density of genes, Alu and long interspersed nuclear element repeats and the different types of single nucleotide polymorphisms on LHGRs also coincide with expectations in true isochores. Potential applications of isochore maps range from the improvement of gene-finding algorithms to the prediction of linkage disequilibrium levels in association studies between marker genes and complex traits. The coordinates for the LHGRs identified in all the contigs longer than 2 Mb in the human genome sequence are available at the online resource on isochore mapping: http://bioinfo2.ugr.es/isochores.

Identification of BING-4 cancer antigen translated from an alternative open reading frame of a gene in the extended MHC class II region using lymphocytes from a patient with a durable complete regression following immunotherapy

Multiple human cancer Ags have been identified, although little is known concerning which would be most effectively used in cancer immunotherapy. To gain insight into the selection of appropriate Ags, the immunologic reactivity of a patient who had a durable complete regression of melanoma metastases was measured. PBMCs were directly cloned using the monoclonal anti-CD3 Ab OKT3 and IL-2 without any bias introduced by previous culture. A lymphocyte clone recognized a previously unknown shared melanoma Ag that was identified as the BING-4 protein encoded in a gene-rich region of the extended class II MHC. The HLA-A2-restricted BING-4 immunodominant peptide was translated from a 10-aa-long alternative open reading frame. In vitro sensitization against this peptide generated lymphocytes reactive against HLA-A2(+) melanomas. Real-time semiquantitative RT-PCR analysis revealed that 8 of 15 melanoma cell lines overexpressed BING-4, and this correlated with recognition by lymphocytes. Overexpression was not found in normal tissues or other tumor types. Thus, BING-4 represents another candidate Ag for possible use in the immunotherapy of patients with melanoma.

The human major histocompatability complex: lessons from the DNA sequence

The entire 3.6-MbpDNA sequence of a human major histocompatibility complex derived from a composite of DNA clones from different haplotypes, was completed in 1999, primarily through the work of four main groups. At that time, it was the longest contiguous human DNA sequence to have been determined. The sequence is of extremely high quality and accuracy. In this review, we discuss how the DNA sequence has facilitated our understanding of the biology and genetics of the major histocompatibility complex. We suggest some ways in which the sequence may be exploited in the future to explore the relationship between the extraordinary polymorphism of the region and its association with both autoimmune and infectious diseases.

Isochore chromosome maps of eukaryotic genomes

Analytical DNA ultracentrifugation revealed that eukaryotic genomes are mosaics of isochores: long DNA segments (>>300 kb on average) relatively homogeneous in G+C. Important genome features are dependent on this isochore structure, e.g. genes are found predominantly in the GC-richest isochore classes. However, no reliable method is available to rigorously partition the genome sequence into relatively homogeneous regions of different composition, thereby revealing the isochore structure of chromosomes at the sequence level. Homogeneous regions are currently ascertained by plain statistics on moving windows of arbitrary length, or simply by eye on G+C plots. On the contrary, the entropic segmentation method is able to divide a DNA sequence into relatively homogeneous, statistically significant domains. An early version of this algorithm only produced domains having an average length far below the typical isochore size. Here we show that an improved segmentation method, specifically intended to determine the most statistically significant partition of the sequence at each scale, is able to identify the boundaries between long homogeneous genome regions displaying the typical features of isochores. The algorithm precisely locates classes II and III of the human major histocompatibility complex region, two well-characterized isochores at the sequence level, the boundary between them being the first isochore boundary experimentally characterized at the sequence level. The analysis is then extended to a collection of human large contigs. The relatively homogeneous regions we find show many of the features (G+C range, relative proportion of isochore classes, size distribution, and relationship with gene density) of the isochores identified through DNA centrifugation. Isochore chromosome maps, with many potential applications in genomics, are then drawn for all the completely sequenced eukaryotic genomes available.

Delineating relative homogeneous G+C domains in DNA sequences

The concept of homogeneity of G+C content is always relative and subjective. This point is emphasized and quantified in this paper using a simple example of one sequence segmented into two subsequences. Whether the sequence is homogeneous or not can be answered by whether the two-subsequence model describes the DNA sequence better than the one-sequence model. There are at least three equivalent ways of looking at the 1-to-2 segmentation: Jensen-Shannon divergence measure, log likelihood ratio test, and model selection using Bayesian information criterion. Once a criterion is chosen, a DNA sequence can be recursively segmented into multiple domains. We use one subjective criterion called segmentation strength based on the Bayesian information criterion. Whether or not a sequence is homogeneous and how many domains it has depend on this criterion. We compare six different genome sequences (yeast S. cerevisiae chromosome III and IV, bacterium M. pneumoniae, human major histocompatibility complex sequence, longest contigs in human chromosome 21 and 22) by recursive segmentations at different strength criteria. Results by recursive segmentation confirm that yeast chromosome IV is more homogeneous than yeast chromosome III, human chromosome 21 is more homogeneous than human chromosome 22, and bacterial genomes may not be homogeneous due to short segments with distinct base compositions. The recursive segmentation also provides a quantitative criterion for identifying isochores in human sequences. Some features of our recursive segmentation, such as the possibility of delineating domain borders accurately, are superior to those of the moving-window approach commonly used in such analyses.

DPB1*8601, a previously unrecognized DPB1 variant in the Caucasoid population

A new, previously unrecognized DPB1* allele, DPB1*8601, was found in a Swedish family. The new allele was carried on the common North European haplotype HLA A1-B8-DR3. Both individuals carrying the new allele were initially typed as clear DPB1*4601,*6601 but after family studies and further typing with allele-specific primers it was concluded that a new allele was present together with the common DPB1*0401. The new allele was investigated by direct sequencing of exon 2 in both forward and reverse directions employing intron primers combined with either an allele-specific sense or anti-sense biotinylated primer for bi-directional sequencing. The new allele is identical to DPB1*1701 in the five first variable regions. In the sixth region, however, DPB1*8601 carries the GGPM motif shared by several common alleles such as DPB1*0201 and 0401and 0402.

Characterization of clustered MHC-linked olfactory receptor genes in human and mouse

Olfactory receptor (OR) loci frequently cluster and are present on most human chromosomes. They are members of the seven transmembrane receptor (7-TM) superfamily and, as such, are part of one of the largest mammalian multigene families, with an estimated copy number of up to 1000 ORs per haploid genome. As their name implies, ORs are known to be involved in the perception of odors and possibly also in other, nonolfaction-related, functions. Here, we report the characterization of ORs that are part of the MHC-linked OR clusters in human and mouse (partial sequence only). These clusters are of particular interest because of their possible involvement in olfaction-driven mate selection. In total, we describe 50 novel OR loci (36 human, 14 murine), making the human MHC-linked cluster the largest sequenced OR cluster in any organism so far. Comparative and phylogenetic analyses confirm the cluster to be MHC-linked but divergent in both species and allow the identification of at least one ortholog that will be useful for future regulatory and functional studies. Quantitative feature analysis shows clear evidence of duplications of blocks of OR genes and reveals the entire cluster to have a genomic environment that is very different from its neighboring regions. Based on in silico transcript analysis, we also present evidence of extensive long-distance splicing in the 5'-untranslated regions and, for the first time, of alternative splicing within the single coding exon of ORs. Taken together with our previous finding that ORs are also polymorphic, the presented data indicate that the expression, function, and evolution of these interesting genes might be more complex than previously thought.

Characterization of Clustered MHC-Linked Olfactory Receptor Genes in Human and Mouse

Olfactory receptor (OR) loci frequently cluster and are present on most human chromosomes. They are members of the seven transmembrane receptor (7-TM) superfamily and, as such, are part of one of the largest mammalian multigene families, with an estimated copy number of up to 1000 ORs per haploid genome. As their name implies, ORs are known to be involved in the perception of odors and possibly also in other, nonolfaction-related, functions. Here, we report the characterization of ORs that are part of the MHC-linked OR clusters in human and mouse (partial sequence only). These clusters are of particular interest because of their possible involvement in olfaction-driven mate selection. In total, we describe 50 novel OR loci (36 human, 14 murine), making the human MHC-linked cluster the largest sequenced OR cluster in any organism so far. Comparative and phylogenetic analyses confirm the cluster to be MHC-linked but divergent in both species and allow the identification of at least one ortholog that will be useful for future regulatory and functional studies. Quantitative feature analysis shows clear evidence of duplications of blocks of OR genes and reveals the entire cluster to have a genomic environment that is very different from its neighboring regions. Based on in silico transcript analysis, we also present evidence of extensive long-distance splicing in the 5′-untranslated regions and, for the first time, of alternative splicing within the single coding exon of ORs. Taken together with our previous finding that ORs are also polymorphic, the presented data indicate that the expression, function, and evolution of these interesting genes might be more complex than previously thought.

[The sequence data described in this paper have been submitted to the EMBL nucleotide data library under accession nos.Z84475, Z98744, Z98745, AL021807, AL021808, AL022723, AL022727,AL031893, AL035402, AL035542, AL050328, AL050339, AL078630, AL096770,AL121944, AL133160, and AL133267.]

Conservation of Mhc class III region synteny between zebrafish and human as determined by radiation hybrid mapping

In the HLA, H2, and other mammalian MHC:, the class I and II loci are separated by the so-called class III region comprised of approximately 60 genes that are functionally and evolutionarily unrelated to the class I/II genes. To explore the origin of this island of unrelated loci in the middle of the MHC: 19 homologues of HLA class III genes, we identified 19 homologues of HLA class III genes as well as 21 additional non-class I/II HLA homologues in the zebrafish and mapped them by testing a panel of 94 zebrafish-hamster radiation hybrid cell lines. Six of the HLA class III and eight of the flanking homologues were found to be linked to the zebrafish class I (but not class II) loci in linkage group 19. The remaining homologous loci were found to be scattered over 14 zebrafish linkage groups. The linkage group 19 contains at least 25 genes (not counting the class I loci) that are also syntenic on human chromosome 6. This gene assembly presumably represents the pre-MHC: that existed before the class I/II genes arose. The pre-MHC: may not have contained the complement and other class III genes involved in immune response.

Compositional bias in DNA

Experimental approaches, as well as computer analysis on genomic sequences, have revealed a large variability in base composition between regions in the same genome or between genomes of different species. In most cases, however, the biological causes of these compositional biases remain unknown. The recent large increase in the availability of completely sequenced genomes can give new insight into evolution processes involved in these compositional biases.

An isochore transition in the NF1 gene region coincides with a switch in the extent of linkage disequilibrium

Whole-genome association studies will be a powerful tool to identify genes responsible for common human diseases. A crucial task for association-mapping studies is the evaluation of the relationship between linkage disequilibrium (LD) and physical distance for the genomic region under study. Since it is known that the extent of LD is nonuniformly distributed throughout the human genome, the required marker density has to be determined specifically for the region under study. These regions may be related to isochores and chromosomal bands, as indicated by earlier cytogenetic findings concerning chiasma distribution in meiosis. Therefore we analyzed the neurofibromatosis type 1 (NF1) gene region on chromosome 17q11.2, which is characterized by a nonuniform LD pattern and an L1-to-H2 isochore transition. Long-range LD within the NF1 gene was found to extend over 200 kb (D' = 0.937) in the L1 isochore, whereas, in the neighboring H2 isochore, no LD is apparent between markers spaced by 26 kb (D' = 0.144). Recombination frequencies derived from the LD are at.00019 (high LD) and.01659 (low LD) per megabase, the latter identical to the average value from segregation analysis. The boundary between these regions coincides precisely with a transition in the GC content of the sequences, with low values (37.2%) in the region with long-range LD and high values (51%) in the other. Our results suggest a correlation between the LD pattern and the isochores, at least in the NF1 region. If this correlation can be generalized, the marker densities required for association studies have to be adjusted to the regional GC content and may be chosen according to the isochores.

Isochores and the evolutionary genomics of vertebrates

The nuclear genomes of vertebrates are mosaics of isochores, very long stretches (>>300kb) of DNA that are homogeneous in base composition and are compositionally correlated with the coding sequences that they embed. Isochores can be partitioned in a small number of families that cover a range of GC levels (GC is the molar ratio of guanine+cytosine in DNA), which is narrow in cold-blooded vertebrates, but broad in warm-blooded vertebrates. This difference is essentially due to the fact that the GC-richest 10-15% of the genomes of the ancestors of mammals and birds underwent two independent compositional transitions characterized by strong increases in GC levels. The similarity of isochore patterns across mammalian orders, on the one hand, and across avian orders, on the other, indicates that these higher GC levels were then maintained, at least since the appearance of ancestors of warm-blooded vertebrates. After a brief review of our current knowledge on the organization of the vertebrate genome, evidence will be presented here in favor of the idea that the generation and maintenance of the GC-richest isochores in the genomes of warm-blooded vertebrates were due to natural selection.

Complete sequence and gene map of a human major histocompatibility complex. The MHC sequencing consortium

Here we report the first complete sequence and gene map of a human major histocompatibility complex (MHC), a region on chromosome 6 which is essential to the immune system. When it was discovered over 50 years ago the region was thought to specify histocompatibility genes, but their nature has been resolved only in the last two decades. Although many of the 224 identified gene loci (128 predicted to be expressed) are still of unknown function, we estimate that about 40% of the expressed genes have immune system function. Over 50% of the MHC has been sequenced twice, in different haplotypes, giving insight into the extraordinary polymorphism and evolution of this region. Several genes, particularly of the MHC class II and III regions, can be traced by sequence similarity and synteny to over 700 million years ago, clearly predating the emergence of the adaptive immune system some 400 million years ago. The sequence is expected to be invaluable for the identification of many common disease loci. In the past, the search for these loci has been hampered by the complexity of high gene density and linkage disequilibrium.