MHC and Evolution in Teleosts

Multiomics Analyses Explore the Immunometabolic Interplay in the Liver of White Crucian Carp (Carassius cuvieri) After Aeromonas veronii Challenge

Aeromonas veronii is one of the predominant pathogenic species that can imperil the survival of farmed fish. However, the interactive networks of immune regulation and metabolic response in A. veronii-infected fish are still unclear. In this investigation, we aimed to explore immunometabolic interplay in white crucian carp (WCC) after the A. veronii challenge. Elevated levels of immune-related genes were observed in various tissues after A. veronii infection, along with the sharp alteration of disease-related enzymatic activities. Besides, decreased levels of antioxidant status were observed in the liver, but most metabolic gene expressions increased dramatically. Multiomics analyses revealed that metabolic products of amino acids, such as formiminoglutamic acid (FIGLU), L-glutamate (L-Glu), and 4-hydroxyhippuric acid, were considered the crucial liver biomarkers in A. veronii-infected WCC. In addition, A. veronii infection may dysregulate endoplasmic reticulum (ER) function to affect the metabolic process of lipids, carbohydrates, and amino acids in the liver of WCC. These results may have a comprehensive implication for understanding immunometabolic response in WCC upon A. veronii infection.

A nanocarrier immersion vaccine encoding surface immunogenic protein confers cross-immunoprotection against Streptococcus agalactiae and Streptococcus iniae infection in tilapia

Streptococcosis is a highly contagious aquatic bacterial disease that poses a significant threat to tilapia. Vaccination is a well-known effective measure to prevent and control fish bacterial diseases. Among the various immunization methods, immersion vaccination is simple and can be widely used in aquaculture. Besides, nanocarrier delivery technology has been reported as an effective solution to improve the immune effect of immersion vaccine. In this study, the surface immunogenic protein (Sip) was proved to be conserved and potential to provide cross-immunoprotection for both Streptococcus agalactiae (S. agalactiae) and Streptococcus iniae (S. iniae) by multiple sequences alignment and Western blotting analysis. On this basis, we expressed and obtained the recombinant protein rSip and connected it with functionalized carbon nanotubes (CNT) to construct the nanocarrier vaccine system CNT-rSip. After immersion immunization, the immune effect of CNT-rSip against above two streptococcus infections was evaluated in tilapia based on some aspects including the serum specific antibody level, non-specific enzyme activities, immune-related genes expression and relative percent survival (RPS) after bacteria challenge. The results showed that compared with control group, CNT-rSip significantly (P < 0.05) increased the serum antibody levels, related enzyme activities including acid phosphatase, alkaline phosphatase, lysozyme and total antioxidant capacity activities, as well as the expression levels of immune-related genes from 2 to 4 weeks post immunization (wpi), and all these indexes peaked at 3 wpi. Besides, the above indexes of CNT-rSip were higher than those of rSip group with different extend during the experiment. Furthermore, the challenge test indicated that CNT-rSip provided cross-immunoprotection against S. agalactiae and S. iniae infection with RPS of 75 % and 72.41 %, respectively, which were much higher than those of other groups. Our study indicated that the nanocarrier immersion vaccine CNT-rSip could significantly improve the antibody titer and confer cross-immuneprotection against S. agalactiae and S. iniae infection in tilapia.

Expression and Immune Characterization of Major Histocompatibility Complex in Paralichthys olivaceus after Antigen Stimulation

The Major histocompatibility complex (Mhc) is an important molecule for antigen presenting and binds to T cell receptors, activating T lymphocytes and triggering specific immune responses. To investigate the role of MhcII in adaptive immunity, in this study, mhcIIα and mhcIIβ of flounder (Paralichthys olivaceus) were cloned, polyclonal antibodies (Abs) against their extracellular regions were produced, respectively, and their distribution on cells and tissues and expression patterns, which varied by antigen stimulation or pathogen infection, were investigated. The results showed that the open reading frame (ORF) of mhcIIα is 708 bp, including 235 amino acids (aa); and the ORF of mhcIIβ is 741 bp, encoding 246aa. The mhcIIα and mhcIIβ were significantly expressed in gills, spleen, and peripheral blood leukocytes (PBLs). Their antibodies could specifically recognize eukaryotic expressed MhcIIα and MhcIIβ. MhcIIα+ and MhcIIβ+ cells were 30.2 ± 2.9% of the percentage in peripheral blood leukocytes. MhcII molecules were co-localized with CD83 and IgM on leukocytes, respectively, but not on CD4+ or CD8+ T lymphocyte subpopulations. The expression of both mhcIIα and mhcIIβ were significantly upregulated in flounder after bacteria and virus challenges. The percentages of MhcII+ cells, MhcII+/CD83+, and MhcII+/IgM+ double-positive cells increased significantly after PHA and ConA stimulation, respectively; they varied significantly in PBLs after polyI:C stimulation, and no variations were found after LPS treatment. In the meantime, variations in MhcII+ cells were consistent with that of CD4+ T lymphocytes. These results suggest that MhcII, mainly expressed in B cells and dendritic cells, play an essential role in antigen presentation, and respond significantly to exogenous antigens and T cell-dependent antigens. These results may provide an important reference for the study of cellular immunity in teleosts.

Major histocompatibility complex class IIα and IIβ of pufferfish (Takifugu obscurus): Identification and functional characterization in response to bacterial challenge

Classical major histocompatibility complex (MHC) class II molecules play an essential role in immune system. In this study, MHC IIα (Pf-MHC IIα) and MHC IIβ (Pf-MHC IIβ) homology genes from pufferfish (Takifugu obscurus) were cloned and their functional characterization in response to bacterial challenge was identified. The nucleotide sequences of the open reading frames (ORFs) of pufferfish Pf-MHC IIα and Pf-MHC IIβ were 708 bp and 750 bp, encoding 235 aa and 249 aa, respectively. The structure of Pf-MHC IIα or Pf-MHC IIβ contained a signal peptide, an α1/β1 domain, an α2/β2 domain, a transmembrane region and a cytoplasmic region. Multiple sequence alignment and phylogenetic analysis showed that Pf-MHC IIα and Pf-MHC IIβ molecules had the highest similarity with Fugu rubripes (Takifugu rubripes). Cellular localization analysis indicated that the distribution of Pf-MHC IIα and Pf-MHC IIβ was in the lymphocyte membrane and cytoplasm. qRT-PCR results showed that Pf-MHC IIα and Pf-MHC IIβ expressed relatively high in skin, gills and gut. In addition, after stimulation challenge in vitro (lipopolysaccharide, or polyinosinic: polycytidylic acid) and in vivo (A. hydrophila), the mRNA expressions of Pf-MHC IIα and Pf-MHC IIβ were significantly up-regulated in lymphocytes and in tissues of skin, gills, gut and head kidney. Moreover, Pf-MHC IIα or Pf-MHC IIβ neutralization reduced the ability of A. hydrophila to induce the expressions of lymphocyte cytokines (TNF-α, IL-1β and IL-10). Overall, it is speculated that Pf-MHC IIα and Pf-MHC IIβ may play an important role in the host response against A. hydrophila in pufferfish.

Exploring the Immunoprotective Potential of a Nanocarrier Immersion Vaccine Encoding Sip against Streptococcus Infection in Tilapia (Oreochromis niloticus)

Tilapia, as one of the fish widely cultured around the world, is suffering severe impact from the streptococcus disease with the deterioration of the breeding environment and the increasing of breeding density, which brings serious economic loss to tilapia farming. In this study, the surface immunogenic protein (Sip) of Streptococcus agalactiae (S. agalactiae) was selected as the potential candidate antigen and connected with bacterial nano cellulose (BNC) to construct the nanocarrier subunit vaccine (BNC-rSip), and the immersion immune effects against S. agalactiae and Streptococcus iniae (S. iniae) in Nile tilapia were evaluated on the basis of the serum antibody level, non-specific enzyme activity, the immune-related gene expression and relative percent survival (RPS). The results indicated that Sip possessed the expected immunogenicity according to the immunoinformatic analysis. Compared with the rSip group, BNC-rSip significantly induced serum antibody production and improved the innate immunity level of tilapia. After challenge, the RPS of BNC-rSip groups were 78.95% (S. agalactiae) and 67.86% (S. iniae), which were both higher than those of rSip groups,31.58% (S. agalactiae) and 35.71% (S. iniae), respectively. Our study indicated that BNC-rSip can induce protective immunity for tilapia through immersion immunization and may be an ideal candidate vaccine for controlling tilapia streptococcal disease.

Characteristics and expression profiles of MHC class Ⅰ molecules in Carassius auratus

Major histocompatibility complex class Ⅰ (MHC Ⅰ) molecules play a vital role in adaptive immune systems in vertebrates by presenting antigens to effector T cells. Understanding the expression profiling of MHC Ⅰ molecules in fish is essential for improving our knowledge of the relationship between microbial infection and adaptive immunity. In this study, we conducted a comprehensive analysis of MHC Ⅰ gene characteristics in Carassius auratus, an important freshwater aquaculture fish in China that is susceptible to Cyprinid herpesvirus 2 (CyHV-2) infection. We identified approximately 20 MHC Ⅰ genes discussed, including U, Z, and L lineage genes. However, only U and Z lineage proteins were identified in the kidney of Carassius auratus using high pH reversed-phase chromatography and mass spectrometry. The L lineage proteins were either not expressed or present at an extremely low level in the kidneys of Carassius auratus. We also used targeted proteomics to analyze changes in protein MHC Ⅰ molecules abundance in healthy and CyHV-2-infected Carassius auratus. We observed that five MHC Ⅰ molecules were upregulated, and Caau-UFA was downregulated in the diseased group. This study is the first to reveal the expression of MHC Ⅰ molecules at a large scale in Cyprinids, which enhances our understanding of fish adaptive immune systems.

Understanding the Pathogenesis of Red Mark Syndrome in Rainbow Trout (Oncorhynchus mykiss) through an Integrated Morphological and Molecular Approach

Red mark syndrome (RMS) is a widespread skin disorder of rainbow trout in freshwater aquaculture, believed to be caused by a Midichloria-like organism (MLO). Here, we aimed to study the pathologic mechanisms at the origin of RMS by analyzing field samples from a recent outbreak through gene expression, MLO PCR, quantitative PCR, and a histopathological scoring system proposed for RMS lesions. Statistical analyses included a One-Way Analysis of Variance (ANOVA) with a Dunnett's multiple comparisons test to assess differences among gene expression groups and a nonparametric Spearman correlation between various categories of skin lesions and PCR results. In short, the results confirmed the presence of a high quantity of 16S gene copy numbers of Midichloria-like organisms in diseased skin tissues. However, the number of Midichloria-like organisms detected was not correlated to the degree of severity of skin disease. Midichloria-like organism DNA was found in the spleen and head kidney. The spleen showed pathologic changes mainly of hyperplastic type, reflecting its direct involvement during infection. The most severe skin lesions were characterized by a high level of inflammatory cytokines sustaining and modulating the severe inflammatory process. IL-1 β, IL-6, IL-10, MHC-II, and TCR were upregulated in severe skin lesions, while IL-10 was highly expressed in moderate to severe ones. In the moderate form, the response was driven to produce immunoglobulins, which appeared crucial in controlling the skin disease's severity. Altogether our results illustrated a complex immune interaction between the host and Midichloria-like organism.

Meta-analysis of major histocompatibility complex (MHC) class IIA reveals polymorphism and positive selection in many vertebrate species

Pathogen-mediated selection and sexual selection are important drivers of evolution. Both processes are known to target genes of the major histocompatibility complex (MHC), a gene family encoding cell-surface proteins that display pathogen peptides to the immune system. The MHC is also a model for understanding processes such as gene duplication and trans-species allele sharing. The class II MHC protein is a heterodimer whose peptide-binding groove is encoded by an MHC-IIA gene and an MHC-IIB gene. However, our literature review found that class II MHC papers on infectious disease or sexual selection included IIA data only 18% and 9% of the time, respectively. To assess whether greater emphasis on MHC-IIA is warranted, we analysed MHC-IIA sequence data from 50 species of vertebrates (fish, amphibians, birds, mammals) to test for polymorphism and positive selection. We found that the number of MHC-IIA alleles within a species was often high, and covaried with sample size and number of MHC-IIA genes assayed. While MHC-IIA variability tended to be lower than that of MHC-IIB, the difference was only ~25%, with ~3 fewer IIA alleles than IIB. Furthermore, the unexpectedly high MHC-IIA variability showed clear signatures of positive selection in most species, and positive selection on MHC-IIA was stronger in fish than in other surveyed vertebrate groups. Our findings underscore that MHC-IIA can be an important target of selection. Future studies should therefore expand the characterization of MHC-IIA at both allelic and genomic scales, and incorporate MHC-IIA into models of fitness consequences of MHC variation.

Protective efficacy evaluation of immunogenic protein AHA_3793 of Aeromonas hydrophila as vaccine candidate for largemouth bass Micropterus salmoides

Aeromonas hydrophila is a Gram-negative pathogen that can infect various fish, including largemouth bass (Micropterus salmoides), which have caused huge economic losses. In present study, largemouth bass anti-A. hydrophila antibodies were produced, then a highly immunogenic outer membrane proteins, AHA_3793, was identified by combined western blotting and mass spectrometry analysis. Moreover, AHA_3793 was expressed, and its immunogenicity was further verified by western blotting. Subsequently, the protective efficacy of AHA_3793 were evaluated in largemouth bass. The results showed that rAHA_3793 could produce a relative percentage survival (RPS) of 61.76% for largemouth bass against A. hydrophila challenge. ELISA analysis showed the specific serum antibodies of largemouth bass against rAHA_3793 and A. hydrophila in vaccinated group in weeks 4 and 5 after immunization were significantly higher than those in control group, which suggested that rAHA_3793 induced production of specific serum antibodies against rAHA_3793 and A. hydrophila. The qRT-PCR analysis showed that expressions of CD4-2 and MHC IIα were also significantly up-regulated after immunization. These results collectively demonstrated that rAHA_3793 could induce a strong humoral immune response of largemouth bass, and then produce high immune protection effects against A. hydrophila infection.

Application of reverse vaccinology to design a multi-epitope subunit vaccine against a new strain of Aeromonas veronii

Background

Aeromonas veronii is one of the most common pathogens of freshwater fishes that cause sepsis and ulcers. There are increasing numbers of cases showing that it is a significant zoonotic and aquatic agent. Epidemiological studies have shown that A. veronii virulence and drug tolerance have both increased over the last few years as a result of epidemiological investigations. Cadaverine reverse transporter (CadB) and maltoporin (LamB protein) contribute to the virulence of A. veronii TH0426. TH0426 strain is currently showing severe cases on fish species, and its resistance against therapeutic has been increasing. Despite these devastating complications, there is still no effective cure or vaccine for this strain of A.veronii.

Results

In this regard, an immunoinformatic method was used to generate an epitope-based vaccine against this pathogen. The immunodominant epitopes were identified using the CadB and LamB protein of A. veronii. The final constructed vaccine sequence was developed to be immunogenic, non-allergenic as well as have better solubility. Molecular dynamic simulation revealed significant binding stability and structural compactness. Finally, using Escherichia coli K12 as a model, codon optimization yielded ideal GC content and a higher CAI value, which was then included in the cloning vector pET2+ (a).

Conclusion

Altogether, our outcomes imply that the proposed peptide vaccine might be a good option for A. veronii TH0426 prophylaxis.

Some thoughts about what non-mammalian jawed vertebrates are telling us about antigen processing and peptide loading of MHC molecules

The major histocompatibility complex (MHC) of mammals encodes highly polymorphic classical class I and class II molecules with crucial roles in immune responses, as well as various nonclassical molecules encoded by the MHC and elsewhere in the genome that have a variety of functions. These MHC molecules are supported by antigen processing and peptide loading pathways which are well-understood in mammals. This review considers what has been learned about the MHC, MHC molecules and the supporting pathways in non-mammalian jawed vertebrates. From the initial understanding from work with the chicken MHC, a great deal of diversity in the structure and function has been found. Are there underlying principles?

•

The genomic organisation of the MHC varies enormously across jawed vertebrates.

•

Total numbers of MHC genes vary among vertebrates, with only a few classical MHC genes.

•

Some nonclassical MHC and classical pathway genes appear earlier than others.

•

Obvious co-evolution within MHC pathways occurs in some species, but not others.

•

The promiscuity of interactions may correlate with differences in genomic organisation.

Confocal Characterization of Intestinal Dendritic Cells from Myxines to Teleosts

Dendritic cells (DCs) are antigen-presenting cells (APCs) that regulate the beginning of adaptive immune responses. The mechanisms of tolerance to antigens moving through the digestive tract are known to be regulated by intestinal DCs. Agnatha and Gnathostoma are descendants of a common ancestor. The Ostracoderms gave rise to Cyclostomes, whereas the Placoderms gave rise to Chondrichthyes. Sarcopterygii and Actinopterygii are two evolutionary lines of bony fishes. Brachiopterygii and Neopterygii descend from the Actinopterygii. From Neopterygii, Holostei and Teleostei evolved. Using immunohistochemistry with TLR-2, Langerin/CD207, and MHC II, this study aimed to characterize intestinal DCs, from myxines to teleosts. The findings reveal that DCs are positive for the antibodies tested, highlighting the presence of DCs and DC-like cells phylogenetically from myxines, for the first time, to teleosts. These findings may aid in improving the level of knowledge about the immune system’s evolution and these sentinel cells, which are crucial to the body’s defense.

Comparative assessment of immunological tolerance in fish with natural immunodeficiency

Natural occurrences of immunodeficiency by definition should lead to compromised immune function. The major histocompatibility complexes (MHC) are key components of the vertebrate adaptive immune system, charged with mediating allorecognition and antigen presentation functions. To this end, the genomic loss of the MHC II pathway in Syngnathus pipefishes raises questions regarding their immunological vigilance and allorecognition capabilities. Utilising allograft and autograft fin-transplants, we compared the allorecognition immune responses of two pipefish species, with (Nerophis ophidion) and without (Syngnathus typhle) a functional MHC II. Transcriptome-wide assessments explored the immunological tolerance and potential compensatory measures occupying the role of the absent MHC II. Visual observations suggested a more acute rejection response in N. ophidion allografts compared with S. typhle allografts. Differentially expressed genes involved in innate immunity, angiogenesis and tissue recovery were identified among transplantees. The intriguing upregulation of the cytotoxic T-cell implicated gzma in S. typhle allografts, suggests a prominent MHC I related response, which may compensate for the MHC II and CD4 loss. MHC I related downregulation in N. ophidion autografts hints at an immunological tolerance related reaction. These findings may indicate alternative measures evolved to cope with the MHC II genomic loss enabling the maintenance of appropriate tolerance levels. This study provides intriguing insights into the immune and tissue recovery mechanisms associated with syngnathid transplantation, and can be a useful reference for future studies focusing on transplantation transcriptomics in non-model systems.

Use of Integrated Core Proteomics, Immuno-Informatics, and In Silico Approaches to Design a Multiepitope Vaccine against Zoonotic Pathogen Edwardsiella tarda

Multidrug-resistant Edwardsiella tarda has been reported as the main causative agent for massive fish mortality. The pathogen is well-known for causing hemorrhagic septicemia in fish and has been linked to gastrointestinal infections in humans. Formalin-inactivated Edwardsiella vaccination has previously been found to be ineffective in aquaculture species. Therefore, based on E. tarda’s integrated core complete sequenced genomes, the study aimed to design a subunit vaccine based on T and B cell epitopes employing immunoinformatics approach. Initially, the top immunodominant and antigenic epitopes were predicted from the core complete sequenced genomes of the E. tarda genome and designed the vaccine by using linkers and adjuvant. In addition, vaccine 3D structure was predicted followed by refinement, and molecular docking was performed for the analysis of interacting residues between vaccines with TLR5, MHC-I, and MHC-II, respectively. The final vaccine constructs demonstrated strong hydrogen bond interactions. Molecular dynamic simulation of vaccine-TLR5 receptor complex showed a stable structural binding and compactness. Furthermore, E. coli used as a model organism for codon optimization proved optimal GC content and CAI value, which were subsequently cloned in vector pET2+ (a). Overall, the findings of the study imply that the designed epitope vaccine might be a good option for prophylaxis for E. tarda.

Intra-specific copy number variation of MHC class II genes in the Siamese fighting fish

Duplicates of genes for major histocompatibility complex (MHC) molecules can be subjected to selection independently and vary markedly in their evolutionary rates, sequence polymorphism, and functional roles. Therefore, without a thorough understanding of their copy number variation (CNV) in the genome, the MHC-dependent fitness consequences within a species could be misinterpreted. Studying the intra-specific CNV of this highly polymorphic gene, however, has long been hindered by the difficulties in assigning alleles to loci and the lack of high-quality genomic data. Here, using the high-quality genome of the Siamese fighting fish (Betta splendens), a model for mate choice studies, and the whole-genome sequencing (WGS) data of 17 Betta species, we achieved locus-specific amplification of their three classical MHC class II genes - DAB1, DAB2, and DAB3. By performing quantitative PCR and depth-of-coverage analysis using the WGS data, we revealed intra-specific CNV at the DAB3 locus. We identified individuals that had two allelic copies (i.e., heterozygous or homozygous) or one allele (i.e., hemizygous) and individuals without this gene. The CNV was due to the deletion of a 20-kb-long genomic region harboring both the DAA3 and DAB3 genes. We further showed that the three DAB genes were under different modes of selection, which also applies to their corresponding DAA genes that share similar pattern of polymorphism. Our study demonstrates a combined approach to study CNV within a species, which is crucial for the understanding of multigene family evolution and the fitness consequences of CNV.

Designing a novel mRNA vaccine against Vibrio harveyi infection in fish: an immunoinformatics approach

Vibrio harveyi belongs to the family Vibrionaceae of class Gammaproteobacteria. Around 12 Vibrio species can cause gastroenteritis (gastrointestinal illness) in humans. A large number of bacterial particles can be found in the infected cells, which may cause death. Despite these devastating complications, there is still no cure or vaccine for the bacteria. As a result, we used an immunoinformatics approach to develop a multi-epitope vaccine against the most pathogenic hemolysin gene of V. harveyi. The immunodominant T- and B-cell epitopes were identified using the hemolysin protein. We developed a vaccine employing three possible epitopes: cytotoxic T-lymphocytes, helper T-lymphocytes, and linear B-lymphocyte epitopes, after thorough testing. The vaccine was developed to be antigenic, immunogenic, and non-allergenic, as well as have a better solubility. Molecular dynamics simulation revealed significant structural stiffness and binding stability. In addition, the immunological simulation generated by computers revealed that the vaccination might elicit immune reactions Escherichia coli K12 as a model, codon optimization yielded ideal GC content and a higher codon adaptation index value, which was then included in the cloning vector pET2+ (a). Altogether, our experiment implies that the proposed peptide vaccine might be a good option for vibriosis prophylaxis.

Why Does Piscirickettsia salmonis Break the Immunological Paradigm in Farmed Salmon? Biological Context to Understand the Relative Control of Piscirickettsiosis

Piscirickettsiosis (SRS) has been the most important infectious disease in Chilean salmon farming since the 1980s. It was one of the first to be described, and to date, it continues to be the main infectious cause of mortality. How can we better understand the epidemiological situation of SRS? The catch-all answer is that the Chilean salmon farming industry must fight year after year against a multifactorial disease, and apparently only the environment in Chile seems to favor the presence and persistence of Piscirickettsia salmonis. This is a fastidious, facultative intracellular bacterium that replicates in the host’s own immune cells and antigen-presenting cells and evades the adaptive cell-mediated immune response, which is why the existing vaccines are not effective in controlling it. Therefore, the Chilean salmon farming industry uses a lot of antibiotics—to control SRS—because otherwise, fish health and welfare would be significantly impaired, and a significantly higher volume of biomass would be lost per year. How can the ever-present risk of negative consequences of antibiotic use in salmon farming be balanced with the productive and economic viability of an animal production industry, as well as with the care of the aquatic environment and public health and with the sustainability of the industry? The answer that is easy, but no less true, is that we must know the enemy and how it interacts with its host. Much knowledge has been generated using this line of inquiry, however it remains insufficient. Considering the state-of-the-art summarized in this review, it can be stated that, from the point of view of fish immunology and vaccinology, we are quite far from reaching an effective and long-term solution for the control of SRS. For this reason, the aim of this critical review is to comprehensively discuss the current knowledge on the interaction between the bacteria and the host to promote the generation of more and better measures for the prevention and control of SRS.

Stress response gene family expansions correlate with invasive potential in teleost fish

The bluegill sunfish Lepomis macrochirus and the closely related redear sunfish Lepomis microlophus have important ecological and recreational value and are widely used for research and aquaculture. While both species have been introduced outside of their native ranges, only the bluegill is considered invasive. Here, we report de novo transcriptome assemblies for these fish as a resource for sunfish biology. Comparative analyses of the transcriptomes revealed an unexpected, bluegill-specific expansion in the HSP70 and HSP90 molecular chaperone gene families. These expansions were not unique to the bluegill as expansions in HSP70s and HSP90s were identified in the genomes of other teleost fish using the NCBI RefSeq database. To determine whether gene family expansions are specific for thermal stress responses, GST and SOD gene families that are associated with oxidative stress responses were also analyzed. Species-specific expansions were also observed for these gene families in distinct fish species. Validating our approach, previously described expansions in the MHC gene family were also identified. Intriguingly, the number of HSP70 paralogs was positively correlated with thermotolerance range for each species, suggesting that these expansions can impact organismal physiology. Furthermore, fish that are considered invasive contained a higher average number of HSP70 paralogs than non-invasive fish. Invasive fish also had higher average numbers of HSP90, MHC and GST paralogs, but not SOD paralogs. Taken together, we propose that expansions in key cellular stress response gene families represent novel genetic signatures that correlate with invasive potential.

Comprehensive transcriptome analysis reveals the effect of feeding rhythm on the immunity and metabolism of Acipenser dabryanus

Acipenser dabryanus is a species endemic to Yangtze River drainage in China and is listed as a critical endangered species on the IUCN Red List. In the present study, the liver and spleen transcriptomes were analyzed by comparing the data of A. dabryanus that experienced nine different feeding rhythms (once a day diurnal, T1; two times a day diurnal, T2; three times a day diurnal, T3; four times a day, T4; five times a day, T5; six times a day, T6; once a day nocturnal, Tn1; two times a day nocturnal, Tn2; and three times a day nocturnal, Tn3). Transcriptome sequencing generated 1,901,236,482 clean reads, encompassing 570.4 Gb of sequence data. The reads were assembled into 287,372 unigenes with an average length of 803 bp and an N50 of 1004 bp. KEGG analysis showed that 1,080, 1,030, and 1216 unigenes were annotated to lipid metabolism, amino acid metabolism and carbohydrate metabolism, respectively, and 2549 unigenes were annotated to the immune system category. Differentially expressed genes (DEGs) between different feeding frequency groups or between nighttime and daytime feeding were obtained and functionally enriched. Importantly, DEGs participating in nutrition metabolism and various immunoregulation pathways and their expression profiles in A. dabryanus were discussed. Interestingly, the majority of key genes related to lipid metabolism or in immunodependent gene families, such as antimicrobial peptides, Toll-like receptors, chemokines, NOD-like receptors, B cell receptors and the major histocompatibility complex, were all significantly upregulated in animals in the T6 group compared to the characteristics of animals in the T2 group that had a normal feeding frequency. In addition, light/dark rhythm also affected the immunity of A. dabryanus, and fish fed at night possessed an improved immune response than fish fed at daytime. Our study suggested that feeding six times a day is optimal for A. dabryanus juvenile growth as it enhances the organism's nutrition metabolism and immune function.

Transcriptome Profiling of Atlantic Salmon (Salmo salar) Parr With Higher and Lower Pathogen Loads Following Piscirickettsia salmonis Infection

Salmonid rickettsial septicemia (SRS), caused by Piscirickettsia salmonis, is one of the most devastating diseases of salmonids. However, the transcriptomic responses of Atlantic salmon (Salmon salar) in freshwater to an EM-90-like isolate have not been explored. Here, we infected Atlantic salmon parr with an EM-90-like isolate and conducted time-course qPCR analyses of pathogen load and four biomarkers (campb, hampa, il8a, tlr5a) of innate immunity on the head kidney samples. Transcript expression of three of these genes (except hampa), as well as pathogen level, peaked at 21 days post-injection (DPI). Multivariate analyses of infected individuals at 21 DPI revealed two infection phenotypes [lower (L-SRS) and higher (H-SRS) infection level]. Five fish from each group (Control, L-SRS, and H-SRS) were selected for transcriptome profiling using a 44K salmonid microarray platform. We identified 1,636 and 3,076 differentially expressed probes (DEPs) in the L-SRS and H-SRS groups compared with the control group, respectively (FDR = 1%). Gene ontology term enrichment analyses of SRS-responsive genes revealed the activation of a large number of innate (e.g. “phagocytosis”, “defense response to bacterium”, “inflammatory response”) and adaptive (e.g. “regulation of T cell activation”, “antigen processing and presentation of exogenous antigen”) immune processes, while a small number of general physiological processes (e.g. “apoptotic process”, development and metabolism relevant) was enriched. Transcriptome results were confirmed by qPCR analyses of 42 microarray-identified transcripts. Furthermore, the comparison of individuals with differing levels of infection (H-SRS vs. L-SRS) generated insights into the biological processes possibly involved in disease resistance or susceptibility. This study demonstrated a low mortality (~30%) EM-90-like infection model and broadened the current understanding of molecular pathways underlying P. salmonis-triggered responses of Atlantic salmon, identifying biomarkers that may assist to diagnose and combat this pathogen.

Histocompatibility and Reproduction: Lessons from the Anglerfish

Reproduction in certain deep-sea anglerfishes involves the permanent attachment of dwarf males to much larger females and fusion of their tissues leading to the establishment of a shared circulatory system. This unusual phenomenon of sexual parasitism enables anglerfishes to maximize reproductive success in the vast and deep oceans, where females and males otherwise rarely meet. An even more surprising phenomenon relates to the observation that joining of genetically disparate male and female anglerfishes does not evoke a strong anti-graft immune rejection response, which occurs in vertebrates following allogeneic parabiosis. Recent studies demonstrated that the evolutionary processes that led to the unique mating strategy of anglerfishes coevolved with genetic changes that resulted in loss of functional genes encoding critical components of the adaptive immune system. These genetic alterations enabled anglerfishes to tolerate the histoincompatible tissue antigens of their mate and prevent the occurrence of reciprocal graft rejection responses. While the exact mechanisms by which anglerfishes defend themselves against pathogens have not yet been deciphered, it is speculated that during evolution, anglerfishes adopted new immune strategies that compensate for the loss of B and T lymphocyte functions and enable them to resist infection by pathogens.

The pink salmon genome: Uncovering the genomic consequences of a two-year life cycle

Pink salmon (Oncorhynchus gorbuscha) adults are the smallest of the five Pacific salmon native to the western Pacific Ocean. Pink salmon are also the most abundant of these species and account for a large proportion of the commercial value of the salmon fishery worldwide. A two-year life history of pink salmon generates temporally isolated populations that spawn either in even-years or odd-years. To uncover the influence of this genetic isolation, reference genome assemblies were generated for each year-class and whole genome re-sequencing data was collected from salmon of both year-classes. The salmon were sampled from six Canadian rivers and one Japanese river. At multiple centromeres we identified peaks of Fst between year-classes that were millions of base-pairs long. The largest Fst peak was also associated with a million base-pair chromosomal polymorphism found in the odd-year genome near a centromere. These Fst peaks may be the result of a centromere drive or a combination of reduced recombination and genetic drift, and they could influence speciation. Other regions of the genome influenced by odd-year and even-year temporal isolation and tentatively under selection were mostly associated with genes related to immune function, organ development/maintenance, and behaviour.

The bowfin genome illuminates the developmental evolution of ray-finned fishes

The bowfin (Amia calva) is a ray-finned fish that possesses a unique suite of ancestral and derived phenotypes, which are key to understanding vertebrate evolution. The phylogenetic position of bowfin as a representative of neopterygian fishes, its archetypical body plan and its unduplicated and slowly evolving genome make bowfin a central species for the genomic exploration of ray-finned fishes. Here we present a chromosome-level genome assembly for bowfin that enables gene-order analyses, settling long-debated neopterygian phylogenetic relationships. We examine chromatin accessibility and gene expression through bowfin development to investigate the evolution of immune, scale, respiratory and fin skeletal systems and identify hundreds of gene-regulatory loci conserved across vertebrates. These resources connect developmental evolution among bony fishes, further highlighting the bowfin's importance for illuminating vertebrate biology and diversity in the genomic era.

The Spleen as a Target to Characterize Immunomodulatory Effects of Down-Stream Processed Cyberlindnera jadinii Yeasts in Atlantic Salmon Exposed to a Dietary Soybean Meal Challenge

Aquaculture feeds have changed dramatically from being largely based on fishmeal (FM) towards increased use of plant protein sources, which could impact the fish's immune response. In order to characterize immunomodulatory properties of novel functional ingredients, this study used four diets, one based on FM, a challenging diet with 40% soybean meal (SBM), and two diets containing 40% SBM with 5% of Cyberlindnera jadinii yeast exposed to different down-stream processing conditions: heat-inactivated (ICJ) or autolysation (ACJ). The immunomodulatory effects of the diets were analyzed in the spleen of Atlantic salmon after 37 days of feeding, using a transcriptomic evaluation by RNA sequencing (RNA-seq) and the detection of specific immunological markers at the protein level through indirect Enzyme-linked Immunosorbent Assay (indirect ELISA). The results showed that SBM (compared to FM) induced a down-regulation of pathways related to ion binding and transport, along with an increase at the protein level of pro-inflammatory cytokines such as tumor necrosis factor alpha (TNFα) and interferon gamma (IFNγ). On the other hand, while ICJ (compared to FM-group) maintain the inflammatory response associated with SBM, with higher levels of TNFα and IFNγ, and with an upregulation of creatine kinase activity and phosphagen metabolic process, the inclusion of ACJ was able to modulate the response of Atlantic salmon compared to fish fed the SBM-diet by the activation of biological pathways related to endocytosis, Pattern recognition receptor (PPRs)-signal transduction and transporter activity. In addition, ACJ was also able to control the pro-inflammatory profile of SBM, increasing Interleukin 10 (IL-10) levels and decreasing TNFα production, triggering an immune response similar to that of fish fed an FM-based diet. Finally, we suggest that the spleen is a good candidate to characterize the immunomodulatory effects of functional ingredients in Atlantic salmon. Moreover, the inclusion of ACJ in fish diets, with the ability to control inflammatory processes, could be considered in the formulation of sustainable salmon feed.

Diversification of CD1 Molecules Shapes Lipid Antigen Selectivity

Molecular studies of host-pathogen evolution have largely focused on the consequences of variation at protein-protein interaction surfaces. The potential for other microbe-associated macromolecules to promote arms race dynamics with host factors remains unclear. The cluster of differentiation 1 (CD1) family of vertebrate cell surface receptors plays a crucial role in adaptive immunity through binding and presentation of lipid antigens to T-cells. Although CD1 proteins present a variety of endogenous and microbial lipids to various T-cell types, they are less diverse within vertebrate populations than the related major histocompatibility complex (MHC) molecules. We discovered that CD1 genes exhibit a high level of divergence between simian primate species, altering predicted lipid-binding properties and T-cell receptor interactions. These findings suggest that lipid-protein conflicts have shaped CD1 genetic variation during primate evolution. Consistent with this hypothesis, multiple primate CD1 family proteins exhibit signatures of repeated positive selection at surfaces impacting antigen presentation, binding pocket morphology, and T-cell receptor accessibility. Using a molecular modeling approach, we observe that interspecies variation as well as single mutations at rapidly-evolving sites in CD1a drastically alter predicted lipid binding and structural features of the T-cell recognition surface. We further show that alterations in both endogenous and microbial lipid-binding affinities influence the ability of CD1a to undergo antigen swapping required for T-cell activation. Together these findings establish lipid-protein interactions as a critical force of host-pathogen conflict and inform potential strategies for lipid-based vaccine development.

Interferon Gamma Induces the Increase of Cell-Surface Markers (CD80/86, CD83 and MHC-II) in Splenocytes From Atlantic Salmon

Type II interferon gamma (IFNγ) is a pleiotropic cytokine capable of modulating the innate and adaptive immune responses which has been widely characterized in several teleost families. In fish, IFNγ stimulates the expression of cytokines and chemokines associated with the pro-inflammatory response and enhances the production of nitrogen and oxygen reactive species in phagocytic cells. This work studied the effect of IFNγ on the expression of cell-surface markers on splenocytes of Atlantic salmon (Salmo salar). In vitro results showed that subpopulations of mononuclear splenocytes cultured for 15 days were capable of increasing gene expression and protein availability of cell-surface markers such as CD80/86, CD83 and MHC II, after being stimulated with recombinant IFNγ. These results were observed for subpopulations with characteristics associated with monocytes (51%), and features that could be related to lymphocytes (46.3%). In addition, a decrease in the expression of zbtb46 was detected in IFNγ-stimulated splenocytes. Finally, the expression of IFNγ and cell-surface markers was assessed in Atlantic salmon under field conditions. In vivo results showed that the expression of ifnγ increased simultaneously with the up-regulation of cd80/86, cd83 and mhcii during a natural outbreak of Piscirickettsia salmonis. Overall, the results obtained in this study allow us to propose IFNγ as a candidate molecule to stimulate the phenotypic progression of a small population of immune cells, which will increase antigen presenting cells markers. Thereby, modulatory strategies using IFNγ may generate a robust and coordinated immune response in fish against pathogens that affect aquaculture.

Induction of foxp3 during the Crosstalk between Antigen Presenting Like-Cells MHCII+CD83+ and Splenocytes CD4+IgM- in Rainbow Trout

Simple Summary

In aquatic biological models, the communication between cells from the immune system remains poorly characterized. In this work, to determine the gene expression of master transcriptional factors that coordinate the polarization of T cells, co-cultures of rainbow trout splenocytes are analyzed after stimulation with Interferon-gamma and/or Piscirickettsia salmonis. The results showed an upregulation of foxp3 compared to the other transcriptional factors, suggesting a potential communication between cells in the spleen, which may induce a Treg phenotype.

Abstract

In fish, the spleen is one of the major immune organs in the animal, and the splenocytes could play a key role in the activation and modulation of the immune response, both innate and adaptive. However, the crosstalk between different types of immune cells in the spleen has been poorly understood. In this work, an in vitro strategy is carried out to obtain and characterize mononuclear splenocytes from rainbow trout, using biomarkers associated with lymphocytes (CD4 and IgM) and antigen-presenting cells (CD83 and MHC II). Using these splenocytes, co-cultures of 24 and 48 h are used to determine the gene expression of master transcriptional factors that coordinate the polarization of T cells (t-bet, gata3, and foxp3). The results show a proportional upregulation of foxp3 (compared to t-bet and gata3) in co-cultures (at 24 h) of IFNγ-induced splenocytes with and without stimulation of Piscirickettsia salmonis proteins. In addition, foxp3 upregulation was established in co-cultures with IFNγ-induced cells and in cells only stimulated previously with P. salmonis proteins at 48 h of co-culture. These results show a potential communication between antigen-presenting-like cells and lymphocyte in the spleen, which could be induced towards a Treg phenotype.

Characterization of novel zebrafish MHC class I U lineage genes and their haplotype

Major histocompatibility complex (MHC) genes are essential for distinguishing between individuals in all jawed vertebrates. Although MHC class I (mhc1) genes in zebrafish comprise distinct haplotypes, not all members of the mhc1 gene family have been fully characterized. In this study, we report the identification of two novel U lineage genes isolated from the WIK strain of zebrafish. These new mhc1 genes, named una and uoa, are located in tandem on chromosome 19 with >70% homology to previously isolated U genes. Sequencing of their neighboring genes revealed that una and uoa form a unique haplotype different from the previously known U lineage haplotypes. Additionally, we determined the expression profiles of U, Z, and L genes in three different tissues. These findings collectively suggest that mhc1 U lineage genes and their haplotypes in zebrafish are more divergent than previously considered, and their expression patterns vary significantly among different tissues.

MHC class I evolution; from Northern pike to salmonids

Background

Salmonids are of major importance both as farmed and wild animals. With the changing environment comes changes in pathogenic pressures so understanding the immune system of all salmonid species is of essence. Major histocompatibility complex (MHC) genes are key players in the adaptive immune system signalling infection to responding T-cells populations. Classical MHC class I (MHCI) genes, defined by high polymorphism, broad expression patterns and peptide binding ability, have a key role in inducing immunity. In salmonids, the fourth whole genome duplication that occurred 94 million years ago has provided salmonids with duplicate MHCI regions, while Northern Pike, a basal sister clade to salmonids, represent a species which has not experienced this whole genome duplication.

Results

Comparing the gene organization and evolution of MHC class I gene sequences in Northern pike versus salmonids displays a complex picture of how many of these genes evolved. Regional salmonid Ia and Ib Z lineage gene duplicates are not orthologs to the Northern pike Z lineage sequences. Instead, salmonids have experienced unique gene duplications in both duplicate regions as well as in the Salmo and Oncorhynchus branch. Species-specific gene duplications are even more pronounced for some L lineage genes.

Conclusions

Although both Northern pike as well as salmonids have expanded their U and Z lineage genes, these gene duplications occurred separately in pike and in salmonids. However, the similarity between these duplications suggest the transposable machinery was present in a common ancestor. The salmonid MHCIa and MHCIb regions were formed during the 94 MYA since the split from pike and before the Oncorhynchus and Salmo branch separated. As seen in tetrapods, the non-classical U lineage genes are diversified duplicates of their classical counterpart. One MHCI lineage, the L lineage, experienced massive species-specific gene duplications after Oncorhynchus and Salmo split approximately 25 MYA. Based on what we currently know about L lineage genes, this large variation in number of L lineage genes also signals a large functional diversity in salmonids.

Fate of MHCII in salmonids following 4WGD

Major histocompatibility complex (MHC) genes are key players in the adaptive immunity providing a defense against invading pathogens. Although the basic structures are similar when comparing mammalian and teleost MHC class II (MHCII) molecules, there are also clear-cut differences. Based on structural requirements, the teleosts non-classical MHCII molecules do not comply with a function similar to the human HLA-DM and HLA-DO, i.e., assisting in peptide loading and editing of classical MHCII molecules. We have previously studied the evolution of teleost class II genes identifying various lineages and tracing their phylogenetic occurrence back to ancient ray-finned fishes. We found no syntenic MHCII regions shared between cyprinids, salmonids, and neoteleosts, suggesting regional instabilities. Salmonids have experienced a unique whole genome duplication 94 million years ago, providing them with the opportunity to experiment with gene duplicates. Many salmonid genomes have recently become available, and here we set out to investigate how MHCII has evolved in salmonids using Northern pike as a diploid sister phyla, that split from the salmonid lineage prior to the fourth whole genome duplication (4WGD) event. We identified 120 MHCII genes in pike and salmonids, ranging from 11 to 20 genes per species analyzed where DB-group genes had the most expansions. Comparing the MHC of Northern pike with that of Atlantic salmon and other salmonids species provides a tale of gene loss, translocations, and genome rearrangements.

Selective Stimulation of Duplicated Atlantic Salmon MHC Pathway Genes by Interferon-Gamma

Induction of cellular immune responses rely on Major histocompatibility complex (MHC) molecules presenting pathogenic peptides to T cells. Peptide processing, transport, loading and editing is a constitutive process in most cell types, but is accelerated upon infection. Recently, an unexpected complexity in the number of functional genes involved in MHC class I peptide cleavage, peptide transport, peptide loading and editing was found in teleosts, originating from the second and third whole genome duplication events. Salmonids have expanded upon this with functional duplicates also from a fourth unique salmonid whole genome duplication. However, little is known about how individual gene duplicates respond in the context of stimulation. Here we set out to investigate how interferon gamma (IFNg) regulates the transcription of immune genes in Atlantic salmon with particular focus on gene duplicates and MHC pathways. We identified a range of response patterns in Atlantic salmon gene duplicates, with upregulation of all duplicates for some genes, like interferon regulatory factor 1 (IRF1) and interferon induced protein 44-like (IFI44.L), but only induction of one or a few duplicates of other genes, such as TAPBP and ERAP2. A master regulator turned out to be the IRF1 and not the enhanceosome as seen in mammals. If IRF1 also collaborates with CIITA and possibly NLRC5 in regulating IFNg induction of MHCI and MHCII expression in Atlantic salmon, as in zebrafish, remains to be established. Altogether, our results show the importance of deciphering between gene duplicates, as they often respond very differently to stimulation and may have different biological functions.

The immunogenetics of sexual parasitism

Sexual parasitism has evolved as a distinctive mode of reproduction among deep-sea anglerfishes. The permanent attachment of males to host females observed in these species represents a form of anatomical joining, which is otherwise unknown in nature. Pronounced modifications to immune facilities are associated with this reproductive trait. The genomes of species with temporarily attaching males lack functional aicda genes that underpin affinity maturation of antibodies. Permanent attachment is associated with additional alterations, culminating in the loss of functional rag genes in some species, abolishing somatic diversification of antigen receptor genes, the hallmark of canonical adaptive immunity. In anglerfishes, coevolution of innate and adaptive immunity has been disentangled, implying that an alternative form of immunity supported the emergence of this evolutionarily successful group of vertebrates.

Extreme genomic volatility characterizes the evolution of the immunoglobulin heavy chain locus in cyprinodontiform fishes

The evolution of the adaptive immune system has provided vertebrates with a uniquely sophisticated immune toolkit, enabling them to mount precise immune responses against a staggeringly diverse range of antigens. Like other vertebrates, teleost fishes possess a complex and functional adaptive immune system; however, our knowledge of the complex antigen-receptor genes underlying its functionality has been restricted to a small number of experimental and agricultural species, preventing systematic investigation into how these crucial gene loci evolve. Here, we analyse the genomic structure of the immunoglobulin heavy chain (IGH) gene loci in the cyprinodontiforms, a diverse and important group of teleosts present in many different habitats across the world. We reconstruct the complete IGH loci of the turquoise killifish (Nothobranchius furzeri) and the southern platyfish (Xiphophorus maculatus) and analyse their in vivo gene expression, revealing the presence of species-specific splice isoforms of transmembrane IGHM. We further characterize the IGH constant regions of 10 additional cyprinodontiform species, including guppy, Amazon molly, mummichog and mangrove killifish. Phylogenetic analysis of these constant regions suggests multiple independent rounds of duplication and deletion of the teleost-specific antibody class IGHZ in the cyprinodontiform lineage, demonstrating the extreme volatility of IGH evolution. Focusing on the cyprinodontiforms as a model taxon for comparative evolutionary immunology, this work provides novel genomic resources for studying adaptive immunity and sheds light on the evolutionary history of the adaptive immune system.

The Mechanism of β2m Molecule-Induced Changes in the Peptide Presentation Profile in a Bony Fish

Contemporary antigen presentation knowledge is based on the existence of a single β2m locus, and a classical MHC class I forms a complex with a peptide (i.e., pMHC-I) to trigger CTL immunity. However, two β2m loci have been found in diploid bony fish; the function of the two β2m molecules is unclear. Here, we determined the variant peptide profiles originating from different products of the β2m loci binding to the same MHC-I molecule and further solved the crystal structures of the two pMHC-I molecules (i.e., pCtid-UAA-β2m-2 and pCtid-UAA-β2m-1-II). Of note, in pCtid-UAA-β2m-2, a unique hydrogen bond network formed in the bottom of the peptide-binding groove (PBG) led to α2-helix drift, ultimately leading to structural changes in the PBG. The mechanism of the change in peptide presentation profiles by β2m molecules is illustrated. The results are also of great significance for antivirus and antitumor functions in cold-blooded vertebrates and even humans.

Major Histocompatibility Complex in Osteichthyes

Based on analysis of available genome sequences, five gene lineages of MHC class I molecules (MHC I-U, -Z, -S, -L and -P) and one gene lineage of MHC class II molecules (MHC II-D) have been identified in Osteichthyes. In the latter lineage, three MHC II molecule sublineages have been identified (MHC II-A, -B and -E). As regards MHC class I molecules in Osteichthyes, it is important to take note of the fact that the lineages U and Z in MHC I genes have been identified in almost all fish species examined so far. Phylogenetic studies into MHC II molecule genes of sublineages A and B suggest that they may be descended from the genes of the sublineage named A/B that have been identified in spotted gar (Lepisosteus oculatus). The sublineage E genes of MHC II molecules, which represent the group of non-polymorphic genes with poor expression in the tissues connected with the immune system, are present in primitive fish, i.e. in paddlefish, sturgeons and spotted gar (Lepisosteus oculatus), as well as in cyprinids (Cyprinidae), Atlantic salmon (Salmo salar), and rainbow trout (Oncorhynchus mykiss). Full elucidation of the details relating to the organisation and functioning of the particular components of the major histocompatibility complex in Osteichthyes can advance the understanding of the evolution of the MHC molecule genes and the immune mechanism.

An Ancient, MHC-Linked, Nonclassical Class I Lineage in Cartilaginous Fish

Cartilaginous fishes, or chondrichthyans, are the oldest jawed vertebrates that have an adaptive immune system based on the MHC and Ig superfamily-based AgR. In this basal group of jawed vertebrates, we identified a third nonclassical MHC class I lineage (UDA), which is present in all species analyzed within the two major cartilaginous subclasses, Holocephali (chimaeras) and Elasmobranchii (sharks, skates, and rays). The deduced amino acid sequences of UDA have eight out of nine typically invariant residues that bind to the N and C termini of bound peptide found in most vertebrae classical class I (UAA); additionally, the other predicted 28 peptide-binding residues are perfectly conserved in all elasmobranch UDA sequences. UDA is distinct from UAA in its differential tissue distribution and its lower expression levels and is mono- or oligomorphic unlike the highly polymorphic UAA UDA has a low copy number in elasmobranchs but is multicopy in the holocephalan spotted ratfish (Hydrolagus colliei). Using a nurse shark (Ginglymostoma cirratum) family, we found that UDA is MHC linked but separable by recombination from the tightly linked cluster of UAA, TAP, and LMP genes, the so-called class I region found in most nonmammalian vertebrates. UDA has predicted structural features that are similar to certain nonclassical class I genes in other vertebrates, and, unlike polymorpic classical class I, we anticipate that it may bind to a conserved set of specialized peptides.

An Illumina approach to MHC typing of Atlantic salmon

The IPD-MHC Database represents the official repository for non-human major histocompatibility complex (MHC) sequences, overseen and supported by the Comparative MHC Nomenclature Committee, providing access to curated MHC data and associated analysis tools. IPD-MHC gathers allelic MHC class I and class II sequences from classical and non-classical MHC loci from various non-human animals including pets, farmed and experimental model animals. So far, Atlantic salmon and rainbow trout are the only teleost fish species with MHC class I and class II sequences present. For the remaining teleost or ray-finned species, data on alleles originating from given classical locus is scarce hampering their inclusion in the database. However, a fast expansion of sequenced genomes opens for identification of classical loci where high-throughput sequencing (HTS) will enable typing of allelic variants in a variety of new teleost or ray-finned species. HTS also opens for large-scale studies of salmonid MHC diversity challenging the current database nomenclature and analysis tools. Here we establish an Illumina approach to identify allelic MHC diversity in Atlantic salmon, using animals from an endangered wild population, and alter the salmonid MHC nomenclature to accommodate the expected sequence expansions.

The Reproducibility of an Inferred Tree and the Diploidization of Gene Segregation after Genome Duplication

We previously introduced a numerical quantity called the stability (Ps) of an inferred tree and showed that for the tree to be reliable this stability as well as the reliability of the tree, which is usually computed as the bootstrap probability (Pb), must be high. However, if genome duplication occurs in a species, a gene family of the genome also duplicates, and for this reason alone some Ps values can be high in a tree of the duplicated gene families. In addition, the topology of the duplicated gene family can be similar to that of the original gene family if such gene families are identifiable. After genome duplication, however, the gene families are often partially deleted or partially duplicated, and the duplicated gene family may not show the same topology as that of the original family. It is therefore necessary to compute the similarity of the topologies of the duplicated and the original gene families. In this paper, we introduce another quantity called the reproducibility (Pr) for measuring the similarity of the two gene families. To show how to compute the Pr values, we first compute the Pb and Ps values for each of the MHC class II α and β chain gene families, which were apparently generated by genome duplication. We then compute the Pr values for the MHC class II α and β chain gene families. The Pr values for the α and β chain gene families are now low, and this suggests that the diploidization of gene segregation has occurred after the genome duplication. Currently higher animals, defined as animals with complex phenotypic characters, generally have a higher genome size, and this increase in genome size appears to have been caused by genome duplication and diploidization of gene segregation after genome duplication.

Not all birds have a single dominantly expressed MHC-I gene: Transcription suggests that siskins have many highly expressed MHC-I genes

Passerine birds belong to the most species rich bird order and are found in a wide range of habitats. The extremely polymorphic adaptive immune system of passerines, identified through their major histocompatibility complex class I genes (MHC-I), may explain some of this extreme radiation. Recent work has shown that passerines have higher numbers of MHC-I gene copies than other birds, but little is currently known about expression and function of these gene copies. Non-passerine birds have a single highly expressed MHC-I gene copy, a pattern that seems unlikely in passerines. We used high-throughput sequencing to study MHC-I alleles in siskins (Spinus spinus) and determined gene expression, phylogenetic relationships and sequence divergence. We verified between six and 16 MHC-I alleles per individual and 97% of these were expressed. Strikingly, up to five alleles per individual had high expression. Out of 88 alleles 18 were putatively non-classical with low sequence divergence and expression, and found in a single phylogenetic cluster. The remaining 70 alleles were classical, with high sequence divergence and variable degrees of expression. Our results contradict the suggestion that birds only have a single dominantly expressed MHC-I gene by demonstrating several highly expressed MHC-I gene copies in a passerine.

Microbial Danger Signals Control Transcriptional Induction of Distinct MHC Class I L Lineage Genes in Atlantic Salmon

Antigen processing and presentation by major histocompatibility complex (MHC) molecules is a cornerstone in vertebrate immunity. Like mammals, teleosts possess both classical MHC class I and multiple families of divergent MHC class I genes. However, while certain mammalian MHC class I-like molecules have proven to be integral in immune regulation against a broad array of pathogens, the biological relevance of the different MHC class I lineages in fish remains elusive. This work focuses on MHC class I L lineage genes and reveals unique regulatory patterns of six genes (Sasa-lia, Sasa-lda, Sasa-lca, Sasa-lga, Sasa-lha, and Sasa-lfa) in antimicrobial immunity of Atlantic salmon (Salmo salar L.). Using two separate in vivo challenge models with different kinetics and immune pathologies combined with in vitro stimulation using viral and bacterial TLR ligands, we show that de novo synthesis of different L lineage genes is distinctly regulated in response to various microbial stimuli. Prior to the onset of classical MHC class I gene expression, lia was rapidly and systemically induced in vivo by the single-stranded (ss) RNA virus salmonid alpha virus 3 (SAV3) but not in response to the intracellular bacterium Piscirickettsia salmonis. In contrast, lga expression was upregulated in response to both viral and bacterial stimuli. A role for distinct MHC class I L-lineage genes in anti-microbial immunity in salmon was further substantiated by a marked upregulation of lia and lga gene expression in response to type I IFNa stimulation in vitro. Comparably, lha showed no transcriptional induction in response to IFNa stimulation but was strongly induced in response to a variety of viral and bacterial TLR ligands. In sharp contrast, lda showed no response to viral or bacterial challenge. Similarly, induction of lca, which is predominantly expressed in primary and secondary lymphoid tissues, was marginal with the exception of a strong and transient upregulation in pancreas following SAV3 challenge Together, these findings suggest that certain Atlantic salmon MHC class I L lineage genes play important and divergent roles in early anti-microbial response and that their regulation, in response to different activation signals, represents a system for selectively promoting the expression of distinct non-classical MHC class I genes in response to different types of immune challenges.

Complete loss of the MHC II pathway in an anglerfish, Lophius piscatorius

Genome studies in fish provide evidence for the adaptability of the vertebrate immune system, revealing alternative immune strategies. The reported absence of the major compatibility complex (MHC) class II pathway components in certain species of pipefish (genus Syngnathus) and cod-like fishes (order Gadiformes) is of particular interest. The MHC II pathway is responsible for immunization and defence against extracellular threats through the presentation of exogenous peptides to T helper cells. Here, we demonstrate the absence of all genes encoding MHC II components (CD4, CD74 A/B, and both classical and non-classical MHC II α/β) in the genome of an anglerfish, Lophius piscatorius, indicating loss of the MHC II pathway. By contrast, it has previously been reported that another anglerfish, Antennarius striatus, retains all MHC II genes, placing the loss of MHC II in the Lophius clade to their most recent common ancestor. In the three taxa where MHC II loss has occurred, the gene loss has been restricted to four or five core MHC II components, suggesting that, in teleosts, only these genes have functions that are restricted to the MHC II pathway.

Molecular cloning, characterization and expression analysis of MHCI and chemokines CXCR3 and CXCR4 gene from freshwater carp, Catla catla

The immune system with large number of molecules protects the host against a plethora of continuously evolving microbes. Major histocompatibility complex (MHC) molecules serve as cardinal elements of the adaptive immune system responsible for the activation of the adaptive immunity in the host. The present study reports MHCI molecule in freshwater carp, Catla catla, and its differential expression in immunologically relevant tissues post-infection with Gram-negative and Gram-positive bacteria. The MHCI sequence of C. catla had 502 bp nucleotides encoding putative 146 amino acids. The phylogenetic analysis exhibited its evolutionary conservation within the Cyprinidae family and formed a different clade with the higher vertebrates. Simultaneously, CXCR3 and CXCR4 chemokines were cloned and characterized for their expression in infected tissues. Analysis of immunologically relevant tissues of the infected fish exhibited an increase of MHCI gene expression and the down-regulation of CXCR3 and CXCR4 chemokines, indicating a tricky interaction between the innate and adaptive immune system. It was found that intestine, skin and spleen played a crucial role in the contribution of the defense activity which instigated the self-immunity. These immune activities can provide useful information to understand the interaction of self and non-self- immune system in freshwater fish, Catla catla.

The key differentially expressed genes and proteins related to immune response in the spleen of pufferfish (Takifugu obscurus) infected by Aeromonas hydrophila

The immune mechanism elicited in pufferfish (Takifugu obscurus) against the invasion of Aeromonas hydrophila is still poorly understood. We examined the spleen of pufferfish at the transcriptome and proteome levels by using Illumina-seq and TMT coupled mass spectrometry after 12 h infection by A. hydrophila, respectively. A total of 2,339 genes (1,512 up-regulated and 827 down-regulated) and 537 (237 up-regulated and 300 down-regulated) proteins were identified. GO and KEGG analyses revealed that the responses to stimulus were the main biological processes, intestinal immune network for IgT production and calcium signaling pathway. Fourteen genes (8 up-regulated and 6 down-regulated) and proteins (5 up-regulated and 9 down-regulated) involved immune responses or signal transduction were validated by qRT-PCR and parallel reaction monitoring to confirm the reliability of the transcriptomic and proteomic analyses, respectively. Moreover, qRT-PCR and flow cytometry were used to detect dynamics of the genes in calcium signaling pathway and changes of concentration of cytoplasm Ca²⁺ in spleen cells within a 72 h challenge. This study provides the findings regarding immune response, especially intestinal immune network for IgT production pathway and calcium signaling pathway at the molecular, protein and cellular in pufferfish after infection by A. hydrophila. These results would provide a new insight and molecular targets into the response to pathogenic infection in pufferfish.

Association of MHC IIA polymorphisms with disease resistance in Aeromonas hydrophila-challenged Nile tilapia

The major histocompatibility complex (MHC) genes show high polymorphisms in vertebrates depending on animal immunity status. Herein, MHC class IIA gene in Aeromonas hydrophila-challenged Nile tilapia was screened for presence of polymorphisms using sequencing. Twelve nucleotides deletion polymorphism was determined with a PCR product size of 267 bp in the resistant fish and 255 bp in the control and susceptible/diseased fish. Additionally, a non-synonymous right frameshift c.712 T > G (P. 238 * > G) SNP was detected at the stop codon (*). SNP-susceptibility association analysis revealed that fish carrying GG genotype and allele G were high susceptible (risk) for A. hydrophila, and had lower immune response as indicated by significant reduction in non-specific immune parameters (total protein, globulin, IgM, phagocytic activity, phagocytic index, and lysosome activity) and mRNA level of MHC IIA, interleukin 1 beta (IL1β), tumor necrosis factor alfa (TNFα), and toll-like receptor 7 (TLR7) in the spleen and head kidney. Thus, G allele could be considered as a risk (recessive or mutant) allele for c. 712 T > G (P. 238 * > G) SNP and so selection of Nile tilapia with protective allele (T) for this SNP could improve the disease resistant of the fish.

Dynamic distribution of formalin-inactivated Edwardsiella tarda in flounder (Paralichthys olivaceus) post intraperitoneal vaccination

In order to investigate the dynamic distribution of antigen in different tissues post vaccination, an absolute real-time quantitative PCR was employed to detect the amount of antigen in flounder (Paralichthys olivaceus) post intraperitoneal (i.p.) injection with three concentrations (10⁷, 10⁸, 10⁹ CFU ml^-1) of formalin-inactivated Edwardsiella tarda bacterin. The results showed that the amount of uptaken antigen quickly increased and then decreased in different tissues. The peak occurred first in the spleen and head kidney at 6-9 h after injection, and in the liver and blood at 9-15 h, then in the gill, intestine and skin at 15-24 h, finally in the muscle at 24-36 h. The amount of antigen was highest in the spleen and head kidney, followed by the blood, liver and gill, and lowest in the intestine, skin and muscle. Among the three concentration groups, the amount of antigen increased with the increasing concentration of the vaccine in the blood, liver, gill, intestine, skin and muscle, except for the spleen and head kidney, in which more antigens were found in the 10⁸ CFU ml^-1 group than that in 10⁹ CFU ml^-1 group. Moreover, IIFA and western blotting was performed to examine the tissue distribution of antigen at 9 h after vaccination with 10⁸ CFU ml^-1 formalin-inactivated E. tarda. The bacteria were mainly observed in the spleen and head kidney, then the liver, gill and blood, and least in the intestine, skin and muscle, which was roughly in accordance with the results of absolute qPCR. Furthermore, the expressions of CD4-1, MHC IIα, CD8α and MHC Iα in different tissues were detected by RT-qPCR, and the expression levels of these genes were highest in the spleen and head kidney, then in the blood, gill, liver, and lowest in the intestine, skin and muscle. All these results provided useful information for dynamic transportation of antigen uptake post vaccination, and also deepened the understanding of immune response to the injection vaccination.

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.

Disentangling the immune response and host-pathogen interactions in Francisella noatunensis infected Atlantic cod

The genetic repertoire underlying teleost immunity has been shown to be highly variable. A rare example is Atlantic cod and its relatives Gadiformes that lacks a hallmark of vertebrate immunity: Major Histocompatibility Complex class II. No immunological studies so far have fully unraveled the functionality of this particular immune system. Through global transcriptomic profiling, we investigate the immune response and host-pathogen interaction of Atlantic cod infected with the facultative intracellular bacterium Francisella noatunensis. We find that Atlantic cod displays an overall classic innate immune response with inflammation, acute-phase proteins and cell recruitment through up-regulation of e.g. IL1B, fibrinogen, cathelicidin, hepcidin and several chemotactic cytokines such as the neutrophil attractants CXCL1 and CXCL8. In terms of adaptive immunity, we observe up-regulation of interferon gamma followed by up-regulation of several MHCI transcripts and genes related to antigen transport and loading. Finally, we find up-regulation of immunoglobulins and down-regulation of T-cell and NK-like cell markers. Our analyses also uncover some contradictory transcriptional findings such as up-regulation of anti-inflammatory IL10 as well as down-regulation of the NADPH oxidase complex and myeloperoxidase. This we interpret as the result of host-pathogen interactions where F. noatunensis modulates the immune response. In summary, our results suggest that Atlantic cod mounts a classic innate immune response as well as a neutrophil-driven response. In terms of adaptive immunity, both endogenous and exogenous antigens are being presented on MHCI and antibody production is likely enabled through direct B-cell stimulation with possible neutrophil help. Collectively, we have obtained novel insight in the orchestration of the Atlantic cod immune system and determined likely targets of F. noatunensis host-pathogen interactions.

Molecular and Structural Characterization of MHC Class II β Genes Reveals High Diversity in the Cold-Adapted Icefish Chionodraco hamatus

This study reports the presence of two distinct MHC class II β genes in the Antarctic icefish Chionodraco hamatus, belonging to the classical (ChhaDAB) and nonclassical (ChhaDBB) evolutionary lineages. By the application of targeted sequencing approach, a remarkable molecular diversity in the exon 2 sequence of the highly expressed gene ChhaDAB has been observed, resulting in an estimate of 92 different variants translated in 87 different peptides from 54 analysed icefish individuals. A highly conservative estimate, based on a 95% sequence identity threshold clustering, translate this variability in 41 different peptide clusters belonging to four different clades and showing the signature of different kinds of selection. In stark contrast, the poorly expressed ChhaDBB gene displayed a very low level of molecular diversity within exon 2, in agreement with expectations for a nonclassical MHC class II β gene.

Generation of virus-specific CD8+ T cells by vaccination with inactivated virus in the intestine of ginbuna crucian carp

Although a previous study using ginbuna crucian carp suggested that cell-mediated immunity can be induced by the oral administration of inactivated viruses, which are exogenous antigens, there is no direct evidence that CD8⁺ cytotoxic T cells (CTLs) in teleost fish are generated by vaccination with exogenous antigens. In the present study, we investigated whether antigen-specific CD8⁺ CTLs in ginbuna crucian carp can be elicited by intestinal immunization with an exogenous antigen without any adjuvant. The IFNγ-1 and T-bet mRNA expressions were up-regulated in intestinal leukocytes following the administration of formalin-inactivated crucian hematopoietic necrosis virus (FI-CHNV), whereas the down-regulation of these genes was observed in kidney leukocytes. Furthermore, an increase in the percentage of proliferating CD8⁺ cells was detected in the posterior portion of the hindgut, suggesting that the virus-specific CTLs are locally generated in this site. In addition, cell-mediated cytotoxicity against CHNV-infected syngeneic cells and the in vivo inhibition of viral replication were induced by immunization with FI-CHNV. Unexpectedly, intraperitoneal immunization with FI-CHNV induced a type I helper T cell (Th1)-response in the intestine, but not in the kidney; however, its effect was slightly lower than that reported after intestinal immunization. These findings suggest that the posterior portion of the intestine is an important site for generating virus-specific CTLs by vaccination with the inactivated vaccine.