MicroRNA expression profiles from HEK293 cells expressing H5N1 avian influenza virus non-structural protein 1.

H5N1 avian influenza poses a serious threat to the poultry industry and human health. Non-structural protein 1 (NS1) plays an important role in the replication and pathogenesis of avian influenza virus (AIV). However, the function of the NS1 gene is still unclear. In this study, illumina genome analyzer iix screening was used to identify the differentially expressed microRNAs (miRNAs) in HEK293 cells expressing H5N1 AIV NS1. There were 13 differentially expressed miRNAs (hsa-miR-17-5p, hsa-miR-221-3p, hsa-miR-22-3p, hsa-miR-31-5p, hsa-miR-20a-5p, hsa-miR-222-3p, hsa-miR-24-3p, hsa-miR-3613-3p, hsa-miR-3178, hsa-miR-4505, hsa-miR-345-3p, hsa-miR-3648, and hsa-miR-455-3p) ( P < 0.01). The qRT-PCR validation results demonstrated that hsa-miR-221-3p, hsa-miR-22-3p, hsa-miR-20a-5p, and hsa-miR-3613-3p were upregulated, while hsa-miR-3178 and hsa-miR-4505 were down-regulated. The softwares targetscan and miranda were further used to predict their target genes, and the gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis results showed that 20 GO terms and 20 KEGG pathways were significantly enriched. Our findings are the first to report expression profiling of miRNA and their functions in H5N1 AIV NS1-expressing HEK293 cells, and pave the way to further elucidating the accurate interaction mechanism between NS1 and virus replication, thus providing brand new insight into the prophylaxis and treatment of H5N1 AIV.

Introduction

H5N1 is a subtype of avian influenza (AI). The zoonosis caused by H5N1 avian influenza virus (AIV) has attracted much attention worldwide because it is zoonotic and highly pathogenic. AIV has a special mechanism of infecting mammalian cells; it was reported that AIV virulence is related directly to the avian innate immune system.[1-3]

Non-structural protein 1 (NS1) is found mainly in the nucleus and is synthesized largely in the early stages of AIV infection. NS1 can inhibit the production of IFN in infected cells and plays an important role in the antiviral effect of AIV on IFN.[4] NS1 is an important regulator in the process of replication and infection of AIV. The latest research has shown that NS1 can interact with annexin A2 (ANXA2) to modulate replication of H5N1 AIV.[5] NS1 protein can induce cellular apoptosis via activating reactive oxygen species (ROS) accumulation and mitochondria-mediated apoptotic signaling in chicken oviduct epithelial cells (COECs).[6],[7] These reports have shown that the molecular mechanisms linking NS1 and host are very complex.

MicroRNA (miRNA) is a class of non-coded single-stranded RNA molecules of about 22 nucleotides encoded by endogenous genes, and are involved in the regulation of post-transcriptional gene expression in animals.[8-11] Administration of miR-2911 or a honeysuckle (HS) decoction dramatically reduced mouse mortality caused by H5N1 infection. These results suggested that miR-2911 is involved directly in the process of H5N1 AIV infection.[12] Since recent research has shown that gga-miR-1644 targets muscle blind-like protein 2 (MBNL2), and gga-miR-6675 targets polymerase basic protein 1 (PB1) to regulate the ability of the H9N2 AIV replication, this study may provide brand new insights into the interaction of avian miRNAs with avian DCs and the inhibition of virus replication.[13]

Here, we aimed to reveal the interaction between NS1 and miRNAs in host cells. We reconstructed stable NS1-EGFP and EGFP cell lines. miRNA-Seq-based screening was performed to identify the differential expression miRNA profiling in NS1-expressing HEK293 cells, qRT-PCR was used to validate them, their target genes were predicted by targetscan and miranda software, and the confirmed miRNAs were assigned functions by gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis.

Materials and methods

Construction of NS1-expressing HEK293 cells

The HEK293 cell line was purchased from the Cell Bank of the Chinese Academy of Science (Shanghai, China) and cultured at 5% CO2, 37°C, in DMEM (Thermo Scientific, San Diego, CA) containing 10% heat-inactivated FBS (Thermo Scientific), supplemented with penicillin (100 U/ml) (Thermo Scientific) and streptomycin (100 µg/ml) (Thermo Scientific). Recombinant plasmids pEGFP-N1 and pNS1-EGFP, which express H5N1 AIV NS1-EGFP fusion protein, were transfected into HEK293 cells with lipofectamine™ 3000, and NS1-EGFP and EGFP stable cell lines were constructed, as identified by fluorescence microscopy (Olympus, Tokyo, Japan) observation. FACSCalibur flow cytometer (Becton-Dickinson, Franklin Lakes, NJ) determination and Western blotting detection were performed according to our previous research.[14]

miRNA array screening

Small RNAs were extracted from the NS1-EGFP and EGFP stable cell lines, named NS and E, respectively. Of each group, three independent samples were prepared; a total of six samples were chosen randomly for an initial genome-wide microarray screening. An mParaflo miRNA microarray assay was performed by an external service provider, as previously described.[15] miRNAs that were significantly differentially expressed (P < 0.01) based on the miRNA array were selected for further qRT-PCR validation and functional enrichment analysis.

qRT-PCR validation for miRNA profiling

miRNA profiling was verified using an miRNA specific-TaqMan® probe (Thermo Scientific), and qRT-PCR was used to validate the differentially expressed miRNAs (P < 0.01). U6 snRNA (Thermo Scientific) was used as an internal control. Three independent repeated experiments were performed as described previously.[16]

Prediction of target genes of the confirmed miRNAs

Two online target prediction software programs, targetscan and miranda, were used to predict target genes of miRNAs (P < 0.01). The prediction results were screened according to the grading criteria of each software. Finally, the data predicted by both softwares were combined and overlaps were calculated.

GO and KEGG functional enrichment analysis of miRNAs

According to our previous study, GO and KEGG were performed to analyze the miRNAs.[15] GO annotations were divided into biological process, cellular components, and molecular function categories; each term corresponds to one category. Pathway analysis helps further understand the biological functions of genes. KEGG is the major public database of pathways; pathways were analyzed by significance enrichment. Using the KEGG pathway as a unit, we applied hypergeometric testing to find out the background of the whole genome. Pathways were significantly enriched in significant differentially expressed genes (DEGs).

Statistical analysis

Data were analyzed with the Student’s t-test and one-way ANOVA. P < 0.05 was considered statistically significant, and P < 0.01 was considered statistically highly significant.

Results

miRNA profilings in NS1-expressing HEK293 cells

NS1-EGFP and EGFP stable cell lines were constucted according to our previous research,[14] and expression of EGFP and NS1-EGFP fusion protein were detected by fluorescence microscopy and Western blotting. The percentage of positive rate and the fluorescent mean value of the stable cells lines were detected by flow cytometry, as shown in Supplementary Figure 1.

The small RNAs of NS and E were extracted from the NS1-EGFP and EGFP stable cell lines. High-throughput based microarray screening was performed. As shown in Figure 1, there were 3 differentially expressed miRNAs (P < 0.01), namely hsa-miR-17-5p, hsa-miR-221-3p, hsa-miR-22-3p, hsa-miR-31-5p, hsa-miR-20a-5p, hsa-miR-222-3p, hsa-miR-24-3p, hsa-miR-3613-3p, hsa-miR-3178, hsa-miR-4505, hsa-miR-345-3p, hsa-miR-3648, and hsa-miR-455-3p. The P values, group E means, group NS means and log2 (NS/E) values of these are shown in Table 1.

Figure 1.

Heatmap of differentially expressed miRNAs in NS1-expressing HEK293 cells (P < 0.01). NS, NS1-EGFP stable cell line; E EGFP stable cell line. NS-1, NS-2, and NS-3 represent the samples of three independently repeated experiments, respectively, as do E1, E2 and E3. The red color part represents upregulated expression, and green color, down-regulated expression.

Table 1.

The significantly differentially expressed miRNAs identified by microarray analysis.

		Group1E	Group2NS
Reporter name	P-Value	Mean	Mean	Log2 (G2/G1)
hsa-miR-17-5p	6.08E-03	3,034	4,955	0.71
hsa-miR-221-3p	7.94E-03	1,654	3,421	1.05
hsa-miR-22-3p	9.68E-03	184	594	1.69
hsa-miR-31-5p	2.61E-03	418	655	0.65
hsa-miR-20a-5p	1.39E-03	3,106	4,606	0.57
hsa-miR-222-3p	2.11E-03	2,625	5,065	0.95
hsa-miR-24-3p	6.86E-03	756	1,016	0.43
hsa-miR-3613-3p	7.19E-03	598	1233	1.04
hsa-miR-3178	4.05E-03	11,642	7,211	–0.69
hsa-miR-4505	7.43E-03	1,124	509	–1.14
hsa-miR-345-5p	4.72E-03	55	30	–0.86
hsa-miR-3648	3.03E-03	76	19	–1.98
hsa-miR-455-3p	2.07E-03	174	76	–1.19

qRT-PCR validation of the 13 differentially expressed miRNAs

To validate the miRNA profiling, miRNA specific TaqMan® probes of the 13 differentially expressed miRNAs were obtained, and qRT-PCR was performed to confirm the miRNA array-based data. U6 snRNA was used as an internal control. As shown in Figure 2, when compared with the E group, the qRT-PCR validation results showed that hsa-miR-221-3p, hsa-miR-22-3p, hsa-miR-20a-5p, and hsa-miR-3613-3p were upregulated, and hsa-miR-3178 and hsa-miR-4505 were down-regulated in NS1-expressing HEK293 cells. The qRT-PCR validation ΔΔCт values of the E and NS groups are shown in Table 2.

Figure 2.

qRT-PCR validation of differentially expressed miRNA profiling obtained from the miRNA microarray assay. Up-regulation of hsa-miR-221-3p, hsa-miR-22-3p, hsa-miR-20a-5p, and hsa-miR-3613-3p was confirmed, as well as down-regulation of hsa-miR-3178 and hsa-miR-4505.

Table 2.

The qRT-PCR validation ΔΔCт value of the E and NS groups.

Reporter name	ΔΔCт
	Group1E	Group2NS
hsa-miR-221-3p	0	1.02
hsa-miR-22-3p	0	1.65
hsa-miR-20a-5p	0	0.60
hsa-miR-3613-3p	0	1.08
hsa-miR-3178	0	–0.62
hsa-miR-4505	0	–1.11

Prediction of target genes of the six differentially expressed miRNAs and GO functional enrichment analysis

The targetscan and miranda software programs were used to predict the target genes of the six differentially expressed miRNAs. In the targetscan algorithm, target genes with a context score percentile of < 50 were removed. In the miranda algorithm, target genes with maximum free energy (Max Energy) greater than –10 were removed. Finally, the overlap of the two software programs was taken to represent the final target genes of the six differentially expressed miRNAs.

As shown in Figure 3, the target genes of the six confirmed miRNAs were assigned to three categories: biological process, cellular component, and molecular function (P < 0.05), and the percentage of the genes of 50 GO terms were determined (Supplementary Table S1).

Figure 3.

GO classifications for the predicted target genes of the six differentially expressed miRNAs.

Further, a functional enrichment analysis of all GO terms was performed. As shown in Figure 4, 20 GO terms were significantly enriched, and the target genes of these are shown in Supplementary Table S2. The enrichment factor indicated the number of differential genes located in the GO relative to the total GO base, and the smaller the P value, the higher the concentration of GO.

Figure 4.

GO functional enrichment analysis of the predicted target genes of the six differentially expressed miRNAs. y-axis: GO terms; x-axis: rich factor. The color of each bubble represents the P value, and bubble size the gene number.

Prediction of target genes of the six differentially expressed miRNAs and KEGG functional enrichment analysis

Target genes were predicted as above. As shown in Figure 5, functional enrichment analysis of the pathway was performed. Twenty KEGG pathways were significantly enriched, the target genes of which are shown in Supplementary Table S3. The enrichment factor indicated the number of differential genes located in the KEGG relative to the total base of the KEGG, and, the smaller the P value, the higher the concentration of KEGG.

Figure 5.

KEGG functional enrichment analysis of the predicted target genes of the six differentially expressed miRNAs. y-axis: pathway terms; x-axis: rich factor. The color of each bubble represents the P value and bubble size the gene number.

Discussion

H5N1 is the most common AIV and represents an important pathogen for zoonosis. New vaccines against highly pathogenic AIV have been strongly discussed worldwide.[17-20] In the process of AIV infection in birds, NS1 proteins mediate the translation and replication of viral proteins by affecting host signaling pathways. There were reports that the poly ADP-ribose polymerase 10 (PARP10) protein was synergistic with H5N1 subtype AIV NS1 protein and co-regulated the replication process of AIV virus.[21]

Genome microarray experiments were used to identify the DEGs in chicken embryo fibroblasts (CEF) stimulated by the highly pathogenic H5N1 AIV, providing new insights into how AIV affects expression of host genes, as well as the host immune response and viral replication.[22] Our study was performed to reveal the molecular mechanism of AIV-infected HEK293 target cells, and to identify how the differentially expressed miRNAs affect the host response and replicate in host cells. Stable E and NS cell lines were constructed, and small RNAs extracted. An initial genome-wide microarray screening was performed to identify the miRNA profiling, and qRT-PCR validation results showed that there were six differentially expressed miRNAs, consistent with the results of the array-based chip. However, there were seven false positive miRNAs that could not be validated by qRT-PCR. The reason for this might be that the transcripts were statistically significant but had low signals (< 500) and the mean values were too low to confirm.

The hsa-miR-221-3p plays an important role in the diagnosis of breast cancer (BC) and other tumors. The interaction between hsa-miRNA-221-3p and its target gene may serve as a molecular marker for the diagnosis of BC, and as a target for subsequent treatment.[23] Long non-coding RNAs (lncRNAs) play an important role in many cancers and tumor pathogenesis, e.g., lncRNA digeorge syndrome critical region gene 5 (DGCR5) inhibits expression of hsa-miR-22-3p, which promotes the development of lung adenocarcinoma (LUAD).[24] LncRNAs and hsa-miR-20a-5p co-target the target gene transcription 1, which mediates hepatocellular carcinoma (HCC) growth, cell cycle progression, and the migration process.[25] In HCC, hsa-miR-3613-3p affects cell proliferation and the cell cycle, and is considered as a biomarker and therapeutic target for liver cancer.[26] The hsa-miRNA-3178 regulates expression of the target gene early growth responsive gene 3 (EGR3), thereby mediating the proliferation, invasion, angiogenesis and apoptosis of HCC tumor endothelial cells (TECs).[27] No research on hsa-miR-4505 has been reported yet.

miRNA microarray chip analysis was performed to investigate the functions of H5N1 AIV NS1. Each sample of total RNA first passed quality control, and the whole chip experiment process, included labeling, chip hybridization, hybridization signal amplification, and image acquisition, was then carried out. Six miRNAs were identified in the array-based screen. Up-regulation of hsa-miR-221-3p, hsa-miR-22-3p, hsa-miR-20a-5p, and hsa-miR-3613-3p, and down-regulation of hsa-miR-3178 and hsa-miR-4505 were validated, mediated by a stable expression of H5N1 AIV NS1 protein in HEK293 cells. Bioinformatic software tools were used to predict the target genes of the six confirmed miRNAs and to determine their functional enrichment. Transfection of miRNA mimics, miRNA inhibitors and 3’UTR dual-luciferase reporter gene experiments will be performed further to identify the target genes of the sex miRNAs, which may help elucidate the regulatory relationship between H5N1 AIV NS1 and the host.

This study is the first to report expression profiling of miRNAs and determining their functions in HEK293 cells stably expressing NS1 protein of H5N1 AIV. However, future research should focus not only on how to elucidate the accurate interaction mechanism between NS1 and virus replication, but also on how to reveal NS1 affecting the host immune response. Our findings provide new insights into possible routes for the prophylaxis and treatment of H5N1 AIV.

Click here for additional data file.

Supplemental material, Supplemental Figure for MicroRNA expression profiles from HEK293 cells expressing H5N1 avian influenza virus non-structural protein 1 by Hanwei Jiao, Zonglin Zheng, Xuehong Shuai, Li Wu, Jixuan Chen, Yichen Luo, Yu Zhao, Hongjun Wang and Qingzhou Huang in Innate Immunity

Click here for additional data file.

Supplemental material, Supplemental Table1 for MicroRNA expression profiles from HEK293 cells expressing H5N1 avian influenza virus non-structural protein 1 by Hanwei Jiao, Zonglin Zheng, Xuehong Shuai, Li Wu, Jixuan Chen, Yichen Luo, Yu Zhao, Hongjun Wang and Qingzhou Huang in Innate Immunity

Click here for additional data file.

Supplemental material, Supplemental Table2 for MicroRNA expression profiles from HEK293 cells expressing H5N1 avian influenza virus non-structural protein 1 by Hanwei Jiao, Zonglin Zheng, Xuehong Shuai, Li Wu, Jixuan Chen, Yichen Luo, Yu Zhao, Hongjun Wang and Qingzhou Huang in Innate Immunity

Click here for additional data file.

Supplemental material, Supplemental Table3 for MicroRNA expression profiles from HEK293 cells expressing H5N1 avian influenza virus non-structural protein 1 by Hanwei Jiao, Zonglin Zheng, Xuehong Shuai, Li Wu, Jixuan Chen, Yichen Luo, Yu Zhao, Hongjun Wang and Qingzhou Huang in Innate Immunity