Next Article in Journal
Rab39a and Rab39b Display Different Intracellular Distribution and Function in Sphingolipids and Phospholipids Transport
Previous Article in Journal
Distinct Dopamine D2 Receptor Antagonists Differentially Impact D2 Receptor Oligomerization
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Investigating the Role of BATF3 in Grass Carp (Ctenopharyngodon idella) Immune Modulation: A Fundamental Functional Analysis

1
State Key Laboratory of Freshwater Ecology and Biotechnology, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan 430072, China
2
University of Chinese Academy of Sciences, Beijing 100049, China
3
Innovative Academy of Seed Design, Chinese Academy of Sciences, Beijing 100101, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2019, 20(7), 1687; https://doi.org/10.3390/ijms20071687
Submission received: 8 March 2019 / Revised: 25 March 2019 / Accepted: 27 March 2019 / Published: 4 April 2019
(This article belongs to the Section Molecular Biology)

Abstract

:
Basic leucine zipper transcription factor ATF-like (BATF)-3, belonging to activator protein 1 (AP-1) superfamily transcription factors, is essential for homeostatic development of CD8α+ classical dendritic cells activating CD8 T-cell responses to intracellular pathogens. In this study, the characteristics and cDNA cloning of the CiBATF3 molecule were described in grass carp (Ctenopharyngodon idella). CiBATF3 had abundant expression in immune-related organizations, including liver, spleen and gill, and grass carp reovirus (GCRV) infection had significantly changed its expression level. After Ctenopharyngodon idella kidney (CIK) cells were challenged with pathogen-associated molecular patterns (PAMPs), polyinosinic:polycytidylic acid (poly(I:C)) stimulation induced higher mRNA levels of CiBATF3 than that of lipopolysaccharide (LPS). Subcellular localization showed that CiBATF3-GFP was entirely distributed throughout cells and nuclear translocation of CiBATF3 was found after poly(I:C) treatment. Additionally, the interaction between CiBATF3 and interleukin 10 (IL-10) was proven by bimolecular fluorescence complementation (BiFC) system. The small interfering RNA (siRNA)-mediated CiBATF3 silencing showed that the mRNA of CiBATF3 and its downstream genes were down-regulated in vitro and in vivo. CiBATF3 played a negative regulatory role in the transcriptional activities of AP-1 and NF-κB reporter gene. In summary, the results may provide valuable information on fundamental functional mechanisms of CiBATF3.

1. Introduction

The basic leucine zipper transcription factor ATF-like (BATF) family consists of BATF, BATF2 and BATF3 [1], belonging to the activator protein 1 (AP-1) family of transcription factors [2]. An alpha-helical bZIP domain containing a leucine zipper motif and a basic DNA binding region is characteristic of members of this family [3,4,5]. In every tissue and cell of the immune system, bZIP proteins are crucial transcriptional regulators by interactions with AP-1-binding elements. Two of the earliest identified bZIP proteins, FBJ osteosarcoma oncogene (FOS) and JUN (AP-1 transcription factor subunit), constituting the heterodimeric transcription factor termed AP-1 [6], can widely regulate gene expression. It is known that the leucine zipper motif is essential for interaction with a bZIP protein or a non-bZIP transcription factor such as interferon regulatory factors (IRF) that modulates target genes [7].
Dendritic cells (DCs) are hematopoietic cells that belong to the antigen-presenting cell (APC) family. Their primary function is to process antigens and present them to T cells to promote immunity to foreign antigens and tolerance to self-antigens [8]. Mouse DCs are subdivided into only three main subtypes: two main lineages of classical DCs (cDCs) (“classical type 1 DCs (cDC1s)” for CD8α+ and CD103+ DCs, and “cDC2s” for CD11b+ and CD172a+ DCs based on their distinct developmental pathways) and plasmacytoid DCs (pDCs) [9]. The lineages of cDC1s include both lymph node (LN)-resident and migratory subsets [10,11]. Lymphoid tissue-resident cDCs differentiate in, and spend their entire lives within, lymphoid tissues [12]. Tissue-migratory cDCs are located in the peripheral lymph nodes (LNs). They refer to nonlymphoid tissue DCs that have migrated to the tissue-draining LNs through the lymphatics as opposed to the blood-borne lymphoid-resident DCs [13]. The resident and migratory cDC1s cells are either CD8α+ (LN resident) or CD103+ (migratory). cDC1 development is mediated by the transcription factors IRF8 and BATF3. Mice lacking IRF8 [14] or BATF3 [15] exhibit a severe defect in the development of cDC1s. In particular, IRF8 mutant and IRF8-deficient mice lack CD8α+ cDCs but have additional immune defects [16,17,18]. In BATF3−/− mice, DCs lack virus specific CD8+ T cell in responses to West Nile virus and also have a defect in cross-presentation [15]. Migratory DCs mediate the exchange of archived antigens between lymphatic endothelial cells (LECs) and antigen-presenting cells (APCs). After vaccination, antigen exchange and cross-presentation are performed by both BATF3-dependent and BATF3-independent DCs [19]. In adult and neonatal mice, the intestinal BATF3-dependent DC is vital for optimal antiviral T-cell responses [20]. Recently, many studies [15,21,22] have shown that BATF3 plays very important roles in promoting expansion of functional CD8+ cDCs to control infection of intracellular pathogens, and via the interaction of the conserved LZ domain with IRF4 or IRF8, this function may be compensated by other members of the BATF family [22].
Up to now, research on BATF3 has mostly focused on mammals rather than fish. In fish, the BATF3 genes have been identified in Danio rerio [23], Salmo salar [24], and Oncorhynchus mykiss [25]. In China, Grass carp is one of the most important economic freshwater aquaculture species, and deeply studying its functional genes [26,27,28], whole-genome [29], transcriptome [30,31,32] and microRNA regulation [33] has provided new data for studying the immune system of fish. In our study, we amplified the full-length cDNA sequence of Ctenopharyngodon idella BATF3 (CiBATF3) and, to a better understanding of the role of CiBATF3 gene in fish immune responses, sequence analysis and functional experiments of CiBATF3 were also completed. The bioinformatics and functions analysis of CiBATF3 will help for studying the transcriptional regulation by BATF3 in fish immune systems in the future.

2. Results

2.1. Analysis of CiBATF3

The full-length cDNA of CiBATF3 is 1198 bp with a 5′-terminal untranslated region (UTR) of 147 bp, a 3′-UTR of 688 bp, two mRNA instability motifs (ATTTA), a polyadenylation signal (AATAA), and an open reading frame of 363 bp encoding for 120 amino acid residues with a predicted molecular mass of 13.884 kDa and theoretical isoelectric point of 8.54. SMART analysis revealed that a typical basic leucine zipper (bZIP, 30–87 aa) domain was discovered in the deduced amino acid of CiBATF3 (Figure S1). No signal peptide was detected by the signalP program. The genomic sequence of CiBATF3 ORF, which contained three exons and two introns, was 1621 bp (Figure 1A) and all the intron–exon boundaries displayed the conserved GT/AG rule. The bZIP domain was composed of a leucine zipper (LZ), a hinge (H) region and a DNA binding domain (DB), which existed in all the BATF3. The amino acid sequences of these domains are highly homologous with six absolutely conserved leucine residues. The N- and C-terminal regions of CiBATF3 have a relatively low sequence homology (Figure 1B). Swiss-model analysis indicated the CiBATF3 protein contained a single α-helix, similar to that of human and zebrafish (Figure 1C). A comparison of homology revealed that the deduced CiBATF3 shares 36.9–85% sequence similarity with BATF3 of other vertebrate, with it being the most similar to that of the Danio rerio (85% sequence similarity), followed by the Ictalurus punctatus (69.8% sequence similarity) (Table 1).
To detect the evolutionary conservation in genomic structure, the syntenic analysis of CiBATF3 genes was performed with neighboring genes from other species. The result show that the upstream genes of higher vertebrates BATF3 had relatively conserved synteny, whereas CiBATF3 had highly conserved synteny with other fish BATF3 genes (Figure 2A). Especially, the syntenic relationship of the BATF3 and ATF3 genes in the genome is highly conserved among the selected species.
To determine the molecular evolutionary relationship, BATF family amino acid sequences from various vertebrate species were used to construct a phylogenetic tree. As shown in Figure 2B, three major clades corresponding to BATF1s, BATF2s and BATF3s were detected. All the fish BATF3 proteins clustered together, of which the CiBATF3 branched most closely to Danio rerio BATF3.

2.2. The mRNA Expression Levels of CiBATF3 in Tested Tissues

The spatial mRNA expression profiles of CiBATF3 in all tested tissues (skin (SK), middle kidney (MK), intestine (I), liver (L), gill (G), heart (H), spleen (S), head kidney (HK), muscle (M) and brain (B)) are shown in Figure 3A. The CiBATF3 had the dominant abundance expression in liver, followed by spleen, gills, heart and head kidney, and was low in skin.

2.3. Analysis of CiBATF3, IRF8 and IL-10 Genes After GCRV Infection In Vivo

The temporal expression of the CiBATF3 gene in six immune tissues after GCRV challenge is shown in Figure 3B–G. In the spleen, intestine and head kidney, CiBATF3 were down-regulated on Day 1, and then were up-regulated on Days 3 and 5. In the middle kidney, the CiBATF3 expression levels were promoted during the entire detection period. Expression of GCRV S6 gene increased in gill and intestine on Days 3 and 5 after injection (Figure 3H), which verified infection had occurred.
To find the antiviral mechanism in grass carp after GCRV infection, the expression levels of IRF8 (Figure 4A–D) and IL-10 (Figure 4E–H), which were regulated by CiBATF3, were detected in four immune tissues (spleen, intestine, head kidney and middle kidney). After GCRV infection, interestingly, the response patterns of IRF8 and IL-10 were up-regulated in the four examined tissues similar to CiBATF3, respectively.

2.4. Time-Course Analysis of CiBATF3 Expression in CIK Cells After LPS Exposure or Poly(I:C) Challenge

Time-course expression of CiBATF3 expression in CIK cells after LPS exposure is shown in Figure 5A. The expression level of CiBATF3 gradually increased from 0 h to 6 h, reached the peak at 6 h (1.54-fold, p < 0.05), and then decreased until 48 h. Upon poly(I:C) stimulation, the CiBATF3 transcripts were induced from 3 h to 48 h (Figure 5B), and reached peak at 36 h, at which time the expression level was 18.36-fold (p < 0.01) compared to the control group (0 h).

2.5. siRNA-Mediated CiBATF3 Silencing Could Affect the Expression of its Downstream Molecules In Vitro and In Vivo

The effect of siRNA-mediated CiBATF3 silencing on the expression of its downstream genes was examined by RT-qPCR analysis, including GATA binding protein 3 (GATA3), interleukin 4 (IL-4), interleukin 10 (IL-10), interferon regulatory factor 8 (IRF8), interleukin 12a (IL-12 p35), interleukin 12b (IL-12 p40) and cellular myelocytomatosis oncogene (c-myc), because it has been reported that these genes are regulated by BATF family members and involved in immune response [7,22,35,36]. In the in vitro assay, siRNA transfection depressed BATF3 mRNA level in CIK cells (Figure 6A). In this condition, the mRNA levels of IL-4, IL-10 and IL-12 p40 displayed a 0.49-fold, 0.13-fold and 0.15-fold decrease, respectively (Figure 6B). However, these genes, including GATA3, IRF8, IL-12 p35 and c-myc, did not change significantly compared to negative the control group.
For in vivo test, the CiBATF3 mRNA were down-regulated 0.17-fold (p < 0.01) and 0.33-fold (p < 0.01) in middle kidney (MK) and head kidney (HK), respectively (Figure 7A), after siRNA silencing of CiBATF3. Most of its downstream genes were significantly down-regulated in middle kidney and head kidney except that IRF8 and IL-12 p35 did not change significantly in head kidney at this condition (Figure 7B–H).

2.6. Cell Transfection and Luciferase Activity Analysis

To investigate the function of CiBATF3, HEK293 cells were co-transfected with empty vector or pCDNA3.1-CiBATF3 together with NF-κB or AP-1 reporter gene. As shown in Figure 8, these results indicated that CiBATF3 can suppress the transcriptional activities of AP-1 and NF-κB reporter gene.

2.7. Subcellular Localization of CiBATF3 in CIK Cells and HEK293 Cells

The subcellular localization of the CiBATF3 protein was observed using confocal microscope in CIK cells or HEK293 cells. In CIK cells, CiBATF3-pEGFP was strongly dispersed throughout the cytoplasm and the nuclear areas, which was consistent with pEGFP-N3 (Figure 9A). In HEK293 cells, CiBATF3-pEGFP was also distributed same as that of CIK cells (Figure 9B).

2.8. The CiBATF3-GFP Nuclear Translocation Induced by Poly(I:C) Stimulation

As shown in Figure 9C, PAMPs induce the relocation of CiBATF3 to varying degrees. In poly(I:C) stimulated-cells, CiBATF3 proteins appeared to have a tendency to concentrate in the nucleus. After LPS stimulation, the protein expression pattern of CiBATF3 did not change.

2.9. Interaction between CiBATF3 and IL-10

The interaction of CiBATF3 and IL-10 in living cells was visualized by the mNeptune-based BiFC system. No fluorescence signal was observed when pBATF3-MN155 or pMC156-IL-10 was transfected into CIK cells alone. However, after co-transfection of pBATF3-MN155 and pMC156-IL-10, the red mNeptune fluorescence signal in the cytoplasm was bright red (Figure 9D). Thus, it confirmed the interaction between CiBATF3 and IL-10 in CIK cells.

3. Discussion

It is known that BATF3 is crucial for regulating numerous cellular processes [7] and immune-specific functions [15,21,37,38]. In fish, BATF3 signaling attracts attention [25]. Here, a BATF family member (CiBATF3) was firstly cloned and identified in C. idella. Multiple sequence alignments, domain structure and phylogenetic analyses strongly showed functional homology of CiBATF3 orthologs in other species. The high expressions of CiBATF3 in immune-related tissues (liver, spleen, gills and head kidney) suggested that it plays critically important roles in the immune system of grass carp. Thus, expression pattern of CiBATF3 gene post GCRV, LPS and poly(I:C) stimulation were further analyzed.
Toll-like receptors (TLRs) [39,40,41,42,43,44] play a key role in the innate immune response in mammals by identifying conserved molecular patterns associated with different microorganisms. TLR4 has been identified as a signaling molecule essential for the responses to LPS [45]. TLR3 can be stimulated by double-stranded RNA [46]. Poly(I:C) is well known as a dsRNA analog. In rainbow trout, Poly(I:C) and R848 stimulation dramatically up-regulate BATF3a and BATF3b. After LPS treatment, BATF3a is also up-regulated to some extent. However, no significant changes in BATF3b are detected in cells treated with LPS and Flagellin [25]. This indicates that BATF3 may be involved in the TLRs signaling pathway. In this study, the mRNA levels of CiBATF3 were induced by poly(I:C) and LPS. However, poly(I:C) can induce higher mRNA levels of CiBATF3 than LPS. Moreover, CiBATF3 proteins localization appeared to have a tendency to concentrate in the nucleus following poly(I:C) stimulation but there was no change after LPS stimulation (Figure 9C). Therefore, these suggest that CiBATF3 may be more sensitive to the TLR3 signaling pathway than TLR4 signaling pathway.
It is reported that BATF family members have the function of regulating the gene expression. In BATF−/− BATF3−/− TH2 cells, the expressions of IL-4 and IL-10 are significantly lost, indicating that BATF and BATF3 act together on the expression of these two cytokine genes [22]. BATF contributes to express the major TH2 cell-associated transcription factor GATA3 and also has direct effects on the genes encoding other TH2 cell effector factors (e.g., cytotoxic T lymphocyte antigen 4 (CTLA4), IL-4 and IL-10). BATF3 also helps regulate these TH2 cell-associated genes [7]. STAT-mediated BATF3 expression is crucial for lymphoma cell survival and promotes MYC activity in classical Hodgkin lymphoma and anaplastic large cell lymphoma, and a new oncogenic axis is recognized in these lymphomas [35]. The evidence for the direct role of BATF3 in regulating IRF8 expression has been discovered. In CD24+ cDCs, the auto-activation of IRF8, dependent on BATF3 may operate at an element +32 kb from the IRF8 transcriptional start site (TSS), within an IRF8 superenhancer [36]. In this study, siRNA-mediated CiBATF3 silencing showed that CiBATF3 can regulate most of the tested gene expression levels.
According to many studies, transcriptional regulator NF-κB can activate downstream anti-apoptotic genes and promote cell survival [47,48,49,50]. AP-1 has been proven to promote cell proliferation and malignant transformation [27,51,52]. In mice, BATF negatively regulates the transactivation of an AP-1 reporter gene in vivo [53]. BATF2 inhibites both AP-1 and NF-κB reporter gene transcriptional activities [54]. In our study, CiBATF3 negatively regulated the transcriptional activities of AP-1 and NF-κB reporter gene. Therefore, we speculate that CiBATF3 protein can inhibit cell proliferation by inhibiting AP-1 activity. Inhibition of NF-κB activity by CiBATF3 is likely to attenuate the effects of anti-apoptosis, thereby enhancing the effects of pro-apoptotic.
In conclusion, the results of this study may provide valuable information for the study of the potential functional mechanism of BATF3 in teleosts.

4. Materials and methods

4.1. Cells, Plasmid and Fish

The HEK293 cells were cultured in high glucose Dulbecco’s modified Eagle’s medium (DMEM; Hyclone, Logan, UT, USA), with 10% FBS, 100 mg/mL streptomycin (Sigma) and 100 IU/mL penicillin (Sigma) under a humidified condition with 5% CO2 at 37 °C. The CIK cells (provided by China Center for Type Culture Collection, Wuhan, China) were cultured under a humidified condition with 5% CO2 at 28 °C, using Medium 199 (Sigma, St. Louis, MO, USA) with 10% FBS, 100 mg/mL streptomycin (Sigma) and 100 IU/mL penicillin (Sigma). The pEGFP-N3 plasmid was previously in our lab. The pMN155 and pMC156 plasmids were given by Professor Zongqiang Cui, Wuhan Institute of Virology, Chinese Academy of Sciences. pMD18-T were obtained from TransGen Biotech (Beijing, China). pRL-TK, pGL3-bacic, pCDNA3.1 (-), pNFқB-TA-luc, and pAP1-TA-luc were purchased from Beyotime (Shanghai, China). Grass carps (3 months old; weight, 10 ± 2 g; length, 7 ± 3 cm) were collected at the GuanQiao Experimental Station, Institute of Hydrobiology, Chinese Academic of Sciences, and acclimatized in aerated freshwater at 28 °C for one week. The fish were fed with a commercial feed (Tong Wei, Chengdu, China) to adapt to the environment until 24 h before the experiments under the same conditions.

4.2. GCRV Infection and Sampling

For the tissue distribution experiment, tissues including liver, skin, intestine, gill, spleen, muscle, heart, head kidney, brain and middle kidney were collected from five healthy grass carp.
The GCRV infection experiment was carried out as described previously [55] with modifications. Five healthy grass carps were injected intraperitoneally with 200 μL PBS (pH 7.4) as control group while each fish of challenge group was injected with an equal volume of GCRV (GD108 strain). The titer of virus detected by RT-qPCR (the special primers are listed in Table 2) was 3.12 × 103 copy/μL. These injected fish were kept under the same conditions as mentioned above with commercial feed during the test. Five individuals were collected at 0, 1, 2, 3, 4, 5 and 6 days post injection (dpi). The head kidney, gill, liver, intestine, middle kidney and spleen were harvested in TRIzol reagent (Invitrogen, Carlsbad, CA, USA) and stored at −80 °C until RNA extraction.
Total RNAs were extracted and the first-strand cDNA was synthesized using ReverTra Ace kit (Toyobo, Osaka, Japan). The reaction was performed at 42 °C for 1 h, and terminated by heating at 95 °C for 5 min.

4.3. Cloning the Full-Length CiBATF3 cDNA

The specific fragments with the coverage of open reading frame (ORF) region were obtained by blasting the CiBATF3 sequences of zebrafish (Accession no. NM_001045392.1) with the C. idella transcriptional database [29]. Then, the ORF sequence was amplified and then the 5ʹ and 3ʹ untranslated regions (UTR) were cloned using the 5′ and 3′ Full RACE Kit (TaKaRa, Kusatsu, Japan) according to the manufacturer’s protocol. The PCR fragments were purified by gel extraction and then were incubated with linearized pMD18-T vector and T4 ligase at 16 °C overnight. The recombinant vectors were transfected to E. coli DH5α (TransGen, Beijing, China), and then monoclonal E. coli were sequenced by Tsing Ke company. Finally, the CiBATF3 full length cDNA was assembled using DNAMAN software.

4.4. Sequence Analysis

The Sequence Manipulation Suite (STS) (http://www.bio-soft.net/sms/) were used to analyze the amino acid sequence and nucleotide sequence of CiBATF3. The CiBATF3 full-length cDNA sequence was aligned to the genomic sequence using BLAST to determine the intron/exon structure of the genomic sequence. The protein domain characteristics were analyzed using the Simple SMART program (http://smart.emblheidelberg.de/) and the CDD tool (http://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi). The multiple sequence alignment was performed with ClustalW2 (http://www.ebi.ac.uk/Tools/msa/clustalw2/) tool. The the amino acid sequence and nucleotide sequence from other species were collected by NCBI. Gene synteny relationship of BATF3 genes were obtained by GnomicusV88.01 [34], and then these genes used as references were compared with grass carp genomic sequence. The signal peptide was predicted using the SignaIP 4.0 Server (http://www.cbs.dtu.dk/services/SignalP/). The tertiary structure was predicted using SWISS-MODEL (http://swissmodel.expasy.org/) and displayed by Swiss-PdbViewer Version 4.1. The phylogenetic tree of CiBATF3 was constructed by MEGA v 5.0 program based on neighbor-joining (NJ) method [56]. The amino acid sequence of CiBATF3 was comparatively analyzed with orthologous sequences by Geneious software to assess the similarity percentages.

4.5. Tissue Distribution of CiBATF3

To understand the expression profile of CiBATF3 in different tissues, quantitative real-time PCR (RT-qPCR) was carried out using iQ™ SYBR Green Supermix (Bio-Rad, Hercules, CA, USA ) on a CFX96TM Real Time Detection System (Bio-Rad, Hercules, CA, USA). Ctenopharyngodon idella β-actin (Accession No. M25013.1) was selected as reference gene. The specific primers are listed in Table 2. Finally, the Ct values for respective reaction were subjected to comparative Ct method (2−ΔΔCT) [57] to calculate the CiBATF3 mRNA levels. Data are presented as the CiBATF3 expression-fold relative to that of β-actin mRNA and are presented as mean ± standard deviation (SD).

4.6. Responses of CiBATF3, IRF8 and IL-10 Genes to GCRV Infection In Vivo

RT-qPCR was employed for examining the responses of CiBATF3, IRF8 and IL-10 following GCRV infection. β-actin (Accession No. M25013.1) was selected as internal control. Relative expression levels were calculated as the ratio of gene expression in the infected grass carp at each time point (1, 3 and 5 days after GCRV infection) relative to that in the uninfected fish (Day 0). To test the effectiveness of GCRV infection, the relative copy numbers of the virus in gill and intestine were examined by RT-qPCR. The expression level of Day 1 was set as the baseline (1.0). Special primers were designed based on the S6 segments of GCRV-II (see Table 2).

4.7. The mRNA Expression Profiles of CiBATF3 Gene in CIK Cells Following Poly(I:C) Challenge, LPS Stimulation

To further investigate the effect on CiBATF3 in vitro post PAMPs challenges, the expression patterns of RT-qPCR were employed. For LPS stimulation, cells were exposed to 10 μg mL−1 LPS (Escherichia coli 055: B5, Sigma, USA) for 1, 2, 6, 12, 24 and 48 h. Untreated cells served as controls. For poly(I:C) challenge, cells were exposed to 20 μg mL−1 poly(I:C) (sigma, USA) for 0, 3, 8, 24, 36, and 48 h. The cells at 0 h served as controls. Three biological replicate cells were collected and lysed in Trizol for subsequent RT-qPCR analysis.

4.8. CiBATF3 Gene Silencing Using siRNA In Vitro and In Vivo

Small interfering RNA (siRNA) that specifically targets CiBATF3 and negative control were synthesized by GenePharma (Shanghai, China). The siRNA sequence information is listed in Table 2. For transfection experiments in vitro, 5 μL CiBATF3 specific siRNA (20 mM) or negative controls were transfected into CIK cells using Lipo6000™ Transfection Reagent (Beyotime). At 24 h post-transfection, the treated and control CIK cells were collected for expression analysis. Three biological replicates were obtained for each group.
The transfection experiments in vivo were performed as described previously [58] with some modifications. Briefly, 10 μL of CiBATF3 specific siRNA (20 mM) or the negative control was mixed with 10 μL of Lipo6000™ Transfection Reagent and 80 μL of PBS to serve as the transfection solution. Grass carp (weight, 10 ± 2 g; length, 7 ± 3 cm) were injected intraperitoneally with 100 μL of above transfection solution. After 24 h transfection, the head kidney and middle kidney tissues from treated or control group were collected for expression analysis using RT-qPCR. There were three biological replications at each time point.

4.9. Dual-Luciferase Activity Assays

The ORF sequence encoding the mature peptide of CiBATF3 was amplified by specific primers (Table 3) with Xho I/EcoR I sites, and then ligated into pCDNA3.1 (-) eukaryotic expression vector. All recombinant vectors were transformed into E. coli DH5a and confirmed by sequencing. Then, positive plasmids were extracted by Endo-free plasmid mini kit I.
The process of luciferase assays was performed as described previously [26]. For the transient transfection assays, HEK293 cells were transfected with 0.25 μg pCDNA3.1-CiBATF3, 0.25 μg reporter gene vectors (pNFқB-TA-luc and pAP1-TA-luc, respectively), 0.025 μg pRL-TK renilla luciferase vector and 1 μL Lipo6000™ Transfection Reagent. At 24 h post- transfection, the cells were collected and lysed. Luciferase activities were measured using the Luciferase Assay System (Promega, Wisconsin, USA). The experiments were repeated at least three times. Differences were considered significant at p < 0.05.

4.10. Subcellular Localization of CiBATF3 in CIK Cells and HEK293 Cells

Specific primer pairs containing XhoI and BamHI restriction enzyme cutting sites (Table 3) were designed to amplify the complete ORF sequence of CiBATF3. After being digested by restriction enzyme, the PCR product was cloned into pEGFP-N3 vectors producing GFP-tagged expression plasmids. Subsequently, CIK cells and HEK293T cells were seeded onto a sterile microscope cover glass and placed in a 6-well cell culture plate prior to transfection. After cells were grown to 70–90% confluence, 5 μg of pEGFP-N3-CiBATF3 vectors were transiently transfected, while empty pEGFP-N3 vector was used as the control. After transfection for 24 h, CIK cells and HEK293T cells were washed three times with PBS and fixed with 4% (v/v) paraformaldehyde for 15 min at room temperature. Then, CIK cells and HEK293T cells were dyed with Hoechst 33342 (Beyotime, China) and observed under the UltraVIEW VOX confocal system (PerkinElmer, Fremont, CA, USA) and a 63× oil immersion objective lens.
To evaluate the influence of PAMPs stimulation on subcellular localization of CiBATF3 proteins, HEK293T cells were transfected with CiBATF3-EGFP fusion plasmid or control plasmid (pEGFP-N3) and seeded into six-well plates. Approximately 24 h later, the cells were treated with 50 mg/mL poly(I:C) and 20 mg/mL LPS, respectively. After PAMPs stimulation, the HEK293T cells were observed using confocal system.

4.11. Verification of the Protein Interaction between CiBATF3 and IL-10

Interleukin-10 (IL-10) is known as an anti-inflammatory cytokine, which plays critical inhibitory roles in broad spectrum of immune responses [59]. It was predicted by online InBio Map™ version 2016_09_12 (https://www.intomics.com/inbio/map.html) that BATF3 interacts with IL-10. The bimolecular fluorescence complementation (BiFC) system is based on the reconstitution of two non-fluorescent fragments of a fluorescent protein, which could be used to determine whether two proteins are interactive [60,61]. In the present study, the mNeptune-based BiFC system was introduced to visualize whether IL-10 interacts with CiBATF3 in CIK cells. Briefly, the ORF sequences of CiBATF3 and IL-10 (amplified by primers in Table 3) were cloned from grass carp and inserted into the pMN155 and pMC156 plasmids, respectively. The final plasmids were named pBATF3-MN155 and pMC156- IL-10, which contained the N-terminal of mNeptune (mNeptune aa 1–155, MN155) and C-terminal of mNeptune (mNeptune aa 156C-terminal, MC156), respectively. Then, plasmids pBATF3-MN155 and pMC156- IL-10 were transfected into CIK cells alone or together. At 24 h post-transfection, cells were observed as described above.

4.12. Statistical Analysis

All the data were analyzed using the SPSS software (version 22.0 USA) and assessed by one-way analysis of variance (ANOVA). All the experiments were repeated at least three times. The significance level was set at p < 0.05 and the highest significance level was set at p < 0.01.

Supplementary Materials

Supplementary Materials can be found at https://www.mdpi.com/1422-0067/20/7/1687/s1. Figure S1. Deduced amino acid sequence and nucleotide sequence of CiBATF3. The start codon (ATG) and the stop codon (TAA/TGA) were boxed. The shaded region is the bZIP domain. The motif associated mRNA instability (ATTTA) is shown as double underscore, and the polyadenylation signal motif (AATAA) is shown as underscore.

Author Contributions

R.H., Y.W., Z.Z. and D.Z. designed the experiments. R.H., L.H., Y.L. and L.L. (Lanjie Liao) supported sampling. D.Z., R.H., P.C. and P.F. contributed to data analysis. D.Z., L.L. (Lifei Luo) and L.C. performed the RT-qPCR validation and graphic analysis. D.Z. and R.H. drafted the manuscript. Y.W. and Z.Z. provided advice on manuscript writing. All authors reviewed the manuscript.

Funding

This work was supported by the Strategic Priority Research Program of the Chinese Academy of Sciences (grant numbers XDA08030203) and the National Natural Science Foundation of China (grant numbers 31572614 and 31721005).

Acknowledgments

We appreciate Zongqiang Cui (Wuhan Institute of Virology, Chinese Academy of Sciences) for giving the pMN155 and pMC156 plasmids.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

BATF3basic leucine zipper transcription factor ATF-like (BATF)-3
STAT3signal transducer and activator of transcription 3
DBDNA binding domain
AP-1activator protein 1
DCdendritic cell
IFNinterferon
ILinterleukin
IRFinterferon regulated factor
CIKCtenopharyngodon idella kidney
LPSlipopolysaccharide
FOSFBJ osteosarcoma oncogene
LZleucine zipper
c-myccellular myelocytomatosis oncogene
PAMPpathogen-associated molecular pattern
Poly(I:C)polyinosinic:polycytidylic acid
bZIPbasic leucine zipper
CDcluster of differentiation
GATA3GATA binding protein 3
IL-4interleukin 4
IL-10interleukin 10
IRF8interferon regulatory factor 8
IL-12 p35interleukin 12a
IL-12 p40interleukin 12b
JUNAP-1 transcription factor subunit
dsRNAdouble-stranded RNA
PGNpeptidoglycan

References

  1. Dorsey, M.J.; Tae, H.J.; Sollenberger, K.G.; Mascarenhas, N.T.; Johansen, L.M.; Taparowsky, E.J. B-ATF: A novel human bZIP protein that associates with members of the AP-1 transcription factor family. Oncogene 1995, 11, 2255–2265. [Google Scholar] [PubMed]
  2. Finn, R.D.; Bateman, A.; Clements, J.; Coggill, P.; Eberhardt, R.Y.; Eddy, S.R.; Heger, A.; Hetherington, K.; Holm, L.; Mistry, J.; et al. Pfam: The protein families database. Nucleic Acids Res. 2014, 42, D222–D230. [Google Scholar] [CrossRef] [PubMed]
  3. Landschulz, W.H.; Johnson, P.F.; McKnight, S.L. The leucine zipper: A hypothetical structure common to a new class of DNA binding proteins. Science 1988, 240, 1759–1764. [Google Scholar] [CrossRef] [PubMed]
  4. Turner, R.; Tjian, R. Leucine repeats and an adjacent DNA binding domain mediate the formation of functional cFos-cJun heterodimers. Science 1989, 243, 1689–1694. [Google Scholar] [CrossRef] [PubMed]
  5. O’Shea, E.K.; Rutkowski, R.; Stafford, W.F., 3rd; Kim, P.S. Preferential heterodimer formation by isolated leucine zippers from fos and jun. Science 1989, 245, 646–648. [Google Scholar] [CrossRef] [PubMed]
  6. Lee, W.; Mitchell, P.; Tjian, R. Purified transcription factor AP-1 interacts with TPA-inducible enhancer elements. Cell 1987, 49, 741–752. [Google Scholar] [CrossRef]
  7. Murphy, T.L.; Tussiwand, R.; Murphy, K.M. Specificity through cooperation: BATF–IRF interactions control immune-regulatory networks. Nat. Rev. Immunol. 2013, 13, 499. [Google Scholar] [CrossRef]
  8. Merad, M.; Manz, M.G. Dendritic cell homeostasis. Blood 2009, 113, 3418–3427. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Guilliams, M.; Ginhoux, F.; Jakubzick, C.; Naik, S.H.; Onai, N.; Schraml, B.U.; Segura, E.; Tussiwand, R.; Yona, S. Dendritic cells, monocytes and macrophages: A unified nomenclature based on ontogeny. Nat. Rev. Immunol. 2014, 14, 571–578. [Google Scholar] [CrossRef]
  10. Worbs, T.; Hammerschmidt, S.I.; Forster, R. Dendritic cell migration in health and disease. Nat. Rev. Immunol. 2017, 17, 30–48. [Google Scholar] [CrossRef] [PubMed]
  11. Durai, V.; Murphy, K.M. Functions of Murine Dendritic Cells. Immunity 2016, 45, 719–736. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Randolph, G.J.; Ochando, J.; Partida-Sanchez, S. Migration of dendritic cell subsets and their precursors. Annu. Rev. Immunol. 2008, 26, 293–316. [Google Scholar] [CrossRef] [PubMed]
  13. Randolph, G.J.; Angeli, V.; Swartz, M.A. Dendritic-cell trafficking to lymph nodes through lymphatic vessels. Nat. Rev. Immunol. 2005, 5, 617–628. [Google Scholar] [CrossRef] [PubMed]
  14. Schiavoni, G.; Mattei, F.; Sestili, P.; Borghi, P.; Venditti, M.; Morse, H.C.; Belardelli, F.; Gabriele, L. ICSBP is essential for the development of mouse type I interferon-producing cells and for the generation and activation of CD8 alpha(+) dendritic cells. J. Exp. Med. 2002, 196, 1415–1425. [Google Scholar] [CrossRef]
  15. Hildner, K.; Edelson, B.T.; Purtha, W.E.; Diamond, M.; Matsushita, H.; Kohyama, M.; Calderon, B.; Schraml, B.U.; Unanue, E.R.; Diamond, M.S.; et al. Batf3 Deficiency Reveals a Critical Role for CD8 alpha(+) Dendritic Cells in Cytotoxic T Cell Immunity. Science 2008, 322, 1097–1100. [Google Scholar] [CrossRef] [PubMed]
  16. Tailor, P.; Tamura, T.; Morse, H.C.; Ozato, K. The BXH2 mutation in IRF8 differentially impairs dendritic cell subset development in the mouse. Blood 2008, 111, 1942–1945. [Google Scholar] [CrossRef] [PubMed]
  17. Tsujimura, H.; Tamura, T.; Gongora, C.; Aliberti, J.; Sousa, C.R.E.; Sher, A.; Ozato, K. ICSBP/IRF-8 retrovirus transduction rescues dendritic cell development in vitro. Blood 2003, 101, 961–969. [Google Scholar] [CrossRef] [PubMed]
  18. Gabriele, L.; Ozato, K. The role of the interferon regulatory factor (IRF) family in dendritic cell development and function. Cytokine Growth Factor Rev. 2007, 18, 503–510. [Google Scholar] [CrossRef]
  19. Kedl, R.M.; Lindsay, R.S.; Finlon, J.M.; Lucas, E.D.; Friedman, R.S.; Tamburini, B.A.J. Migratory dendritic cells acquire and present lymphatic endothelial cell-archived antigens during lymph node contraction. Nat. Commun. 2017, 8. [Google Scholar] [CrossRef] [PubMed]
  20. Sun, T.; Rojas, O.L.; Li, C.; Ward, L.A.; Philpott, D.J.; Gommerman, J.L. Intestinal Batf3-dependent dendritic cells are required for optimal antiviral T-cell responses in adult and neonatal mice. Mucosal Immunol. 2017, 10, 775. [Google Scholar] [CrossRef]
  21. Edelson, B.T.; Wumesh, K.C.; Juang, R.; Kohyama, M.; Benoit, L.A.; Klekotka, P.A.; Moon, C.; Albring, J.C.; Ise, W.; Michael, D.G.; et al. Peripheral CD103(+) dendritic cells form a unified subset developmentally related to CD8 alpha(+) conventional dendritic cells. J. Exp. Med. 2010, 207, 823–836. [Google Scholar] [CrossRef] [PubMed]
  22. Tussiwand, R.; Lee, W.L.; Murphy, T.L.; Mashayekhi, M.; Wumesh, K.C.; Albring, J.C.; Satpathy, A.T.; Rotondo, J.A.; Edelson, B.T.; Kretzer, N.M.; et al. Compensatory dendritic cell development mediated by BATF-IRF interactions. Nature 2012, 490, 502–507. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Pasquier, J.; Cabau, C.; Nguyen, T.; Jouanno, E.; Severac, D.; Braasch, I.; Journot, L.; Pontarotti, P.; Klopp, C.; Postlethwait, J.H.; et al. Gene evolution and gene expression after whole genome duplication in fish: The PhyloFish database. BMC Genom. 2016, 17, 368. [Google Scholar] [CrossRef] [PubMed]
  24. Leong, J.S.; Jantzen, S.G.; Schalburg, K.R.V.; Cooper, G.A.; Messmer, A.M.; Liao, N.Y.; Munro, S.; Moore, R.; Holt, R.A.; Jones, S.J. Salmo salar and Esox lucius full-length cDNA sequences reveal changes in evolutionary pressures on a post-tetraploidization genome. BMC Genom. 2010, 11, 1–17. [Google Scholar] [CrossRef] [PubMed]
  25. Wang, J.; Liu, M.; Wu, Y.; Yoon, S.; Alnabulsi, A.; Liu, F.; Fernández-Álvarez, C.; Wang, T.; Holland, J.W.; Secombes, C.J.; et al. Immune-modulation of two BATF3 paralogues in rainbow trout Oncorhynchus mykiss. Mol. Immunol. 2018, 99, 104–114. [Google Scholar] [CrossRef] [PubMed]
  26. Chu, P.; He, L.; Zhu, D.; Chen, L.; Huang, R.; Liao, L.; Li, Y.; Zhu, Z.; Wang, Y. Identification, characterisation and preliminary functional analysis of IRAK-M in grass carp (Ctenopharyngodon idella). Fish Shellfish Immunol. 2019, 84, 312–321. [Google Scholar] [CrossRef] [PubMed]
  27. Xiong, L.; He, L.; Luo, L.; Li, Y.; Liao, L.; Huang, R.; Zhu, Z.; Wang, Y. Global and Complement Gene-Specific DNA Methylation in Grass Carp after Grass Carp Reovirus (GCRV) Infection. Int. J. Mol. Sci. 2018, 19, 1110. [Google Scholar] [CrossRef] [PubMed]
  28. Chu, P.; He, L.; Li, Y.; Huang, R.; Liao, L.; Li, Y.; Zhu, Z.; Wang, Y. Molecular cloning and functional characterisation of NLRX1 in grass carp (Ctenopharyngodon idella). Fish Shellfish Immunol. 2018, 81, 276–283. [Google Scholar] [CrossRef]
  29. Wang, Y.; Lu, Y.; Zhang, Y.; Ning, Z.; Li, Y.; Zhao, Q.; Lu, H.; Huang, R.; Xia, X.; Feng, Q.; et al. The draft genome of the grass carp (Ctenopharyngodon idellus) provides insights into its evolution and vegetarian adaptation. Nat. Genet. 2015, 47, 625–631. [Google Scholar] [CrossRef] [PubMed]
  30. Li, G.X.; Zhao, Y.L.; Wang, J.; Liu, B.Z.; Sun, X.L.; Guo, S.; Feng, J.X. Transcriptome profiling of developing spleen tissue and discovery of immune-related genes in grass carp (Ctenopharyngodon idella). Fish Shellfish Immunol. 2017, 60, 400–410. [Google Scholar] [CrossRef]
  31. Chen, J.; Li, C.; Huang, R.; Du, F.; Liao, L.; Zhu, Z.; Wang, Y. Transcriptome analysis of head kidney in grass carp and discovery of immune-related genes. BMC Vet. Res. 2012, 8, 108. [Google Scholar] [CrossRef] [PubMed]
  32. Wei, L.L.; Liu, Y.; Zhong, S.; Wu, H.; Ruan, J.; Liu, M.; Zhou, Q.; Zhong, Q. Transcriptome analysis of grass carp provides insights into the immune-related genes and pathways in response to MC-LR induction. Aquaculture 2018, 488. [Google Scholar] [CrossRef]
  33. He, L.; Zhang, A.; Chu, P.; Li, Y.; Huang, R.; Liao, L.; Zhu, Z.; Wang, Y. Deep Illumina sequencing reveals conserved and novel microRNAs in grass carp in response to grass carp reovirus infection. BMC Genom. 2017, 18, 195. [Google Scholar] [CrossRef]
  34. Muffato, M.; Louis, A.; Poisnel, C.E.; Roest Crollius, H. Genomicus: A database and a browser to study gene synteny in modern and ancestral genomes. Bioinformatics 2010, 26, 1119–1121. [Google Scholar] [CrossRef] [PubMed]
  35. Lollies, A.; Hartmann, S.; Schneider, M.; Bracht, T.; Weiss, A.L.; Arnolds, J.; Klein-Hitpass, L.; Sitek, B.; Hansmann, M.L.; Kuppers, R.; et al. An oncogenic axis of STAT-mediated BATF3 upregulation causing MYC activity in classical Hodgkin lymphoma and anaplastic large cell lymphoma. Leukemia 2018, 32, 92–101. [Google Scholar] [CrossRef]
  36. Grajales-Reyes, G.E.; Iwata, A.; Albring, J.; Wu, X.; Tussiwand, R.; Kc, W.; Kretzer, N.M.; Briseno, C.G.; Durai, V.; Bagadia, P.; et al. Batf3 maintains autoactivation of Irf8 for commitment of a CD8alpha(+) conventional DC clonogenic progenitor. Nat. Immunol. 2015, 16, 708–717. [Google Scholar] [CrossRef] [PubMed]
  37. Schraml, B.U.; Hildner, K.; Ise, W.; Lee, W.L.; Smith, W.A.E.; Solomon, B.; Sahota, G.; Sim, J.; Mukasa, R.; Cemerski, S.; et al. The AP-1 transcription factor Batf controls T(H)17 differentiation. Nature 2009, 460, 405–409. [Google Scholar] [CrossRef] [PubMed]
  38. Ise, W.; Kohyama, M.; Schraml, B.U.; Zhang, T.T.; Schwer, B.; Basu, U.; Alt, F.W.; Tang, J.; Oltz, E.M.; Murphy, T.L.; et al. The transcription factor BATF controls the global regulators of class-switch recombination in both B cells and T cells. Nat. Immunol. 2011, 12, 536–543. [Google Scholar] [CrossRef]
  39. Medzhitov, R.; Preston-Hurlburt, P.; Janeway, C.A., Jr. A human homologue of the Drosophila Toll protein signals activation of adaptive immunity. Nature 1997, 388, 394–397. [Google Scholar] [CrossRef] [PubMed]
  40. Rock, F.L.; Hardiman, G.; Timans, J.C.; Kastelein, R.A.; Bazan, J.F. A family of human receptors structurally related to Drosophila Toll. Proc. Natl. Acad. Sci. USA 1998, 95, 588–593. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Takeuchi, O.; Kawai, T.; Sanjo, H.; Copeland, N.G.; Gilbert, D.J.; Jenkins, N.A.; Takeda, K.; Akira, S. TLR6: A novel member of an expanding toll-like receptor family. Gene 1999, 231, 59–65. [Google Scholar] [CrossRef]
  42. Chuang, T.H.; Ulevitch, R.J. Cloning and characterization of a sub-family of human toll-like receptors: hTLR7, hTLR8 and hTLR9. Eur. Cytokine Netw. 2000, 11, 372–378. [Google Scholar]
  43. Hemmi, H.; Takeuchi, O.; Kawai, T.; Kaisho, T.; Sato, S.; Sanjo, H.; Matsumoto, M.; Hoshino, K.; Wagner, H.; Takeda, K.; et al. A Toll-like receptor recognizes bacterial DNA. Nature 2000, 408, 740–745. [Google Scholar] [CrossRef] [PubMed]
  44. Zhang, D.K.; Zhang, G.L.; Hayden, M.S.; Greenblatt, M.B.; Bussey, C.; Flavell, R.A.; Ghosh, S. A toll-like receptor that prevents infection by uropathogenic bacteria. Science 2004, 303, 1522–1526. [Google Scholar] [CrossRef] [PubMed]
  45. Poltorak, A.; He, X.L.; Smirnova, I.; Liu, M.Y.; Van Huffel, C.; Du, X.; Birdwell, D.; Alejos, E.; Silva, M.; Galanos, C.; et al. Defective LPS signaling in C3H/HeJ and C57BL/10ScCr mice: Mutations in Tlr4 gene. Science 1998, 282, 2085–2088. [Google Scholar] [CrossRef] [PubMed]
  46. Alexopoulou, L.; Holt, A.C.; Medzhitov, R.; Flavell, R.A. Recognition of double-stranded RNA and activation of NF-kappaB by Toll-like receptor 3. Nature 2001, 413, 732–738. [Google Scholar] [CrossRef]
  47. Beg, A.A.; Baltimore, D. An essential role for NF-kappaB in preventing TNF-alpha-induced cell death. Science 1996, 274, 782–784. [Google Scholar] [CrossRef]
  48. Hayden, M.S.; Ghosh, S. Signaling to NF-kappaB. Genes Dev. 2004, 18, 2195–2224. [Google Scholar] [CrossRef]
  49. Liu, Z.G.; Hsu, H.L.; Goeddel, D.V.; Karin, M. Dissection of TNF receptor 1 effector functions: JNK activation is not linked to apoptosis while NF-kappa B activation prevents cell death. Cell 1996, 87, 565–576. [Google Scholar] [CrossRef]
  50. Angel, P.; Karin, M. The Role of Jun, Fos and the Ap-1 Complex in Cell-Proliferation and Transformation. Biochim. Biophys. Acta 1991, 1072, 129–157. [Google Scholar] [CrossRef]
  51. Blank, V. Small Maf proteins in mammalian gene control: Mere dimerization partners or dynamic transcriptional regulators? J. Mol. Biol. 2008, 376, 913–925. [Google Scholar] [CrossRef] [PubMed]
  52. Echlin, D.; Tae, H.N.; Taparowsky, E. B-ATF functions as a negative regulator of AP-1 mediated transcription and blocks cellular transformation by Ras and FOS. Oncogene 2000, 19, 1752–1763. [Google Scholar] [CrossRef] [PubMed]
  53. Williams, K.L.; Nanda, I.; Lyons, G.E.; Kuo, C.T.; Schmid, M.; Leiden, J.M.; Kaplan, M.H.; Taparowsky, E.J. Characterization of murine BATF: A negative regulator of activator protein-1 activity in the thymus. Eur. J. Immunol. 2015, 31, 1620–1627. [Google Scholar] [CrossRef]
  54. Zhao, L.U.; Peng, Z.S.; Jing, N.; Jia, H.T.; Wei, D. BATF2/SARI Induces Tumor Cell Apoptosis by Inhibiting p53-dependent NF-κB Activity. Chin. J. Biochem. Mol. Biol. 2011. [Google Scholar] [CrossRef]
  55. Zhu, D.; Huang, R.; Chen, L.; Fu, P.; Luo, L.; He, L.; Li, Y.; Liao, L.; Zhu, Z.; Wang, Y. Cloning and characterization of the LEF/TCF gene family in grass carp (Ctenopharyngodon idella) and their expression profiles in response to grass carp reovirus infection. Fish Shellfish Immunol. 2019, 86, 335–346. [Google Scholar] [CrossRef] [PubMed]
  56. Tamura, K.; Peterson, D.; Peterson, N.; Stecher, G.; Nei, M.; Kumar, S. MEGA5: Molecular Evolutionary Genetics Analysis Using Maximum Likelihood, Evolutionary Distance, and Maximum Parsimony Methods. Mol. Biol. Evol. 2011, 28, 2731–2739. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(T)(-Delta Delta C) method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
  58. Cui, Y.; Jiang, L.; Xing, R.; Wang, Z.; Wang, Z.; Shao, Y.; Zhang, W.; Zhao, X.; Li, C. Cloning, expression analysis and functional characterization of an interleukin-1 receptor-associated kinase 4 from Apostichopus japonicus. Mol. Immunol. 2018, 101, 479–487. [Google Scholar] [CrossRef]
  59. Saraiva, M.; O’Garra, A. The regulation of IL-10 production by immune cells. Nat. Rev. Immunol. 2010, 10, 170–181. [Google Scholar] [CrossRef] [Green Version]
  60. Han, Y.; Wang, S.; Zhang, Z.; Ma, X.; Li, W.; Zhang, X.; Deng, J.; Wei, H.; Li, Z.; Zhang, X.E.; et al. In vivo imaging of protein-protein and RNA-protein interactions using novel far-red fluorescence complementation systems. Nucleic Acids Res. 2014, 42, e103. [Google Scholar] [CrossRef]
  61. Li, S.; Chen, M.; Xiong, Q.; Zhang, J.; Cui, Z.; Ge, F. Characterization of the Translationally Controlled Tumor Protein (TCTP) Interactome Reveals Novel Binding Partners in Human Cancer Cells. J. Proteome Res. 2016, 15, 3741–3751. [Google Scholar] [CrossRef] [PubMed]
Figure 1. (A) Prediction of domain and gene structure of CiBATF3 by the online website SMART. The lines represent introns; boxes represent exons; and numbers represent the base pair length of each region; and bZIP represent basic leucine zipper domain. (B) Amino acid sequence alignment of CiBATF3 with other species BATF3. The amino acid sequences were aligned using the ClustalW 2.1 program. The core domains were framed with square, including leucine zipper (LZ) regions, hinge (H) and DNA binding domain (DB). Conserved leucine residues are highlighted in gray. Below the alignment, identity (*), strong similarity (:) and weak similarity (.) are represented. (C) The spatial structure of BATF3 protein from human, zebrafish and grass carp were predicted by SWISS-MODEL program. Green, α- helices.
Figure 1. (A) Prediction of domain and gene structure of CiBATF3 by the online website SMART. The lines represent introns; boxes represent exons; and numbers represent the base pair length of each region; and bZIP represent basic leucine zipper domain. (B) Amino acid sequence alignment of CiBATF3 with other species BATF3. The amino acid sequences were aligned using the ClustalW 2.1 program. The core domains were framed with square, including leucine zipper (LZ) regions, hinge (H) and DNA binding domain (DB). Conserved leucine residues are highlighted in gray. Below the alignment, identity (*), strong similarity (:) and weak similarity (.) are represented. (C) The spatial structure of BATF3 protein from human, zebrafish and grass carp were predicted by SWISS-MODEL program. Green, α- helices.
Ijms 20 01687 g001
Figure 2. (A) Syntenic analysis of BATF3 gene in C. idella with different species. Gene synteny relationship of BATF3 genes were obtained by GnomicusV88.01 [34], and then these genes used as references were compared with grass carp genomic sequence. The syntenic analysis of BATF3 genes was performed with neighboring genes from other species. The arrows indicate the direction of transcription. Homologous genes are displayed in the same color. (B) The phylogenetic tree of BATF family genes. The phylogenetic tree of BATF3 was constructed by MEGA v 5.0 program based on neighbor-joining (NJ) method. Abbreviations of the species presented in the figure are as follows: Mus musculus (Mm), Oryzias latipes (Ol), Bos taurus (Bt), Takifugu rubripes (Tr), Gallus gallus (Gg), Xenopus tropicalis (Xt), Danio rerio (Dr), Oncorhynchus mykiss (Om), Homo sapiens (Hs), Poecilia formosa (Pf), Salmo salar (Ss), Latimeria chalumnae (Lc), Lepisosteus oculatus (Lo), Oreochromis niloticus (On), Esox lucius (El), Ictalurus punctatus (Ip), Cynoglossus semilaevis (Cs). The accession numbers are given after the species names in the tree.
Figure 2. (A) Syntenic analysis of BATF3 gene in C. idella with different species. Gene synteny relationship of BATF3 genes were obtained by GnomicusV88.01 [34], and then these genes used as references were compared with grass carp genomic sequence. The syntenic analysis of BATF3 genes was performed with neighboring genes from other species. The arrows indicate the direction of transcription. Homologous genes are displayed in the same color. (B) The phylogenetic tree of BATF family genes. The phylogenetic tree of BATF3 was constructed by MEGA v 5.0 program based on neighbor-joining (NJ) method. Abbreviations of the species presented in the figure are as follows: Mus musculus (Mm), Oryzias latipes (Ol), Bos taurus (Bt), Takifugu rubripes (Tr), Gallus gallus (Gg), Xenopus tropicalis (Xt), Danio rerio (Dr), Oncorhynchus mykiss (Om), Homo sapiens (Hs), Poecilia formosa (Pf), Salmo salar (Ss), Latimeria chalumnae (Lc), Lepisosteus oculatus (Lo), Oreochromis niloticus (On), Esox lucius (El), Ictalurus punctatus (Ip), Cynoglossus semilaevis (Cs). The accession numbers are given after the species names in the tree.
Ijms 20 01687 g002
Figure 3. (A) Tissue distribution of CiBATF3; (BG) the expression patterns of CiBATF3 in immune-related tissues after GCRV challenge; and (H) the expression of the S6 gene in gill and intestine. The expression level of the control group was set as the baseline (1.0). The results are presented as the mean ± SD (n = 5). The abbreviations of 10 examined tissues. L, liver; S, spleen; G, gill; HK, head kidney; B, brain; MK, middle kidney; H, heart; M, muscle; SK, skin; I, intestine. * = p ≤ 0.05, ** = p ≤ 0.01.
Figure 3. (A) Tissue distribution of CiBATF3; (BG) the expression patterns of CiBATF3 in immune-related tissues after GCRV challenge; and (H) the expression of the S6 gene in gill and intestine. The expression level of the control group was set as the baseline (1.0). The results are presented as the mean ± SD (n = 5). The abbreviations of 10 examined tissues. L, liver; S, spleen; G, gill; HK, head kidney; B, brain; MK, middle kidney; H, heart; M, muscle; SK, skin; I, intestine. * = p ≤ 0.05, ** = p ≤ 0.01.
Ijms 20 01687 g003
Figure 4. The expression levels of IRF8 (AD) and IL-10 (EH) in four immune tissues (spleen, intestine, head kidney and middle kidney) after GCRV infection. The results are presented as the mean ± SD (n = 5). The expression level of the control group was set as the baseline (1.0) and were normalized to the expression in the untreated groups (Day 0). * = p ≤ 0.05.
Figure 4. The expression levels of IRF8 (AD) and IL-10 (EH) in four immune tissues (spleen, intestine, head kidney and middle kidney) after GCRV infection. The results are presented as the mean ± SD (n = 5). The expression level of the control group was set as the baseline (1.0) and were normalized to the expression in the untreated groups (Day 0). * = p ≤ 0.05.
Ijms 20 01687 g004
Figure 5. The expression patterns of CiBATF3 in CIK cells after LPS (A) or poly(I:C) (B) challenge. CIK cells were stimulated with poly(I:C) (20 μg mL−1) and LPS (10 μg mL−1), respectively. Untreated cells served as controls. The expression level of the control group was set as the baseline (1.0) and the results are presented as the mean ± SD (n = 3). * = p ≤ 0.05, ** = p ≤ 0.01.
Figure 5. The expression patterns of CiBATF3 in CIK cells after LPS (A) or poly(I:C) (B) challenge. CIK cells were stimulated with poly(I:C) (20 μg mL−1) and LPS (10 μg mL−1), respectively. Untreated cells served as controls. The expression level of the control group was set as the baseline (1.0) and the results are presented as the mean ± SD (n = 3). * = p ≤ 0.05, ** = p ≤ 0.01.
Ijms 20 01687 g005
Figure 6. Gene expression profiles in CiBATF3 deficient CIK cells generated by siRNA. (A) Three siRNAs targeted to different sites of CiBATF3 were used to infected CIK cells, and siRNA-BATF3-2 has the highest efficiency in CiBATF3 gene silencing. (B) The expression of CiBATF3, GATA3, IL-4, IL-10 IRF8, IL-12 p35, IL-12 p40 and c-myc at 24 h were assayed by RT-qPCR after infected with siRNA-BATF3-2. The expression level of the control group was set as the baseline (1.0) and the results are presented as the mean ± SD (n = 3). siRNA-scramble: control siRNA infection; siRNA-BATF3: CiBATF3 siRNA infection. * = p ≤ 0.05, ** = p ≤ 0.01.
Figure 6. Gene expression profiles in CiBATF3 deficient CIK cells generated by siRNA. (A) Three siRNAs targeted to different sites of CiBATF3 were used to infected CIK cells, and siRNA-BATF3-2 has the highest efficiency in CiBATF3 gene silencing. (B) The expression of CiBATF3, GATA3, IL-4, IL-10 IRF8, IL-12 p35, IL-12 p40 and c-myc at 24 h were assayed by RT-qPCR after infected with siRNA-BATF3-2. The expression level of the control group was set as the baseline (1.0) and the results are presented as the mean ± SD (n = 3). siRNA-scramble: control siRNA infection; siRNA-BATF3: CiBATF3 siRNA infection. * = p ≤ 0.05, ** = p ≤ 0.01.
Ijms 20 01687 g006
Figure 7. Gene expression profiling in head kidney (HK) and middle kidney (MK) at 24 h after injected siRNA-CiBATF3 in vivo. The expression of CiBATF3 (A), GATA3 (B), IL-4 (C), IL-10 (D), IRF8 (E), IL-12 p35 (F), IL-12 p40 (G) and c-myc (H) at 24 h were assayed by RT-qPCR. The expression level of the control group was set as the baseline (1.0) and the results are presented as the mean ± SD (n = 3). siRNA-scramble: control siRNA infection; siBATF3: CiBATF3 siRNA infection. * = p ≤ 0.05, ** = p ≤ 0.01.
Figure 7. Gene expression profiling in head kidney (HK) and middle kidney (MK) at 24 h after injected siRNA-CiBATF3 in vivo. The expression of CiBATF3 (A), GATA3 (B), IL-4 (C), IL-10 (D), IRF8 (E), IL-12 p35 (F), IL-12 p40 (G) and c-myc (H) at 24 h were assayed by RT-qPCR. The expression level of the control group was set as the baseline (1.0) and the results are presented as the mean ± SD (n = 3). siRNA-scramble: control siRNA infection; siBATF3: CiBATF3 siRNA infection. * = p ≤ 0.05, ** = p ≤ 0.01.
Ijms 20 01687 g007
Figure 8. Impact of CiBATF3 on AP-1 and NF-κB reporter gene. HEK293 cells were transfected with 0.25 μg pCDNA3.1-CiBATF3, 0.25 μg reporter gene vectors (pNFқB-TA-luc and pAP1-TA-luc, respectively), 0.025 μg pRL-TK renilla luciferase vector and 1 μL Lipo6000™ Transfection Reagent. At 24 h post-transfection, the cells were collected and lysed. Luciferase activities were measured using the Luciferase Assay System (Promega). The results are presented as the mean ± SD (n = 3). Differences were considered significant at p < 0.01.
Figure 8. Impact of CiBATF3 on AP-1 and NF-κB reporter gene. HEK293 cells were transfected with 0.25 μg pCDNA3.1-CiBATF3, 0.25 μg reporter gene vectors (pNFқB-TA-luc and pAP1-TA-luc, respectively), 0.025 μg pRL-TK renilla luciferase vector and 1 μL Lipo6000™ Transfection Reagent. At 24 h post-transfection, the cells were collected and lysed. Luciferase activities were measured using the Luciferase Assay System (Promega). The results are presented as the mean ± SD (n = 3). Differences were considered significant at p < 0.01.
Ijms 20 01687 g008
Figure 9. Subcellular localization of CiBATF3 in: CIK cells (A); and HEK293 cells (B). (C) CiBATF3-GFP nuclear translocation induced by PAMPs stimulation. Green fluorescence shows the distribution of EGFP or EGFP-tagged proteins and blue fluorescence shows the nucleus that was stained by Hoechst 33342 under a 60× oil immersion objective lens (scale bars, 10 μm and 25 μm). (D) Imaging of the protein–protein interaction by using far-red mNeptune-based BiFC in HEK293 cells. The fluorescence of the mNeptune channel was red. The images were acquired using fluorescence microscopy and a 40× oil immersion objective lens (scale bar, 100 μm).
Figure 9. Subcellular localization of CiBATF3 in: CIK cells (A); and HEK293 cells (B). (C) CiBATF3-GFP nuclear translocation induced by PAMPs stimulation. Green fluorescence shows the distribution of EGFP or EGFP-tagged proteins and blue fluorescence shows the nucleus that was stained by Hoechst 33342 under a 60× oil immersion objective lens (scale bars, 10 μm and 25 μm). (D) Imaging of the protein–protein interaction by using far-red mNeptune-based BiFC in HEK293 cells. The fluorescence of the mNeptune channel was red. The images were acquired using fluorescence microscopy and a 40× oil immersion objective lens (scale bar, 100 μm).
Ijms 20 01687 g009
Table 1. Amino acid sequence similarity of CiBATF3 with other species by Geneious.
Table 1. Amino acid sequence similarity of CiBATF3 with other species by Geneious.
The Species Name and Accession NumbersPercent of Identity
Percent of identity123456789101112131415161718
1. CiBATF38569.850.850.851.653.6505044.548.741.744.143.336.940.339.746.9
2. DrBATF3_NP_001038857.18567.25252.85253.151.251.244.547.940.242.444.136.941.140.545.3
3. IpBATF3_XP_017310165.169.867.249.25053.454.651.651.65048.745.14341.537.341.642.147.2
4. OmBATF3a_APH08552.150.85249.298.489.480.5696948.853.741.741.733.931.532.63142.2
5. OmBATF3_XP_021454099.150.852.85098.489.480.5696947.952.940.94034.632.333.331.742.2
6. SsBATF3_XP_014060214.151.65253.489.489.488.573.673.649.654.544.142.53735.436.435.746.9
7. SsBATF3_NP_001135082.153.653.154.680.580.588.569.768.951.452.344.443.136.837.537.837.447.5
8. OmBATF3b_APH08553.15051.251.6696973.669.796.95056.346.641.639.835.937.736.248.1
9. SsBATF3_XP_014000176.15051.251.6696973.668.996.949.254.8454039.136.636.235.448.1
10. PfBATF3_XP_007574812.144.544.55048.847.949.651.45049.270.661.35338.134.437.538.738.9
11. AoBATF3_XP_023123541.148.747.948.753.752.954.552.356.354.870.66354.234.933.338.339.240.2
12. OlBATF3_XP_011472010.141.740.245.141.740.944.144.446.64561.3635235.132.840.63842
13. CsBATF3_XP_016889000.144.142.44341.74042.543.141.6405354.25235.230.537.635.541.8
14. XtBATF3_XP_004914833.143.344.141.533.934.63736.839.839.138.134.935.135.261.561.256.653.1
15. GgBATF3_XP_419428.436.936.937.331.532.335.437.535.936.634.433.332.830.561.561.557.754.1
16. HsBATF3_NP_061134.140.341.141.632.633.336.437.837.736.237.538.340.637.661.261.580.357.5
17. MmBATF3_NP_084336.139.740.542.13131.735.737.436.235.438.739.23835.556.657.780.355.2
18. LoBATF3_XP_015194606.146.945.347.242.242.246.947.548.148.138.940.24241.853.154.157.555.2
Abbreviations of the species presented in the figure are as follows: Ctenopharyngodon idella (Ci), Amphiprion ocellaris (Ao), Mus musculus (Mm), Oryzias latipes (Ol), Bos taurus (Bt), Gallus gallus (Gg), Xenopus tropicalis (Xt), Danio rerio (Dr), Oncorhynchus mykiss (Om), Homo sapiens (Hs), Poecilia formosa (Pf), Salmo salar (Ss), Lepisosteus oculatus (Lo), Ictalurus punctatus (Ip), Cynoglossus semilaevis (Cs). The accession numbers are given after the species names in the table.
Table 2. Primers used for cloning and RT-qPCR and the siRNA sequence information.
Table 2. Primers used for cloning and RT-qPCR and the siRNA sequence information.
PrimersSequences (5′—3′)Purpose
BATF3-5′RoutGCACTGTGCCTGTGGACCTT5′ RACE
BATF3-5′RinACGCCTCGTGCAACTCGTCA
BATF3-3′RoutAGTGATGCTCCAGCTTTACGGT3′ RACE
BATF3-3′RinACCGAGTTGCTGCCCAGAGA
BATF3-FATGTCACTTTTCAATGCGACAAGTAAcDNA cloning
BATF3-RTCAGATGTGAATGTCTTGTGGCACTGT
qBATF3-FAGTGATGCTCCAGCTTTACGGTRT-qPCR
qBATF3-RACGCCTCGTGCAACTCGTCA
qIL10-FTATTAAACGAGAACGTGCAACAGAA
qIL10-RTCCCGCTTGAGATCTTGAAATATACT
qIRF8-FCAGAGGAGGAACAGAAGTTGGGTAA
qIRF8-RACGCTTCAGGATGCCCATGTA
q β-actin-FTCGGTATGGGACAGAAGGAC
q β-actin-RGACCAGAGGCATACAGGGAC
qIL4-FCTCAGGTGAAGCCCTTTGCC
qIL4-RACTGGATGTTCCTCTGAAGCTGTAA
qIL-12 p35-FAGGCTCGGATGATTCCTTACA
qIL-12 p35-RTCACACTGGGCTGGTAGGAG
qIL-12 p40-FGGAGAAGTCTACGAAGGGCAA
qIL-12 p40-FGTGTGTGGTTTAGGTAGGAGCC
q c-myc-FGAGCGAAGACATTTGGAAGA
q c-myc-RTGATGAAGGACTGGGAGTAG
qGATA-3-FTACGAGGAGGACAAAGAGT
qGATA-3-RGTAAGTGGCGATGGGATGGT
qS6-FAGCGCAGCAGGCAATTACTATCTGCRV RT-qPCR
qS6-RATCTGCTGGTAATGCGGAACG
Negative controlUUCUCCGAACGUGUCACGUTTsiRNA
siRNA-BATF3-1GGCGUAUGAGUGUCUGGAATT
siRNA-BATF3-2GAGGAACAGCAACGCUUAATT
siRNA-BATF3-3GCACAGUGCCACAAGACAUTT
Table 3. Primers used in plasmid construction. The different restriction sites are shown in lowercase letters and the names are represented above the respective nucleotides.
Table 3. Primers used in plasmid construction. The different restriction sites are shown in lowercase letters and the names are represented above the respective nucleotides.
PrimersSequences (5′–3′)Purpose
L-BATF3-F   XhoI
CGctcgagGTATGTCACTTTTCAATGCGACAAGTAA
Subcellular localization
L-BATF3-R   BamHI
CGggatccGATGTGAATGTCTTGTGGCACTGT
BF-BATF3-F    XhoI
CCCctcgagCTATGTCACTTTTCAATGCGACAAGTAA
BiFC Analysis
BF-BATF3-R    EcoRI
CAAgaattcGAACTCCCGCCACCTCCACTCCCG
CCACCTCCGATGTGAATGTCTTGTGGCACTGT
BF-IL10-F   EcoRI
CGgaattcTGGGAGGTGGCGGGAGTGGAGGTGGCGG
GAGTATGATTTTCTCTAGAGTCATCTTTTCTGC
BF-IL10-R    KpnI
CGCggtaccGTGCTTTTCTCTCTTTGATGCCAG

Share and Cite

MDPI and ACS Style

Zhu, D.; Huang, R.; Fu, P.; Chen, L.; Luo, L.; Chu, P.; He, L.; Li, Y.; Liao, L.; Zhu, Z.; et al. Investigating the Role of BATF3 in Grass Carp (Ctenopharyngodon idella) Immune Modulation: A Fundamental Functional Analysis. Int. J. Mol. Sci. 2019, 20, 1687. https://doi.org/10.3390/ijms20071687

AMA Style

Zhu D, Huang R, Fu P, Chen L, Luo L, Chu P, He L, Li Y, Liao L, Zhu Z, et al. Investigating the Role of BATF3 in Grass Carp (Ctenopharyngodon idella) Immune Modulation: A Fundamental Functional Analysis. International Journal of Molecular Sciences. 2019; 20(7):1687. https://doi.org/10.3390/ijms20071687

Chicago/Turabian Style

Zhu, Denghui, Rong Huang, Peipei Fu, Liangming Chen, Lifei Luo, Pengfei Chu, Libo He, Yongming Li, Lanjie Liao, Zuoyan Zhu, and et al. 2019. "Investigating the Role of BATF3 in Grass Carp (Ctenopharyngodon idella) Immune Modulation: A Fundamental Functional Analysis" International Journal of Molecular Sciences 20, no. 7: 1687. https://doi.org/10.3390/ijms20071687

APA Style

Zhu, D., Huang, R., Fu, P., Chen, L., Luo, L., Chu, P., He, L., Li, Y., Liao, L., Zhu, Z., & Wang, Y. (2019). Investigating the Role of BATF3 in Grass Carp (Ctenopharyngodon idella) Immune Modulation: A Fundamental Functional Analysis. International Journal of Molecular Sciences, 20(7), 1687. https://doi.org/10.3390/ijms20071687

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop