MHC class IIalpha gene polymorphism and its association with resistance/susceptibility to Vibrio anguillarum in Japanese flounder (Paralichthys olivaceus)

Characterization, expression, and polymorphism of MHC II α and MHC II β in Sichuan taimen (Hucho bleekeri)

The major histocompatibility complex (MHC) is involved in antigen presentation and plays an essential role in regulating immune function. In the present study, we identified two MHC class II genes and investigated their potential roles in Hucho bleekeri. The MHC II α and MHC II β of H. bleekeri had typical leading peptides, extracellular domains, connecting peptides, transmembrane region, and cytoplasmic region. Amino acid sequence comparison revealed that MHC II of H. bleekeri has high homology with other vertebrates, among which homology with salmonid fish was the highest. Phylogenetic analysis showed that H. bleekeri MHC II clustered with salmonid fish; moreover they clustered with orthologous genes of other fish, whereas mammalian MHC II clustered into a separate branch. Tissue distribution analysis revealed MHC II was widely expressed in all tested tissues, with both MHC II α and MHC II β highly expressed in the spleen, gill, kidney, and hindgut. After lipopolysaccharide (LPS) and polyinosinic:polycytidylic acid (poly(I:C)) stimulation, the expression of MHC II in the head kidney and spleen of H. bleekeri was significantly upregulated. Compared with MHC II α, MHC II β acted faster in response to the stimulation. Polymorphism analysis of MHC II revealed that all the different alleles belonged to the same major type, and very limited polymorphisms were found in H. bleekeri MHC II α and II β. Selection pressure analysis showed signs of weak and non-significant positive selection in the MHC II α and MHC II β extracellular region. Our study reveals the potential role of MHC II in the immune response of H. bleekeri and provides a reference for studying the evolutionary model of teleost MHC II.

Genome-Wide Association Study of Resistance to Largemouth Bass Ranavirus (LMBV) in Micropterus salmoides

The disease caused by Largemouth bass ranavirus (LMBV) is one of the most severe viral diseases in largemouth bass (Micropterus salmoides). It is crucial to evaluate the genetic resistance of largemouth bass to LMBV and develop markers for disease-resistance breeding. In this study, 100 individuals (45 resistant and 55 susceptible) were sequenced and evaluated for resistance to LMBV and a total of 2,579,770 variant sites (SNPs-single-nucleotide polymorphisms (SNPs) and insertions-deletions (InDels)) were identified. A total of 2348 SNPs-InDels and 1018 putative candidate genes associated with LMBV resistance were identified by genome-wide association analyses (GWAS). Furthermore, GO and KEGG analyses revealed that the 10 candidate genes (MHC II, p38 MAPK, AMPK, SGK1, FOXO3, FOXO6, S1PR1, IL7R, RBL2, and GADD45) were related to intestinal immune network for IgA production pathway and FoxO signaling pathway. The acquisition of candidate genes related to resistance will help to explore the molecular mechanism of resistance to LMBV in largemouth bass. The potential polymorphic markers identified in this study are important molecular markers for disease resistance breeding in largemouth bass.

Genome-Wide Identification of Trachinotus ovatus Antimicrobial Peptides and Their Immune Response against Two Pathogen Challenges

Golden pompano, Trachinotus ovatus, as a highly nutritious commercially valuable marine fish, has become one of the preferred species for many fish farmers due to its rapid growth, wide adaptability, and ease of feeding and management. However, with the expansion of aquaculture scale, bacterial and parasitic diseases have also become major threats to the golden pompano industry. This study, based on comparative genomics, shows the possibility of preferential evolution of freshwater fish over marine fish by analyzing the phylogenetic relationships and divergence times of 14 marine fish and freshwater fish. Furthermore, we identified antimicrobial peptide genes from 14 species at the genomic level and found that the number of putative antimicrobial peptides may be related to species evolution. Subsequently, we classified the 341 identified AMPs from golden pompano into 38 categories based on the classification provided by the APD3. Among them, TCP represented the highest proportion, accounting for 23.2% of the total, followed by scolopendin, lectin, chemokine, BPTI, and histone-derived peptides. At the same time, the distribution of AMPs in chromosomes varied with type, and covariance analysis showed the frequency of its repeat events. Enrichment analysis and PPI indicated that AMP was mainly concentrated in pathways associated with disease immunity. In addition, our transcriptomic data measured the expression of putative AMPs of golden pompano in 12 normal tissues, as well as in the liver, spleen, and kidney infected with Streptococcus agalactiae and skin infected with Cryptocaryon irritans. As the infection with S. agalactiae and C. irritans progressed, we observed tissue specificity in the number and types of responsive AMPs. Positive selection of AMP genes may participate in the immune response through the MAPK signaling pathway. The genome-wide identification of antimicrobial peptides in the golden pompano provided a complete database of potential AMPs that can contribute to further understanding the immune mechanisms in pathogens. AMPs were expected to replace traditional antibiotics and be developed into targeted drugs against specific bacterial and parasitic pathogens for more precise and effective treatment to improve aquaculture production.

Construction of an alanine dehydrogenase gene deletion strain for vaccine development against Nocardia seriolae in hybrid snakehead (Channa maculata ♀ × Channa argus ♂)

Nocardia seriolae is the main pathogen of fish nocardiosis. In our previous study, alanine dehydrogenase was identified as a potential virulence factor of N. seriolae. On the basis of this fact, the alanine dehydrogenase gene of N. seriolae (NsAld) was knocked out to establish the strain ΔNsAld for vaccine development against fish nocardiosis in this study. The LD₅₀ of strain ΔNsAld was 3.90 × 10⁵ CFU/fish, higher than that of wild strain (5.28 × 10⁴ CFU/fish) significantly (p < 0.05). When the strain ΔNsAld was used as a live vaccine to immunize hybrid snakehead (Channa maculata ♀ × Channa argus ♂) at 2.47 × 10⁵ CFU/fish by intraperitoneal injection, the non-specific immune indexes (LZM, CAT, AKP, ACP and SOD activities), specific antibody (IgM) titers and several immune-related genes (CD4, CD8α, IL-1β, MHCIα, MHCIIα and TNFα) were up-regulated in different tissues, indicating that this vaccine could induce humoral and cell-mediated immune responses. Furthermore, the relative percentage survival (RPS) of ΔNsAld vaccine was calculated as 76.48% after wild N. seriolae challenge. All these results suggest that the strain ΔNsAld could be a potential candidate for live vaccine development to control fish nocardiosis in aquaculture.

Development and Evaluation of Recombinant B-Cell Multi-Epitopes of PDHA1 and GAPDH as Subunit Vaccines against Streptococcus iniae Infection in Flounder (Paralichthys olivaceus)

Streptococcus iniae is a severe Gram-positive pathogen that can infect a wide range of freshwater and marine fish species. In continuation of our earlier studies on the development of S. iniae vaccine candidates, pyruvate dehydrogenase E1 subunit alpha (PDHA1) and glyceraldehyde-3-phosphate dehydrogenase (GAPDH) were highly efficacious in protecting flounder (Paralichthys olivaceus) against S. iniae. In the present study, to investigate the potential of multi-epitope vaccination strategy to prevent flounder against S. iniae infection, the liner B-cell epitopes of PDHA1 and GAPDH proteins were predicted using a bioinformatics approach and were identified by immunoassay, and recombinant B-cell multi-epitopes of PDHA1 and GAPDH (rMEPIP and rMEPIG) containing immunodominant epitope-concentrated domains were expressed in Escherichia coli BL21 (DE3) and were used as a subunit vaccine to immunize healthy flounder, while recombinant PDHA1 (rPDHA1), GAPDH (rGAPDH) and formalin-inactivated S. iniae (FKC) served as controls. Then, the immunoprotection efficacy of rMEPIP and rMEPIG was evaluated by determining the percentages of CD4-1+, CD4-2+, CD8β+ T lymphocytes and surface-IgM-positive (sIgM+) lymphocytes in peripheral blood leucocytes (PBLs), spleen leucocytes (SPLs) and head kidney leucocytes (HKLs), as well as total IgM, specific IgM, and relative percentage survival (RPS) post immunization, respectively. It was found that fish immunized with rPDHA1, rGAPDH, rMEPIP, rMEPIG and FKC showed significant increases in sIgM+, CD4-1+, CD4-2+, and CD8β+ lymphocytes and production of total IgM and specific IgM against S. iniae or recombinant proteins rPDHA1 and rGAPDH, which indicated the activation of humoral and cellular immune responses after vaccination. Moreover, RPS rate of the multi-epitope vaccine rMEPIP and rMEPIG groups reached 74.07% and 77.78%, higher than that of rPDHA1 and rGAPDH (62.96% and 66.67%) and KFC (48.15%). These results demonstrated that B-cell multi-epitope protein vaccination, rMEPIP and rMEPIG, could give a better protective effect against S. iniae infection, which provided a promising strategy to design the efficient vaccine in teleost fish.

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.

Induction of attenuated Nocardia seriolae and their use as live vaccine trials against fish nocardiosis

Nocardia seriolae, a Gram-positive facultative intercellular pathogen, has been identified as the causative agent of fish nocardiosis, causing substantial mortality and morbidity of a wide range of fish species. Looking into that fact, the effective vaccine against this pathogen is urgently needed to control the significant losses in aquaculture practices. In order to induct attenuated strains for developing the potential live vaccines, the mutagenic N. seriolae strain S-250 and U-20 were obtained from wild-type strain ZJ0503 through continuous passaging and ultraviolet (UV) irradiation, respectively. Additionally, the biological characteristic, virulence, stability, mediating immune response and supplying protective efficacy to hybrid snakehead of the S-250 and U-20 strains were determined in the present study. The results showed that U-20 strain displayed dramatic changes in morphological characteristic and significant decreased in the virulence to hybrid snakehead, while that of S-250 strain had no obvious different in comparison to ZJ0503 strain. When hybrid snakehead were intraperitoneally injected with ZJ0503, S-250 and U-20 strains at their respective sub-clinical dosage, the non-specific immunity parameters (serum LYZ, POD, ACP, AKP and SOD activities), specific antibody (IgM) titers production and immune-related genes (CC1, CC2, IL-1β, IL-8, TNFα, IFNγ, MHCIα, MHCIIα, CD4, CD8α, TCRα and TCRβ) expression were up-regulated, indicating that they were able to trigger humoral and cell-mediated immune responses. Furthermore, the protective efficacy in hybrid snakehead after vaccination with ZJ0503, S-250 and U-20 strains, in terms of relative percentage survival (RPS), were 28.85%, 56.89% and 89.65% respectively. Taken together, two attenuated N. seriolae strains S-250 and U-20 were obtained successfully and they could elicit strong immune response and supply protective efficacy to hybrid snakehead against N. seriolae, which suggested that these two attenuated strains were the potential candidates for live vaccine development to control fish nocardiosis in aquaculture.

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.

Transcriptome analysis and protein-protein interaction in resistant and susceptible families of Japanese flounder (Paralichthys olivaceus) to understand the mechanism against Edwardsiella tarda

Edwardsiella tarda is one of the most harmful bacterial pathogens for aquaculture flatfish. After artificial infection of 47 Japanese flounder (Paralichthys olivaceus) families, resistant and susceptible families were identified in this study. High-throughput sequencing was performed on the liver transcriptome of uninfected groups (PoRU and PoSU) and infected groups (PoRC and PoSC). Through assembly and annotation, a total of 3012 and 1386 differentially expressed genes (DEGs) were identified in PoRU vs. PoSU and PoRC vs. PoSC. The significant enrichment pathways between PoRU and PoSU were mainly in metabolic and biosynthesis pathways. A total of thirty dominant enrichment pathways between PoRC and PoSC mainly focused on some immune-related pathways, including the hematopoietic cell lineage, cytokine-cytokine receptor interaction, complement and coagulation cascades, antigen processing and presentation, the intestinal immune network for immunoglobulin A (IgA) production and T/B cell receptor signaling pathway. Under the protein-protein interaction (PPI) analysis, hub genes, including CD molecules, complement component factors and chemokines, were identified in the network, and 16 core genes were differentially expressed in resistant and sustainable families in quantitative polymerase chain reaction (qPCR) validation. This study represents the first transcriptome analysis based on resistant and susceptible families and provides resistant genes to understand the potential molecular mechanisms of antibacterial function in marine fish. The results obtained in this study provide crucial information on gene markers for resistant breeding of Japanese flounder.

Immunogenicity study of OmpU subunit vaccine against Vibrio mimicus in yellow catfish, Pelteobagrus fulvidraco

The outer membrane protein U (OmpU) is a conserved outer membrane protein in a variety of pathogenic Vibrio species and has been considered as a vital protective antigen for vaccine development. Vibrio mimicus (V. mimicus) is the pathogen causing ascites disease in aquatic animals. In this study, the prokaryotically expressed and purified His-tagged OmpU of V. mimicus (His-OmpU) was used as a subunit vaccine. The formalin inactivated V. mimicus, purified His tag (His-tag), and PBS were used as controls. The vaccinated yellow catfish were challenged with V. mimicus at 28 days post-vaccination, and the results showed that the His-OmpU and inactivated V. mimicus groups exhibited much higher survival rates than the His-tag and PBS groups. To fully understand the underlying mechanism, we detected the expression levels of several immune-related genes in the spleen of fish at 28 days post-vaccination and 24 h post-challenge. The results showed that most of the detected immune-related genes were significantly upregulated in His-OmpU and inactivated V. mimicus groups. In addition, we performed the serum bactericidal activity assay, and the results showed that the serum from His-OmpU and inactivated V. mimicus groups exhibited much stronger bactericidal activity against V. mimicus than those of His-tag and PBS groups. Finally, the serum agglutination antibody was detected, and the antibody could be detected in His-OmpU and inactivated V. mimicus groups with the antibody titers increasing along with the time post-vaccination, but not in His-tag or PBS group. Our data reveal that the recombinant OmpU elicits potent protective immune response and is an effective vaccine candidate against V. mimicus in yellow catfish.

A SNP of the hemocyanin gene (LvHc) is a marker for high growth and ammonia-tolerance in Pacific white shrimp Litopenaeus vannamei

Expression levels of hemocyanin (LvHc), activating transcription factor 4 (LvAtf4), glutathione S-transferase (LvGst), caspase 2 (LvCasp2) and anti-lipopolysaccharide factor (LvAlf) were examined in the hepatopancreas of Pacific white shrimp Litopenaeus vannamei juveniles exposed to a lethal concentration of ammonia-N (32.15 mg/l). The expression levels of all transcripts except LvAlf were significantly greater (P < 0.05) in tolerant shrimp (Lv-AT; N = 30) that survived up to 72 h post treatment (hpt) than in susceptible shrimp (Lv-AS24 and Lv-AS72; N = 45 and 15), that died within 24 h or between 24 and 72 hpt, respectively. Subsequently, effects of non-lethal concentrations of ammonia-N (control, 10 and 20 mg/l) on the expression of LvHc in juvenile shrimp were examined. Compared to the control, expression levels of LvHc transcripts in hemocytes and the hepatopancreas of tested shrimp changed after exposure to ammonia-N. One SNP (C > T₅₄₅) was found in the LvHc₃₂₂ gene segment. Real-time PCR amplification of specific alleles (real-time PASA) was developed for detection of C > T₅₄₅ genotypes. Juveniles in the lethal exposure test that carried a C/T₅₄₅ genotype showed a greater average body weight and total length (8.46 ± 0.36 g and 10.05 ± 0.16 cm) than those with a C/C₅₄₅ genotype (7.48 ± 0.31 g and 9.60 ± 0.13 cm) (P < 0.05). Similar results were found in the second generation (G2) of a growth-improved stock (3 and 4 families of BIOTEC-G2-L1 and BIOTEC-G2-L2) and in commercially farmed shrimp (2 groups). Accordingly, expression levels and SNP of LvHc can serve as markers for selection high growth performance in ammonia-tolerant L. vannamei.

Identification of a Growth-Associated Single Nucleotide Polymorphism (SNP) in Cyclin C of the Giant Tiger Shrimp Penaeus monodon

The full-length cDNA of cyclin C of the giant tiger shrimp Penaeus monodon (PmCyC) was isolated by RACE-PCR. It was 1443 bp in length containing an open reading frame (ORF) of 804 bp and 267 deduced amino acids. Tissue distribution analysis indicated that PmCyC was more abundantly expressed in ovaries and testes than other tissues of female and male juveniles (P < 0.05). A pair of primers was designed, and an amplification product of 403 bp containing an intron of 123 bp was obtained. Polymorphism of amplified PmCyC gene segments of the 5th (3-month-old G5, N = 30) and 7th (5-month-old G7, N = 18) generations of domesticated juveniles was analyzed. Four conserved SNPs (T>C₁₃₄, T>C₁₈₈, G>A₃₇₉, and T>C₃₈₂) were found within the examined sequences. A TaqMan genotyping assay was developed for detection of a T>C₁₃₄ SNP. Association analysis indicated that this SNP displayed significant association with body weight (P < 4.2e-10) and total length (P < 2e-09) of the examined G7 P. monodon (N = 419) with an allele substitution effect of 5.02 ± 0.78 g and 1.41 ± 0.19 cm, respectively. Juveniles with C/C₁₃₄ (22.80 ± 2.51 g and 12.97 ± 0.53 cm, N = 19) and T/C₁₃₄ (20.41 ± 0.93 g and 12.77 ± 0.21 cm, N = 129) genotypes exhibited a significantly greater average body weight and total length than those with a T/T₁₃₄ genotype (14.72 ± 0.53 g and 11.37 ± 0.13 cm, N = 271) (P < 0.05).

Recombinant outer membrane protein C of Aeromonas salmonicida subsp. masoucida, a potential vaccine candidate for rainbow trout (Oncorhynchus mykiss)

Aeromonas salmonicida subsp. masoucida (ASM) is an important bacterial pathogen of salmonid fish, which can cause huge economic losses to the fish farming industry. In order to screen effective vaccine candidate proteins, four outer membrane proteins of ASM, including OmpA, OmpC, OmpK and OmpW, were selected and recombinantly expressed in Escherichia coli. The result of western blotting showed that these four recombinant proteins could be recognized by rainbow trout anti-ASM antibodies. The immune protective effects of the four rOMPs were also investigated, and the relative percentage survival (RPS) of rOmpA, rOmpC, rOmpK and rOmpW were 71.1%, 81.6%, 55.3% and 42.1%, respectively. The RPS of rOmpC was significantly higher than the other three rOMPs, so the immune responses of rainbow trout induced by rOmpC were further investigated. The results showed that vaccination with rOmpC could significantly induced the production of specific serum antibodies and proliferation of sIg + lymphocytes in peripheral blood. Meanwhile, RT-qPCR analysis showed that rOmpC could significantly enhance the expression of the MHC-II, TCR, CD4, CD8, IL-8 and IgM genes compared with the BSA immunized group. These results demonstrated that rOmpC could induce strong humoral immune response in rainbow trout and provided effective immune protection against ASM challenge, which indicated that OmpC is a promising vaccine candidate against Aeromonas salmonicida infection.

Molecular polymorphism and expression of MHC I α, II α, II β and II invariant chain in the critically endangered Dabry's sturgeon (Acipenser dabryanus)

The major histocompatibility complex (MHC) is a key player in the regulation of immune responses through presenting foreign antigens to T lymphocytes. In this study, three MHC genes, namely, MHC I α, II α, II β and the II invariant chain (Ii), were identified and characterized in the critically endangered Dabry's sturgeon (Acipenser dabryanus). A tissue distribution study showed that the MHC and Ii transcripts were widely expressed in various tissues. The highest expression levels of MHC I α, II α and Ii were found in the gill, while MHC II β was primarily expressed in the spleen. Challenge of A. dabryanus with a pathogenic bacterium in vivo resulted in significant upregulation of both MHC and Ii expression, indicating potential roles of these genes in immune response. Phylogenetic analysis showed that A. dabryanus MHC grouped with other teleost MHC genes and sequences from Polyodon spathula and A. dabryanus had an intermingling of alleles. According to the split time between paddlefishes and sturgeons, this result indicated that trans-species MHC lineages in Chondrostei were much older than those in tetrapods. The molecular polymorphisms of the complete open reading frame regions of the MHC genes were analysed in several A. dabryanus individuals. MHC II α and II β were highly polymorphic in different individuals, while MHC I α was more conserved. The ratio of non-synonymous substitution occurred at a significantly higher frequency than synonymous substitution in peptide-binding regions (PBR) of MHC II α and II β, demonstrating the existence of positive selection at peptide-binding sites. Our study suggested potential roles of the MHC chains in immune response to pathogen microbial infection, and the numerous alleles identified in this study will help further genetic management and molecular marker-assisted selective breeding programmes in A. dabryanus.

T and B lymphocytes immune responses in flounder (Paralichthys olivaceus) induced by two forms of outer membrane protein K from Vibrio anguillarum: Subunit vaccine and DNA vaccine

To further elucidate the roles of T and B lymphocytes in fish subunit and DNA candidate vaccines for immunisation, the immune responses of T and B lymphocytes to recombinant protein (rOmpK) and plasmid OmpK (pOmpK) from Vibrio anguillarum plus cyclosporine A (CsA) were investigated in flounder (Paralichthys olivaceus). The results showed that in the rOmpK-immunised groups, the percentages of CD4-1⁺ and CD4-2⁺ T (PCD4-1⁺ and PCD4-2⁺ T) lymphocytes significantly increased to a peak on days 5 or 7. The percentages of IgM⁺ B (PIgM⁺ B) lymphocytes and specific antibodies markedly increased to a peak at weeks 4 or 5. The nine immune-related genes were significantly up-regulated and the expression levels of CD4-1, CD4-2 and MHC II genes were higher than that of CD8α, CD8β and MHC I genes. The CD4⁺ T lymphocytes, IgM⁺ B lymphocytes, and specific antibodies were significantly inhibited by CsA. Therefore, the responses of CD4⁺ T lymphocytes influenced the responses of the B lymphocytes and antibodies. In the pOmpK-immunised groups, the PCD4-1⁺, PCD4-2⁺, and PCD8β⁺ T lymphocytes significantly increased to a peak on days 11 or 14, days 9 or 11, and days 7 or 9, respectively. The PIgM⁺ B lymphocytes and specific antibodies significantly increased to a peak at weeks 5 or 6. Immune related genes upregulated, and CD4⁺ and CD8⁺ T lymphocytes, IgM⁺ B lymphocytes and specific antibodies all suppressed by CsA, suggesting that the responses of T lymphocytes subpopulations influenced B lymphocytes and antibodies responses. Therefore, the subpopulations of T lymphocytes played an important role in the immune responses induced by subunit and DNA candidate vaccines of OmpK and regulated the immune responses of B lymphocytes in flounder.

Identification of common antigens of three pathogenic Nocardia species and development of DNA vaccine against fish nocardiosis

Fish nocardiosis is a chronic granulomatous bacterial disease and three pathogens have been reported so far, including Nocardia asteroids, N. seriolae and N. salmonicida. However, the absence of antigen markers is a bottleneck for developing effective vaccines against fish nocardiosis. In this study, the antigenicity of whole-cell protein of these three pathogenic Nocardia species were profiled by immunoproteomic analysis and 7 common immunogenic proteins were identified as follows: molecular chaperone DnaK (DnaK), molecular chaperone GroEL (GroEL), 30 S ribosomal protein S1 (RpsA), TerD family protein (TerD), FHA domain-containing protein (FHA), 50 S ribosomal protein L7/L12 (RplL) and PspA/IM30 family protein (PspA). Furthermore, the DNA vaccine encoding FHA gene against fish nocardiosis was developed and its efficacy was investigated in hybrid snakehead. The results suggested that it needed at least 7 d to transport pcDNA-FHA DNA vaccine from injected muscle to head kidney, spleen and liver and stimulate host's immune system for later protection. In addition, non-specific immunity paraments (serum lysozyme (LYZ), peroxidase (POD), acid phosphatase (ACP), alkaline phosphatase (AKP) and superoxide dismutase (SOD) activities), specific antibody (IgM) titers production and immune-related genes (MHCIα, MHCIIα, CD4, CD8α, IL-1β and TNFα) were used to evaluate the immune response induced in pcDNA-FHA vaccinated hybrid snakehead, it proved that all these mentioned immune activities were significantly enhanced after immunization. The results also showed hybrid snakehead vaccinated with pcDNA-FHA had higher survival rate (79.33%) compared with the controls after challenge with N. seriolae, indicating that the pcDNA-FHA DNA vaccine can supply immune protection against N. seriolae infection. Taken together, this study may warrant further development of these common immunogenic proteins as the antigens for vaccine or diagnosis and facilitate the prevention and treatment of fish nocardiosis.

Comparison of protective efficacy between two DNA vaccines encoding DnaK and GroEL against fish nocardiosis

Fish nocardiosis is a chronic granulomatous bacterial disease mainly caused by three pathogenic bacteria, including Nocardia seriolae, N. asteroids and N. salmonicida. Molecular chaperone DnaK and GroEL were identified to be the common antigens of the three pathogenic Nocardia species in our previous studies. To evaluate the immune protective effect of two DNA vaccines encoding DnaK or GroEL against fish nocardiosis, hybrid snakehead were vaccinated and the immune responses induced by these two vaccines were comparatively analyzed. The results suggested it needed at least 7 d to transport DnaK or GroEL gene from injected muscle to head kidney, spleen and liver and stimulate host's immune system for later protection after immunization by DNA vaccines. Additionally, non-specific immunity parameters (serum lysozyme (LYZ), peroxidase (POD), acid phosphatase (ACP), alkaline phosphatase (AKP) and superoxide dismutase (SOD) activities), specific antibody (IgM) production and immune-related genes (MHCIα, MHCIIα, CD4, CD8α, IL-1β and TNFα) were used to evaluate the immune responses induced in vaccinated hybrid snakehead. It proved that all the above-mentioned immune activities were significantly enhanced after immunization with these two DNA vaccines. The protective efficacy of pcDNA-DnaK and pcDNA-GroEL DNA vaccines, in terms of relative percentage survival (RPS), were 53.01% and 80.71% respectively. It demonstrated that these two DNA vaccines could increase the survival rate of hybrid snakehead against fish nocardiosis, albeit with variations in immunoprotective effects. Taken together, these results indicated that both pcDNA-DnaK and pcDNA-GroEL DNA vaccines could boost the innate, humoral and cellular immune response in hybrid snakehead and show highly protective efficacy against fish nocardiosis, suggesting that DnaK and GroEL were promising vaccine candidates. These findings will promote the development of DNA vaccines against fish nocardiosis in aquaculture.

Divergence between genes but limited allelic polymorphism in two MHC class II A genes in Leach's storm-petrels Oceanodroma leucorhoa

The major histocompatibility complex (MHC) is critical to host-pathogen interactions. Class II MHC is a heterodimer, with α and β subunits encoded by different genes. The peptide-binding groove is formed by the first domain of both subunits (α₁ and β₁), but studies of class II variation or natural selection focus primarily on the β subunit and II B genes. We explored MHC II A in Leach's storm-petrel, a seabird with two expressed, polymorphic II B genes. We found two II A genes, Ocle-DAA and Ocle-DBA, in contrast to the single II A gene in chicken and duck. In exon 2 which encodes the α₁ domain, the storm-petrel II A genes differed strongly from each other but showed little within-gene polymorphism in 30 individuals: just one Ocle-DAA allele, and three Ocle-DBA alleles differing from each other by single non-synonymous substitutions. In a comparable sample, the two II B genes had nine markedly diverged alleles each. Differences between the α₁ domains of Ocle-DAA and Ocle-DBA showed signatures of positive selection, but mainly at non-peptide-binding site (PBS) positions. In contrast, positive selection within and between the II B genes corresponded to putative PBS codons. Phylogenetic analysis of the conserved α₂ domain did not reveal deep or well-supported lineages of II A genes in birds, in contrast to the pronounced differentiation of DQA, DPA, and DRA isotypes in mammals. This uncertain homology complicates efforts to compare levels of functional variation and modes of evolution of II A genes across taxa.

Genetic parameters for resistance against Vibrio anguillarum in turbot Scophthalmus maximus

Genetic parameters for resistance to Vibrio anguillarum in Scophthalmus maximus were estimated using three different statistical models. Data were recorded from an experimental infection performed on 2,400 individuals from 30 full-sib groups. Cross-sectional linear model and cross-sectional threshold probit model were used to analyse the test-period survival; the cross-sectional threshold logit models were used to analyse the test-day survival. The heritability values estimated by cross-sectional linear model (CSL), cross-sectional threshold (probit) model (THRp) and cross-sectional threshold (logit) model (THRl) were 0.202 ± 0.101, 0.296 ± 0.168 and 0.110 ± 0.023, respectively. The correlation coefficients between the full-sib families' estimated breeding values (EBVs) based on these three models were over 0.993. An almost identical ranking of families was generated using any of these models. Accuracy of selection using CSL, THRp and THRl models was 0.783, 0.789 and 0.801, respectively. Accuracy of selection based on the THRl model was higher than that based on CSL and THRp models.

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.

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.

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.

The effects of IL-1β, IL-8, G-CSF and TNF-α as molecular adjuvant on the immune response to an E. tarda subunit vaccine in flounder (Paralichthys olivaceus)

Cytokines play vital roles in mounting immune responses and activating host defense network. In this study, the expression plasmid pcDNA3.1 (pcN3) encoding four flounder (Paralichthys olivaceus) cytokines including IL-1β, TNF-α, IL-8 or G-CSF (pcIL-1β, pcTNF-α, pcIL-8 and pcG-CSF) were successfully constructed, and their adjuvant potential on an Edwardsiella tarda (E. tarda) subunit vaccine OmpV (rOmpV) were comparatively analyzed in vaccinated flounder model. Results revealed that flounder vaccinated with rOmpV plus pcIL-1β, pcIL-8 or pcG-CSF produced the relative percent survivals (RPS) of 71%, 65% and 49% respectively, which were higher than that in flounder vaccinated with rOmpV plus pcTNF-α (39%) or pcN3 (36%, the control group). Immunological analysis showed that: (1) except pcTNF-α, higher levels of anti-E. tarda serum antibodies and sIg + lymphocytes in spleen, head kidney and peripheral blood were significantly enhanced by pcIL-1β, pcIL-8 or pcG-CSF, however, pcIL-8 and pcIL-1β enhanced higher levels of sIg + lymphocytes and anti-E. tarda antibodies than pcG-CSF; (2) pcTNF-α could promote the up-regulation of genes participated in cellular immunity (MHCIα, IFN-γ, CD8α and CD8β), pcIL-1β could enhance the expression of genes related to humoral immunity (CD4-1, CD4-2, MHCIIα and IgM), and all the detected genes were augmented by pcIL-8 and pcG-CSF; Among the four cytokines, pcIL-8 and pcIL-1β could strengthen the highest levels of genes participated in cellular immunity and humoral immunity, respectively. These results demonstrated that pcIL-8 and pcIL-1β could enhance stronger cellular and/or humoral immunity induced by rOmpV than pcG-CSF and pcTNF-α, and evoked higher RPS against E. tarda challenge in flounder, which indicated that pcIL-8 and pcIL-1β are promising adjuvants of vaccines in controlling E. tarda infection.

Construction of a Vibrio alginolyticus hopPmaJ (hop) mutant and evaluation of its potential as a live attenuated vaccine in orange-spotted grouper (Epinephelus coioides)

Vibrio alginolyticus, a bacterial pathogen in fish and humans, expresses a type III secretion system (T3SS) that is critical for pathogen virulence and disease development. However, little is known about the associated effectors (T3SEs) and their physiological role. In this study, the T3SE gene hopPmaJ (hop) was cloned from V. alginolyticus wild-type strain HY9901 and the mutant strain HY9901Δhop was constructed by the in-frame deletion method. The results showed that the deduced amino acid sequence of V. alginolyticus HopPmaJ shared 78-98% homology with other Vibrio spp. In addition, the HY9901Δhop mutant showed an attenuated swarming phenotype and a 2600-fold decrease in the virulence to grouper. However, the HY9901Δhop mutant showed no difference in morphology, growth, biofilm formation and ECPase activity. Finally, grouper vaccinated via intraperitoneal (IP) injection with HY9901Δhop induced a high antibody titer with a relative percent survival (RPS) value of 84% after challenging with the wild-type HY9901. Real-time PCR assays showed that vaccination with HY9901Δhop enhanced the expression of immune-related genes, including MHC-Iα, MHC-IIα, IgM, and IL-1β after vaccination, indicating that it is able to induce humoral and cell-mediated immune response in grouper. These results demonstrate that the HY9901Δhop mutant could be used as an effective live vaccine to combat V. alginolyticus in grouper.

Genomic structure, expression pattern and polymorphisms of GILT in golden pompano Trachinotus ovatus (Linnaeus 1758)

The interferon-g-inducible lysosomal thiol reductase (GILT) plays a significant character in the processing and presentation of MHC class II restricted antigen (Ag) by catalyzing disulfide bond reduction in mammals. To explore the function of GILT in the immune system of fish, we cloned a GILT gene homologue from Trachinotus ovatus, the full-length cDNA of GILT, which consisted of 2, 747 bp with a 771 bp open reading frame, encoding a protein of 256 amino acids. Moreover, similar to other species GILT gene, 7 exons and 6 introns were identified in T. ovatus, the deduced protein also possessed a representative characteristic of known GILT proteins. The result of real-time quantitative PCR showed that GILT mRNA was dramatically expressed in immune-associated tissues, such as spleen (p < 0.01) and kidney (p < 0.05). Bacterial challenge revealed that GILT mRNA level remarkably up-regulation in liver, spleen, kidney and intestine after induction with Photobacterium damsela. Furthermore, based on cloned sequences and genome BLAST, only one SNP site (ToGILT-S1-g.148C>G) was identified, and the allele C was significantly associated with high-susceptibility (HS) group, nevertheless, the allele G was dramatically associated with high-resistance (HR) group, indicating potential application for disease resistant breeding selection in T. ovatus.

Identification of SNP markers associated with tolerance to ammonia toxicity by selective genotyping from de novo assembled transcriptome in Litopenaeus vannamei

The high concentration of ammonia from deteriorated aquaculture environments and the intensive culture system could increase the susceptibility to pathogens and even cause high mortality in Litopenaeus vannamei. In addition, we have revealed that the ammonia-tolerant shrimp also have high disease resistance in L. vannamei. In the present study, in order to identify SNP markers associated with tolerance to ammonia toxicity, we developed and characterized SNPs from our previous transcriptome sequencing data of ammonia-stressed and control groups, and a marker-trait association analysis was performed for marker-assisted selection (MAS) to increase production in L. vannamei. A total of 318,919 SNPs were identified from the transcriptome sequences, and 25,772 SNPs were found from the 1826 ammonia-responsive genes with functional annotation. We selected 49 SNPs from 26 ammonia-responsive genes that had strong homologies to known genes in the shrimp and probably involved in immune function as candidate markers for genotyping, among which 39 SNPs were polymorphic for further marker-trait association analysis with the ammonia-tolerant (AT) and ammonia-sensitive (AS) groups. Finally, 12 out of the 49 SNP markers were identified to be associated with ammonia tolerance, containing 10 loci with significantly different allele frequencies and 10 loci with significantly different genotyping frequencies between the AT and AS groups. Among the associated markers, the G allele of TSP-1 (the first locus from the thrombospondin gene), the A allele of TSP-3, and the C allele of XBP1-5 (the fifth locus from X-box binding protein 1) only presented in the AT groups, but they were absent from the AS groups, which would be the preference of the MAS for the ammonia-tolerant shrimp. In addition, when the 12 associated SNP markers were used for analysis, the genetic diversity of the AT groups were significantly higher than that of the AS groups, but when the 39 loci were used there was no difference. This is the first report for the markers associated with ammonia tolerance in this species, indirectly with disease resistance, which provided important potential for genetic selection to increase survival rate and production in shrimp farming.

Transcriptome profiling provides gene resources for understanding gill immune responses in Japanese flounder (Paralichthys olivaceus) challenged with Edwardsiella tarda

Marine organisms are commonly under threats from various pathogens. Edwardsiella tarda is one of the fish pathogens that seriously infect cultured and wild fish species. Bacteremia caused by E. tarda can be a fatal disease in humans. Fish gill is a mucosa-associated lymphoid tissue that directly contacted with sea water. Generating gill transcriptomic resources that challenged by E. tarda is crucial for understanding the molecular mechanisms underlying gill immune responses. In this study, we performed transcriptome profiling of gene expression in Japanese flounder gills (Paralichthys olivaceus) challenged by E. tarda with different stress duration. An average of 40 million clean reads per library were obtained, of which approximately 83.2% were successfully mapped to the reference genome. 456 and 1037 differential expressed genes (DEGs) were identified at 8 h and 48 h post-injection, respectively. Gene annotation analysis and protein-protein interaction networks were conducted to obtain the key interaction relationships of immune-related DEGs during pathogens infection. 24 hub genes with multiple protein-protein interaction relationships or involved in multiple KEGG signaling pathways were discovered and validated by qRT-PCR. These hub genes mainly participated in Leukocyte transendothelial migration signaling pathway, B cell receptor signaling pathway, Wnt signaling pathway and Apoptosis signaling pathway. This study represents the first gill transcriptomic analysis based on protein-protein interaction networks in fish and provides valuable gene resources for understanding the fish gill immunity, which can pave the way to understand the molecular mechanisms of immune responses with E. tarda infection.

Transcriptome profiling based on protein-protein interaction networks provides a core set of genes for understanding blood immune response mechanisms against Edwardsiella tarda infection in Japanese flounder (Paralichthys olivaceus)

Marine organisms are commonly under threat from various pathogens. Edwardsiella tarda is one of the fish pathogens that can infect both cultured and wild fish species. E. tarda can also infect other vertebrates, including amphibians, reptiles, and mammals. Bacteremia caused by E. tarda can be a fatal disease in humans. Blood acts as a pipeline for the fish immune system. Generating blood transcriptomic resources from fish challenged by E. tarda is crucial for understanding molecular mechanisms underlying blood immune response process. In this study, we performed transcriptome-wide gene expression profiling of Japanese flounder (Paralichthys olivaceus) challenged by 8 and 48 h E. tarda stress. An average of 37 million clean reads per library was obtained, and approximately 85.6% of these reads were successfully mapped to the reference genome. In addition, 808 and 1265 differential expression genes (DEGs) were found at 8 and 48 h post-injection, respectively. Gene Ontology (GO) functional enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were conducted to search immune-related DEGs. A protein-protein interaction network was constructed to obtain the interaction relationship of immune genes during pathogens stress. Based on KEGG and protein association networks analysis, 30 hub genes were discovered and validated by quantitative RT-PCR. This study represents the first transcriptome analysis based on protein-protein interaction networks in fish and provides us with valuable gene resources for the research of fish blood immunity, which can significantly assist us to further understand the molecular mechanisms of humans and other vertebrates against E. tarda.

Outer membrane protein A: An immunogenic protein induces highly protective efficacy against Vibrio ichthyoenteri

Vibrio ichthyoenteri was an important causative agent of bacterial enteritis in flounder (Paralichthys olivaceus). Outer membrane protein A (OmpA) of Gram-negative pathogen was a major cell surface antigen. In the present study, OmpA of V. ichthyoenteri was recombinantly expressed in Escherichia coli, and the immunogenicity of OmpA was identified by western blotting using flounder anti-rOmpA and anti-V. ichthyoenteri antibodies. The vaccine potential of rOmpA was tested in a flounder model, and a high relative percentage of survival rate was obtained with 73.1% after challenge with V. ichthyoenteri. Meanwhile, the immune response of flounder induced by rOmpA was also investigated, and the results showed that the sIg + lymphocytes in blood, spleen, and pronephros significantly proliferated, and the peak levels occurred at week 4 after immunization. Moreover, rOmpA could induce higher levels of specific serum antibodies than the control group after immunization, and the peak level occurred at week 5 after immunization. Meanwhile, qRT-PCR analysis showed that the expressions of CD4-1, CD8α, IL-1β, IFN-γ, MHCIα and MHCIIα genes were significantly up-regulated after immunization with rOmpA. Taking together, these results demonstrated that rOmpA could evoke highly protective effects against V. ichthyoenteri challenge and induce strong immune response of flounder, which indicated that OmpA was a promising vaccine candidate.

Sequencing-based network analysis provides a core set of gene resource for understanding kidney immune response against Edwardsiella tarda infection in Japanese flounder

Marine organisms are under a frequent threat from various pathogens. Edwardsiella tarda is one of the major fish pathogens infecting both cultured and wild fish species. It can also infect a variety of other vertebrates, including amphibians, reptiles, and mammals, and bacteremia caused by E. tarda can be fatal in humans. The kidney is the largest lymphoid organ in fish, and generating kidney transcriptomic information under different stresses is crucial for understanding molecular mechanisms underlying the immune responses in the kidneys. In this study, we performed transcriptome-wide gene expression profiling of the Japanese flounder (Paralichthys olivaceus) challenged by 8 and 48 h of E. tarda infection. An average of 40 million clean reads per library was obtained, and approximately 81.6% of these reads were successfully mapped to the reference genome. In addition, 1319 and 4439 differentially expressed genes (DEGs) were found at 8 and 48 h post-injection, respectively. Gene Ontology (GO) functional enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were performed to search immune-related DEGs. A protein-protein interaction network was constructed to ascertain the relationship between interacting immune genes during pathogen-induced stress. Based on the KEGG and protein association networks analysis, 24 hub genes were discovered and validated by qRT-PCR. To our knowledge, this study is the first to represent the kidney transcriptome analysis based on protein-protein interaction networks in fish. Our results provide valuable gene resources for further research on kidney immune response in fish, which can significantly improve our understanding of the molecular mechanisms underlying the immune response to E. tarda in humans and other vertebrates.

Immunogenicity and efficacy of DNA vaccine encoding antigenic AcfA via addition of the molecular adjuvant Myd88 against Vibrio alginolyticus in Epinephelus coioides

DNA vaccines had been widely used against microbial infection in animals. The use of molecular adjuvants to improve the immunogenicity of DNA vaccines has been increasingly studied in recent years. MyD88 is one of the adapter molecules to activate the signaling cascades and produces inflammatory mediators, and its immunological role and adjuvant potential which had been proved in mammals were rarely reported in fish species. In this study, plasmid pcMyD88 was constructed and the capacity of MyD88 as molecular adjuvant was explored by co-injecting with a DNA vaccine encoding AcfA against Vibrio alginolyticus infection in orange spotted grouper. The results suggested that it needed at least 7 days to transported DNA vaccine pcacfA or molecular adjuvant pcMyD88 from the injected muscle to kidney and spleens and stimulate host's immune system for later protection. The co-injection of pcMyD88 with DNA vaccine pcacfA could increase significantly specific antibody levels and the expression levels of the immune-related genes including MHCIα, MHCIIα, CD4, CD8α, IL-1β and TNFα. Furthermore, pcMyD88 enhanced the immunoprotection of pcacfA against V. alginolyticus infection, with the significantly higher RPS of 83.3% in pcMyD88 + pcacfA group compared with that of pcacfA alone (73.3%) at challenging test of 10 weeks post vaccination. Together, these results clearly demonstrate that MyD88 is an effective adjuvant for the DNA vaccine pcacfA in orange spotted grouper.

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.

Studying the Genetics of Resistance to CyHV-3 Disease Using Introgression from Feral to Cultured Common Carp Strains

Sustainability and further development of aquaculture production are constantly challenged by outbreaks of fish diseases, which are difficult to prevent or control. Developing fish strains that are genetically resistant to a disease is a cost-effective and a sustainable solution to address this challenge. To do so, heritable genetic variation in disease resistance should be identified and combined together with other desirable production traits. Aquaculture of common carp has suffered substantial losses from the infectious disease caused by the cyprinid herpes virus type 3 (CyHV-3) virus and the global spread of outbreaks indicates that many cultured strains are susceptible. In this research, CyHV-3 resistance from the feral strain “Amur Sassan” was successfully introgressed into two susceptible cultured strains up to the first backcross (BC₁) generation. Variation in resistance of families from F₁ and BC₁ generations was significantly greater compared to that among families of any of the susceptible parental lines, a good starting point for a family selection program. Considerable additive genetic variation was found for CyHV-3 resistance. This phenotype was transferable between generations with contributions to resistance from both the resistant feral and the susceptible cultured strains. Reduced scale coverage (mirror phenotype) is desirable and common in cultured strains, but so far, cultured mirror carp strains were found to be susceptible. Here, using BC₁ families ranging from susceptible to resistant, no differences in resistance levels between fully scaled and mirror full-sib groups were found, indicating that CyHV-3 resistance was successfully combined with the desirable mirror phenotype. In addition, the CyHV-3 viral load in tissues throughout the infection of susceptible and resistant fish was followed. Although resistant fish get infected, viral loads in tissues of these fish are significantly lesser than in those of susceptible fish, allowing them to survive the disease. Taken together, in this study we have laid the foundation for breeding CyHV-3-resistant strains and started to address the mechanisms underlying the phenotypic differences in resistance to this disease.

Molecular characterization of MHC class IIB genes of sympatric Neotropical cichlids

BACKGROUND

The Major Histocompatibility Complex (MHC) is a key component of the adaptive immune system of all vertebrates and consists of the most polymorphic genes known to date. Due to this complexity, however, MHC remains to be characterized in many species including any Neotropical cichlid fish. Neotropical crater lake cichlids are ideal models to study evolutionary processes as they display one of the most convincing examples of sympatric and repeated parallel radiation events within and among isolated crater lakes.

RESULTS

Here, we characterized the genes of MHC class IIB chain of the Midas cichlid species complex (Amphilophus cf. citrinellus) including fish from five lakes in Nicaragua. We designed 19 new specific primers anchored in a stepwise fashion in order to detect all alleles present. We obtained 866 genomic DNA (gDNA) sequences from thirteen individuals and 756 additional sequences from complementary DNA (cDNA) of seven of those individuals. We identified 69 distinct alleles with up to 25 alleles per individual. We also found considerable intron length variation and mismatches of alleles detected in cDNA and gDNA suggesting that some loci have undergone pseudogenization. Lastly, we created a model of protein structure homology for each allele and identified their key structural components.

CONCLUSIONS

Overall, the Midas cichlid has one of the most diverse repertoires of MHC class IIB genes known, which could serve as a powerful tool to elucidate the process of divergent radiations, colonization and speciation in sympatry.

Characterization and functional analysis of a novel C1q-domain-containing protein in Japanese flounder (Paralichthys olivaceus)

The complement system is important in the innate immune response. C1q-domain-containing proteins have multiple functions and occur extensively in invertebrates and vertebrates. In this study, PoC1ql3 encoding a C1q-domain-containing protein in the Japanese flounder was identified. The 266-amino-acid polypeptide encoded, PoC1ql3, shares high sequence and structural similarity with orthologues in other fish and mammals. PoC1ql3 is abundantly expressed in the brain, but less in the blood, gills, and liver. Transcripts of PoC1ql3 were down-regulated in the spleen and liver 6-24 h after bacterial infection, but were significantly up-regulated after 48 h. Full-length PoC1ql3 (C1ql3-full) and its gC1q domain (C1ql3-part) were both exerted anti-Edwardsiella tarda activity. C1ql3-part bound to lipopolysaccharide and peptidoglycan, and exerted antibacterial effects against E. tarda in vivo, suggesting that C1ql3 functions as a pathogen-recognition receptor. Therefore, PoC1ql3 functions in the innate immune system, which would facilitate the investigation of the immune system in Japanese flounder.

Edwardsiella tarda Outer Membrane Protein C: An Immunogenic Protein Induces Highly Protective Effects in Flounder (Paralichthys olivaceus) against Edwardsiellosis

Outer membrane protein C of Edwardsiella tarda is a major cell surface antigen and it was identified to be an immunogenic protein by Western blot using flounder (Paralichthys olivaceus) anti-recombinant OmpC (rOmpC), and anti-E. tarda antibodies. rOmpC tested the immune protective effect against E. tarda challenge in a flounder model and produced a relative percentage of survival rate of 85%. The immune response of flounder induced by rOmpC was investigated, and the results showed that: (1) the levels of specific serum antibodies induced by rOmpC were significantly higher than the control group after the second week after immunization, and the peak level occurred at week five after immunization; (2) rOmpC could induce the proliferation of sIg+ lymphocytes, and the peak levels of sIg+ lymphocytes in blood, spleen, and pronephros occurred at 4–5 weeks after immunization; and (3) the MHCIIα, CD4-1, IL-1β, IL-6 and TNF-α genes were significantly induced after being injected with rOmpC. Taken together, these results demonstrated that rOmpC could evoke highly protective effects against E. tarda challenge and induce strong innate immune response and humoral immune response of flounder, which indicated that OmpC was a promising vaccine candidate against E. tarda infection.

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.

Characterization of 40 full-length MHC class IIA functional alleles in miiuy croaker: Polymorphism and positive selection

The major histocompatibility complex is a highly polymorphic gene superfamily in vertebrates that plays an important role in adaptive immune response. In the present study, we identified 40 full-length miiuy croaker MHC class IIA (Mimi-DAA) functional alleles from 26 miiuy croaker individuals and found that the alleles encode 30 amino acid sequences. A high level of polymorphism in Mimi-DAA was detected in miiuy croaker. The rate of non-synonymous substitutions (d(N)) occurred at a significantly higher frequency than that of synonymous substitutions (d(S)) in the peptide-binding region (PBR) and non-PBR. This result suggests that balancing selection maintains polymorphisms at the Mimi-DAA locus. Phylogenetic analysis based on the full-length sequences showed that the Mimi-DAA alleles clustered into three groups. However, the phylogenetic tree constructed using the exon 2 sequences indicated that the Mimi-DAA alleles clustered into two groups. A total of 22 positively selected sites were identified on the Mimi-DAA alleles after testing for positive selection, and five sites were predicted to be associated with the binding of peptide antigen, suggesting that a few selected residues may play a significant role in immune function.

Genome-wide SNP identification for the construction of a high-resolution genetic map of Japanese flounder (Paralichthys olivaceus): applications to QTL mapping of Vibrio anguillarum disease resistance and comparative genomic analysis

High-resolution genetic maps are essential for fine mapping of complex traits, genome assembly, and comparative genomic analysis. Single-nucleotide polymorphisms (SNPs) are the primary molecular markers used for genetic map construction. In this study, we identified 13,362 SNPs evenly distributed across the Japanese flounder (Paralichthys olivaceus) genome. Of these SNPs, 12,712 high-confidence SNPs were subjected to high-throughput genotyping and assigned to 24 consensus linkage groups (LGs). The total length of the genetic linkage map was 3,497.29 cM with an average distance of 0.47 cM between loci, thereby representing the densest genetic map currently reported for Japanese flounder. Nine positive quantitative trait loci (QTLs) forming two main clusters for Vibrio anguillarum disease resistance were detected. All QTLs could explain 5.1-8.38% of the total phenotypic variation. Synteny analysis of the QTL regions on the genome assembly revealed 12 immune-related genes, among them 4 genes strongly associated with V. anguillarum disease resistance. In addition, 246 genome assembly scaffolds with an average size of 21.79 Mb were anchored onto the LGs; these scaffolds, comprising 522.99 Mb, represented 95.78% of assembled genomic sequences. The mapped assembly scaffolds in Japanese flounder were used for genome synteny analyses against zebrafish (Danio rerio) and medaka (Oryzias latipes). Flounder and medaka were found to possess almost one-to-one synteny, whereas flounder and zebrafish exhibited a multi-syntenic correspondence. The newly developed high-resolution genetic map, which will facilitate QTL mapping, scaffold assembly, and genome synteny analysis of Japanese flounder, marks a milestone in the ongoing genome project for this species.

De Novo assembly of the Japanese flounder (Paralichthys olivaceus) spleen transcriptome to identify putative genes involved in immunity

Background

Japanese flounder (Paralichthys olivaceus) is an economically important marine fish in Asia and has suffered from disease outbreaks caused by various pathogens, which requires more information for immune relevant genes on genome background. However, genomic and transcriptomic data for Japanese flounder remain scarce, which limits studies on the immune system of this species. In this study, we characterized the Japanese flounder spleen transcriptome using an Illumina paired-end sequencing platform to identify putative genes involved in immunity.

Methodology/Principal Findings

A cDNA library from the spleen of P. olivaceus was constructed and randomly sequenced using an Illumina technique. The removal of low quality reads generated 12,196,968 trimmed reads, which assembled into 96,627 unigenes. A total of 21,391 unigenes (22.14%) were annotated in the NCBI Nr database, and only 1.1% of the BLASTx top-hits matched P. olivaceus protein sequences. Approximately 12,503 (58.45%) unigenes were categorized into three Gene Ontology groups, 19,547 (91.38%) were classified into 26 Cluster of Orthologous Groups, and 10,649 (49.78%) were assigned to six Kyoto Encyclopedia of Genes and Genomes pathways. Furthermore, 40,928 putative simple sequence repeats and 47, 362 putative single nucleotide polymorphisms were identified. Importantly, we identified 1,563 putative immune-associated unigenes that mapped to 15 immune signaling pathways.

Conclusions/Significance

The P. olivaceus transciptome data provides a rich source to discover and identify new genes, and the immune-relevant sequences identified here will facilitate our understanding of the mechanisms involved in the immune response. Furthermore, the plentiful potential SSRs and SNPs found in this study are important resources with respect to future development of a linkage map or marker assisted breeding programs for the flounder.

Genetic and genomic analyses for economically important traits and their applications in molecular breeding of cultured fish

The traits of cultured fish must continually be genetically improved to supply high-quality animal protein for human consumption. Economically important fish traits are controlled by multiple gene quantitative trait loci (QTL), most of which have minor effects, but a few genes may have major effects useful for molecular breeding. In this review, we chose relevant studies on some of the most intensively cultured fish and concisely summarize progress on identifying and verifying QTLs for such traits as growth, disease and stress resistance and sex in recent decades. The potential applications of these major-effect genes and their associated markers in marker-assisted selection and molecular breeding, as well as future research directions are also discussed. These genetic and genomic analyses will be valuable for elucidating the mechanisms modulating economically important traits and to establish more effective molecular breeding techniques in fish.

A genome scan for quantitative trait loci associated with Vibrio anguillarum infection resistance in Japanese flounder (Paralichthys olivaceus) by bulked segregant analysis

A recent genetic linkage map was employed to detect quantitative trait loci (QTLs) associated with Vibrio anguillarum resistance in Japanese flounder. An F1 family established and challenged with V. anguillarum in 2009 was used for QTL mapping. Of the 221 simple sequence repeat (SSR) markers used to detect polymorphisms in the parents of F1, 170 were confirmed to be polymorphic. The average distance between the markers was 10.6 cM. Equal amounts of genomic DNA from 15 fry that died early and from 15 survivors were pooled separately to constitute susceptible bulk and resistance bulk DNA. Bulked segregant analysis and QTL mapping were combined to detect candidate SSR markers and regions associated with the disease. A genome scan identified four polymorphic SSR markers, two of which were significantly different between susceptible and resistance bulk (P=0.008). These two markers were located in linkage group (LG) 7; therefore, all the SSR markers in LG7 were genotyped in all the challenged fry by single marker analysis. Using two different models, 11-17 SSR markers were detected with different levels of significance. To confirm the associations of these markers with the disease, composite interval mapping was employed to genotype all the challenged individuals. One and three QTLs, which explained more than 60 % of the phenotypic variance, were detected by the two models. Two of the QTLs were located at 48.6 cM. The common QTL may therefore be a major candidate region for disease resistance against V. anguillarum infection.

Genetic polymorphism of major histocompatibility complex class IIB alleles and pathogen resistance in the giant spiny frog Quasipaa spinosa

Major histocompatibility complex (MHC) genes are candidates for determining disease susceptibility due to their pivotal role in both innate and adaptive immune responses. Accordingly, the association between the genetic variation of MHC genes and the pathogen resistance has been investigated in numerous vertebrates. To date, however, little is reported in amphibians. In this study, we investigate the genetic variation at the MHC class IIB gene in the giant spiny frog Quasipaa spinosa, which has high commercial value in China. The full length of MHC class IIB cDNA was cloned from Q. spinosa by homology cloning and rapid amplification of cDNA end-polymerase chain reaction (RACE-PCR). Two MHC class IIB loci were identified in Q. spinosa. We also developed PCR primers for a portion of the second exon of the MHC class IIB gene. A total of 26 MHC class IIB alleles were identified. The dN rate was significantly higher than the dS rate in the putative peptide-binding region, thereby proving the positive selection hypothesis. In addition, individuals intraperitoneally injected with Aeromonas hydrophila were used to study the association between MHC class IIB alleles and pathogen resistance/susceptibility, to explore the specific alleles in balancing selection. Eighty frogs were used after exposure to A. hydrophila infection. Nine alleles were used to study the association between the alleles and disease resistance. Two alleles, namely, Pasa-DAB(∗)1301 and Pasa-DAB(∗)0901, were significantly associated with resistance against A. hydrophila. This study provides valuable information on the structure of the MHC class IIB gene and confirms the association between MHC class IIB gene alleles and disease resistance to bacterial infection in Q. spinosa. Moreover, pathogen resistance-related MHC markers can be used for the selective breeding of the giant spiny frog.