Allelic variation of ovine MHC class II DQA1 and DQA2 genes

In the present study we characterize allelic variation of polymorphic OLA-DQA1 and OLA-DQA2 genes in sheep. To achieve this, PCR primers were designed to independently amplify the second exons of OLA-DQA1 and OLA-DQA2 genes. Single strand conformation polymorphism (SSCP) gel analyses reveals that there are at least 12 distinct OLA-DQA2 sequences, 10 of which have been characterized by sequencing. Six distinct OLA-DQA1 alleles have been sequenced in sheep and we can detect at least seven DQA1 alleles, including a null allele, by SSCP analysis. The second exon of the OLA-DQA2 gene is more polymorphic than the equivalent region of the OLA-DQA1 gene. Thirty-two per cent of nucleotide and 49% of amino acid sites showed variation at the DQA2 locus, compared to 20% of nucleotide and 33% of amino acid sites for DQA1 . Phylogenetic analysis of DQA sequences from a number of species show that sheep DQA1 sequences group together and are more similar to bovine DQA1 sequences than to sheep DQA2 alleles. The majority of OLA-DQA2 sequences are on the same main branch of the phylogenetic tree as bovine DQA2 sequences. However, three sheep DQA2 sequences have a tendency to group with putative bovine DQA3 sequences rather than to other ovine DQA2 alleles. A variety of SSCP gel conditions were tried in order to develop a typing system for the OLA-DQA2 gene. We describe a set of PCR and SSCP conditions which distinguish between all known OLA-DQA2 alleles.     

Gene duplications and sequence polymorphism of bovine class II DQB genes

The genetic diversity of bovine class II DQB genes was investigated by polymerase chain reaction amplification and DNA sequencing. The first domain exon was amplified from genomic DNA samples representing 14 class II haplotypes, defined by restriction fragment length polymorphism (RFLP) analysis. The presence of a polymorphism in the copy number of DQB genes was confirmed since two DQB sequences were isolated from certain haplotypes. Four subtypes of bovine DQB genes were found. DQB1 is the major type and was found in almost all haplotypes. DQB2 is very similar DQB1 but was found only in the duplicated haplotypes DQ9 to 12. DQB3 and DQB4 are two quite divergent genes only present in certain duplicated haplotypes. The bovine DQB complexity thus resembles that in the human DRB region. Bovine DQB genes were found to be highly polymorphic as ten DQB1 alleles and four DQB2 alleles were identified. The observed sequence polymorphism correlated well with previously defined DQB RFLPs. Bovine and human DQB alleles show striking similarities at the amino acid level. In contrast, the frequency of silent substitutions is much higher in comparisons of DQB alleles between species than within species ruling out the possibility that any of the contemporary DQB alleles have been maintained since the divergence of humans and cattle. The frequency of silent substitutions between DQB alleles was markedly lower in cattle than in humans, in agreement with a previous comparison of human and bovine DRB alleles.     

Polymorphism of the IGHA gene in sheep

Genetic variation in immunoglobulin A, the most abundant immunoglobulin in mammalian cells, has not been reported in ruminants. In this study, variation in the immunoglobulin heavy alpha chain constant gene (IGHA) of sheep was investigated by amplification of a fragment that included the hinge coding sequence, followed by single-strand conformational polymorphism (SSCP) analysis and DNA sequencing. Three novel sequences, each characterized by unique SSCP banding patterns, were identified. One or two sequences were detected in individual sheep and all the sequences identified shared high homology to the published ovine and bovine IGHA sequences, suggesting that these sequences represent allelic variants of the IGHA gene in sheep. Sequence alignment showed that these sequences differed mainly in the 3′ end of exon 1 and in the coding sequence of the hinge region. There was either a deletion or an insertion of two codons in the hinge coding region in these allelic variants. Codon usage in the hinge coding region was quite different from that in the non-hinge coding regions of the gene, suggesting different evolution of the IGHA hinge sequence. Three novel amino acid sequences of ovine IGHA were also predicted, and variation in these sequences might not only affect antigen recognition but also susceptibility to cleavage by bacterial or parasitic proteases. Nucleotide sequence data reported in this paper have been submitted to the NCBI GenBank nucleotide sequence database and have been assigned the accession nos. AY956424–AY956426.     

The nucleotide sequence of bovine MHC class IIDQB andDRB genes

The nucleotide sequences of most of the exons and parts of the introns of twoBoLA-DQB genes and twoBoLA-DRB genes have been determined. The structure of these genes is very similar to that of human major histocompatibility complex (MHC) class II genes. The twoDQB genes probably represent true alleles. Based on the exons sequenced, bothDQB genes and one of theDRB genes seem to be functional. The otherDRB gene is a pseudogene; stopcodons are found in the exons encoding the second and transmembrane domain and, furthermore, a 2 base pair (bp) deletion has occured in the leader exon which places the initiation start codon out of frame. Also in this pseudogene, an almost perfect inverted repeat of 200 bp is found flanking the exon encoding the first domain, which might have been the result of a duplication/inversion event. The sequences presented in this paper do not contain any repetitions. Therefore, DNA fragments containing these sequences can be used as homologous bovine probes in restriction fragment length polymorphism (RFLP) analysis to study disease association in 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 M30002–M30014. Address correspondence and offprint requests to: M. A. M. Groenen.     

Diversity and duplication of DQB and DRB-like genes of the MHC in baleen whales (suborder: Mysticeti)

The molecular diversity and phylogenetic relationships of two class II genes of the baleen whale major histocompatibility complex were investigated and compared to toothed whales and out-groups. Amplification of the DQB exon 2 provided sequences showing high within-species and between-species nucleotide diversity and uninterrupted reading frames consistent with functional class II loci found in related mammals (e.g., ruminants). Cloning of amplified products indicated gene duplication in the humpback whale and triplication in the southern right whale, with average nucleotide diversity of 5.9 and 6.3%, respectively, for alleles of each species. Significantly higher nonsynonymous divergence at sites coding for peptide binding (32% for humpback and 40% for southern right) suggested that these loci were subject to positive (overdominant) selection. A population survey of humpback whales detected 23 alleles, differing by up to 21% of their inferred amino acid sequences. Amplification of the DRB exon 2 resulted in two groups of sequences. One was most similar to the DRB3 of the cow and present in all whales screened to date, including toothed whales. The second was most similar to the DRB2 of the cow and was found only in the bowhead and right whales. Both loci showed low diversity among species and apparent loss of function or altered function including interruption of reading frames. Finally, comparison of inferred protein sequence of the DRB3-like locus suggested convergence with the DQB, perhaps resulting from intergenic conversion or recombination.Electronic Supplementary Material Supplementary material is available for this article at     

Identification of Allelic Polymorphism in the Ovine Leptin Gene

Leptin is an adipocyte-derived hormone/cytokine that influences the physiological control of numerous biological functions and links nutritional status with both neuroendocrine and immune functions. In livestock, variation in the leptin (LEP) gene has been characterized in cattle and pig, but it has not been reported in sheep. In this study, variation in the exon 3 coding sequence of the ovine LEP gene was investigated by polymerase chain reaction–single-strand conformational polymorphism (PCR–SSCP) analysis and DNA sequencing. Five novel SSCP patterns, representing five different sequences, were identified under a combination of two different electrophoresis conditions. Either one or two different sequences were detected in individual sheep and all the sequences identified shared high homology with the LEP sequences from a variety of species, suggesting that these sequences represent alleles of the ovine LEP gene. Four single nucleotide polymorphisms (SNPs) were detected, and three of these resulted in amino acid changes. Variation detected here might have an impact on leptin activity and function.     

Characterization of the horse (<Emphasis Type="Italic">Equus caballus</Emphasis>)<Emphasis Type="Italic"> IGHA</Emphasis> gene

Nucleotide sequences of the immunoglobulin constant heavy chain genes of the horse have been described for IGHM, IGHG and IGHE genes, but not for IGHA. Here, we provide the nucleotide sequence of the genomic IGHA gene of the horse (Equus caballus), including its secretion region and the transmembrane exon. The equine IGHA gene shows the typical structure of a mammalian IGHA gene, with only three exons, separated by two introns of similar size. The hinge exon is located at the 5 end of the CH2 exon and encodes a hinge region of 11 amino acids, which contains five proline residues. The coding nucleotide sequence of the secreted form of the equine IGHA gene shares around 72% identity with the human IGHA1 and IGHA2 genes, as well as the bovine, ovine, porcine and canine IGHA genes, without distinct preference for any of these species. The same species also cluster together in a phylogenetic tree of the IGHA coding regions of various mammals, whereas rodent, rabbit, marsupial and monotreme IGHA genes each build a separate cluster.The nucleotide sequences reported in this paper have been assigned the EMBL/GenBank accession numbers AY247966 and AY351982     

Isolation, characterization and evolution of ovine major histocompatibility complex class II DRA and DQA genes

Four full-length ovine major histocompatibility complex (MHC) class II A cDNA clones coding for new alleles of DRA, DQA1 and DQA2 genes were isolated from two ovine Λgt10 cDNA libraries. The derived amino acid sequences of these clones resemble class II A molecules from other species in both size and structure. Restriction fragment length polymorphism analysis, using an Ovar-DRA probe on DNA from Merino and Romney sheep revealed only limited polymorphism in contrast to the high levels of polymorphism revealed by Ovar-DQA probes. Comparison of the predicted amino acid sequences for the three ovine A genes with class II A genes from five other species revealed that the most variable region of the molecule is the signal peptide. Although virtually every amino acid site shows variation, within or between species, there are some blocks of highly conserved residues. Within gene comparisons of nucleotide differences reveal that the greatest number of changes is found between the alleles of Ovar-DQA1 and -DQA2 genes and the least between Ovar-DRA1 alleles. Phylogenetic analysis of class IIA sequences from several species place DRA and DQA genes on two distinct branches, with Ovar-DRA1 and BOLA-DRA, and Ovar-DQA1 and BOLA-DQA being most similar on their respective branches.     

Sequence diversity and molecular evolution of the leukotoxin (lktA) gene in bovine and ovine strains of Mannheimia (Pasteurella) haemolytica

The molecular evolution of the leukotoxin structural gene (lktA) of Mannheimia (Pasteurella) haemolytica was investigated by nucleotide sequence comparison of lktA in 31 bovine and ovine strains representing the various evolutionary lineages and serotypes of the species. Eight major allelic variants (1.4 to 15.7% nucleotide divergence) were identified; these have mosaic structures of varying degrees of complexity reflecting a history of horizontal gene transfer and extensive intragenic recombination. The presence of identical alleles in strains of different genetic backgrounds suggests that assortative (entire gene) recombination has also contributed to strain diversification in M. haemolytica. Five allelic variants occur only in ovine strains and consist of recombinant segments derived from as many as four different sources. Four of these alleles consist of DNA (52.8 to 96.7%) derived from the lktA gene of the two related species Mannheimia glucosida and Pasteurella trehalosi, and four contain recombinant segments derived from an allele that is associated exclusively with bovine or bovine-like serotype A2 strains. The two major lineages of ovine serotype A2 strains possess lktA alleles that have very different evolutionary histories and encode divergent leukotoxins (5.3% amino acid divergence), but both contain segments derived from the bovine allele. Homologous segments of donor and recipient alleles are identical or nearly identical, indicating that the recombination events are relatively recent and probably postdate the domestication of cattle and sheep. Our findings suggest that host switching of bovine strains from cattle to sheep, together with inter- and intraspecies recombinational exchanges, has played an important role in generating leukotoxin diversity in ovine strains. In contrast, there is limited allelic diversity of lktA in bovine strains, suggesting that transmission of strains from sheep to cattle has been less important in leukotoxin evolution.     

Class II genes of miniature swine

Genomic clones corresponding to class II genes of theSLA ^c haplotype of miniature swine have been isolated and characterized. These genes have been grouped into seven non-overlapping clusters on the basis of restriction mapping. Ordering of exons within each cluster was accomplished by hybridization of Southern blots of restriction fragments with exon-specific probes. The two clusters (clusters 2 and 3) encoding theDRB andDQB genes were identified on the basis of hybridization with locus-specific 3 untranslated cDNA probes. Cluster 4 contained exons of bothDOB andDQB genes, the basis for which remains to be determined. The remaining four clusters (1, 5, 6, 7) were identified as containingDP, DR, andDO coding sequences, respectively, on the basis of sequence analysis. The porcine class II region appears very similar to that of man in number and nature of the class II genes identified and in the intron/exon organization of corresponding genes.The nucleotide sequence data reported in this paper have been submitted to the GenBank nucleotide sequence database and have been assigned the association number M29944. Address correspondence and offprint requests to: C. LeGuern.     

Diversity of the glycine/tyrosine-rich keratin-associated protein 6 gene (<Emphasis Type="Italic">KAP6</Emphasis>) family in sheep

Keratin-associated proteins (KAPs) are structural components of wool and variation in them may affect wool characteristics. In this study, we used PCR-SSCP to analyse the ovine KAP6 family which encodes glycine and tyrosine-rich KAPs. Five unique PCR-SSCP patterns were detected in the 250 sheep investigated. Between two and five patterns were observed in individual sheep and none with only one pattern was detected. This suggests the amplicons were heterogeneous and derived from more than one locus. To analyse these heterogeneous PCR amplicons, a sequencing approach using SSCP to separate individual amplified sequences, was developed. Using this approach, five DNA sequences (A–E) representing five unique PCR-SSCP patterns were obtained. D was identical to a published ovine KAP6-1 sequence (GenBank accession no. M95719), whereas the others were novel, but the closest homology was with KAP6 sequences from human, sheep, goats and cattle. The five ovine KAP6 sequences could be assigned into three distinct groups. B and D were identical to each other, with the exception of a 57-bp deletion/insertion and a single nucleotide polymorphism (SNP) in the 3′-UTR region. These appear to be allelic variants of ovine KAP6-1. A and C could form another group, as they were similar to each other (with only one synonymous SNP), but different to the other sequences. This group appears to be related to a sheep KAP6 amino acid sequence, and represent allelic variation at another KAP6 locus (designated KAP6-2). The remaining sequence E did not show high sequence homology with either the KAP6-1 or KAP6-2 sequences, but exhibited homology with a bovine KAP6-3 sequence, with the exception of a deletion/insertion of 30 nucleotides. This suggests that E represents ovine KAP6-3. This sequence was detected in only 11% of the sheep investigated, suggesting either a KAP6-3 null allele, or failure to amplify allleles. These results suggest that ovine KAP6 is a complex gene family, that is not only comprised multiple loci, but that is also polymorphic.     

Haplotypic Diversity Within the Ovine Calpastatin <Emphasis Type="Italic">(CAST)</Emphasis> Gene

Calpastatin (CAST) is a protein inhibitor that acts specifically on calpains and plays a regulatory role in postmortem beef tenderization and muscle proteolysis. Polymorphisms in the bovine CAST gene have been associated with meat tenderness, but little is known about how the ovine CAST gene may affect sheep meat quality traits. In this study, we selected two parts of the ovine CAST gene that have been previously reported to be polymorphic (region 1—part of intron 5 and exon 6, and region 2—part of intron 12), to investigate haplotype diversity across an extended region of the ovine gene. First, we developed a simple and efficient polymerase chain reaction-single-strand conformational polymorphism (PCR-SSCP) method for genotyping region 2, which allowed the detection of a novel allele as well as the three previously reported alleles. Next, we genotyped both regions 1 and 2 of the ovine CAST gene from a large number of sheep to determine the haplotypes present. Nine different haplotypes were found across this extended region of the ovine CAST gene and four haplotypes were identified that suggested historical recombination events within this gene. Haplotypes are typically more informative than single nucleotide polymorphisms (SNPs) for analyzing associations between genes and complex production traits, such as meat tenderness, but the potential for intragenic recombination within the ovine CAST may make finding associations challenging.     

Partial nucleotide sequence of a novel bovine major histocompatibility complex class II β-chain gene, BoLA-DIB

A bovine genomic clone that hybridized to HLA-DQ beta cDNA was isolated and fragments containing the beta 1, beta 2 and transmembrane (TM) exons subcloned. The nucleotide sequences of the exons and flanking intron regions were determined. Comparisons of these exon nucleotide sequences and derived amino acid sequences to human class II beta-chain sequences showed that this gene is only 77% identical to HLA-DQ beta and about 75% identical to bovine DQ beta-like genes. The exon sequences were more divergent from other class II beta-chain genes. However, structural features such as conserved cysteines and regions of amino acids strongly suggest this to be a class II beta-chain gene. When exon-containing fragments were used as hybridization probes on Southern blots of bovine genomic DNA digested with Eco RI or Pvu II, each exon hybridized to a single band. Based on these results we have referred to this gene as a novel bovine class II beta-chain gene, BoLA-DIB.     

Genetic Variation at Exon 2 of the MHC Class II DQB Locus in Blue Whale (Balaenoptera musculus) from the Gulf of California

The genes of the Major Histocompatibility Complex (MHC) play an important role in the vertebrate immune response and are among the most polymorphic genes known in vertebrates. In some marine mammals, MHC genes have been shown to be characterized by low levels of polymorphism compared to terrestrial taxa; this reduction in variation is often explained as a result of lower pathogen pressures in marine habitats. To determine if this same reduction in variation applies to the migratory population of blue whales (Balaenoptera musculus) that occurs in the Gulf of California, we genotyped a 172 bp fragment of exon 2 of the MHC Class II DQB locus for 80 members of this population. Twenty-two putatively functional DQB allotypes were identified, all of which were homologous with DQB sequences from other cetacean species. Up to 5 putative alleles per individual were identified, suggesting that gene duplication has occurred at this locus. Rates of non-synonymous to synonymous substitutions (ω) and maximum likelihood analyses of models of nucleotide variation provided potential evidence of ongoing positive selection at this exon. Phylogenetic analyses of DQB alleles from B. musculus and 16 other species of cetaceans revealed trans-specific conservation of MHC variants, suggesting that selection has acted on this locus over prolonged periods of time. Collectively our findings reveal that immunogenic variation in blue whales is comparable to that in terrestrial mammals, thereby providing no evidence that marine taxa are subject to reduced pathogen-induced selective pressures.     

Polymorphism of the ovine keratin-associated protein 1-4 gene (<Emphasis Type="Italic">KRTAP1-4</Emphasis>)

Keratin-associated proteins (KAPs) are one of the main structural components of the wool fibre and form a semi-rigid matrix in which the keratin intermediate filaments are embedded. Variation in the KAP genes may affect the structure of KAPs and hence wool characteristics. In this study, we used PCR-SSCP to analyse ovine KRTAP1-4 (previously B2D), a gene encoding a member of the KAP1-x family. Nine different PCR-SSCP patterns were detected in the 320 sheep that were analysed. Either one or a combination of two patterns was observed for each sheep, which was consistent with these sheep being either homozygous or heterozygous for this gene. DNA sequencing revealed that these patterns represent nine different DNA sequences. All of these sequences were unique, but shared a high homology with the published ovine KRTAP1-4 sequence, suggesting that these sequences represent allelic variants of KRTAP1-4. There were a total of 14 single nucleotide polymorphisms (SNPs) identified and these SNPs tended to be clustered in two regions. Of the 13 SNPs found in the coding region, nine were non-synonymous SNPs and would result in amino acid changes. The variation detected here may have an impact on the structure of KAP1-4 and hence affect wool traits.     

Sequence diversity of MHC class-II <Emphasis Type="BoldItalic">DRB</Emphasis> gene in gazelles (<Emphasis Type="BoldItalic">Gazella subgutturosa</Emphasis>) raised in Sanliurfa of Turkey

In this study, we aimed to assess the sequence diversity of major histocompatibility complex (MHC) class-II DRB gene at exon 2 in gazelles raised in Sanliurfa Province of Turkey. Twenty DNA samples isolated from gazelles (Gazella subgutturosa) were used for sequencing exon 2 of MHC class-II DRB gene. Target region was amplified by polymerase chain reaction (PCR) and their products were directly sequenced. Nine of these 20 samples yielded unambiguously readable sequences. Three of the nine samples were homozygotes and each showed different sequences. A 262-bp sequence obtained from the three homozygote samples were submitted to GenBank (accession numbers: KC309405, KC309406 and KC309407). Using an allele specific PCR, we detected 10 additional haplotypes. Among 13 haplotypes, 45 nucleotide positions were polymorphic and most of the polymorphic nucleotide positions localized at peptide-binding region (PBR). Rates of nonsynonymous substitutions were significantly higher than synonymous substitutions at PBR. Phylogenetic analysis of the haplotypes showed that 10 haplotypes of the gazelles were clustered together while three were clustered with ovine and bovine haplotypes. The results indicated that at least 13 haplotypes at exon 2 of MHC class-II DRB gene were showing high degree of nucleotide and amino acid diversity, and certain haplotypes of G. subgutturosa were more similar to haplotypes from sheep or cattle than to each other. Rates of synonymous and nonsynonymous substitutions suggested that positive selection was a driving force for diversity at this locus in G. subgutturosa.     

The sheep CD1 gene family contains at least four CD1B homologues

We identified four cDNA sequences encoding sheep homologues of the CD1 molecule. The sheep sequences were selected from λgt11 thymocyte cDNA libraries by hybridization with a humanCD1C probe and a homologous sheep probe. TheSCD1B-42 andSCD1A25 sequences encode complete CD1 molecules. The third sequence,SCD1B-52, which is closely related toSCD1B-42 and may be an allele, has the sequence encoding the α3 region precisely deleted. The fourth sequence,SCD1T10, is truncated at the 5′ end. All four sequences are related to the humanCD1B and domestic rabbitCD1B-like sequences at both nucleotide and amino acid level. Comparison of the derivedCD1 amino acid sequences with the sequence of major histocompatibility complex class I molecules showed that the sheep CD1 molecules, like human CD1 molecules, lack most of the conserved class I residues known to be involved in interaction with 132-microglobulin and the CD8 molecule. They do not contain the peptide docking residues involved in anchoring peptides in the peptide binding groove of class I molecules. Southern hybridization of sheep DNA with a sheepCD1 exon 4/ga3 probe showed that the sheep genome encodes at least sevenCD1 genes. The implications of these analyses for CD1 function are discussed. The nucleotide sequence data reported in this paper have been submitted to the EMBL/GenBank nucleotide sequence databases and have been assigned the accession numbers Z36890 (SCD1A25), Z36891 (SCD1B-42), Z36892 (SCD1B-52), and X90567 (SCDIT10)     

Sequence analysis of HLA class II domains: characterization of the DQw3 family of DQB genes

HLA class II allelic variants within the DQw3-related family of genes carry distinct allo-specificities and have been implicated in specific HLA-disease associations, such as insulin-dependent diabetes mellitus. To investigate the nucleotide variations which characterize DQw3 genes, we applied a novel cDNA cloning strategy that uses a single-stranded vector/primer system to facilitate DNA sequencing of allelically variable gene families. Using a DQB-specific primer sequence and M13 bacteriophage as the cloning vector, direct cloning and sequencing of multiple DQB genes was performed without the need for second strand synthesis or for subcloning. Sequence analysis from eight lymphoblastoid cell lines selected to represent different ethnic backgrounds revealed three DQw3-related DQB genes, DQB3.1, 3.2, and 3.3, corresponding to the newly designated HLA-DQw7, w8, and w9 specificities, respectively. An unusual Pro-Pro couplet at codons 55–56 is characteristic of all DQw3-positive sequences and may be contributing to the broad DQw3 allospecificity. Comparisons among ethnically disparate DQw3-related sequences showed no additional expressed or silent nucleotide substitutions among these DQB alleles. Thus, polymorphism within the DQw3 family of genes appears to be extremely limited, with a paucity of nucleotide variations accumulated by evolutionary distance.