Innate Immune and Inflammatory Responses to Respiratory VirusesView this Special Issue
Research Article | Open Access
MicroRNA Expression Profile of Whole Blood Is Altered in Adenovirus-Infected Pneumonia Children
Human adenovirus (Adv) infection is responsible for most community-acquired pneumonia in infants and children, which results in significant morbidity and mortality in children every year. MicroRNAs (miRNAs) are associated with viral replication and host immune response. Knowing the miRNA expression profile will help understand the role of miRNAs in modulating the host response to adenovirus infection and possibly improve the diagnosis of adenovirus-infected pneumonia. In our study, total RNA extracted from whole blood of adenovirus-infected pneumonia children and healthy controls were analyzed by small RNA deep sequencing. Expression profiles of whole blood microRNAs were altered and distinctly different in adenovirus-infected children. The top 3 upregulated miRNA (hsa-miR-127-3p, hsa-miR-493-5p, and hsa-miR-409-3p) were identified in adenovirus-infected children and provided a clear distinction between infected and healthy individuals. Potential host target genes were predicated and validated by qRT-PCR to study the impact of microRNAs on the host genes. Most of the target genes were involved in the MAPK signaling pathway and innate immune response. These highly upregulated microRNAs may have crucial roles in Adv pathogenesis and are potential biomarkers for adenovirus-infected pneumonia.
Human adenovirus (Adv) infection is responsible for most community-acquired pneumonia in infants and children [1, 2]. Adv causes infections for 5–10% of upper and lower respiratory tract infections in children, which results in pneumonia and nearly 1.3 million deaths of children every year [3, 4]. The fatality rates for untreated severe pneumonia or disseminated disease caused by Adv may even exceed to 50% [5, 6]. There are no efficacious antiviral drugs for Adv treatment until now. Also, the traditional diagnosis of Adv infection is limited. Therefore, to discover the interaction between the virus and its host will help us to find novel treatment and diagnosis for Adv infection.
Human Advs are nonenveloped double-stranded DNA viruses and belong to the Adenoviridae family [7, 8]. Human Advs are divided into seven subgroups including 53 serotypes based on immunologic and biological characteristics. Adv replicates efficiently in human cells and triggers an innate immune response such as inflammatory response in the host cells. Also, viral infection has proved to have a great impact on cellular small RNA expression and gene expression [9–11]. Adv infects the host cells through binding to different cellular receptors such as coxsackievirus and adenovirus receptor (CAR) . During Adv DNA replication, host cellular proteins such as nuclear factor I and POU2F1 are used by Adv [13, 14]. In turn, the host will trigger an innate immune response against Adv infection. However, the details of Adv and host interaction still remains poorly defined.
Small RNAs are important regulators that modulate development, proliferation, differentiation, and apoptosis of organisms [15, 16]. Small RNAs include microRNA (miRNA), siRNA, tRFs, piRNA, and rasiRNAs, which regulate gene expression in a wide range of processes such as viral replication and host immune response. miRNAs are the most well-studied small RNAs during the latest decades. miRNAs are very important regulators that modulate transcriptome changes . miRNAs regulate gene expression in a wide range of physiological and pathological processes such as in immune response and viral replication . Although miRNAs have been examined in Adv type 3-infected human laryngeal epithelial cells and Adv type 2-infected human lung fibroblast cells [19, 20], there is no research about small RNA profiling in whole blood of Adv-infected pneumonia children. In our study, we sought to present the different miRNA profiles between Adv-infected pneumonia children and healthy controls, identify candidate diagnostic biomarkers for pneumonia with Adv infection in children, and examine the role of miRNAs in host defense response in Adv-infected children.
2. Materials and Methods
The whole blood samples used in the study were obtained from Guangzhou Women and Children’s Medical Center. Children diagnosed with human Adv pneumonia were included in the study. The diagnosis of human Adv pneumonia was considered certain when it was associated with the following criteria: (1) lower respiratory and/or systemic symptoms, (2) lung infiltration on chest radiography or computed tomography (CT) scan, and (3) positive results for human Adv IgM antibody in sera and/or human Adv DNA by PCR in throat swabs and/or bronchoalveolar lavage (BAL) fluid. A total of 33 samples from patients and 33 samples from healthy volunteers were used in the study. The ages of all patients (male or female) and healthy volunteers (male or female) range from one year to three years. The study was approved by the Ethics Committee at Guangzhou Women and Children’s Medical Center (number 2014121815), and written informed consent was obtained from all guardians.
2.2. RNA Extraction and Small RNA Sequencing
The blood samples of the patients and volunteers were collected in anticoagulant tubes. The total RNA was isolated using RiboPure™ Blood RNA Isolation Kit (Ambion, USA) according to the manufacturer’s protocol. The extracted small RNAs were treated with DNase before deep sequencing. RNA concentration was determined using a NanoDrop ND1000 system (Thermo Fisher Scientific, South San Francisco, CA), and their integrity was verified using an Agilent 2100 Bioanalyzer (Agilent Technologies, USA). Small RNA deep sequencing was performed as previously described . A total of 3 samples from patients and 3 samples from healthy volunteers were used for small RNA deep sequencing in the study. The clinical characteristics of the 3 patients and 3 healthy volunteers are summarized in Table 1.
aFever duration from onset to normothermia. bData extracted from the first test for the children on admission. cSpecific IgM antibodies against HAdV pneumoniae were detected in 2 mL of acute phase (on admission) patient serum using a commercial ELISA kit. dHAdV pneumonia DNA was detected by real-time PCR in throat swabs or bronchoalveolar lavage fluid (BALF). eJudged by chest radiograph or CT scan in whole course of the patients.
2.3. Cell and Virus Culture
Human primary lung fibroblasts (IMR-90) and human 293T cells were grown in Dulbecco’s modified Eagle’s medium (DMEM) (HyClone) supplemented with 10% fetal bovine serum (Invitrogen), streptomycin, and penicillin (Invitrogen). All virus infections were carried out in serum-free medium for 1 h, followed by addition of saved complete medium.
miRNA inhibitors (GenePharma Co., China) were transfected into IMR-90 cells with Lipofectamine RNAiMAX (Invitrogen) at final concentrations of 100 nM according to the manufacturer’s protocol.
2.5. Virus Growth Assay
IMR-90 cells transfected with miRNA inhibitors were infected with HAdV5 (isolated from patients) at an MOI of 10 in serum-free medium. Virus titers were determined 72 h after infection by plaque assays performed on 293T cells [22, 23].
2.6. Quantitative Real-Time PCR (qRT-PCR) Analysis
A qRT-PCR experiment was performed using the Power SYBR Green PCR Master Mix. Each reaction was performed in a 10 μL volume system containing 0.5 μL of cDNA, 0.5 μL of each primer, 5 μL of Power SYBR Green PCR Master Mix, and 3.5 μL of ddH2O. The reactions were incubated in a 96-well plate at 95°C for 10 minutes, followed by 39 cycles of 95°C for 15 seconds and 62°C for 1 minute. For miRNA quantitation, 10 ng of total RNA was reverse-transcribed using specific stem-loop primers. U6 was used as an endogenous control.
2.7. Data Analysis
miRNA cluster and family information from miRBase (miRBase 20, http://www.mirbase.org/) was used to annotate the cluster of miRNAs. miRNA target genes were predicted based on two software: miRDB software (http://mirdb.org/) and TargetScan software (http://www.targetscan.org/).
2.8. Statistical Analysis
Data were analyzed using GraphPad Prism 6.0 software (La Jolla, CA, USA). The two-tailed Student’s -test was used to determine the significance of statistical data between two experimental groups. Data were considered significant at , , and .
3.1. Different Expressions of MicroRNAs in Adv-Infected Children vs. Healthy Children
To study the impact of Adv infection on cellular small RNA expression in pneumonia children, deep sequencing of small RNAs was performed in our study. Accordingly, we found an apparent small RNA peak at 21–24 nt for miRNAs (Figure 1(a)). When we analyzed the differently expressed small RNAs, 118 miRNAs were found differently expressed in Adv-infected children vs. healthy children in a volcano plot (Figure 1(b)). We mapped the clean reads from each group to the known miRNA sequences and identified 908 miRNAs in Adv-infected children versus healthy controls (Figure 1(c)).
Furthermore, hierarchical cluster analysis of differentially expressed miRNA is shown in Figure 2. Specially, 77 differentially expressed miRNAs in the 6 samples passed our fold-change filter (log2 fold ), among which 20 miRNAs that have high expression (reads up to 1000 in Adv samples) and showed significant different expressions were selected for further analysis (Table 2).
3.2. Validation of Differentially Expressed miRNAs
To confirm the differential expression of miRNAs in Adv patients vs. healthy controls, we performed qRT-PCR assays in our study. The result showed that hsa-miR-127-3p, hsa-miR-379-5p, hsa-miR-493-5p, hsa-miR-409-3p, hsa-miR-99b-5p, hsa-miR-370-3p, and hsa-miR-381-3p were upregulated in whole blood samples from 5 Adv-infected children vs. 5 healthy controls (Figure 3(a)), while hsa-miR-101-3p, hsa-miR-150-5p, hsa-miR-29a-3p, and hsa-miR-342-3p were downregulated in whole blood samples from 5 Adv-infected children vs. 5 healthy controls (Figure 3(a)), which was comparable with our sequencing data. To identify candidate diagnostic miRNA biomarkers, we focused on the upregulated miRNAs. More samples (15 Adv-infected children vs. 15 healthy controls) were collected to verify the expression of the upregulated miRNAs (Figure 3(b)). From the result, we found that hsa-miR-127-3p, hsa-miR-493-5p, and hsa-miR-409-3p were significantly increased. Also, after depleting these miRNAs with miRNA inhibitors, the viral replications were significantly decreased (Figure 3(c)). These findings imply that our selected miRNAs may reflect the infection of Adv, and such miRNAs can likely serve as biomarker candidates for Adv-infected patients.
3.3. Predict Target Genes of Differentially Expressed miRNAs
To study the biological significance of miRNAs, we then predicted the computational target genes of miRNAs. We focused on the targets of verified miRNAs with distinct expression profiles. Go enrichment of the predicted target genes of the miRNAs showed that the target genes were mainly involved in cellular process and molecular function (Figure 4(a)). In particular, most of the target genes participated in the MAPK signaling pathway and Ras signaling pathway (Figure 4(b)). The top five predicted target genes of hsa-miR-127-3p, hsa-miR-493-5p, or hsa-miR-409-3p are listed in Table 3 by the highest miRNA target score from two miRNA predictive software.
The target genes were selected for further validation with qRT-PCR, and we found that the mRNA expression of 8 genes (PSMB5, ITGA6, MYCBP2, TCF7L2, UBE2V2, HIPK1, UBE2D2, and KANSL1) were downregulated in Adv-infected patients compared to healthy controls (Figure 5). Most of them are transcription factors or factors involved in the ubiquitin pathway. The downregulated mRNAs may indicate the mechanism of Adv infection and Adv-induced pneumonia.
Human Advs are common causative pathogens of acute respiratory infections in children. The treatment of human Adv infections is limited because prospective, randomized therapeutic trials have not been done. Therefore, it is very important to discover the mechanism of Adv-infected pneumonia and search the biomarkers for Adv-infected pneumonia in children. The regulatory potential of miRNA is well defined, and the different profiles of miRNA expression are the result of diverse diseases including viral infections. In our study, we found that the miRNA profile in whole blood of Adv-infected children was different from that of healthy children. Blood samples of Adv-infected children reflect the associated pathology of Adv infection and thus provide a better understanding of the disease.
Profiling of miRNA expression from Adv-infected blood samples identified a cluster of 118 miRNAs significantly altered. The altered blood miRNA profile was similar to cells infected with Adv reported previously [19, 20], indicating that the different expressed miRNAs identified in our study could be taken as the candidate diagnostic biomarkers for pneumonia with Adv infection in children. Among these different expressed miRNAs, hsa-let-7e-5p was reported to involve in the replication of influenza infection . hsa-miR-127-3p was reported to affect the Epstein-Barr virus-associated lymphoma through targeting the PTEN-AKT-mTOR pathway . Altogether, these miRNAs play important roles in innate immune response or viral replication, thus affecting the outcome of the disease.
Through miRNA target gene GO analysis, we found that most target genes of different expressions of miRNAs were involved in the MAPK signaling pathway. The MAPK signaling pathway is activated by Toll-like receptors , which play important roles in innate immune response against viral infection. Also, the activation of the MAPK pathway will result in the activation of the NF-κB signal pathway, which will stimulate the production of inflammatory cytokines and MMPs  and finally lead to pneumonia. Thus, the altered miRNA expression profile of whole blood from Adv-infected children partly reflected the mechanism of Adv-infected pneumonia.
With a developmental framework or disease process, miRNAs exhibit dynamic expression patterns. In our study, we characterized the miRNA expression profile of Adv-infected children using deep sequencing analysis and identified that many miRNAs were differently expressed in Adv-infected children when compared with healthy children. Those with the greatest differences were chosen for further verification. In particular, the expressions of hsa-miR-127-3p, hsa-miR-493-5p, and hsa-miR-409-3p from 20 Adv-infected children were significantly higher than those from 20 healthy controls, indicating that these miRNAs could be taken as good diagnostic biomarkers for Adv-infected pneumonia.
To further explore the possible molecular mechanisms of the differently expressed miRNAs in Adv-infected children, we predicted the possible target genes of hsa-miR-127-3p, hsa-miR-493-5p, and hsa-miR-409-3p and found that most of them are transcription factors or factors involved in the ubiquitin pathway. Especially, after verifying the predicted target genes with qRT-PCR, 8 genes (PSMB5, ITGA6, MYCBP2, TCF7L2, UBE2V2, HIPK1, UBE2D2, and KANSL1) were found significantly downregulated in samples from Adv-infected children. In particular, MYCBP2 is an E3 ubiquitin protein ligase [28, 29], which regulated the cAMP and mTOR signaling pathway. The mTOR signaling pathway plays a critical role in effector T cell function, the downregulation of which will result in impaired cell cytolysis and incapability of virus elimination. UBE2V2 is an ubiquitin-conjugating enzyme and is involved in the differentiation of monocytes, which can produce proinflammatory cytokines like MIP-1β [30, 31]. Thus, the downregulated UBE2V2 will result in abnormal differentiation of monocytes and more proinflammatory will be produced, which in turn aggravates pneumonia. UBE2D2 is essential for the activation of MAVS and RIG-I in response to viral infection [32–34]. The downregulation of UBE2D2 will result in abnormal activation of MAVS and RIG-I signal and virus clearance. These findings imply that the miRNA expression profile changed in Adv-infected children resulting in different transcriptome profiles, which reflects the mechanism of Adv replication and the formation Adv-infected pneumonia.
In summary, we identified the 3 most markedly differently expressed miRNAs in whole blood from Adv-infected children, which can be taken as biomarkers for Adv-infected pneumonia. Simultaneously, based on target gene prediction and qRT-PCR analysis, we found that genes MYCBP2, UBE2V2, and UBE2D2 may play important roles in viral replication and Adv-induced pneumonia. However, additional studies are necessary to clarify their roles in these processes, which will provide a physiological basis for the treatment of Adv-infected pneumonia.
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare no conflict of interest.
Feng Huang and Junsong Zhang contributed equally to this work.
This study was supported by the National Natural Science Foundation of China (81701990, 81702327, and 81601759), the National Science and Technology Major Project (2018ZX10101004003001), and the General Projects of Guangzhou Scientific Research (201707010183).
- C. M. Nascimento-Carvalho, “Etiology of childhood community acquired pneumonia and its implications for vaccination,” The Brazilian Journal of Infectious Diseases, vol. 5, no. 2, pp. 87–97, 2001.
- S. J. Kim, K. Kim, S. B. Park, D. J. Hong, and B. W. Jhun, “Outcomes of early administration of cidofovir in non-immunocompromised patients with severe adenovirus pneumonia,” PLoS One, vol. 10, no. 4, article e0122642, 2015.
- Z. A. Bhutta, J. K. Das, N. Walker et al., “Interventions to address deaths from childhood pneumonia and diarrhoea equitably: what works and at what cost?” The Lancet, vol. 381, no. 9875, pp. 1417–1429, 2013.
- A. Agweyu, M. Kibore, L. Digolo et al., “Prevalence and correlates of treatment failure among Kenyan children hospitalised with severe community-acquired pneumonia: a prospective study of the clinical effectiveness of WHO pneumonia case management guidelines,” Tropical Medicine & International Health, vol. 19, no. 11, pp. 1310–1320, 2014.
- J. P. Lynch III and A. E. Kajon, “Adenovirus: epidemiology, global spread of novel serotypes, and advances in treatment and prevention,” Seminars in Respiratory and Critical Care Medicine, vol. 37, no. 4, pp. 586–602, 2016.
- U. Sandkovsky, L. Vargas, and D. F. Florescu, “Adenovirus: current epidemiology and emerging approaches to prevention and treatment,” Current Infectious Disease Reports, vol. 16, no. 8, p. 416, 2014.
- M. G. Ison and R. T. Hayden, “Adenovirus,” Microbiology Spectrum, vol. 4, 2016.
- J. P. Lynch, M. Fishbein, and M. Echavarria, “Adenovirus,” Seminars in Respiratory and Critical Care Medicine, vol. 32, no. 4, pp. 494–511, 2011.
- K. L. Wu, X. Zhang, J. Zhang et al., “Inhibition of hepatitis B virus gene expression by single and dual small interfering RNA treatment,” Virus Research, vol. 112, no. 1-2, pp. 100–107, 2005.
- J. Yuan, P. K. M. Cheung, H. M. Zhang, D. Chau, and D. Yang, “Inhibition of coxsackievirus B3 replication by small interfering RNAs requires perfect sequence match in the central region of the viral positive strand,” Journal of Virology, vol. 79, no. 4, pp. 2151–2159, 2005.
- J. Huang, F. Wang, E. Argyris et al., “Cellular microRNAs contribute to HIV-1 latency in resting primary CD4+ T lymphocytes,” Nature Medicine, vol. 13, no. 10, pp. 1241–1247, 2007.
- N. Arnberg, “Adenovirus receptors: implications for targeting of viral vectors,” Trends in Pharmacological Sciences, vol. 33, no. 8, pp. 442–448, 2012.
- Z. C. Hartman, E. P. Black, and A. Amalfitano, “Adenoviral infection induces a multi-faceted innate cellular immune response that is mediated by the toll-like receptor pathway in A549 cells,” Virology, vol. 358, no. 2, pp. 357–372, 2007.
- J. S. Yoon, H. H. Kim, Y. Lee, and J. S. Lee, “Cytokine induction by respiratory syncytial virus and adenovirus in bronchial epithelial cells,” Pediatric Pulmonology, vol. 42, no. 3, pp. 277–282, 2007.
- A. Heintz-Buschart, D. Yusuf, A. Kaysen et al., “Small RNA profiling of low biomass samples: identification and removal of contaminants,” BMC Biology, vol. 16, no. 1, p. 52, 2018.
- W. Ahmed and M. Razzaq, “RETRACTED: A small non-coding RNA AbcR2 regulate gntR transcription factor that modulate the intracellular survival of Brucella melitensis,” Microbial Pathogenesis, vol. 118, pp. 118–125, 2018.
- D. P. Bartel, “MicroRNAs: genomics, biogenesis, mechanism, and function,” Cell, vol. 116, no. 2, pp. 281–297, 2004.
- D. Piedade and J. M. Azevedo-Pereira, “MicroRNAs as important players in host-adenovirus interactions,” Frontiers in Microbiology, vol. 8, article 1324, 2017.
- Y. Qi, J. Tu, L. Cui et al., “High-throughput sequencing of microRNAs in adenovirus type 3 infected human laryngeal epithelial cells,” Journal of Biomedicine & Biotechnology, vol. 2010, Article ID 915980, 8 pages, 2010.
- H. Zhao, M. Chen, C. Tellgren-Roth, and U. Pettersson, “Fluctuating expression of microRNAs in adenovirus infected cells,” Virology, vol. 478, pp. 99–111, 2015.
- F. Zhong, N. Zhou, K. Wu et al., “A SnoRNA-derived piRNA interacts with human interleukin-4 pre-mRNA and induces its decay in nuclear exosomes,” Nucleic Acids Research, vol. 43, no. 21, pp. 10474–10491, 2015.
- J. R. Frost, M. Mendez, A. M. Soriano et al., “Adenovirus 5 E1A-mediated suppression of p 53 via FUBP1,” Journal of Virology, vol. 92, no. 14, 2018.
- J. Zhang, F. Huang, L. Tan et al., “Host protein Moloney leukemia virus 10 (MOV10) acts as a restriction factor of influenza a virus by inhibiting the nuclear import of the viral nucleoprotein,” Journal of Virology, vol. 90, no. 8, pp. 3966–3980, 2016.
- P. A. Tambyah, S. Sepramaniam, J. Mohamed Ali et al., “MicroRNAs in circulation are altered in response to influenza a virus infection in humans,” PLoS One, vol. 8, no. 10, article e76811, 2013.
- A. Onnis, M. Navari, G. Antonicelli et al., “Epstein-Barr nuclear antigen 1 induces expression of the cellular microRNA hsa-miR-127 and impairing B-cell differentiation in EBV-infected memory B cells. New insights into the pathogenesis of Burkitt lymphoma,” Blood Cancer Journal, vol. 2, no. 8, p. e84, 2012.
- Q.-W. Wang, Y. Su, J.-T. Sheng et al., “Anti-influenza A virus activity of rhein through regulating oxidative stress, TLR4, Akt, MAPK, and NF-κB signal pathways,” PLoS One, vol. 13, no. 1, article e0191793, 2018.
- W. Kang, L. Shang, T. Wang, H. Liu, and S. Ge, “Rho-kinase inhibitor Y-27632 downregulates LPS-induced IL-6 and IL-8 production via blocking p38 MAPK and NF-κB pathways in human gingival fibroblasts,” Journal of Periodontology, vol. 89, no. 7, pp. 883–893, 2018.
- G. James, B. Key, and A. Beverdam, “The E3 ubiquitin ligase Mycbp 2 genetically interacts with Robo 2 to modulate axon guidance in the mouse olfactory system,” Brain Structure & Function, vol. 219, no. 3, pp. 861–874, 2014.
- S. Holland and K. Scholich, “Regulation of neuronal functions by the E3-ubiquitinligase protein associated with MYC (MYCBP2),” Communicative & Integrative Biology, vol. 4, no. 5, pp. 513–515, 2014.
- Y. Zhao, M. J. C. Long, Y. Wang, S. Zhang, and Y. Aye, “Ube2V2 is a Rosetta stone bridging redox and ubiquitin codes, coordinating DNA damage responses,” ACS Central Science, vol. 4, no. 2, pp. 246–259, 2018.
- L. J. Simpson and J. E. Sale, “UBE2V2 (MMS2) is not required for effective immunoglobulin gene conversion or DNA damage tolerance in DT40,” DNA Repair (Amst), vol. 4, no. 4, pp. 503–510, 2005.
- C. Polge, N. Koulmann, A. Claustre et al., “UBE2D2 is not involved in MuRF1-dependent muscle wasting during hindlimb suspension,” The International Journal of Biochemistry & Cell Biology, vol. 79, pp. 488–493, 2016.
- S. Geisler, S. Vollmer, S. Golombek, and P. J. Kahle, “The ubiquitin-conjugating enzymes UBE2N, UBE2L3 and UBE2D2/3 are essential for Parkin-dependent mitophagy,” Journal of Cell Science, vol. 127, no. 15, pp. 3280–3293, 2014.
- M.-H. Chiang, L.-F. Chen, and H. Chen, “Ubiquitin-conjugating enzyme UBE2D2 is responsible for FBXW2 (F-box and WD repeat domain containing 2)-mediated human GCM1 (glial cell missing homolog 1) ubiquitination and degradation,” Biology of Reproduction, vol. 79, no. 5, pp. 914–920, 2008.
Copyright © 2018 Feng Huang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.