ZnO NPs treatment induces an extensive transcriptional response
The high throughput RNA-sequencing was performed on the 143B and MG-63 OS cells treated with or without ZnO NPs. Then, total RNA were extracted from different groups of cells and pooled for RNA-seq. All of the sequencing reads were aligned to the NCBI human reference genome. To observe the gene expression patterns, cluster analysis of all the DEGs based on the log10 (FPKM+1) was conducted with heat maps for the 143B cells (Figure 1A) and MG-63 cells (Figure 1B). In 143B cells, 2544 significantly DEGs were identified in ZnO NPs treated OS cells compared with control, which consisted of 1628 up- and 916 down-regulated genes.In MG-63 cells, the results obtained 2768 significantly DEGs in two groups treated like 143B cell, Among the 2768 DEGs, the expression of 1756 were increased upon ZnO NPs intervention, and 1012 genes were down-regulated (Figure 1C). Of all the genes that differently regulated between the two cells, 928 overlapping genes existed between these two cellular lines (Figure 1D). These data demonstrated that ZnO NPs treatment induces an extensive transcriptional response in OS cells. The full list of DEGs is shown in Supplementary Table S1 and S2.
Functional gene ontology (GO) classification analysis of DEGs upon ZnO NPs treatment
To analysis DEGs and obtain the biological functional influence of ZnO NPs treatment on OS cells, we employed GO analysis to evaluate the cellular components, biological processes, and molecular functions correlated with DEGs between non-ZnO NPs and ZnO NPs-treated groups. The functions of three structured vocabularies were significant differences between the two cell lines (Figure 2A, 2B). In 143B cell line, we found functions of biological processes (BP) focused on tissue development- and interleukin-7-mediated signaling pathway. In MG-63 cells, however, functions mainly focused on the innate immune pathway and cell adhesion. In the cellular components (CC) aspect, functions concentrated on extracellular space in both cellular lines. In specific cells, nucleosome assembly and Wnt signalosome were focus on 143B cells, whereas integrin complex and chromatoid body were specific located in MG-63 cells. In the last section of molecular function (MF), the nucleic acid binding related pathways were enriched in both cell lines. However, hormone activity was enriched in 143B cell line while specific enrichment in MAP kinase tyrosine phosphatase pathway in MG-63 cell line.These data suggested that ZnO NPs treatment triggers multiples responses in OS cell functions.
Functional KEGG: Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of DEGs upon ZnO NPs treatment
To get further understanding of the metabolic and signal transduction pathways involved in ZnO NPs treatment, KEGG analysis was conducted between non-ZnO NPs and ZnO NPs-treated groups in two OS cell lines. In both 143B and MG-63 cells, total of seven pathways were significantly enriched in both cell lines, including MAPK, Toll-like receptor (TLRs) and NF-kappa B signaling pathways which were related to tumor growth and inflammation functions (Figure 3A, 3B). Moreover, the analysis displayed that multiple pathways were different in these two cellular lines. The TNF and p53 signaling pathways were enriched in 143B cells (Figure 3A). The PI3K-AKT and RIG-like receptor signaling pathways were specificly abundant in MG-63 cells (Figure 3B). Interestingly, the result also revealed the alcoholism regulated pathway was enriched in 143B cells. However, in MG-63 cells, the changed microRNA related pathways were enriched. These KEGG terms exhibited some pathway similarity in both cell lines as well as certain cell-specific pathways.
PPI networks of DEGs between non-ZnO NPs and ZnO NPs-treated OS cellular lines
The protein–protein interaction (PPI) analysis provided a novel approach for the discovery of novel gene interaction and identify neighbors of known mediators. To investigate the PPI network in ZnO NPs treated OS cell lines, we used mapping tool of NetworkAnalyst3.0)/Cytoscape to conduct PPI analysis (Figure 4A) with the top 20 genes classified as hub genes with the greatest differences in expression of each cellular lines. In the PPI network, the most significantly proteins were HMOX1 and MAFB in 143B cell line. On the other hand, CXCL10 and CXCL11 proteins were the most significantly regulated hub proteins in MG-63 cells (Figure 4B). These proteins are involved in multiple aspects of biological events in OS cells during ZnO NPs treatment.
Verification of the expression level of DEGs by RT-qPCR
To verify the DEGs data of RNA-seq analysis, the RT-qPCR assay was conducted to verify the expression levels of 8 genes displayed the most significant changes upon ZnO NPs treatment in two cellular lines (Table 1). As shown in Figure 5A, RT-qPCR verified that MT1M, HSPA6, RASD1 and CLDN6 were significantly increased while HOXA6, SERTM2 and TNFRSF10D were significantly downregulated in ZnO NPs-treated 143B cells. In Figure 5B, HSPA6, XIRP1, C9orf152 and HMOX1 were significantly up-regulated while TNFSF10, CXCL10, SLFN12L and CHL1 obviously downregulated. Taken together, these results indicated that data of most DEGs obtained from RNA-seq are well identified by RT-qPCR detection.
Verification of ZnO NPs stimulation of pathways by Western Blot
Based on the KEGG results from RNA-seq, multiple signaling pathways were influenced by ZnO NPs treatment in both 143B and MG-63 cellular lines. We choose to verify pathways important for cancer biological, cancer treatment and drug resistance aspects, including MAPK, Toll-like receptor, NF-κB pathway. As shown in Figure 6A, the MAPK signaling pathway displayed overall repression, as the expression of p-ERK, p-p38 MAPK and p-c-Jun were all downregulated to some degree upon ZnO NPs treatment. Next, TLR4, MyD88 and TRAF6 which are the key proteins of the Toll-like receptor pathway were significantly upregulated upon ZnO NPs treatment in both cellular lines (Figure 6B). Similarly, the NF-κB pathway showed a significant inhibitory characterized by p-IκBα and p-p65 upon ZnO NPs treatment (Figure 6C). These above changes resembled responses to multiple aspects of signaling pathway in ZnO NPs-treated OS cells, which further disturbed cellular functions.