Molecular cloning, genomic structure, polymorphism and expression analysis of major histocompatibility complex class IIA and IIB genes of half-smooth tongue sole (Cynoglossus semilaevis)

High water temperature-mediated immune gene expression of olive flounder, Paralichthys olivaceus according to pre-stimulation at high temperatures

Increased ocean temperature due to global warming affects the health and immunity of fish. In this study, juvenile Paralichthys olivaceus were exposed to high temperature after pre-heat (Acute: Acute heat shock at 32 °C, AH-S: Acquired heat shock at 28 °C & short recovery (2h) and heat shock at 32 °C, AH-L: acquired heat shock at 28 °C and long recovery (2 days), AH-LS: acquired heat shock at 28 °C & long (2 days) + short (2h) recovery). Heat shock after pre-heat significantly upregulated various immune-related genes, including interleukin 8 (IL-8), c-type lysozyme (c-lys), immunoglobulin M (IgM), Toll-like receptor 3 (tlr3), major histocompatibility complex IIα (mhcIIα) and cluster of differentiation 8α (cd8α) in the liver and brain of P. olivaceus. This study showed pre-exposure to high temperatures below the critical temperature can activate fish immunity and increase tolerance to high temperatures.

Comprehensive transcriptomics and proteomics analysis of Carassius auratus gills in response to Aeromonas hydrophila

As one of the mucosal barriers, fish gills represent the first line of defense against pathogen infection. However, the exact mechanism of gill mucosal immune response to bacterial infection still needs further investigation in fish. Here, to investigate pathological changes and molecular mechanisms of the mucosal immune response in the gills of crucian carp (Carassius auratus) challenged by Aeromonas hydrophila, the transcriptomics and proteomics were performed by using multi-omics analyses of RNA-seq coupled with iTRAQ techniques. The results demonstrated gill immune response were mostly related to the activation of complement and coagulation cascades, antigen processing and presentation, phagosome, NOD-like receptor (NLR) and nuclear factor κB (NFκB) signaling pathway. Selected 21 immune-related DEGs (ie., Clam, nfyal, snrpf, acin1b, psme, sf3b5, rbm8a, rbm25, prpf18, g3bp2, snrpd3l, tecrem-2, cfl-A, C7, lysC, ddx5, hsp90, α-2M, C9, C3 and slc4a1a) were verified for their immune roles in the A. hydrophila infection via using qRT-PCR assay. Meanwhile, some complement (C3, C7, C9, CFD, DF and FH) and antigen presenting (HSP90, MHC Ⅱ, CALR, CANX and PSME) proteins were significantly participated in the process of defense against infections in gill tissues, and protein-protein interaction (PPI) network displayed the immune signaling pathways and interactions among these DEPs. The correlation analysis indicated that the iTRAQ and qRT-PCR results was significantly correlated (Pearson's correlation coefficient = 0.70, p < 0.01). To our knowledge, the transcriptomics and proteomics of gills firstly identified by multi-omics analyses contribute to understanding on the molecular mechanisms of local mucosal immunity in cyprinid species.

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.

Effects of autochthonous strains mixture on gut microbiota and metabolic profile in cobia (Rachycentron canadum)

The fish immune system is a topic or subject that offers a unique understanding of defensive system evolution in vertebrate heredity. While gut microbiota plays several roles in fish: well-being, promoting health and growth, resistance to bacterial invasion, regulation of energy absorption, and lipid metabolism. However, studies on fish gut microbiota face practical challenges due to the large number of fish varieties, fluctuating environmental conditions, and differences in feeding habits. This study was carried out to evaluate the impacts of supplemented three autochthonous strains, Bacillus sp. RCS1, Pantoea agglomerans RCS2, and Bacillus cereus RCS3 mixture diet on cobia fish (Rachycentron canadum). Also, chromatography, mass spectrometry and high throughput sequencing were combined to explore composition and metabolite profile of gut microbiota in juvenile cobia fed with supplemented diet. In the trial group, juvenile cobia received diets supplemented with 1 × 10¹² CFU mL^-1 autochthonous strains for ten weeks and a control diet without supplementation. Juvenile cobia receiving diets supplementation exhibited significantly improved growth than those without additives (control). Haematological indices, such as red blood cells, white blood cells, corpuscular haemoglobin concentration, mean corpuscular volume, haemoglobin, and mean corpuscular haemoglobin, were higher in the supplemented group. Similarly, digestive enzymes (trypsin, lipase, amylase, pepsin and cellulose, activities) activities were higher in supplemented diet with an indigenous isolates mixture. Serum biochemical parameters albumin, globulin, and total protein were significantly higher, while triglyceride, alanine aminotransferase, aspartate aminotransferase, alkaline phosphatase, and cholesterol showed no significant difference. On the other hand, glucose was significantly (P < 0.05) higher in the group without supplementation. On gene expression in the midgut, Immunoglobulin, Colony-stimulating factor receptor 1, major histocompatibility complex 1 were up-regulated by native isolates while T cell receptor beta, and Major histocompatibility complex 2 showed no significant difference. Gut bacterial composition was altered in fish receiving supplemented diet with autochthonous strains. Metabolomics also revealed that some metabolic pathways were considerably enriched in fish fed with supplemented diet; pathway analysis based on Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment revealed that differentially expressed metabolites were involved in galactose metabolism, tryptophan metabolism, carbohydrate digestion and absorption, purine metabolism, and ABC transporters. Functional analysis of bacterial community showed that differences in enriched metabolic pathways generally comprised carbohydrate and its metabolites, nucleotide and its metabolites, amino acid and its metabolites, heterocyclic compounds, and tryptamines, cholines, pigments. The current investigation results showed that autochthonous strains mixture has significantly enhanced the growth, survival, and innate and adaptive immunities of juvenile cobia.

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

Background

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

Results

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

Conclusion

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

Molecular cloning and immune characterization of CIITA in artificially challenged Chinese tongue sole (Cynoglossus semilaevis) with Vibrio harveyi

In mammals, Class II, major histocompatibility complex (MHC II) transactivator (CIITA) recognizes microbial pathogens and triggers immune responses. In Chinese tongue sole Cynoglossus semilaevis, Cs-CIITA was prevalently expressed in various tissues. Cs-CIITA, Cs-MHC IIA and Cs-MHC IIB were expressed significantly higher in skin in susceptible families infected with Vibrio harveyi, while higher expression of Cs-CIITA and Cs-MHC IIB was examined in liver in resistant families. In addition, the three genes were up-regulated in gill, skin, intestine, liver, spleen and kidney at 48 h or 72 h after V. harveyi infection. Furthermore, the three genes were co-expressed in the epithelial mucous cells of gill, skin, and intestine. Knockdown of Cs-CIITA regulates the expression of other inflammation-related genes, including CD40, IL-1β, IL-8, RelB, NFκB, and Myd88. These results suggest that CIITA functions in the inflammatory responses of C. semilaevis against V. harveyi, via MHC II transcriptional regulation.

Role of major histocompatibility complex II antigen-presentation pathway genes in orange-spotted grouper infected with Cryptocaryon irritans

Cryptocaryon irritans, a pathogen model for fish mucosal immunity, causes skin mucosal and systematic humoral immune response. Where and how MHC II antigen presentation occurs in fish infected with C. irritans remain unknown. In this study, the full-length cDNA of the grouper cysteine protease CTSS was cloned. The expression distributions of six genes (CTSB, CTSL, CTSS, GILT, MHC IIA and MHC IIB) involved in MHC II antigen presentation pathway were tested. These genes were highly expressed in systematic immune tissues and skin and gill mucosal-associated immune tissues. All six genes were upregulated in skin at most time points. Five genes expected CTSS was upregulated in spleen at most time points. CTSB, CTSL and MHC IIA were upregulated in the gill and head kidney at some time points. These results indicate that the presentation of MHC II antigen intensively occurred in local infected skin and gill. Spleen, not head kidney, had the most extensive systematic antigen presentation. In skin, six genes most likely peaked at day 2, earlier than in spleen (5-7 days), marking an earlier skin antibody peak than any recorded in serum previously. This significant and earlier mucosal antigen presentation indicates that specific immune response occurs in local mucosal tissues.

Identification and application of piwi-interacting RNAs from seminal plasma exosomes in Cynoglossus semilaevis

BACKGROUND

Piwi-interacting RNAs (piRNAs) have been linked to epigenetic and post-transcriptional gene silencing of retrotransposons in germ line cells, particularly in spermatogenesis. Exosomes are important mediators of vesicle transport, and the piRNAs in exosomes might play an important role in cell communication and signal pathway regulation. Moreover, exosomic piRNAs are promising biomarkers for disease diagnosis and physiological status indication. We used Cynoglossus semilaevis because of its commercial value and its sexual dimorphism, particularly the sex reversed "pseudomales" who have a female karyotype, produce sperm, and copulate with normal females to produce viable offspring.

RESULTS

To determine whether piRNAs from fish germ line cells have similar features, seminal plasma exosomes from half-smooth tongue sole, C. semilaevis, were identified, and their small RNAs were sequenced and analysed. We identified six signature piRNAs as biomarkers in exosomes of seminal plasma from males and pseudomale C. semilaevis. Bioinformatic analysis showed that all six signatures were sex-related, and four were DNA methylation-related and transposition-related piRNAs. Their expression profiles were verified using real-time quantitative PCR. The expression of the signature piRNAs was markedly higher in males than in pseudomales. The signature piRNAs could be exploited as male-specific biomarkers in this fish.

CONCLUSIONS

These signatures provide an effective tool to explore the regulatory mechanism of sex development in C. semilaevis and may provide guidance for future research on the function of piRNAs in the generative mechanism of sex reversed "pseudomales" in C. semilaevis.

Expression of immune genes in Indian major carp, Catla catla challenged with Flavobacterium columnare

Columnaris disease, caused by Flavobacterium columnare, is one of the important bacterial diseases responsible for large-scale mortalities in numerous freshwater fishes globally. This disease can cause up to 100% mortality within 24 h of infection and is considered to be a cause of concern for aquaculture industry. Despite being a serious disease, scarce information is available regarding host-pathogen interaction, particularly the modulation of different immune genes in response to F. columnare infection. Therefore, in the present study, an attempt has been made to study expression of important immune regulatory genes, namely IL-1β, iNOS, INF-γ, IL-10, TGF-β, C3, MHC-I and MHC-II in gills and kidney of Catla catla following experimental infection with F. columnare. The expression analysis of immune genes revealed that transcript levels of IL-1β, iNOS, IL-10, TGF-β, C3 and MHC-I were significantly up-regulated (p < 0.05) in both the organs of the infected catla. IFN-γ and MHC-II were up-regulated in gills of infected catla whereas, both the genes showed down-regulation in kidney. The results indicate that important immune genes of C. catla are modulated following infection with F. columnare. The knowledge thus generated will strengthen the understanding of molecular pathogenesis of F. columnare in Indian major carp C. catla.

Transcriptome analysis in the skin of Carassius auratus challenged with Aeromonas hydrophila

Skin plays an important role in the innate immune responses of fish, particularly towards bacterial infection. To understand the molecular mechanism of mucosal immunity of fish during bacterial challenge, a de novo transcriptome assembly of crucian carp Carassius auratus skin upon Aeromonas hydrophila infection was performed, the latter with Illumina Hiseq 2000 platform. A total of 118111 unigenes were generated and of these, 9693 and 8580 genes were differentially expressed at 6 and 12 h post-infection, respectively. The validity of the transcriptome results of eleven representative genes was verified by quantitative real-time PCR (qRT-PCR) analysis. A comparison with the transcriptome profiling of zebrafish skin to A. hydrophila with regards to the mucosal immune responses revealed similarities in the complement system, chemokines, heat shock proteins and the acute-phase response. GO and KEGG enrichment pathway analyses displayed the significant immune responses included TLR, MAPK, JAK-STAT, phagosome and three infection-related pathways (ie., Salmonella, Vibrio cholerae and pathogenic Escherichia coli) in skin. To our knowledge, this study is the first to describe the transcriptome analysis of C. auratus skin during A. hydrophila infection. The outcome of this study contributed to the understanding of the mucosal defense mechanisms in cyprinid species.

Analysis of apolipoprotein multigene family in spotted sea bass (Lateolabrax maculatus) and their expression profiles in response to Vibrio harveyi infection

Apolipoproteins (Apos), which are the protein components of plasma lipoproteins, play important roles in lipid transport in vertebrates. It has been demonstrated that in teleosts, several Apos display antimicrobial activity and play crucial roles in innate immunity. Despite their importance, apo genes have not been systematically characterized in many aquaculture fish species. In our study, a complete set of 23 apo genes was identified and annotated from spotted sea bass (Lateolabrax maculatus). Phylogenetic and homology analyses provided evidence for their annotation and evolutionary relationships. To investigate their potential roles in the immune response, the expression patterns of 23 apo genes were determined in the liver and intestine by qRT-PCR after Vibrio harveyi infection. After infection, a total of 20 differentially expressed apo genes were observed, and their expression profiles varied among the genes and tissues. 5 apo genes (apoA1, apoA4a.1, apoC2, apoF and apoO) were dramatically induced or suppressed (log2 fold change >4, P < 0.05), suggesting their involvement in the immune response of spotted sea bass. Our study provides a valuable foundation for future studies aimed at uncovering the specific roles of each apo gene during bacterial infection in spotted sea bass and other teleost species.

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.

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.

Transcriptome profiling and expression analysis of immune responsive genes in the liver of Golden mahseer (Tor putitora) challenged with Aeromonas hydrophila

Transcriptome profiling has been used to decipher the novel mechanisms behind immune responses of the fishes. However, the molecular mechanism underlining immune response in mahseer is not studied so far. Fishes are greatly affected by bacterial pathogens such as Aeromonas hydrophila. In this study, transcriptome response of golden mahseer (Tor putitora) infected with A. hydrophila was examined using paired end Illumina sequencing of liver tissue to understand the immune response of the fish. The de novo assembly generated 61,042 unigenes ranging from 200 to 9322 bp in length and an average length of 463 bp. The gene ontology annotations resulted a total of 131,826 term assignments to the annotated transcriptome including 60,846 (46.16%) allocations from the biological process; 21,603 (16.39%) from molecular function and 49,377 (37.46%) from cellular components. Differential gene expression analysis of the transcriptome data from challenged and control group revealed 1104 upregulated and 1304 down-regulated unigenes. The differentially expressed genes were mainly involved in the pathways including cell surface receptor signaling, TH1 and TH2 cell differentiation, pathogen recognition, and immune system process/defense response especially complement cascade. Twelve unigenes including ankyrin, serum amyloid, hsp4b, STAT3, complement factor c3 and c7 were validated using qPCR and found differentially expressed in accordance with in silico expression analysis. The results obtained in this study will provide the first and crucial information on the molecular mechanism of mahseer fishes against bacterial infection.

The influence of concentration of inactivated Edwardsiella tarda bacterin and immersion time on antigen uptake and expression of immune-related genes in Japanese flounder (Paralichthys olivaceus)

Our previous work has demonstrated that the immune response of Japanese flounder was associated with the concentration of formalin-inactivated Edwardsiella tarda and immersion time. In order to further investigate the influence of immersion vaccine dose and bath time on the antigen uptake, formalin-killed Edwardsiella tarda bacterin was prepared and adjusted to four concentrations (10⁹, 10⁸, 10⁷, 10⁶ cfu ml^-1) for 30, 60 and 90 min immersion in Japanese flounder model, respectively. Absolute quantitative real-time PCR was employed to examine the bacterin uptake in gill, skin, spleen and kidney at 3 and 6 h post vaccination. The results showed that the antigen uptaken in gills and skin were significant higher than spleen and kidney, and the antigen amounts in gill and skin both declined from 3 to 6 h, whereas the antigen amounts in spleen and kidney gradually increased. Significant higher antigen amounts were detected in 10⁹-30, 10⁹-60, 10⁸-60, 10⁸-90 and 10⁸-90 groups than other groups (P < 0.05), especially the 10⁸-60min group displayed the highest antigen uptaken. Meanwhile, the expression profiles of antigen recognization and presentation genes (MHCⅡα, TcRα, CD4-1), immunoglobulins (IgM, IgT), inflammatory cytokines (IL-1β, IL-6), heat shock protein 70 (HSP70) and c-type lysozyme were analyzed using real-time PCR. On the whole, the transcription levels of the eight genes exhibited to be higher in 10⁷-90, 10⁸ and 10⁹ cfu ml^-1 groups than other groups (P < 0.05), especially the 10⁸-60 group displayed the highest up-regulation. These results demonstrated that immersion with formalin-inactivated E. tarda, especially under 10⁸-60 min condition could efficiently enhance the antigen uptake and the expression of immune-related genes, which provided evidences for an enhanced vaccination effects under an optimized combination of vaccine dose and immersion time.

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.

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.

Complement factor I from flatfish half-smooth tongue (Cynoglossus semilaevis) exhibited anti-microbial activities

Complement factor I (Cfi) is a soluble serine protease which plays a crucial role in the modulation of complement cascades. In the presence of substrate modulating cofactors (such as complement factor H, C4bp, CR1, etc), Cfi cleaves and inactivates C3b and C4b, thereby controlling the complement-mediated processes. In this study, we sequenced and characterized Cfi gene from Cynoglossus Semilaevis (designated as CsCfi) for the first time. The full-length cDNA of CsCfi was 2230 bp in length, including a 98 bp 5'-untranslated region (UTR), a 164 bp 3'-UTR and a 1968 bp open reading frame (ORF). It encoded a polypeptide of 656 amino acids, with a molecular mass of 72.28 kDa and an isoelectric point of 7.71. A signal peptide was defined at N-terminus, resulting in a 626-residue mature protein. Multiple sequence alignment revealed that Cfi proteins were well conserved with the typical modular architecture and identical active sites throughout the vertebrates, which suggested the conserved function of Cfi. Phylogenetic analysis indicated that CsCfi and the homologous Cfi sequences from teleosts clustered into a clade, separating from another clade from the cartilaginous fish and other vertebrates. Tissue expression profile analysis by quantitative real-time PCR (qRT-PCR) showed that CsCfi mRNA constitutively expressed in all tested tissues, with the predominant expression in liver and the lowest in stomach. Temporal expression levels of CsCfi after challenging with Vibrio anguillarum showed different expression patterns in intestine, spleen, skin, blood, head kidney and liver. The recombinant CsCfi (rCsCfi) protein showed broad-spectrum antimicrobial activities against the Gram-positive bacteria Staphylococcus aureus and the Gram-negative bacteria Escherichia coli, Pseudomonas aeruginosa and Shewanella putrefaciens. The research revealed that CsCfi plays an important role in C. Semilaevis immunity.

Transcriptome analysis revealed changes of multiple genes involved in immunity in Cynoglossus semilaevis during Vibrio anguillarum infection

Half-smooth tongue sole (Cynoglossus semilaevis) is one of the most valuable marine aquatic species in Northern China. Given to the rapid development of aquaculture industry, the C. semilaevis was subjected to disease-causing bacteria Vibrio anguillarum. It therefore is indispensable and urgent to understand the mechanism of C. semilaevis host defense against V. anguillarum infection. In the present study, the extensively analysis at the transcriptome level for V. Anguillarum disease in tongue sole was carried out. In total, 94,716 high quality contigs were generated from 75,884,572 clean reads in three libraries (HOSG, NOSG, and CG). 22,746 unigenes were identified when compared with SwissProt, an NR protein database and NT nucleotide database. 954 genes exhibiting the differentially expression at least one pair of comparison in all three libraries were identified. GO enrichment for these genes revealed gene response to biotic stimulus, immune system regulation, and immune response and cytokine production. Further, the pathways such as complement and coagulation cascades and Vibrio cholerae infection pathways were enriched in defensing of pathogen. Besides, 13,428 SSRs and 118,239 SNPs were detected in tongue sole, providing further support for genetic variation and marker-assisted selection in future. In summary, this study identifies several putative immune pathways and candidate genes deserving further investigation in the context of development of therapeutic regimens and lays the foundation for selecting resistant lines of C. semilaevis against V. anguillarum.

sghC1q, a novel C1q family member from half-smooth tongue sole (Cynoglossus semilaevis): identification, expression and analysis of antibacterial and antiviral activities

The C1q family includes many proteins that contain a globular (gC1q) domain, and this family is widely conserved from bacteria to mammals. The family is divided into three subgroups: C1q, C1q-like and ghC1q. In this study, a novel C1q family member, sghC1q, was cloned and identified from Cynoglossus semilaevis (named CssghC1q). The full-length CssghC1q cDNA spans 905 bp, including an open reading frame (ORF) of 768 bp, a 5'-untranslated region (UTR) of 25 bp and a 3'-UTR of 112 bp. The ORF encodes a putative protein of 255 amino acids (aa) with a deduced molecular weight of 28 kDa. The predicted protein contains a signal peptide (aa 1-19), a coiled-coil region (aa 61-102) and a globular C1q (gC1q) domain (aa 117-255). Protein sequence alignment indicated that the C-terminus of CssghC1q is highly conserved across several species. Phylogenetic analysis indicated that CssghC1q is most closely related to Maylandia zebra C1q-like-2-like. The CssghC1q genomic sequence spanned 1562 bp, with three exons and two introns. CssghC1q is constitutively expressed in all evaluated tissues, with the highest expression in the liver and the weakest in the heart. After a challenge with Vibrio anguillarum, CssghC1q transcript levels exhibited distinct time-dependent response patterns in the blood, head kidney, skin, spleen, intestine and liver. Recombinant CssghC1q protein exhibited antimicrobial activities against Gram-negative bacteria, Gram-positive bacteria and viruses. The minimum inhibitory concentration (MIC) values against Vibrio harveyi, Vibrio anguillarum, Pseudomonas aeruginosa and Staphylococcus aureus were 0.043 mg/mL, 0.087 mg/mL, 0.174 mg/mL and 0.025 mg/mL, respectively. A low concentration (0.06 mg/mL) of CssghC1q showed significant antiviral activity in vitro against nervous necrosis virus (NNV). These results suggest that CssghC1q plays a vital role in immune defense against bacteria and viruses.

Transcriptome analysis of the Trachinotus ovatus: identification of reproduction, growth and immune-related genes and microsatellite markers

Background

The Trachinotus ovatus (Teleostei, Carangidae) is an economically important marine fish species in the world. However, the lack of genomic information regarding this species limits our understanding of the genetics and biological mechanisms in Trachinotus ovatus. In this study, high throughput transcriptome sequencing was used to obtain comprehensive genomic information in Trachinotus ovatus.

Principal Findings

Transcriptome sequencing was performed by using Illumina paired-end sequencing technology. The 98,534,862 high quality reads were yielded, and were de novo assembled into 156,094 unigenes with an average sequence length of 1179 bp. Transcriptome annotation revealed that 75,586 and 67,923 unigenes were functionally annotated in the NCBI non-redundant database and Swiss-Prot protein database, respectively. Functional analysis demonstrated that 67,923 unigenes were grouped into 25 Cluster of Orthologous Groups (COG) functional categories, 37,976 unigenes were clustered into 61 Gene Ontology (GO) terms, and 38,172 unigenes were assigned to 275 different Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Based on the transcriptome dataset, a large number of unigenes associated with reproduction, growth and immunity were identified. Furthermore, a total number of 38,794 simple sequence repeats (SSRs) were discovered and 16 polymorphic loci were characterized in Trachinotus ovatus.

Conclusion/Significance

The present study is the first transcriptome analysis of a fish species belonging to the genus Trachinotus and provides a valuable genomic resource for novel gene discovery, gene expression and regulation studies, and the identification of genetic markers in Trachinotus ovatus and the other fish of the genus Trachinotus.

Molecular cloning and expression analysis of major histocompatibility complex class I, IIA and IIB genes of blunt snout bream (Megalobrama amblycephala)

Major histocompatibility complex (MHC) plays an important role in the immune response of vertebrates. In this study, we isolated MHC class IIA and IIB genes from blunt snout bream (Megalobrama amblycephala) by rapid amplification of cDNA ends polymerase chain reaction (RACE-PCR). In order to study the function of the MHC genes in M. amblycephala, tissue distribution and immune response of the MHC genes to bacterial challenge were analyzed. All the characteristic features of MHC class II chain structure could be identified in the deduced amino sequences of MHC IIA and IIB, including the leader peptide, α1/β1 and α2/β2 domains, connecting peptide and transmembrane and cytoplasmic regions, as well as conserved cysteines and N-glycosylation site. The deduced amino acid sequence of the MHC IIA and IIB molecules shared from 48% to 88% and from 65% to 77% similarity with those of other teleosts, respectively. Quantitative real-time PCR (qRT-PCR) demonstrated that MHC I and II genes were ubiquitously expressed in ten tissues, with high level in immune related tissues, including kidney, intestine, gill and spleen. Challenge of M. amblycephala with the extracellular pathogen, Aeromonas hydrophila, resulted in a significant increase in the expression of MHC I, MHC IIA and IIB mRNA within 72 h after infection in gill, kidney, intestine and liver, followed by a recovery to normal level after 120 h. The changes of expression levels for MHC IIA and IIB in most tissues were significantly higher than that of MHC I in the corresponding tissues at most time points (P<0.05). These results demonstrated the MHC genes played an important role in response to bacterial infection in M. amblycephala; however, MHC class I and II genes showed different functional activity, which need be further investigated in teleost.

Molecular cloning, genomic structure, polymorphism and expression analysis of major histocompatibility complex class IIA gene of swamp eel Monopterus albus

Major histocompatibility complex (MHC) class II molecules play an important role in the immune response of vertebrates. In this paper, full-length MHC IIA cDNA was isolated from swamp eel (Monopterus albus) by rapid amplification of cDNA ends PCR. The genomic structure, molecular polymorphism, tissue distribution, and immune response of the MHC IIA gene to bacterial challenge were investigated. The full-length cDNA (GenBank accession No.: KC616308) is 1,509 bp in length including an 83 bp-long 5' untranslated region (UTR) and a 709 bp-long 3' UTR, which encoded a 238 amino acids protein. In the 2,339 bp-long MHC IIA genomic DNA, four exons and three introns were identified. Sequence comparison exhibited that the deduced amino acid sequence shared 27.1-66.3% identity with those of other species. Seven alleles were identified from five healthy individuals. Number of alleles per individual diversified from two to five. Five different 5' UTR sequences and two different 3' UTR sequences from one individual may infer the existence of five loci at least. Real-time quantitative PCR demonstrated that swamp eel MHC IIA transcripts were ubiquitously expressed in ten tissues, but the expression level was distinctly different. Significant changes were observed in liver, spleen, kidney and intestine after challenged with pathogenic bacteria Aeromonas hydrophilia.

MH II-DAB gene expression in grass carp Ctenopharyngodon idella (Valenciennes) after infection with the ciliate parasite, Ichthyophthirius multifiliis

The grass carp, Ctenopharyngodon idella (Valenciennes), is one of the most extensively aquacultured freshwater fish in China. However, because of the lack of effective control measures and the high-density culture environment, considerable economic losses are caused by infection of C. idella with the parasitic ciliate, Ichthyophthirius multifiliis. The major histocompatibility (MH) DAB gene belongs to antigen-presented genes in the class II genomic region, which is associated with parasite resistance. To understand the relationship of the DAB gene with I. multifiliis infection in grass carp, the expression profiles of MH II-DAB were studied in tissues using real-time quantitative polymerase chain reaction. The results showed that expression of the MH II-DAB gene was up-regulated in head kidney after I. multifiliis infection, and the expression peak appeared earlier in the study (case) group than in the control group. The obvious up-regulation peak of MH II-DAB gene was found at days 2 and 4 in skin; at 12 h to day 4 in spleen; at 12 h and days 1 and 6 in gill; and at day 10 in blood, whereas the MH II-DAB gene was down-regulated in liver and intestines after I. multifiliis infection. These results have implications for better understanding C. idella resistance to I. multifiliis infection.

Molecular cloning and expression analysis of major histocompatibility complex class IIB gene of the Whitespotted bambooshark (Chiloscyllium plagiosum)

Major histocompatibility complex (MHC) plays an important role in the immune response to antigenic peptides in vertebrates. In this study, the full length of MHC IIB cDNA was isolated from the Whitespotted bambooshark (Chiloscyllium plagiosum) by homology cloning, and the rapid amplification of cDNA ends polymerase chain reaction. As a result, the MHC IIB cDNA is 1,407 bp, which contains an open reading frame (ORF) of 831 bp encoding a protein of 276 amino acids. Furthermore, seven alleles of the complete MHC IIB ORF were detected and the variable sites were mainly located in the immunoglobulin-like (β2) region. Tissue distribution analysis showed that MHC IIB can be detected in all the ten tissues examined, with the highest expression in the spleen and gill. Challenge of C. plagiosum with the pathogenic bacteria, Vibrio harveyi, resulted in significant changes in the expression of MHC IIB mRNA in the three immune-related tissues (gill, liver and spleen). These results show that the MHC IIB plays an important role in response to bacterial infection in elasmobranches.

Major histocompatibility complex class IIA and IIB genes of Nile tilapia Oreochromis niloticus: genomic structure, molecular polymorphism and expression patterns

Major histocompatibility complex (MHC) is a large genomic region characterized by extremely high polymorphism, and it plays an important role in the immune response of vertebrates. In the present study, we isolated MHC class II genes from Nile tilapia in order to investigate the immune mechanism in tilapia and develop better strategies for disease prevention. Moreover, we cloned the full-length cDNA sequences of MHC IIA and IIB from Nile tilapia by the RACE approach. In addition, the genomic structure, molecular polymorphism and expression patterns of MHC II genes in Nile tilapia were also examined. Compared with that of other teleosts, Nile tilapia MHC class IIA contained four exons and three introns. The deduced amino acid sequence of the MHC IIA molecule shared 25.4-64.5% similarity with those of other teleosts and mammals. Six exons and five introns were identified from Nile tilapia MHC IIB, and the deduced amino acid sequence shared 26.9-74.7% similarity with those of other teleosts and mammals. All the characteristic features of MHC class II chain structure could be identified in the deduced sequences of MHC IIA and IIB molecules, including the leader peptide, α1/β1 and α2/β2 domains, connecting peptide and transmembrane and cytoplasmic regions, as well as conserved cysteines and N-glycosylation site. A total of 12 MHC IIA alleles were identified from six individuals. Four alleles originating from a single individual suggested that at least four MHC IIA loci existed. Moreover, 10 MHC IIB alleles were identified, among which four were detected in a single individual, suggesting that at least four MHC IIB loci existed. The expression of MHC IIA and IIB at the mRNA level in 10 types of normal tissues was determined using quantitative real-time PCR analysis. The highest expression level was detected in stomach and gill, whereas the lowest expression was detected in muscle and brain. Furthermore, MHC IIA and IIB were probably two candidate immune molecules involved in the resistance against streptococcosis, because their expression was significantly up-regulated in gill, kidney, intestine and spleen after the intraperitoneal injection of Streptococcus agalactiae.

Advances in research of fish immune-relevant genes: a comparative overview of innate and adaptive immunity in teleosts

Fish is considered to be an important model in comparative immunology studies because it is a representative population of lower vertebrates serving as an essential link to early vertebrate evolution. Fish immune-relevant genes have received considerable attention due to its role in improving understanding of both fish immunology and the evolution of immune systems. In this review, we discuss the current understanding of teleost immune-relevant genes for both innate and adaptive immunity, including pattern recognition receptors, antimicrobial peptides, complement molecules, lectins, interferons and signaling factors, inflammatory cytokines, chemokines, adaptive immunity relevant cytokines and negative regulators, major histocompatibility complexes, immunoglobulins, and costimulatory molecules. The implications of these factors on the evolutionary history of immune systems were discussed and a perspective outline of innate and adaptive immunity of teleost fish was described. This review may provide clues on the evolution of the essential defense system in vertebrates.

Identification and expression analysis of major histocompatibility complex IIB gene in orange-spotted grouper Epinephelus coioides

In this study, complementary DNA (cDNA) and DNA sequences of major histocompatibility complex (MHC) class IIB genes (mhcIIB) were cloned from orange-spotted grouper Epinephelus coioides. The gene structure of E. coioides mhcIIB consists of five exons and four introns, and its deduced amino acid sequence length is 249 amino acids, including a signal peptide, a peptide-binding region, an IGC1 domain, a transmembrane region and a cytoplasmic tail. A phylogenetic study showed that E. coioides mhcIIB shared 32.0-79.1% identity with those of other teleosts and mammals. Real-time reverse transcriptase (RT)-PCR was performed to detect the class IIB gene expression in eight different tissues. To characterize the relationship between E. coioides mhcIIB gene and pathogens, in vivo and in vitro studies were performed. Challenge of Cryptocaryon irritans revealed that class IIB genes were down-regulated after 24 and 48 h of challenge, and their expression was later restored at 72 h. Stimulation of isolated E. coioides leukocytes with lipopolysaccharide (LPS) and polyinosinic:polycytidylic acid (PolyI:C) significantly increased peripheral blood and spleen mhcIIB expression, while head kidney mhcIIB expression remained constant.

Major histocompatibility complex class IIA and IIB genes of the spotted halibut Verasper variegatus: genomic structure, molecular polymorphism, and expression analysis

The major histocompatibility complex (MHC) is a large genomic region characterized by extremely high polymorphism and its association with resistance/susceptibility to disease in vertebrates. In this study, the full lengths of MHC IIA and IIB cDNA were obtained from spotted halibut (Verasper variegates) by expressed sequence tag (EST) and rapid amplification of cDNA ends (RACE) approaches. The genomic structure, molecular polymorphism, and expression patterns were examined to study MHC II gene functions in fish. As in other teleosts, the genomic structure of the spotted halibut MHC IIA contained 4 exons and 3 introns. The deduced amino acid sequence of the class IIA molecule shared 28-79% similarity with those of teleosts and mammals. Nine class IIA alleles were identified from five individuals. Three alleles originating from a single individual suggested the existence of at least two class IIA loci in the genome. Six exons and 5 introns were identified from spotted halibut MHC IIB, and the deduced amino acid sequence shared 33-79% similarity with those of teleosts and mammals. Twelve alleles were identified, among which five were observed in a single individual, which suggested at least three class IIB loci. Quantitative real-time PCR analysis revealed the presence of class IIA and IIB transcripts in nine normal tissues with high expression level in kidney and gill. Furthermore, MHC IIA and IIB are probably two candidates of immune molecules involved in the acute-phase response in spotted halibut, because their transcriptional levels were significantly up-regulated in blood and liver after bacterial challenge.

MHC polymorphism and disease resistance to Vibrio anguillarum in 8 families of half-smooth tongue sole (Cynoglossus semilaevis)

Background

Genes in the major histocompatibility complex (MHC) have a critical role in both the innate and adaptive immune responses because of their involvement in presenting foreign peptides to T cells. However, the nature has remained largely unknown.

Results

We examined the genetic variation in MHC class IIB in half-smooth tongue sole (Cynoglossus semilaevis) after challenge with vibrio anguillarum. Two thousand and four hundred fry from 12 half-smooth tongue sole families were challenged with Vibrio anguillarum. To determine any association between alleles and resistance or susceptibility to V. anguillarum, 160 individuals from four high-resistance (HR, < 40.55% mortality) families and four low-resistance (LR, > 73.27% mortality) families were selected for MHC IIB exon2 gene sequence analysis. The MHC IIB exon2 genes of tongue sole displayed a high level of polymorphism and were discovered at least four loci. Meanwhile, the d_N/d_S [the ratio of non-synonymous (d_N) substitutions to synonymous (d_S) substitutions] in the peptide-binding region (PBR) was higher than that in the non-peptide-binding region (non-PBR). Eighty-eight alleles were discovered among 160 individuals, and 13 out of 88 alleles were used to analyze the distribution pattern between the resistant and susceptible families. Certain alleles presented in HR and LR with a different frequency, while other alleles were discovered in only the HR or LR families, not both. Five alleles, Cyse-DBB*6501, Cyse-DBB*4002, Cyse-DBB*6102, Cyse-DBB*5601 and Cyse-DBB*2801, were found to be associated with susceptibility to V. anguillarum with a frequency of 1.25%, 1.25%, 1.25%, 1.25% and 2.5% in the HR families, and 35%, 33.75%, 27.5%, 16.25%, 15% in the LR families (p < 0.01, 0.01, 0.01, 0.01, 0.01), respectively. Four alleles, Cyse-DBB*3301, Cyse-DBB*4701, Cyse-DBB*6801 and Cyse-DBB*5901, were found to be associated with resistance to V. anguillarum, with a frequency of 13.75%, 11.25%, 11.25%, 8.75% in the HR families and 1.25%, 1.25%, 1.25%, 1.25% and 1.25% in the LR families (p < 0.01, 0.05, 0.05 and p = 0.064), respectively.

Conclusions

Elucidation of the role of MHC II B genes in half-smooth tongue sole should prove to be helpful to the in-depth development of marker-assisted selective breeding in half-smooth tongue sole.

Characterization of the major histocompatibility complex class II genes in miiuy croaker

Major histocompatibility complex (MHC) has a central role in the adaptive immune system by presenting foreign peptide to the T-cell receptor. In order to study the molecular function and genomic characteristic of class II genes in teleost, the full lengths of MHC class IIA and IIB cDNA and genomic sequence were cloned from miiuy croaker (Miichthys miiuy). As in other teleost, four exons and three introns were identified in miiuy croaker class IIA gene; but the difference is that six exons and five introns were identified in the miiuy croaker class IIB gene. The deduced amino acid sequence of class IIA and class IIB had 26.3–85.7% and 11.0–88.8% identity with those of mammal and teleost, respectively. Real-time quantitative RT-PCR demonstrated that the MHC class IIA and IIB were ubiquitously expressed in ten normal tissues; expression levels of MHC genes were found first upregulated and then downregulated, and finally by a recovery to normal level throughout the pathogenic bacteria infection process. In addition, we report on the underlying mechanism that maintains sequences diversity among many fish species. A series of site-model tests implemented in the CODEML program revealed that positive Darwinian selection is likely the cause of the molecular evolution in the fish MHC class II genes.

Molecular cloning, genomic structure and expression analysis of major histocompatibility complex class Iα gene of half-smooth tongue sole (Cynoglossus semilaevis)

Major histocompatibility complex (MHC) has a central role in the adaptive immune system by presenting foreign peptide to the T-cell receptor. The full length of MHC class Iα cDNA was cloned from half-smooth tongue sole by homology cloning and rapid amplification of cDNA ends polymerase chain reaction (RACE-PCR), genomic organization and expression of MHC Iα were examined to study the function of MHC gene in fish. The domain structure feature and antigen-binding motifs of other teleost and mammals MHC are conserved in the half-smooth tongue sole MHC Iα gene. The deduced amino acid sequence of half-smooth tongue sole MHC Iα (GenBank accession no. FJ372720) had 12.1-61.8% identity with those of human and other fish. Eight exons and seven introns were identified in MHC Iα gene. Real-time quantitative PCR demonstrated that MHC Iα gene was ubiquitously expressed in normal tissues, while that in Vibrio anguillarum infected fish was significantly increased in intestines and decreased in spleen and liver from 24 to 72 h after infection, followed by a recovery to normal level after 96 h.

Allelic polymorphism, gene duplication and balancing selection of the MHC class II DAB gene of Cynoglossus semilaevis (Cynoglossidae)

Major histocompatibility complex (MHC) genes play an important role in the immune response of vertebrates. Allelic polymorphism and evolutionary mechanism of MHC genes have been investigated in many mammals, but much less is known in teleosts. We examined the polymorphism, gene duplication and balancing selection of the MHC class II DAB gene of the half-smooth tongue sole (Cynoglossus semilaevis); 23 alleles were found in this species. Gene duplication manifested as three to six distinct sequences at each domain in the same individuals. Non-synonymous substitutions occurred at a significantly higher frequency than synonymous substitutions in the PBR domain, suggesting balancing selection for maintaining polymorphisms at the MHC II DAB locus. Many positive selection sites were found to act very intensely on antigen-binding sites of MHC class II DAB gene.

Gene duplication and evidence for balancing selection acting on MHC class II DAA gene of the half-smooth tongue sole (Cynoglossus semilaevis)

Allelic polymorphism and evolution mechanism of major histocompatibility complex (MHC) genes has been investigated in many mammals, however, much less is known in teleost. In order to investigate the mechanisms creating and maintaining variability at the MHC class II DAA locus, we examined the polymorphism, gene duplication and balancing selection of MHC class II DAA gene of the half-smooth tongue sole (Cynoglossus semilaevis). We described 33 alleles in the C. semilaevis, recombination and gene duplication seems to play more important roles in the origin of new alleles. The rate of non-synonymous substitutions (d(N)) occurred at a significantly higher frequency than that of synonymous substitutions (d(S)) in peptide-binding region (PBR) and non-PBR, suggesting balancing selection for maintaining polymorphisms at the MHC II DAA locus. Many positive selection sites were found to act very intensively on antigen-binding sites. Our founding suggests a snapshot in an evolutionary process of MHC-DAA gene evolution of the C. semilaevis.