Gene duplication and evidence for balancing selection acting on MHC class II DAA gene of the half-smooth tongue sole (Cynoglossus semilaevis)

Allelic polymorphism and evolution mechanism of major histocompatibility complex (MHC) genes has been investigated in many mammals, however, much less is known in teleost. In order to investigate the mechanisms creating and maintaining variability at the MHC class II DAA locus, we examined the polymorphism, gene duplication and balancing selection of MHC class II DAA gene of the half-smooth tongue sole (Cynoglossus semilaevis). We described 33 alleles in the C. semilaevis, recombination and gene duplication seems to play more important roles in the origin of new alleles. The rate of non-synonymous substitutions (d(N)) occurred at a significantly higher frequency than that of synonymous substitutions (d(S)) in peptide-binding region (PBR) and non-PBR, suggesting balancing selection for maintaining polymorphisms at the MHC II DAA locus. Many positive selection sites were found to act very intensively on antigen-binding sites. Our founding suggests a snapshot in an evolutionary process of MHC-DAA gene evolution of the C. semilaevis.     

Contrasting patterns of selection acting on MHC class I and class II DRB genes in the Alpine marmot (Marmota marmota)

The major histocompatibility complex (MHC) genes code for proteins that play a critical role in the immune system response. The MHC genes are among the most polymorphic genes in vertebrates, presumably due to balancing selection. The two MHC classes appear to differ in the rate of evolution, but the reasons for this variation are not well understood. Here, we investigate the level of polymorphism and the evolution of sequences that code for the peptide-binding regions of MHC class I and class II DRB genes in the Alpine marmot (Marmota marmota). We found evidence for four expressed MHC class I loci and two expressed MHC class II loci. MHC genes in marmots were characterized by low polymorphism, as one to eight alleles per putative locus were detected in 38 individuals from three French Alps populations. The generally limited degree of polymorphism, which was more pronounced in class I genes, is likely due to bottleneck the populations undergone. Additionally, gene duplication within each class might have compensated for the loss of polymorphism at particular loci. The two gene classes showed different patterns of evolution. The most polymorphic of the putative loci, Mama-DRB1, showed clear evidence of historical positive selection for amino acid replacements. However, no signal of positive selection was evident in the MHC class I genes. These contrasting patterns of sequence evolution may reflect differences in selection pressures acting on class I and class II genes.     

Allelic variation, balancing selection and positive selected sites detected from MHC class Iα gene of olive flounder

Major histocompatibility complex (MHC) genes play an important role in the immune response of vertebrates. Allelic polymorphism and pattern of evolution in MHC genes has been investigated in many mammals, however, much less is known in teleost. In the present study, we have investigated complete MHC Iα gene consists of 7 exons and 6 introns in Olive flounder (Paralichthys olivaceus). Genetic variation in the MHC class Iα gene was also tested in flounder. In 32 individuals, a total of 62 alleles were detected from exon 2 of MHC class Iα gene. The rate of non-synonymous substitutions (d ( N )) occurred at a significantly higher frequency than that of synonymous substitutions (d ( S )) in PBR and non-PBR, suggesting that balancing selection for maintaining polymorphisms at the MHC Iα locus. Many positive selection sites were found to act very intensively on antigen binding sites. Our founding suggests a snapshot in an evolutionary process of MHC Iα gene evolution of the P. olivaceus.     

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.     

Structure,diversity and evolutionary patterns of expressed MHC class II<Emphasis Type="Italic">B</Emphasis> genes in chub (<Emphasis Type="Italic">Squalius cephalus</Emphasis>), a cyprinid fish species from Europe

The polymorphism of exon 2 of the DAB genes (major histocompatibility complex [MHC] class IIB) was investigated for the first time in the freshwater cyprinid fish species, Squalius cephalus, in the wide range of its distribution in Europe. We identified 111 different MHC class IIB variants in 15 chub populations distributed from Finland to Spain. The sequence analysis showed that many structurally important amino acid sites that were conserved among tetrapods were also conserved in chub. The analysis of recombination indicated that it does not play an important role in producing and maintaining the variation of DAB genes analyzed in the present study. The exon 2 was shown to be subjected to intense positive selection. Phylogenetic analysis and sequence identities suggest the presence of two class IIB loci (DAB1-like and DAB3-like) in chub. Nevertheless, the presence of three DAB3-like sequence variants in several individuals indicates the duplication of the DAB3 gene. A contrasting selection pattern was found in DAB1-like and DAB3-like genes, which suggests the potential functional differences between these genes. Some DAB sequence variants were shared among the populations of different mtDNA lineages. The phylogenetic analyses did not confirm any biogeographical pattern of the genetic structure of MHC IIB in chub, which is in line with balancing selection and trans-species polymorphism in MHC genes. Nevertheless, cluster analysis based on the presence/absence of DAB sequence variants in the populations showed the phylogeophraphical pattern corresponding to the mtDNA lineages, which indicates that neutral selection can partially explain the MHC IIB evolution in chub.     

Sequence variability at three MHC loci of finless porpoises (<Emphasis Type="Italic">Neophocaena phocaenoides</Emphasis>)

Major histocompatibility complex (MHC) class II DQB and DRA genes and class I gene of finless porpoises (Neophocaena phocaenoides) were investigated by single-strand conformation polymorphism and sequence analysis. The DRA, DQB, and MHC-I loci each contained 5, 14, and 34 unique sequences, respectively, and considerable sequence variation was found at the MHC-I and DQB loci. Gene duplication was manifested as three to five distinct sequences at each of the DQB and MHC-I loci from some individuals, and these sequences at each of the two loci separately clustered into four groups (cluster A, B, C, and D) based on the phylogenetic trees. Phylogenetic reconstruction revealed a trans-species pattern of evolution. Relatively high rates of non-synonymous (dN) vs synonymous (dS) substitution in the peptide-binding region (PBR) suggested balancing selection for maintaining polymorphisms at the MHC-I and DQB loci. In contrast, one single locus with little sequence variation was detected in the DRA gene, and no non-synonymous substitutions in the PBR indicated no balancing selection on this gene.     

Spatial pattern of MHC class II variation in the great snipe (Gallinago media)

The genes of the major histocompatibility complex (MHC) code for proteins involved in antigen recognition and triggering of the adaptive immune response, and are therefore likely to be under selection from parasites. These selection regimes may vary in space and time. Here we report a strong geographical structure in MHC class II B genes of a migrating bird, the great snipe (Gallinago media). Genetic differentiation in the MHC between two ecologically distinct distributional regions (Scandinavian mountain populations vs. East European lowland populations) was still present after statistically controlling for the effect of selectively neutral variation (microsatellites) using partial Mantel tests. This suggests a role for selection in generating this spatial structure and that it represents local adaptation to different environments. Differentiation between populations within the two regions was negligible. Overall, we found a high number of MHC alleles (50, from 175 individuals). This, together with a tendency for a higher rate of nonsynonymous than synonymous substitutions in the peptide binding sites, and high Tajima's D in certain regions of the gene, suggests a history of balancing selection. MHC variation is often thought to be maintained by some form of balancing selection, but the nature of this selection remains unclear. Our results support the hypothesis that spatial variation in selection regimes contributes to the high polymorphism.     

A comparison of variability and population structure for major histocompatibility complex and microsatellite loci in California coastal steelhead (Oncorhynchus mykiss Walbaum)

The major histocompatibility complex (MHC) contains genes integral to immune response in vertebrates. MHC genes have been shown to be under selection in a number of vertebrate taxa, making them intriguing for population genetic studies. We have conducted a survey of genetic variation in an MHC class II gene for steelhead trout from 24 sites in coastal California and compared this variation to that observed at 16 presumably neutral microsatellite loci. A high amount of allelic variation was observed at the MHC when compared to previously published studies on other Pacific salmonids. Elevated nonsynonymous substitutions, relative to synonymous substitutions, were detected at the MHC gene, indicating the signature of historical balancing selection. The MHC data were tested for correlations to and deviations from the patterns found with the microsatellite data. Estimates of allelic richness for the MHC gene and for the microsatellites were positively correlated, as were estimates of population differentiation (F(ST)). An analysis for F(ST) outliers indicates that the MHC locus has an elevated F(ST) relative to the neutral expectation, although a significant result was found for only one particular geographical subgroup. Relatively uniform allele frequency distributions were detected in four populations, although this finding may be partially due to recent population bottlenecks. These results indicate that, at the scale studied here, drift and migration play a major role in the observed geographical variability of MHC genes in steelhead, and that contemporary selection is relatively weak and difficult to detect.     

Balancing selection,sexual selection and geographic structure in MHC genes of Great Snipe

Signatures of balancing selection are often found when investigating the extremely polymorphic regions of major histocompatibility complex (MHC) genes, and it is generally accepted that selective forces maintain this polymorphism. However, the exact nature of the selection is controversial. Theoretical studies have mainly focused on overdominance and/or frequency dependent selection while laboratory studies have emphasised the role of mate choice. Empirical field data, on the other hand, have been relatively scarce. Previously we have found that geographic structure in MHC class II genes of the Great Snipe (Gallinago media) is too pronounced to be explained by neutral forces alone. Here we test the hypothesis that sexual selection on MHC alleles may be influencing this geographic structure between mountain and lowland populations. We found evidence of balancing selection acting on MHC genes in the form of a higher rate of amino-acid changing substitutions compared to silent substitutions in the peptide binding regions. Not only natural selection but also sexual selection may influence MHC polymorphism in this bird because certain MHC alleles have been found to be associated with higher male mating success. Contrary to predictions from negative frequency dependent selection, males carrying locally rare alleles did not have a mating advantage. Instead, the mating success of alleles in a mountain population was positively correlated to their relative frequency in the mountains compared to the lowlands, implying that locally adapted MHC alleles may also be favoured by sexual selection.     

Non-neutral evolution of the major histocompatibility complex class II gene DRB1 in the sac-winged bat Saccopteryx bilineata

The immune genes of the major histocompatibility complex (MHC) are classical examples for high levels of genetic diversity and non-neutral evolution. This is particularly true for the regions containing the antigen-binding sites as, for instance, in the exon 2 of the MHC class II gene DRB. We surveyed, for the first time in the order Chiroptera, the genetic diversity within this exon in the sac-winged bat Saccopteryx bilineata. We detected 11 alleles among 85 bats, of which 79 were sampled in one population. Pairwise comparisons revealed that interallelic sequence differences ranged between 3 and 22%, although nucleotide substitutions were not evenly distributed along the exon sequence. This was most probably the result of intragenic recombination. High levels of sequence divergence and significantly more nonsynonymous than synonymous substitutions (d(N)/d(S)>1) suggest long-term balancing selection. Thus, the data are consistent with the hypothesis that recombination gives rise to new alleles at the DRB locus of the sac-winged bat, and these are maintained in the population through balancing selection. In this respect, the sac-winged bat closely resembles other mammalian species.     

中国地方鸡种MHC B-LBⅡ新等位基因的遗传多态性研究

为研究鸡MHC B-LBⅡ基因的遗传多态性,首先在8个中国地方鸡种(藏鸡、仙居鸡、北京油鸡、固始鸡、斗鸡、丝羽乌骨鸡、白耳鸡和狼山鸡)B-LBⅡ基因第二外显子扩增了一长度为 175 bp 的 DNA 片段并进行 SSCP 基因型分析;在8 个地方鸡种共 467 个个体中检测到 37 个 PCR-SSCP 基因型;从被检样品中筛选出不同基因型的个体,并在其 B-LBⅡ基因组中扩增了一个包括其第二外显子和第二内含子在内长度为374 bp的片段,通过克隆和测序获得了该片段的核苷酸序列。经序列分析,在前述地方鸡种被筛选出的 30 个无血缘关系的个体中发现了 31 个 B-LBⅡ新等位基因,并参照哺乳动物 MHC II 类 B 等位基因命名规则进行了命名。对这 31 个 B-LBⅡ新等位基因长度为 374 bp 的 DNA 片段进行比对表明,在其第二外显子序列上共有 68 个多态性变异位点,其中简约性信息位点 51 个,单变异位点 17 个,具有丰富的遗传多态性。在这些多态性变异位点中,出现在遗传密码子第一和第二位上的碱基替换率分别为 36.76% 和 35.29%。等位基因序列间的相似性估测为 90.6%-99.5%;B-LBⅡ基因第二外显子的错义替换率和同义替换率分别为 14.64±2.67%和 2.92±0.94%。结果表明,B-LBⅡ基因的丰富遗传多态性主要是由基因重组和平衡选择效应所引起的。对 B-LBⅡ等位基因第二外显子所编码的 B-LBⅡ分子β1 结构域氨基酸序列比对发现,31 个 B-LBⅡ新等位基因属于 26 个等位基因主型;在β1结构域氨基酸序列的 33个变异位点上,存在 6 个同义替换和 27 个错义替换。分析认为,那些发生在多肽结合位点上的氨基酸错义替换与鸡 MHC B-LBⅡ分子的免疫特异性有关。该结果可为鸡的抗病育种研究提供分子生物学依据。     

Historical and contemporary selection of teleost MHC genes: did we leave the past behind?

The extreme polymorphism of antigen‐presenting genes of the major histocompatibility complex (MHC) has spurred intense research unparalleled for any other gene family. This applies also to teleosts where sequence information is available for 3559 MHC class I and class II allelic variants from 137 species. This review summarizes current knowledge on the origin and maintenance of diversity at classical MHC loci. Most studies identified positive selection (i.e. elevated rates of non‐synonymous over synonymous substitutions, dN/dS) as a sign of balancing selection. A meta‐analysis on nine species with sufficient numbers of class I and class II sequences revealed that recombination rate and intensity of positive selection were positively correlated, suggesting that recombination and gene conversion played a significant role in shaping the allelic repertoire. Processes that create diversity over long timescales need to be complemented by contemporary balancing selection, either through overdominance or frequency‐dependent selection, in order to explain the high allelic diversity observed today. While some evidence for overdominance exists for a few taxa (mainly salmonids) by correlating parasite infection data or survival to MHC genotypes, field or experimental data on negative frequency‐dependent selection are lacking altogether, even though some fish species are particularly suitable as model systems. Theoretical predictions suggest that negative frequency‐dependent selection is necessary to maintain the existing polymorphism. Hence, future empirical studies should focus on detecting signals that differentiate between mechanisms of contemporary selection rather than repeatedly showing historical selection events.     

Polymorphism and signature of selection in the MHC class I genes of the three-spined stickleback Gasterosteus aculeatus

The role and intensity of positive selection maintaining the polymorphism of major histocompatibility complex (MHC) class I genes in the three-spined stickleback Gasterosteus aculeatus was investigated. The highly polymorphic set of MHC class I genes found was organized in a single linkage group. Between 5 and 14 sequence variants per individual were identified by single-stranded conformation polymorphism (SSCP) analysis. Segregation analysis studied in 10 three-spined stickleback families followed the expected pattern of Mendelian inheritance. The gamete fusion in three-spined stickleback thus seems to be random with respect to the MHC class I genes. The DNA sequence analyses showed that the expressed MHC class I loci are under strong selection pressure, possibly mediated by parasites. Codons that were revealed to be under positive selection are potentially important in antigen binding. MHC class I sequences did not form significant supported clusters within a phylogenetic tree. Analogous to MHC class II genes, it was not possible to assign the class I sequences to a specific locus, suggesting that the class I genes may have been generated by recent gene duplication.     

Population fragmentation and major histocompatibility complex variation in the spotted suslik,Spermophilus suslicus

The fragmentation of populations typically enhances depletion of genetic variation, but highly polymorphic major histocompatibility complex (MHC) genes are thought to be under balancing selection and therefore retain polymorphism despite population bottlenecks. In this study, we investigate MHC DRB (class II) exon 2 variation in 14 spotted suslik populations from two regions differing in their degree of habitat fragmentation and gene flow. We found 16 alleles that segregated in a sample of 248 individuals. The alleles were highly divergent and revealed the hallmark signs of positive selection acting on them in the past, showing a significant excess of nonsynonymous substitutions. This excess was concentrated in putative antigen‐binding sites, which suggests that past selection was driven by pathogens. MHC diversity was significantly lower in fragmented western populations than in the eastern populations, characterized by significant gene flow. In contrast to neutral variation, amova did not reveal genetic differentiation between the two regions. This may indicate similar selective pressures shaping MHC variation in both regions until the recent past. However, MHC allelic richness within a population was correlated with that for microsatellites. F_ST outlier analyses have shown that population differentiation at DRB was neither higher nor lower than expected under neutrality. The results suggest that selection on MHC is not strong enough to counteract drift that results from recent fragmentation of spotted suslik populations.     

Extensive polymorphism and geographical variation at a positively selected MHC class II B gene of the lesser kestrel (Falco naumanni)

Understanding the selective forces that shape genetic variation in natural populations remains a high priority in evolutionary biology. Genes at the major histocompatibility complex (MHC) have become excellent models for the investigation of adaptive variation and natural selection because of their crucial role in fighting off pathogens. Here we present one of the first data sets examining patterns of MHC variation in wild populations of a bird of prey, the lesser kestrel, Falco naumanni . We report extensive polymorphism at the second exon of a putatively functional MHC class II gene, Fana- DAB*1. Overall, 103 alleles were isolated from 121 individuals sampled from Spain to Kazakhstan. Bayesian inference of diversifying selection suggests that several amino acid sites may have experienced strong positive selection (ω = 4.02 per codon). The analysis also suggests a prominent role of recombination in generating and maintaining MHC diversity (ρ = 4 Nc  = 0.389 per codon, θ = 0.017 per codon). Both the Fana -DAB*1 locus and a set of eight polymorphic microsatellite markers revealed an isolation-by-distance pattern across the Western Palaearctic ( r  = 0.67; P  = 0.01 and r  = 0.50; P  = 0.04, respectively). Nonetheless, geographical variation at the MHC contrasts with relatively uniform distributions in the frequencies of microsatellite alleles. In addition, we found lower fixation rates in the MHC than those predicted by genetic drift after controlling for neutral mitochondrial sequences. Our results therefore underscore the role of balancing selection as well as spatial variations in parasite-mediated selection regimes in shaping MHC diversity when gene flow is limited.     

Evolution and trans-species polymorphism of MHC class IIbeta genes in cyprinid fish

The polymorphism of DAB genes encoding MHC IIbeta was investigated in 11 cyprinid species from central Europe. The species belonged to four subfamilies: Cyprininae, Tincinae, Gobioninae and Leuciscinae. Two paralogous groups of sequences, DAB1 and DAB3, were recognised according to the similarity of their nucleotide and amino-acid sequences and from phylogenetic analyses using either partial exon 2 or partial exon 3 sequences. A high allelic variability among species was found for exon 2, indicating extensive MHC polymorphism. Time divergence estimation supports the separation of DAB1 and DAB3 groups predating the separation into fish subfamilies, and a cyprinid origin of the DAB genes. Phylogenetic trees using exon 2 support the hypothesis of trans-species polymorphism, which appears to be limited to the subfamily level, i.e. the presence of sequences from different species in the same allelic group was more often recognised within subfamilies Cyprininae and Leuciscinae than between them. Phylogenetic trees using exon 3 reflect the phylogenetic patterns previously found for Cyprinidae systematics. Specific nucleotides and amino-acids in exon 3 that separate both subfamilies, as well as the species within the Cyprininae subfamily were observed. A lack of segregation in leuciscin species was recognised and the alleles of different leuciscin species tend to share similar motifs in exon 3. This could be explained by the ancient and complicated dispersion history of Cyprininae and the radiation of Leuciscinae. The effects of selective pressures were investigated: (1) within species, (2) among lineages, and (3) among sites. From intraspecific analyses, exon 2 sequences were identified as the targets of diversifying selection, whilst the evolution of exon 3 seems to be under the influence of purifying selection. The analyses among lineages indicate positive selection in many branches when using exon 2, therefore confirming trans-species polymorphism, whilst the DAB lineages of exon 3 are potentially submitted to purifying selection to some extent. Moreover, our results suggest the secondary acquisition of function of DAB1 group after duplication. The analyses among sites reveal that exon 2 exhibits sites under positive selection mostly corresponding to the putative PBR sites involved in the alpha-helix structure of the protein.     

Marsupial MHC class II beta genes are not orthologous to the eutherian beta gene families

The major histocompatibility complex (MHC) class II DRB, DQB, DPB, and DOB gene clusters are shared by different eutherian orders. Such an orthologous relationship is not seen between the beta genes of birds and eutherians. A high degree of uncertainty surrounds the evolutionary relationship of marsupial class II beta sequences with eutherian beta gene families. In particular, it has been suggested that marsupials utilize the DRB gene cluster. A cDNA encoding an MHC class II beta molecule was isolated from a brushtail possum mesenteric lymph node cDNA library. This clone is most similar to Macropus rufogriseus DBB. Our analysis suggests that all known marsupial beta-chain genes, excluding DMB, fall into two separate clades, which are distinct from the eutherian DRB, DQB, DPB, or DOB gene clusters. We recommend that the DAB and DBB nomenclature be reinstated. DAB and DBB orthologs are not present in eutherians. It appears that the marsupial and eutherian lineages have retained different gene clusters following gene duplication events early in mammalian evolution.     

Sequence diversity of the MHC DRB gene in the Eurasian beaver (Castor fiber)

Major histocompatibility complex (MHC) genes, coding molecules which play an important role in immune response, are the most polymorphic genes known in vertebrates. However, MHC polymorphism in some species is limited. MHC monomorphism at several MHC class I and II loci was previously reported for two neighbouring northern European populations of the Eurasian beaver (Castor fiber) and reduced selection for polymorphism has been hypothesized. Here, we analysed a partial sequence of the second exon of the MHC II DRB locus from seven relict European and Asian beaver populations. We detected 10 unique alleles among 76 beavers analysed. Only a western Siberian population was polymorphic, with four alleles detected in 10 individuals. Each of the remaining populations was fixed for a different allele. Sequences showed considerable divergence, suggesting the long persistence of allelic lineages. A significant excess of nonsynonymous substitutions was detected at the antigen binding sites, indicating that sequence evolution of beaver DRB was driven by positive selection. Current MHC monomorphism in the majority of populations may be the result of the superimposition of the recent bottleneck on pre-existing genetic structure resulting from population subdivision and differential pathogen pressure.     

Two patterns of variation among MHC class I loci in Tuatara (Sphenodon punctatus)

The genes of the major histocompatibility complex (MHC) are a central component of the immune system in vertebrates and have become important markers of functional, fitness-related genetic variation. We have investigated the evolutionary processes that generate diversity at MHC class I genes in a large population of an archaic reptile species, the tuatara (Sphenodon punctatus), found on Stephens Island, Cook Strait, New Zealand. We identified at least 2 highly polymorphic (UA type) loci and one locus (UZ) exhibiting low polymorphism. The UZ locus is characterized by low nucleotide diversity and weak balancing selection and may be either a nonclassical class I gene or a pseudogene. In contrast, the UA-type alleles have high nucleotide diversity and show evidence of balancing selection at putative peptide-binding sites. Twenty-one different UA-type genotypes were identified among 26 individuals, suggesting that the Stephens Island population has high levels of MHC class I variation. UA-type allelic diversity is generated by a mixture of point mutation and gene conversion. As has been found in birds and fish, gene conversion obscures the genealogical relationships among alleles and prevents the assignment of alleles to loci. Our results suggest that the molecular mechanisms that underpin MHC evolution in nonmammals make locus-specific amplification impossible in some species.