The genomic sequence and comparative analysis of the rat major histocompatibility complex

Interferon-triggered transcriptional cascades in the oligodendroglial lineage: a comparison of induction of MHC class II antigen between oligodendroglial progenitor cells and mature oligodendrocytes

Interferon-gamma induces major histocompatibility complex class II (MHC-II) in proliferating oligodendroglial progenitor cells (OPC), but to a much lesser extent in mature oligodendrocytes. Interferon-beta has virtually no effects on MHC-II induction even in OPC. Interferon-gamma-mediated transcriptional induction of CIITA, a critical regulator of MHC-II induction, was 6-fold lower in mature oligodendrocytes than in OPC, and entirely dependent on promoter IV, suggesting that the transcriptional activity of promoter IV is down-regulated after differentiation. The distinct difference in MHC-II induction between interferon-gamma and interferon-beta is attributed to transient interferon-beta-mediated activation of STAT1-IRF1 signaling compared to the sustained interferon-gamma-mediated activation.

Plasma from stored packed red blood cells and MHC class I antibodies causes acute lung injury in a 2-event in vivo rat model

Transfusion-related acute lung injury (TRALI) is the leading cause of transfusion death. We hypothesize that TRALI requires 2 events: (1) the clinical condition of the patient and (2) the infusion of antibodies against MHC class I antigens or the plasma from stored blood. A 2-event rat model was developed with saline (NS) or endotoxin (LPS) as the first event and the infusion of plasma from packed red blood cells (PRBCs) or antibodies (OX18 and OX27) against MHC class I antigens as the second event. ALI was determined by Evans blue dye leak from the plasma to the bronchoalveolar lavage fluid (BALF), protein and CINC-1 concentrations in the BALF, and the lung histology. NS-treated rats did not evidence ALI with any second events, and LPS did not cause ALI. LPS-treated animals demonstrated ALI in response to plasma from stored PRBCs, both prestorage leukoreduced and unmodified, and to OX18 and OX27, all in a concentration-dependent fashion. ALI was neutrophil (PMN) dependent, and OX18/OX27 localized to the PMN surface in vivo and primed the oxidase of rat PMNs. We conclude that TRALI is the result of 2 events with the second events consisting of the plasma from stored blood and antibodies that prime PMNs.

An informative set of SSLP markers and genomic profiles in the rat MHC, the RT1 complex

Almost 10,000 single nucleotide polymorphisms (SNPs) had been identified in the RT1 complex, the major histocompatibility complex of the rat, but less than approximately 0.5% have been characterized. In the context of the incomplete characterization of most SNPs, simple sequence length polymorphism (SSLP) marker development is still valuable for understanding the involvement of genes in the RT1 in controlling disease susceptibility, since SSLPs are user-friendly and cost-effective genetic markers in rat genome analysis. In this study, we developed a set of 67 SSLP markers, including 57 novel markers, to cover the entire RT1 complex and then created genetic profiles across 67 rat strains. These markers are located almost every 50 kb in the RT1 complex and show comparable polymorphism; the average number of alleles was 8.04 +/- 3.44 and the average polymorphic rate was 71 +/- 23%. Interestingly, markers failing to amplify polymerase chain reaction products were highly observed in all strains except for BN/SsNHsd, which suggests the existence of highly variable genomic sequences or genomic rearrangements in the RT1 region across rat strains. Based on the phylogenic tree and individual genotyping data, we identified 28 SSLP marker haplotypes in the RT1 region that roughly consisted of three genomic regions. These findings provided new insight into the genomic organization of the RT1 complex and we recognized the need of additional RT1 genome sequences in different strains. Owing to the accuracy and ease of determination, PCR-based SSLP genotyping could replace serological typing in genetic analyses and characterization of rat major histocompatibility.

Evolution of the opossum major histocompatibility complex: evidence for diverse alternative splice patterns and low polymorphism among class I genes

The opossum major histocompatibility complex (MHC) shares a similar organization with that of non-mammals while containing a diverse set of class I genes more like that of eutherian (placental) mammals. There are 11 class I loci in the opossum MHC region, seven of which are known to be transcribed. The previously described Monodelphis domestica (Modo)-UA1 and Modo-UG display characteristics consistent with their being classical and non-classical class I genes, respectively. Here we describe the characteristics of the remaining five transcribed class I loci (Modo-UE, -UK, -UI, -UJ and -UM). All five genes have peptide-binding grooves with low or no polymorphism, contain unpaired cysteines with the potential to produce homodimer formation and display genomic organizational features that would be unusual for classical class I loci. In addition, Modo-UJ and -UM were expressed in alternatively spliced mRNA forms, including a potentially soluble isoform of Modo-UJ. Thus, the MHC region of the opossum contains a single class I gene that is clearly classical and six other class I genes each with its own unique characteristics that probably perform roles other than or in addition to antigen presentation.

Both MHC and non-MHC genes regulate development of experimental neuropathic pain in rats

We have previously demonstrated that differences in neuropathic pain-like behaviors after sciatic nerve injury genetically maps to the major histocompatibility complex (MHC) in rats carrying RT1(c) or RT1(av1) haplotypes on the Piebald Virol Glaxo (PVG) background. In order to further explore the genetic contribution to neuropathic pain, we have here examined the MHC-congenic rat strains PVG-RT1(n) and PVG-RT1(av1) and the inbred strains PVG (RT1(c)) and Brown-Norway (BN; RT1(n)). All studied strains developed mechanical hypersensitivity (allodynia-like behavior) of the hind paw after photochemically induced sciatic nerve injury. However, the PVG-RT1(n) and PVG strains displayed significantly more allodynia than PVG-RT1(av1) and BN rats. In addition, the BN strain demonstrated an elevated threshold for the baseline response. The results demonstrate that both MHC and non-MHC genes influence experimental neuropathic pain in rats and also suggest that allelic variation contained in the RT1(av1) haplotype on the PVG background protects against neuropathic pain.

Contrasting genetic effects of major histocompatibility complex on ischemic peripheral nerve and spinal cord injury in female rats

We have recently shown that the major histocompatibility complex (MHC) exerts a regulatory influence on the development of neuropathic pain-like behaviors after partial sciatic nerve injury in male rats. In the present study, we assessed the role of the MHC in peripheral nerve injury-induced pain as well as central pain following spinal cord injury in female rats using the following inbred strains: Dark Agouti (DA; RT1(av1)), Piebald Virol Glaxo (PVG; RT1(c)) and in the MHC-congenic strain PVG-RT1(av1). In line with our previous observation in male rats, PVG-RT1(c) displayed more severe allodynia compared to the strains carrying the RT1(av1) haplotype (PVG-RT1(av1) and DA-RT1(av1)) following sciatic nerve injury in female rats. However, the MHC did not seem to impact the development of allodynia following spinal cord injury since the two congenic strains PVG-RT1(c) and PVG-RT1(av) did not differ after spinal cord injury. Interestingly, the DA-RT1(av1) strain displayed significantly more severe allodynia than both PVG strains and this difference was not explained by the extent of spinal cord injury. These results suggest that there are differences in the genetic mechanisms for neuropathic pain development following peripheral or central nervous system injury, both in regarding to the role of the MHC complex as well as non-MHC genes.

Impact of MHC mismatch and freezing on bone graft incorporation: an experimental study in rats

Cortical bone graft failure develops for poorly defined reasons, and the effects of the immune responses on the incorporation of an allograft are less clear. In a rat model of tibial allotransplantation, we have studied biometric and histological changes of the graft and the humoral immune response against it. We have also compared fresh with prefrozen grafts to study putative effects of freezing on the healing of the graft and the immune response against it. Fresh and frozen cortical bone grafts matched or mismatched for major histocompatibility complex antigens (syngeneic and allogeneic grafts) were implanted in an 8-mm segmental defect in the tibia. The construct was stabilized with intramedullary nailing. Incorporation of the graft was assessed with use of conventional radiography, micro computed tomography (CT(, biomechanical testing and histological examination. The immune response was evaluated by monitoring distribution of leukocytes in the blood and by measuring antibodies in a tailor-made fluorescence activating cell scanning (FACS( analysis. We found that the fresh syngeneic grafts were well integrated after 8 weeks with intact bone cells. In the fresh allogeneic grafts, all cells were dead with radiological signs of resorption, and mechanical testing indicated failure of incorporation. The frozen grafts showed poorer overall reconstruction than the fresh syngeneic grafts, but the incorporation was better than the fresh allogeneic grafts. A measurable alloantibody response was only detected after fresh allografting. The combined results suggest that freezing of bone allograft impedes the antibody response against major histocompatibility complex (MHC( antigens and improves incorporation, but frozen allografts still perform poorer than do frozen syngeneic grafts.

Genetic analysis of neuropathic pain-like behavior following peripheral nerve injury suggests a role of the major histocompatibility complex in development of allodynia

Neuropathic pain is a common consequence of damage to the nervous system. We here report a genetic analysis of development of neuropathic pain-like behaviors after unilateral photochemically-induced ischemic sciatic nerve injury in a panel of inbred rat strains known to display different susceptibility to autoimmune neuroinflammation. Pain behavior was initially characterized in Dark-Agouti (DA; RT1(av1)), Piebald Virol Glaxo (PVG; RT1(c)), and in the major histocompatibility complex (MHC)-congenic strain PVG-RT1(av1). All strains developed mechanical hypersensitivity (allodynia) following nerve injury. However, the extent and duration of allodynia varied significantly among the strains, with PVG displaying more severe allodynia compared to DA rats. Interestingly, the response of PVG-RT1(avRT1) was similar to that of DA, suggesting regulation by the MHC locus. This notion was subsequently confirmed in an F2 cohort derived from crossing of the PVG and PVG-RT1(av1)strains, where allodynia was reduced in homozygous or heterozygous carriers of the RT1(av1) allele in comparison to rats homozygous for the RT1(c) allele. These results indicate that certain allelic variants of the MHC could influence susceptibility to develop and maintain neuropathic pain-like behavior following peripheral nerve injury in rats.

Coevolutionary relationship between helminth diversity and MHC class II polymorphism in rodents

Parasite-mediated selection on major histocompatibility complex (MHC) genes has mainly been explored at the intraspecific level, although many molecular studies have revealed trans-species polymorphism. Interspecific patterns of MHC diversity might reveal factors responsible for the long-term evolution of MHC polymorphism. We hypothesize that host taxa harbouring high parasite diversity should exhibit high levels of MHC genetic diversity. We test this assumption using data on rodent species and their helminth parasites compiled from the literature. Controlling for similarity due to common descent, we present evidence indicating that high helminth species richness in rodent species is associated with increased MHC class II polymorphism. Our results are consistent with the idea that parasites sharing a long-term coevolutionary history with their hosts are the agents of selection explaining MHC polymorphism.

Major histocompatibility complex (Mhc) class Ib gene duplications, organization and expression patterns in mouse strain C57BL/6

BACKGROUND

The mouse has more than 30 Major histocompatibility complex (Mhc) class Ib genes, most of which exist in the H2 region of chromosome 17 in distinct gene clusters. Although recent progress in Mhc research has revealed the unique roles of several Mhc class Ib genes in the immune and non-immune systems, the functions of many class Ib genes have still to be elucidated. To better understand the roles of class Ib molecules, we have characterized their gene duplication, organization and expression patterns within the H2 region of the mouse strain C57BL/6.

RESULTS

The genomic organization of the H2-Q, -T and -M regions was analyzed and 21 transcribed Mhc class Ib genes were identified within these regions. Dot-plot and phylogenetic analyses implied that the genes were generated by monogenic and/or multigenic duplicated events. To investigate the adult tissue, embryonic and placental expressions of these genes, we performed RT-PCR gene expression profiling using gene-specific primers. Both tissue-wide and tissue-specific gene expression patterns were obtained that suggest that the variations in the gene expression may depend on the genomic location of the duplicated genes as well as locus specific mechanisms. The genes located in the H2-T region at the centromeric end of the cluster were expressed more widely than those at the telomeric end, which showed tissue-restricted expression in spite of nucleotide sequence similarities among gene paralogs.

CONCLUSION

Duplicated Mhc class Ib genes located in the H2-Q, -T and -M regions are differentially expressed in a variety of developing and adult tissues. Our findings form the basis for further functional validation studies of the Mhc class Ib gene expression profiles in specific tissues, such as the brain. The duplicated gene expression results in combination with the genome analysis suggest the possibility of long-range regulation of H2-T gene expression and/or important, but as yet unidentified nucleotide changes in the promoter or enhancer regions of the genes. Since the Mhc genomic region has diversified among mouse strains, it should be a useful model region for comparative analyses of the relationships between duplicated gene organization, evolution and the regulation of expression patterns.

Structure and expression of MHC class Ib genes of the central M region in rat and mouse: M4, M5, and M6

The M region at the telomeric end of the murine major histocompatibility complex (MHC) contains class I genes that are highly conserved in rat and mouse. We have sequenced a cosmid clone of the LEW rat strain (RT1 haplotype) containing three class I genes, RT1.M6-1, RT1.M4, and RT1.M5. The sequences of allelic genes of the BN strain (RT1n haplotype) were obtained either from cDNAs or genomic clones. For the coding parts of the genes few differences were found between the two RT1 haplotypes. In LEW, however, only RT1.M5 and RT1.M6 have open reading frames; whereas in BN all three genes were intact. In line with the findings in BN, transcription was found for all three rat genes in several tissues from strain Sprague Dawley. Protein expression in transfectants could be demonstrated for RT1.M6-1 using the monoclonal antibody OX18. By sequencing of transcripts obtained by RT-PCR, a second, transcribed M6 gene, RT1.M6-2, was discovered, which maps next to RT1.M6-1 outside of the region covered by the cosmid. In addition, alternatively spliced forms for RT1.M5 and RT1.M6 were detected. Of the orthologous mouse genes, H2-M4, H2-M5, and H2-M6, only H2-M5 has an open reading frame. Other important differences between the corresponding parts of the M region of the two species are insertion of long LINE repeats, duplication of RT1.M6, and the inversion of RT1.M5 in the rat. This demonstrates substantial evolutionary dynamics in this region despite conservation of the class I gene sequences themselves.

Disruption and pseudoautosomal localization of the major histocompatibility complex in monotremes

BACKGROUND

The monotremes, represented by the duck-billed platypus and the echidnas, are the most divergent species within mammals, featuring a flamboyant mix of reptilian, mammalian and specialized characteristics. To understand the evolution of the mammalian major histocompatibility complex (MHC), the analysis of the monotreme genome is vital.

RESULTS

We characterized several MHC containing bacterial artificial chromosome clones from platypus (Ornithorhynchus anatinus) and the short-beaked echidna (Tachyglossus aculeatus) and mapped them onto chromosomes. We discovered that the MHC of monotremes is not contiguous and locates within pseudoautosomal regions of two pairs of their sex chromosomes. The analysis revealed an MHC core region with class I and class II genes on platypus and echidna X3/Y3. Echidna X4/Y4 and platypus Y4/X5 showed synteny to the human distal class III region and beyond. We discovered an intron-containing class I pseudogene on platypus Y4/X5 at a genomic location equivalent to the human HLA-B,C region, suggesting ancestral synteny of the monotreme MHC. Analysis of male meioses from platypus and echidna showed that MHC chromosomes occupy different positions in the meiotic chains of either species.

CONCLUSION

Molecular and cytogenetic analyses reveal new insights into the evolution of the mammalian MHC and the multiple sex chromosome system of monotremes. In addition, our data establish the first homology link between chicken microchromosomes and the smallest chromosomes in the monotreme karyotype. Our results further suggest that segments of the monotreme MHC that now reside on separate chromosomes must once have been syntenic and that the complex sex chromosome system of monotremes is dynamic and still evolving.

Density-related changes in selection pattern for major histocompatibility complex genes in fluctuating populations of voles

Host-pathogen interactions are of particular interest in studies of the interplay between population dynamics and natural selection. The major histocompatibility complex (MHC) genes of demographically fluctuating species are highly suitable markers for such studies, because they are involved in initiating the immune response against pathogens and display a high level of adaptive genetic variation. We investigated whether two MHC class II genes (DQA1, DRB) were subjected to contemporary selection during increases in the density of fossorial water vole (Arvicola terrestris) populations, by comparing the neutral genetic structure of seven populations with that estimated from MHC genes. Tests for heterozygosity excess indicated that DQA1 was subject to intense balancing selection. No such selection operated on neutral markers. This pattern of selection became more marked with increasing abundance. In the low-abundance phase, when populations were geographically isolated, both overall differentiation and isolation-by-distance were more marked for MHC genes than for neutral markers. Model-based simulations identified DQA1 as an outlier (i.e. under selection) in a single population, suggesting the action of local selection in fragmented populations. The differences between MHC and neutral markers gradually disappeared with increasing effective migration between sites. In the high-abundance year, DQA1 displayed significantly lower levels of overall differentiation than the neutral markers. This gene therefore displayed stronger homogenization than observed under drift and migration alone. The observed signs of selection were much weaker for DRB. Spatial and temporal fluctuations in parasite pressure and locus-specific selection are probably the most plausible mechanisms underlying the observed changes in selection pattern during the demographic cycle.

Comparative genomic structure of human, dog, and cat MHC: HLA, DLA, and FLA

Comparisons of the genomic structure of 3 mammalian major histocompatibility complexes (MHCs), human HLA, canine DLA, and feline FLA revealed remarkable structural differences between HLA and the other 2 MHCs. The 4.6-Mb HLA sequence was compared with the 3.9-Mb DLA sequence from 2 supercontigs generated by 7x whole-genome shotgun assembly and 3.3-Mb FLA draft sequence. For FLA, we confirm that 1) feline FLA was split into 2 pieces within the TRIM (member of the tripartite motif) gene family found in human HLA, 2) class II, III, and I regions were placed in the pericentromeric region of the long arm of chromosome B2, and 3) the remaining FLA was located in subtelomeric region of the short arm of chromosome B2. The exact same chromosome break was found in canine DLA structure, where class II, III, and I regions were placed in a pericentromeric region of chromosome 12 whereas the remaining region was located in a subtelomeric region of chromosome 35, suggesting that this chromosome break occurred once before the split of felid and canid more than 55 million years ago. However, significant differences were found in the content of genes in both pericentromeric and subtelomeric regions in DLA and FLA, the gene number, and amplicon structure of class I genes plus 2 other class I genes found on 2 additional chromosomes; canine chromosomes 7 and 18 suggest the dynamic nature in the evolution of MHC class I genes.

Association between the DQA MHC class II gene and Puumala virus infection in Myodes glareolus, the bank vole

Puumala virus (PUUV) is a hantavirus specifically harboured by the bank vole, Myodes (earlier Clethrionomys) glareolus. It causes a mild form of hemorrhagic fever with renal syndrome (HFRS) in humans, called Nephropathia epidemica (NE). The clinical severity of NE is variable among patients and depends on their major histocompatibility complex (MHC) genetic background. In this study we investigated the potential role of class II MHC gene polymorphism in the susceptibility/resistance to PUUV in the wild reservoir M. glareolus. We performed an association study between the exon 2 of the DQA gene and PUUV antibodies considering a natural population of bank voles. Because immune gene polymorphism is likely to be driven by multiple parasites in the wild, we also screened bank voles for other potential viral and parasitic infections. We used multivariate analyses to explore DQA polymorphism/PUUV associations while considering the potential antagonist and/or synergistic effects of the whole parasite community. Our study suggests links between class II MHC characteristics and viral infections including PUUV and Cowpox virus. Several alleles are likely to be involved in the susceptibility or in the resistance of bank voles to these infections. Alternatively, heterozygosity does not seem to be associated with PUUV or any other parasite infections. This result thus provides no evidence in favour of the hypothesis of selection through overdominance. Finally this multivariate approach reveals a strong antagonism between ectoparasitic mites and PUUV, suggesting direct or indirect immunogenetic links between infections by these parasites. Other datasets are now required to confirm these results and to test whether the associations vary in space and/or time.

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.

Adenovirus subversion of immune surveillance, apoptotic and growth regulatory pathways: a model for tumorigenesis

The adenovirus system provides a novel model for evaluating the roles of multiple factors involved in tumour progression. In common with other DNA tumour viruses, adenovirus employs a variety of strategies to evade immune surveillance and perturbs cellular apoptotic and growth regulatory pathways to ensure efficient replication of progeny virions. Such subversion of cellular networks is also found in tumour cells. The mechanism behind the avoidance of immune surveillance and the extent of cellular network interference achieved by adenovirus is still being uncovered and is predicted to have ramifications for the design of cancer therapeutics.

A BAC clone-based physical map of ovine major histocompatibility complex

An ovine bacterial artificial chromosome (BAC) library containing 190,000 BAC clones was constructed and subsequently screened to construct a BAC-based physical map for the ovine major histocompatibility complex (MHC). Two hundred thirty-three BAC clones were selected by 84 overgo probes designed on human, mouse, and swine MHC sequence homologies. Ninety-four clones were ordered by DNA fingerprinting to form contigs I, II, and III that correspond to ovine MHC class I-class III, class IIa, and class IIb. The minimum tiling paths of contigs I, II, and III are 15, 4, and 4 BAC clones, spanning approximately 1900, 400, and 300 kb, respectively. The order and orientation of most BAC clones in each contig were confirmed by BAC-end sequencing. An open gap exists between class IIa and class III. This work helps to provide a foundation for detailed study of ovine MHC genes and of evolution of MHCs in mammals.

The genomic sequence and analysis of the swine major histocompatibility complex

We describe the generation and analysis of an integrated sequence map of a 2.4-Mb region of pig chromosome 7, comprising the classical class I region, the extended and classical class II regions, and the class III region of the major histocompatibility complex (MHC), also known as swine leukocyte antigen (SLA) complex. We have identified and manually annotated 151 loci, of which 121 are known genes (predicted to be functional), 18 are pseudogenes, 8 are novel CDS loci, 3 are novel transcripts, and 1 is a putative gene. Nearly all of these loci have homologues in other mammalian genomes but orthologues could be identified with confidence for only 123 genes. The 28 genes (including all the SLA class I genes) for which unambiguous orthology to genes within the human reference MHC could not be established are of particular interest with respect to porcine-specific MHC function and evolution. We have compared the porcine MHC to other mammalian MHC regions and identified the differences between them. In comparison to the human MHC, the main differences include the absence of HLA-A and other class I-like loci, the absence of HLA-DP-like loci, and the separation of the extended and classical class II regions from the rest of the MHC by insertion of the centromere. We show that the centromere insertion has occurred within a cluster of BTNL genes located at the boundary of the class II and III regions, which might have resulted in the loss of an orthologue to human C6orf10 from this region.

A novel discoidin domain receptor 1 (Ddr1) transcript is expressed in postmeiotic germ cells of the rat testis depending on the major histocompatibility complex haplotype

The Ddr1 gene encoding the discoidin domain receptor 1 (DDR1), a member of a small subfamily of receptor tyrosine kinases, is known to be involved in differentiation, proliferation, and cell adhesion. The extracellular discoidin domain is responsible for the binding of the ligand collagen. As the human homologue, the rat Ddr1 gene consists of 17 exons and is located in the major histocompatibility complex, the RT1 complex in rats. A novel testis-specific Ddr1 transcript of 3.5 kb is described here which is expressed specifically in the postmeiotic germ cells of the rat testis. The exons 1 to 5 are missing in this transcript and the putative protein would lack the discoidin domain and parts of the stalk region. The expression level of both, the full-length 4.3 kb and the novel 3.5 kb Ddr1 transcript, is dependent on the RT1 haplotype. In the RT1(avl) haplotype, carried by DA and LEW.1AV1 rats, the 3.5 kb Ddr1 transcript is completely missing. This might be explained by the lack of four nucleotides GGGC in the RT1(avl) haplotype, which appear to contribute to a SP1 binding site in intron 5 of the Ddr1 gene in the presumed testis-specific alternative promoter region of the 3.5 kb Ddr1 transcript. In addition, two novel exons in the 5'-untranslated region of the Ddr1 gene were found that give rise to further alternative Ddr1 transcripts. Interestingly, the 3.5 kb Ddr1 transcript is not only expressed in a cell type-specific manner in postmeiotic germ cells but also controlled by the RT1 haplotype.

Duplication, balancing selection and trans-species evolution explain the high levels of polymorphism of the DQA MHC class II gene in voles (Arvicolinae)

Major histocompatibility complex (MHC) genes play important role in host-parasite interactions and parasites are crucial factors influencing the population dynamics of hosts. We described the structure and diversity of exon 2 of the MHC class II DQA gene in three species of voles (Arvicolinae) exhibiting regular multi-annual fluctuations of population density and analysed the processes leading to the observed MHC polymorphism. By using cloning-sequencing methodology and capillary electrophoresis-single strand conformation polymorphism, we described seven sequences in the water, eight in the common, and seven in the bank voles coming from an area of 70 km(2) around the Nozeroy canton in the Jura Mountains (Franche Comté, France). All exon 2 sequences translate to give unique amino acid sequences and positive selection was found to act very intensively on antigen binding sites. We documented the presence of recombination at vole DQA region but its importance in generating allelic polymorphism seems to be relatively limited. For the first time within rodents, we documented the duplication of the DQA gene in all three species with both copies being transcriptionally active. Phylogenetic analysis of allelic sequences revealed extensive trans-species polymorphism within the subfamily although no alleles were shared between species in our data set. We discuss possible role of parasites in forming the recent polymorphism pattern of the DQA locus in voles.

The orthology of HLA-E and H2-Qa1 is hidden by their concerted evolution with other MHC class I molecules

BACKGROUND

Whether MHC molecules undergo concerted evolution or not has been the subject of a long-standing debate.

RESULTS

By comparing sequences of eight functional homologues of HLA-E from primates and rodents with those of MHC class Ia molecules from the same eight species, we find that different portions of MHC class I molecules undergo different patterns of evolution. By focusing our analyses sequentially on these various portions, we have obtained clear evidence for concerted evolution of MHC class I molecules, suggesting the occurrence of extensive interallelic and intergenic exchanges. Intra-species homogenisation of sequences is particularly noticeable at the level of exon 4, which codes for the alpha3 domain, but our results suggest that homogenisation also concerns certain residues of the alpha1-alpha2 codomain that lie outside the antigen recognition site.

CONCLUSION

A model is presented in which Darwinian selective pressures due to pathogens could, at the same time, favour diversification of MHC class Ia molecules and promote concerted evolution of separate loci by spreading advantageous motifs arising by mutations in individual MHC molecules to other alleles and to other loci of the MHC region. This would also allow MHC molecules to co-evolve with the proteins with which they interact to fulfil their functions of antigen presentation and regulation of NK cell activity. One of the raisons d'être of the MHC may therefore be to favour at the same time both diversification of MHC class Ia molecules and homogenisation of the whole pool of MHC class I molecules (Ia and Ib) involved in antigen presentation.

REVIEWERS

This article was reviewed by Stephan Beck, Lutz Walter and Pierre Pontarotti.

Additional human β2-microglobulin curbs HLA–B27 misfolding and promotes arthritis and spondylitis without colitis in male HLA–B27–transgenic rats

OBJECTIVE

Ankylosing spondylitis and related spondylarthritides are associated with HLA-B27, and also with intestinal inflammation, by unknown mechanisms. The folded HLA-B27 molecule is a trimer of heavy chain, beta2-microglobulin (beta2m), and short peptide. However, B27 heavy chain has an unusual propensity to misfold and trigger the unfolded protein response (UPR). This study was initiated to test the hypothesis that B27 misfolding plays a role in the pathogenesis of spondylarthritis.

METHODS

Rats with high transgene copy numbers of HLA-B27 heavy chain together with human beta2m (Hubeta2m) spontaneously develop colitis, peripheral arthritis, and occasional spondylitis, and those with lower transgene copy numbers remain healthy. We crossed disease-prone and healthy HLA-B27/Hubeta2m-transgenic rat lines with a healthy line, 283-2, carrying only the Hubeta2m transgene. HLA-B27 assembly was assessed by pulse-chase analysis of B27 molecules, and UPR triggering was assessed by measuring BiP/Grp78 messenger RNA (mRNA) in splenic concanavalin A blasts. Surface expression of B27 and Hubeta2m was determined by flow cytometry. Disease manifestations were identified by clinical observation, histology, and measurement of cytokine mRNA.

RESULTS

The extra Hubeta2m from the 283-2 line significantly reduced B27 misfolding and UPR triggering. Unexpectedly, however, F1 male offspring of the healthy 21-3 line crossed with the 283-2 line showed a high prevalence, severity, and duration of arthritis and spondylitis, in the absence of colitis. The arthropathy showed many features characteristic of human spondylarthritis.

CONCLUSION

These results suggest that B27 misfolding is associated with intestinal inflammation, but that neither B27 misfolding nor intestinal inflammation is critical to the development of B27-associated arthropathy.

Genomic sequence analysis of the 238-kb swine segment with a cluster of TRIM and olfactory receptor genes located, but with no class I genes, at the distal end of the SLA class I region

Continuous genomic sequence has been previously determined for the swine leukocyte antigen (SLA) class I region from the TNF gene cluster at the border between the major histocompatibility complex (MHC) class III and class I regions to the UBD gene at the telomeric end of the classical class I gene cluster (SLA-1 to SLA-5, SLA-9, SLA-11). To complete the genomic sequence of the entire SLA class I genomic region, we have analyzed the genomic sequences of two BAC clones carrying a continuous 237,633-bp-long segment spanning from the TRIM15 gene to the UBD gene located on the telomeric side of the classical SLA class I gene cluster. Fifteen non-class I genes, including the zinc finger and the tripartite motif (TRIM) ring-finger-related family genes and olfactory receptor genes, were identified in the 238-kilobase (kb) segment, and their location in the segment was similar to their apparent human homologs. In contrast, a human segment (alpha block) spanning about 375 kb from the gene ETF1P1 and from the HLA-J to HLA-F genes was absent from the 238-kb swine segment. We conclude that the gene organization of the MHC non-class I genes located in the telomeric side of the classical SLA class I gene cluster is remarkably similar between the swine and the human segments, although the swine lacks a 375-kb segment corresponding to the human alpha block.

Association of the truncating splice site mutation in BTNL2 with multiple sclerosis is secondary to HLA-DRB1*15

The major histocompatibility complex human leukocyte antigen (HLA)-DRB1*15 (DR2) haplotype is strongly associated with risk of multiple sclerosis (MS). The primary susceptibility has been localized to only approximately 200 kb encompassing the HLA-DR and -DQ loci. Further dissection of disease association with this region is demanding because of the high levels of linkage disequilibrium (LD). Recently, evidence was obtained for the involvement of a gene, potentially encoding an immune co-receptor, in another DR2-associated inflammatory condition, sarcoidosis. The implicated gene, BTNL2, is adjacent to DR and is in strong LD with HLA-DRB1. This fact, combined with a sequence relationship between BTNL2 and myelin oligodendrocyte glycoprotein, an autoantigen associated with MS, makes the gene an attractive candidate. To determine whether BTNL2 contributes to MS, we genotyped 1136 well-characterized MS families from the UK and the USA, as well as an African-American case-control data set, making this among the largest genetic studies in MS. Family-based and case-control association studies were performed for the BTNL2 and HLA-DRB1 loci. In all family data sets, the protein-truncating allele of BTNL2, implicated in sarcoidosis, was significantly over-transmitted to cases (combined data sets: global P=2.4x10(-11)). Given that the protein-truncating allele of BTNL2 virtually always occurred with DRB1*15, an effect could only be tested in DRB1*15-negative individuals or pedigrees. However, despite adequate power to detect an independent association, no difference in transmission of BTNL2 alleles or genotypes was observed in DRB1*15-negative individuals with MS. Conditional logistic regression modeling also strongly supported the conclusion that BTNL2 does not confer additional disease risk. The association of BTNL2 with MS observed in the African-American data set was also secondary to the primary DRB1*15 association.

HLA genomics in the third millennium

The MHC region contains several unique characteristics that set it apart as the most important region in the vertebrate genome in relation to disease. Recent data fit with the long-held view that the polymorphism of this region is driven by resistance to infection, although this is not yet proven. Interestingly, the MHC gene complex is associated with most, if not all, of the common autoimmune conditions. It has been difficult to identify the precise MHC genes associated with infection and autoimmunity, mainly because of the strong linkage disequilibrium over the region. Over the past few years, tools have been developed in an attempt to overcome these problems, including multiple fully sequenced MHC haplotypes, which have led to high-density hapmaps. In conjunction with large well-documented patient/control groups and sophisticated statistical methods these advances are starting to provide a comprehensive view of the genetics of the HLA region and disease susceptibility.

Comparative genomics of natural killer cell receptor gene clusters

Many receptors on natural killer (NK) cells recognize major histocompatibility complex class I molecules in order to monitor unhealthy tissues, such as cells infected with viruses, and some tumors. Genes encoding families of NK receptors and related sequences are organized into two main clusters in humans: the natural killer complex on Chromosome 12p13.1, which encodes C-type lectin molecules, and the leukocyte receptor complex on Chromosome 19q13.4, which encodes immunoglobulin superfamily molecules. The composition of these gene clusters differs markedly between closely related species, providing evidence for rapid, lineage-specific expansions or contractions of sets of loci. The choice of NK receptor genes is polarized in the two species most studied, mouse and human. In mouse, the C-type lectin-related Ly49 gene family predominates. Conversely, the single Ly49 sequence is a pseudogene in humans, and the immunoglobulin superfamily KIR gene family is extensive. These different gene sets encode proteins that are comparable in function and genetic diversity, even though they have undergone species-specific expansions. Understanding the biological significance of this curious situation may be aided by studying which NK receptor genes are used in other vertebrates, especially in relation to species-specific differences in genes for major histocompatibility complex class I molecules.

Unprecedented intraspecific diversity of the MHC class I region of a teleost medaka, Oryzias latipes

The major histocompatibility complex (MHC) is present at a single chromosomal locus of all jawed vertebrate analyzed so far, from sharks to mammals, except for teleosts whose orthologs of the mammalian MHC-encoded genes are dispersed at several chromosomal loci. Even in teleosts, several class IA genes and those genes directly involved in class I antigen presentation preserve their linkage, defining the teleost MHC class I region. We determined the complete nucleotide sequence of the MHC class I region of the inbred HNI strain of medaka, Oryzias latipes (northern Japan population-derived), from four overlapping bacterial artificial chromosome (BAC) clones spanning 540,982 bp, and compared it with the published sequence of the corresponding region of the inbred Hd-rR strain of medaka (425,935 bp, southern Japan population-derived) as the first extensive study of intraspecies polymorphisms of the ectotherm MHC regions. A segment of about 100 kb in the middle of the compared sequences encompassing two class Ia genes and two immunoproteasome subunit genes, PSMB8 and PSMB10, was so divergent between these two inbred strains that a reliable sequence alignment could not be made. The rest of the compared region (about 320 kb) showed a fair correspondence, and an approximately 96% nucleotide identity was observed upon gap-free segmental alignment. These results indicate that the medaka MHC class I region contains an approximately 100-kb polymorphic core, which is most probably evolving adaptively by accumulation of point mutations and extensive genetic rearrangements such as insertions, deletions, and duplications.

The evolutionary ecology of the major histocompatibility complex

The major histocompatibility complex (MHC) has become a paradigm for how selection can act to maintain adaptively important genetic diversity in natural populations. Here, we review the contribution of studies on the MHC in non-model species to our understanding of how selection affects MHC diversity, emphasising how ecological and ethological processes influence the tempo and mode of evolution at the MHC, and conversely, how variability at the MHC affects individual fitness, population dynamics and viability. We focus on three main areas: the types of information that have been used to detect the action of selection on MHC genes; the relative contributions of parasite-mediated and sexual selection on the maintenance of MHC diversity; and possible future lines of research that may help resolve some of the unanswered issues associated with MHC evolution.

Structure of a pheromone receptor-associated MHC molecule with an open and empty groove

Neurons in the murine vomeronasal organ (VNO) express a family of class Ib major histocompatibility complex (MHC) proteins (M10s) that interact with the V2R class of VNO receptors. This interaction may play a direct role in the detection of pheromonal cues that initiate reproductive and territorial behaviors. The crystal structure of M10.5, an M10 family member, is similar to that of classical MHC molecules. However, the M10.5 counterpart of the MHC peptide-binding groove is open and unoccupied, revealing the first structure of an empty class I MHC molecule. Similar to empty MHC molecules, but unlike peptide-filled MHC proteins and non-peptide–binding MHC homologs, M10.5 is thermally unstable, suggesting that its groove is normally occupied. However, M10.5 does not bind endogenous peptides when expressed in mammalian cells or when offered a mixture of class I–binding peptides. The F pocket side of the M10.5 groove is open, suggesting that ligands larger than 8–10-mer class I–binding peptides could fit by extending out of the groove. Moreover, variable residues point up from the groove helices, rather than toward the groove as in classical MHC structures. These data suggest that M10s are unlikely to provide specific recognition of class I MHC–binding peptides, but are consistent with binding to other ligands, including proteins such as the V2Rs.

MHC-like protein M10.5 is expressed in the vomeronasal organ. The structure does not bind endogenous class I-binding peptides, but is thought to interact with a larger V2R pheromone receptor.

Two structurally related rat Ly49 receptors with opposing functions (Ly49 stimulatory receptor 5 and Ly49 inhibitory receptor 5) recognize nonclassical MHC class Ib-encoded target ligands

The Ly49 family of lectin-like receptors in rodents includes both stimulatory and inhibitory members. Although NK alloreactivity in mice is regulated primarily by inhibitory Ly49 receptors, in rats activating Ly49 receptors are equally important. Previous studies have suggested that activating rat Ly49 receptors are triggered by polymorphic ligands encoded within the nonclassical class Ib region of the rat MHC, RT1-CE/N/M, while inhibitory Ly49 receptors bind to widely expressed classical class Ia molecules encoded from the RT1-A region. To further investigate rat Ly49-mediated regulation of NK alloreactivity, we report in this study the identification and characterization of two novel paired Ly49 receptors that we have termed Ly49 inhibitory receptor 5 (Ly49i5) and Ly49 stimulatory receptor 5 (Ly49s5). Using a new mAb (mAb Fly5), we showed that Ly49i5 is an inhibitory receptor that recognizes ligands encoded within the class Ib region of the u and l haplotypes, while the structurally related Ly49s5 is an activating receptor that recognizes class Ib ligands of the u haplotype. Ly49s5 is functionally expressed in the high NK-alloresponder PVG strain, but not in the low alloresponder BN strain, in which it is a pseudogene. Ly49s5 is hence not responsible for the striking anti-u NK alloresponse previously described in BN rats (haplotype n), which results from repeated alloimmunizations with u haplotype cells. The present studies support the notion of a complex regulation of rat NK alloreactivity by activating and inhibitory Ly49 members, which may be highly homologous in the extracellular region and bind similar class Ib-encoded target ligands.

The genes and gene organization of the Ly49 region of the rat natural killer cell gene complex

We here report the cDNA sequences of 11 new rat Ly49 genes with full and three with incomplete open reading frames. Although obtained from different inbred rat strains, these as well as six previously published cDNA represent non-allelic genes matching different loci in the Brown Norway (BN) rat genome, which is predicted to contain 34 Ly49 loci distributed over the distal part of the NK cell gene complex. Some of the cloned genes appear to be mutated to non-function in the BN genome, which harbors additional genes with full open reading frames, suggesting at least 26 non-allelic functional Ly49 genes in the rat. Of the encoded receptors, 13 are predicted to be inhibitory, eight to be activating, whereas five may be both ('bifunctional'). Phylogenetic analysis bears evidence of a highly dynamic genetic region, in which only the most distally localized Ly49 gene has a clear-cut mouse ortholog. In phylograms, the majority of the genes cluster into three subgroups with the genes mapping together, defining three chromosomal regions that seem to have undergone recent expansions. When comparing the lectin-like domains, the receptors form smaller subgroups, most containing at least one inhibitory and one activating or 'bifunctional' receptor, where close sequence similarities suggest recent homogenization events.

Genetics of the central MHC

The MHC, primarily known for its antigen-presenting class I and II molecules, harbours, within a central segment of less than 1 Mb, a dense collection of genes involved in various biological functions. Although MHC I and MHC II are principal players of adaptive immunity, several loci within this central (still called class III) MHC region encode members of the innate immune system. These include the long known factors of the complement system--potentially inhibitory and triggering natural killer receptors as well as stress proteins. Whether this physical proximity is fortuitous or functionally advantageous is an important question for the future of MHC genetics.

The Rat Expresses Two Complement Factor C4 Proteins, but Only One Isotype Is Expressed in the Liver

The complement component C4 is well known for its complex genetics in human and mouse where it is part of a tandemly duplicated module. For the rat, no such information had been available until recently. A C4 gene duplication could be identified also in the rat, but the duplicated module maps approximately 200 kb centromerically from the canonical C4-1 gene. In this study, we present the genomic organization of the two C4 gene-containing modules and the expression of the two C4 genes in the rat (Rattus norvegicus). The duplicated module contains an intact C4 gene as well as Cyp21 and Stk19 pseudogenes. Quantitative mRNA expression analyses revealed that both C4 genes are transcribed in various organs and tissues, but displaying ample differences of C4-1 and C4-2 expression. Most notably, C4-2 is not expressed in the liver. At variance to the mouse, the expression of the rat C4 genes does not exhibit any sex dependency. By using two-dimensional gel electrophoresis and mass spectrometry, products of both C4 genes could be identified in rat serum samples. These two rat C4 isotypes are nearly identical, but differ in a functionally important amino acid residue that is known to influence the functional properties of the C4 isotypes in human.

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.

The feline major histocompatibility complex is rearranged by an inversion with a breakpoint in the distal class I region

In order to determine the genomic organization of the major histocompatibility complex (MHC) of the domestic cat (Felis catus), DNA probes for 61 markers were designed from human MHC reference sequences and used to construct feline MHC BAC contig map spanning ARE1 in the class II region to the olfactory receptor complex in the extended class I region. Selected BAC clones were then used to identify feline-specific probes for the three regions of the mammalian MHC (class II-class III-class I) for radiation hybrid mapping and fluorescent in situ hybridization to refine the organization of the domestic cat MHC. The results not only confirmed that the p-arm of domestic cat B2 is inverted relative to human Chromosome 6, but also demonstrated that one inversion breakpoint localized to the distal segment of the MHC class I between TRIM39 and TRIM26. The inversion thus disjoined the approximately 2.85 Mb of MHC containing class II-class III-class I (proximal region) from the approximately 0.50 Mb of MHC class I/extended class I region, such that TRIM39 is adjacent to the Chromosome B2 centromere and TRIM26 is adjacent to the B2 telomere in the domestic cat.

Considerable haplotypic diversity in the RT1-CE class I gene region of the rat major histocompatibility complex

The major histocompatibility complex (MHC) class I region extending between the Bat1 and Pou5f1 genes shows considerable genomic plasticity in mouse and rhesus macaque but not in human haplotypes. In the rat, this region is known as the RT1-CE region. The recently published rat MHC sequence gave rise to a complete set of class I gene sequences in a single MHC haplotype, namely the RT1(n) haplotype of the widely used BN inbred strain. To study the degree of genetic diversity, we compared the RT1-CE region-derived class I genes of the RT1(n) haplotype with class I sequences of other rat haplotypes. By using phylogenetic tree analyses, we obtained evidence for extensive "presence and absence" polymorphisms of single loci and even small subfamilies of class I genes in the rat. Alleles of RT1-CE region class I genes could also be identified, but the rate of allelic nucleotide substitutions appeared rather low, indicating that the diversity in the RT1-CE region is mainly based on genomic plasticity.

Open reading frame sequencing and structure-based alignment of polypeptides encoded by RT1-Bb, RT1-Ba, RT1-Db, and RT1-Da alleles

MHC class II genes are major genetic components in rats developing autoimmunity. The majority of rat MHC class II sequencing has focused on exon 2, which forms the first external domain. Sequence of the complete open reading frame for rat MHC class II haplotypes and structure-based alignment is lacking. Herein, the complete open reading frame for RT1-Bbeta, RT1-Balpha, RT1-Dbeta, and RT1-Dalpha was sequenced from ten different rat strains, covering eight serological haplotypes, namely a, b, c, d, k, l, n, and u. Each serological haplotype was unique at the nucleotide level of the sequenced RT1-B/D region. Within individual genes, the number of alleles identified was seven, seven, six, and three and the degree of amino-acid polymorphism between allotypes for each gene was 22%, 16%, 19%, and 0.4% for RT1-Bbeta, RT1-Balpha, RT1-Dbeta, and RT1-Dalpha, respectively. The extent and distribution of amino-acid polymorphism was comparable with mouse and human MHC class II. Structure-based alignment identified the beta65-66 deletion, the beta84a insertion, the alpha9a insertion, and the alpha1a-1c insertion in RT1-B previously described for H2-A. Rat allele-specific deletions were found at RT1-Balpha76 and RT1-Dbeta90-92. The mature RT1-Dbeta polypeptide was one amino acid longer than HLA-DRB1 due to the position of the predicted signal peptide cleavage site. These data are important to a comprehensive understanding of MHC class II structure-function and for mechanistic studies of rat models of autoimmunity.