Characterisation and In Silico Analysis of Interleukin-4 cDNA of Nilgai (Boselaphus tragocamelus) and Indian Buffalo (Bubalus bubalis)
Interleukin-4 (IL-4) produced from Th2 cells modulates both innate and adaptive immune responses. It is a common belief that wild animals possess better immunity against diseases than domestic and laboratory animals; however, the immune system of wild animals is not fully explored yet. Therefore, a comparative study was designed to explore the wildlife immunity through characterisation of IL-4 cDNA of nilgai, a wild ruminant, and Indian buffalo, a domestic ruminant. Total RNA was extracted from peripheral blood mononuclear cells of nilgai and Indian buffalo and reverse transcribed into cDNA. Respective cDNA was further cloned and sequenced. Sequences were analysed in silico and compared with their homologues available at GenBank. The deduced 135 amino acid protein of nilgai IL-4 is 95.6% similar to that of Indian buffalo. N-linked glycosylation sequence, leader sequence, Cysteine residues in the signal peptide region, and 3′ UTR of IL-4 were found to be conserved across species. Six nonsynonymous nucleotide substitutions were found in Indian buffalo compared to nilgai amino acid sequence. Tertiary structure of this protein in both species was modeled, and it was found that this protein falls under 4-helical cytokines superfamily and short chain cytokine family. Phylogenetic analysis revealed a single cluster of ruminants including both nilgai and Indian buffalo that was placed distinct from other nonruminant mammals.
The discoveries of Interleukin-1 (IL-1) and IL-2 led to a better understanding of the effects of ILs, and till now more than 40 cytokines are discovered with specific functions . Interleukin-4 (IL-4) is one of the extensively studied cytokines which induces specific functions in wide range of immune cells defining its pleotropic character . IL-4 was identified originally as a B cell growth factor-1 in mice  and was subsequently shown to modulate other cellular interactions of immune response . It is the primary cytokine which promotes the development of Th2 effector cells and antagonises the activity of interferon gamma (IFN-) induced development of Th1 cells [5, 6]. Upon activation by IL-4, Th2 cells subsequently produce additional IL-4. These cytokines act synergistically with IL-5 to either activate IgE producing B cells or induce isotype switching and enhance IgE mediated responses in allergy and asthma [7–9]. IL-21 that was discovered recently is homologous to IL-4 in its ability to modulate both innate and adaptive immune responses . Wide diversity of IL-4 activity reported to date suggests that it is a key regulator in humoral and adaptive immunity.
The gene encoding IL-4 is found in chromosomes 11, 5, and 7 in mouse , human , and cattle , respectively. In mouse and human, the gene comprises 4 exons spanning 6 kb and 10 kb, respectively [14, 15]. This cytokine was initially cloned and characterised in mouse and human [16–18]. Further exploration was carried out by characterising it in domestic animals like dog , cat , camel , horse , pig , and so forth. In addition, IL-4 of some ruminants like cow , African buffalo , sheep , and goat  was also previously cloned and identified. In addition, IL-4 has been identified and reported in chimpanzee  and bottle-nosed dolphin . Wild animals are presumed to possess stronger immune system as compared to their domestic counterparts. Due to difference in habitat/environment, the immune function of wild animals could be different from that of laboratory bred/domestic animals . A comparison of sequence encoding IL-4 among various wild and domestic species could explain the difference, if any, in structure and function with respect to this cytokine. Indian buffaloes are the centre of dairy industry and are fast replacing indigenous cattle in contribution towards total milk production of India. In India, nilgai is sympatric with one domestic ruminant, that is, Indian buffalo. The present study reports the characterisation of IL-4 of nilgai (Boselephus tragocamelus) and Indian buffalo (Bubalus bubalis) as a model for the comparison of wild versus domestic ruminants of Bovidae family and their phylogenetic lineage.
2. Materials and Methods
2.1. Sample Collection and RNA Isolation
Total RNA was isolated from peripheral blood mononuclear cells (PBMs). Blood was obtained aseptically by jugular puncture from nilgai maintained in semicaptivity at Deer Park, Indian Veterinary Research Institute (IVRI), Izatnagar, and Indian buffalo from slaughter house, Bareilly. PBM cells were extracted using Histopaque 1077 (Sigma, USA) density gradient centrifugation following a method previously described  and stimulated with Concanavalin A (Con A) at the concentration of 10 μg/mL for 20 h at 37°C in a humidified incubator with 5% CO2. Total RNA of both the samples was isolated using Trizol LS reagent (Life Technologies, New York, NY) following the manufacturer’s instructions.
2.2. cDNA Synthesis and Amplification
Two respective first strands of cDNA were synthesized at 37°C from two RNA samples by using oligo dT primers (Promega, Madison, WI). Nilgai and Indian buffalo IL-4 genes were amplified from their respective cDNA using specific oligonucleotide primers (Forward 5′-TAATGGGTCTCACCTACCAG-3′ and Reverse 5′-TTCAGCTTCAACACTTGGAG-3′) designed based on the sequence of cattle (Accession NM_173921.2). The oligonucleotide primers were designed using OLIGO 4.0 software (USA). The IL-4 specific cDNAs were amplified using sequence specific primers (50 pmol/μL) 1.0 μL each; Template cDNA 1.0 μL; dNTPs (10 mM) 1.0 μL; 10X Taq polymerase buffer 5 μL; 25 mM MgCl2 3 μL; Taq DNA polymerase (MBI Fermentas, 5 U/μL) 1.0 μL; and nuclease free water making final reaction mixture of volume 50 μL. PCR amplification program followed was: 95°C for 5 min, 35 repeated cycles of 1 min denaturation at 94°C, 1 min annealing at 60°C and 1 min extension at 72°C, and one cycle of final extension at 72°C for 10 min. The PCR amplified product was analysed on 1% agarose gel containing ethidium bromide along with DNA molecular weight marker.
2.3. cDNA Cloning and Sequencing
The amplified products were purified from the agarose gel using Gel extraction Kit (Qiagen, Germany). Nilgai IL-4 PCR product was cloned into pTZ57R/T vector (MBI Fermentas, MD) and buffalo amplified product using pGEMT-Easy (Promega, Madison, USA) vector following the manufacturers’ protocol and further screened by blue white screening. The recombinant plasmids were characterized by restriction enzymes NotI, PstI, NcoI, and EcoRI (MBI Fermentas, MD) and by PCR using gene specific primers predicted from cattle IL-4 sequence.
2.4. Sequencing and Analysis
The characterized plasmids were sequenced using T7 and SP6 universal primer using ABI PRISM 377 Version 3.0 DNA sequencer (Applied Biosystem, Foster city, CA). The nucleotide sequences of both insert IL-4 were first BLAST analyzed (http://www.ncbi.nlm.nih.gov/) and further submitted to GenBank. Multiple sequence alignment was carried out with IL-4 gene sequences of nilgai and Indian buffalo with its homologues from other species like cattle (Bos taurus) (GenBank Accession no. NM_173921), African buffalo (Syncerus caffer) (EU000421), goat (Capra hircus) (U34273), sheep (Ovis aries) (M96845), pig (Sus scrofa) (JF906512), camel (Camelus dromedarius) (HM051106), red deer (Cervus elaphus) (L07081), giraffe (Giraffa camelopardalis) (EU000423), bison (Bison bonasus) (EU000422), llama (Lama glama) (AB107648), dog (Canis lupus familiaris) (NM_001003159), cat (Felis catus) (NM_001043339), and bottle-nosed dolphin (Tursiops truncatus) (AB020732).
Amino acid sequences were predicted using DNA Star software (Lasergene). Nucleotide and deduced amino acid sequence were aligned to predict phylograms using Mega 5.1 software . Nilgai and Indian buffalo IL-4 protein structure was predicted using PHYRE2 software (Protein Homology/analog Y Recognition Engine; http://www.sbg.bio.ic.ac.uk/phyre2). The N-glycosylation sites were predicted using HIV sequence database (http://www.hiv.lanl.gov). Leader peptide cleavage site was predicted using SignalP 4.1 server (http://www.cbs.dtu.dk) .
The concentration of RNA was measured using UV spectrometer, and the purity and integrity were checked by analyzing the ratio of optical density (OD) at 260 and 280 nm. The ratios of OD260/OD280 in total RNA from nilgai and Indian buffalo were found to be 1.83 and 1.85, respectively. Amplification of cDNA through PCR was confirmed through agarose gel electrophoresis which gave a product size 417 bp in both the cases. Purified PCR product of respective species was cloned, and the recombinant plasmid was characterized by restriction analysis (Figure 1, Lanes 2–5) and sequencing. Recombinant plasmid was linearised by digesting with BamH1 (Figure 1, Lane 1); insert was released from vector using NotI enzyme (Figure 1, Lane 2). The presence of insert was confirmed by digesting recombinant plasmid with PstI which yielded a single band of around 200 bp (Figure 1, Lane 3). Digestion with NcoI yielded product of approximately 70 bp (Figure 1, Lane 4), and EcoRI yielded two fragments around 134 and 270 bp (Figure 1, Lane 5). Since the results of characterization of buffalo and nilgai plasmids in agarose gel electrophoresis were identical, the results of nilgai IL-4 are shown in Figure 1.
The sequences encoding nilgai and Indian buffalo IL-4 cDNA were assigned GenBank accession numbers AY939910 and AY293620, respectively. The full length cDNA of nilgai and Indian buffalo IL-4 contained open reading frame (ORF) of size 408 bp each encoding a protein of 135 amino acids. The molecular weight and isoelectric point of protein were predicted to be 15.039 and 8.785 kDa, respectively, in case of nilgai; however, in case of Indian buffalo these values were found to be 15.158 and 8.940 kDa, respectively.
A comparison of deduced amino acid sequences is provided in Figure 2. Multiple sequence alignment revealed conserved leader sequence of 24 amino acids and a potential N-linked glycosylation site, that is, Asn-Thr-Thr at positions 64–66. Cys (C) residues at positions 13 and 17 in the signal peptide were found to be conserved except in dog and cat for the 17th position. One significant difference is the presence of Gly at position 123 in case of nilgai IL-4 protein, whereas the Arg 123 residue was conserved across other species.
Comparative nonsynonymous nucleotide substitutions leading to change in amino acid at different positions of various species as compared to nilgai are given in Table 1.
Phylogenetic analysis based on the nucleotide and predicted amino acid (Figure 3) sequences of IL-4 revealed that the ruminants formed a single cluster indicating their recent divergence from other mammalian species. The phylogenetic tree showed that nilgai IL-4 was more related to cattle, Indian buffalo, African buffalo, and bison than to other species sequence included in the comparison.
From the deduced amino acid sequences, similar tertiary structure was predicted for both Indian buffalo and nilgai IL-4 proteins (Figure 4).
IL-4 is one of the key cytokines in Th2 mediated immune responses, which has been shown to regulate the responses of many other cytokines like IL-1, interferon-gamma, and tumor necrosis factor-alpha. Several reports are there regarding human, murine, and domestic animal IL-4 with scanty reports on wild ruminants. In the present study, nilgai and Indian buffalo IL-4 cDNA was sequenced, and amino acid sequences were predicted for the precursor of the protein in both these species that were further compared with other sequences available in the database. Upon alignment of nilgai IL-4 ORF region sequences with its homologues revealed maximum similarity with cattle (98.3%) followed by Indian buffalo (97.5%), sheep (95.6%), goat (93.9%), camel (85.3%), dolphin (84.3%), pig (81.6%), cat (73.9%), and dog (72.2%). The deduced amino acid sequences showed that nilgai IL-4 is highly similar to cattle with 97.8% match followed by 95.6% with Indian buffalo, 90.4% with goat, and 91.9% with sheep.
One striking difference is the presence of Gly at position 123 in case of nilgai IL-4 protein instead of conserved Arg 123 residue across other species. At this position, the arginine which is a polar amino acid changed to a nonpolar neutral amino acid glycine. This could lead to minor changes in the folding of protein as evident from the predicted tertiary structure (Figure 4) and cause observed difference in isoelectric point that may influence the activity of this cytokine in both the species.
Nilgai and Indian buffalo amino acid sequence revealed a leader sequence of 24 amino acids, and their mature peptides are predicted to be of molecular weight 12.44 kDa and 12.56 kDa, respectively. Similar finding was earlier reported in cattle .
A potential N-linked glycosylation sequence Asn-Thr-Thr (positions 64, 65, and 66) is found to be conserved in all species. Similar findings were also reported in cattle  and in human IL-4 . In addition to earlier reported glycosylation site additional N-glycosylation sites were identified upon in silico analysis (Table 1). But whether glycosylations occur in all these sites is not yet to be established.
Sequence analysis also revealed that the Cys (C) residues at positions 13 and 17 in the signal peptide were found conserved in all the species except in dog and cat for 17th position. Similar result was reported on comparison of ovine and bovine IL-4 sequence . N-linked glycosylation sites and Cys residues were found to be located in the same position in all species. This observation suggests that this region is highly conserved in all the species and may play an important role in determining tertiary structure and functional integrity of the cytokine.
It was observed that 3′ UTR of IL-4 gene contains A+T rich stretches which include both tandem repeats of TAAT or ATTTA and also the polyadenylation signal sequence. Similar observations were earlier reported in cattle and human [18, 24].
Findings of phylogram that ruminants form a cluster, and nilgai IL-4 is evolutionarily closer to buffalo and cattle than other mammals studied, were also corroborated previously on different cytokines of nilgai, that is, in IL-2  and IL-18 .
A comparative analysis on nonsynonymous nucleotide substitutions leading to change in amino acid at different positions of various species as compared to nilgai is given in Table 1. In spite of six variations in predicted amino acid sequences, the tertiary structure predicted for both Indian buffalo and nilgai IL-4 proteins was nearly the same (Figure 4). It is evident that both these modeled proteins fall under 4-helical cytokines superfamily and short-chain cytokine family. Since 81% of the amino acid sequences submitted have been modeled with 100% confidence by the single highest scoring template, few alterations in amino acid between two species did not result in change in the predicted structure of the protein.
This comparison of nilgai and Indian buffalo IL-4 precursors will be useful to correlate the molecular aspect of immunity in wild and domestic ruminants.
Conflict of Interests
The authors declare that they have no conflict of interests.
The authors are thankful to the Central Zoo Authority, New Delhi, for providing financial help and to the Director, Joint Director (Research), and Joint Director (Academic), Indian Veterinary Research Institute, Izatnagar, for providing necessary facilities to carry out the work.
G. Shubinsky and M. Schlesinger, “Kinetics of the pleiotropic effect of interleukin 4 on the surface properties of human B-Lymphoma cells,” Leukemia and Lymphoma, vol. 15, no. 3-4, pp. 333–340, 1994.View at: Google Scholar
M. Howard, J. Farrar, M. Hilfiker et al., “Identification of a T cell-derived B cell growth factor distinct from interleukin 2,” Journal of Experimental Medicine, vol. 155, no. 3, pp. 914–923, 1982.View at: Google Scholar
W. E. Paul, “Interleukin 4: signalling mechanisms and control of T cell differentiation,” CIBA Foundation Symposia, no. 204, pp. 208–219, 1997.View at: Google Scholar
P. D'Eustachio, M. Brown, C. Watson, and W. E. Paul, “The IL-4 gene maps to chromosome 11, near the gene encoding IL-3,” Journal of Immunology, vol. 141, no. 9, pp. 3067–3071, 1988.View at: Google Scholar
M. M. Le Beau, R. S. Lemons, R. Espinosa III, R. A. Larson, N. Arai, and J. D. Rowley, “Interleukin-4 and interleukin-5 map to human chromosome 5 in a region encoding growth factors and receptors and are deleted in myeloid leukemias with a del(5q),” Blood, vol. 73, no. 3, pp. 647–650, 1989.View at: Google Scholar
J. Buitkamp, F. W. Schwaiger, S. Solinas-Toldo, R. Fries, and J. T. Epplen, “The bovine interleukin-4 gene: genomic organization, localization, and evolution,” Mammalian Genome, vol. 6, no. 5, pp. 350–356, 1995.View at: Google Scholar
T. Otsuka, D. Villaret, T. Yokota et al., “Structural analysis of the mouse chromosomal gene encoding interleukin 4 which expresses B cell, T cell and mast cell stimulating activities,” Nucleic Acids Research, vol. 15, no. 1, pp. 333–344, 1987.View at: Google Scholar
N. Arai, D. Nomura, D. Villaret et al., “Complete nucleotide sequence of the chromosomal gene for human IL-4 and its expression,” Journal of Immunology, vol. 142, no. 1, pp. 274–282, 1989.View at: Google Scholar
F. Lee, T. Yokota, T. Otsuka et al., “Isolation and characterization of a mouse interleukin cDNA clone that expresses B-cell stimulatory factor 1 activities and T-cell- and mast-cell-stimulating activities,” Proceedings of the National Academy of Sciences of the United States of America, vol. 83, no. 7, pp. 2061–2065, 1986.View at: Google Scholar
Y. Noma, P. Sideras, T. Naito et al., “Cloning of cDNA encoding the murine IgG1 induction factor by a novel strategy using SP6 promoter,” Nature, vol. 319, no. 6055, pp. 640–646, 1986.View at: Google Scholar
T. Yokota, T. Otsuka, T. Mosmann et al., “Isolation and characterization of a human interleukin cDNA clone, homologous to mouse B-cell stimulatory factor 1, that expresses B-cell- and T-cell-stimulating activities,” Proceedings of the National Academy of Sciences of the United States of America, vol. 83, no. 16, pp. 5894–5898, 1986.View at: Google Scholar
R. Harley, C. R. Helps, D. A. Harbour, T. J. Gruffydd-Jones, and M. J. Day, “Cytokine mRNA expression in lesions in cats with chronic gingivostomatitis,” Clinical and Diagnostic Laboratory Immunology, vol. 6, no. 4, pp. 471–478, 1999.View at: Google Scholar
G. Nagarajan, S. K. Swami, S. K. Ghorui, K. M. L. Pathak, R. K. Singh, and N. V. Patil, “Cloning and sequence analysis of IL-2, IL-4 and IFN-γ from Indian Dromedary camels (Camelus dromedarius),” Research in Veterinary Science, vol. 92, no. 3, pp. 420–426, 2012.View at: Publisher Site | Google Scholar
C. R. Engwerda and R. M. Sandeman, “The isolation and sequence of sheep interleukin 4,” DNA Sequence, vol. 3, no. 2, pp. 111–113, 1992.View at: Google Scholar
I. Gautherot, N. Burdin, D. Seguin, L. Aujame, and R. Sodoyer, “Cloning of interleukin-4 delta2 splice variant (IL-4δ2) in chimpanzee and cynomolgus macaque: phylogenetic analysis of δ2 splice variant appearance, and implications for the study of IL-4-driven immune processes,” Immunogenetics, vol. 54, no. 9, pp. 635–644, 2002.View at: Publisher Site | Google Scholar
Y. Inoue, T. Itou, T. Sakai, and T. Oike, “Cloning and sequencing of a bottle-nosed dolphin (Tursiops truncatus) interleukin-4-encoding cDNA,” Journal of Veterinary Medical Science, vol. 61, no. 6, pp. 693–696, 1999.View at: Google Scholar
P. E. Baker and K. F. Knoblock, “Bovine costimulator. I. Production kinetics, partial purification, and quantification in serum-free Iscove's medium,” Veterinary Immunology and Immunopathology, vol. 3, no. 4, pp. 365–379, 1982.View at: Google Scholar
K. Tamura, D. Peterson, N. Peterson, G. Stecher, M. Nei, and S. Kumar, “MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods,” Molecular Biology and Evolution, vol. 28, no. 10, pp. 2731–2739, 2011.View at: Publisher Site | Google Scholar
D. K. Das, M. Saini, D. Swarup, and P. K. Gupta, “Comparison of nucleotide and amino acids sequence of Nilgai (Boselaphus tragocamelus) interleukin-18 (IL-18) with other ruminants,” Indian Journal of Biotechnology, vol. 7, no. 2, pp. 195–199, 2008.View at: Google Scholar