Polymorphism at the ovine major histocompatibility complex class II loci

Southern hybridization analysis of the ovine major histocompatibility complex (MHC) ( MhcOvar ) class II region, using sheep-specific probes for the DQA1, DQA2, DQB and DRA loci, has revealed extensive polymorphism. DQA1 and DQAP had eight and 16 alleles respectively, DQB had six and DRA had three alleles. Little information was derived from the DRB locus owing to extensive cross-hybridization between the DRB probe and the DQB locus. Differences in allele frequency between breeds were revealed. At the DQA1 locus a null allele (DQA1-N) was observed with a frequency of between 27% and 45%, making this the most common DQA1 allele in all breeds examined. The frequency of DQA1-N homozygotes was between 11% and 18%, raising questions as to the functional significance of the DQA1 gene. Linkage analysis between the DQA1, DQA2, DQB and DRA loci did not reveal any recombination.     

白鱀豚MHC基因类DQB1座位第二外元的序列变异分析

测定了 4 5个克隆的白豚 (Lipotesvexillifer)MHCⅡ类基因DQB座位第二外元 172bp的核苷酸序列 ,共获得 15种序列 ,发现了 2 2个变异位点。核苷酸的非同义替换明显多于同义替换 ,并造成了 15个氨基酸的改变。氨基酸的替换趋于集中在假定的与抗原的选择性识别相关的位点附近。白豚DQB基因的核苷酸和氨基酸序列与文献报道的白鲸 (Delphinapterusleucas)和一角鲸 (Monodonmonoceros)DQB1序列具有较高的同源性。氨基酸序列不具备人及其它一些灵长类动物DQB2基因所共有的基序 (Motif) ,而与牛DQB1基因的基序相近 ,说明本研究得到的白豚MHC序列应属于类DQB1基因。同一个体出现了多种序列的情况 ,提示白豚的DQB基因可能存在着座位重复。白豚的类DQB1座位的序列中存在多种基序的不同组合 ,推测是由于基因转换造成的.     

Evolutionary Analysis for Functional Divergence of the Toll-Like Receptor Gene Family and Altered Functional Constraints

The Toll-like receptor (TLR) gene family consists of type 1 transmembrane receptors, which play essential roles in both innate immunity and adaptive immune response by ligand recognition and signal transduction. Using all available vertebrate TLR protein sequences, we inferred the phylogenetic tree and then characterized critical amino acid residues for functional divergence by detecting altered functional constraints after gene duplications. We found that the extracellular domain of TLR genes showed higher functional divergence than that of the cytoplasmic domain, particularly in the region between leucine-rich repeat (LRR) 10 and LRR 15 of TLR 4. Our finding supports the concept that sequence evolution in the extracellular domain may be responsible for the broad diversity of TLR ligand-binding affinity, providing a testable hypothesis for potential targets that could be verified by further experimentation.     

Identification of three duplicated Spin genes in medaka (Oryzias latipes)

Gene and genomic duplications are very important and frequent events in fish evolution, and the divergence of duplicated genes in sequences and functions is a focus of research on gene evolution. Here, we report the identification and characterization of three duplicated Spindlin (Spin) genes from medaka (Oryzias latipes): OlSpinA, OlSpinB, and OlSpinC. Molecular cloning, genomic DNA Blast analysis and phylogenetic relationship analysis demonstrated that the three duplicated OlSpin genes should belong to gene duplication. Furthermore, Western blot analysis revealed significant expression differences of the three OlSpins among different tissues and during embryogenesis in medaka, and suggested that sequence and functional divergence might have occurred in evolution among them.     

Limited polymorphism at major histocompatibility complex (MHC) loci in the Swedish moose A. alces

The Swedish moose was analysed for genetic variability at major histocompatibility complex (MHC) class I and class II DQA, DQB and DRB loci using restriction fragment length polymorphism (RFLP) and single strand conformation polymorphism (SSCP) techniques. Both methods revealed limited amounts of polymorphism. Since the SSCP analysis concerned an expressed DRB gene it can be concluded that the level of functional MHC class II polymorphism, at least at the DRB locus, is low in Swedish moose. DNA fingerprinting was used to determine if the unusual pattern of low MHC variability could be explained by a low degree of genome-wide genetic diversity. Hybridizations with two minisatellite probes gave similarity indices somewhat higher than the average for other natural population, but the data suggest that the low MHC variability cannot be explained by a recent population bottleneck. However, since minisatellite sequences evolve more rapidly than MHC sequences, the low levels of MHC diversity may be attributed to a bottleneck of more ancient origin. The selection pressure for MHC variability in moose may also be reduced and we discuss the possibility that its solitary life style may reduce lateral transmission of pathogens in the population.     

The major histocompatibility complex of sheep (OLA) and two minor loci

Among 11 lymphocyte factors defined in sheep, 9 are the products of multiple alleles at 2 closely linked loci: OLA-A and OLA-B. A tenth factor is the product of a gene at a third locus: OL-X probably on the same chromosome, but in this case very distant from OLA. The last factor is the product of a gene at a fourth locus: OL-Z, independent of OLA-A and B.     

Gene conversion vs point mutation in generating variability at the antigen recognition site of major histocompatibility complex loci

We have prepared a [³²p]-labeled oligonucleotide probe carrying a ureido (-NH-CO-NH₂) function at its 3-terminus. This labeled oligomer was used to study polycondensations of urea and formaldehyde and of various phenols and formaldehyde in aqueous solution. The formation of formaldehyde copolymers attached to the amido-function of the probe was monitored by gel electrophoresis. Our results are generally in agreement with those obtained using conventional techniques. Our method is suitable for monitoring potentially prebiotic polycondensation reactions involving formaldehyde.Abbreviations HBA 4-hydroxybenzoic acid, Na⁺ salt - HBzOH 2-hydroxybenzyl alcohol - HBSA 4-hydroxybenzenesulfonic acid, Na⁺ salt Correspondence to: L.E. Orgel     

Family‐assisted inference of the genetic architecture of major histocompatibility complex variation

With their direct link to individual fitness, genes of the major histocompatibility complex (MHC) are a popular system to study the evolution of adaptive genetic diversity. However, owing to the highly dynamic evolution of the MHC region, the isolation, characterization and genotyping of MHC genes remain a major challenge. While high‐throughput sequencing technologies now provide unprecedented resolution of the high allelic diversity observed at the MHC, in many species, it remains unclear (i) how alleles are distributed among MHC loci, (ii) whether MHC loci are linked or segregate independently and (iii) how much copy number variation (CNV) can be observed for MHC genes in natural populations. Here, we show that the study of allele segregation patterns within families can provide significant insights in this context. We sequenced two MHC class I (MHC‐I) loci in 1267 European barn owls (Tyto alba), including 590 offspring from 130 families using Illumina MiSeq technology. Coupled with a high per‐individual sequencing coverage (~3000×), the study of allele segregation patterns within families provided information on three aspects of the architecture of MHC‐I variation in barn owls: (i) extensive sharing of alleles among loci, (ii) strong linkage of MHC‐I loci indicating tandem architecture and (iii) the presence of CNV in the barn owl MHC‐I. We conclude that the additional information that can be gained from high‐coverage amplicon sequencing by investigating allele segregation patterns in families not only helps improving the accuracy of MHC genotyping, but also contributes towards enhanced analyses in the context of MHC evolutionary ecology.     

Duplication and population dynamics shape historic patterns of selection and genetic variation at the major histocompatibility complex in rodents

Genetic variation at the major histocompatibility complex (MHC) is vitally important for wildlife populations to respond to pathogen threats. As natural populations can fluctuate greatly in size, a key issue concerns how population cycles and bottlenecks that could reduce genetic diversity will influence MHC genes. Using 454 sequencing, we characterized genetic diversity at the DRB Class II locus in montane voles (Microtus montanus), a North American rodent that regularly undergoes high‐amplitude fluctuations in population size. We tested for evidence of historic balancing selection, recombination, and gene duplication to identify mechanisms maintaining allelic diversity. Counter to our expectations, we found strong evidence of purifying selection acting on the DRB locus in montane voles. We speculate that the interplay between population fluctuations and gene duplication might be responsible for the weak evidence of historic balancing selection and strong evidence of purifying selection detected. To further explore this idea, we conducted a phylogenetically controlled comparative analysis across 16 rodent species with varying demographic histories and MHC duplication events (based on the maximum number of alleles detected per individual). On the basis of phylogenetic generalized linear model‐averaging, we found evidence that the estimated number of duplicated loci was positively related to allelic diversity and, surprisingly, to the strength of purifying selection at the DRB locus. Our analyses also revealed that species that had undergone population bottlenecks had lower allelic richness than stable species. This study highlights the need to consider demographic history and genetic structure alongside patterns of natural selection to understand resulting patterns of genetic variation at the MHC.     

Cloning, localization, and structure of new members of the butyrophilin gene family in the juxta-telomeric region of the major histocompatibility complex

 New members of the butyrophilin (BT) gene family have been identified by cDNA and genomic cloning. Six genes are described: BT2.1, 2.2, 2.3, and BT3.1, 3.2, and 3.3. BT2, BT3, and BT form three distinct subfamilies sharing about 95% amino acid identity at the intra subfamily level and 50% identity at the interfamily level. All the BT2 and BT3 subfamily members map close to BT in the juxta-telomeric region of the major histocompatibility complex. The BT2 members have the canonical structural organization of BT, i.e., two immunoglobulin domains followed by a transmembrane anchor and a B30.2 intracytoplasmic domain. In the BT3 subfamily, only BT3.3 has the structural organization of BT. The two other genes, BT3.1 and BT3.2, code for putative protein without the B30.2 domain. In the case of BT3.2, this is due to an Alu insertion in the B30.2 coding exon, leading to a new polyadenylation site. Received: 21 April 1997 / Revised: 13 June 1997     

Parallel signatures of selection at genomic islands of divergence and the major histocompatibility complex in ecotypes of sockeye salmon across Alaska

Understanding the genetic mechanisms that facilitate adaptive radiation is an important component of evolutionary biology. Here, we genotyped 82 neutral SNPs, seven SNPs in islands of divergence identified in a previous study (island SNPs), and a region of the major histocompatibility complex (MHC) in 32 populations of sockeye salmon to investigate whether conserved genes and genomic regions are involved in adaptive radiation. Populations representing three ecotypes were sampled from seven drainages with differing habitats and colonization histories spanning a range of 2,000 km. We found strong signatures of parallel selection across drainages at the island SNPs and MHC, suggesting that the same loci undergo divergent selection during adaptive radiation. However, patterns of differentiation at most island SNPs and the MHC were not associated with ecotypes, suggesting that these loci are responding differently to a mosaic of selective pressures. Our study provides some of the first evidence that conserved genomic islands may be involved in adaptive divergence of salmon populations. Additionally, our data provide further support for the hypothesis that sockeye salmon inhabiting rivers unconnected to lakes harbour similar genetic diversity across large distances, are likely the ancestral form of the species, and have repeatedly recolonized lake systems as they have become available after glacial recession. Finally, our results highlight the value and importance of validating outlier loci by screening additional populations and regions, a practice that will hopefully become more common in the future.     

Gene conversion of major histocompatibility complex genes is associated with CpG-rich regions

 We examined 32 DNA sequences of mouse and human major histocompatibility complex (MHC) genes believed to have been subjected to gene conversion events. All regions of the mouse H2 genes as well as the human HLA genes which have been implied to be involved in gene conversion events had elevated levels of CpG dinucleotides, whereas the rest of the genes showed extensive CpG suppression. Mouse MHC genes which have been suspected but not directly implied to be involved in gene conversion events also showed elevated levels of CpG dinucleotides. Moreover, both mouse and human MHC genes which have never been suspected of undergoing gene conversion had low levels of CpG throughout the genes. These results indicate that high CpG levels are correlated with gene conversion rather than with polymorphism, as non-polymorphic genes that have been implicated as gene conversion donors also have elevated levels of CpG dimers in the involved regions, whereas polymorphic genes which have never been considered to undergo gene conversion events have a low level of CpG dinucleotides. We also studied the methylation pattern of CpG dimers in the Abk gene by restriction enzyme digestion of mouse testis DNA followed by Southern blot and hybridization to an Abk-specific probe. The examined CpG dimers in prepubescent mice, where the latest germline stages are spermatogonia, leptene, or pachytene, are respectively non-methylated. Accordingly, the CpG dimers appear to be non-methylated in germline DNA from the testis of prepubescent mice, where gene conversions have been reported to occur. Received: 6 May 1998 / Revised: 2 September 1998     

Diversity at two genetic loci associated with the major histocompatibility complex in the golden snub-nosed monkey (Rhinopithecus roxellana)

Habitat loss via human activity has fragmented populations of the golden snub-nosed monkey (Rhinopithecus roxellana), and thus affected patterns of gene flow. We investigated in-depth a single troop in the Qinling Mountains, central China, two major histocompatibility complex (MHC) class II loci, DQA1 and DQB1, and compared the resulting data with data from troops from the wider Qinling Mountains region. We found that a novel DQB1 allele was only present in the study troop and relatively few divergent alleles at the DQA1 and DQB1 loci were present compared with the wider population. The inbreeding coefficient (F_is) at the MHC region was lower than previous measurements, which may have reflected different sampling strategies. However, R. roxellana has relatively high diversity in MHC genes, even though it has probably experienced serious past population bottlenecks and reduced gene flow between populations. We also found that some alleles present in the wider population had been lost in the study troop, and suggest that conservation management strategies be implemented to increase gene flow between troops in order to increase genetic variation.     

Balancing selection and heterozygote advantage in major histocompatibility complex loci of the bottlenecked Finnish wolf population

Maintaining effective immune response is an essential factor in the survival of small populations. One of the most important immune gene regions is the highly polymorphic major histocompatibility complex (MHC). We investigated how a population bottleneck and recovery have influenced the diversity and selection in three MHC class II loci, DLA‐DRB1, DLA‐DQA1 and DLA‐DQB1, in the Finnish wolf population. We studied the larger Russian Karelian wolf population for comparison and used 17 microsatellite markers as reference loci. The Finnish and Karelian wolf populations did not differ substantially in their MHC diversities ( = 0.047, P = 0.377), but differed in neutral microsatellite diversities ( = 0.148, P = 0.008). MHC allele frequency distributions in the Finnish population were more even than expected under neutrality, implying balancing selection. In addition, an excess of nonsynonymous compared to synonymous polymorphisms indicated historical balancing selection. We also studied association between helminth (Trichinella spp. and Echinococcus canadensis) prevalence and MHC diversity at allele and SNP level. MHC‐heterozygous wolves were less often infected by Trichinella spp. and carriers of specific MHC alleles, SNP haplotypes and SNP alleles had less helminth infections. The associated SNP haplotypes and alleles were shared by different MHC alleles, which emphasizes the necessity of single‐nucleotide‐level association studies also in MHC. Here, we show that strong balancing selection has had similar effect on MHC diversities in the Finnish and Russian Karelian wolf populations despite significant genetic differentiation at neutral markers and small population size in the Finnish population.     

Convergent evolution of major histocompatibility complex molecules in humans and New World monkeys

  In both Old World and New World monkeys Mhc-DRB sequences have been found which resemble human DRB1*03 and DRB3 genes in their second exon. The resemblance is shared sequence motifs and clustering of the genes or the encoded proteins in phylogenetic trees. This similarity could be due to common ancestry, convergence at the molecular level, or chance. To test which of these three explanations applies, we sequenced segments of New World monkey and macaque genes which encompass the entire second exon and large parts of both flanking introns. The test strongly supports the monophyly of New World monkey DRB intron sequences. The phylogenies of introns 1 and 2 from DRB1*03-like and DRB3-like genes are congruent, but both are incongruent with the exon 2-based phylogeny. The matching of intron 1- and intron 2-based phylogenies with each other suggests that reciprocal recombination has not played a major role in exon 2 evolution. Statistical comparisons of exon 2 from different DRB1*03 and DRB3 lineages indicate that it was neither gene conversion (descent), nor chance, but molecular convergence that has shaped their characteristic motifs. The demonstration of convergence in anthropoid Mhc-DRB genes has implications for the classification, age, and mechanism of generation of DRB allelic lineages. Received: 30 August 1999 / Revised: 19 October 1999     

Sequence variation at the major histocompatibility complex DRB loci in beluga (Delphinapterus leucas) and narwhal (Monodon monoceros)

 The variation at loci with similarity to DRB class II major histocompatibility complex loci was assessed in 313 beluga collected from 13 sampling locations across North America, and 11 narwhal collected in the Canadian high Arctic. Variation was assessed by amplification of exon 2, which codes for the peptide binding region, via the polymerase chain reaction, followed by either cloning and DNA sequencing or single-stranded conformation polymorphism analysis. Two DRB loci were identified in beluga: DRB1, a polymorphic locus, and, DRB2, a monomorphic locus. Eight alleles representing five distinct lineages (based on sequence similarity) were found at the beluga DRB1 locus. Although the relative number of alleles is low when compared with terrestrial mammals, the amino acid variation found among the lineages is moderate. At the DRB1 locus, the average number of nonsynonymous substitutions per site is greater than the average number of synonymous substitutions per site (0.0806 : 0.0207, respectively;P<0.01). Most of the 31 amino acid substitutions do not conserve the physiochemical properties of the residue, and 21 of these are located at positions implicated as forming pockets responsible for the selective binding of foreign peptide side chains. Only DRB1 variation was examined in 11 narwhal, revealing a low amount of variation. These data are consistent with an important role for the DRB1 locus in the cellular immune response of beluga. In addition, the ratio of nonsynonymous to synonymous substitutions is similar to that among primate alleles, arguing against a reduction in the balancing selection pressure in the marine environment. Two hypotheses may explain the modest amount of Mhc variation when compared with terrestrial mammals: small population sizes at speciation or a reduced neutral substitution rate in cetaceans. Received: 15 July 1997 / Revised: 24 March 1998     

Comparative analysis of population structure across environments and geographical scales at major histocompatibility complex and microsatellite loci in Atlantic salmon (Salmo salar)

Evidence of selection acting on major histocompatibility complex (MHC) genes has been illustrated with the analysis of their nucleotide sequences and allele frequency distribution. Comparing the patterns of population differentiation at neutral markers and MHC genes in the wild may provide further insights about the relative role of selection and neutrality in shaping their diversity. In this study, we combine both methods to assess the role of selection on a MHC gene in Atlantic salmon. We compare variation at a MHC class II B locus and microsatellites among 14 samples from seven different rivers and seven subpopulations within a single river system covering a variety of habitats and different geographical scales. We show that diversifying selection is acting on the sites involved in antigen presentation and that balancing selection maintains a high level of polymorphism within populations. Despite important differences in habitat type, the comparison of the population structure at MHC and microsatellites on large geographical scales reveals a correlation between patterns of differentiation, indicating that drift and migration have been more important than selection in shaping population differentiation at the MHC locus. In contrast, strong discrepancies between patterns of population differentiation at the two types of markers provides support for the role of selection in shaping population structure within rivers. Together, these results confirm that natural selection is influencing MHC gene diversity in wild Atlantic salmon although neutral forces may also be important in their evolution.     

Evolutionary analysis for functional divergence of Jak protein kinase domains and tissue-specific genes

Jak (Janus kinase) is a nonreceptor tyrosine kinase, which plays important roles in signal transduction pathways. The unique feature of Jak is that, in addition to a fully functional tyrosine kinase domain (JH1), Jak possesses a pseudokinase domain (JH2). Although JH2 lost its catalytic function, experimental evidence has shown that this domain may have acquired some new but unknown functions. This apparent functional divergence after the (internal) domain duplication may result in dramatic changes of selective constraints at some sites. We conducted a data analysis to test this hypothesis. Our result shows that shifted selective constraints (or shifted evolutionary rates) between the JH1 and the JH2 domains are statistically significant. Predicted amino acid sites by posterior analysis can be classified into two groups: very conserved in JH1 but highly variable in JH2, and vice versa. Moreover, we have studied the evolutionary pattern of four tissue-specific genes, Jak1, Jak2, Jak3, and Tyk2, which were generated in the early stages of vertebrates. We found that after the (first) gene duplication, site-specific rate shifts between Jak2/Jak3 and Jak1/Tyk are significant, presumably as a consequence of functional divergence among these genes. The implication of our study for functional genomics is discussed.     

Patterns of nucleotide substitutions inferred from the phylogenies of the class I major histocompatibility complex genes

Summary Patterns of nucleotide substitutions in human major histocompatibility complex (MHC) class I genes were estimated by using phylogenetic trees of DNA sequences. The pattern is defined as a set of 12 parameters, each of which represents the relative frequency of substitutions from a particular nucleotide to another. The pattern at the antigen recognition sites (ARS) in functional MHC genes was remarkably different from that at the remaining coding region (non-ARS). In particular, the proportion of transitions among all the nucleotide substitutions (P _s) was extremely low at the third codon positions of ARS. In the HLA-A genes, P _s at the third codon positions was only 6% in ARS, whereas it was 69% in non-ARS. In HLA-B, the corresponding values were 30% in ARS and 80% in non-ARS, respectively. On the other hand, P _s in a class I pseudogene (HLA-H) was 57%, which was in good agreement with P _s in other pseudogenes. Because pseudogenes are selectively neutral, the pattern in pseudogenes is regarded as the pattern of spontaneous substitution mutations. In general, the pattern in functional genes that are subject to selective forces deviates from the pattern in pseudogenes. At the third codon positions in coding regions, transitions scarcely cause amino acid replacements, whereas about half of transversions do cause replacements. Accordingly, P _s at the third codon positions decreases if amino acid replacements are accelerated by natural selection but increases if amino acids are conserved by functional constraint. Our observations imply that the ARS region is subject to natural selection favoring amino acid replacements, whereas the non-ARS region is subject to functional constraint. Offprint requests to: T. Gojobori