Molecular cloning and analysis of common carp CcIgZ3
Identification of the constant region of CcIgZ3 cDNA in Cyprinus carpio.
The obtained sequence of full-length CcIgZ3 cDNA was 2144 bp with a 5’-UTR of 112 bp and a 3’-UTR of 382 bp with a putative typical polyadenylation signal sequence AATAAA located upstream of the poly (A) tail. The deduced CcIgZ3 amino acid sequence contained 549 aa and spanned V domain, four constant domains and a secreted tail (Fig. 1.). The analysis of the IgT sequence showed the presence of a putative 20 amino acids signal peptide. Concerning the V domain of the CcIgZ3 sequences the IMGT unique numbering has been used. The entire C-domain could be divided into four CH domains.
Multiple amino acid sequence alignment
The alignment of common carp CcIgZ3 with zebrafish IgZ1 and IgZ2, common carp IgZ1 and IgZ2, grass carp IgZ and IgZ2 showed that it was comprised of four Ig-like constant domains (CH1, CH2, CH3, CH4) and a secretory tail. The conserved cysteine residues for disulfide bond formation and tryptophan residues for the folding of IgSF domain were present in each CH domain of CcIgZ3 (Fig. 2.). Five N-linked glycosylation sites were predicted to be present in CcIgZ3 (in CH2, CH3 and CH4) (Fig. 1.). The number and distribution of putative N-glycosylation sites including the sequons NXS or NXT for each domain are varied in different teleost species. CH1 domain could possess none to three putative N-glycosylation sites. There was no N-linked glycosylation site in CcIgZ3 CH1 domain, which is different from the grass carp (three N-linked glycosylation sites in CH1) and zebrafish (two N-linked glycosylation sites in CH1). The first N-linked glycosylation site of CcIgZ3 existed in the CH2 domain and it is at the very beginning of CH2. CcIgZ3 presented two N-linked glycosylation sites in CH3 and CH4 domains respectively, which is similar to the grass carp[5].
The percentage of identity between the amino acid sequences of all the IgZ subclasses of common carp and that of grass carp, zebrafish and rainbow trout were calculated considering every single domain. The results showed that the value of identity obtained ranged between 37.5% and 74.76% for all the three IgZ subclasses of common carp, 11.59% and 70.21% for the three IgZ subclasses of common carp with grass carp IgZ/IgZ2, 24.18% and 64.13% for the three IgZ subclasses of common carp with zebrafish IgZ1/IgZ2, 21% and 40.82% for the three IgZ subclasses of common carp with rainbow trout IgT1 ~ IgT3. Interestingly, the identity between CcIgZ3 CH1 domain and grass carp IgZ CH2 domain (71.26%) was higher than that between CcIgZ3 CH1 domain and grass carp IgZ CH1 domain (11.69%), the other two IgZ of common carp also showed this characteristic (41.67% vs 13.1 in IgZ1, 66.67% vs 17.07% in IgZ2).
Phylogenetic analysis
The phylogenetic analysis of immunoglobulins in vertebrates was performed to reveal the relationship of common carp CcIgZ3 with other immunoglobulins. The amino acid sequences from different species including mammals, birds, amphibians and fish were aligned by Clustal X2 program. Based on the alignment, we constructed an NJ tree supported by 1000 bootstrap replications. It showed that common carp CcIgZ3 are grouped with grass carp IgZ and in the same branch as IgZ/T, further demonstrated that CcIgZ3 is a homology of IgZ (Fig. 3).
Basal expression comparison of CcIgZ3 with other Ig molecules of common carp
To investigate the tissue-dependent expression pattern, we performed qRT-PCR analysis using gene-specific primers for all the three IgZ and IgM of common carp in normal adult common carp tissues including liver, spleen, head kidney, gill, skin, gonad, brain, muscle, blood, foregut, hindgut and oral epithelium. The expression of the CcIgZ3, IgZ1, IgZ2 and IgM were found in all the tested tissues and IgM was the most abundant Ig in all the tissues. The expression of CcIgZ3 was found to be highest in the head kidney, gill, and gonad, followed by spleen, hindgut, oral epithelium, liver, brain, muscle, foregut, and blood, but at a very low level in skin. IgZ1 and IgZ2 demonstrated similar expression patterns with CcIgZ3 in most of the detected tissues except that the lowest transcription level of IgZ1 and IgZ2 were found in blood (Fig. 4). To compare the expression of the four IgH transcripts in leukocytes, we performed qRT-PCR analysis using gene-specific primers for all the three IgZ subclasses and IgM of common carp in isolated leukocytes from peripheral blood cells and spleen. The result showed that significant difference in CcIgZ3 expression existed between PBL and leukocytes from the spleen. In contrast, no significant difference in IgM, IgZ1 and IgZ2 expression existed between PBL and leukocytes of the spleen (Fig. 5).
Basal expression comparison of CcIgZ3 with other Ig molecules of common carp during the different developmental stage
When testing for constitutive expression of the four immunoglobulin genes during the early development stage of common carp by qRT-PCR, expression was recorded from 1 to 65 days post fertilization (dpf). (Fig. 6 and Fig. 7). Expression of CcIgZ3 and the other three immunoglobulin genes were detected from 6 days post fertilization (dpf). IgM was expressed highest and was the dominant immunoglobulin isotype during the early development stages tested. (Fig. 6). Furthermore, to investigate the tissue expression pattern of the IgH transcripts, tissues including gill, spleen, hindgut and liver from common carp at 65 days post fertilization (dpf) were used to evaluate the expression level of CcIgZ3, IgM, IgZ1, and IgZ2 respectively. The expression of IgM was still higher than that of the other three subclasses and it was the only isotype that could be detected in all the tested tissues. The gene expression of CcIgZ3, IgZ1 and IgZ2 exhibited varied tissue expression preference. CcIgZ3 was found in spleen, hindgut and gill but not in liver, however, IgZ1 was expressed at a higher level in liver but its expression level in gill was too low to be detected. IgZ2 was primarily detected in gill and spleen but hardly detected in liver and hindgut (Fig. 7).
Tissue location of CcIgZ3 and IgM mRNA-expressing cells at early development stage of common carp
As both of the CcIgZ3 and IgM expression were detected in spleen, gill and hindgut from common carp at 65 dpf, we further detected the CcIgZ3 mRNA-expressing cells and IgM mRNA-expressing cells in tissues by in situ hybridization. Positive expression of both IgM and CcIgZ3 were seen in splenocytes (Fig. 8F and 8H), gill epithelial cells (Fig. 8J and 8L) and cells in lamina propria of hindgut (Fig. 8B and 8D). The staining was specific for both IgM and CcIgZ3 as only the antisense probes (Fig. 8B, 8D, 8F, 8H, 8J, 8L) and not the sense probes (Fig. 8A, 8C, 8E, 8G, 8I, 8K.) hybridized.
Organ-specific CcIgZ3 and IgM expression after Aeromonas hydrophila challenge by different routes
The qRT-PCR results showed that both CcIgZ3 and IgM mRNA expressions were up-regulated compared to the blank control in all tested tissues in the injection and immersion group from 3 days to 21 days except for CcIgZ3 in liver in injection group (0.46-fold at 3 days) (Fig. 9F). Compared with the control, the peak value of CcIgZ3 expression was 86.78-fold in spleen (Fig. 9D), 31-fold in hindgut (Fig. 9H), 13.61-fold in head kidney (Fig. 9B), and 6.39-fold in liver (Fig. 9F) in the immersion group, and 65.13-fold in spleen (Fig. 9D), 10.96-fold in hindgut (Fig. 9H), 26.91-fold in head kidney (Fig. 9B), and 1.21-fold in liver (Fig. 9F) in the injection group respectively. The peak expression of CcIgZ3 in immersion group, in contrast to injection group was comparatively higher in liver (5.28-fold), hindgut (2.83-fold), and spleen (1.33-fold), but lower in head kidney (0.51-fold). However, the peak value of IgM expression was 13.79-fold in head kidney (Fig. 9A), 0.57-fold in spleen (Fig. 9C), 6.5 in liver (Fig. 9E) and 4.07 in hindgut (Fig. 9G) in immersion group and 3.09-fold in head kidney (Fig. 9A), 5.27-fold in spleen (Fig. 9C), 5.85-fold in liver (Fig. 9E), and 5.8-fold in hindgut (Fig. 9G) in injection group. The peak expression of IgM in immersion group, in contrast to injection group were comparatively higher in head kidney (4.46-fod), spleen (1.25-fold), and liver (1.11-fold), but lower in hindgut (0.7-fold). The results of IgM displayed different tissue expression pattern and lower fold change from that of CcIgZ3.
The peak time of CcIgZ3 and IgM were most at 3–7 days in spleen, liver, and hindgut in both immersion group and injection group. However, the peak time of IgM in head kidney was much later in immersion group (14 days) than in injection group (1 day). The peak time of CcIgZ3 in head kidney was also later in immersion group (14 days) than in injection group (3 days). In hindgut, the peak time of CcIgZ3 appeared later in immersion group (3 days) than in injection group (1 day). The immersion challenge in these tissues showed a slower immune response mediated by CcIgZ3 and IgM than injection challenge (Fig. 9).
As for the response in hindgut, the up-regulation of CcIgZ3 expression was significantly higher in the immersion group than the injection group from 3 days to 21 days post challenge, while the expression of IgM had no significant difference from 3 days to 21 days in the hindgut. The value of CcIgZ3 expression in hindgut was 1.73-fold at 1 day, 31-fold at 3 days, 11.53-fold at 7 days, 5.97-fold at 14 days, and 3.76-fold at 21 days in immersion group, and 10.96-fold at 1 day, 1.69-fold at 3 days, 1.13-fold at 7 days, 1.06-fold at 14 days, and 1.09-fold at 21 days in injection group. The ratio of expression of CcIgZ3 in immersion group in contrast to injection group at different time point were 0.15-fold at 1 day, 18.34-fold at 3 days, 10.2-fold at 7 days, 15.63-fold at 14 days, and 3.45-fold at 21 days (Fig. 9H). The value of IgM expression in hindgut was 0.75-fold at 1 day, 2.93-fold at 3 days, 4.07-fold at 7 days, 2.07-fold at 14 days, and 1.73-fold at 21 days in immersion group, and 2.94-fold at 1 day, 2.13-fold at 3 days, 5.80-fold at 7 days, 3-fold at 14 days, and 1.62-fold at 21 days in injection group. The ratio of expression of IgM in immersion group in contrast to injection group at different time points was 0.26-fold at 1 day, 1.38-fold at 3 days, 0.7-fold at 7 days, 0.69-fold at 14 days, and 1.07-fold at 21 days (Fig. 9G). The results showed that Aeromonas hydrophila challenge triggered a stronger CcIgZ3 immune response than IgM in hindgut, especially by the immersion route.