Transspecies dimorphic allelic lineages of the proteasome subunit beta-type 8 gene (PSMB8) in the teleost genus Oryzias

Coevolution between MHC class I and Antigen Processing Genes in salamanders

Proteins encoded by antigen-processing genes (APGs) provide major histocompatibility complex (MHC) class I (MHC-I) with antigenic peptides. In mammals, polymorphic multigenic MHC-I family is served by monomorphic APGs, whereas in certain nonmammalian species both MHC-I and APGs are polymorphic and coevolve within stable haplotypes. Coevolution was suggested as an ancestral gnathostome feature, presumably enabling only a single highly expressed classical MHC-I gene. In this view coevolution, while optimizing some aspects of adaptive immunity, would also limit its flexibility by preventing the expansion of classical MHC-I into a multigene family. However, some nonmammalian taxa, such as salamanders, have multiple highly expressed MHC-I genes, suggesting either that coevolution is relaxed or that it does not prevent the establishment of multigene MHC-I. To distinguish between these two alternatives, we use salamanders (30 species from 16 genera representing six families) to test, within a comparative framework, a major prediction of the coevolution hypothesis: the positive correlation between MHC-I and APG diversity. We found that MHC-I diversity explained both within-individual and species-wide diversity of two APGs, TAP1 and TAP2, supporting their coevolution with MHC-I, whereas no consistent effect was detected for the other three APGs (PSMB8, PSMB9, and TAPBP). Our results imply that although coevolution occurs in salamanders, it does not preclude the expansion of the MHC-I gene family. Contrary to the previous suggestions, nonmammalian vertebrates thus may be able to accommodate diverse selection pressures with flexibility granted by rapid expansion or contraction of the MHC-I family, while retaining the benefits of coevolution between MHC-I and TAPs.

Molecular Evolution of Antigen-Processing Genes in Salamanders: Do They Coevolve with MHC Class I Genes?

Proteins encoded by antigen-processing genes (APGs) prepare antigens for presentation by the major histocompatibility complex class I (MHC I) molecules. Coevolution between APGs and MHC I genes has been proposed as the ancestral gnathostome condition. The hypothesis predicts a single highly expressed MHC I gene and tight linkage between APGs and MHC I. In addition, APGs should evolve under positive selection, a consequence of the adaptive evolution in MHC I. The presence of multiple highly expressed MHC I genes in some teleosts, birds, and urodeles appears incompatible with the coevolution hypothesis. Here, we use urodele amphibians to test two key expectations derived from the coevolution hypothesis: 1) the linkage between APGs and MHC I was studied in Lissotriton newts and 2) the evidence for adaptive evolution in APGs was assessed using 42 urodele species comprising 21 genera from seven families. We demonstrated that five APGs (PSMB8, PSMB9, TAP1, TAP2, and TAPBP) are tightly linked (<0.5 cM) to MHC I. Although all APGs showed some codons under episodic positive selection, we did not find a pervasive signal of positive selection expected under the coevolution hypothesis. Gene duplications, putative gene losses, and divergent allelic lineages detected in some APGs demonstrate considerable evolutionary dynamics of APGs in salamanders. Overall, our results indicate that if coevolution between APGs and MHC I occurred in urodeles, it would be more complex than envisaged in the original formulation of the hypothesis.

Whole genome duplications have provided teleosts with many roads to peptide loaded MHC class I molecules

Background

In sharks, chickens, rats, frogs, medaka and zebrafish there is haplotypic variation in MHC class I and closely linked genes involved in antigen processing, peptide translocation and peptide loading. At least in chicken, such MHCIa haplotypes of MHCIa, TAP2 and Tapasin are shown to influence the repertoire of pathogen epitopes being presented to CD8+ T-cells with subsequent effect on cell-mediated immune responses.

Results

Examining MHCI haplotype variation in Atlantic salmon using transcriptome and genome resources we found little evidence for polymorphism in antigen processing genes closely linked to the classical MHCIa genes. Looking at other genes involved in MHCI assembly and antigen processing we found retention of functional gene duplicates originating from the second vertebrate genome duplication event providing cyprinids, salmonids, and neoteleosts with the potential of several different peptide-loading complexes. One of these gene duplications has also been retained in the tetrapod lineage with orthologs in frogs, birds and opossum.

Conclusion

We postulate that the unique salmonid whole genome duplication (SGD) is responsible for eliminating haplotypic content in the paralog MHCIa regions possibly due to frequent recombination and reorganization events at early stages after the SGD. In return, multiple rounds of whole genome duplications has provided Atlantic salmon, other teleosts and even lower vertebrates with alternative peptide loading complexes. How this affects antigen presentation remains to be established.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1138-9) contains supplementary material, which is available to authorized users.

Trans-Species Polymorphism in Immune Genes: General Pattern or MHC-Restricted Phenomenon?

Immunity exhibits extraordinarily high levels of variation. Evolution of the immune system in response to host-pathogen interactions in particular ecological contexts appears to be frequently associated with diversifying selection increasing the genetic variability. Many studies have documented that immunologically relevant polymorphism observed today may be tens of millions years old and may predate the emergence of present species. This pattern can be explained by the concept of trans-species polymorphism (TSP) predicting the maintenance and sharing of favourable functionally important alleles of immune-related genes between species due to ongoing balancing selection. Despite the generality of this concept explaining the long-lasting adaptive variation inherited from ancestors, current research in TSP has vastly focused only on major histocompatibility complex (MHC). In this review we summarise the evidence available on TSP in human and animal immune genes to reveal that TSP is not a MHC-specific evolutionary pattern. Further research should clearly pay more attention to the investigation of TSP in innate immune genes and especially pattern recognition receptors which are promising candidates for this type of evolution. More effort should also be made to distinguish TSP from convergent evolution and adaptive introgression. Identification of balanced TSP variants may represent an accurate approach in evolutionary medicine to recognise disease-resistance alleles.

Evolution of dimorphisms of the proteasome subunit beta type 8 gene (PSMB8) in basal ray-finned fish

The proteasome subunit beta type 8 (PSMB8) gene encodes a catalytic subunit of immunoproteasome that plays a central role in the processing of antigenic peptides presented by major histocompatibility complex class I molecules. The A- and F-type alleles defined by the 31st amino acid residue determining cleaving specificity have been identified from ray-finned fish, amphibia, and reptiles. These two types show extremely long-term trans-species polymorphism in Polypteriformes, Cypriniformes, and Salmoniformes, suggesting the presence of very ancient lineages termed A and F. To elucidate the evolution of the PSMB8 dimorphism in basal ray-finned fish, we analyzed Pantodon buchholzi (Osteoglossiformes), seven species of Anguilliformes, and Hypomesus nipponensis (Osmeriformes). Both A and F lineage sequences were identified from P. buchholzi and H. nipponensis, confirming that these two lineages have been conserved by basal ray-finned fish. However, both the A- and F-type alleles found in Anguilliformes species belonged to the F lineage irrespective of their types. This apparently suggests that the A lineage was lost in the common ancestor of Anguilliformes, and recovery of the A type within the F lineage occurred in Anguilliformes. The apparent loss of the F lineage and recovery of the F type within the A lineage have already been reported from tetrapods and higher teleosts. However, this is the first report on the reverse situation and reveals the dynamic evolution of the PSMB8 dimorphism.

Dimorphisms of the proteasome subunit beta type 8 gene (PSMB8) of ectothermic tetrapods originated in multiple independent evolutionary events

The proteasome subunit beta type 8 gene (PSMB8) encodes one of the beta subunits of the immunoproteasome responsible for the generation of peptides presented by major histocompatibility complex class I molecules. Dimorphic alleles of the PSMB8 gene, termed A and F types, based on the deduced 31st amino acid residue of the mature protein have been reported from various vertebrates. Phylogenetic analysis revealed the presence of dichotomous ancient lineages, one comprising the F-type PSMB8 of basal ray-finned fishes, and the other comprising the A-type PSMB8 of these animals and both the F- and A-type PSMB8 of Xenopus and acanthopterygians, indicating that evolutionary history of the PSMB8 dimorphism was not straightforward. We analyzed the PSMB8 gene of five reptile and one amphibian species and found both the A and F types from all six. Phylogenetic analysis indicated that the PSMB8 F type was apparently regenerated from the PSMB8 A type at least five times independently during tetrapod evolution. Genomic typing of wild individuals of geckos and newts indicated that the frequencies of the A- and F-type alleles are not highly biased in these species. Phylogenetic analysis of each exon of the reptile PSMB8 gene suggested interallelic sequence homogenization as a possible evolutionary mechanism for the apparent recurrent regeneration of PSMB8 dimorphism in tetrapods. An extremely strong balancing selection acting on PSMB8 dimorphism was implicated in an unprecedented pattern of allele evolution.

Where are the wormy mice? A reexamination of hybrid parasitism in the European house mouse hybrid zone

Wormy mice in a hybrid zone have been interpreted as evidence of low hybrid fitness, such that parasites contribute to species separation. However, because of its natural heterogeneity, observations of parasite load must be numerous with good field area coverage. We sampled 689 mice from 107 localities across the Bavaria-Bohemia region of the European house mouse hybrid zone and calculated their hybrid indices using 1401 diagnostic single nucleotide polymorphisms (SNPs). We tested whether hybrids have greater or lesser diversity and load of parasite helminths than additive expectations, performing load analyses on the four most common taxa. We found hybrids have significantly reduced diversity and load of each of the commonest helminths; rarer helminths further support reduced load. Although within-locality comparisons have little power, randomization tests show the repeated pattern is unlikely to be due to local parasite heterogeneity, and simulations show a patch of low parasite diversity is unlikely to fall by chance just so in the field area, such that it produces the observed effects. Our data therefore contradict the idea that helminths reduce hybrid fitness through increased load. We discuss a vicariant Red Queen model that implies immune genes tracking parasites will escape Dobzhansky-Muller incompatibilities, generating hybrid variants untargeted by parasites.

Long-lived dichotomous lineages of the proteasome subunit beta type 8 (PSMB8) gene surviving more than 500 million years as alleles or paralogs

On an evolutionary time scale, polymorphic alleles are believed to have a short life, persisting at most tens of millions of years even under long-term balancing selection. Here, we report highly diverged trans-species dimorphism of the proteasome subunit beta type 8 (PSMB8) gene, which encodes a catalytic subunit of the immunoproteasome responsible for the generation of peptides presented by major histocompatibility complex (MHC) class I molecules, in lower teleosts including Cypriniformes (zebrafish and loach) and Salmoniformes (trout and salmon), whose last common ancestor dates to 300 Ma. Moreover, phylogenetic analyses indicated that these dimorphic alleles share lineages with two shark paralogous genes, suggesting that these two lineages have been maintained for more than 500 My either as alleles or as paralogs, and that conversion between alleles and paralogs has occurred at least once during vertebrate evolution. Two lineages termed PSMB8A and PSMB8F show an A(31)F substitution that would probably affect their cleaving specificity, and whereas the PSMB8A lineage has been retained by all analyzed jawed vertebrates, the PSMB8F lineage has been lost by most jawed vertebrates except for cartilaginous fish and basal teleosts. However, a possible functional equivalent of the PSMB8F lineage has been revived as alleles within the PSMB8A lineage at least twice during vertebrate evolution in the amphibian Xenopus and teleostean Oryzias species. Dynamic evolution of the PSMB8 polymorphism through long-term persistence, loss, and regaining of dimorphism and conversion between alleles and paralogs implies the presence of strong selective pressure for functional polymorphism of this gene.

Highly divergent dimorphic alleles of the proteasome subunit beta type-8 (PSMB8) gene of the bichir Polypterus senegalus: implication for evolution of the PSMB8 gene of jawed vertebrates

The proteasome subunit beta type-8 (PSMB8) gene encodes a catalytic subunit of the immunoproteasome, which is involved in the generation of peptides presented by MHC class I molecules. To date, highly diverged dichotomous alleles of PSMB8 have been reported in Oryzias species (actinopterygian teleosts) and Xenopus species (sarcopterygian amphibians). These dimorphic alleles share a similar substitution (A/V(31)F/Y) at the 31st position of the mature protein, which is most probably involved in formation of the S1 pocket. This substitution likely confers different cleavage specificities on the dimorphic PSMB8s. In addition, two paralogous PSMB8 genes possessing the A and F residues at the 31st position have been reported in sharks. Phylogenetic analysis indicated that the two types of PSMB8 of Oryzias, Xenopus, and sharks arose by independent evolutionary events. Here, we identified another pair of dimorphic alleles of PSMB8, which have the A and F residues at the 31st position of the mature protein, from bichir, Polypterus senegalus, a basal actinopterygian. The sequences of the mature proteins-encoding region of the dimorphic alleles of bichir PSMB8, the A and F types, showed only 72.7% and 77.5% identities at the nucleotide and the deduced amino acid levels, respectively. Their intronic sequences show almost no similarity, indicating that the dimorphic alleles of bichir PSMB8 have a very ancient origin. However, phylogenetic analysis showed that the dimorphisms of PSMB8 of bichir, Xenopus, and Oryzias arose by independent evolutionary events, suggesting the presence of a strong selective pressure for possessing the dimorphism.