Diversity and evolution of the highly polymorphic tandem repeat LEI0258 in the chicken MHC-B region

The chicken major histocompatibility complex (MHC) is located on the microchromosome 16 and is described as the most variable region in the genome. The genes of the MHC play a central role in the immune system. Particularly, genes encoding proteins involved in the antigen presentation to T cells. Therefore, describing the genetic polymorphism of this region is crucial in understanding host–pathogen interactions. The tandem repeat LEI0258 is located within the core area of the B region of the chicken MHC (MHC-B region) and its genotypes correlate with serology. This marker was used to provide a picture of the worldwide diversity of the chicken MHC-B region and to categorize chicken MHC haplotypes. More than 1,600 animals from 80 different populations or lines of chickens from Africa, Asia, and Europe, including wild fowl species, were genotyped at the LEI0258 locus. Fifty novel alleles were described after sequencing. The resulting 79 alleles were classified into 12 clusters, based on the SNPs and indels found within the sequences flanking the repeats. Furthermore, hypotheses were formulated on the evolutionary dynamics of the region. This study constitutes the largest variability report for the chicken MHC and establishes a framework for future diversity or association studies.     

A Naturally Occurring Variant of Porcine Mx1 Associated with Increased Susceptibility to Influenza Virus In Vitro

Mx1 has been implicated in resistance to the influenza virus. We have now identified four alleles of the Mxl gene in domesticated breeds of pigs. Two of the alleles encode deletion variants (a 3-bp deletion in exon 13 and an 11-bp deletion in exon 14), which might be expected to interfere with Mx activity. The porcine Mxl genes corresponding to wild type, the 3-bp deletion mutant, and the 11-bp deletion mutant were cloned and expressed in NIH3T3 cells, and the antiviral activity for influenza virus was assayed. Virus yield was observed to be 10–100-fold greater with the 11-bp deletion allele than that for wild type and the 3-bp deletion alleles. The results suggest that the 11-bp deletion type is lacking antiviral activity able to contribute to the interference of influenza virus replication.     

MHC class I variation in a natural blue tit population (Cyanistes caeruleus)

The major histocompatibility complex (MHC) is central to the vertebrate immune system and its highly polymorphic genes are considered to influence several life-history traits of individuals. To characterize the MHC in a natural population of blue tits (Cyanistes caeruleus) we investigated the class I exon 3 diversity of more than 900 individuals. We designed two pairs of motif-specific primers that reliably amplify independent subsets of MHC alleles. Applying denaturing gradient gel electrophoresis (DGGE) we obtained 48 independently inherited units of unique band patterns (DGGE-haplogroups), which were validated in a segregation analysis within 105 families. In a second approach, we extensively sequenced 6 unrelated individuals to confirm that DGGE-haplogroup composition reflects individual allelic variation. The highest number of different DGGE-haplogroups in a single individual corresponded in 19 MHC exon 3 sequences, suggesting a minimum of 10 amplified MHC class I loci in the blue tit. In total, we identified 50 unique functional and 3 non-functional sequences. Functional sequences showed high levels of recombination and strong positive selection in the antigen binding region, whereas nucleotide diversity was comparatively low in the range of all passerine species. Finally, in a phylogenetic comparison of passerine MHC class I exon 3 sequences we discuss conflicting evolutionary signals possibly due to recent gene duplication, recombination events and concerted evolution. Our results indicate that the described method is suitable to effectively explore the MHC diversity and its ecological impacts in blue tits in future studies.     

Characterization of MHC class II B polymorphism in bottlenecked New Zealand saddlebacks reveals low levels of genetic diversity

The major histocompatibility complex (MHC) is integral to the vertebrate adaptive immune system. Characterizing diversity at functional MHC genes is invaluable for elucidating patterns of adaptive variation in wild populations, and is particularly interesting in species of conservation concern, which may suffer from reduced genetic diversity and compromised disease resilience. Here, we use next generation sequencing to investigate MHC class II B (MHCIIB) diversity in two sister taxa of New Zealand birds: South Island saddleback (SIS), Philesturnus carunculatus, and North Island saddleback (NIS), Philesturnus rufusater. These two species represent a passerine family outside the more extensively studied Passerida infraorder, and both have experienced historic bottlenecks. We examined exon 2 sequence data from populations that represent the majority of genetic diversity remaining in each species. A high level of locus co-amplification was detected, with from 1 to 4 and 3 to 12 putative alleles per individual for South and North Island birds, respectively. We found strong evidence for historic balancing selection in peptide-binding regions of putative alleles, and we identified a cluster combining non-classical loci and pseudogene sequences from both species, although no sequences were shared between the species. Fewer total alleles and fewer alleles per bird in SIS may be a consequence of their more severe bottleneck history; however, overall nucleotide diversity was similar between the species. Our characterization of MHCIIB diversity in two closely related species of New Zealand saddlebacks provides an important step in understanding the mechanisms shaping MHC diversity in wild, bottlenecked populations.     

Selection, trans-species polymorphism, and locus identification of major histocompatibility complex class IIβ alleles of New World ranid frogs

Genes encoded by the major histocompatibility complex (MHC) play key roles in the vertebrate immune system. However, our understanding of the evolutionary processes and underlying genetic mechanisms shaping these genes is limited in many taxa, including amphibians, a group currently impacted by emerging infectious diseases. To further elucidate the evolution of the MHC in frogs (anurans) and develop tools for population genetics, we surveyed allelic diversity of the MHC class II β1 domain in both genomic and complementary DNA of seven New World species in the genus Rana (Lithobates). To assign locus affiliation to our alleles, we used a “gene walking” technique to obtain intron 2 sequences that flanked MHC class IIβ exon 2. Two distinct intron sequences were recovered, suggesting the presence of at least two class IIβ loci in Rana. We designed a primer pair that successfully amplified an orthologous locus from all seven Rana species. In total, we recovered 13 alleles and documented trans-species polymorphism for four of the alleles. We also found quantitative evidence of selection acting on amino acid residues that are putatively involved in peptide binding and structural stability of the β1 domain of anurans. Our results indicated that primer mismatch can result in polymerase chain reaction (PCR) bias, which influences the number of alleles that are recovered. Using a single locus may minimize PCR bias caused by primer mismatch, and the gene walking technique was an effective approach for generating single-copy orthologous markers necessary for future studies of MHC allelic variation in natural amphibian populations.     

Frequent segmental sequence exchanges and rapid gene duplication characterize the<Emphasis Type="Italic"> MHC</Emphasis> class I genes in lemurs

Major histocompatibility complex (MHC) class I genes have complicated and profound evolutionary histories. To reconstruct and better understand their histories, partial class I genes (exon 2–intron 2–exon 3) were sequenced in a sampling of prosimians (Strepsirhini, Primates). In total, we detected 117 different sequences from 36 Malagasy prosimians (lemurs) and 1 non-Malagasy prosimian (galago) representing 4 families, 7 genera, and 13 species. Unlike the MHC class II genes (MHC-DRB), MHC class I genes show a generally genus-specific mode of evolution in lemurs. Additionally, no prosimian class I loci were found to be orthologous to HLA genes, even at highly conserved loci (such as HLA-E, HLA-F). Phylogenetic analysis indicates that nucleotide diversity among loci was very small and the persistence time of the polymorphisms was short, suggesting that the origin of the lemur MHC class I genes detected in this study was relatively recent. The evolutionary mode of these genes is similar to that of classic HLA genes, HLA-A, HLA-B, and HLA-C, in terms of their recent origin and rarity of pseudogenes, and differs from them with respect to the degree of gene duplications. From the viewpoint of MHC genes evolution, some interlocus sequence exchanges were apparently observed in the lemur lineage upon phylogenetic and amino acid motif analyses. This is also in contrast to the evolutionary mode of HLA genes, where intralocus exchanges have certainly occurred but few interlocus exchanges have taken place. Consequently, the gene conversion model for explaining the generation of the MHC diversity among different loci can be thought to play more important roles in the evolution of lemur MHC class I genes than in that of HLA genes.Electronic Supplementary Material Supplementary material is available in the online version of this article at     

Detection of a common BOLA-DRB3 deletion by sequence-specific oligonucleotide typing

The bovine MHC class II BoLA-DRB3*2A has an amino acid deletion of unknown function at codon 65 in the second exon, which codes for the antigen-binding site. Sequence-specific oligonucleotides were designed based on published nucleotide sequences on BoLA-DRB3 alleles, and used to detect this deletion in 51 Hereford cattle. Probes 65+ and 65- detect the presence or absence of codon 65 respectively. Oligonucleotide probes were labelled with Digoxigenin (DIG), hybridized to dot blots of BOLA-DRB3 exon 2 polymerase chain reaction (PCR) product, and detected by chemiluminescence. Of the 51 animals screened, two were homozygous and 11 were heterozygous for the deletion at codon 65. The methodology described here provides the necessary tools to screen rapidly for this deletion in a large number of animals in order to study its effect on antigen binding and immune response.     

兴国红鲤MHCⅡ类基因多态性、表达及其与鱼体抗病力关系的分析

对41尾感病和22尾抗病兴国红鲤的308个有效克隆进行测序,获得171条不同的MHC Ⅱ类基因编码序列,分属26个不同的等位基因,其中Cyca-DXA24Cyca-DXA36为新发现的13个等位基因。MHCⅡ类基因片段的长度为624 bp,包括第14个外显子,分别编码信号肽、1和2结构域及连接肽/跨膜区。1结构域的变异明显大于2结构域,表现在1结构域中核苷酸和氨基酸变异位点比例（55.16%和79.76%）明显高于2结构域的变异位点比例（45.96%和68.42%）。1结构域的PBR区的非同义碱基替换率（dN）与同义碱基替换率（dS）的比值 （=dN/dS）为5.742,远远高于非抗原结合位点（non-PBR）及2结构域的0.755、0.592,揭示兴国红鲤MHC Ⅱ类基因的1结构域在进化过程中受到正向选择作用。等位基因Cyca-DXA24 （P0.01）与兴国红鲤对嗜水气单胞菌的抗性相关,等位基因Cyca-DXA3 （P0.05）、Cyca-DXA4 （P0.01）、Cyca-DXA6 （P0.05）、Cyca-DXA33 （P0.05）与兴国红鲤对嗜水气单胞菌的易感性相关。荧光定量PCR结果表明,MHCⅡ类基因在健康兴国红鲤的肾、肝、鳃等10个组织均能普遍表达。人工感染嗜水气单胞菌后,肾、肝、脾3个组织中的MHCⅡ类基因的表达量均发生了不同程度的变化,表明MHCⅡ类分子在兴国红鲤的免疫反应中起到重要作用。       

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.     

The genetic constitutions of complementary genes Pp and Pb determine the purple color variation in pericarps with cyanidin-3-O-glucoside depositions in black rice

The purple pericarp color in rice was controlled by two dominant complementary genes, Pb and Pp. Crossing black rice ‘Heugnambyeo’ variants with three varieties of white pericarp rice gave a segregation ratio of 9 purple: 3 brown: 4 white. The Pp genes were segregated by homozygous PpPp alleles for the dark purple pericarps, heterozygous Pppp alleles for the medium and mixed purple pericarps, and homozygous pppp alleles for either brown or white pericarps with a 1 PpPp: 2 Pppp: 1 pppp segregation ratio, indicating that the Pp allele in rice is incompletely dominant to the recessive pp allele. Among the purple seeds, the amount of cyanidin-3-O-glucoside was higher in the dark purple seeds (Pb_PpPp) than in the medium purple seeds (Pb_Pppp). Moreover, no cyanidin-3-glucoside was detected in brown (Pb_pppp) or white pericarp seeds (pbpbpppp). These findings indicated that the level of cyanidin-3-glucoside was determined by the copy number of the Pp allele. Further genotype investigation of the F₃ progeny demonstrated that the dominant Pb allele was present in either purple or brown pericarp. A 2-bp (GT) deletion from the DNA sequences of the dominant and functional Pb was found in the same DNA sequences of the recessive and non-functional pb allele. These findings suggested that the presence of at least a dominant Pb allele was an essential factor for color development in rice pericarps. In conclusion, the Pp allele in rice is incompletely dominant to the recessive pp allele; thus, the number of dominant Pp alleles determines the concentration of cyanidin-3-O-glucoside in black rice.     

A new major histocompatibility complex class I b gene expressed in the mouse blastocyst and placenta

 Because of the role major histocompatibility complex (MHC) class I b molecules may play during mouse embryonic development, we thought it would be interesting to search for additional MHC class I b molecules that might be expressed in preimplantation embryos, and in particular in the trophoblastic lineage. We therefore screened a mouse preimplantation blastocyst cDNA library for MHC class I sequences. This search led to the identification and characterization of a new MHC class I b gene, blastocyst MHC. Sequences identical to the exons and 3′ untranslated region of this gene have been found in many laboratory mouse strains, as well as in the related mouse species Mus spreciligus. The presence of this gene in mouse strains of different MHC class I haplotypes argues that blastocyst MHC is a unique, newly-described gene rather than a new allele of a previously described mouse MHC class I gene. Blastocyst MHC has the structure of an MHC class I b gene, with the six exons characteristic of T-region genes. It is linked to H2-D. The amino acid sequence encoded by this gene maintains all the features of a functional antigen-presentation domain. The blastocyst MHC gene, like the human class I b gene HLA-G, is expressed at the blastocyst stage and in the placenta, and may be the mouse analog for HLA-G. Received: 31 May 1996 / Revised: 19 August 1996     

MHC class I typing in a songbird with numerous loci and high polymorphism using motif-specific PCR and DGGE

The major histocompatibility complex (MHC) has a central role in the specific immune defence of vertebrates. Exon 3 of MHC class I genes encodes the domain that binds and presents peptides from pathogens that trigger immune reactions. Here we develop a fast population screening method for detecting genetic variation in the MHC class I genes of birds. We found evidence of at least 15 exon 3 sequences in the investigated great reed warbler individual. The organisation of the great reed warbler MHC class I genes suggested that a locus-specific screening protocol is impractical due to the high similarity between alleles across loci, including the introns flanking exon 3. Therefore, we used motif-specific PCR to amplify two subsets of alleles (exon 3 sequences) that were separated with by DGGE. The motif-specific primers amplify a substantial proportion of the transcribed class I alleles (2-12 alleles per individual) from as many as six class I loci. Although not exhaustive, this gives a reliable estimate of the class I variation. The method is highly repeatable and more sensitive in detecting genetic variation than the RFLP method. The motif-specific primers also allow us to avoid screening pseudogenes. In our study population of great reed warblers, we found a high level of genetic variation in MHC class I, and no less than 234 DGGE genotypes were detected among 248 screened individuals.     

Copy Number Variants in the Kallikrein Gene Cluster

The kallikrein gene family (KLK1-KLK15) is the largest contiguous group of protease genes within the human genome and is associated with both risk and outcome of cancer and other diseases. We searched for copy number variants in all KLK genes using quantitative PCR analysis and analysis of inheritance patterns of single nucleotide polymorphisms. Two deletions were identified: one 2235-bp deletion in KLK9 present in 1.2% of alleles, and one 3394-bp deletion in KLK15 present in 4.0% of alleles. Each deletion eliminated one complete exon and created out-of-frame coding that eliminated the catalytic triad of the resulting truncated gene product, which therefore likely is a non-functional protein. Deletion breakpoints identified by DNA sequencing located the KLK9 deletion breakpoint to a long interspersed element (LINE) repeated sequence, while the deletion in KLK15 is located in a single copy sequence. To search for an association between each deletion and risk of prostate cancer (PC), we analyzed a cohort of 667 biopsied men (266 PC cases and 401 men with no evidence of PC at biopsy) using short deletion-specific PCR assays. There was no association between evidence of PC in this cohort and the presence of either gene deletion. Haplotyping revealed a single origin of each deletion, with most recent common ancestor estimates of 3000-8000 and 6000-14 000 years for the deletions in KLK9 and KLK15, respectively. The presence of the deletions on the same haplotypes in 1000 Genomes data of both European and African populations indicate an early origin of both deletions. The old age in combination with homozygous presence of loss-of-function variants suggests that some kallikrein-related peptidases have non-essential functions.     

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.     

Current perspectives on the intensity of natural selection of MHC loci

Polymorphism of genes in the major histocompatibility complex (MHC) is believed to be maintained by balancing selection. However, direct evidence of selection has proven difficult to demonstrate. In 1994, Satta and colleagues estimated the selection intensity of the human MHC (human leukocyte antigen (HLA)) loci; however, at that time the number of HLA sequences was limited. By comparing five different methods, this study demonstrated the best way to calculate the selection coefficient, through a computer simulation study. Since the study, many HLA nucleotide sequences have been made available. Our new analysis takes advantage of these newly available sequences and compares new estimates with those of the previous study. Generally, our new results are consistent with those of the 1994 study. Our results show that, even after 20 years of exhaustive sequencing of human HLA, the number of dominant HLA alleles, on which our original estimate of selection intensity depended, appears to be conserved. Indeed, according to the frequency distribution for each HLA allele, most sequences in the database were minor or private alleles; therefore, we conclude that the selection intensities of HLA loci are at most 4.4 % even though the HLA is the prominent example on which the natural selection has been operating.     

Resequencing Reveals Different Domestication Rate for BADH1 and BADH2 in Rice (Oryza sativa)

BADH1 and BADH2 are two homologous genes, encoding betaine aldehyde dehydrogenase in rice. In the present study, we scanned BADHs sequences of 295 rice cultivars, and 10 wild rice accessions to determine the polymorphisms, gene functions and domestication of these two genes. A total of 16 alleles for BADH1 and 10 alleles for BADH2 were detected in transcribed region of cultivars and wild species. Association study showed that BADH1 has significant correlation with salt tolerance in rice during germination stage, the SNP (T/A) in exon 4 is highly correlated with salt tolerance index (STI) (P<10⁻⁴). While, BADH2 was only responsible for rice fragrance, of which two BADH2 alleles (8 bp deletion in exon 7 and C/T SNP in exon 13) explain 97% of aroma variation in our germplasm. Theses indicate that there are no overlapping functions between the two homologous genes. In addition, a large LD block was detected in BADH2 region, however, there was no large LD blocks in a 4-Mb region of BADH1. We found that BADH2 region only showed significant bias in Tajima’s D value from the balance. Extended haplotype homozygosity study revealed fragrant accessions had a large LD block that extended around the mutation site (8 bp deletion in exon 7) of BADH2, while both of the BADH1 alleles (T/A in exon 4) did not show large extended LD block. All these results suggested that BADH2 was domesticated during rice evolution, while BADH1 was not selected by human beings.     

Diversity of Extended HLA-DRB1 Haplotypes in the Finnish Population

The Major Histocompatibility Complex (MHC, 6p21) codes for traditional HLA and other host response related genes. The polymorphic HLA-DRB1 gene in MHC Class II has been associated with several complex diseases. In this study we focus on MHC haplotype structures in the Finnish population. We explore the variability of extended HLA-DRB1 haplotypes in relation to the other traditional HLA genes and a selected group of MHC class III genes. A total of 150 healthy Finnish individuals were included in the study. Subjects were genotyped for HLA alleles (HLA-A, -B, -DRB1, -DQB1, and -DPB1). The polymorphism of TNF, LTA, C4, BTNL2 and HLA-DRA genes was studied with 74 SNPs (single nucleotide polymorphism). The C4A and C4B gene copy numbers and a 2-bp silencing insertion at exon 29 in C4A gene were analysed with quantitative genomic realtime-PCR. The allele frequencies for each locus were calculated and haplotypes were constructed using both the traditional HLA alleles and SNP blocks. The most frequent Finnish A∼B∼DR -haplotype, uncommon in elsewhere in Europe, was A*03∼B*35∼DRB1*01∶01. The second most common haplotype was a common European ancestral haplotype AH 8.1 (A*01∼B*08∼DRB1*03∶01). Extended haplotypes containing HLA-B, TNF block, C4 and HLA-DPB1 strongly increased the number of HLA-DRB1 haplotypes showing variability in the extended HLA-DRB1 haplotype structures. On the contrary, BTNL2 block and HLA-DQB1 were more conserved showing linkage with the HLA-DRB1 alleles. We show that the use of HLA-DRB1 haplotypes rather than single HLA-DRB1 alleles is advantageous when studying the polymorphisms and LD patters of the MHC region. For disease association studies the HLA-DRB1 haplotypes with various MHC markers allows us to cluster haplotypes with functionally important gene variants such as C4 deficiency and cytokines TNF and LTA, and provides hypotheses for further assessment. Our study corroborates the importance of studying population-specific MHC haplotypes.     

Molecular characterization,balancing selection,and genomic organization of the tree shrew (Tupaia belangeri) MHC class I gene

The major histocompatibility complex (MHC) class I genes play a pivotal role in the adaptive immune response among vertebrates. Accordingly, in numerous mammals the genomic structure and molecular characterization of MHC class I genes have been thoroughly investigated. To date, however, little is known about these genes in tree shrews, despite the increasingly popularity of its usage as an animal model. To address this deficiency, we analyzed the structure and characteristic of the tree shrew MHC class I genes (Tube-MHC I) and performed a comparative gene analysis of the tree shrew and other mammal species. We found that the full-length cDNA sequence of the tree shrew MHC class I is 1074 bp in length. The deduced peptide is composed of 357 amino acids containing a leader peptide, an α1 and α2 domain, an α3 domain, a transmembrane domain and a cytoplasmic domain. Among these peptides, the cysteines, CD8⁺ interaction and N-glycosylation sites are all well conserved. Furthermore, the genomic sequence of the tree shrew MHC class I gene was identified to be 3180 bp in length, containing 8 exons and 7 introns. In 21 MHC class I sequences, we conducted an extensive study of nucleotide substitutions. The results indicated that in the peptide binding region (PBR) the rate of non-synonymous substitutions (dN) to synonymous substitutions (dS) was greater than 1, suggesting balancing selection at the PBR. These findings provide valuable contributions in furthering our understanding of the structure, molecular polymorphism, and function of the MHC class I genes in tree shrews, further improving their utility as an animal model in biomedical research.