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

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

Segmental exchange between MHC class I genes in a higher primate: recombination in the gorilla between the ancestor of a human non-functional gene and an A locus gene

Classical human major histocompatibility complex (MHC) class I molecules are the products of highly diverse gene loci. It has been suggested that segmental exchange may play a role in the generation of diversity at the antigen recognition site of MHC class I molecules. Here we present the cloning, sequencing and expression of two gorilla A locus cDNAs. One of these cDNAs shows remarkable similarity to the non-functional HLA-AR locus gene (5.4-LBF) only in exon 2. The remainder of the cDNA, however, is most closely related to other classical higher primate A locus genes. This suggests that a segmental exchange may have occurred between the ancestor of the non-functional HLA-AR gene and a classical gorilla A locus gene. Furthermore, the recombination event resulting in Gogo-A3 has affected its antigen recognition site. These data, therefore, demonstrate that segmental exchange can generate diversity at the antigen recognition sites of primate MHC class I molecules and suggest that non-functional genes can contribute to the generation of diversity of classical MHC class I genes.The nucleotide sequence data reported in this paper have been submitted to the GenBank nucleotide sequence data base and have been assigned the accession numbers X54375 (Gogo-A3) and X54376 (Gogo-A4). Address correspondence and offprint requests to: D. I. Watkins.     

Unusually limited nucleotide sequence variation of the expressed major histocompatibility complex class I genes of a New World primate species (Saguinus oedipus)

Although major histocompatibility complex (MHC) class I molecules are, as a rule, highly polymorphic in mammalian species, those of the New World primate Saguinus oedipus (cotton-top tamarin) exhibit limited polymorphism. We have cloned and sequenced twelve MHC class I cDNAs from this species. Since cloned cotton-top tamarin cell lines express three to six MHC class I molecules, this species must have at least three functional MHC class I loci. There was, however, no evidence of locus-specific substitutions in the tamarin cDNAs. Unlike all other species studied, tamarin MHC class I cDNAs displayed limited nucleotide sequence variation. The sequence similarity between the two most divergent tamarin cDNAs was 95%. To ensure that the polymerase chain reaction (PCR) primers employed in these studies had amplified all of the tamarins' expressed MHC class I genes, we used another set of primers to amplify only exons 2 and 3 from RNA and DNA. PCR of genomic DNA resulted in the amplification of six distinct clones, of which only three were well expressed. Two of these nonexpressed genes were pseudogenes and the other was a nonclassical gene. Southern blot analysis demonstrated that the tamarin has 8–11 MHC class I genes, suggesting we had indeed cloned the majority of these genes. Cotton-top tamarins are, therefore, unique among mammalian species studied to date in that they express MHC class I molecules with limited nucleotide sequence variation.The nucleotide sequence data reported in this paper have been submitted to the GenBank nucleotide sequence database and have been assigned the accession numbers M38403-15.     

The MHC E locus in macaques is polymorphic and is conserved between macaques and humans

Although the functions of the molecules encoded by the classical MHC class I loci are well defined, no function has been ascribed to the molecules encoded by the non-classical MHC class I loci. To investigate the evolution and conservation of the non-classical loci, we cloned and sequenced HLA-E homologues in macaques. We isolated four E locus alleles from five rhesus monkeys and two E locus alleles from one cynomolgus monkey, which indicated that the E locus in macaques is polymorphic. We also compared the rate of nucleotide substitution in the second intron of the macaque and human E locus alleles with that of exons two and three. The rate of nucleotide substitution was significantly higher in the introns, which suggested that the E locus has evolved under selective pressure. Additionally, comparison of the rates of synonymous and non-synonymous substitutions in the peptide binding region versus the remainder of the molecule suggested that the codons encoding the amino acids in the peptide binding region had been conserved in macaques and humans over the 36 million years since macaques and humans last shared a common ancestor.The nucleotide sequence data reported in this paper have been submitted to the GenBank nucleotide sequence database and have been assigned the accession numbers UO2976–UO2981     

The dominant MHC class I gene is adjacent to the polymorphic<Emphasis Type="Italic"> TAP2</Emphasis> gene in the duck,<Emphasis Type="Italic"> Anas platyrhynchos</Emphasis>

We are investigating the expression and linkage of major histocompatibility complex (MHC) class I genes in the duck (Anas platyrhynchos) with a view toward understanding the susceptibility of ducks to two medically important viruses: influenza A and hepatitis B. In mammals, there are multiple MHC class I loci, and alleles at a locus are polymorphic and co-dominantly expressed. In contrast, in lower vertebrates the expression of one locus predominates. Southern-blot analysis and amplification of genomic sequences suggested that ducks have at least four loci encoding MHC class I. To identify expressed MHC genes, we constructed an unamplified cDNA library from the spleen of a single duck and screened for MHC class I. We sequenced 44 positive clones and identified four MHC class I sequences, each sharing approximately 85% nucleotide identity. Allele-specific oligonucleotide hybridization to a Northern blot indicated that only two of these sequences were abundantly expressed. In chickens, the dominantly expressed MHC class I gene lies adjacent to the transporter of antigen processing (TAP2) gene. To investigate whether this organization is also found in ducks, we cloned the gene encoding TAP2 from the cDNA library. PCR amplification from genomic DNA allowed us to determine that the dominantly expressed MHC class I gene was adjacent to TAP2. Furthermore, we amplified two alleles of the TAP2 gene from this duck that have significant and clustered amino acid differences that may influence the peptides transported. This organization has implications for the ability of ducks to eliminate viral pathogens.The nucleotide sequence data reported in this paper have been submitted to the GenBank nucleotide sequence database and have been assigned the accession numbers AY294416–22     

Variation in the number of expressed MHC genes in different cattle class I haplotypes

 Analysis of cattle major histocompatibility complex (MHC) (BoLA) class I gene expression using serological and biochemical methods has demonstrated a high level of polymorphism. However, analysis of class I cDNA sequences has failed to produce conclusive evidence concerning the number and nature of expressed genes. Such information is essential for detailed studies of cattle immune responses, and to increase our understanding of the mechanisms of MHC evolution. In this study a selective breeding programme has been used to generate a number of MHC homozygous cattle expressing common serologically defined class I specificities. Detailed analysis of five class I haplotypes was carried out, with transcribed class I genes identified and characterized by cDNA cloning, sequence analysis, and transfection/expression studies. Surface expression of the gene products (on lymphocytes) was confirmed using monoclonal antibodies of defined BoLA specificity. Phylogenetic analysis of available transcribed cattle MHC class I sequences revealed complex evolutionary relationships including possible evidence for recombination. The study of individual haplotypes suggests that certain groupings of related sequences may correlate with loci, but overall it was not possible to define the origin of individual alleles using this approach. The most striking finding of this study is that none of the cattle class I genes is consistently expressed, and that in contrast to human, haplotypes differ from one another in both the number and composition of expressed classical class I genes. Received: 15 February 1999 / Revised: 23 June 1999     

Association between BoLA and subclinical bovine leukemia virus infection in a herd of Holstein-Friesian cows

The role of the bovine major histocompatibility system (BoLA) in subclinical bovine leukemia virus (BLV) infection was investigated in a herd of Holstein-Friesian cows (n=240). The BoLA W8.1 allele was negatively associated with the presence of antibodies to the major BLV envelope glycoprotein, BLV-gp51 (corrected P<0.001, relative risk =0.31). These results suggest that a BoLA-linked gene(s) may influence the early spread of BLV infection. Since B cells are the primary target of BLV infection, we then determined the relationship between BoLA-A locus phenotypes and B-cell numbers in peripheral blood of seropositive and seronegative cows. There were no significant differences between BoLA-A alleles for any hematological parameter in seronegative cows. Seropositive cows with the W12.1 allele had significantly greater absolute numbers of lymphocytes per microliter and B cells per microliter than did seropositive cows with other BoLA-A phenotypes (P<0.01, respectively). The average effect associated with the W12.1 allele in BLV-infected cows was an increase of 2010 B cells per microliter of whole blood relative to BLV-infected cows with other BoLA-A phenotypes. These results demonstrate that susceptibility to the polyclonal expansion of BLV-infected B lymphocytes is associated with the W12.1 allele in Holstein-Friesian cattle. Compared with results of a previous study in a herd of Shorthorn cattle, it appears that resistance and susceptibility to subclinical progression of BLV infection are associated with different BoLA-A locus alleles in different cattle breeds.Abbreviations used in this paper AGID agar gel immunodiffusion - BLV bovine leukemia virus - BoLA bovine lymphocyte antigen - EBL enzootic bovine leukosis - HLA human leukocyte antigen - MHC major histocompatibility complex - PL persistent lymphocytosis     

Identification of expressed bovine class I MHC genes at two loci and demonstration of physical linkage

A cDNA library prepared from lymphocytes of a cow (E98), homozygous at major histocompatibility complex (MHC) loci (BoLA phenotype w10, KN104), was screened with a bovine MHC class I probe. Of the cDNA clones isolated, two, (2.1 and 5.1) were selected and showed divergence at both 5 and 3 termini. E98 DNA was digested with rare-cutter enzymes (Sfi I, Mlu I, Not I, and Cla I) and fragments were size-separated by field inversion gel electrophoresis (FIGE). Hybridization with an entire class I cDNA probe revealed multiple fragments generated by each enzyme. When the 3 untranslated regions (UT) of 2.1 and 5.1 were used as probes, only one fragment was revealed in each digested sample, showing locus specificity of these probes in cattle. Further, DNA of transfected mouse fibroblasts L4 (expressing KN104) and L10 (expressing w10) hybridized to the 3UT regions of clones 2.1 and 5.1, respectively, Northern blot analysis of the mRNA of the L4 and L10 transfected cells provided further evidence that the cDNA clones 2.1 and 5.1 code for the BoLA-KN104 and BoLA-w10 class I molecules respectively, and thus these represent the products of two different genes. A long range physical mapping of the BoLA-w10 and KN104 genes was performed using FIGE analysis of DNA of and homozygous and an heterozygous animal. This analysis revealed that the BoLA-w10 and KN104 genes are separated by not more than 210 kilobases (kb) and that they are components of a multigene family spanning 1550 kb. As the] w10 gene is at the BoLA-A locus we assign the KN104 gene to a B locus.     

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 ∼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 nucleotide sequence data of HNI MHC class I region reported in this paper have been submitted to the DDBJ/EMBL/GenBank and were assigned the accession number AB183488.     

Description of bovine major histocompatibility complex class IIa haplotypes using parthenogenetic embryo‐derived cells

In this study, we report an approach to characterize individual BoLA haplotypes using cells from parthenogenetic bovine embryos derived from slaughterhouse ovaries. Eight of the 15 parthenogenetic embryos so obtained had not undergone meiotic recombination on the BoLA region and were suitable to describe BoLA haplotypes. Detailed analysis of the BoLA class IIa region identified seven different class IIa haplotypes, including six not previously described and two new alleles of BoLA‐DQA and one BoLA‐DQB. Our method provided reliable sources of homozygous DNA to describe BoLA haplotypes.     

Independent gene duplications,not concerted evolution,explain relationships among class I MHC genes of murine rodents

It has been claimed that class I MHC loci are homogenized within species by frequent events of interlocus genetic exchange (concerted evolution). Evidence for this process includes the fact that certain rat class I loci (including RT1.A) located centromeric to class II and class III are more similar to each other than to the mouse K locus (also centromeric to class II/class III). However, a phylogenetic analysis showed that the rat RT1.A locus is in fact orthologous to the mouse K1 pseudogene (also centromeric to class II/class III). Thus, two independent events of translocation of genes centromeric to class II/class III have occurred in the history of the murine rodents, at least one of which (involving the ancestor of RT1.A and K1) occurred prior to the divergence of rat and mouse. It was also found that the rat nonclassical class I gene RT.BM1 is orthologous to the mouse nonclassical gene 37 ^d. These results argue that intelocus genetic exchange does not occur at a rate sufficient to cause within-species homogenization of class I MHC loci.     

Identification of major histocompatibility complex genes in the guppy,Poecilia reticulata

The guppy, Poecilia reticulata, a teleostean fish of the order Cyprinodontiformes, has been used extensively in studies of host-parasite interactions, courtship behavior, and mating preference, as well as in ecological and evolutionary genetics. A related species was among the first poikilotherm vertebrates to be used in the study of histocompatibility genes. All these studies could benefit from the identification and characterization of the guppy major histocompatibility complex (Mhc) genes. Here, both class I and class II genes of the guppy are described. The number of expressed loci, as determined by representation of clones in a cDNA library, sequencing, and Southern blot analysis, may be low in both Mhc classes: combined evidence suggests that there may be one expressed class II locus only and one or two expressed class I loci. The variability of aquaristic guppy stocks is very low: only three and two genes have been detected at the class I and class II loci, respectively, in the stocks examined. This genetic paucity is most likely the consequence of breeding practices employed by aquarists and commercial establishments. Limited sampling of wild guppy populations revealed extensive Mhc polymorphism at loci of both classes in nature. Comparison of guppy Mhc sequences with those of other vertebrates has revealed the existence of a set of insertions/deletions which can be used as characters in cladistic analysis to infer phylogenetic relationships among vertebrate taxa and the Mhc genes themselves. These indels are particularly frequent in the regions coding for the loops of 1 and 2 domains of class I proteins.The nucleotide sequence data reported in this publication have been submitted to the EMBL nucleotide sequence database and have been assigned the accession numbers Z54076-Z54095     

Evolution of the major histocompatibility complex: a hundred-fold amplification of MHC class I genes in the African pigmy mouse Nannomys setulosus

The genome of the African murine rodent Nannomys setulosus was found to harbor several thousand major histocompatibility complex (MHC) class I genes instead of the 30–40 genes found in conventional laboratory mice, which are mostly of Mus musculus domesticus origin. Other genes of N. setulosus, either functionally or physically linked to class I genes, are not amplified. Amplified genes derive from as few as three ancestors and amplification has likely occured after the divergence of the two Nannomys species, N. setulosus and N. minutoides, which took place about three million years ago. Amplified genes are mostly pseudogenes. Statistical analysis of dinucleotide frequencies leads us to propose that inactivation of the genes has occured through the repeat induced mutation process, a possible newcomer in the evolution of the MHC.The nucleotide sequence data reported in this paper have been submitted to the EMBL nucleotide data library and have been assigned the accession numbers X61685 and X61686. Correspondence to : G. Gachelin.     

Expressed Peromyscus maniculatus (Pema) MHC class I genes: evolutionary implications and the identification of a gene encoding a Qa1-like antigen

To gain insight into the evolution of rodent major histocompatibility complex (MHC) class I genes and identify important (conserved) nonclassical class I (class Ib) gene products and residues in these proteins, sixPeromyscus maniculatus MHC (Pema) class I cDNA clones were isolated and sequenced. FivePema class I cDNAs appeared most similar to mouse and rat classical class I (class Ia) genes. One exhibited highest similarity to anH2 class Ib gene,H2-T23 (encoding the Qa1 antigen). Phylogenetic trees constructed withPema, RT1, andH2 class I sequences suggested that the lineages of some rodent class Ib genes (e.g.,T23 andT24) originated prior toMus andPeromyscus speciation [>50 million years (My) ago]. Sequences of four Qa1-like proteins from three species permitted the identification of ten Qa1-specific amino acids. On the basis of molecular modeling, three residues showed the potential to interact with T-cell receptors and three residues (all corresponding to polymorphic positions among H2 class Ia proteins) were predicted to influence antigen binding. The recognition of mouse Qa1 proteins by a subset of T-cells in influenced by a locus,Qdm, which encodes the H2-D leader peptide. One of thePema class I cDNA clones classified asH2-K, D/L-like (class Ia) is predicted to encode an identical peptide, implying that an antigen binding protein (Qa1) and the antigen to which it binds (the product ofQdm) has been conserved for over 50 My. The nucleotide sequence data reported in this paper have been submitted to the GenBank nucleotide sequence database and have been assigned the accession numbers U12822 (Pm13), U12885 (Pm41), U12886 (Pm52), U12887 (Pm62), U16846 (Pm11), and U16847 (Pm53)     

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.     

Analysis of genomic and expressed major histocompatibility class Ia and class II genes in a hexaploid Lake Tana African ‘large’ barb individual (<Emphasis Type="Italic">Barbus intermedius</Emphasis>)

Expression of too many co-dominant major histocompatibility complex (MHC) alleles is thought to be detrimental to proper functioning of the immune system. Polyploidy of the genome will increase the number of expressed MHC genes unless they are prone to a silencing mechanism. In polyploid Xenopus species, the number of MHC class I and II genes has been physically reduced, as it does not increase with higher ploidy genomes. In the zebrafish some class IIB loci have been silenced, as only two genomically bona fide loci, DAA/DAB and DEA/DEB, have been described. Earlier studies indicated a reduction in the number of genomic and expressed class II MHC genes in a hexaploid African large barb. This prompted us to study the number of MHC genes present in the genome of an African large barb individual (Barbus intermedius) in relation to those expressed, adopting the following strategy. Full-length cDNA sequences were generated from mRNA and compared with partial genomic class Ia and II sequences generated by PCR using the same primer set. In addition, we performed Southern hybridizations to obtain a verification of the number of class I and IIB genes. Our study revealed three ₂-microglobulin, five class Ia, four class IIA, and four class IIB genes at the genomic level, which were shown to be expressed in the hexaploid barb individual. The class Ia and class II data indicate that the ploidy status does not correlate with the presence and expression of these MHC genes.     

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

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

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

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

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