MHC polymorphism and disease resistance in Atlantic salmon (Salmo salar); facing pathogens with single expressed major histocompatibility class I and class II loci

Understanding the evolution of immune genes in jawed vertebrates

Driven by co-evolution with pathogens, host immunity continuously adapts to optimize defence against pathogens within a given environment. Recent advances in genetics, genomics and transcriptomics have enabled a more detailed investigation into how immunogenetic variation shapes the diversity of immune responses seen across domestic and wild animal species. However, a deeper understanding of the diverse molecular mechanisms that shape immunity within and among species is still needed to gain insight into-and generate evolutionary hypotheses on-the ultimate drivers of immunological differences. Here, we discuss current advances in our understanding of molecular evolution underpinning jawed vertebrate immunity. First, we introduce the immunome concept, a framework for characterizing genes involved in immune defence from a comparative perspective, then we outline how immune genes of interest can be identified. Second, we focus on how different selection modes are observed acting across groups of immune genes and propose hypotheses to explain these differences. We then provide an overview of the approaches used so far to study the evolutionary heterogeneity of immune genes on macro and microevolutionary scales. Finally, we discuss some of the current evidence as to how specific pathogens affect the evolution of different groups of immune genes. This review results from the collective discussion on the current key challenges in evolutionary immunology conducted at the ESEB 2021 Online Satellite Symposium: Molecular evolution of the vertebrate immune system, from the lab to natural populations.

Functional immune diversity in reindeer reveals a high Arctic population at risk

Climate changes the geographic range of both species as well as pathogens, causing a potential increase in the vulnerability of populations or species with limited genetic diversity. With advances in high throughput sequencing (HTS) technologies, we can now define functional expressed genetic diversity of wild species at a larger scale and identify populations at risk. Previous studies have used genomic DNA to define major histocompatibility complex (MHC) class II diversity in reindeer. Varying numbers of expressed genes found in many ungulates strongly argues for using cDNA in MHC typing strategies to ensure that diversity estimates relate to functional genes. We have used available reindeer genomes to identify candidate genes and established an HTS approach to define expressed MHC class I and class II diversity. To capture a broad diversity we included samples from wild reindeer from Southern Norway, semi-domesticated reindeer from Northern Norway and reindeer from the high Artic archipelago Svalbard. Our data show a medium MHC diversity in semi-domesticated and wild Norwegian mainland reindeer, and low MHC diversity reindeer in Svalbard reindeer. The low immune diversity in Svalbard reindeer provides a potential risk if the pathogenic pressure changes in response to altered environmental conditions due to climate change, or increased human-related activity.

In silico prediction of B and T cell epitopes of infectious salmon anemia virus proteins and molecular modeling of T cell epitopes to salmon major histocompatibility complex (MHC) class I

Infectious salmon anemia (ISA) can be devastating in farmed Atlantic salmon (Salmo salar). The disease can evolve into epidemics if it is not contained and controlled. ISA epidemics were seen in Norway in the early 1990s and Chile in 2007-2009. Consequently, there is an urgent need to develop a vaccine to prevent or treat the infection. In this study, an immunoinformatic approach was employed to predict 32 lineal B-cell epitopes based on antigenicity and surface accessibility prediction for ISAV fusion (F), hemagglutinin-esterase (HE), and matrix (M) proteins. On the other hand, twelve conformational B-cell epitopes were also predicted. We further identified six antigenic cytotoxic T lymphocyte (CTL) epitopes and investigated the binding interactions with five salmon MHC-I proteins after docking the peptides to the binding groove of the MHC-I proteins. Our results showed that all the predicted epitopes could bind to salmon MHC-I with high negative ΔG values with medium to high binding affinities. Hence, the predicted epitopes have a high potential of being recognized by Atlantic salmon MHC-I to elicit a CD8⁺ T cell response in salmon. The predicted and analyzed B and T cell antigenic epitopes in this work might present an initial set of peptides for future vaccine development against ISAV. The ability to model and predict these interactions will ultimately lead to the ability to predict potential binding for MHCs and epitopes that were not studied previously. As current knowledge of salmon MHC specificity is limited, studying and modeling interactions in the peptide/MHC complex is a key to resolving unknown epitope specificity.

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.

MHC class I-α population differentiation in a commercial fish, the European sea bass (Dicentrarchus labrax)

Identifying population structuring in highly fecund marine species with high dispersal rates is challenging, but critical for conservation and stock delimitation for fisheries management. European sea bass (Dicentrarchus labrax) is a commercial species of fisheries and aquaculture relevance whose stocks are declining in the North Atlantic, despite management measures to protect them and identifying their fine population structure is needed for managing their exploitation. As for other marine fishes, neutral genetic markers indicate that eastern Atlantic sea bass form a panmictic population and is currently managed as arbitrarily divided stocks. The genes of the major histocompatibility complex (MHC) are key components of the adaptive immune system and ideal candidates to assess fine structuring arising from local selective pressures. We used Illumina sequencing to characterise allelic composition and signatures of selection at the MHC class I-α region of six D. labrax populations across the Atlantic range. We found high allelic diversity driven by positive selection, corresponding to moderate supertype diversity, with 131 alleles clustering into four to eight supertypes, depending on the Bayesian information criterion threshold applied, and a mean number of 13 alleles per individual. Alleles could not be assigned to particular loci, but private alleles allowed us to detect regional genetic structuring not found previously using neutral markers. Our results suggest that MHC markers can be used to detect cryptic population structuring in marine species where neutral markers fail to identify differentiation. This is particularly critical for fisheries management, and of potential use for selective breeding or identifying escapes from sea farms.

Functional analysis of a novel MHC-Iα genotype in orange-spotted grouper: Effects on Singapore grouper iridovirus (SGIV) replication and apoptosis

The classical major histocompatibility complex class I (MHC-Ⅰ) molecule plays a key role in vertebrate immune response for its important functions in antigen presentation and immune regulation. MHC pathway is closely related to many diseases involving autoimmunity, antigen intrusion and inflammation. However, rare literatures about the effect of MHC-I on fish cells apoptosis were reported. In this study, a novel type of MHC-Ⅰα genotype from orange-spotted grouper (named EcMHC-ⅠA*01) were cloned and characterized. It shared a 77% identity to its Epinephelus coioides MHC-Iα homology that has been uploaded to NCBI (ACZ97571.1). Molecular characterization analysis showed that EcMHC-ⅠA*01 encodes a 357-amino-acid protein, containing a signal peptide，α1，α2，α3, Cytoplasmic (Cyt) and Transmembrane (TM) domains. Tissue expression pattern showed that EcMHC-ⅠA*01 was extensively distributed in twelve selected tissues, with higher expression in the gill, intestine and skin. The expression of EcMHC-ⅠA*01 in grouper liver and spleen tissues were significantly induced by different stimuli (Zymosan A, LPS, Ploy I:C, RGNNV and SGIV). Comparing with the EcMHC-ⅠA*01 expression levels induced by Zymosan A, Ploy I:C and RGNNV, the effects induced by SGIV and LPS were more significant. Subcellular localization analysis showed that EcMHC-ⅠA*01 localizes throughout the cytoplasm appeared both diffuse and focal intracellular expression pattern. Overexpression of EcMHC-ⅠA*01 inhibited the CPE progression, the mRNA expression of the SGIV related genes (MCP, LITAF, ICP-18 and VP19) and the protein expression of MCP. Meanwhile, qRT-PCR result showed that EcMHC-ⅠA*01 overexpression upregulated the expression of interferon signaling molecules (IFN-γ, ISG56, MDA5 and MXI) and inflammatory cytokines (IL-1β, IL-6, TNF-α and TRAF6). In addition, our results showed that overexpression of EcMHC-ⅠA*01 promoted the apoptosis of normal fathead minnow (FHM) cells as well as the apoptosis of FHM cells induced by SGIV. However, there was no significant change in the activity of caspase 3 between control group and EcMHC-ⅠA*01 overexpression group, suggesting that EcMHC-ⅠA*01-induced apoptosis may not depend on the caspase 3 pathway. Taken together, these data in our study provide new insights into the role of MHC-I in antiviral immune response and apoptosis in fish.

Comparison of Gene Expression Between Resistant and Susceptible Families Against VPAHPND and Identification of Biomarkers Used for Resistance Evaluation in Litopenaeus vannamei

Acute hepatopancreatic necrosis disease (AHPND) has caused a heavy loss to shrimp aquaculture since its outbreak. Vibrio parahaemolyticus (VP_AHPND) is regarded as one of the main pathogens that caused AHPND in the Pacific white shrimp Litopenaeus vannamei. In order to learn more about the mechanism of resistance to AHPND, the resistant and susceptible shrimp families were obtained through genetic breeding, and comparative transcriptome approach was used to analyze the gene expression patterns between resistant and susceptible families. A total of 95 families were subjected to VP_AHPND challenge test, and significant variations in the resistance of these families were observed. Three pairs of resistant and susceptible families were selected for transcriptome sequencing. A total of 489 differentially expressed genes (DEGs) that presented in at least two pairwise comparisons were screened, including 196 DEGs highly expressed in the susceptible families and 293 DEGs in the resistant families. Among these DEGs, 16 genes demonstrated significant difference in all three pairwise comparisons. Gene set enrichment analysis (GSEA) of all 27,331 expressed genes indicated that some energy metabolism processes were enriched in the resistant families, while signal transduction and immune system were enriched in the susceptible families. A total of 32 DEGs were further confirmed in the offspring of the detected families, among which 19 genes were successfully verified. The identified genes in this study will be useful for clarifying the genetic mechanism of shrimp resistance against Vibrio and will further provide molecular markers for evaluating the disease resistance of shrimp in the breeding program.

Whole-genome association study searching for QTL for Aeromonas salmonicida resistance in rainbow trout

Aeromonas salmonicida subsp. salmonicida, the causative agent of furunculosis, has extensive negative effects on wild and farmed salmonids worldwide. Vaccination induces some protection under certain conditions but disease outbreaks occur even in vaccinated fish. Therefore, alternative disease control approaches are required to ensure the sustainable expansion of rainbow trout aquaculture. Selective breeding can be applied to enhance host resistance to pathogens. The present work used genome-wide association study (GWAS) to identify quantitative trait loci (QTL) associated with A. salmonicida resistance in rainbow trout. A total 798 rainbow trout exposed to A. salmonicida by bath challenge revealed 614 susceptible and 138 resistant fish. Genotyping was conducted using the 57 K single nucleotide polymorphism (SNP) array and the GWAS was performed for survival and time to death phenotypes. We identified a QTL on chromosome 16 and located positional candidate genes in the proximity of the most significant SNPs. In addition, samples from exposed fish were examined for expression of 24 immune-relevant genes indicating a systematic immune response to the infection. The present work demonstrated that resistance to A. salmonicida is moderately heritable with oligogenic architecture. These result will be useful for the future breeding programs for improving the natural resistance of rainbow trout against furunculosis.

Challenges and Solutions to Viral Diseases of Finfish in Marine Aquaculture

Aquaculture is the fastest food-producing sector in the world, accounting for one-third of global food production. As is the case with all intensive farming systems, increase in infectious diseases has adversely impacted the growth of marine fish farming worldwide. Viral diseases cause high economic losses in marine aquaculture. We provide an overview of the major challenges limiting the control and prevention of viral diseases in marine fish farming, as well as highlight potential solutions. The major challenges include increase in the number of emerging viral diseases, wild reservoirs, migratory species, anthropogenic activities, limitations in diagnostic tools and expertise, transportation of virus contaminated ballast water, and international trade. The proposed solutions to these problems include developing biosecurity policies at global and national levels, implementation of biosecurity measures, vaccine development, use of antiviral drugs and probiotics to combat viral infections, selective breeding of disease-resistant fish, use of improved diagnostic tools, disease surveillance, as well as promoting the use of good husbandry and management practices. A multifaceted approach combining several control strategies would provide more effective long-lasting solutions to reduction in viral infections in marine aquaculture than using a single disease control approach like vaccination alone.

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.

Transcriptome analysis reveals immune pathways underlying resistance in the common carp Cyprinus carpio against the oomycete Aphanomyces invadans

Infection with Aphanomyces invadans is a serious fish disease with major global impacts. Despite affecting over 160 fish species, some of the species like the common carp Cyprinus carpio are resistant to A. invadans infection. In the present study, we investigated the transcriptomes of head kidney of common carp experimentally infected with A. invadans. In time course analysis, 5288 genes were found to be differentially expressed (DEGs), of which 731 were involved in 21 immune pathways. The analysis of immune-related DEGs suggested that efficient processing and presentation of A. invadans antigens, enhanced phagocytosis, recognition of pathogen-associated molecular patterns, and increased recruitment of leukocytes to the sites of infection contribute to resistance of common carp against A. invadans. Herein, we provide a systematic understanding of the disease resistance mechanisms in common carp at molecular level as a valuable resource for developing disease management strategies for this devastating fish-pathogenic oomycete.

Salmonid Antibacterial Immunity: An Aquaculture Perspective

The aquaculture industry is continuously threatened by infectious diseases, including those of bacterial origin. Regardless of the disease burden, aquaculture is already the main method for producing fish protein, having displaced capture fisheries. One attractive sector within this industry is the culture of salmonids, which are (a) uniquely under pressure due to overfishing and (b) the most valuable finfish per unit of weight. There are still knowledge gaps in the understanding of fish immunity, leading to vaccines that are not as effective as in terrestrial species, thus a common method to combat bacterial disease outbreaks is the use of antibiotics. Though effective, this method increases both the prevalence and risk of generating antibiotic-resistant bacteria. To facilitate vaccine design and/or alternative treatment efforts, a deeper understanding of the teleost immune system is essential. This review highlights the current state of teleost antibacterial immunity in the context of salmonid aquaculture. Additionally, the success of current techniques/methods used to combat bacterial diseases in salmonid aquaculture will be addressed. Filling the immunology knowledge gaps highlighted here will assist in reducing aquaculture losses in the future.

Transcriptome sequencing and histology reveal dosage compensation in the liver of triploid pre-smolt Atlantic salmon

Triploid Atlantic salmon (Salmo salar L.) is seen as one of the best solutions to solve key issues in the salmon farming industry, such as the impact of escapees on wild stocks and pre-harvest sexual maturation. However, the effects of triploidy on salmon smoltification are poorly understood at the molecular level, even though smoltification is a very sensitive period that has a major influence on survival rate and performance of farmed salmon. In this study, we have compared the liver transcriptomes of diploid and triploid Atlantic salmon at three ontogeny stages: fry, parr and smolt. In diploid fish, a total of 2,655 genes were differentially expressed between fry and parr, whereas 506 genes had significantly different transcript levels between parr and smolts. In triploids, 1,507 and 974 genes were differentially expressed between fry and parr, and between parr and smolts, respectively. Most of these genes were down-regulated and 34 genes were differentially expressed between ploidies at the same stage. In both ploidy groups, the top differentially expressed genes with ontogeny stage belonged to common functional categories that can be related to smoltification. Nucleotide and energy metabolism were significantly down-regulated in fry when compared to parr, while immune system processes were significantly down-regulated in parr when compared to smolts. The close resemblance of enriched biological processes and pathways between ploidy groups suggests that triploidy is regulated by genome dosage compensation in Atlantic salmon. Histological analysis revealed that areas of vacuolization (steatosis) were present only in fry and parr stages, in contrast to a compact cellular histology with glycogen granules after smoltification. There was no significant difference in vacuolization between ploidy groups at the fry stage but the liver of diploid parr had a 33.5% higher vacuolization area compared to their triploid counterparts. Taken together, our data provide novel insights into the changes that occur at the molecular and histological level in the liver of both diploid and triploid Atlantic salmon prior to and during smoltification.

QTLs Associated with Resistance to Cardiomyopathy Syndrome in Atlantic Salmon

Cardiomyopathy syndrome (CMS) caused by piscine myocarditis virus is a major disease affecting the Norwegian Atlantic salmon industry. Three different populations of Atlantic salmon from the Mowi breeding program were used in this study. The first 2 populations (population 1 and 2) were naturally infected in a field outbreak, while the third population (population 3) went through a controlled challenged test. The aim of the study was to estimate the heritability, the genetic correlation between populations and perform genome-wide association analysis for resistance to this disease. Survival data from population 1 and 2 and heart atrium histology score data from population 3 was analyzed. A total of 571, 4312, and 901 fish from population 1, 2, and 3, respectively were genotyped with a noncommercial 55,735 Affymetrix marker panel. Genomic heritability ranged from 0.12 to 0.46 and the highest estimate was obtained from the challenge test dataset. The genetic correlation between populations was moderate (0.51–0.61). Two chromosomal regions (SSA27 and SSA12) contained single nucleotide polymorphisms associated with resistance to CMS. The highest association signal (P = 6.9751 × 10⁻²⁷) was found on chromosome 27. Four genes with functional roles affecting viral resistance (magi1, pi4kb, bnip2, and ha1f) were found to map closely to the identified quantitative trait loci (QTLs). In conclusion, genetic variation for resistance to CMS was observed in all 3 populations. Two important quantitative trait loci were detected which together explain half of the total genetic variance, suggesting strong potential application for marker-assisted selection and genomic predictions to improve CMS resistance.

Combining Multiple Approaches and Models to Dissect the Genetic Architecture of Resistance to Infections in Fish

Infectious diseases represent a major threat for the sustainable development of fish farming. Efficient vaccines are not available against all diseases, and growing antibiotics resistance limits the use of antimicrobial drugs in aquaculture. It is therefore important to understand the basis of fish natural resistance to infections to help genetic selection and to develop new approaches against infectious diseases. However, the identification of the main mechanisms determining the resistance or susceptibility of a host to a pathogenic microbe is challenging, integrating the complexity of the variation of host genetics, the variability of pathogens, and their capacity of fast evolution and adaptation. Multiple approaches have been used for this purpose: (i) genetic approaches, QTL (quantitative trait loci) mapping or GWAS (genome-wide association study) analysis, to dissect the genetic architecture of disease resistance, and (ii) transcriptomics and functional assays to link the genetic constitution of a fish to the molecular mechanisms involved in its interactions with pathogens. To date, many studies in a wide range of fish species have investigated the genetic determinism of resistance to many diseases using QTL mapping or GWAS analyses. A few of these studies pointed mainly toward adaptive mechanisms of resistance/susceptibility to infections; others pointed toward innate or intrinsic mechanisms. However, in the majority of studies, underlying mechanisms remain unknown. By comparing gene expression profiles between resistant and susceptible genetic backgrounds, transcriptomics studies have contributed to build a framework of gene pathways determining fish responsiveness to a number of pathogens. Adding functional assays to expression and genetic approaches has led to a better understanding of resistance mechanisms in some cases. The development of knock-out approaches will complement these analyses and help to validate putative candidate genes critical for resistance to infections. In this review, we highlight fish isogenic lines as a unique biological material to unravel the complexity of host response to different pathogens. In the future, combining multiple approaches will lead to a better understanding of the dynamics of interaction between the pathogen and the host immune response, and contribute to the identification of potential targets of selection for improved resistance.

Genome-Wide Association Study Confirms Previous Findings of Major Loci Affecting Resistance to Piscine myocarditis virus in Atlantic Salmon (Salmo salar L.)

Cardiomyopathy syndrome is a viral disease of Atlantic salmon, mostly affecting fish during the late stages of production, resulting in significant losses to the industry. It has been shown that resistance to this disease has a strong genetic component, with quantitative trait loci (QTL) on chromosomes 27 (Ssa27) and Ssa12 to explain most of the additive genetic variance. Here, by analysing animals from a different year-class and a different population, we further aimed to confirm and narrow down the locations of these QTL. The data support the existence of the two QTL and suggest that the causative mutation on Ssa27 is most likely within the 10–10.5 Mbp segment of this chromosome. This region contains a cluster of major histocompatibility complex class I (MHC I) genes with the most strongly associated marker mapped to one of these loci. On Ssa12, the data confirmed the previous finding that the location of the causative mutation is within the 61.3 to 61.7 Mbp region. This segment contains several immune-related genes, but of particular interest are genes related to MHC II. Together, these findings highlight the likely key role of MHC genes in Atlantic salmon following infection with Piscine myocarditis virus (PMCV) and their potential impact on influencing the trajectory of this disease.

Linking Virus Discovery to Immune Responses Visualized during Zebrafish Infections

The discovery of new viruses currently outpaces our capacity for experimental examination of infection biology. To better couple virus discovery with immunology, we genetically modified zebrafish to visually report on virus infections. After generating a strain that expresses green fluorescent protein (GFP) under an interferon-stimulated gene promoter, we repeatedly observed transgenic larvae spontaneously expressing GFP days after hatching. RNA sequencing comparisons of co-housed GFP-positive and GFP-negative zebrafish revealed a naturally occurring picornavirus that induced a canonical interferon-mediated response and hundreds of antiviral defense genes not observed following immunostimulatory treatments or experimental infections with other viruses. Among the many genes induced by picornavirus infection was a large set encoding guanosine triphosphatase (GTPase) of immunity-associated proteins (GIMAPs). The GIMAP gene family is massively expanded in fish genomes and may also play a crucial role in antiviral responses in mammals, including humans. We subsequently detected zebrafish picornavirus in publicly available sequencing data from seemingly asymptomatic zebrafish in many research institutes and found that it altered gene expression in a previous study of zebrafish development. Experiments revealed a horizontal mode of virus transmission, highlighting a system for studying the spread of picornavirus infections within and between individuals. Our study describes a naturally occurring picornavirus that elicits strong antiviral responses in zebrafish and provides new strategies for simultaneously discovering viruses and their impact on vertebrate hosts.

Single-Nucleotide Polymorphisms (SNP) Mining and Their Effect on the Tridimensional Protein Structure Prediction in a Set of Immunity-Related Expressed Sequence Tags (EST) in Atlantic Salmon (Salmo salar)

Single-nucleotide polymorphisms (SNPs) are single genetic code variations considered one of the most common forms of nucleotide modifications. Such SNPs can be located in genes associated to immune response and, therefore, they may have direct implications over the phenotype of susceptibility to infections affecting the productive sector. In this study, a set of immune-related genes (cc motif chemokine 19 precursor [ccl19], integrin β2 (itβ2, also named cd18), glutathione transferase omega-1 [gsto-1], heat shock 70 KDa protein [hsp70], major histocompatibility complex class I [mhc-I]) were analyzed to identify SNPs by data mining. These genes were chosen based on their previously reported expression on infectious pancreatic necrosis virus (IPNV)-infected Atlantic salmon phenotype. The available EST sequences for these genes were obtained from the Unigene database. Twenty-eight SNPs were found in the genes evaluated and identified most of them as transition base changes. The effect of the SNPs located on the 5'-untranslated region (UTR) or 3'-UTR upon transcription factor binding sites and alternative splicing regulatory motifs was assessed and ranked with a low-medium predicted FASTSNP score risk. Synonymous SNPs were found on itβ2 (c.2275G > A), gsto-1 (c.558G > A), and hsp70 (c.1950C > T) with low FASTSNP predicted score risk. The difference in the relative synonymous codon usage (RSCU) value between the variant codons and the wild-type codon (ΔRSCU) showed one negative (hsp70 c.1950C > T) and two positive ΔRSCU values (itβ2 c.2275G > A; gsto-1 c.558G > A), suggesting that these synonymous SNPs (sSNPs) may be associated to differences in the local rate of elongation. Nonsynonymous SNPs (nsSNPs) in the gsto-1 translatable gene region were ranked, using SIFT and POLYPHEN web-tools, with the second highest (c.205A > G; c484T > C) and the highest (c.499T > C; c.769A > C) predicted score risk possible. Using homology modeling to predict the effect of these nonsynonymous SNPs, the most relevant nucleotide changes for gsto-1 were observed for the nsSNPs c.205A > G, c484T > C, and c.769A > C. Molecular dynamics was assessed to analyze if these GSTO-1 variants have significant differences in their conformational dynamics, suggesting these SNPs could have allosteric effects modulating its catalysis. Altogether, these results suggest that candidate SNPs identified may play a crucial potential role in the immune response of Atlantic salmon.

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.

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.

Early Developmental Stress Affects Subsequent Gene Expression Response to an Acute Stress in Atlantic Salmon: An Approach for Creating Robust Fish for Aquaculture?

Stress during early life has potential to program and alter the response to stressful events and metabolism in later life. Repeated short exposure of Atlantic salmon to cold water and air during embryonic (E), post-hatch (PH) or both phases of development (EPH) has been shown to alter the methylome and transcriptome and to affect growth performance during later life compared to untreated controls (CO). The aim of this study was to investigate how the transcriptome of these fish responds to subsequent acute stress at the start feeding stage, and to describe methylation differences that might steer these changes. EPH treated fish showed the strongest down-regulation of corticotropin releasing factor 1, up-regulation of glucocorticoid receptor and 3-oxo-5-alpha-steroid 4-dehydrogenase 2 gene expression and a suppressed cortisol response 3 hr after the acute stress, differences that could influence hormesis and be affecting how EPH fish cope and recover from the stress event. Growth hormone 2 and insulin-like growth factor 1 were more strongly down-regulated following acute stress in EPH treated fish relative to E, PH and CO fish. This indicates switching away from growth toward coping with stress following stressful events in EPH fish. Genes implicated in immune function such as major histocompatibility class 1A, T-cell receptor and toll-like receptor also responded to acute stress differently in EPH treated fish, indicating that repeated stresses during early life may affect robustness. Differential DNA methylation was detected in regions mapping <500 bases from genes differentially responding to acute stress suggesting the involvement of epigenetic mechanisms. Stress treatments applied during early development therefore have potential as a husbandry tool for boosting the productivity of aquaculture by affecting how fish respond to stresses at critical stages of production.

MHC Class IIB gene polymorphisms associated with resistance/susceptibility to Streptococcus agalactiae in Nile tilapia Oreochromis niloticus

Genetic variation in the major histocompatibility complex (MHC) Class IIB was tested in Nile tilapia Oreochromis niloticus, and the association between the MHC IIB alleles and disease resistance was also studied. F3 fry offspring (n = 1200) from 12 full-sib families were challenged with Streptococcus agalactiae, which caused significantly different mortalities in different Nile tilapia families (11.00-81.10%). Twenty fry (F1) from each of the 12 families were selected to study the polymorphisms of the MHC Class IIB gene using PCR followed by cloning and sequencing methods. The results showed that the size of the amplified fragment was 770-797 bp. Thirty-seven sequences from 240 individuals revealed 22 different alleles, which belonged to 9 major allele types. Up to 63.58% of nucleotide positions were variable, while the proportion of the amino acid variable positions was up to 68.73%. According to the survival rate of offspring (F3) from 12 full-sib families, we deduced that the alleles Orni-DAB*0107, Orni-DAB*0201 and Orni-DAB*0302 were highly associated with resistance to S. agalactiae, while the allele Orni-DAB*0701 was associated with susceptibility to S. agalactiae. In addition, our previous study found that the allele Orni-DAB*0201 was more frequently distributed in the disease-resistant groups. Therefore, the allele Orni-DAB*0201 could be used as an S. agalactiae resistance-related MHC marker in molecular marker-assisted selective breeding programs for S. agalactiae-resistant Nile tilapia.

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.

Reproductive Strategy Inferred from Major Histocompatibility Complex-Based Inter-Individual, Sperm-Egg, and Mother-Fetus Recognitions in Giant Pandas (Ailuropoda melanoleuca)

Few major histocompatibility complex (MHC)-based mate choice studies include all MHC genes at the inter-individual, sperm-egg, and mother-fetus recognition levels. We tested three hypotheses of female mate choice in a 17-year study of the giant panda (Ailuropoda melanoleuca) while using ten functional MHC loci (four MHC class I loci: Aime-C, Aime-F, Aime-I, and Aime-L; six MHC class II loci: Aime-DRA, Aime-DRB3, Aime-DQA1, Aime-DQA2, Aime-DQB1, and Aime-DQB2); five super haplotypes (SuHa, SuHaI, SuHaII, DQ, and DR); and, seven microsatellites. We found female choice for heterozygosity at Aime-C, Aime-I, and DQ and for disassortative mate choice at Aime-C, DQ, and DR at the inter-individual recognition level. High mating success occurred in MHC-dissimilar mating pairs. No significant results were found based on any microsatellite parameters, suggesting that MHCs were the mate choice target and there were no signs of inbreeding avoidance. Our results indicate Aime-DQA1- and Aime-DQA2-associated disassortative selection at the sperm-egg recognition level and a possible Aime-C- and Aime-I-associated assortative maternal immune tolerance mechanism. The MHC genes were of differential importance at the different recognition levels, so all of the functional MHC genes should be included when studying MHC-dependent reproductive mechanisms.

MHC class IIA polymorphisms and their association with resistance-susceptibility to Streptococcus agalactiae in Nile tilapia, Oreochromis niloticus

The association between major histocompatibility complex (MHC) class IIA polymorphisms and the severity of infection by Streptococcus agalactiae was investigated using 40 susceptible and 40 resistant individuals of Nile tilapia Oreochromis niloticus. Twenty-five alleles were identified from 80 individuals, which belong to 22 major allele types. High polymorphism of mhcIIa gene and at least two loci were discovered in O. niloticus. In peptide-binding region (PBR) and non-PBR, the ratio of nonsynonymous substitution (dN) to synonymous substitution (dS) was 1.294 (>1) and 1.240 (>1), suggesting that the loci are evolving under positive balancing selection. Association analysis showed that the allele, orni-daa*0501, was significantly associated with resistance to S. agalactiae, while the alleles, orni-daa*1101, orni-daa*1301, orni-daa*1401 and orni-daa*1201, were associated with susceptibility to S. agalactiae. To confirm these correlations, another independent challenge experiment was performed in the Huizhou population of the O. niloticus. The frequency distribution showed that the orni-daa*1101 allele was significantly more frequent in the Huizhou-Susceptible group (HZ-SG) than in the Huizhou-Resistant group (HZ-RG) (P < 0.05), which was consistent with the first challenge. However, orni-daa*0501 did not present in HZ-SG and HZ-RG and the distribution frequencies of the orni-daa*1201, orni-daa*1301 and orni-daa*1401 alleles were not significantly more frequent in HZ-SG than in HZ-RG. These results indicate that the orni-daa*1101 allele confers susceptibility to S. agalactia infection. These results suggest that the diversity of exon 2 of mcaIIa alleles could be used to explore the association between disease susceptibility or resistance and the multiformity of mcaIIa and to achieve the molecular-assisted selection of O. niloticus with enhanced disease resistance.

Parallel evolution of gene classes, but not genes: Evidence from Hawai'ian honeycreeper populations exposed to avian malaria

Adaptation in nature is ubiquitous, yet characterizing its genomic basis is difficult because population demographics cause correlations with nonadaptive loci. Introduction events provide opportunities to observe adaptation over known spatial and temporal scales, facilitating the identification of genes involved in adaptation. The pathogen causing avian malaria, Plasmodium relictum, was introduced to Hawai'i in the 1930s and elicited extinctions and precipitous population declines in native honeycreepers. After a sharp initial population decline, the Hawai'i 'amakihi (Chlorodrepanis virens) has evolved tolerance to the parasite at low elevations where P. relictum exists, and can sustain infection without major fitness consequences. High-elevation, unexposed populations of 'amakihi display little to no tolerance. To explore the genomic basis of adaptation to P. relictum in low-elevation 'amakihi, we genotyped 125 'amakihi from the island of Hawai'i via hybridization capture to 40,000 oligonucleotide baits containing SNPs and used the reference 'amakihi genome to identify genes potentially under selection from malaria. We tested for outlier loci between low- and high-elevation population pairs and identified loci with signatures of selection within low-elevation populations. In some cases, genes commonly involved in the immune response (e.g., major histocompatibility complex) were associated with malaria presence in the population. We also detected several novel candidate loci that may be implicated in surviving malaria infection (e.g., beta-defensin, glycoproteins and interleukin-related genes). Our results suggest that rapid adaptation to pathogens may occur through changes in different immune genes, but in the same classes of genes, across populations.

MHC class IIα polymorphism and its association with resistance/susceptibility to Vibrio harveyi in golden pompano (Trachinotus ovatus)

The major histocompatibility complex (MHC) plays an important role in the vertebrate immune response to antigenic peptides, and it is essential for recognizing foreign pathogens in organisms. In this study, MHC class IIα (Trov-MHC IIα) from the golden pompano (Trachinotus ovatus) was first cloned and identified. The gene structure of Trov-MHC IIα was contained four exons and three introns. High levels of polymorphism were found in the exon 2 of Trov-MHC IIα. A total of 29 different MHC class IIα alleles with high polymorphism were identified from 80 individuals. The ratio of non-synonymous substitutions (dN) to synonymous substitutions (dS) was 3.157 (>1) in the peptide binding regions (PBRs) of Trov-MHC IIα, suggesting positive balancing selection. Six alleles were selected to analyze the association between alleles and resistance/susceptibility to Vibrio harveyi in golden pompano. The results showed that Trov-DAA*6401 and Trov-DAA*6702 alleles were associated with the resistance to V. harveyi in golden pompano, while alleles Trov-DAA*6304 and Trov-DAA*7301 were associated with the susceptibility to V. harveyi in golden pompano. This study confirmed the association between alleles of MHC class IIα and disease resistance, and also detected some alleles which might be correlated with high V. harveyi-resistance. These disease resistance-related MHC alleles could be used as potential genetic markers for molecular marker-assisted selective breeding in the golden pompano.

Disentangling the mechanisms of mate choice in a captive koala population

Successful captive breeding programs are crucial to the long-term survival of many threatened species. However, pair incompatibility (breeding failure) limits sustainability of many captive populations. Understanding whether the drivers of this incompatibility are behavioral, genetic, or a combination of both, is crucial to improving breeding programs. We used 28 years of pairing data from the San Diego Zoo koala colony, plus genetic analyses using both major histocompatibility complex (MHC)-linked and non-MHC-linked microsatellite markers, to show that both genetic and non-genetic factors can influence mating success. Male age was reconfirmed to be a contributing factor to the likelihood of a koala pair copulating. This trend could also be related to a pair's age difference, which was highly correlated with male age in our dataset. Familiarity was reconfirmed to increase the probability of a successful copulation. Our data provided evidence that females select mates based on MHC and genome-wide similarity. Male heterozygosity at MHC class II loci was associated with both pre- and post-copulatory female choice. Genome-wide similarity, and similarity at the MHC class II DAB locus, were also associated with female choice at the post-copulatory level. Finally, certain MHC-linked alleles were associated with either increased or decreased mating success. We predict that utilizing a variety of behavioral and MHC-dependent mate choice mechanisms improves female fitness through increased reproductive success. This study highlights the complexity of mate choice mechanisms in a species, and the importance of ascertaining mate choice mechanisms to improve the success of captive breeding programs.

Interpopulation Variation in the Atlantic Salmon Microbiome Reflects Environmental and Genetic Diversity

The microbiome has a crucial influence on host phenotype and is of broad interest to ecological and evolutionary research. Yet, the extent of variation that occurs in the microbiome within and between populations is unclear. We characterized the skin and gut microbiomes of seven populations of juvenile Atlantic salmon (Salmo salar) inhabiting a diverse range of environments, including hatchery-reared and wild populations. We found shared skin operational taxonomic units (OTUs) across all populations and core gut microbiota for all wild fish, but the diversity and structure of both skin and gut microbial communities were distinct between populations. There was a marked difference between the gut microbiomes of wild and captive fish. Hatchery-reared fish had lower intestinal microbial diversity, lacked core microbiota found in wild fish, and showed altered community structure and function. Skin and gut microbiomes were also less varied within captive populations, reflecting more uniform artificial rearing conditions. The surrounding water influenced the microbiome of the gut and, especially, the skin, but could not explain the degree of variation observed between populations. For both gut and skin, we found that there was greater difference in microbiome structures between more genetically distinct fish populations, and that population genetic diversity was positively correlated with microbiome diversity. However, diet is likely to be the major factor contributing to the large differences in gut microbiota between wild and captive fish. Our results highlight the scope of interpopulation variation in the Atlantic salmon microbiome and offer insights into the deterministic factors contributing to microbiome diversity and structure.IMPORTANCE Variation in the microbiome has a fundamental influence on host health, ecology, and evolution, but the scope and basis of this variation are not fully understood. We identified considerable variation in skin and gut microbial communities between seven wild and captive populations of Atlantic salmon, reflecting divergent environmental conditions and fish genetic diversity. In particular, we found very pronounced differences in the intestinal microbiomes of wild and hatchery-reared fish, likely reflecting differences in diet. Our results offer an insight into how the microbiome potentially contributes to the generation of local adaptations in this species and how domestication alters intestinal microbial communities, highlighting future research directions in these areas.

Interpopulation Variation in the Atlantic Salmon Microbiome Reflects Environmental and Genetic Diversity

Variation in the microbiome has a fundamental influence on host health, ecology, and evolution, but the scope and basis of this variation are not fully understood. We identified considerable variation in skin and gut microbial communities between seven wild and captive populations of Atlantic salmon, reflecting divergent environmental conditions and fish genetic diversity. In particular, we found very pronounced differences in the intestinal microbiomes of wild and hatchery-reared fish, likely reflecting differences in diet. Our results offer an insight into how the microbiome potentially contributes to the generation of local adaptations in this species and how domestication alters intestinal microbial communities, highlighting future research directions in these areas.

The microbiome has a crucial influence on host phenotype and is of broad interest to ecological and evolutionary research. Yet, the extent of variation that occurs in the microbiome within and between populations is unclear. We characterized the skin and gut microbiomes of seven populations of juvenile Atlantic salmon (Salmo salar) inhabiting a diverse range of environments, including hatchery-reared and wild populations. We found shared skin operational taxonomic units (OTUs) across all populations and core gut microbiota for all wild fish, but the diversity and structure of both skin and gut microbial communities were distinct between populations. There was a marked difference between the gut microbiomes of wild and captive fish. Hatchery-reared fish had lower intestinal microbial diversity, lacked core microbiota found in wild fish, and showed altered community structure and function. Skin and gut microbiomes were also less varied within captive populations, reflecting more uniform artificial rearing conditions. The surrounding water influenced the microbiome of the gut and, especially, the skin, but could not explain the degree of variation observed between populations. For both gut and skin, we found that there was greater difference in microbiome structures between more genetically distinct fish populations, and that population genetic diversity was positively correlated with microbiome diversity. However, diet is likely to be the major factor contributing to the large differences in gut microbiota between wild and captive fish. Our results highlight the scope of interpopulation variation in the Atlantic salmon microbiome and offer insights into the deterministic factors contributing to microbiome diversity and structure.

IMPORTANCE Variation in the microbiome has a fundamental influence on host health, ecology, and evolution, but the scope and basis of this variation are not fully understood. We identified considerable variation in skin and gut microbial communities between seven wild and captive populations of Atlantic salmon, reflecting divergent environmental conditions and fish genetic diversity. In particular, we found very pronounced differences in the intestinal microbiomes of wild and hatchery-reared fish, likely reflecting differences in diet. Our results offer an insight into how the microbiome potentially contributes to the generation of local adaptations in this species and how domestication alters intestinal microbial communities, highlighting future research directions in these areas.

Teleosts Genomics: Progress and Prospects in Disease Prevention and Control

Genome wide studies based on conventional molecular tools and upcoming omics technologies are beginning to gain functional applications in the control and prevention of diseases in teleosts fish. Herein, we provide insights into current progress and prospects in the use genomics studies for the control and prevention of fish diseases. Metagenomics has emerged to be an important tool used to identify emerging infectious diseases for the timely design of rational disease control strategies, determining microbial compositions in different aquatic environments used for fish farming and the use of host microbiota to monitor the health status of fish. Expounding the use of antimicrobial peptides (AMPs) as therapeutic agents against different pathogens as well as elucidating their role in tissue regeneration is another vital aspect of genomics studies that had taken precedent in recent years. In vaccine development, prospects made include the identification of highly immunogenic proteins for use in recombinant vaccine designs as well as identifying gene signatures that correlate with protective immunity for use as benchmarks in optimizing vaccine efficacy. Progress in quantitative trait loci (QTL) mapping is beginning to yield considerable success in identifying resistant traits against some of the highly infectious diseases that have previously ravaged the aquaculture industry. Altogether, the synopsis put forth shows that genomics studies are beginning to yield positive contribution in the prevention and control of fish diseases in aquaculture.

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.

Generalists and Specialists: A New View of How MHC Class I Molecules Fight Infectious Pathogens

In comparison with the major histocompatibility complexes (MHCs) of typical mammals, the chicken MHC is simple and compact with a single dominantly expressed class I molecule that can determine the immune response. In addition to providing useful information for the poultry industry and allowing insights into the evolution of the adaptive immune system, the simplicity of the chicken MHC has allowed the discovery of phenomena that are more difficult to discern in the more complicated mammalian systems. This review discusses the new concept that poorly expressed promiscuous class I alleles act as generalists to protect against a wide variety of infectious pathogens, while highly expressed fastidious class I alleles can act as specialists to protect against new and dangerous pathogens.

A broad overview of classical MHC I expression and bound peptides reveals an inverse correlation between repertoire breadth and cell-surface expression in some chicken and human alleles.

Several chicken class I alleles with wide peptide-binding repertoires (promiscuity) are associated with resistance to a variety of common diseases.

Conversely, a narrow peptide-binding repertoire (fastidiousness) in some human HLA-B alleles is associated with resistance to HIV progression.

Cell-surface expression of some classical class I alleles depends on the regulation of translocation to the cell surface rather than of transcription or translation. MHC translocation is influenced by peptide translocation in chickens and by tapasin interaction in humans.

Expression of immune genes in skin of channel catfish immunized with live theronts of Ichthyophthirius multifiliis

The objective of this study was to evaluate differential expression of innate and adaptive immune genes, including immunoglobulin, immune cell receptor, cytokine, inflammatory protein, toll-like receptors (TLR) and recombination-activating gene (RAG) in skin from channel catfish, Ictalurus punctatus after immunization with live theronts of Ichthyophthirius multifiliis (Ich) by intraperitoneal injection. The immunized catfish showed significantly higher survival rate (95%) than those of mock-immunized control fish (0% survival) after the theront challenge. The gene expression of innate immune system, such as cytokines (IL-1β type a, IL-1β type b, IFN-γ, TGF1-β and TNF-α) and inflammatory proteins (NF-kB and iNOS 2), showed significant upregulation at day 1 (D1) post-immunization. Expression of TLR genes exhibited a rapid increase from hour 4 (h4) to D10 post-immunization. Genes of the adaptive response, such as the cell receptor MHC I, CD8⁺ , CD4⁺ and TCR-α, showed upregulation at D1, D6 and D10. The TCR-β expression increased rapidly at h4 and remained upregulated until D10. Immunoglobulin IgM upregulation was detected from h4 until D2 while IgD expression was increased from D1 until D10. Rapid upregulation of innate and adaptive immune genes in skin of catfish following live theront vaccination was demonstrated in this study ultimately resulting in significant protection against Ich infection.

Genes of the major histocompatibility complex highlight interactions of the innate and adaptive immune system

Background

A well-functioning immune defence is crucial for fitness, but our knowledge about the immune system and its complex interactions is still limited. Major histocompatibility complex (MHC) molecules are involved in T-cell mediated adaptive immune responses, but MHC is also highly upregulated during the initial innate immune response. The aim of our study was therefore to determine to what extent the highly polymorphic MHC is involved in interactions of the innate and adaptive immune defence and if specific functional MHC alleles (FA) or heterozygosity at the MHC are more important.

Methods

To do this we used captive house sparrows (Passer domesticus) to survey MHC diversity and immune function controlling for several environmental factors. MHC class I alleles were identified using parallel amplicon sequencing and to mirror immune function, several immunological tests that correspond to the innate and adaptive immunity were conducted.

Results

Our results reveal that MHC was linked to all immune tests, highlighting its importance for the immune defence. While all innate responses were associated with one single FA, adaptive responses (cell-mediated and humoral) were associated with several different alleles.

Discussion

We found that repeated injections of an antibody in nestlings and adults were linked to different FA and hence might affect different areas of the immune system. Also, individuals with a higher number of different FA produced a smaller secondary response, indicating a disadvantage of having numerous MHC alleles. These results demonstrate the complexity of the immune system in relation to the MHC and lay the foundation for other studies to further investigate this topic.

Identification, polymorphism and expression of MHC class Iα in golden pompano, Trachinotus ovatus

The classical major histocompatibility complex class I (MHC I) plays a vital role in the immune system. In this study, we cloned and identified golden pompano (Trachinotus ovatus) MHC Iα (Trov-MHC Iα), which encodes 351 amino acid residues including a leader peptide, α1, α2, α3 domain, a transmembrane region and a cytoplasmic domain. Twenty six different sequences, which encoded various numbers of amino acid residues ranging from 348 to 354, were obtained from 12 individuals. Highly genetic polymorphism was found in the Trov-MHC Iα, especially in the α1 and α2 domains. Meanwhile, in the α1 and α2 domains, 21 positive selected positions were revealed by site models, indicating the diversity of Trov-MHC Iα may be mainly generated by positive selection. Moreover, quantitative real-time reverse transcription PCR and western blotting analyses demonstrated that Trov-MHC Iα was ubiquitously expressed in the nine tested tissues and more highly expressed in intestine, head kidney, gill, and spleen. In the head kidney and spleen, Trov-MHC Iα was significantly upregulated under LPS or poly I:C stimulation. The results of this study provide valuable insight into molecular polymorphism, evolutionary mechanism, expression and function of MHC Iα in the immune system of golden pompano.

The polymorphism of chicken-type lysozyme gene in Japanese flounder (Paralichthys olivaceus) and its association with resistance/susceptibility to Listonella anguillarum

Lysozyme is a crucially spread hydrolase in organisms that can defend against bacterial infection in innate immunity. In this study, we successfully sequenced the coding region of chicken-type lysozyme gene (PoLysC) in Paralichthys olivaceus and identified nine single nucleotide polymorphisms (SNPs). We then amplified the 2500 bp promoter region of lysozyme and identified the eight sites of polymorphisms. All SNPs were genotyped between susceptible and resistance groups after Listonella anguillarum challenge. One of these SNP sites in the codon of PoLysC was genotyped and determined to be a significant marker by analyzing its distribution in the susceptible and resistant groups. As a nonsynonymous mutation, the frequency of 140G/C genotype in the resistant group was higher (67.74%) than that in the susceptible group (32.26%). The linkage between SNP140 and polymorphisms in the promoter region was also studied. Results revealed that the frequency of haplotype CC_-536/CC_-1200/GG₁₄₀ in the resistance group was significantly higher than that in the susceptible group. The quantitative expression of lysozyme gene in the resistant group was also higher than that in the susceptible group. This finding indicated that the linkage between polymorphism -536 and -1200 sites in promoter and SNP140 in codon sequence was associated with the resistance of P. olivaceus to L. anguillarum. All these results suggest that the mutations in promoter and coding region were related to changes in PoLysC for resisting L. anguillarum. The haplotype CC_-536/CC_-1200/GG₁₄₀ was a potential marker and can thus be applied to selective breeding for the disease resistance of P. olivaceus.

Expression of immune genes in systemic and mucosal immune tissues of channel catfish vaccinated with live theronts of Ichthyophthirius multifiliis

Ichthyophthiriasis caused by Ichthyophthirius multifiliis (Ich) has a worldwide distribution and affects most freshwater fishes. Fish surviving natural infection and/or immunized with Ich develop strong innate and adaptive immune responses. However, there is a lack of the knowledge regarding immune gene expression patterns in systemic and mucosal immune tissues, and how immune genes interact and lead to innate and adaptive immune protection against Ich infection in fish. The objective of this study was to investigate the expression of innate and adaptive immune-related genes in systemic (liver, spleen) and mucosal (gill, intestine) tissues of channel catfish over time following vaccination with live Ich theronts. The vaccinated fish showed significantly higher antibody titers and survival (95%) than those of mock immunized fish. Expression of IgM and IgD heavy chain genes exhibited a rapid increase from 4 h (h4) to 2 days (d2) post-vaccination in systemic immune tissues. Immune cell receptor genes (CD4, CD8-α, MHC I, MHC II β, TcR-α, and TcR-β) were more highly upregulated and remained upregulated for longer duration in systemic tissues than in mucosal tissues of the vaccinated fish. The cytokine genes IL-1βa and IFN-γ were rapidly upregulated in both systemic and mucosal tissues of vaccinated fish, with peak expression from h4 to d1 post-vaccination. Toll-like receptor genes TLR-1 and TLR-9 showed relatively stable upregulation in the gill of immunized fish following vaccination. Results of this study revealed the molecular immune responses in mucosal and systemic tissues of vaccinated fish and demonstrated that Ich vaccination resulted in innate and adaptive immune responses against Ich infection.

Temporal and spatial instability in neutral and adaptive (MHC) genetic variation in marginal salmon populations

The role of marginal populations for the long-term maintenance of species’ genetic diversity and evolutionary potential is particularly timely in view of the range shifts caused by climate change. The Centre-Periphery hypothesis predicts that marginal populations should bear reduced genetic diversity and have low evolutionary potential. We analysed temporal stability at neutral microsatellite and adaptive MHC genetic variation over five decades in four marginal Atlantic salmon populations located at the southern limit of the species’ distribution with a complicated demographic history, which includes stocking with foreign and native salmon for at least 2 decades. We found a temporal increase in neutral genetic variation, as well as temporal instability in population structuring, highlighting the importance of temporal analyses in studies that examine the genetic diversity of peripheral populations at the margins of the species’ range, particularly in face of climate change.

Genomic structure and expression pattern of MHC IIα and IIβ genes reveal an unusual immune trait in lined seahorse Hippocampus erectus

The major histocompatibility complex (MHC) genes are crucial in the adaptive immune system, and the gene duplication of MHC in animals can generally result in immune flexibility. In this study, we found that the lined seahorse (Hippocampus erectus) has only one gene copy number (GCN) of MHC IIα and IIβ, which is different from that in other teleosts. Together with the lack of spleen and gut-associated lymphatic tissue (GALT), the seahorse may be referred to as having a partial but natural "immunodeficiency". Highly variable amino acid residues were found in the IIα and IIβ domains, especially in the α1 and β1 domains with 9.62% and 8.43% allelic variation, respectively. Site models revealed seven and ten positively selected positions in the α1 and β1 domains, respectively. Real-time PCR experiments showed high expression levels of the MHC II genes in intestine (In), gill (Gi) and trunk kidney (TK) and medium in muscle (Mu) and brood pouch (BP), and the expression levels were significantly up-regulated after bacterial infection. Specially, relative higher expression level of both MHC IIα and IIβ was found in Mu and BP when compared with other fish species, in which MHC II is expressed negligibly in Mu. These results indicate that apart from TK, Gi and In, MU and BP play an important role in the immune response against pathogens in the seahorse. In conclusion, high allelic variation and strong positive selection in PBR and relative higher expression in MU and BP are speculated to partly compensate for the immunodeficiency.

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.

Does density influence relative growth performance of farm, wild and F1 hybrid Atlantic salmon in semi-natural and hatchery common garden conditions?

The conditions encountered by Atlantic salmon, Salmo salar L., in aquaculture are markedly different from the natural environment. Typically, farmed salmon experience much higher densities than wild individuals, and may therefore have adapted to living in high densities. Previous studies have demonstrated that farmed salmon typically outgrow wild salmon by large ratios in the hatchery, but these differences are much less pronounced in the wild. Such divergence in growth may be explained partly by the offspring of wild salmon experiencing higher stress and thus lower growth when compared under high-density farming conditions. Here, growth of farmed, wild and F1 hybrid salmon was studied at contrasting densities within a hatchery and semi-natural environment. Farmed salmon significantly outgrew hybrid and wild salmon in all treatments. Importantly, however, the reaction norms were similar across treatments for all groups. Thus, this study was unable to find evidence that the offspring of farmed salmon have adapted more readily to higher fish densities than wild salmon as a result of domestication. It is suggested that the substantially higher growth rate of farmed salmon observed in the hatchery compared with wild individuals may not solely be caused by differences in their ability to grow in high-density hatchery scenarios.

Identification of 48 full-length MHC-DAB functional alleles in miiuy croaker and evidence for positive selection

Major histocompatibility complex (MHC) molecules play a vital role in the immune response and are a highly polymorphic gene superfamily in vertebrates. As the molecular marker associated with polymorphism and disease susceptibility/resistance, the polymorphism of MHC genes has been investigated in many tetrapods and teleosts. Most studies were focused on the polymorphism of the second exon, which encodes the peptide-binding region (PBR) in the α1- or β1-domain, but few studies have examined the full-length coding region. To comprehensive investigate the polymorphism of MHC gene, we identified 48 full-length miiuy croaker (Miichthys miiuy) MHC class IIB (Mimi-DAB) functional alleles from 26 miiuy croaker individuals. All of the alleles encode 34 amino acid sequences, and a high level of polymorphism was detected in Mimi-DAB alleles. The rate of non-synonymous substitutions (dN) occurred at a significantly higher frequency than that of synonymous substitutions (dS) in the PBR, and this result suggests that balancing selection maintains polymorphisms at the Mimi-DAB locus. Phylogenetic analysis based on the full-length and exon 2 sequences of Mimi-DAB alleles both showed that the Mimi-DAB alleles were clustered into two major groups. A total of 19 positive selected sites were identified on the Mimi-DAB alleles after testing for positive selection, and 14 sites were predicted to be associated with antigen-binding sites, which suggests that most of selected sites are significant for disease resistance. The polymorphism of Mimi-DAB alleles provides an important resource for analyzing the association between the polymorphism of MHC gene and disease susceptibility/resistance, and for researching the molecular selective breeding of miiuy croaker with enhanced disease resistance.