IF: 0.957
Cite Score:

Gene Profile Expression Related to Type I Interferons in HT-29 Cells Exposed to Cryptosporidium parvum


Seyede Manizhe Heidarnejadi 1 , 2 , 3 , Abdollah Rafiei 1 , 2 , * , Manoochehr Makvandi ORCID 1 , 4 , Majid Pirestani 5 , Jasem Saki ORCID 6 , Ataallah Ghadiri 7


1 Infectious and Tropical Diseases Research Center, Health Research Institute, Ahvaz Jundishapur University of Medical Sciences, Ahvaz, IR Iran

2 Department of Medical Parasitology, School of Medicine, Ahvaz Jundishapur University of Medical Sciences, Ahvaz, IR Iran

3 Shoshtar Faculty of Medical Sciences, Shoshtar, IR Iran

4 Virology Department, Infectious and Tropical Disease Research Center, Health Research Institute, Ahvaz Jundishapur University of Medical Sciences, Ahvaz, IR Iran

5 Department of Parasitology and Entomology, Faculty of Medical Sciences, Tarbiat Modares University, Tehran, IR Iran

6 Department of Medical Parasitology, School of Medicine, Ahvaz Jundishapur University of Medical Sciences, Ahvaz, IR Iran

7 Department of Immunology, School of Medicine, Ahvaz Jundishapur University of Medical Sciences, Ahvaz, IR Iran


Jundishapur Journal of Microbiology: 11 (7); e63071
Published Online: July 2, 2018
Article Type: Research Article
Received: October 27, 2017
Revised: June 13, 2018
Accepted: June 17, 2018




Background: Cryptosporidium parvum may contribute to upregulation or downregulation of host cellular genes among which type I (α, β) and type II (γ) interferons play key roles to eliminate infectious agents.

Objectives: The current study aimed at evaluating the expressed genes related to human type I interferon response in HT-29 cell line after exposure to C. parvum for six and 24 hours.

Methods: Subsequently, the overexpression and under expression of 84 human genes related to type I interferons were investigated using RT2Profiler™ PCR (polymerase chain reaction) Array.

Results: Four top overexpressed genes including IL10, SH2D1A, MX1, and HLA-A after six hours exposure, and 10 top overexpressed genes as IL15, DDX58, CXCL10, NMI, MYD88, STAT3, IFNAR2, IFIH1, CASP1, and TLR3 were observed after 24 hours exposure. Five underexpressed genes such as TIMP1, TYK2, IRF2, PML and IRF5 were monitored for six hours.

Conclusions: The current study findings revealed that the overexpressed genes IL15, TIMP1, and SH2D1A may have an important role to inhibit the invasion of C. parvum. Also, the overexpressed genes, namely SH2D1A, MX1, and NMI, may have antiviral properties while TIMP1 may have anticancer properties. Further, the pertinent results demonstrated that the type I interferons and the relevant genes had significant effects on stimulating innate immune system against C. parvum.


HT29 Cells Type I Interferon Response Gene Profile Expression Cryptosporidium parvum Iowa Strain

Copyright © 2018, Author(s). This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial 4.0 International License (http://creativecommons.org/licenses/by-nc/4.0/) which permits copy and redistribute the material just in noncommercial usages, provided the original work is properly cited
1. Background

Cryptosporidium is a coccidian protozoan parasite, which belongs to the phylum Apicomplexa and causes cryptosporidiosis. The spread of infection hinges on several factors such as environmental situation, host growing old, immune system status, geographic locality, and contact with infected humans and/or animals (1). Epithelial cells contribute primary and quick defense in opposition to cryptosporidiosis, and also activate immune effector cells to the infectivity position. The attachment of Cryptosporidium parvum to apical cell surface can stimulate host cell signal pathways and hence modify cell function (2). The development of C. parvum infection in intestinal epithelial cells is effectively controlled by both type I and II interferons (IFNs) (3). Type I IFNs bring cell-inherent antimicrobial states in infected and adjacent cells that confine the increase of infectious agents, mainly viral pathogens.

The investigation on type I IFNs in immunity to parasitic protozoa is limited; however, the role of type I IFNs showed that they reduce the severity of INF of Leishmania major, and that they can act in a straight line on host cells to induce parasite inactivation. Some studies showed that the development of C. parvum infection in intestinal epithelial cells is effectively controlled both by type I and type II IFNs. The enterocytes can produce type I IFNs in rejoinder to C. parvum infection and lead to the inhibition of parasite development (3). The cultivation of C. parvum in tissue culture provided a new approach to recognize different gene expressions and immunological responses to different parasite proteins. Initially, primary epithelial cell culture (PECs) can stimulate different cellular genes after being exposed to C. parvum (4). Additionally, the continuous cell lines HCT8 and HT-29 showed good support for C. parvum (5). The implementation of advanced molecular tool can assess the gene expressions and their relative functions. The application of RT2Profiler™ PCR Array is a golden standard, which can measure low levels of gene expression (6).

2. Objectives

Since very few publications with focus on gene expression profile in the cells exposed to C. parvum oocysts can be found, the current study aimed at evaluating overexpressed and underexpressed genes related to type I INFs in HT29 cells infected by C. parvum using RT2 Profiler™ PCR Array.

3. Methods
3.1. Ethics Statement

The current study was granted by research deputy (Grant no: 94147) and ethically approved by the Ethics Committee of Ahvaz Jundishapur University of Medical Sciences.

3.2. Parasite Provision and Infection of Epithelial Cells

Cryptosporidium parvum oocysts (Iowa strain, BTF Company from Sydney –Australia) were kept in a solution of 0.01% Tween 20, 100 U/mL penicillin, and 100 mg/mL gentamicin, and stored at 4°C. To determine the infectivity of C. parvum, two-day-old neonatal rats (Razi Vaccine and Serum Research Institute, Karaj, Iran) were inoculated orally with 100,000 - 120,000 C. parvum oocysts (7). The shedding of the oocysts in the feces of neonatal rats were observed after eight days following the infection and confirmed with Ziehl–Neelsen staining.

3.3. Treatment of Oocysts

Oocysts were washed with 0.15 M/L phosphate-buffered saline (PBS; pH 7.2), centrifuged (3000 g for 15 minutes at 4°C) and treated with acidic H2O (pH 2.5; 20 minutes, vortexed at five minutes interval at 37°C). The oocysts were centrifuged (3000 g for four minutes at 20°C), re-suspended in parasite-maintenance Dulbecco's modified eagle's medium (DMEM; Gibco, Auckland, New Zealand) containing 4.5 g/L d-glucose, 0.58 g/L l-glutamine, 3.7 g/L sodium bicarbonate, 0.20 g/L bovine bile, 0.004 g/L folic acid, 0.001 g/L 4-aminobenzoic acid, 0.004 g/L d-calcium pantothenate, 0.88 g/L ascorbic acid, 1% heat-inactivated FBS, 2.4 g/L HEPES (4-(2-hydroxyethyl)-1-piperazineethanesulfonic acid), 100,000 U/L penicillin, 100,000 mg/L streptomycin, and 250 mg/L amphotericin B, adjusted to pH 7.4), and incubated for three hours at 37°C. The supernatant was removed, fresh DMEM was added, and oocysts were kept at 37°C for 90 minutes. The final dose of oocysts was adjusted 5 000/well prior to inoculation and cell culture (8).

3.4. In Vitro Cultivation
3.4.1. Preparation and Infection of Host Cells

The cell line HT-29 was purchased from Iranian Biological Research Centre, Tehran. The HT-29 cells were harvested in 25-cm2 flasks, followed by subculture in 24-well plates, supplemented with DMEM containing 10% fatal bovine serum (FBS; Gibco, Auckland, New Zealand), 100 U/mL penicillin, 100µg/mL streptomycin, and 25 µg/mL amphotericin B incubated with 5% CO2 at 37°C for 48 hours, and confluent monolayer cells were developed. For 24-well plates, the supernatant DMEM of the wells were removed and washed once with PBS (pH = 7.2). Next, 100 µL DMEM containing 5000 oocysts were inoculated onto each well, and the plates were incubated with 5% CO2 at 37°C for six hours. The 24-well plates were also inoculated with 100 µL DMEM containing 5000 oocysts/cm2 and incubated with 5% CO2 at 37°C for 24 hours. A 24-well plate was inoculated with PBS as control (mock). The HT-29 cell lines were exposed to C. parvum for 6 and 24 hours, and then, the cells were lysed with buffer RLT (RNeasy mini kit, Qiagen).

3.5. RNA Isolation and RT2 Profiler™ PCR Array

The total RNA was extracted for 6 and 24 hours separately using RNeasy mini kit (Qiagen), as per the manufacturer’s instructions. The cDNA was prepared using RT2 first standard kit (Qiagen). The qPCR Array was performed using the RT2 Profiler™ PCR Array of SABiosciences (http://www.sabiosciences.com/ArrayList.php) according to the manufacturer’s instruction. The triplicate test was conducted on each sample. Briefly, 2700 µL reaction mixture containing 1350 µL 2x RT2 SYBR Green qPCR Mastermix, 102µL cDNA, and 1248 µL RNAase-free water were prepared. The 25 μL of reaction mixture was added to each 96-well custom RT2 Profiler™ PCR Arrays and was subjected to real-time PCR (Applied Biosystems, Foster City, CA) with the following program: 95°C, 10 minutes, one cycle, 40 cycle, 15 seconds, 95°C, followed by 60°C for one minute (fluorescence data collection). The threshold cycle (Ct) for each well was calculated as described in web-based PCR Array Data Analysis Software (www.SABiosciences.com/pcrarraydataanalysis.php). A triplicate test was conducted to evaluate gene expression.

The gene expressions of the each sample were appraised on 84 genes (Table 1) related to human type I INFs (PAHS-016ZA). In addition, five housekeeping genes (ACTB, B2M, GAPDH, HPRT1, PRLPO), one genomic DNA control (GDC), three reversed transcription control (RTC), and three positive PCR control were used in the assay. The following overexpression and underexpression values were determined by fold change assessment. Fold-change (2^ [-delta delta Ct]) is the normalized gene expression (2^(-delta Ct)) in the test sample divided by the normalized gene expression (2^ [-delta Ct]) in the control sample. Fold-regulation represents fold-change results in a biologically significant way. Fold-change values > 1 indicate a positive or upregulation; the fold-regulation is equal to the fold-change. Fold-change values > 1 indicate a negative or down-regulation, and the fold-regulation is the negative inverse of the fold-change.

Table 1. The Profile of Human Genes Related to Type I Interferon by RT² Profiler™ PCR Array
Gene SymbolDescription
ADARAdenosine deaminase, RNA-specific
BAG3BCL2-associated athanogene 3
BST2Bone marrow stromal cell antigen 2
CASP1Caspase 1, apoptosis-related cysteine peptidase (interleukin 1, beta, convertase)
CAV1Caveolin 1, caveolae protein, 22 kDa
CCL2Chemokine (C-C motif) ligand 2
CCL5Chemokine (C-C motif) ligand 5
CD70CD70 molecule
CD80CD80 molecule
CD86CD86 molecule
CDKN1BCyclin-dependent kinase inhibitor 1B (p27, Kip1)
CIITAClass II, major histocompatibility complex, transactivator
CRPC-reactive protein, pentraxin-related
CXCL10Chemokine(C-X-C motif) ligand 10
DDX58DEAD(Asp-Glu-Ala-Asp) box polypeptide 58
EIF2AK2Eukaryotic translation initiation factor 2-alpha kinase 2
GBP1Guanylate binding protein 1, interferon-inducible
HLA-AMajor histocompatibility complex, class I, A
HLA-BMajor histocompatibility complex, class I, B
HLA-EMajor histocompatibility complex, class I, E
HLA-GMajor histocompatibility complex, class I, G
IFI16Interferon, gamma-inducible protein 16
IFI27Interferon, alpha-inducible protein 27
IFI30Interferon, gamma-inducible protein 30
IFI6Interferon, alpha-inducible protein 6
IFIH1Interferon induced with helicase C domain 1
IFIT1Interferon-induced protein with tetratricopeptide repeats 1
IFIT2Interferon-induced protein with tetratricopeptide repeats 2
IFIT3Interferon-induced protein with tetratricopeptide repeats 3
IFITM1Interferon induced transmembrane protein 1(9–27)
IFITM2Interferon induced transmembrane protein 2(1–8D)
IFIInterferon induced transmembrane protein 3
IFNA1Interferon, alpha 1
IFNA2Interferon, alpha 2
IFNA4Interferon, alpha 4
IFNAR1Interferon (alpha, beta and omega) receptor 1
IFNAR2Interferon (alpha, beta and omega) receptor 2
IFNB1Interferon, beta 1, fibroblast
IFNEInterferon, epsilon
IFNW1Interferon, omega 1
IL10Interleukin 10
IL15Interleukin 15
IL6Interleukin 6 (interferon, beta 2)
IRF1Interferon regulatory factor 1
IRF2Interferon regulatory factor 2
IRF3Interferon regulatory factor 3
IRF5Interferon regulatory factor 5
IRF7Interferon regulatory factor 7
IRF9Interferon regulatory factor 9
ISG15ubiquitin-like modifier
ISG20Interferon stimulated exonuclease gene 20kDa
JAK1Janus kinase 1
JAK2Janus kinase 2
MALMal, T-cell differentiation protein
METMet proto-oncogene(hepatocyte growth factor receptor)
MNDAMyeloid cell nuclear differentiation antigen
MX1Myxovirus (influenza virus) resistance 1, interferon-inducible protein p78(mouse)
MX2Myxovirus (influenza virus) resistance 2(mouse)
MYD88Myeloid differentiation primary response gene(88)
NMIN-myc (and STAT) interactor NM_ NOS2 Nitric oxide synthase 2, inducible
OAS12'-5'-oligoadenylate synthetase 1, 40/46kDa
OAS22'-5'-oligoadenylate synthetase 2, 69/71kDa
PMLPromyelocytic leukemia
PRKCZProtein kinase C, zeta
PSME2Proteasome(prosome, macropain) activator subunit 2(PA28 beta)
SH2D1ASH2 domain containing 1A
SHBSrc homology 2 domain containing adaptor protein B
SOCS1Suppressor of cytokine signaling 1
STAT1Signal transducer and activator of transcription 1, 91kDa
STAT2Signal transducer and activator of transcription 2, 113kDa
STAT3Signal transducer and activator of transcription 3 (acute-phase response factor)
TAP1Transporter 1, ATP-binding cassette, sub-family B (MDR/TAP)
TICAM1Toll-like receptor adaptor molecule 1
TIMP1TIMP metallopeptidase inhibitor 1
TLR3Toll-like receptor 3
TLR7Toll-like receptor 7
TLR8Toll-like receptor 8 Toll-like receptor
TMEM173Transmembrane protein 173
TNFSF10Tumor necrosis factor (ligand) superfamily, member 10
TRAF3TNF receptor-associated factor 3
TYK2Tyrosine kinase 2
GAPDHGlyceraldehyde-3-phosphate dehydrogenase
HPRT1Hypoxanthine phosphoribosyltransferase 1
RPLP0Ribosomal protein, large, P0
HGDCHuman genomic DNA contamination
RTCReverse transcription control
RTCReverse transcription control
RTCReverse transcription control
PPCPositive PCR control
PPCPositive PCR control
PPCPositive PCR control
3.6. Statistical Analysis

The results of overexpression and underexpression of each gene related to type I INFs were compared with mock cells using online Software (www.SABiosciences.com/pcrarraydataanalysis.php). In vitro, cellular assays were done in triplicate. Data were analyzed using student t test. Based on fold change, the overexpression of genes and underexpression of genes were determined for six and 24 hours after infection. P value < 0.05 was considered significant. The P values are calculated based on a Student t test of the replicate 2^(-delta Ct) values for each gene in the control and treatment groups. The automatic selection from housekeeping genes with genes (HKG) panel with ACTB and GAPDH was conducted to evaluate the normalization of the test samples for six and 24 hours. This method automatically selects an optimal set of internal control/housekeeping/normalization genes for the analysis from the available HKG panel on the PCR array.

4. Results
4.1. Inoculate Cryptosporidium parvum Iowa Strain

The post-infected neonatal rat with C. parvum showed symptoms such as weakness, severe diarrhea, dehydration, and volume depletion three to four days after oral inoculation. The microscopy presence of Cryptosporidium oocysts in feces of the infected neonatal rats was stained by Ziehl-Neelsen. The amount of 522 ng/µL of the extracted RNA from exposed cells to C. parvum was used for RT² PCR™ Array.

4.2. Gene Expression Pattern of Human Type I INF Response at RT² Profiler™ PCR Array

Based on technical triplicate data, fold change for the significant upregulation of the expressed gene was defined (fold-change > 1.0: P < 0.05) for six and 24 hours and fold change for significant downregulated gene was described (fold-change < 1.0; P < 0.05) following the infection. It should be noted that the GSEA analysis confirmed all the overexpressed genes related to type I INFs pathway.

4.3. Overexpressed, Underexpressed, and Unchanged Genes Six Hours Post-Infection

Four top overexpressed genes related to type I INFs namely HLA-A, IL-10, MX1, and SH2D1A were selected (Table 2 and Figure 1) and out of 60 downexpressed genes, five genes including IRF2, IRF5, PML, TIMP1, and TYK2 were chosen for further investigation (Table 3 and Figure 1). The remaining 55 underexpressed genes were as follows: CAV1, TAP1, TLR9, CD80, JAK1, IFI30, IFI27, IFITM2, SHB, OAS2, IFI6, IFNAR2, BST2, ADAR, STAT2, MET, PRKCZ, IRF9, IFNAR1, MX2, TICAM1, SOCS1, JAK2, IFI16, TRAF3, STAT1, BAG3, TLR3, ISG15, TNFSF10, IRF3, IFNA2, HLA-G, CXCL10, ISG20, STAT3, CASP1, IRF7, MYD88, NMI, IL15, IFITM3, IFNE, IFNA1, IRF1, IFITM1, IFNW1, IFIT1, IL6, IFIT3, IFIH1, NOS2, IFIT2, CIITA, and CCL5. Twenty genes residual unchanged were as follows: CCL2, CD70, CD86, CDKN1B, CRP, DDX58, EIF2AK2, GBP1, HLA-B, HLA-E, IFNA4, IFNB1, MAL, MNDA, OAS1, PSME2, TLR7, TLR8, TMEM173, and VEGFA.

Table 2. Overexpressed Genes for 6 h (Fold-change > 1; P < 0.05)
Overexpressed GenesFold-change 6 hP-value h
Table 3. Underexpressed Genes for 6 h (Fold-change < 1; P < 0.05)
Under expressed GenesFold change 6 hoursP value hours
The scatter pattern shows overexpressed and downexpressed genes in infected HT-29 cells with Cryptosporidium parvum oocysts for 6 hours.
Figure 1. The scatter pattern shows overexpressed and downexpressed genes in infected HT-29 cells with Cryptosporidium parvum oocysts for 6 hours.
4.4. Overexpressed, Underexpressed, and Unchanged Genes 24 Hours Post-Infection

Sixty-seven genes were overexpressed among which 10 top genes such as CASP1, CXCL10, DDX58, IFIH1, IFNAR2, IL15, MYD88, NMI, STAT3, and TLR3 were selected for further studies (Table 4, Figure 2). The following 57 overexpressed genes were as follows: IFNE, IFNB1, CAV1, EIF2AK2, CDKN1B, ISG15, CCL5, ISG20, BAG3, PRKCZ, GBP1, IFNAR1, IFITM3, TIMP1, STAT1, HLA-A, VEGFA, JAK1, IFIT2, IFIT3, ADAR, TRAF3, IFNA1, OAS2, HLA-E, TMEM173, IRF2, IFITM2, IFI16, IFI27, SOCS1, MX2, PML, TYK2, SHB, OAS1, MX1, IFI30, CIITA, BST2, NOS2, IRF7, IRF1, CD70, TNFSF10, STAT2, TICAM1, IRF9, IRF5, IFNW1, IFI6, IFNA4, TAP1, CCL2, MNDA, IFNA2, and IL6. Three underexpressed genes IFITM1, PSME2, and IRF3 were observed in post-infected cells. The remaining 14 unchanged genes related to type I INFs were diagnosed in the infected cells as follows: TLR7, TLR8, TLR9, SH2D1A, JAK2, MAL, MET, IL10, IFIT1, HLA-B, HLA-G, CRP, CD80, and CD86. The cluster gram pattern showed that the overexpressed and underexpressed genes were related to human type I INF after six and 24 hours exposure (Figure 3 A & B).

Table 4. Overexpressed Genes for 24 h (Fold-change > 1; P < 0.05)
Overexpressed GeneFold-change 24 hP-value 24 h
The scatter pattern shows overexpressed and downexpressed genes for 24 hours after infection mock cell (control) and stimulated cell of HT-29 with Cryptosporidium parvum oocysts.
Figure 2. The scatter pattern shows overexpressed and downexpressed genes for 24 hours after infection mock cell (control) and stimulated cell of HT-29 with Cryptosporidium parvum oocysts.
The clustergram pattern shows overexpressed and underexpressed genes related to human type I interferon response in infected cells and mock cells for 6 (A) and 24 hours(B).
Figure 3. The clustergram pattern shows overexpressed and underexpressed genes related to human type I interferon response in infected cells and mock cells for 6 (A) and 24 hours(B).
5. Discussion

It is important to understand whether the reaction of exposed cells to C. parvum cells are resistant or the disease develops. Therefore, it is important to evaluate the innate immune genes related to type I INF response. Limited information is available to describe immune response expression of genes in cells exposed to C. parvum. Therefore, the interest of the current study was to find out overexpression and underexpression of related genes to human type I INFs (α and β) in HT-29 cells exposed to C. parvum for six and 24 hours. Based on fold-change outcomes, four top overexpressed genes with fold-change > 1.0 (P < 0.05) including HLA-A and MX1 genes were evaluated in infected cells for six hours. Birte Pantenburg et al., reported that the HLA-A and -B were involved in presenting the C. parvum antigens to CD8 + T-lymphocytes to eliminate the parasite (8). The current study found that the HLA-A gene was overexpressed in infected cells after six hours exposure. This finding was in agreement with what Birte Pantenburg et al., reported.

IFN-α and -β induce MX1 gene expression that encodes a GTPase activity and acts as a potent inhibitor for replication of many RNA viruses and Salmonella typhimurium (9-11). The current study data showed that MX1 gene was over expressed in the exposed cell six hours post-exposure. It was thought-provoking to express the MX1 gene as an antiviral factor induced by a C. parvum protein, which needs further investigation. Expression of SH2D1A gene synthesizes protein 2B4, which can interact with SLAM (signaling lymphocyte-activating molecule) to monitor T-cell and B-cell response to Epstein-Barr virus (EBV) infection (12-14). The current study findings revealed that SH2D1A gene was overexpressed in vitro for six hours. It is necessary to investigate the characterization of antigen or protein that induces SH2D1A gene expression against C. parvum infection. Wyatt et al., found that C. parvum epitopes (44 to > 200 kD) could prompt overexpression of IL-10 gene in intestinal epithelial cells. They concluded that the expression of IL-10 gene may have an important role in developing parasite infection (15). Based on the current study data, IL-10 gene was overexpressed in the infected cells.

The following genes, TIMP1, TYK2, IRF2, PML, and IRF5 were under expressed in the infected cells for six hours. The expression of TIMP1 (tissue inhibitors of metalloproteinases) gene inhibits expression of MMPs (matrix metalloproteinases) gene, which is implicated in developing cancer, neurological, and infectious pathologies (16). Overexpression of TIMP1 gene is implicated in clearance of parasites infection (16, 17). According to the current study findings, the TIMP1 gene was underexpressed for six hours, but overexpressed in the infected cells for 24 hours. Thus, it was concluded that the overexpression of TIMP1 and SH2D1A genes might be an important key for C. parvum to become self-limited infection.

The expression of TYK2 gene plays a significant role in the intracellular signaling of type I IFNs using phosphorylation and activation of STAT proteins (18). The expression of TYK2 gene plays a vital role in the synthesis of INF-β. The stimulation of INF-β production via overexpression of TYK2 gene is detected in the infected cells by L. major. The current study observed that TYK2 gene was downexpressed for six hours, but overexpressed for 24 hours following the infection. In contrast, Piyal Ganguli stated that leishmaniosis stimulated the production of INF-β, which in turn can upregulate the secretion of RAP1 and SOCS3 proteins within the T-cells that potentially inhibits MAPK and JAK-STAT signaling pathway via TYK2 mediated pathway (19).

The expression of IRF2 (interferon regulatory factor 2) gene upregulated type I INF associated with signaling pathway and resulted in the synthesis of IL-12 in cells infected by L. major (20, 21). IRF2 protein acts multiple functions such as transcription factor, IRF2, transcriptional activator of histone H4, IFN-stimulated regulatory element (ISRE)-responsive genes and requires optimal maturation of natural killer (NK) cells (22). The current study recorded the IRF2 gene as underexpressed for six hours, but it was overexpressed for 24 hours. Thus, it is very important to investigate the characterization of a protein of C. parvum to induce IRF2 gene expression.

Promyelocytic leukemia (PML) gene produces a protein, acts as a tumour suppress and prevents cells from growing and dividing. Promyelocytic leukemia protein is also a cofactor to induce apoptosis. In the current study, the PML gene was downexpressed in the infected cells for six and 24 hours. TLR 7, TLR 9 via signaling pathway of INF type I induced IRF5 gene expression, which activates pro-inflammatory cytokines such as IL6 and TNF (tumor necrosis factor) (23, 24). The expression of IRF5 gene was detected in the infected cells by L. donovani (25). The current study data showed that IRF5 gene was downexpressed in the infected cells for six and 24 hours. The current study selected 10 top overexpressed genes with high enrichment including IL15, DDX58, CXCL10, NMI, MYD88, STAT3, IFNAR2, IFIH1, CASP1, and TLR3 genes relevant to types I IFNs in post-infected cells for 24 hours. IL-15 (mRNA is overexpressed in cells of secretary tissues, dendritic cells (DCs), epithelial cells, bone marrow stromal cells, and fibroblasts) and its functions stimulate the proliferation of T-lymphocytes and important role of IL-15 to control infection (26-28). The current study results indicated that IL-15 gene was overexpressed in the infected cell for 24 hours; to the authors` best knowledge, an overexpression of IL-15 is reported in infected cells by Cryptosporidiosis.

The expression of DDX58 gene can promote immunological responses (29). Alvarez et al., reported that DDX58 gene was downexpressed in lamb liver tissue cells infected by Fasiola hepatica (30). In contrast, the current study found that DDX58 gene was overexpressed in the infected cells for 24 hours. CXCL10 is a Chemokine and can attract the immune cells at the site of infection, while in the patients with AIDS, due to lack of immune cells, CXCL10 is enabled to recruit immune inflammatory cells causing immune pathogenesis. Carneiro et al., reported that the genes related to IFN-γ such as CXCL10 were upregulated when HPs and LPs (high or low IFN-γ producers) following motivation of peripheral blood cells were exposed to L. braziliense, which was in line with the current study findings (31). The expression of CXCL10 gene via stimulation of IFN-δ gene was detected in intestinal epithelial cells infected by C. parvum (32). Type I IFN and IFN-γ stimulated the expression of N-myc interactor (Nmi) gene (33). NMI is a cytokine with a multiple function including enhancement IFN-γ-induced STAT-dependent transcription and increased interleukin-2 (IL-2)-driven STAT5-dependent relative transcription (34). NMI protein acts as antiviral activities against prototype foamy virus (35), Sendai virus (33), and foot and mouth disease virus (FMDV) (36). To the authors` best knowledge, parasite induced NMI protein is not reported yet. In the current study, the NMI gene was overexpressed in the infected cells for 24 hours.

MyD88 gene encodes a protein and implicates in innate immunity. Toll-like receptor 4 induced type I INF via MYD88 (37). Lacroix-Lamande et al., injected polyC to neonatal mice infected with C. parvum, and found enhancement signalling MYD88 pathway, which resulted in resistance of mice intestine cells to parasite oocysts infection (38, 39). Rogers et al., reported that mice lacking MyD88 were more at risk to infection than the wild-type mice (38). In the current study, MYD88 gene was overexpressed for 24 hours indicating an important key role in activation of type I INF. The TLR3 protein performs multiple functions such as pathogen recognition, activation of innate immunity, recognition of dsRNA associated with viral infection, and inducing NF-Кβ and type I INF production (40). The current study observed that the TLR-3 gene was overexpressed for 24 hours in vitro, which was in agreement with the findings of Koblansky et al., (41) but in contrast with those of the study by Beiting (42). STAT3 is a transcription factor that can be activated by interferon ligand that implicates in apoptosis and cell growth processes (43). Kyeong et al. showed that STAT3 gene was downexpressed in Caco-2 cells infected with Entamoeba histolytica, but it was in contrast with the current study results. The STAT3 gene was overexpressed for 24 hours (44).

Caspase-1, which can activate two inflammatory cytokines and initiate the pro-apoptosis program, results in cell death (45). Liu et al., reported that about 20% of infected cells led to apoptosis for 48 hours, but the majority of cells were protected from apoptosis by C. parvum infection. They found high levels of caspase-3, -4, -6, and -7 activation and low levels of caspase-2, -8, and -9 activation, while no caspase-1 and -10 genes were expressed at a 48-hour course of infection. Cryptosporidium parvum suppressed staurosporine-induced apoptosis and initiated caspase-3/7 activation of intestinal epithelial cells via upregulation of survivin, which sustains parasite infection (46). The current study found that the caspase-1 gene was overexpressed for 24 hours. In the current study, caspase-2-10 gene expressions were not evaluated, even though it was necessary for further study.

IFNAR2 (interferon-alpha/beta receptor beta) protein bands activate the receptor stimulation Janus protein kinases, which in turn activate STAT1 and STAT2. It subsequently produces type I IFN and pro-inflammatory cytokines, in that way up regulating a family of IFN-stimulated genes (ISGs) that exert pleiotropic inhibitory effects on viral reproduction in the neighboring cells (47, 48). The overexpression of IFNAR2 gene is yet to be reported in the infected cell by C. parvum. The current study found that the IFNAR2 gene was overexpressed in the infected cells for 24 hours, even though its functions should be investigated.

IFIH1 (interferon-induced helicase 1) recognizes dsRNA of many viruses and acts as mediator for antiviral response (49). The IFIH1 gene was overexpressed, but its function requires further evaluation. The cells infected by C. parvum result in overexpressed genes related to subtype 1 interferon including synthesis of MX1, SH2D1A, TLR3, NMI, and IFIH1 proteins, which may have antiviral properties and are promising new candidate biological drugs against some viral infections. Han et al., reported that anti-viral cytokine interferon-β (IFN-β) was produced after IMs (inflammatory monocytes) exposed to Toxoplasma gondii (50).

6. Conclusions

In conclusion, the overexpression of TIMP1 gene had an important role in the clearance of C. parvum. The overexpression of TIMP1, SH2D1A, and MYD88 genes may lead C. parvum to become a self-limited infection. The overexpressed TIMP1 may have anticancer properties. The overexpression of TYK2 and MYD88 genes synthetizing type I INF may play an important role in innate immunity. It is important to investigate the induction of production the protein and overexpression of IRF2 gene with multiple functions related to subtype I INF. The CXCL10 gene function is involved in the attraction of immune effector cells to the site of infection, the overexpressed CXCL10 was observed. The current study identified the upregulated genes, which can prevent the invasion of parasite to cells. Finally, the identification of all unknown induced overexpressed and underexpressed genes by C. parvum requires thorough investigation. In the light of aforementioned overexpressed genes related to type I IFN as a part of innate immunity may play an important role against C. parvum infection.

1 Farthing MJG. Clinical Aspects of Human Cryptosporidiosis. In: Petry F, editor. Cryptosporidiosis and Microsporidiosis. 6. Karger; 2000. p. 50-74. doi: 10.1159/000060368.
2 Zhou R, Gong AY, Eischeid AN, Chen XM. miR-27b targets KSRP to coordinate TLR4-mediated epithelial defense against Cryptosporidium parvum infection. PLoS Pathog. 2012;8(5). e1002702. doi: 10.1371/journal.ppat.1002702. [PubMed: 22615562]. [PubMed Central: PMC3355088].
3 Barakat FM, McDonald V, Foster GR, Tovey MG, Korbel DS. Cryptosporidium parvum infection rapidly induces a protective innate immune response involving type I interferon. J Infect Dis. 2009;200(10):1548-55. doi: 10.1086/644601. [PubMed: 19821721].
4 Castellanos-Gonzalez A, Cabada MM, Nichols J, Gomez G, White ACJ. Human primary intestinal epithelial cells as an improved in vitro model for Cryptosporidium parvum infection. Infect Immun. 2013;81(6):1996-2001. doi: 10.1128/IAI.01131-12. [PubMed: 23509153]. [PubMed Central: PMC3676030].
5 Karanis P, Aldeyarbi HM. Evolution of Cryptosporidium in vitro culture. Int J Parasitol. 2011;41(12):1231-42. doi: 10.1016/j.ijpara.2011.08.001. [PubMed: 21889507].
6 Park WD, Stegall MD. A meta-analysis of kidney microarray datasets: investigation of cytokine gene detection and correlation with rt-PCR and detection thresholds. BMC Genomics. 2007;8:88. doi: 10.1186/1471-2164-8-88. [PubMed: 17397532]. [PubMed Central: PMC1852103].
7 Meloni BP, Thompson RC. Simplified methods for obtaining purified oocysts from mice and for growing Cryptosporidium parvum in vitro. J Parasitol. 1996;82(5):757-62. [PubMed: 8885885].
8 Pantenburg B, Castellanos-Gonzalez A, Dann SM, Connelly RL, Lewis DE, Ward HD, et al. Human CD8(+) T cells clear Cryptosporidium parvum from infected intestinal epithelial cells. Am J Trop Med Hyg. 2010;82(4):600-7. doi: 10.4269/ajtmh.2010.09-0590. [PubMed: 20348507]. [PubMed Central: PMC2844566].
9 Haller O, Arnheiter H, Lindenmann J, Gresser I. Host gene influences sensitivity to interferon action selectively for influenza virus. Nature. 1980;283(5748):660-2. [PubMed: 7354853].
10 Haller O, Frese M, Kochs G. Mx proteins: mediators of innate resistance to RNA viruses. Rev Sci Tech. 1998;17(1):220-30. [PubMed: 9638812].
11 Shtrichman R, Heithoff DM, Mahan MJ, Samuel CE. Tissue selectivity of interferon-stimulated gene expression in mice infected with Dam(+) versus Dam(-) Salmonella enterica serovar Typhimurium strains. Infect Immun. 2002;70(10):5579-88. [PubMed: 12228285]. [PubMed Central: PMC128359].
12 Nagy N, Maeda A, Bandobashi K, Kis LL, Nishikawa J, Trivedi P, et al. SH2D1A expression in Burkitt lymphoma cells is restricted to EBV positive group I lines and is downregulated in parallel with immunoblastic transformation. Int J Cancer. 2002;100(4):433-40. doi: 10.1002/ijc.10498. [PubMed: 12115526].
13 Sumegi J, Huang D, Lanyi A, Davis JD, Seemayer TA, Maeda A, et al. Correlation of mutations of the SH2D1A gene and epstein-barr virus infection with clinical phenotype and outcome in X-linked lymphoproliferative disease. Blood. 2000;96(9):3118-25. [PubMed: 11049992].
14 Chuang HC, Lay JD, Hsieh WC, Wang HC, Chang Y, Chuang SE, et al. Epstein-Barr virus LMP1 inhibits the expression of SAP gene and upregulates Th1 cytokines in the pathogenesis of hemophagocytic syndrome. Blood. 2005;106(9):3090-6. doi: 10.1182/blood-2005-04-1406. [PubMed: 16002423].
15 Wyatt CR, Barrett WJ, Brackett EJ, Schaefer DA, Riggs MW. Association of IL-10 expression by mucosal lymphocytes with increased expression of Cryptosporidium parvum epitopes in infected epithelium. J Parasitol. 2002;88(2):281-6. doi: 10.1645/0022-3395(2002)088[0281:AOIEBM]2.0.CO;2. [PubMed: 12053998].
16 Clark RT, Nance JP, Noor S, Wilson EH. T-cell production of matrix metalloproteinases and inhibition of parasite clearance by TIMP-1 during chronic Toxoplasma infection in the brain. ASN Neuro. 2011;3(1). e00049. doi: 10.1042/AN20100027. [PubMed: 21434872]. [PubMed Central: PMC3024837].
17 Nasser JA, Falavigna A, Ferraz F, Duigou G, Bruce J. Transcription analysis of TIMP-1 and NM23-H1 genes in glioma cell invasion. Arq Neuropsiquiatr. 2006;64(3B):774-80. [PubMed: 17057884].
18 Babon JJ, Lucet IS, Murphy JM, Nicola NA, Varghese LN. The molecular regulation of Janus kinase (JAK) activation. Biochem J. 2014;462(1):1-13. doi: 10.1042/BJ20140712. [PubMed: 25057888]. [PubMed Central: PMC4112375].
19 Ganguli P, Chowdhury S, Chowdhury S, Sarkar RR. Identification of Th1/Th2 regulatory switch to promote healing response during leishmaniasis: a computational approach. EURASIP J Bioinform Syst Biol. 2015;2015(1):13. doi: 10.1186/s13637-015-0032-7. [PubMed: 26660865]. [PubMed Central: PMC4666900].
20 Favila MA, Geraci NS, Zeng E, Harker B, Condon D, Cotton RN, et al. Human dendritic cells exhibit a pronounced type I IFN signature following Leishmania major infection that is required for IL-12 induction. J Immunol. 2014;192(12):5863-72. doi: 10.4049/jimmunol.1203230. [PubMed: 24808365]. [PubMed Central: PMC4052223].
21 McDowell MA, Marovich M, Lira R, Braun M, Sacks D. Leishmania priming of human dendritic cells for CD40 ligand-induced interleukin-12p70 secretion is strain and species dependent. Infect Immun. 2002;70(8):3994-4001. [PubMed: 12117904]. [PubMed Central: PMC128119].
22 Lohoff M, Duncan GS, Ferrick D, Mittrucker HW, Bischof S, Prechtl S, et al. Deficiency in the transcription factor interferon regulatory factor (IRF)-2 leads to severely compromised development of natural killer and T helper type 1 cells. J Exp Med. 2000;192(3):325-36. [PubMed: 10934221]. [PubMed Central: PMC2193225].
23 Hammami A, Charpentier T, Smans M, Stager S. IRF-5-Mediated Inflammation Limits CD8+ T Cell Expansion by Inducing HIF-1alpha and Impairing Dendritic Cell Functions during Leishmania Infection. PLoS Pathog. 2015;11(6). e1004938. doi: 10.1371/journal.ppat.1004938. [PubMed: 26046638]. [PubMed Central: PMC4457842].
24 Joshi T, Rodriguez S, Perovic V, Cockburn IA, Stager S. B7-H1 blockade increases survival of dysfunctional CD8(+) T cells and confers protection against Leishmania donovani infections. PLoS Pathog. 2009;5(5). e1000431. doi: 10.1371/journal.ppat.1000431. [PubMed: 19436710]. [PubMed Central: PMC2674929].
25 Schoenemeyer A, Barnes BJ, Mancl ME, Latz E, Goutagny N, Pitha PM, et al. The interferon regulatory factor, IRF5, is a central mediator of toll-like receptor 7 signaling. J Biol Chem. 2005;280(17):17005-12. doi: 10.1074/jbc.M412584200. [PubMed: 15695821].
26 Dann SM, Wang HC, Gambarin KJ, Actor JK, Robinson P, Lewis DE, et al. Interleukin-15 activates human natural killer cells to clear the intestinal protozoan cryptosporidium. J Infect Dis. 2005;192(7):1294-302. doi: 10.1086/444393. [PubMed: 16136475].
27 Budagian V, Bulanova E, Paus R, Bulfone-Paus S. IL-15/IL-15 receptor biology: a guided tour through an expanding universe. Cytokine Growth Factor Rev. 2006;17(4):259-80. doi: 10.1016/j.cytogfr.2006.05.001. [PubMed: 16815076].
28 Kennedy MK, Park LPR. Interleukin-15. Thomson; 1998.
29 Imaizumi T, Hatakeyama M, Yamashita K, Yoshida H, Ishikawa A, Taima K, et al. Interferon-gamma induces retinoic acid-inducible gene-I in endothelial cells. Endothelium. 2004;11(3-4):169-73. doi: 10.1080/10623320490512156. [PubMed: 15370293].
30 Alvarez Rojas CA, Ansell BR, Hall RS, Gasser RB, Young ND, Jex AR, et al. Transcriptional analysis identifies key genes involved in metabolism, fibrosis/tissue repair and the immune response against Fasciola hepatica in sheep liver. Parasit Vectors. 2015;8:124. doi: 10.1186/s13071-015-0715-7. [PubMed: 25885344]. [PubMed Central: PMC4382932].
31 Carneiro MW, Fukutani KF, Andrade BB, Curvelo RP, Cristal JR, Carvalho AM, et al. Gene Expression Profile of High IFN-gamma Producers Stimulated with Leishmania braziliensis Identifies Genes Associated with Cutaneous Leishmaniasis. PLoS Negl Trop Dis. 2016;10(11). e0005116. doi: 10.1371/journal.pntd.0005116. [PubMed: 27870860]. [PubMed Central: PMC5117592].
32 Wang HC, Dann SM, Okhuysen PC, Lewis DE, Chappell CL, Adler DG, et al. High levels of CXCL10 are produced by intestinal epithelial cells in AIDS patients with active cryptosporidiosis but not after reconstitution of immunity. Infect Immun. 2007;75(1):481-7. doi: 10.1128/IAI.01237-06. [PubMed: 17043107]. [PubMed Central: PMC1828373].
33 Wang J, Yang B, Hu Y, Zheng Y, Zhou H, Wang Y, et al. Negative regulation of Nmi on virus-triggered type I IFN production by targeting IRF7. J Immunol. 2013;191(6):3393-9. doi: 10.4049/jimmunol.1300740. [PubMed: 23956435].
34 Zhu M, John S, Berg M, Leonard WJ. Functional association of Nmi with Stat5 and Stat1 in IL-2- and IFNgamma-mediated signaling. Cell. 1999;96(1):121-30. [PubMed: 9989503].
35 Hu X, Yang W, Liu R, Geng Y, Qiao W, Tan J. N-Myc interactor inhibits prototype foamy virus by sequestering viral Tas protein in the cytoplasm. J Virol. 2014;88(12):7036-44. doi: 10.1128/JVI.00799-14. [PubMed: 24719420]. [PubMed Central: PMC4054356].
36 Wang J, Wang Y, Liu J, Ding L, Zhang Q, Li X, et al. A critical role of N-myc and STAT interactor (Nmi) in foot-and-mouth disease virus (FMDV) 2C-induced apoptosis. Virus Res. 2012;170(1-2):59-65. doi: 10.1016/j.virusres.2012.08.018. [PubMed: 22974759].
37 Yamamoto M, Sato S, Hemmi H, Hoshino K, Kaisho T, Sanjo H, et al. Role of adaptor TRIF in the MyD88-independent toll-like receptor signaling pathway. Science. 2003;301(5633):640-3. doi: 10.1126/science.1087262. [PubMed: 12855817].
38 Rogers KA, Rogers AB, Leav BA, Sanchez A, Vannier E, Uematsu S, et al. MyD88-dependent pathways mediate resistance to Cryptosporidium parvum infection in mice. Infect Immun. 2006;74(1):549-56. doi: 10.1128/IAI.74.1.549-556.2006. [PubMed: 16369011]. [PubMed Central: PMC1346622].
39 Lacroix-Lamande S, Guesdon W, Drouet F, Potiron L, Lantier L, Laurent F. The gut flora is required for the control of intestinal infection by poly(I:C) administration in neonates. Gut Microbes. 2014;5(4):533-40. doi: 10.4161/gmic.29154. [PubMed: 24918602].
40 Alexopoulou L, Holt AC, Medzhitov R, Flavell RA. Recognition of double-stranded RNA and activation of NF-kappaB by Toll-like receptor 3. Nature. 2001;413(6857):732-8. doi: 10.1038/35099560. [PubMed: 11607032].
41 Koblansky AA, Jankovic D, Oh H, Hieny S, Sungnak W, Mathur R, et al. Recognition of profilin by Toll-like receptor 12 is critical for host resistance to Toxoplasma gondii. Immunity. 2013;38(1):119-30. doi: 10.1016/j.immuni.2012.09.016. [PubMed: 23246311]. [PubMed Central: PMC3601573].
42 Beiting DP, Peixoto L, Akopyants NS, Beverley SM, Wherry EJ, Christian DA, et al. Differential induction of TLR3-dependent innate immune signaling by closely related parasite species. PLoS One. 2014;9(2). e88398. doi: 10.1371/journal.pone.0088398. [PubMed: 24505488]. [PubMed Central: PMC3914978].
43 Yuan ZL, Guan YJ, Wang L, Wei W, Kane AB, Chin YE. Central role of the threonine residue within the p+1 loop of receptor tyrosine kinase in STAT3 constitutive phosphorylation in metastatic cancer cells. Mol Cell Biol. 2004;24(21):9390-400. doi: 10.1128/MCB.24.21.9390-9400.2004. [PubMed: 15485908]. [PubMed Central: PMC522220].
44 Kim KA, Min A, Lee YA, Shin MH. Degradation of the transcription factors NF-kappaB, STAT3, and STAT5 is involved in Entamoeba histolytica-induced cell death in Caco-2 colonic epithelial cells. Korean J Parasitol. 2014;52(5):459-69. doi: 10.3347/kjp.2014.52.5.459. [PubMed: 25352693]. [PubMed Central: PMC4210727].
45 Mariathasan S, Newton K, Monack DM, Vucic D, French DM, Lee WP, et al. Differential activation of the inflammasome by caspase-1 adaptors ASC and Ipaf. Nature. 2004;430(6996):213-8. doi: 10.1038/nature02664. [PubMed: 15190255].
46 Liu J, Enomoto S, Lancto CA, Abrahamsen MS, Rutherford MS. Inhibition of apoptosis in Cryptosporidium parvum-infected intestinal epithelial cells is dependent on survivin. Infect Immun. 2008;76(8):3784-92. doi: 10.1128/IAI.00308-08. [PubMed: 18519556]. [PubMed Central: PMC2493198].
47 Akira S, Uematsu S, Takeuchi O. Pathogen recognition and innate immunity. Cell. 2006;124(4):783-801. doi: 10.1016/j.cell.2006.02.015. [PubMed: 16497588].
48 Beutler B, Eidenschenk C, Crozat K, Imler JL, Takeuchi O, Hoffmann JA, et al. Genetic analysis of resistance to viral infection. Nat Rev Immunol. 2007;7(10):753-66. doi: 10.1038/nri2174. [PubMed: 17893693].
49 Satoh T, Kato H, Kumagai Y, Yoneyama M, Sato S, Matsushita K, et al. LGP2 is a positive regulator of RIG-I- and MDA5-mediated antiviral responses. Proc Natl Acad Sci U S A. 2010;107(4):1512-7. doi: 10.1073/pnas.0912986107. [PubMed: 20080593]. [PubMed Central: PMC2824407].
50 Han SJ, Melichar HJ, Coombes JL, Chan SW, Koshy AA, Boothroyd JC, et al. Internalization and TLR-dependent type I interferon production by monocytes in response to Toxoplasma gondii. Immunol Cell Biol. 2014;92(10):872-81. doi: 10.1038/icb.2014.70. [PubMed: 25155465]. [PubMed Central: PMC4245188].