Micro-RNAs in regenerating lungs: an integrative systems biology analysis of murine influenza pneumonia
BMC Genomics volume 15, Article number: 587 (2014)
Tissue regeneration in the lungs is gaining increasing interest as a potential influenza management strategy. In this study, we explored the role of microRNAs, short non-coding RNAs involved in post-transcriptional regulation, during pulmonary regeneration after influenza infection.
We profiled miRNA and mRNA expression levels following lung injury and tissue regeneration using a murine influenza pneumonia model. BALB/c mice were infected with a sub-lethal dose of influenza A/PR/8(H1N1) virus, and their lungs were harvested at 7 and 15 days post-infection to evaluate the expression of ~300 miRNAs along with ~36,000 genes using microarrays. A global network was constructed between differentially expressed miRNAs and their potential target genes with particular focus on the pulmonary repair and regeneration processes to elucidate the regulatory role of miRNAs in the lung repair pathways. The miRNA arrays revealed a global down-regulation of miRNAs. TargetScan analyses also revealed specific miRNAs highly involved in targeting relevant gene functions in repair such as miR-290 and miR-505 at 7 dpi; and let-7, miR-21 and miR-30 at 15 dpi.
The significantly differentially regulated miRNAs are implicated in the activation or suppression of cellular proliferation and stem cell maintenance, which are required during the repair of the damaged lungs. These findings provide opportunities in the development of novel repair strategies in influenza-induced pulmonary injury.
Vulnerability of lungs to influenza virus infection
Our lungs are constantly in contact with the external environment and exposed to airborne debris, pathogens and chemicals, which can inflict damage acutely or chronically . Pathogens such as influenza virus can cause acute lung injury (ALI) through acute viral pneumonia . ALI induced by influenza may also lead to complications that severely impair pulmonary functions, and not only affect an individual’s quality of life, but can also be life-threatening if left untreated [2–5]. Furthermore, with the emergence of pandemic H1N1-2009, avian H5N1 and H7N9 influenza strains, and seasonal H1N1 and H3N2 strains, which cause varying severity of ALI, there is a need for novel and timely intervention strategies for post-influenza lung repair to mitigate permanent tissue damage and its sequelae. Among available options, promoting pulmonary tissue regeneration represents an effective approach against influenza-induced lung injury [6, 7].
Regenerative medicine as a novel influenza management strategy
Influenza infection of the lungs constitutes an important source of severe inflammatory damage to the lung architecture via respiratory burst of the innate immune response that cannot be generalized to other types of lung injury [5, 8]. The injury can give rise to complications and/or chronic damage if not managed properly, particularly in susceptible groups such as the elderly and allergic individuals with diminished repair functions to cope with these infections [3, 4]. Being a common respiratory pathogen, influenza has therefore been extensively investigated to elucidate its infection kinetics and pathogenicity by clinical studies and animal models [9, 10]. Hence, with better understanding of its kinetics, novel strategies can be formulated based on the phases of the influenza infection process. The repair stage of influenza pneumonia is gaining increasing interest given that complications can be prevented if the damaged lung architecture is restored appropriately. Studies on the recovery period of pulmonary influenza infection have indicated changes in proliferative and differentiation traits of lung epithelial cells [11, 12], suggesting potential regulation of these changes to facilitate the repair of damaged lungs.
MicroRNAs as potential agents in post-influenza lung repair strategies
MicroRNAs (miRNAs) are short RNA nucleotides involved in post-transcriptional regulation, and are of major interest in view of their ease of manipulation and delivery into host cells due to their relatively small sizes [13–15]. In addition, they regulate various disease processes as well as normal physiologic functions, including embryonic lung development [13–20]. The recent success of miRNA-based drug discovery  further fuelled our interest to screen for miRNAs that mediate repair of influenza-afflicted lungs. We hypothesized that miRNAs play major roles in priming pulmonary tissues for repair and regeneration following influenza pneumonia, similar to how miRNAs are intimately involved in embryonic lung development [16, 17].
This study aimed to investigate miRNAs that are differentially regulated during the repair phases of influenza using an established murine model of influenza H1N1 pneumonia, drawing parallels to seasonal H1N1 infection that may culminate in pulmonary complications in susceptible individuals. Another objective was to attain complete systems-level analysis to elucidate potential miRNAs that are crucial for lung repair in vivo. These findings can improve our understanding of influenza-induced lung repair responses, and also identify miRNA candidates for further validation as critical players in pulmonary regeneration.
Infected murine lungs exhibit similar degrees of damage, while miRNA and mRNA microarray expression profiles are largely homogenous at specific time-points
The lungs on 7 days post-infection (dpi) were scored using a modified damage scoring system  to ascertain that the mice suffered similar extents of damage, while lungs on 15 dpi were scored based on the area of recovery. Lungs that displayed similar degrees of injury within each individual time-point at 7 and 15 dpi were subjected to miRNA and mRNA microarrays (Table 1). Pearson correlation was computed between all samples to compare the homogeneity of miRNA and mRNA expression within each animal group, which showed that miRNA and mRNA expression profiles were generally highly correlated. For subsequent analyses, only the probe-set intensity data that were relatively homogenous between replicates of the specific infected groups at 7 and 15 dpi were included. Thus, consistent miRNA profiles from 3 mice at 7 dpi and 4 mice at 15 dpi, together with consistent mRNA profiles from 4 mice each at 7 and 15 dpi were analyzed (Additional file 1: Figure S1).
Cell types, degree of cell proliferation and repair in the lungs vastly differ between 7 and 15 dpi
Compared to uninfected control, histopathologic analyses of infected lungs (Figure 1a) revealed more damaged areas, and a greater number of infiltrating polymorphs at 7 dpi, suggesting predominant innate immune responses against the virus in the lungs. This finding is congruent with the viral NS1 expression, which was still detected in lungs at 7 dpi (Figure 1b). Most areas of the lungs exhibited widening interstitia and septae, with the presence of fibrin and red blood cells. In contrast, at 15 dpi, while some features of broncho-alveolar damage remained (i.e. apoptotic debris, desquamation), viral expression was absent in the lungs. Furthermore, the most striking change was the appearance of new epithelial cells at 15 dpi. Also, the subjective increase in all types of infiltrating cells on 15 dpi is likely to account for the apparent increase in area of lung damage compared with 7 dpi.Immunohistochemistry (IHC) analyses with proliferating cell nuclear antigen (PCNA) and surfactant protein-C (SP-C) staining were also conducted to ascertain the degree of cellular proliferation, DNA synthesis and repair in the lungs (specifically for alveolar type 2 or AT2 pneumocytes). Minimal PCNA-positive cells were present in the control mice at both days, while SP-C-positive AT2 pneumocytes were evenly distributed across the lungs (Figure 1c). Lungs of infected mice at 7 dpi displayed areas without SP-C-positive cells likely attributed to virus-induced damage. Increased PCNA-positive cells were observed in the infected lungs. However, there was only a small percentage of PCNA and SP-C double-positive cells at 7 dpi, which was not significant compared to control uninfected lungs (Figure 1d, e). On 15 dpi, multiple areas deprived of SP-C-positive cells were likely attributed to pneumocytes damaged by infection and innate immune responses. Relatively more PCNA-positive cells as well as PCNA and SP-C double-positive cells were observed in the lungs at 15 dpi, indicating enhanced cell proliferation, DNA synthesis and DNA repair of the pneumocytes.
General down-regulation of miRNAs and up-regulation of genes at repair phases on both 7 and 15 dpi
The miRNA array data on both 7 and 15 dpi revealed a global down-regulation of miRNAs (at false discovery rate or FDR of 5%). On 7 dpi, 3 miRNAs were up-regulated and 22 miRNAs were down-regulated out of 336 miRNAs analyzed. However, on 15 dpi, 12 miRNAs were up-regulated and 114 down-regulated out of 314 miRNAs analyzed (Figure 2a). Table 2 outlines some of the top significant miRNAs on both days. This landscape suggests that miRNAs, being gene regulatory factors via translational repression, were down-regulated to yield augmented gene expression in response to H1N1 infection. Indeed, the mRNA expression patterns were skewed towards up-regulation in the infected groups, with 591 up-regulated genes and 246 down-regulated genes at 7 dpi (out of 36,991 genes analyzed); and 1,768 up-regulated genes and 744 down-regulated genes at 15 dpi (out of 35,577 genes analyzed). In comparison, infected mice revealed ~5 times more differentially expressed (DE) miRNAs (at least 1.3× fold change) and ~3 times more DE genes (at least 1.5× fold change) on 15 dpi compared to 7 dpi. This indicates more significant alterations in expression at 15 dpi, presumably due to complete viral clearance. There was almost no overlap between the DE miRNAs and genes between the two time-points, implying that the biological processes were largely distinct due to the different conditions of the infected lungs at each time-point (Figure 2b).
Real-time RT-PCR data of representative significantly DE miRNAs and genes indicated good consistency with the microarray expression data (Figure 3). Thus, both miR-335 and miR-582 were down-regulated at 7 dpi; while miR-21 was up-regulated, miR-34b and miR-542 were down-regulated at 15 dpi. For gene expression at 7 dpi, genes such as ANGPTL4 and PLAT were up-regulated, while genes such as SOX4 and CTHRC1 were down-regulated. At 15 dpi, genes such as EREG and GDF6 were up-regulated, while genes such as ESM1 and FGF1 were down-regulated. The directional trends of changes in all the DE genes were consistent with the array data, thus validating the accuracy of both miRNA and gene microarrays.
Overall global gene networks by Ingenuity pathway analysis (IPA)
Table 3 shows the top 5 IPA networks of the DE genes for each time-point. At 7 dpi, functions and pathways such as immune cell proliferation, inflammation, and innate viral responses were significantly enriched (P < 0.05), while repair functions were relatively sparse. The number of functions and pathways was much lower at 7 dpi likely due to viral presence and immediate innate immune responses to combat the virus, thus masking other functions. On the other hand, following viral clearance by 15 dpi, innate immune responses were likely to be largely replaced by more influenza-specific adaptive responses that are more tightly controlled and regulated. Thus, the pathway enrichment shifted from immune responses to largely cellular responses of repair and proliferation, such as DNA damage response, cellular proliferation, DNA replication and cell cycle activities (P < 0.05).
Repair pathway selection and gene ontology (GO) analyses
To focus on miRNA regulatory activities in the relevant biological functions, we applied the IPA function tree to manually select genes that were annotated in repair, development and regeneration functions. Thus, 227 and 876 genes were selected for further analyses for days 7 and 15 (out of a total of 837 and 2,512), respectively. In this selection, the DE genes were mainly involved in cell death, cell survival, cell cycle, cell proliferation, DNA replication, and developmental functions.The GO term analyses applied after the initial manual selection revealed different processes occurring at each time-point (Figure 4). At 7 dpi, only a relatively small number of 53 GO terms were identified to be directly involved in repair, mainly in tissue development, cell cycle, cell death, and wound healing, implying minimal repair and developmental functions during active viral clearance. In contrast, 254 GO terms were found to be associated directly with repair and development at 15 dpi. Many of them were predominantly involved mainly in cell cycle, cell death, tissue development and growth factors, which corroborates the repair and replacement of damaged lung tissue occurring at 15 dpi.
TargetScan reveals the close relationship between miRNAs and selected repair-related genes
To evaluate the extent to which miRNAs were involved with gene regulation in repair and regeneration, TargetScan analyses of the DE miRNAs targeting the selected genes related to repair and regeneration were applied after the GO analyses. Thus, the DE genes were further narrowed down to 61 up- and 43 down-regulated genes at 7 dpi, targeted by the corresponding 2 up- and 10 down-regulated miRNAs. At 15 dpi, there were 556 up- and 193 down-regulated repair-related genes targeted by the corresponding 12 up- and 91 down-regulated miRNAs. When the DE miRNA and repair-associated genes were paired by TargetScan, we thus observed 12 miRNAs at 7 dpi and 103 miRNAs at 15 dpi highly involved in targeting repair-associated genes that were differentially regulated at both time-points. Figure 5 shows the major targets of the altered miRNAs grouped into their respective GO terms. The high frequency of repair GO terms targeting genes in cell cycle, differentiation, wound healing, tissue development and angiogenesis, implied that the repair and regeneration of lung tissues following influenza pneumonia were tightly modulated by miRNAs.
Selection of candidates having high association with lung repair functions
To screen for the final list of potential candidate miRNA regulators of repair and regeneration, we applied other criteria (Figure 6a). The rationale of this additional filter is that miRNAs changes are usually subtle, and therefore selecting those significantly DE miRNAs with the highest or lowest raw abundance in the control samples would likely culminate in a larger magnitude of change and effect in response to infection. These additional criteria narrowed the selection down to 3 miRNAs at 7 dpi (miR-290, miR-1940, miR-505) which were marked to be crucial in the early repair phase. At 15 dpi, 17 miRNAs (let-7b,c, miR-10a, miR-21, miR-25, miR-26a, miR-29c, miR-30a,b,c,d, miR-99a, miR-103, miR-151, miR-195, and miR-200b,c) were identified to be important in the late repair phase. Most of these miRNAs were also found to target genes in lung developmental GO terms, rendering them the most likely regulators of the transcriptome in influenza-induced pulmonary repair (Figure 6b).
Possible functions of high-scoring miRNAs involved in targeting the lung repair pathway
Of the 20 shortlisted miRNAs, their putative functions associated with lung repair were inferred through literature searches. (Additional file 2: Table S4) indicates that the DE miRNAs on days 7 and 15 largely differ in their functions. The miRNAs at 7 dpi were engaged in endothelial and fibroblast proliferation (wound healing), whereas miRNAs at 15 dpi were more broadly involved in cell proliferation, apoptosis, tumor suppression and DNA repair, consistent with the functional annotation. The full list of targets of these shortlisted miRNAs is provided in Additional file 3: Table S1.
miRNAs are believed to play crucial roles during the repair of the lungs following injury. We therefore identified miRNAs that are differentially regulated in a time-dependent manner during the early and later recovery phases of the lungs using a murine influenza infection model. The recovery time-points were selected based on previously documented kinetics of lung-repairing cell types initially present on day 7 (early phase) and that peak on day 15 (late phase) [10, 12]. Our in vivo investigation unveiled potential candidate miRNAs strongly associated with pulmonary repair and regeneration after acute inflammatory lung injury due to influenza pneumonia.
Marked expression of repair, proliferation and differentiation genes after viral clearance at 15 dpi
At 7 dpi, we observed highly specific differential expression of mainly inflammatory and viral innate immune response-related genes in the lungs. This likely reflects the specific commitment of damaged lungs to combat viral infection rather than for repair functions. In addition, the respiratory burst of innate immune responses also causes oxidative stress and DNA damage, culminating in greater genetic instability and apoptosis which hamper cellular proliferation for replacing damaged cells [24–26]. However, day 15 corresponds to the later phase of recovery, characterized by viral clearance e.g. by lymphocytes in the lungs to execute more specific adaptive responses. The emergence of new immature epithelial cells as indicated by histologic analyses as well as increased PCNA-positive cells also suggest re-epithelialization of lung tissues following viral clearance. This was evident from the patches of PCNA-positive cells within areas devoid of SP-C-positive cells, which may reflect epithelial progenitor cells or lung stem cells prior to differentiating into new pneumocytes to replace the damaged ones. Regeneration was also illustrated by the significantly increased double-positive cells lining the edges of these SP-C-deprived areas. Therefore, more efficient repair processes would be expected at 15 dpi, with less oxidative stress when the respiratory burst subsided in the lungs. Alveolar epithelial cells could thus be replaced more readily since repair would less likely be compromised by DNA damage and mutation .
Potential regulatory roles of specific miRNAs in lung repair and development
The number of DE miRNAs was directly proportional to the number of DE genes involved in repair and development. Most DE miRNAs directly target repair-associated genes, especially at 15 dpi when extensive repair occurred. The close association between miRNAs and lung development has been documented [13, 19, 27–29]. The differentially regulated miRNAs were unique with almost no overlap between the two time-points, indicating that miRNA expression changed transiently according to the temporal gene expression requirements in the lungs. The top 20 miRNA candidates were identified, together with their roles and targets during the lung repair phases following influenza infection. Some of these miRNAs are highly involved in development, regeneration, and cancer progression due to their apparent roles in proliferation and tumor suppression. This underscores that tissue repair entails essentially tightly controlled proliferation as opposed to uncontrolled proliferation in cancer which involves similar factors .
We therefore postulate that the increase of certain stem cell-associated miRNAs (e.g. miR-290) promotes the stemness of cells during influenza-induced lung repair [30–32]. The let-7 miRNAs are implicated in developmental functions and cellular proliferation, and also regulate and inhibit proliferation in a cancer context. Therefore, their transiently diminished expression during repair may facilitate the proliferation of new tissues [33–35]. The miR-200 family is often enriched in epithelial tissues, and their decreased expression is linked to epithelial-mesenchymal transition and cell stemness. Hence, their reduced expression may suggest enhanced plasticity of epithelial cells in repairing lungs to replace damaged areas . Lowered miR-200 expression may also augment endothelial development to repair the lung vasculature . Interestingly, Li et al  found that miR-200a is differentially expressed in influenza virus infection, and that miRNA expression profiles vary depending on the virus strains. miR-21 targets proliferation-suppressing factors, and its elevated expression during repair coincides with increased proliferation in repairing lungs . Notably, increased miR-21 may cause detrimental effects due to their targets in the TGF-β pathway, suggesting that miR-21 knockdown may prevent lung complications such as fibrosis . On the other hand, some other miRNAs are possibly involved in the regulation and surveillance of repair to prevent uncontrolled cell proliferation. Thus, members of the miR-30 family were significantly down-regulated so that expression of its main target p53 could be suitably elevated to counteract the higher proliferation in recovering lung tissues, which are more prone to DNA damage and mutation in the presence of increased DNA synthesis . Hence, miR-30 appears to act as a tumor suppressor, with its subdued expression facilitating proliferation, but concurrently activating the negative feedback loop of p53, thus showcasing the intricate roles that miRNAs play in pulmonary damage and repair .
The functions of the physiologically important miRNAs during tissue repair may reflect an intermediate state between organismal development and cancer progression. Transient and controlled changes in expression of proliferation activators and inhibitors may lead to a state akin to development in contrast to cancer progression, as highlighted by the vastly different DE miRNA profiles between days 7 and 15.
miRNAs and their potential applications in lung tissue damage and repair
This study revealed that the regulation of miRNAs was unique during different phases of repair after influenza infection. The regulatory roles of these miRNAs in tissue repair can be tested for their veracity in improving lung repair and the overall outcome of influenza pneumonia. For example, due to the nature of miRNA targeting, some of these miRNAs (such as miR-21 and let-7) may even serve dual roles of limiting damage and accelerating repair as they possess anti-inflammatory properties besides regulating pulmonary repair activities, i.e. they may minimize damage from inflammation secondary to infection [43, 44]. Therefore, future studies are warranted to further characterize the potential miRNA candidates in vitro and in vivo, in order to determine their detailed roles in ALI and repair following influenza pneumonia.
We anticipate that many of the biological processes observed during influenza lung damage and repair may also be relevant to other types of viral infection and lung injury, and will therefore contribute to our basic understanding of tissue repair of damaged lungs. Ultimately, this integrative systems biology approach highlights the significant roles of the miRNAs in acute inflammatory lung injury and regeneration, and revealed certain miRNAs that may be instrumental in the repair and reconstruction of the lung architecture post-infection.
Virus strain and infection of animals
All animal experiments with 6–8 week old female BALB/c mice were approved by the Institutional Animal Care and Use Committee (IACUC), National University of Singapore (protocol 050/11), and fulfilled all listed criteria in the ARRIVE guidelines (Additional file 4). Mice were housed in ABSL-2 facilities in ventilated cages according to IACUC guidelines. Animals were divided into 4 groups each consisting of 4 mice, i.e. control uninfected and infected groups at 7 and 15 dpi. The two time-points were chosen based on the kinetics of influenza infection, i.e. when repair of primary viral damage commences and certain repair cell types initially appear at 7 dpi, whereas repair of inflammatory damage begins after viral clearance and when repair cell types peak at 15 dpi [10, 12]. Mice were anesthetized with a mixture of 7.5 mg/ml ketamine and 0.1 mg/ml medetomidine, and infected intra-tracheally with a sub-lethal dose of 23 plaque-forming units of influenza A/Puerto Rico/8/1934(H1N1) or PR8 virus (in 50 μl volume) to ensure that all mice recovered from the infection. Control mice were given phosphate-buffered saline (PBS) intra-tracheally. Mice were monitored daily for weight loss from infection and weight gain from recovery. Upon reaching the stipulated time-point, mice were euthanized in a CO2 chamber, and the lungs harvested for histology and RNA extraction.
Murine lung tissue processing and extraction of miRNAs and mRNAs
Approximately 10% of each lobe of lung was snap-frozen at -80°C to preserve RNAs. The lungs were then homogenized in Trizol solution (Qiagen, Hilden, Germany) using a gentleMACS tissue dissociator (Miltenyi Biotec, Bergisch Gladbach, Germany), and total RNA was immediately extracted using the miRNeasy miRNA extraction kit (Qiagen). The isolated total RNA samples containing miRNAs were then subjected to miRNA and gene microarray analyses.
The remaining 90% of each lung was fixed in 4% formalin, dehydrated in ascending ethanol concentrations, embedded in paraffin, sectioned into 4-μm slices, and stained with hematoxylin and eosin (H&E). Ideally, the lungs could be inflated prior to fixation. Notwithstanding this, highly satisfactory assessment of the infected areas was achieved by an experienced pulmonary pathologist to include: inflammatory cell infiltration (polymorphs and lymphocytes); septae, hyaline and fibrin; hemorrhage; edema (for tissue damage); and new progenitor cells and re-epithelization (for tissue repair and replacement). This evaluation was then correlated with the molecular analyses. Scoring of damage to ascertain the homogeneity of the infection was also conducted for mouse lungs at 7 dpi based on a modified scoring system, i.e.% infected × [alveolar hemorrhage + 2(alveolar inflitrates) + 3(fibrin) + alveolar septal congestion] . The scoring was carried out based on the worst-affected area of the lungs, which then yielded the percentage of the lung area showing similar scores. Since mice at 15 dpi were recovering and could not be scored with the same system, the percentage of lung area affected by the infection was obtained instead to ascertain the homogeneity of damage.
Paraffin-embedded blocks of lung tissues were also sectioned into 5-μm slices for IHC analyses. The slides were dewaxed with two washes of xylene, rehydrated with two washes of 100% ethanol, and one wash each of 90%, 70% and 50% ethanol, followed by two washes of deionized water. Slides were then subjected to antigen retrieval in boiling sodium citrate buffer with Tween 20. Slides were further permeated with 0.025% of Triton X-100 for 10 min, before being blocked with 3% BSA for 2 h. Primary antibodies for PCNA (SC-9857, Santa Cruz Biotechnology, Dallas, TX) and SP-C (SC-13979, Santa Cruz Biotechnology) at 100× dilution in 3% BSA were used for co-staining dividing and repairing cells and AT2 pneumocytes, respectively. Following overnight incubation at 4°C in the dark, the slides were washed with Tris-buffered saline thrice for 5 min each, and stained with secondary antibody containing Alexa Fluor 488 for PCNA and Alexa Fluor 568 for SP-C (Molecular Probes, Grand Island, NY) at room temperature for 1 h in the dark. Slides were then washed and mounted in colloidal gold Antifade mounting medium with DAPI (Molecular Probes). Slides were then concurrently scanned with the high-resolution MIRAX MIDI system equipped with fluorescence illumination (Carl Zeiss, Jena, Germany) for the complete lung image to quantify dividing AT2 pneumocytes. Stained sections were also visualized with a FX-1000 confocal microscope at higher magnification.
Image analyses for quantification of dividing AT2 pneumocytes in the lungs
An automated computer algorithm was developed to count the number of total cells, PCNA-positive cells, and PCNA-positive AT2 pneumocytes. Nuclei were segmented by converting the DAPI fluorescence channel into black-and-white using manual thresholding of DAPI intensities, and removing areas less than 3 pixels in size. The number of nuclei was counted as the number of total cells. The PCNA fluorescence channel was next converted to black-and-white using manual thresholding of PCNA intensities. PCNA-positive nuclei were identified by selecting co-localized areas of nuclei and PCNA channel, and performing dilating-hole filling-eroding using “disk” operator with 3 pixels in size. The number of PCNA-positive nuclei was counted as the number of PCNA-positive cells. Finally, the SP-C fluorescence channel was converted into black-and-white using manual thresholding of SP-C intensities, and areas less than 2 pixels in size were ignored. The surrounding areas of PCNA-positive nuclei (1 pixel in radius) were then checked for SP-C signal in the SP-C channel. The number of PCNA-positive nuclei with surrounding SP-C signal was counted as the number of PCNA-positive AT2 cells. All image processing and computation algorithms were implemented using MATLAB with image processing toolbox (MathWorks, Natick, MA). The MATLAB codes are available upon request.
miRNA and gene microarray analyses
Extracted total RNAs containing miRNAs were subjected to quality control using a bio-analyzer to ensure proper RNA quality before the microarray experiments. Total RNA samples from 4 mice per group were then hybridized individually onto the respective microarray platforms (Agilent, Santa Clara, CA). The miRNA microarray was a Mouse miRNA, 8 × 60 K format platform (AMADID 38112), while the gene microarray was a SurePrint G3 Mouse GE, 8 × 60 K, 1 color format platform (AMADID 028005). A total of 32 arrays (i.e. 16 miRNA arrays and 16 mRNA arrays) were hybridized to total RNA samples extracted from the lungs of 4 infected mice and 4 control mice for each 7 and 15 dpi. The expression values of the array probes were calculated following the manufacturer’s standard array protocol. Only probes above the expression threshold of 20 (raw signal intensity) were used to derive expression values.
The expression data for the filtered probes of both miRNA (1,179 probes) and gene (55,681 probes) microarrays were normalized using quantile normalization, and were subjected to DE analyses using QPROT software (an extension of the QSPEC software for generic expression data with missing values). DE genes were identified using the Z-statistics threshold that controls the FDR at 5% . The DE miRNAs and genes were further analyzed below to shortlist candidates associated with active roles in pulmonary repair and regeneration.
Real-time RT-PCR validation of miRNA and gene expression
A total of 5 miRNAs and 20 genes that were significantly regulated at both 7 and 15 dpi were selected for real-time RT-PCR to validate their differential expression based on the microarray data. A stem-loop primer miRNA Universal TaqMan RT-PCR system (Applied Biosystems, Foster City, CA) was utilized to validate miRNAs, and normalized using murine SNORD68 small nucleolar RNA [46, 47]. The miRNA stem-loop primers, specific forward primers and universal reverse primer are listed in Additional file 5: Table S2. Real-time RT-PCR assays to validate 20 genes were performed using reverse transcription with random hexamers (Promega, Madison, WI), the SYBR Green Master mix and LightCycler real-time PCR system (Roche, Penzberg, Germany), normalized against the RPL13a gene (Additional file 6: Table S3).
Developmental pathway and TargetScan analyses
The DE genes identified above were filtered for repair-associated genes by the following criteria. Functional selection for repair association was performed manually using the function tree of the IPA software (Qiagen), where genes under developmental and repair functional groups were identified. These genes were then filtered again with the hypergeometric test and FDR control (the Benjamini-Hochberg procedure), and GO analysis . A concurrent miRNA TargetScan analysis at a PCT score of 0.5 was performed using the GeneSpring GX software (Agilent), and only miRNAs that target the filtered genes from the previous steps were selected. The selected genes (grouped as GO terms at P < 0.05) and miRNAs were then visualized in a network with Cytoscape software (v3.0.1).
Further refinement of miRNA candidate selection
To further narrow down the selected miRNA candidates, an additional filter system of miRNA raw abundance in control samples was taken into account. Therefore, a final set of miRNAs was selected based on their statistical significance, targeting related GO term members of gene analysis, highest or lowest raw abundance in control samples where the fold changes were deemed more significant. The selected miRNAs were then classified based on their expression changes and their implicated functions, and also whether they target genes directly associated with lung development and regeneration.
FDR calculation and GO enrichment analyses of the array dataset were performed using R statistics software for selection of significant DE miRNAs and genes. Student’s one sample t-test was performed on real-time RT-PCR verification and dividing AT2 pneumocyte quantification data with GraphPad and SPSS software. Results were represented as mean ± standard deviation (SD).
Acute lung injury
Proliferating cell nuclear antigen
Alveolar type 2
False discovery rate
Ingenuity pathway analysis
Institutional animal care and use committee
Hematoxylin and eosin
van der Vliet A, Cross CE: Oxidants, nitrosants, and the lung. Am J Med. 2000, 109: 398-421.
Hendrickson CM, Matthay MA: Viral pathogens and acute lung injury: investigations inspired by the SARS epidemic and the 2009 H1N1 influenza pandemic. Semin Respir Crit Care Med. 2013, 34: 475-486.
Rothberg MB, Haessler SD, Brown RB: Complications of viral influenza. Am J Med. 2008, 121: 258-264.
Rothberg MB, Haessler SD: Complications of seasonal and pandemic influenza. Crit Care Med. 2010, 38: e91-97.
Narasaraju T, Yang E, Samy RP, Ng HH, Poh WP, Liew AA, Phoon MC, van Rooijen N, Chow VT: Excessive neutrophils and neutrophil extracellular traps contribute to acute lung injury of influenza pneumonitis. Am J Pathol. 2011, 179: 199-210.
Warburton D, El-Hashash A, Carraro G, Tiozzo C, Sala F, Rogers O, Langhe SD, Kemp PJ, Riccardi D, Torday J, Bellusci S, Shi W, Lubkin SR, Jesudason E: Lung organogenesis. Curr Top Dev Biol. 2010, 90: 73-158.
Warburton D, Olver BE: Coordination of genetic, epigenetic, and environmental factors in lung development, injury, and repair. Chest. 1997, 111: 119S-122S.
White MR, Crouch E, Vesona J, Tacken PJ, Batenburg JJ, Leth-Larsen R, Holmskov U, Hartshorn KL: Respiratory innate immune proteins differentially modulate the neutrophil respiratory burst response to influenza A virus. Am J Physiol Lung Cell Mol Physiol. 2005, 289: L606-L616.
Smith AM, Perelson AS: Influenza A virus infection kinetics: quantitative data and models. Wiley Interdiscip Rev Syst Biol Med. 2011, 3: 429-445.
Yin L, Xu S, Cheng J, Zheng D, Limmon GV, Leung NH, Rajapakse JC, Chow VT, Chen J, Yu H: Spatiotemporal quantification of cell dynamics in the lung following influenza virus infection. J Biomed Opt. 2013, 18: 046001-
Zheng D, Limmon GV, Yin L, Leung NH, Yu H, Chow VT, Chen J: Regeneration of alveolar type I and II cells from Scgb1a1-expressing cells following severe pulmonary damage induced by bleomycin and influenza. PLoS One. 2012, 7: e48451-
Zheng D, Limmon GV, Yin L, Leung NH, Yu H, Chow VT, Chen J: A cellular pathway involved in Clara cell to alveolar type II cell differentiation after severe lung injury. PLoS One. 2013, 8: e71028-
Khoshgoo N, Kholdebarin R, Iwasiow BM, Keijzer R: MicroRNAs and lung development. Pediatr Pulmonol. 2013, 48: 317-323.
Sen CK: MicroRNAs as new maestro conducting the expanding symphony orchestra of regenerative and reparative medicine. Physiol Genomics. 2011, 43: 517-520.
Jackson AL, Levin AA: Developing microRNA therapeutics: approaching the unique complexities. Nucleic Acid Ther. 2012, 22: 213-225.
Lu Y, Okubo T, Rawlins E, Hogan BL: Epithelial progenitor cells of the embryonic lung and the role of microRNAs in their proliferation. Proc Am Thorac Soc. 2008, 5: 300-304.
Lu Y, Thomson JM, Wong HY, Hammond SM, Hogan BL: Transgenic over-expression of the microRNA miR-17-92 cluster promotes proliferation and inhibits differentiation of lung epithelial progenitor cells. Dev Biol. 2007, 310: 442-453.
Dakhlallah D, Batte K, Wang Y, Cantemir-Stone CZ, Yan P, Nuovo G, Mikhail A, Hitchcock CL, Wright VP, Nana-Sinkam SP, Piper MG, Marsh CB: Epigenetic regulation of miR-17 ~ 92 contributes to the pathogenesis of pulmonary fibrosis. Am J Respir Crit Care Med. 2013, 187: 397-405.
Sayed D, Abdellatif M: MicroRNAs in development and disease. Physiol Rev. 2011, 91: 827-887.
Li Y, Kowdley KV: MicroRNAs in common human diseases. Genomics Proteomics Bioinformatics. 2012, 10: 246-253.
Hu J, Xu Y, Hao J, Wang S, Li C, Meng S: MiR-122 in hepatic function and liver diseases. Protein Cell. 2012, 3: 364-371.
Zhou T, Garcia JG, Zhang W: Integrating microRNAs into a system biology approach to acute lung injury. Transl Res. 2011, 157: 180-190.
Matute-Bello G, Winn RK, Jonas M, Chi EY, Martin TR, Liles WC: Fas (CD95) induces alveolar epithelial cell apoptosis in vivo: implications for acute pulmonary inflammation. Am J Pathol. 2001, 158: 153-161.
Reliene R, Fischer E, Schiestl RH: Effect of N-acetyl cysteine on oxidative DNA damage and the frequency of DNA deletions in atm-deficient mice. Cancer Res. 2004, 64: 5148-5153.
Bishop AJ, Kosaras B, Carls N, Sidman RL, Schiestl RH: Susceptibility of proliferating cells to benzo[a]pyrene-induced homologous recombination in mice. Carcinogenesis. 2001, 22: 641-649.
Aguilera A, Garcia-Muse T: Causes of genome instability. Annu Rev Genet. 2013, 47: 1-32.
Dong J, Jiang G, Asmann YW, Tomaszek S, Jen J, Kislinger T, Wigle DA: MicroRNA networks in mouse lung organogenesis. PLoS One. 2010, 5: e10854-
Mujahid S, Logvinenko T, Volpe MV, Nielsen HC: miRNA regulated pathways in late stage murine lung development. BMC Dev Biol. 2013, 13: 13-
Tomankova T, Petrek M, Kriegova E: Involvement of microRNAs in physiological and pathological processes in the lung. Respir Res. 2010, 11: 159-
Houbaviy HB, Murray MF, Sharp PA: Embryonic stem cell-specific MicroRNAs. Dev Cell. 2003, 5: 351-358.
Kaspi H, Chapnik E, Levy M, Beck G, Hornstein E, Soen Y: miR-290-295 regulate embryonic stem cell differentiation propensities by repressing Pax6. Stem Cells. 2013, 31: 2266-2272.
Kuppusamy KT, Sperber H, Ruohola-Baker H: MicroRNA regulation and role in stem cell maintenance, cardiac differentiation and hypertrophy. Curr Mol Med. 2013, 13: 757-764.
Cimadamore F, Amador-Arjona A, Chen C, Huang CT, Terskikh AV: SOX2-LIN28/let-7 pathway regulates proliferation and neurogenesis in neural precursors. Proc Natl Acad Sci U S A. 2013, 110: E3017-3026.
Barh D, Malhotra R, Ravi B, Sindhurani P: MicroRNA let-7: an emerging next-generation cancer therapeutic. Curr Oncol. 2010, 17: 70-80.
Zhao B, Han H, Chen J, Zhang Z, Li S, Fang F, Zheng Q, Ma Y, Zhang J, Wu N, Yang Y: MicroRNA let-7c inhibits migration and invasion of human non-small cell lung cancer by targeting ITGB3 and MAP4K3. Cancer Lett. 2014, 342: 43-51.
Creighton CJ, Gibbons DL, Kurie JM: The role of epithelial-mesenchymal transition programming in invasion and metastasis: a clinical perspective. Cancer Manag Res. 2013, 5: 187-195.
Chan YC, Roy S, Khanna S, Sen CK: Downregulation of endothelial microRNA-200b supports cutaneous wound angiogenesis by desilencing GATA binding protein 2 and vascular endothelial growth factor receptor 2. Arterioscler Thromb Vasc Biol. 2012, 32: 1372-1382.
Li Y, Chan EY, Li J, Ni C, Peng X, Rozsenzweig E, Tumpey TM, Katze MG: MicroRNA expression and virulence in pandemic influenza virus-infected mice. J Virol. 2010, 84: 3023-3032.
Zhong Z, Dong Z, Yang L, Gong Z: miR-21 induces cell cycle at S phase and modulates cell proliferation by down-regulating hMSH2 in lung cancer. J Cancer Res Clin Oncol. 2012, 138: 1781-1788.
Vaporidi K, Vergadi E, Kaniaris E, Hatziapostolou M, Lagoudaki E, Georgopoulos D, Zapol WM, Bloch KD, Iliopoulos D: Pulmonary microRNA profiling in a mouse model of ventilator-induced lung injury. Am J Physiol Lung Cell Mol Physiol. 2012, 303: L199-207.
Li J, Donath S, Li Y, Qin D, Prabhakar BS, Li P: miR-30 regulates mitochondrial fission through targeting p53 and the dynamin-related protein-1 pathway. PLoS Genet. 2010, 6: e1000795-
Kao CJ, Martiniez A, Shi XB, Yang J, Evans CP, Dobi A, Devere White RW, Kung HJ: miR-30 as a tumor suppressor connects EGF/Src signal to ERG and EMT. Oncogene. 2014, 33: 2495-2503.
Xue X, Xia W, Wenzhong H: A modeled dynamic regulatory network of NF-kappaB and IL-6 mediated by miRNA. Biosystems. 2013, 114: 214-218.
Kumar M, Ahmad T, Sharma A, Mabalirajan U, Kulshreshtha A, Agrawal A, Ghosh B: Let-7 microRNA-mediated regulation of IL-13 and allergic airway inflammation. J Allergy Clin Immunol. 2011, 128: 1077-1085.
Choi H, Fermin D, Nesvizhskii AI: Significance analysis of spectral count data in label-free shotgun proteomics. Mol Cell Proteomics. 2008, 7: 2373-2385.
Chen C, Ridzon DA, Broomer AJ, Zhou Z, Lee DH, Nguyen JT, Barbisin M, Xu NL, Mahuvakar VR, Andersen MR, Lao KQ, Livak KJ, Guegler KJ: Real-time quantification of microRNAs by stem-loop RT-PCR. Nucleic Acids Res. 2005, 33: e179-
Jung U, Jiang X, Kaufmann SH, Patzel V: A universal TaqMan-based RT-PCR protocol for cost-efficient detection of small noncoding RNA. RNA. 2013, 19: 1864-1873.
Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B. 1995, 57: 289-300.
The authors are grateful to the National Research Foundation Singapore, Singapore-MIT Alliance in Research and Technology, and National University of Singapore for funding this project (research grant C182000036531). We also thank MC Phoon, SH Lau, MX Gu, D Zheng, AN Moorthy, and Connie Foo for their contributions and assistance.
The authors declare that they have no competing interests.
KST helped to design the study, carried out the animal experiments and other analyses, and drafted the manuscript. HC performed the microarray data analyses and helped to draft the manuscript. XJ and VP designed the real-time RT-PCR primers and experiments for the verification. LY performed the immunofluorescence quantification. JES is the pathologist who conducted all histopathologic analyses. BPE participated in the study design and coordination. VTC conceived the study, participated in its design and coordination, and finalized the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Figure S1: Scatter plot showing between-sample correlation of probe-set intensity in infected groups (miRNA and mRNA expression at both 7 and 15 dpi). (TIFF 209 KB)
Additional file 2: Table S4: miRNAs of higher significance based on selection criteria. Implicated functions of 20 shortlisted DE miRNAs. (DOCX 45 KB)
Additional file 3: Table S1: Significantly DE mRNAs targeted by the top 20 shortlisted miRNAs (with directional changes). (XLSX 28 KB)
Additional file 5: Table S2: Sequences (5′-3′) of stem-loop RT primers, real-time specific forward primers and real-time universal reverse primer for miRNA RT-PCR validation. (DOCX 12 KB)
Additional file 6: Table S3: Sequences (5′–3′) of forward and reverse primers for real-time quantitative RT-PCR of selected genes. (DOCX 13 KB)
Authors’ original submitted files for images
About this article
Cite this article
Tan, K.S., Choi, H., Jiang, X. et al. Micro-RNAs in regenerating lungs: an integrative systems biology analysis of murine influenza pneumonia. BMC Genomics 15, 587 (2014). https://doi.org/10.1186/1471-2164-15-587