Mhc class I genes of swordtail fishes, Xiphophorus: variation in the number of loci and existence of ancient gene families

AmpliSAS and AmpliHLA: Web Server and Local Tools for MHC Typing of Non-model Species and Human Using NGS Data

AmpliSAS and AmpliHLA are tools for automatic genotyping of MHC genes from high-throughput sequencing data. AmpliSAS is designed specifically to analyze amplicon sequencing data from non-model species and it is able to perform de novo genotyping without any previous knowledge of the reference alleles. AmpliHLA is a human specific version; it performs HLA typing by comparing sequenced variants against human reference alleles from the IMGT/HLA database. Both tools are available in AmpliSAT web-server as well as scripts for local/server installation. Here we describe the installation and deployment of AmpliSAS and AmpliHLA Perl scripts and dependencies on a local or a server computer. We will show how to run them in the command line using as examples four genotyping protocols: the first two use amplicon sequencing data to genotype the MHC genes of a passerine bird and human respectively; the third and fourth present the HLA typing of a human cell line starting from RNA and exome sequencing data respectively.

The MHC Class Ia Genes in Chenfu's Treefrog (Zhangixalus chenfui) Evolved via Gene Duplication, Recombination, and Selection

Simple Summary

Amphibians, the first terrestrial vertebrates, provide materials for adaptive evolutionary studies, such as the evolution of the major histocompatibility complex (MHC). To date, various MHC evolutionary mechanisms have been identified in frogs, but more research is needed to determine the evolutionary mechanisms of the frog MHC. The main purpose of this study was to evaluate polymorphisms in the MHC class Ia genes of the Chenfu’s Treefrog. The MHC class Ia genes of the Chenfu’s Treefrog have high polymorphism. The mechanisms responsible for the formation of the polymorphisms include gene duplication, recombination, and selection.

Abstract

The molecular mechanisms underlying the evolution of adaptive immunity-related proteins can be deduced by a thorough examination of the major histocompatibility complex (MHC). Currently, in vertebrates, there is a relatively large amount of research on MHCs in mammals and birds. However, research related to amphibian MHC genes and knowledge about the evolutionary patterns is limited. This study aimed to isolate the MHC class I genes from Chenfu’s Treefrog (Zhangixalus chenfui) and reveal the underlying evolutionary processes. A total of 23 alleles spanning the coding region of MHC class Ia genes were identified in 13 individual samples. Multiple approaches were used to test and identify recombination from the 23 alleles. Amphibian MHC class Ia alleles, from NCBI, were used to construct the phylogenetic relationships in MEGA. Additionally, the partition strategy was adopted to construct phylogenetic relationships using MrBayes and MEGA. The sites of positive selection were identified by FEL, PAML, and MEME. In Chenfu’s Treefrog, we found that: (1) recombination usually takes place between whole exons of MHC class Ia genes; (2) there are at least 3 loci for MHC class Ia, and (3) the diversity of genes in MHC class Ia can be attributed to recombination, gene duplication, and positive selection. We characterized the evolutionary mechanisms underlying MHC class Ia genes in Chenfu’s Treefrog, and in so doing, broadened the knowledge of amphibian MHC systems.

AmpliSAS and AmpliHLA: Web Server Tools for MHC Typing of Non-Model Species and Human Using NGS Data

AmpliSAS and AmpliHLA are web server tools for automatic genotyping of MHC genes from high-throughput sequencing data. AmpliSAS is designed specifically to analyze amplicon sequencing data from non-model species and it is able to perform de-novo genotyping without any previous knowledge of the reference alleles. AmpliHLA is a human-specific version, it performs HLA typing by comparing sequenced variants against human reference alleles from the IMGT/HLA database. Here we describe four genotyping protocols: the first two use amplicon sequencing data to genotype the MHC genes of a passerine bird and human respectively; the third and fourth present the HLA typing of a human cell line starting from RNA and exome sequencing data respectively.

Spatial patterns of immunogenetic and neutral variation underscore the conservation value of small, isolated American badger populations

Small and isolated populations often exhibit low genetic diversity due to drift and inbreeding, but may simultaneously harbour adaptive variation. We investigate spatial distributions of immunogenetic variation in American badger subspecies (Taxidea taxus), as a proxy for evaluating their evolutionary potential across the northern extent of the species' range. We compared genetic structure of 20 microsatellites and the major histocompatibility complex (MHC DRB exon 2) to evaluate whether small, isolated populations show low adaptive polymorphism relative to large and well-connected populations. Our results suggest that gene flow plays a prominent role in shaping MHC polymorphism across large spatial scales, while the interplay between gene flow and selection was stronger towards the northern peripheries. The similarity of MHC alleles within subspecies relative to their neutral genetic differentiation suggests that adaptive divergence among subspecies can be maintained despite ongoing gene flow along subspecies boundaries. Neutral genetic diversity was low in small relative to large populations, but MHC diversity within individuals was high in small populations. Despite reduced neutral genetic variation, small and isolated populations harbour functional variation that likely contribute to the species evolutionary potential at the northern range. Our findings suggest that conservation approaches should focus on managing adaptive variation across the species range rather than protecting subspecies per se.

Copy number variation and genetic diversity of MHC Class IIb alleles in an alien population of Xenopus laevis

Xenopus laevis (the African clawed frog), which originated through hybridisation and whole genome duplication, has been used as a model for genetics and development for many years, but surprisingly little is known about immune gene variation in natural populations. The purpose of this study was to use an isolated population of X. laevis that was introduced to Wales, UK in the past 50 years to investigate how variation at the MHC compares to that at other loci, following a severe population bottleneck. Among 18 individuals, we found nine alleles based on exon 2 sequences of the Class IIb region (which includes the peptide binding region). Individuals carried from one to three of the loci identified from previous laboratory studies. Genetic variation was an order of magnitude higher at the MHC compared with three single-copy nuclear genes, but all loci showed high levels of heterozygosity and nucleotide diversity and there was not an excess of homozygosity or decrease in diversity over time that would suggest extensive inbreeding in the introduced population. Tajima’s D was positive for all loci, which is consistent with a bottleneck. Moreover, comparison with published sequences identified the source of the introduced population as the Western Cape region of South Africa, where most commercial suppliers have obtained their stocks. These factors suggest that despite founding by potentially already inbred individuals, the alien population in Wales has maintained substantial genetic variation at both adaptively important and neutral genes.

The MHC class I genes of zebrafish

Major histocompatibility complex (MHC) molecules play a central role in the immune response and in the recognition of non-self. Found in all jawed vertebrate species, including zebrafish and other teleosts, MHC genes are considered the most polymorphic of all genes. In this review we focus on the multi-faceted diversity of zebrafish MHC class I genes, which are classified into three sequence lineages: U, Z, and L. We examine the polygenic, polymorphic, and haplotypic diversity of the zebrafish MHC class I genes, discussing known and postulated functional differences between the different class I lineages. In addition, we provide the first comprehensive nomenclature for the L lineage genes in zebrafish, encompassing at least 15 genes, and characterize their sequence properties. Finally, we discuss how recent findings have shed new light on the remarkably diverse MHC loci of this species.

Evolution by selection, recombination, and gene duplication in MHC class I genes of two Rhacophoridae species

BACKGROUND

Comparison of major histocompatibility complex (MHC) genes across vertebrate species can reveal molecular mechanisms underlying the evolution of adaptive immunity-related proteins. As the first terrestrial tetrapods, amphibians deserve special attention because of their exposure to probably increased spectrum of microorganisms compared with ancestral aquatic fishes. Knowledge regarding the evolutionary patterns and mechanisms associated with amphibian MHC genes remains limited. The goal of the present study was to isolate MHC class I genes from two Rhacophoridae species (Rhacophorus omeimontis and Polypedates megacephalus) and examine their evolution.

RESULTS

We identified 27 MHC class I alleles spanning the region from exon 2 to 4 in 38 tree frogs. The available evidence suggests that these 27 sequences all belong to classical MHC class I (MHC Ia) genes. Although several anuran species only display one MHC class Ia locus, at least two or three loci were observed in P. megacephalus and R. omeimontis, indicating that the number of MHC class Ia loci varies among anuran species. Recombination events, which mainly involve the entire exons, played an important role in shaping the genetic diversity of the 27 MHC class Ia alleles. In addition, signals of positive selection were found in Rhacophoridae MHC class Ia genes. Amino acid sites strongly suggested by program to be under positive selection basically accorded with the putative antigen binding sites deduced from crystal structure of human HLA. Phylogenetic relationships among MHC class I alleles revealed the presence of trans-species polymorphisms.

CONCLUSIONS

In the two Rhacophoridae species (1) there are two or three MHC class Ia loci; (2) recombination mainly occurs between the entire exons of MHC class Ia genes; (3) balancing selection, gene duplication and recombination all contribute to the diversity of MHC class Ia genes. These findings broaden our knowledge on the evolution of amphibian MHC systems.

Molecular cloning and characterization of sea bass (Dicentrarchus labrax, L.) MHC class I heavy chain and β2-microglobulin

In this work, the gene and cDNA of sea bass (Dicentrarchus labrax) β2-microglobulin (Dila-β2m) and several cDNAs of MHC class I heavy chain (Dila-UA) were characterized. While Dila-β2m is single-copy, numerous Dila-UA transcripts were identified per individual with variability at the peptide-binding domain (PBD), but also with unexpected diversity from the connective peptide (CP) through the 3' untranslated region (UTR). Phylogenetic analysis segregates Dila-β2m and Dila-UA into each subfamily cluster, placing them in the fish class and branching Dila-MHC-I with lineage U. The α1 domains resemble those of the recently proposed L1 trans-species lineage. Although no Dila-specific α1, α2 or α3 sub-lineages could be observed, two highly distinct sub-lineages were identified at the CP/TM/CYT regions. The three-dimensional homology model of sea bass MHC-I complex is consistent with other characterized vertebrate structures. Furthermore, basal tissue-specific expression profiles were determined for both molecules, and expression of β2m was evaluated after poly I:C stimulus. Results suggest these molecules are orthologues of other β2m and teleost classical MHC-I and their basic structure is evolutionarily conserved, providing relevant information for further studies on antigen presentation in this fish species.

The research of W.E. Mayer (1953-2012): a spectrum of immune systems

Over a period of some 20 years, Werner Eugen Mayer played a significant role in establishing a framework for molecular studies of Mhc genes in multiple vertebrates. His work largely concerned gene isolation, sequencing, and related bioinformatic analyses both for the Mhc and for immune system genes of about 200 species, ranging from apes, monkeys, rodents, and marsupials, through to birds, bony fishes, and lampreys. In addition to his exploration of diverse Mhc genes, Werner is remembered for playing a critical role in the development of two important insights into the evolution of immune systems. His was among the first published DNA sequence-based descriptions of trans-species evolution of Mhc alleles, including the first description of the long-lived polymorphisms shared by humans and chimpanzees. This research opened the way for using Mhc polymorphisms in demographic analyses. The second important insight in which he played a prominent role involved the characterization of immune cells and their expressed genes in the lamprey, a jawless vertebrate. His findings helped to indicate the considerable degree to which extant immune mechanisms were co-opted in the creation of the adaptive immune system of jawed vertebrates.

Retained orthologous relationships of the MHC Class I genes during euteleost evolution

Major histocompatibility complex (MHC) class I molecules play a pivotal role in immune defense system, presenting the antigen peptides to cytotoxic CD8+ T lymphocytes. Most vertebrates possess multiple MHC class I loci, but the analysis of their evolutionary relationships between distantly related species has difficulties because genetic events such as gene duplication, deletion, recombination, and/or conversion have occurred frequently in these genes. Human MHC class I genes have been conserved only within the primates for up to 46-66 My. Here, we performed comprehensive analysis of the MHC class I genes of the medaka fish, Oryzias latipes, and found that they could be classified into four groups of ancient origin. In phylogenetic analysis using these genes and the classical and nonclassical class I genes of other teleost fishes, three extracellular domains of the class I genes showed quite different evolutionary histories. The α1 domains generated four deeply diverged lineages corresponding to four medaka class I groups with high bootstrap values. These lineages were shared with salmonid and/or other acanthopterygian class I genes, unveiling the orthologous relationships between the classical MHC class I genes of medaka and salmonids, which diverged approximately 260 Ma. This suggested that the lineages must have diverged in the early days of the euteleost evolution and have been maintained for a long time in their genome. In contrast, the α3 domains clustered by species or fish groups, regardless of classical or nonclassical gene types, suggesting that this domain was homogenized in each species during prolonged evolution, possibly retaining the potential for CD8 binding even in the nonclassical genes. On the other hand, the α2 domains formed no apparent clusters with the α1 lineages or with species, suggesting that they were diversified partly by interlocus gene conversion, and that the α1 and α2 domains evolved separately. Such evolutionary mode is characteristic to the teleost MHC class I genes and might have contributed to the long-term conservation of the α1 domain.

Evolutionary analysis of two classical MHC class I loci of the medaka fish, Oryzias latipes: haplotype-specific genomic diversity, locus-specific polymorphisms, and interlocus homogenization

The major histocompatibility complex (MHC) region of the teleost medaka (Oryzias latipes) contains two classical class I loci, UAA and UBA, whereas most lower vertebrates possess or express a single locus. To elucidate the allelic diversification and evolutionary relationships of these loci, we compared the BAC-based complete genomic sequences of the MHC class I region of three medaka strains and the PCR-based cDNA sequences of two more strains and two wild individuals, representing nine haplotypes. These were derived from two geographically distinct medaka populations isolated for four to five million years. Comparison of the genomic sequences showed a marked diversity in the region encompassing UAA and UBA even between the strains derived from the same population, and also showed an ancient divergence of these loci. cDNA analysis indicated that the peptide-binding domains of both UAA and UBA are highly polymorphic and that most of the polymorphisms were established in a locus-specific manner before the divergence of the two populations. Interallelic recombination between exons 2 and 3 encoding these domains was observed. The second intron of the UAA genes contains a highly conserved region with a palindromic sequence, suggesting that this region contributed to the recombination events. In contrast, the alpha3 domain is extremely homogenized not only within each locus but also between UAA and UBA regardless of populations. Two lineages of the transmembrane and cytoplasmic regions are also shared by UAA and UBA, suggesting that these two loci evolved with intimate genetic interaction through gene conversion or unequal crossing over.

Resistance of common carp (Cyprinus carpio L.) to Cyprinid herpesvirus-3 is influenced by major histocompatibility (MH) class II B gene polymorphism

The role of MH class II B (Cyca-DAB1-like) genes in resistance of common carp (Cyprinus carpio L.) to Cyprinid herpesvirus-3 (CyHV-3), also known as koi herpesvirus (KHV) was analysed. The material consisted of 934 fish from six carp crosses. Fish were challenged with CyHV-3 at an age of 7 and 10 months. During challenge experiments the peak of mortality caused by CyHV-3 was observed at days 8-12 p.i. and the overall cumulative mortality reached 79.9%. Among six Cyca-DAB1-like genotypes, revealed by PCR-RF-SSCP analysis, one genotype (E) was found associated with higher resistance to CyHV-3. Three other genotypes (B, H and J) could be linked to higher susceptibility to CyHV-3. Analysis of the alleles that compose the Cyca-DAB1-like genotypes linked one particular allele (Cyca-DAB1*05) to significantly increased, and two alleles (Cyca-DAB1*02 and Cyca-DAB1*06) to significantly decreased resistance to CyHV-3. Our data indicate that MH class II B genes could be used as potential genetic markers in breeding of common carp for resistance to this virus.

Characterisation of grass carp (Ctenopharyngodon idellus) MHC class I domain lineages

In order to characterise grass carp MHC class I (Ctid-MHC I) sequences, 26 Ctid-MHC I genes were cloned from 12 individuals and their alpha domain lineages were analysed. Simultaneously, a quantitative reverse transcription-polymerase chain reaction (Q-RT-PCR) assay was developed to detect Ctid-MHC I tissue-specific expression. The results suggested that Ctid-MHC I could be divided into eight lineages (Ctid-NA-Ctid-NH). Based on whether they contained the motif of eight key amino acids (YYRTKWYY), Ctid-MHC I lineages were divided into two groups [Ctid-MHC I (8(+)) and Ctid-MHC I (8(-))]. The expression analysis showed that the Ctid-MHC I locus/loci appeared in the kidney, gill, intestine, heart, spleen, liver, and brain. A GenBank homology BLAST was performed independently with each alpha domain, and Ctid-MHC I alpha1, alpha2, and alpha3 were categorised into two (V and IX), five (II, IV-VII), and four (IV-VII) domain lineages, respectively. Based on the alphabetic labelling system created in our earlier studies, one alpha1 (IX), four alpha2 (IV-VII), and unique alpha3 (V-VII) domain lineages were observed in grass carp and across the teleostean species.

Comparative genomics of major histocompatibility complexes

The major histocompatibility complex (MHC) is a gene dense region found in all jawed vertebrates examined to date. The MHC contains a high percentage of immune genes, in particular genes involved in antigen presentation, which are generally highly polymorphic. The region plays an important role in disease resistance. The clustering of MHC genes could be advantageous for co-evolution or regulation, and its study in many species is desirable. Even though some linkage of MHC genes is apparent in all gnathostomes, the genomic organization can differ greatly by species, suggesting rapid evolution of MHC genes after divergence from a common ancestor. Previous reviews of comparative MHC organization have been written when relatively fragmentary sequence and mapping data were available on many species. This review compares maps of MHC gene orders in commonly studied species, where extensive sequencing has been performed.

The cattle major histocompatibility complex: is it unique?

Major histocompatibility complex (MHC) class I genes encode highly polymorphic molecules that are expressed on virtually every cell type, and have been identified in all but the most primitive vertebrates. They play a number of crucial roles in the immune response to infectious disease. Most information regarding MHC genes has been generated from humans and mice but, because of the great variability found in the MHC system, it is not always possible to extrapolate from these to other species. Many strategies have evolved to maximise the ability of the MHC to protect individuals and populations against pathogens. Cattle MHC class I genes exhibit a number of unusual features. Evidence from mapping studies, haplotype and phylogenetic analyses suggests the presence of six classical class I loci, in contrast to the more usual two or three, and these are expressed in various combinations of one, two or three on different haplotypes. Although it remains difficult to assign alleles to loci, it appears that none of the loci are expressed on all haplotypes. There is currently limited information relating to polymorphism, but various approaches suggest diversity is high, and may vary between breeds/populations. Functional consequences of variable MHC haplotype composition are discussed. Identifying unique features of the MHC in cattle will lead to new insights into evolution of the immune system.

Major histocompatibility complex diversity influences parasite resistance and innate immunity in sticklebacks

Proteins of the major histocompatibility complex (MHC) play a central role in the presentation of antigens to the adaptive immune system. The MHC also influences the odour-based choice of mates in humans and several animal taxa. It has recently been shown that female three-spined sticklebacks (Gasterosteus aculeatus) aim at a moderately high MHC diversity in their offspring when choosing a mate. Do they optimize the immune systems of their offspring? Using three-spined sticklebacks that varied in their individual numbers of MHC class IIB molecules, we tested, experimentally, whether allelic diversity at the MHC influences parasite resistance and immune parameters. We found that sticklebacks with low MHC diversity suffered more from parasite infection after experimental exposure to Schistocephalus solidus tapeworms and Glugea anomala microsporidians. They also showed the highest proportion of granulocytes and the strongest respiratory burst reaction, which are correlates of innate immunity. This indicates a strong activity of the innate immune system after challenge by parasites when MHC diversity is suboptimal. Individuals with very high allelic diversity at the MHC seemed inferior to those with moderately high diversity. Such a pattern is consistent with theoretical expectations of an optimal balance between the number of recognizable antigens and self-tolerance.

Clustering of C-Type Lectin Natural Killer Receptor-Like Loci in the Bony Fish Oreochromis niloticus

The genome of the cichlid (teleost) fish Oreochromis niloticus contains a set of genes which encode group V C-type lectin proteins homologous to the mammalian NKG2/CD94 family of natural killer (NK) cell receptors. To determine the genomic organization of these killer cell-like receptor (KLR) genes, an O. niloticus BAC library was screened with a cDNA probe derived previously from an expressed sequence tag of the related cichlid species Paralabidochromis chilotes. Four distinct KLR-bearing BAC clones were analysed, three of which could be assembled into a contig. One of the clones was sequenced in its entirety, whereas the others were partially sequenced to identify the KLR loci borne by them. Altogether, 28 distinct KLR loci were identified, of which at least 26 occupy a single chromosomal region, the KLR complex. One half of the loci appear to be occupied by pseudogenes. Compared to the human NK cell receptor complex, the Oreochromis KLR complex is more compact and, apart from transposons, appears to contain only KLR loci. The gene density of the complex is one KLR locus per 18 kb of sequence. All the KLR loci constituting the complex are derived from a single most recent common ancestor, which is estimated to have existed 7.7 million years ago. The 180 kb of the determined sequence is a mosaic of blocks of similar segments reflecting a complex history of duplications, deletions and rearrangements. The transposons found in the sequenced part belong to the TC1, Xena, CR1 and TX1 families.

A novel functional class I lineage in zebrafish (Danio rerio), carp (Cyprinus carpio), and large barbus (Barbus intermedius) showing an unusual conservation of the peptide binding domains

Species from all major jawed vertebrate taxa possess linked polymorphic class I and II genes located in an MHC. The bony fish are exceptional with class I and II genes located on different linkage groups. Zebrafish (Danio rerio), common carp (Cyprinus carpio), and barbus (Barbus intermedius) represent highly divergent cyprinid genera. The genera Danio and Cyprinus diverged 50 million years ago, while Cyprinus and Barbus separated 30 million years ago. In this study, we report the first complete protein-coding class I ZE lineage cDNA sequences with high similarity between the three cyprinid species. Two unique complete protein-coding cDNA sequences were isolated in zebrafish, Dare-ZE*0101 and Dare-ZE*0102, one in common carp, Cyca-ZE*0101, and six in barbus, Bain-ZE*0101, Bain-ZE*0102, Bain-ZE*0201, Bain-ZE*0301, Bain-ZE*0401, and Bain-ZE*0402. Deduced amino acid sequences indicate that these sequences encode bonafide class I proteins. In addition, the presence of conserved potential peptide anchoring residues, exon-intron organization, ubiquitous expression, and polymorphism generated by positive selection on putative peptide binding residues support a classical nature of class I ZE lineage genes. Phylogenetic analyses revealed clustering of the ZE lineage clade with nonclassical cyprinid class I Z lineage clade away from classical cyprinid class I genes, suggesting a common ancestor of these nonclassical genes as observed for mammalian class I genes. Data strongly support the classical nature of these ZE lineage genes that evolved in a trans-species fashion with lineages being maintained for up to 100 million years as estimated by divergence time calculations.