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

Parallel Selection in Domesticated Atlantic Salmon from Divergent Founders Including on Whole-Genome Duplication-derived Homeologous Regions

Domestication and artificial selection for desirable traits have driven significant phenotypic changes and left detectable genomic footprints in farmed animals. Since the 1960s, intensive breeding has led to the rapid domestication of Atlantic salmon (Salmo salar), with multiple independent events that make it a valuable model for studying early domestication stages and the parallel evolution of populations of different origins subjected to similar selection pressures. Some aquatic species, including Atlantic salmon, have undergone whole-genome duplication (WGD), raising the possibility that genetic redundancy resulting from WGD has contributed to adaptation in captive environments, as seen in plants. Here, we examined the genomic responses to domestication in Atlantic salmon, focusing on potential signatures of parallel selection, including those associated with WGD. Candidate genomic regions under selection were identified by comparing whole-genome sequences from aquaculture and wild populations across 2 independently domesticated lineages (Western Norway and North America) using a genome-wide scan that combined 3 statistical methods: allele frequencies (FST), site frequency (Tajima's D), and haplotype differentiation (XP-EHH). These analyses revealed shared selective sweeps on identical SNPs in major histocompatibility complex (MHC) genes across aquaculture populations. This suggests that a combination of long-term balancing selection and recent human-induced selection has shaped MHC gene evolution in domesticated salmon. Additionally, we observed selective sweeps on a small number of gene pairs in homeologous regions originating from WGD, offering insights into how historical genome duplication events may intersect with recent selection pressures in aquaculture species.

The Evolution of Extreme Genetic Variability in a Parasite-Resistance Complex

Genomic regions that play a role in parasite defense are often found to be highly variable, with the major histocompatibility complex serving as an iconic example. Single nucleotide polymorphisms may represent only a small portion of this variability, with Indel polymorphisms and copy number variation further contributing. In extreme cases, haplotypes may no longer be recognized as orthologous. Understanding the evolution of such highly divergent regions is challenging because the most extreme variation is not visible using reference-assisted genomic approaches. Here we analyze the case of the Pasteuria Resistance Complex in the crustacean Daphnia magna, a defense complex in the host against the common and virulent bacterium Pasteuria ramosa. Two haplotypes of this region have been previously described, with parts of it being nonhomologous, and the region has been shown to be under balancing selection. Using pan-genome analysis and tree reconciliation methods to explore the evolution of the Pasteuria Resistance Complex and its characteristics within and between species of Daphnia and other Cladoceran species, our analysis revealed a remarkable diversity in this region even among host species, with many nonhomologous hyper-divergent haplotypes. The Pasteuria Resistance Complex is characterized by extensive duplication and losses of Fucosyltransferase (FuT) and Galactosyltransferase (GalT) genes that are believed to play a role in parasite defense. The Pasteuria Resistance Complex region can be traced back to common ancestors over 250 million years. The unique combination of an ancient resistance complex and a dynamic, hyper-divergent genomic environment presents a fascinating opportunity to investigate the role of such regions in the evolution and long-term maintenance of resistance polymorphisms. Our findings offer valuable insights into the evolutionary forces shaping disease resistance and adaptation, not only in the genus Daphnia, but potentially across the entire Cladocera class.

Convergent copy number increase of genes associated with freshwater colonization in fishes

Copy number variation (CNV) can cause phenotypic changes. However, in contrast to amino acid substitutions and cis-regulatory changes, little is known about the functional categories of genes in which CNV is important for adaptation to novel environments. It is also unclear whether the same genes repeatedly change the copy numbers for adapting to similar environments. Here, we investigate CNV associated with freshwater colonization in fishes, which was observed multiple times across different lineages. Using 48 ray-finned fishes across diverse orders, we identified 23 genes whose copy number increases were associated with freshwater colonization. These genes showed enrichment for peptide receptor activity, hexosyltransferase activity and unsaturated fatty acid metabolism. We further revealed that three of the genes showed copy number increases in freshwater populations compared to marine ancestral populations of the stickleback genus Gasterosteus. These results indicate that copy number increases of genes involved in fatty acid metabolism (FADS2), immune function (PSMB8a) and thyroid hormone metabolism (UGT2) may be important for freshwater colonization at both the inter-order macroevolutionary scale and at the intra-genus microevolutionary scale. Further analysis across diverse taxa will help to understand the role of CNV in the adaptation to novel environments. This article is part of the theme issue 'Genetic basis of adaptation and speciation: from loci to causative mutations'.

Gut-derived Shewanella induces the differentially expressed proteins in leukocytes of Lampetra japonica

Lampreys, one of the most basal jawless vertebrates, are an excellent animal model for investigating vertebrate evolution, embryonic development, and the origin of adaptive immunity. Gut-derived Shewanella strain was isolated and then proved to induce adaptive immunity response in lampreys. Using Shewanella as the antigen, the effect of gut-derived Shewanella on lamprey leukocyte proteome was investigated via label-free liquid chromatography-tandem mass spectrometry for quantitative proteomics analysis. Twenty-five differentially expressed proteins in lamprey leukocytes were identified with significant differences. The differentially expressed proteins were associated with several biological processes. Among these proteins, the signal transducer and activator of transcription 3 (STAT3) was significantly upregulated in leukocytes after Shewanella immunization, indicating that lamprey STAT3 (L-STAT3) was involved in Shewanella-lamprey interactions. Expression pattern analysis revealed that L-STAT3 was mainly distributed in the cytoplasm and upregulated in other tissues after Shewanella immunization. Moreover, L-STAT3 overexpression could promote HEK-293 T and HeLa cell proliferation. However, the functions of L-STAT3 in the adaptive immune response of lamprey induced by gut-derived Shewanella remain to be explored. Overall, the identification of leukocyte proteins involved in Shewanella-lamprey interactions provides important information for understanding the variable lymphocyte receptor-based adaptive immune signal pathways in lampreys. SIGNIFICANCE: Lampreys are considered to be an excellent animal model for studying the origin and development of adaptive immune systems in vertebrates. Lampreys use variable lymphocyte receptors (VLRs) in recognizing antigens. However, the understanding of the VLR-based adaptive immune signal pathways in lampreys remains unclear. Intestinal bacteria could regulate the development of host immune systems. The attempts of inducing lamprey leukocyte differentially expressed proteins using the gut bacterial as the antigen will supply an promising avenue to explore the molecular mechanism of the intestinal bacteria interaction with it's host. Also, the identification of differentially expressed proteins involved in interactions between gut-derived Shewanella and lamprey will supply clues for understanding the variable lymphocyte receptor-based adaptive immune signal pathways in lampreys.

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.

Major Histocompatibility Complex (MHC) Genes and Disease Resistance in Fish

Fascinating about classical major histocompatibility complex (MHC) molecules is their polymorphism. The present study is a review and discussion of the fish MHC situation. The basic pattern of MHC variation in fish is similar to mammals, with MHC class I versus class II, and polymorphic classical versus nonpolymorphic nonclassical. However, in many or all teleost fishes, important differences with mammalian or human MHC were observed: (1) The allelic/haplotype diversification levels of classical MHC class I tend to be much higher than in mammals and involve structural positions within but also outside the peptide binding groove; (2) Teleost fish classical MHC class I and class II loci are not linked. The present article summarizes previous studies that performed quantitative trait loci (QTL) analysis for mapping differences in teleost fish disease resistance, and discusses them from MHC point of view. Overall, those QTL studies suggest the possible importance of genomic regions including classical MHC class II and nonclassical MHC class I genes, whereas similar observations were not made for the genomic regions with the highly diversified classical MHC class I alleles. It must be concluded that despite decades of knowing MHC polymorphism in jawed vertebrate species including fish, firm conclusions (as opposed to appealing hypotheses) on the reasons for MHC polymorphism cannot be made, and that the types of polymorphism observed in fish may not be explained by disease-resistance models alone.

Origin and evolution of the specialized forms of proteasomes involved in antigen presentation

Proteasomes are a multi-subunit protease complex that produces peptides bound by major histocompatibility complex (MHC) class I molecules. Phylogenetic studies indicate that two specialized forms of proteasomes, immunoproteasomes and thymoproteasomes, and the proteasome activator PA28αβ emerged in a common ancestor of jawed vertebrates which acquired adaptive immunity based on the MHC, T cell receptors, and B cell receptors ~ 500 million years ago. Comparative genomics studies now provide strong evidence that the genes coding for the immunoproteasome subunits emerged by genome-wide duplication. On the other hand, the gene encoding the thymoproteasome subunit β5t emerged by tandem duplication from the gene coding for the β5 subunit. Strikingly, birds lack immunoproteasomes, thymoproteasomes, and the proteasome activator PA28αβ, raising an interesting question of whether they have evolved any compensatory mechanisms.

Mechanisms of Fish Macrophage Antimicrobial Immunity

Overcrowding conditions and temperatures shifts regularly manifest in large-scale infections of farmed fish, resulting in economic losses for the global aquaculture industries. Increased understanding of the functional mechanisms of fish antimicrobial host defenses is an important step forward in prevention of pathogen-induced morbidity and mortality in aquaculture setting. Like other vertebrates, macrophage-lineage cells are integral to fish immune responses and for this reason, much of the recent fish immunology research has focused on fish macrophage biology. These studies have revealed notable similarities as well as striking differences in the molecular strategies by which fish and higher vertebrates control their respective macrophage polarization and functionality. In this review, we address the current understanding of the biological mechanisms of teleost macrophage functional heterogeneity and immunity, focusing on the key cytokine regulators that control fish macrophage development and their antimicrobial armamentarium.

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.

Alternative haplotypes of antigen processing genes in zebrafish diverged early in vertebrate evolution

Antigen processing and presentation genes found within the MHC are among the most highly polymorphic genes of vertebrate genomes, providing populations with diverse immune responses to a wide array of pathogens. Here, we describe transcriptome, exome, and whole-genome sequencing of clonal zebrafish, uncovering the most extensive diversity within the antigen processing and presentation genes of any species yet examined. Our CG2 clonal zebrafish assembly provides genomic context within a remarkably divergent haplotype of the core MHC region on chromosome 19 for six expressed genes not found in the zebrafish reference genome: mhc1uga, proteasome-β 9b (psmb9b), psmb8f, and previously unknown genes psmb13b, tap2d, and tap2e We identify ancient lineages for Psmb13 within a proteasome branch previously thought to be monomorphic and provide evidence of substantial lineage diversity within each of three major trifurcations of catalytic-type proteasome subunits in vertebrates: Psmb5/Psmb8/Psmb11, Psmb6/Psmb9/Psmb12, and Psmb7/Psmb10/Psmb13. Strikingly, nearby tap2 and MHC class I genes also retain ancient sequence lineages, indicating that alternative lineages may have been preserved throughout the entire MHC pathway since early diversification of the adaptive immune system ∼500 Mya. Furthermore, polymorphisms within the three MHC pathway steps (antigen cleavage, transport, and presentation) are each predicted to alter peptide specificity. Lastly, comparative analysis shows that antigen processing gene diversity is far more extensive than previously realized (with ancient coelacanth psmb8 lineages, shark psmb13, and tap2t and psmb10 outside the teleost MHC), implying distinct immune functions and conserved roles in shaping MHC pathway evolution throughout vertebrates.

Uncommon functional properties of the first piscine 26S proteasome from the Antarctic notothenioid Trematomus bernacchii

The antioxidant defense mechanisms have a great impact on the life of Antarctic organisms. The present study could represent the first evidence of a direct involvement of the 26S proteasome in the antioxidant defense systems of fish adapted to cold.

Protein homoeostasis is a fundamental process allowing the preservation of functional proteins and it has a great impact on the life of the Antarctic organisms. However, the effect of low temperatures on protein turnover is poorly understood and the cold-adaptation of the degradation machinery remains an unresolved issue. As the 26S proteasome represents the main proteolytic system devoted to the controlled degradation of intracellular proteins, the purpose of the present study was to investigate the functions of this complex in the notothenioid Trematomus bernacchii, in order to better understand its role in the physiology of Antarctic fish. To this aim, we purified and characterized the 26S proteasome from T. bernacchii and isolated the cDNAs codifying seven of the 14 subunits belonging to the proteasome 20S core particle. Results provided evidences of the high resistance of the piscine 26S proteasome to oxidative agents and of its ‘uncommon’ ability to efficiently hydrolyse oxidized bovine serum albumin (BSA), suggesting that this enzymatic complex could play a key role in the antioxidant defense systems in fish inhabiting permanently cold marine environments. These unique properties were also reflected by the 3D model analysis, which revealed a higher structural stability of the piscine complex respect to the murine template. Finally, a comparative analysis, performed in a variety of tissues collected from T. bernacchii and the temperate fish Dicentrarchus labrax, showed a lower protein retention in the cold-adapted fish, possibly due to a better efficiency of its degradation machinery.

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.

Multiple divergent haplotypes express completely distinct sets of class I MHC genes in zebrafish

The zebrafish is an important animal model for stem cell biology, cancer, and immunology research. Histocompatibility represents a key intersection of these disciplines; however, histocompatibility in zebrafish remains poorly understood. We examined a set of diverse zebrafish class I major histocompatibility complex (MHC) genes that segregate with specific haplotypes at chromosome 19, and for which donor-recipient matching has been shown to improve engraftment after hematopoietic transplantation. Using flanking gene polymorphisms, we identified six distinct chromosome 19 haplotypes. We describe several novel class I U lineage genes and characterize their sequence properties, expression, and haplotype distribution. Altogether, ten full-length zebrafish class I genes were analyzed, mhc1uba through mhc1uka. Expression data and sequence properties indicate that most are candidate classical genes. Several substitutions in putative peptide anchor residues, often shared with deduced MHC molecules from additional teleost species, suggest flexibility in antigen binding. All ten zebrafish class I genes were uniquely assigned among the six haplotypes, with dominant or codominant expression of one to three genes per haplotype. Interestingly, while the divergent MHC haplotypes display variable gene copy number and content, the different genes appear to have ancient origin, with extremely high levels of sequence diversity. Furthermore, haplotype variability extends beyond the MHC genes to include divergent forms of psmb8. The many disparate haplotypes at this locus therefore represent a remarkable form of genomic region configuration polymorphism. Defining the functional MHC genes within these divergent class I haplotypes in zebrafish will provide an important foundation for future studies in immunology and transplantation.

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.

Effects of genotype and dietary fish oil replacement with vegetable oil on the intestinal transcriptome and proteome of Atlantic salmon (Salmo salar)

Background

Expansion of aquaculture requires alternative feeds and breeding strategies to reduce dependency on fish oil (FO) and better utilization of dietary vegetable oil (VO). Despite the central role of intestine in maintaining body homeostasis and health, its molecular response to replacement of dietary FO by VO has been little investigated. This study employed transcriptomic and proteomic analyses to study effects of dietary VO in two family groups of Atlantic salmon selected for flesh lipid content, 'Lean' or 'Fat'.

Results

Metabolism, particularly of lipid and energy, was the functional category most affected by diet. Important effects were also measured in ribosomal proteins and signalling. The long-chain polyunsaturated fatty acid (LC-PUFA) biosynthesis pathway, assessed by fatty acid composition and gene expression, was influenced by genotype. Intestinal tissue contents of docosahexaenoic acid were equivalent in Lean salmon fed either a FO or VO diet and expression of LC-PUFA biosynthesis genes was up-regulated in VO-fed fish in Fat salmon. Dietary VO increased lipogenesis in Lean fish, assessed by expression of FAS, while no effect was observed on β-oxidation although transcripts of the mitochondrial respiratory chain were down-regulated, suggesting less active energetic metabolism in fish fed VO. In contrast, dietary VO up-regulated genes and proteins involved in detoxification, antioxidant defence and apoptosis, which could be associated with higher levels of polycyclic aromatic hydrocarbons in this diet. Regarding genotype, the following pathways were identified as being differentially affected: proteasomal proteolysis, response to oxidative and cellular stress (xenobiotic and oxidant metabolism and heat shock proteins), apoptosis and structural proteins particularly associated with tissue contractile properties. Genotype effects were accentuated by dietary VO.

Conclusions

Intestinal metabolism was affected by diet and genotype. Lean fish may have higher responsiveness to low dietary n-3 LC-PUFA, up-regulating the biosynthetic pathway when fed dietary VO. As global aquaculture searches for alternative oils for feeds, this study alerts to the potential of VO introducing contaminants and demonstrates the detoxifying role of intestine. Finally, data indicate genotype-specific responses in the intestinal transcriptome and proteome to dietary VO, including possibly structural properties of the intestinal layer and defence against cellular stress, with Lean fish being more susceptible to diet-induced oxidative stress.