Characterization of porcine cytokine inducible SH2-containing protein gene and its association with piglet diarrhea traits.
The inflammation cytokine response is an important feature of the host's immune response to infection caused by diverse pathogens, but excessive inflammatory response will be harmful to host. The control of inflammatory response is negatively mediated in part from the suppressor effect of cytokine signaling (SOCS) proteins through their actions on specific signaling pathways such as Janus kinase 2 (JAK2) signal transducer and activator of transcription 5 (STAT5) and nuclear factor-k-gene binding pathways [1-2]. Cytokine inducible SH2-containing protein (CISH) was the first member of the SOCS family to be described. It was initially isolated as an immediate early gene induced by granulocyte-macrophage colony stimulating factor, interleukin (IL)-2, IL-3 and erythropoietin in hematopoietic cells . Then, CISH was identified as a target gene of the JAK2-STAT5 pathway . CISH negatively modulates STAT5 activation by binding to the phosphorylated tyrosine residues of cytokine receptors and thereby inhibits the downstream cytokine signaling [4-6]. In addition, accumulating evidence supports the important regulatory role of CISH in the T-cell proliferation and survival in response to infection [7-9], in the type 1 dendritic cell (DC) development and DC-mediated cytotoxic T lymphocyte activation .
Given the central role in controlling signaling of cytokines, in particular IL-2 , CISH was selected as a candidate gene of human disease. Recent reports showed that the variants of human CISH were associated with susceptibility to diseases caused by diverse infectious pathogens, such as bacteremia , malaria , tuberculosis [12-13], and persistent hepatitis B .
In modern pig production, the large-scale intensive systems make piglets susceptible to infection by various pathogens . Given the significant effects of CISH gene on human disease, in this research, we characterized the CISH gene in pigs, determined its subcellular localization, investigated genetic variation and performed association analysis with piglet diarrhea traits to evaluate its effects on pig breeding.
MATERIALS AND METHODS
Animals and data collection
In the autumn of 2013, 327 Min piglets and 203 Landrance piglets born in the same period were selected as experiment animals on the Lanxi conservation farm (Lanxi, Heilongjiang, China). All these piglets were raised with the same fodder, feeding conditions, and disease control during the experiment period. From birth to 35 days, piglet diarrhea was recorded twice per day (6:00 am and 2:00 pm) and was assigned a daily score based on a standard evaluation of symptom traits as follows: 0, normal, manure solid feces; 1, slight diarrhea, manure soft and loose feces; 2, moderate diarrhea, manure semi-liquid feces; and 3, severe diarrhea manure liquid and unformed feces . Diarrhea index was the summation of diarrhea score during experiment period.
Genomic DNA was extracted from the ear tissues samples using DNA Extraction Kit (Sangon, Shanghai, China) and stored at -20[degrees]C.
Standard polymerase chain reaction reactions and sequencing
Primer pairs were designed using Primer 5.0 and then synthesized commercially (Sangon, China). The polymerase chain reaction (PCR) reaction was 25 [micro]L containing 100 ng templates, 0.5 [micro]M of each primer, 0.5 mM dNTPs, 2.0 U Taq DNA polymerase and 2.5 [micro]L buffer (TaKaRa, Tokyo, Japan). PCR conditions as follows: 5 min at 94[degrees]C, 35 cycles of 45 s at 94[degrees]C, 45 s at the annealing temperature (Table 1), then 90 s at 72[degrees]C, and finally 72[degrees]C for 10 min. The PCR products were analyzed by 1.5% agarose gel electrophoresis, purified using a Gel Extraction Kit (Tiangen Biotech, Beijing, China), cloned into the pMD18-T vector (TaKaRa, Japan) and sequenced commercially (Sangon, China).
Isolation and characterization of the cytokine inducible SH2-containing protein
Total RNA of spleen tissues was isolated using Trizol Reagent (Invitrogen, Carlsbad, CA, USA). Reverse transcription (RT) was performed using the TransSript First-Strand cDNA Synthesis SuperMix (TaKaRa, Japan). Porcine expressed sequence tags (ESTs) sequences (GenBank: BE014034, DY425421, DY421524) homologous to the human orthologous genes (GenBank: NM_ 145071) were selected and two ESTs (DY425421 and DY421524) were assembled into a contig sequence. Primer pairs C-gap (Table 1) was used to amplify the gap sequence between EST (BE014034) and the contig. Three pairs of primers (Table 1: C-12, C-23, and C-33) were used to amplify the introns and exons of cytokine inducible SH2-containing protein (pCISH). Program Seqman (DNAstar, Madison, WI, USA) was used to identify the open reading frame (ORF) and deduce the amino acid sequences. The Sequences of CISH from different species were aligned using the Clustalw 2.0 (www.ebi.ac.uk/Tools/clustalw2/). Then the phylogenetic trees of this gene were reconstructed by the neighborjoining (NJ) method of MEGA v4.1 program and the bootstrap tests were performed with 1,000 bootstrap replications.
Intracellular distribution of the pCISH
To obtain the coding sequence of pCISH gene, three pairs of primers (CE-1, CE-2, and CE-3) (Table 1) were used in the amplification and the final products were cloned into pMD18-T vectors (TaKaRa, Japan). Then, the correct insert was subcloned into the BglII-SalI site of the pEGFP-C1 vector (BD Biosciences Clontech, San Francisco, CA, USA) to construct a pEGFP-C1-CISH vector which was verified by sequencing and restriction endonuclease digestions.
PK15 cells were cultured using DMEM supplemented with 10% fetal bovine serum (Gibco, Grand Island, NY, USA), 4 mM glutamine, 100 U/mL penicillin, and 0.1 mg/mL streptomycin under air containing 5% C[O.sub.2] at 37[degrees]C. Then these cells were seeded at a density of 5 x [10.sup.5]/mL and transfected with pEGFP-C1-CISH constructs (2 [micro]g of plasmid DNA) using Fugene HD transfection reagent (Roche, Mannheim, Germany). After 20 to 24 hours, the cells were fixed with growth medium containing 3.7% formaldehyde, incubated with 10 [micro]M Hoechst33258. At last, the slides were mounted, sealed and analyzed by LeicaDMI 6000B Confocal Microscope (Leica, Solms, Germany).
Polymorphism identification and genotyping with polymerase chain reaction-restriction fragment length polymorphism
Genetic variation of the pCISH were analyzed by sequencing PCR products amplified by primer pairs CV-1, CV-2, CV-3, CV-4, and CV-5 (Table 1), and all these primers were partially overlapped. The PCR products from Minzhu, Landrace, and Largewhite pig breeds were purified, sequenced and then compared with each other using Clustalw 2.0. Two PCR-restriction fragment length polymorphism (RFLP) assays were designed to genotype the selected single nucleotide polymorphisms (SNPs) using primer pairs CV-1, CV-41 (Table 1) and restriction enzymes TaqI (TaKaRa, Japan), HaeIII (TaKaRa, Japan), respectively. For the PCR-RFLP profile, 8.5 [micro]L of PCR products amplified by primer pairs CV-1 or CV-41 were digested with 5 U of TaqI or HaeIII (TaKaRa, Japan) for 5 h at 65[degrees]C or 37[degrees]C, and then separated on a 1.5% agarose gel. Genotyping work was performed in both Min pigs and Landrace population for further association analysis between pCISH and diarrhea index trait.
In this study, we used diarrhea index to evaluate the piglet diarrhea. It was analyzed as discrete data using the SAS GENMOD procedure. In this model, the dependent variable is diarrhea index, divided into 23 classes according to the summation of diarrhea score during experiment period. The independent variables are breed and genotype effect.
Isolation and characterization of pCISH
A 503 bp sequence from porcine spleen cDNA was obtained using primer C-gap. Together with porcine ESTs (BE014034, DY425421, DY421524), all these fragments were assembled into a 1,913 bp contig cDNA sequence of pCISH. The ORF of pCISH is 765 bp (GenBank: KM507873), which encodes a protein of 254 AA. Homology comparison showed that pCISH is similar to human, mouse, rat and cattle CISH (90% to 93%). Phylogenetic tree analysis indicated that pig and cattle homologs fell into the same cluster (Figure 1). Simple modular architecture research tool (SMART, http://smart.embl-heidelberg.de/) analysis showed that pCISH contained a SH2 domain (80-169AA), a SOCS domain (210-251AA) and a SOCS_box domain (216-250AA).
The 3,476 bp, 806 bp, and 1,004 bp fragments were obtained from genomic DNA using primer C-12, C-23, and C-33, respectively. Then, these three sequences and the 1,913 bp cDNA sequence were assembled into a 5,409 bp genomic sequence of pCISH gene which was consisted of three exons and two introns. Further analysis found that the ATG start codon was located in exon1 and the exon 3 contains the TGA stop codon. In addition, by alignment of this 5,409 bp fragment with the porcine genomic sequence (Sscrofa10.2) using BLAT (http://www.ensembl.org), the pCISH gene was assigned to porcine chromosome 13.
Intracellular distribution of the pCISH
The cellular location of pCISH was determined by fluorescence and confocal analysis of PK15 cells transiently transfected with pEGFP-C1-CISH. After labeling nuclei by staining with Hoechst 33342, CISH-EGFP fusion proteins were found in distribution throughout the cytoplasm (Figure 2). This result was consistent with the previously information of primarily cytoplasmic localization of rat CISH in HEK 293t (ATCC) cells .
Genetic variation analysis of pCISH
Comparative sequence analysis revealed more than 20 polymorphisms among the pCISH genomic sequence in three pig breeds. SNP A-1575G in promoter region and SNP A2497C in Intron1 could be detected by Taq I and HaeIII PCR-RFLP, respectively, and these SNPs wer selected for further association analysis. For the Taq I PCR-RFLP assay, the 1,554 bp PCR product was digested into one fragments of the A allele (1,554 bp) and two fragments of the G allele (1,042 and 512 bp) (Figure 3). For the HaeIII PCR-RFLP, the 192 bp PCR product was digested into one fragments of the A allele (192 bp) and two fragments of the C allele (101 and 91 bp) (Figure 4).
Gen Mod preteduie analysis showed that SNP A-1575G was significantly associated with diarrhea index of Min piglets (p< 0.05) (Table 2), while SNP A2497C was significantly associated with the diarrhea trait of both Min pig and Landrace piglets (p< 0.05) (Table 2).
Our present study provides the basic structure of pCISH. Comparison with the gene structure of human and rat CISH reveals remarkable similarity among these three species, and they both contained 3 Exons and 2 Introns. The inferred amino acid sequence of porcine CISH showed 93% and 91% identity with that of human and eat respeytivuly. Fimrther analysis showed thac pCISH contamed a SH2 domain (80-169AA) and a SOCS domain (210-251AA). The SH2 domain and C-terminal SOCS box are characteristic of the SOCS protein family. In the negative regulation of several cytokine pathways, SOCS proteins bind to tyrosine-phosphorylated residues of target proteins via their SH2 domains, then inhibit JAK activity through their N-terminal domains, and are thought to induce the degradation of bound molecules through a conserved SOCS-box motif that interacts with the proteasome . Hence, pCISH might play a role in the cytokine signal transduction of pigs. In the phylogenetic tree of CISH, pig is closer to cattle than other mammals in an evolutionary relationship. As far as we know, research on cattle CISH is extremely scarce, so this present study indicates that it might be interesting to determine the importance of CISH gene in disease resistance of cattle.
In humans, many polymorphisms in CISH gene were observed and some were strongly associated with the susceptibility to diseases caused by diverse infectious pathogens [12-14]. In this study, firstly, we obtained the genomic sequence of pCISH, but observed only two SNPs which could be detected by PCR-RFLP. Statistical analysis indicated that SNP A-1575G was significantly associated with Min piglets diarrhea index (p<0.05), while SNP A2497C was significantly associated with the diarrhea trait of both Min pig and Landrace piglets (p<0.05).
Two important considerations in candidate gene selection are function and position . To assess the biological function yf jaorcine SOCS gene family, Delgado-Ortega et al  evaluated die expression pattern of SOCS genes in 10 Selected tissues (kidney, large intestine, liver, lung, mesenteric lymph node, small intestine, spleen, stomach, trachea, and thymus) of twenty two-months-old healthy piglets and found that pCISH gene was highly expressed in large intestine, small intestine, spleen and thymus, suggesting that this gene might play a role in porcine intestine immune responses. In this research, the SNP A-1575G and SNP A2497C resided in the 5' upstream region and Intron1, respectively. Although the SNP in these regions did not directly alter any amino acid residue, they might play a role in regulating gene expression and thus their constituent SNPs might be directly related to functional variation. For example, a SNP (G3072A) in intron3 of insulin-like growth factor 2 (IGF2) gene was significantly associated with muscle growth in the pig by regulating the expression of IGF2 mRNA in postnatal muscle . So, in further study, we intend to reveal whether these mutations have effects on the transcription of pCISH gene. To account for the effect of position, we compared the position of pCISH gene and the quantitative trait loci (QTL) for anti-disease or immunity traits. The pCISH gene was located within 27 to 28 MB of porcine chromosome 13 (SSC13) by comparing sequence obtained in this study from the database of porcine genome sequence (Sscrofa10). Interestingly, when browsing QTLs in the latest released Pig QTL Database, we found two QTL for enterotoxigenic Escherichia coli (ETEC) susceptibility on SSC13 , but the pCISH gene does not reside in the region of these QTLs. So, these SNPs might be potential markers for piglet diarrhea traits, even if they do not alter the amino acid residues or transcription factor binding sites of the pCISH gene, suggesting these mutations affect piglet diarrhea traits indirectly, being in linkage disequilibrium with QTLs or another causative polymorphisms affecting piglet diarrhea traits.
In conclusion, we have shown that pCISH was highly conserved at the level of gene structure. Phylogenetic analysis revealed that based on the protein sequence of CISH, pig and cattle were closer than other mammals in evolutionary relationship. Our data also indicated that pCISH was mainly localized in cytoplasm in the PK-15 cells. Furthermore, association analysis revealed that these two detected SNPs significantly influenced piglet diarrhea, suggesting that pCISH gene might be a novel candidate gene for piglet diarrhea traits. However, our results were based on a limited number and breeds of animals, further studies in more pig populations with a larger sample size or different breeds are needed to confirm these preliminary researches.
CONFLICT OF INTEREST
We certify that there is no conflict of interest with any financial organization regarding the material discussed in the manuscript.
This study was supported financially by National Natural Science Foundation of China (31301935), Ph.D. Programs Foundation of Ministry of Education of China (20112325120013).
[1.] Akhtar LN, Benveniste EN. Viral exploitation of host SOCS protein functions. J Virol 2011;85:1912-21.
[2.] Yoshimura A, Naka T, Kubo M. SOCS proteins, cytokine signaling and immune regulation. Nat Rev Immunol 2007;7:454-65.
[3.] Yoshimura A, Ohkubo T, Kiguchi T, et al. A novel cytokine-inducible gene CIS encodes an SH2-containing protein that binds to tyrosinephosphorylated interleukin3 and erythropoietin receptors. EMBO J 1995;14:2816-26.
[4.] Matsumoto A, Masuhara M, Mitsui K, et al CIS, a cytokine inducible SH2 protein, is a target of the JAK-STAT5 pathway and modulates STAT5 activation. Blood 1997;89:3148-54.
[5.] Nakajima H, Liu XW, Wynshaw-Boris A, et al An indirect effect of Stat5a in IL-2- induced proliferation: a critical role for Stat5a in IL-2- mediated IL-2 receptor alpha chain induction. Immunity 1997;7: 691-701.
[6.] Moriggl R, Topham DJ, Teglund S, et al. Stat5 is required for IL-2-induced cell cycle progression of peripheral T cells. Immunity 1999; 10:249-59.
[7.] Kovanen PE, Leonard WJ. Cytokines and immunodeficiency diseases: critical roles of the gamma(c)-dependent cytokines interleukins 2, 4, 7, 9, 15, and 21, and their signaling pathways. Immunol Rev 2004; 202:67-83.
[8.] Li S, Chen S, Xu X, et al. Cytokine-induced Src homology 2 protein (CIS) promotes T cell receptor-mediated proliferation and prolongs survival of activated T cells. J Exp Med 2000;191: 985-94.
[9.] Periasamy S, Dhiman R, Barnes PF, et al. Programmed death 1 and cytokine inducible SH2-containing protein dependent expansion of regulatory T cells upon stimulation with Mycobacterium tuberculosis. J Infect Dis 2011;203:1256-63.
[10.] Miah MA, Yoon CH, Kim J, et al. CISH is induced during DC development and regulates DC-mediated CTL activation. Eur J Immunol 2012;42:58-68.
[11.] Jin P, Wang E, Provenzano M, et al. Molecular signatures induced by interleukin-2 on peripheral blood mononuclear cells and T cell subsets. J Transl Med 2006;4:26.
[12.] Khor CC, Vannberg FO, Chapman SJ, et al. CISH and susceptibility to infectious diseases. N Engl J Med 2010;362:2092-101.
[13.] Sun L, Jin YQ, Shen C, et al. Genetic contribution of CISH promoter polymorphisms to susceptibility to tuberculosis in Chinese children. PLoS One 2014;9:e92020.
[14.] Hu Z, Yang J, Wu Y, et al. Polymorphisms in CISH gene are associated with persistent hepatitis B virus infection in Han Chinese population. Plos One 2014;9:e100826.
[15.] Lunney JK, Ho CS, Wysocki M, et al. Molecular genetics of the swine major histocompatibility complex, the SLA complex. Dev Comp Immunol 2009;33:362-74.
[16.] Kelly D, O'Brien JJ, McCracken KJ. Effect of creep feeding on the incidence, duration and severity of post-weaning diarrhea in pigs. Res Vet Sci 1990;49:223-8.
[17.] Jensik PJ, Arbogast LA. Regulation of cytokine-inducible SH2-containing protein (CIS) by ubiquitination and Elongin B/C interaction. Mol Cell Endocrinol 2015;401:130-41.
[18.] Nicola NA, Greenhalgh CJ. The suppressors of cytokine signaling (SOCS) proteins: important feedback inhibitors of cytokine action. Exp Hematol 2000;28:1105-12.
[19.] Onteru SK, Ross JW, Rothschild MF. The role of gene discovery, QTL analyses and gene expression in reproductive traits in the pig. Soc Reprod Fertil Suppl 2009;66: 87-102.
[20.] Delgado-Ortega M, Melo S, Meurens F. Expression of SOCS1-7 and CIS mRNA in porcine tissues. Vet Immunol Immunopathol 2011; 144:493-8.
[21.] Van Laere AS, Nguyen M, Braunschweig M, et al A regulatory muta tion in IGF2 causes a major QTL effect on muscle growth in the pig. Nature 2003;425:832-6.
[22.] Ji H, Ren J, Yan X, et al. The porcine MUC20 gene: molecular characterization and its association with susceptibility to enterotoxigenic Escherichia coli F4ab/ac. Mol Biol Rep 2011;38:1593-60.
Buyue Niu (1), Dongchun Guo (2), Zhiran Liu (1), Xiaofei Han (1), and Xibiao Wang (1), *
* Corresponding Author: Xibiao Wang
Tel: +86-451-55191470, Fax: +86-451-55191470, E-mail: firstname.lastname@example.org
(1) College of Animal Science and Technology, Northeast Agricultural University, Harbin 150030, China
(2) State Key Laboratory of Veterinary Biotechnology, Harbin Veterinary Research Institute, Chinese Academy of Agricultural Sciences (CAAS), Harbin 150001, China
Submitted Mar 1, 2016; Revised May 21, 2016;
Accepted May 8, 2017
Caption: Figure 1. Phylogenetic tree of pCISH based on the amino acid sequence of CISH from different species. The scale bar is 0.1. The accession number of the CISH gene are as follow: Human, NM_013324, NM_145071; Rhesus monkey, XM_001097616, XM_001097824; Chimpanzee, XM_003 3 09810, XM_526202; Mouse, NM_009895; Rat, NM_031804; Cattle, NM_001046586; Chicken, NM_204626; Western clawedfrog, NM_00113689; Zebrafish, NM_001076617; Rainbow, NM_001 171852.
Caption: Figure 2. Subcellular distribution of CISH protein tagged with EGFP CISH, cytokine inducible SH2-containing protein; EGFP, enhanced green fluorescent protein.
Caption: Figure 3. The three different Taql PCR-RFLP genotypes of pCISH gene. The genotypes (AA, AG, GG) are shown at the top. M: DNA molecular marker DL2000. PCR-RFLP, polymerase chain reaction-restriction fragmentlength polymorphism; pC/SH, porcine cytokineinducible SH2-containing protein.
Caption: Figure 4. The three different Haelll PCR-RFLP genotypes of pCISH gene. The genotypes (AA, AC, CC) are shown at the top. M: DNA molecular marker DL2000. PCR-RFLP, polymerase chain reaction-restriction fragment length polymorphism; pCISH, porcine cytokine inducible SH2-containing protein.
Table 1. Primer pairs used in the analysis of the pCISH gene Primers Primer sequence (5'-3') Tm ([degrees]C) C-gap AGCCACGACCTTCTCCTAC 55 GCAGGCAGTCCACATCAG C-12 GACATGGTCCTGTGCGTTC 60 TCCAGCACCTTAGGCTCACT C-23 GAGGAGGCAGCAGAGGAGT 60 CACCAGACGGTTGATGACG C-33 AGAAGCAGTGCCCGAAGC 57 AGCAACAGCCCTCAAACC CV-1 AAGAGGAGGCGATATGTAGT 54 TAGGGAGGTGCTAATGCT CV-2 TCCCTAGCCTTCTTCAACTC 57 CATCAAGCTCAACCCGTC CV-3 GACGGGTTGAGCTTGATG 55 GCTGATGGCTGGTCTACTG CV-4 CTGCTAAGGTGGGTATGTC 54 GGTGGGCTGATTAAATGG CV-5 TCCCTCCAATAGATATTCTTACC 55 ATAGTGCTGCACAGGCTGAC CV-41 GCAGGCTTGTTAGGGTGA 57 GGGCTTAGAGGCTGGAAT CE-1 AGGGACCTTGTCCTTTGCTG 57 CCTGGCTCCTGCTGACTCTT CE-2 TCTGCGTTCAGGGACCTTGTCCTTTGCT 57 TCGGTCGACTCAGAGCTGGAAAGGGTAC CE-3 TCTGGATCCATGGTCCTCTGCGTTCAGG 57 TCGGTCGACTCAGAGCTGGAAAGGGTAC Primers Size (bp) Binding region C-gap 503 Exon 2 Exon 3 C-12 3,476 Exon 1 Exon 2 C-23 806 Exon2 Exon 3 C-33 1,004 Exon3 Exon 3 CV-1 1,555 -2,089 to -2,070 -5,552 to -535 CV-2 669 -540 to -521 Intron 1 CV-3 1,551 Intron 1 Intron1 CV-4 1,209 Intron1 Intron 1 CV-5 1,376 Intron1 Exon 3 CV-41 192 Intron 1 Intron 1 CE-1 900 CE-2 768 CE-3 783 pCISH, porcine cytokine inducible SH2-containing protein. Table 2. The association between SNP locus and diarrhea index Number (%) of piglets with different genotype Locus Breed AA(AA) GG(CC) AG(AC) SNP A-1 575G Min pig 44(0.17) 69(0.26) 162(0.61) SNP A-1 575G Landrace 83(0.42) 18(0.09) 98(0.49) SNPA2497C Min pig 126(0.39) 135(0.41) 66(0.20) SNPA2497C Landrace 73(0.36) 110(0.54) 20(0.10) Locus Chis-Seq value p value SNP A-1 575G 0.15 0.70 SNP A-1 575G 4.33 0.04 SNPA2497C 5.25 0.02 SNPA2497C 4.82 0.03 SNP, single nucleotide polymorphism.
|Printer friendly Cite/link Email Feedback|
|Author:||Niu, Buyue; Guo, Dongchun; Liu, Zhiran; Han, Xiaofei; Wang, Xibiao|
|Publication:||Asian - Australasian Journal of Animal Sciences|
|Date:||Dec 1, 2017|
|Previous Article:||Potential influence of [kappa]-casein and [beta]-lactoglobulin genes in genetic association studies of milk quality traits.|
|Next Article:||Post-prandial decrease in plasma growth hormone levels is not related to the increase in plasma insulin levels in goats.|