Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Open Access

Peer-reviewed

Research Article

MicroRNome Analysis Unravels the Molecular Basis of SARS Infection in Bronchoalveolar Stem Cells

  • Bibekanand Mallick ,

    Contributed equally to this work with: Bibekanand Mallick, Zhumur Ghosh

    Affiliation Computational Biology Group, Indian Association for the Cultivation of Science, Calcutta, India

    ⨯
  • Zhumur Ghosh ,

    Contributed equally to this work with: Bibekanand Mallick, Zhumur Ghosh

    zhumurg@stanford.edu

    Affiliation Department of Medicine and Department of Radiology, Stanford University School of Medicine, Stanford, California, United States of America

    ⨯
  • Jayprokas Chakrabarti

    Affiliation Computational Biology Group, Indian Association for the Cultivation of Science, Calcutta, India

    ⨯

Abstract

Severe acute respiratory syndrome (SARS), caused by the coronavirus SARS-CoV, is an acute infectious disease with significant mortality. A typical clinical feature associated with SARS is pulmonary fibrosis and associated lung failure. In the aftermath of the SARS epidemic, although significant progress towards understanding the underlying molecular mechanism of the infection has been made, a large gap still remains in our knowledge regarding how SARS-CoV interacts with the host cell at the onset of infection. The rapidly changing viral genome adds another variable to this equation. We have focused on a novel concept of microRNA (miRNA)–mediated host–virus interactions in bronchoalveolar stem cells (BASCs) at the onset of infection by correlating the “BASC–microRNome” with their targets within BASCs and viral genome. This work encompasses miRNA array data analysis, target prediction, and miRNA–mRNA enrichment analysis and develops a complex interaction map among disease-related factors, miRNAs, and BASCs in SARS pathway, which will provide some clues for diagnostic markers to view an overall interplay leading to disease progression. Our observation reveals the BASCs (Sca-1+ CD34+ CD45- Pecam-), a subset of Oct-4+ ACE2+ epithelial colony cells at the broncho-alveolar duct junction, to be the prime target cells of SARS-CoV infection. Upregulated BASC miRNAs-17*, -574-5p, and -214 are co-opted by SARS-CoV to suppress its own replication and evade immune elimination until successful transmission takes place. Viral Nucleocapsid and Spike protein targets seem to co-opt downregulated miR-223 and miR-98 respectively within BASCs to control the various stages of BASC differentiation, activation of inflammatory chemokines, and downregulation of ACE2. All these effectively accounts for a successful viral transmission and replication within BASCs causing continued deterioration of lung tissues and apparent loss of capacity for lung repair. Overall, this investigation reveals another mode of exploitation of cellular miRNA machinery by virus to their own advantage.

Citation: Mallick B, Ghosh Z, Chakrabarti J (2009) MicroRNome Analysis Unravels the Molecular Basis of SARS Infection in Bronchoalveolar Stem Cells. PLoS ONE 4(11): e7837. https://doi.org/10.1371/journal.pone.0007837

Editor: I. King Jordan, Georgia Institute of Technology, United States of America

Received: August 17, 2009; Accepted: October 12, 2009; Published: November 13, 2009

Copyright: © 2009 Mallick et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Funding: This work is funded by Department of Science and Technology, Government of India. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Severe Acute Respiratory Syndrome (SARS) is a new fulminant atypical pneumonia which emerged as a regional and global threat in 2002–2003 with a high mortality rate resulting from acute lung failure [1]. The disease causing agent has been identified as a novel coronavirus termed as SARS-associated coronavirus (SARS-CoV) [2], [3]. The SARS-CoV is an enveloped virus containing a single stranded, positive-sense RNA genome which encodes 14 putative open reading frames encoding 28 potential proteins [4], [5]. These include four structural proteins, spike (S) glycoprotein, matrix (M) protein, small envelope (E) protein, and nucleocapsid (N) protein [4]. These proteins have various roles in aiding the virus to enter the host and spread infection. While the incidence of new cases of SARS waned in 2003–2004, many aspects of SARS disease pathogenesis and host-pathogen interactions remain unsolved. Limited pathologic studies reveal that the major site of SARS-CoV infection and morbidity is the respiratory tract. The target organ of SARS is mainly lungs [6]. Among the various animal models that have been used to study the pathogenesis of SARS-CoV infection, the monkey model mimics the clinical course of SARS to a certain degree [7]. The cellular tropism of SARS-CoV in mouse lung has also been investigated by Ling et al., 2006 [8]. But still today, very little is known regarding the mode of SARS-CoV interaction with host cells at the onset of infection in the lungs at a molecular level, and also the cell types in which the primary viral infection and replication occurs.

The lung is an extremely complex, conditionally renewing organ which contains anatomically and functionally distinct epithelial stem cell populations which reside in distinct anatomical locations. The basal cells [9], Clara cells [10] and type-II pneumocytes [11] are the candidate stem/progenitor cells which can repair the injured lungs and contribute to local needs in times of tissue damage. Recently, Kim et al. have isolated a regional pulmonary stem cell population termed as bronchoalveolar stem cells (BASCs) residing at the bronchoalveolar duct junction of adult lungs [12]. These have been identified and characterized as CD34+ Sca-1+ CD45- PE-CAM− cells expressing both cytoplasmic Clara cell secretion protein (CCSP) and prosurfactant protein-C proteins, which are markers for Clara cells and type-II pneumocytes respectively [12]. We have searched for the evidence to establish BASC as the prime target cells of SARS infection initiation and replication.

To date, miRNA-mediated RNA interference is reported to be an essential tool to understand the regulatory pathways at molecular level that underlie infection biology [13], [14]. They not only participate in executive decisions but also perform much of the grunt work to micromanage protein output. Inspite of being such potential regulatory elements, miRNAs have so far missed attention in the pathogenesis of SARS-CoV infection.

MiRNAs are small (19–25 nucleotides) endogenous noncoding RNAs that have been shown to influence the abundance and translational efficiency of cognate mRNAs [15], [16]. Viruses, which typically employ many components of the host gene expression machinery, also encode miRNAs. On the contrary, miRNA-biogenesis pathway poses some serious problem for RNA viruses and a group of DNA viruses (poxviruses) to encode miRNAs [14]. SARS-CoV being an RNA virus is assumed not to encode its own miRNAs. But in due course of evolution, the virus might have developed certain highly sophisticated molecular mechanisms to exploit the cellular biosynthetic machinery of host cells and elude the cellular defense mechanisms. It is most likely that in certain situations SARS-CoV uses the cellular miRNA transcripts in order to foster their own agenda.

Given the proposed involvement of BASCs in SARS disease [8], elucidation of the miRNA-mediated regulatory mechanisms responsible for making BASCs the prime targets of SARS infection might provide new avenues to explore the underlying mechanisms of SARS infection.

Onset of SARS infection

The entry of SARS-CoV into cells is mediated through interaction between spike (S) glycoprotein of the virus and angiotensin-converting enzyme 2 (ACE2), the primary receptor of SARS-CoV on the host cell [17]. ACE2 is a homolog of angiotensin-converting enzyme (ACE), which plays an important role in the renin–angiotensin system for blood pressure homeostasis. Recent studies of Li et al. [18] have shown a protective effect of ACE2 against experimental lung fibrosis through its ability to degrade local tissue angiotensin II (ANG II) in response to bleomycin. Furthermore, in the pathogenesis of lung fibrosis, they have proposed the involvement of the down regulation of ACE2 as an integral component of the sequence of events leading to lung collagen deposition. The molecular mechanisms responsible for the loss of ACE2 in lung fibrosis are currently under investigation. A typical clinical feature associated with SARS is also pulmonary fibrosis and associated lung failure whose underlying mechanism still remains elusive [19]. Investigations leading to the underlying molecular mechanisms responsible for the loss of ACE2 in pulmonary fibrosis might provide a clue to the reason behind fibrosis associated lung injury in SARS.

Several reports suggest that there is a close interaction between SARS-CoV and respiratory epithelia which play an important role in the genesis of SARS [20]. Jia et al. [21] have investigated interactions between SARS-CoV and human airway epithelia and indicated that the state of cell differentiation and ACE2 expression levels are important determinants of the susceptibility of human airway epithelia to infection. To date, the molecular switches regulating the state of cellular differentiation of an infected cell and ACE2 expression have not been studied properly. In this work, we have proposed the role of miRNA in modulating the expression of ACE2 along with the stage specific cellular differentiation at the site of infection initiation. Further we have tried to focus on the role of SARS-CoV which dictates these events for its successful transmission and replication.

The POU-homeodomain transcription factor Oct-4 (Pou5f1) plays a central role in self-renewal, pluripotency, and lineage commitment. Alterations in Oct-4 expression promote differentiation and leads to the specification of ectodermal, endodermal or mesodermal primitive progenitors [22]. A small number of Oct-4-expressing cells have been observed at the bronchoalveolar junction of the neonatal lung [8]. We have investigated miRNA-mediated modulatory role on Oct-4+ BASCs and have correlated their contribution to SARS infection.

Our analysis establishes a complex interplay between miRNAs, BASCs and certain essential factors related to the disease which might provide a better understanding of SARS pathogenesis. The multifarious relationships shared by genes related to the disease-system pathway, modulatory effects of miRNAs within the prime target cells of infection in response to indefinite cues impose the need of an interaction map. This study develops a complex interaction map between disease-related factors, miRNAs and BASCs in SARS disease pathway, which will provide some clues for diagnostic markers to view an overall interplay leading to disease progression. Further, our investigation towards investigating the therapeutic potential of miRNA-mediated RNAi as an effective antiviral agent against SARS might unravel some of the efficient therapeutic measures against the deadly disease in future.

Results and Discussion

We have investigated the cellular tropism of SARS-CoV (Table 1) in the BASCs present in lungs and attempted to establish the role of BASC-miRNAs towards understanding the pathogenesis of SARS in lung fibrosis. We describe the role of BASCs as the target cells of SARS infection and show how a repertoire of BASC miRNAs act as molecular switches to contribute to this host-pathogen interaction.

thumbnail
Download:
Table 1. SARS-CoV strains considered in our study.

https://doi.org/10.1371/journal.pone.0007837.t001

Primary site of SARS infection

Which is the primary site of SARS infection and how do these corona virus interacts with those cells to enter the host and eventually invades the host immune system to spread itself - is a pertinent question which prevails and is a debated issue. In order to throw light on this very important question, we looked upon at the SARS+ cell markers and BASC markers (specific to each of its developmental stages from multipotent state to fully differentiated state).

It has been reported that SARS-infected cells are unlikely to be pneumocytes (both type I or type II) or cytokeratin+ epithelial cells and it is quite obvious since SARS-CoV cannot infect or replicate in a fully differentiated cell [23]. Further, SARS+ cells are distinct from cells expressing the macrophage/monocyte specific marker CD68. But SARS+ cells are found to express the functional receptor ACE2 as well as the stem/progenitor cell marker CD34 and Oct-4 [23]. Chen et al have also shown that SARS+ cells are CD45−. Hence it has been confirmed that the SARS+ cells in the infected lung were a subset of putative stem/progenitor cells expressing CD34, Oct-4 and ACE2. Further, it is unlikely for SARS+ cells to produce differentiated cell markers. The presence of SARS+ cells within the multiple cell types at the bronchiolar lining layers makes it difficult to isolate SARS+ cells. Some of the SARS+ cells remains overlapped with adjacent cytokeratin+ cells which might be mistakenly interpreted as colocalized cells. In such a situation BASC (Sca-1+ CD34+ CD45− Pecam−) in terminal bronchioles located exclusively at the broncho-alveolar duct junction (BADJ) stands a good selection to be investigated as the prime target cells of SARS-CoV infection.

Characterizing BASC at BADJ as the prime target cell of SARS infection initiation

A percentage of the Sca-1+ CD4− Pecam− cell populations at BADJ are CD34+. This Sca-1+ CD34+ CD45− Pecam− population is enriched for BASCs [12] became evident from the fact that it contained no ciliated cells. In addition to clonal colony formation, these cells exhibited extensive self-renewal in culture. They also had a greater capacity for differentiation compared to other lung epithelial cells. Sca-1+ CD34+ CD45− Pecam− BASC cultures further confirmed the multilineage differentiation capacity of BASCs. Within seven to ten days they show up positive for differentiated cell markers as they get differentiated into Clara-like cells (CCA+ SP-C−), AT2-like cells (SP-C+ CCA) and AT1-like cells (AQ5+).

Within the Sca-1+ CD4− Pecam− cells the other percentage of cells showed CD34− and contained ciliated cells. The Oct-4+ colony cells which are Sca-1+ CD34− SSEA-1+ cytokeratin+ were shown to be succeptible to SARS-CoV infection by Ling et al. [8]. These might be a subset of Sca-1+ CD4− Pecam− CD34− population (Figure 1). Later studies have shown that cytokeratin+ cells do not express the functional receptor ACE2 [23]. Further, SARS+ cells do not express cytokeratin and are CD34+ SSEA-1−. Phagocytosis is a common mode of virus entry within surrounding cells in SARS infection [24]. Hence, presence of SARS antigen within the subset of Oct-4+ colony cells which are Sca-1+ CD34− SSEA-1+ cytokeratin+ is most likely due to phagocytosis. Thereby, it is evident that these Oct-4 cells (Sca-1+ CD34− SSEA-1+ cytokeratin+) are not the prime target cells for SARS-CoV infection.

thumbnail
Download:
Figure 1. Lower respiratory tract focusing on bronchoalveolar duct junction (BADJ) of the lung.

Venn diagram represents an enlarged section of the BADJ showing up the different subsets of the distinct cellular population (Oct-4+ ACE2+) in this region of the lung. A subset of this cellular population is enriched in BASCs identified as Sca-1+ CD45− Pecam− CD34+, the other subset being Sca-1+ CD45− Pecam− CD34−. On differentiation Sca-1+ CD45− Pecam− CD34+ cells show positive for CCA and SP-C; Sca-1+ CD45− Pecam− CD34− cells show positive for CCA, SP-C, and ciliated cells.

https://doi.org/10.1371/journal.pone.0007837.g001

On the contrary, the BASCs (Sca-1+ CD34+ CD45− Pecam−) must be another subset of these Oct-4+ ACE2+ epithelial colony cells (Figure 1) which provides a more favorable environment for SARS-CoV entry and replication. Absence of ciliated cells and CD34+ makes them the most probable targets of SARS-CoV at the onset of infection.

miRNAs as modulators of BASC differentiation

Since SARS-CoV cannot infect or replicate in a fully differentiated cell, it is indeed essential for the virus to control the differentiation stages in an infected BASC cells so that they do not reach a fully differentiated state until the virus has undergone successful transmission and replication within this primary site of infection. Hence, it is important to monitor the expression levels of the different developmental stage specific markers of SARS infected BASC and also that of the corresponding regulatory switches that control their level of expression.

MiRNAs are widespread agents of post-transcriptional gene silencing and have been strongly linked with stem cells [25], [26]. They exhibit a high degree of stage- and tissue-specificity, and therefore it is important to scan the profiling data for those miRNAs that operate during the narrow windows of development of a BASC in order to understand their regulatory impact on the expression levels of essential marker genes and transcription factors within the cell. In order to model the regulatory mechanisms within a SARS+ cell, we have tried to understand the underlying regulatory mechanisms within a normal BASC cells mediated by the molecular switches- the miRNAs.

It has been established that the virus exploits the host-miRNA milieu to foster their own agenda of evading the host immune system, thereby establish a consistent infection and continue replication.

mh–miRNAs targeting host factors

We undertook an intensive transcriptome-wide search for candidate mh-miRNAs (BASC-miRNAs which are homologous between mouse-human whole genomes) (Table 2) targeting the developmental stage specific transcriptional factors and marker genes of BASCs. Table 3 shows the list of the miRNAs targeting the potential stage specific factors of BASCs. Further, we have also elucidated the mh-miRNAs targeting ACE2. The target factors have further been mapped with the expression levels of their corresponding miRNAs within BASCs which shows their integrated correlation contributing towards each developmental stage of BASCs in a healthy individual (Figure 2). SARS-CoV definitely interferes with these miRNA mediated regulation to dictate the cell fate and continue its replication. Hence as a future advancement of this work, it will be worthy designing benchwork to detect the correlation of these miRNA-mRNA expressions within a SARS-CoV infected cell.

thumbnail
Download:
Figure 2. Fold change of the differentially expressed host miRNAs in BASC targeting host and virulent viral factors [61].

https://doi.org/10.1371/journal.pone.0007837.g002

thumbnail
Download:
Table 2. Significantly expressed mh–miRNAs in bronchio-alveolar stem cells.

https://doi.org/10.1371/journal.pone.0007837.t002

thumbnail
Download:
Table 3. Predicted miRNA–mRNA pairs within the host system.

https://doi.org/10.1371/journal.pone.0007837.t003

Host miRNA targets in viral genome

We have searched for the host miRNA targets within the virulent genes of SARS-CoV. Table 4 includes the BASC-miRNA candidates targeting the most important virulent proteins responsible for viral infection. Interestingly, miRNAs-17*, -574-5p and -214 targets all the four viral virulent proteins viz. S, N, M, E and orf1a. miR-148a has its target in ORF1a, E, S and M. miR-223 and miR-98 has an exclusive correlation with their targets within the 3′ UTRs of N and S respectively (Figure 2).

thumbnail
Download:
Table 4. Host miRNAs targeting virulent protein coding mRNAs of SARS-CoV.

https://doi.org/10.1371/journal.pone.0007837.t004

Normally presence of such host miRNA targets in SARS-CoV genome indicates that the corresponding miRNAs are a part of the host's innate antiviral defense. On the contrary, the RNA virus co-opt host-miRNAs to suppress their own replication to evade immune elimination and establish a strong infection. It has been found that human infecting single-stranded RNA viruses are enriched for target sequences of human miRNAs [42]. Or else rapidly evolving RNA viruses might evolve targets that match host miRNAs in order to increase their host specificity. It is very likely that the rapidly evolving virus take advantage of slowly evolving host miRNAs to increase their survival in the host, analogous to the manner in which viruses take advantage of other cellular processes. Co-option of host miRNAs by SARS-CoV to suppress their own replication will allow it to be in a latent state and escape the host immune system at the initial time point of transmission and infection. Once successful transmission within the host system takes place, antiviral miRNA target sequences in the virus can be expected to mutate rapidly to maximize mismatches and thus minimize the impact of ‘antiviral’ miRNAs. Rapid mutation in SARS-CoV is supported by the fact that SARS-CoV lack proof reading activity in their polymerases. RNA-dependent polymerase in SARS-CoV misincorporates 8.26 bases per million [43].

While investigating the seed sequence conservation of the viral-specific host miRNAs (Figure 3) we found that there is a strong bias towards G, Y(pYrimidine) and R (PuRine) bases at the 3rd, 4th and 5th position respectively of the seed region. This implies that a single mutation within the complementary regions of this conserved 3-mer motif within the viral targets will disrupt the miRNA seed-target match and will enable the virus to escape the antiviral effect of host-miRNAs after successful transmission within the host. Such host-miRNA sequence conservation helps the virus to mutate the target regions with more precision and at ease so as to escape the effect of all the 4 miRNAs together.

thumbnail
Download:
Figure 3. Sequence alignment between the common viral specific host miRNA mature sequences from 5′ to 3′ end.

The seed region (position 2 to 7) and the conserved positions (3rd, 4th, and 5th) within the seed region are highlighted within boxes.

https://doi.org/10.1371/journal.pone.0007837.g003

The miRNA expression profiles of three miRNAs-17*, -574-5p and -214 targeting S, N, Orf1a, M and E shows that they are upregulated within uninfected BASCs. Presumably at the onset of SARS infection the expression of these targets is suppressed which is likely to control viral replication and consequently help the virus to evade the host immune system. But the targets for downregulated miR-148a within ORF1a, M, S and E might escape repression and maintain replication at a lower rate. Further, the exclusive miRNAs viz. miR-223 and -98 for N and S protein respectively are highly downregulated (8.7 fold and 5.6 fold respectively) in BASCs. N and S proteins are important for viral transmission in host. This provides a clue towards the fact that the S protein and N protein (having targets corresponding to these miRNAs) of SARS-CoV takes advantage of this and escapes miRNA-mediated repression which is a rescue to the virus for effective transmission at the initial stage of viral infection.

Key miRNA–mRNA pairs correlating with virus entry, replication, and host–virus interaction

ACE2-S interaction.

It has been shown that proinflammatory and fibrogenic cytokine pathways are activated within the first 24–48 hr following pulmonary insult resulting from SARS-CoV infection [44]. High initial levels of these activities are associated with persistent pulmonary damage and increased risks of subsequent pulmonary fibrosis and poor outcome in diffuse alveolar damage. Further, it has been reported that ACE2 mRNA, protein and enzymatic activity is severely decreased in pulmonary fibrosis resulting from SARS-CoV infection [17], [18]. Studies have also been performed which support that differential regulation of host cell mRNA transcription and protein synthesis occurs in a SARS-CoV infected cells [45].

SARS-CoV S protein is a multifunctional protein which plays pivotal roles in the biology and pathogenesis of SARS-CoV. It has been shown that S protein mediates viral infection by binding to cellular receptor ACE2 and thus induces membrane fusion. SARS-CoV infections and the Spike protein of the SARS-CoV reduce ACE2 expression [46]. Further a recent work of Imai et al. showed that on injecting SARS-CoV Spike into mice worsens acute lung failure in vivo, which can be attenuated by blocking the renin–angiotensin pathway. This suggests that the S-protein mediated activation of the pulmonary RAS influences the pathogenesis of ALI/ARDS and SARS [47]. But the molecular mechanism by which S protein down regulates ACE2 in lung cells following SARS injury is not well addressed till to date.

In our work we observed that the miRNAs- miR-186,-93 and -200b targeting ACE2 are found to be highly down regulated in normal BASCs thereby showing that they do not interfere with the expression level of ACE2 and thus prevents lung fibrogenesis by limiting the local accumulation of the profibrotic peptide ANG II in an uninfected state. Considering the expression of S-protein, we find that its target corresponding to highly downregulated miR-98 probably escapes miRNA-mediated host RNAi defense at the onset of infection and interferes with ACE2 expression. Figure 4 shows how the virus takes advantage of the microRNA-mediated complex interplay between the host and viral factors within BASCs. The question still remains regarding the molecular mechanism by which this S-protein represses the expression of ACE2.

thumbnail
Download:
Figure 4. Modelled miRNA–mediated host–virus interaction within BASC in SARS-CoV pathogenesis.

Host miRNAs are central to viral pathogenesis. This figure illustrates how the host miRNA and SARS-CoV interactions can explain the features of SARS-CoV pathogenesis such as (A) entry of the virus in an undifferentiated BASC, (B) causing acute infection during the proliferative stage and onset of BASC differentiation, and (C) complete clearance of virus in a fully differentiated BASC.

https://doi.org/10.1371/journal.pone.0007837.g004

S-proteins have been shown to act as transcription factor activating the promoters of GRP94/78 genes inducing ER stress in a SARS-CoV infected cell [48]. The role of S protein as a transcription factor has obliged us to think that they might play a role to activate the miRNA promoters of miR-186,-93 and -200b targeting ACE2 to enhance their expression resulting in a reduced expression of ACE2 in an infected host. Whether this decrease in ACE2 expression is a result of miRNA mediated gene silencing controlled by S-protein will be an interesting and compelling question for further study.

Oct-4-N interaction.

Oct-4 plays a crucial role towards regulating the chromatin structure in a state consistent with self-renewal and later facilitating the expression of genes that keeps the cell poised to respond to cues that lead to differentiation. Hence this is one of the essential factors which control the state of differentiation in BASCs. Furthermore, it has a significant role in apoptosis. miR-146b-5p is seen to target Oct-4 in normal BASCs. The downregulated expression of this miRNA explains for the uninterrupted expression of Oct-4 resulting its normal function in an uninfected cell. We have harnessed the high possibility of a complicated crosstalk between the ability of N-protein to control the expression of Oct-4 in an infected cell.

N protein is known to be the most abundantly expressed protein of SARS-CoV. The most unique and significant property of it is revealed in its ability to act as a sequence specific DNA binding factor. It has been reported to bind NFkB response element of COX2 promoter and enhance COX2 gene expression [49] which effectively is beneficial for its existence within the host. This reveals the power of the N-protein to manipulate the entire gene expression programme of the infected cell. The common miRNAs, miR-214,-574-5p and -17* having targets in N protein are upregulated in BASCs. But miR-223, the unique miRNA targeting N-protein is highly downregulated. This signifies that the N-protein targets corresponding to miR-223 is going to escape miRNA-mediated downregulation at the onset of infection (See Figure 4). Further it seems that N-protein plays a role in interfering with the expression modulation of Oct-4 by binding to the promoter of miR-146-5p and controls the different stages of BASC differentiation as well as prevents apoptosis of SARS infected cells.

Immunity-related issues—induction of inflammatory chemokines in SARS-CoV–infected cells might be an advantage for the virus to replicate

All respiratory virus disease results from two concurrent pathological components: ongoing virus replication and the resulting inflammatory response. Even if antivirals clearly inhibit virus replication, the biochemical and cellular inflammatory responses to the initial infection related events continue despite diminished virus titer [50]. While acute inflammatory responses are generally beneficial in nature and have been shown to limit virus replication in situ, prolonged, uncontrolled inflammation has been recognized as a significant component contributing to the pathological sequelae to the morbidity and mortality of severe acute respiratory syndrome coronavirus infection [44]. A histopathologic analysis of SARS-CoV infected lung in mice showed that there has been an induction of significant changes in the expression levels of the chemokine receptor CCR1 which acts as receptors for CCL3 and CCL5. The expression levels of these inflammatory chemokines were increased by viral infection of the lung [33]. It has been shown that CCR1 play an essential role in the pathogenesis of pulmonary fibrosis in lung-injury models [34]. Further, previous work suggests that antagonism of signaling via the chemokine receptor CCR1 is a potent strategy that attenuates leukocyte recruitment in response to severe respiratory virus infection [50].

Although there has been a subsequent increase in inflammatory cytokines related to innate immune response in response to viral infection, there has been a lack of immunoregulatory cytokines (T or B lymphocytes) in the lung during infection which confirms the fact that the expression of viral replication proteins remains at a quiescent state so as to avoid viral clearance by the activation of immunoregulatory cytokines by adaptive immune responses at the onset of the infection. On the contrary, activation of inflammatory chemokines might be acting as a beneficial phenomenon for the replication of various viruses [51], [52].

Emerging studies suggest that miRNAs provide an added layer in orchestrating immune responses [53]. miRNAs function in shaping immunity by regulating the repertoire of genes expressed in immune cells and the magnitude and duration of responses to particular pathogens [54]. Hence we investigated the underlying miRNA induced molecular mechanisms responsible for these changes in the lung and we found that miR-223 has its target in the 3′-UTR of CCR1. Previously, miR-223 was experimentally identified in chorioamnionitis-related inflammation and its predicted target genes include several genes involved in inflammation and immune responses. Further, miR-223 is known to promote granulocyte differentiation and is a part of a regulatory loop that involves C/EBP and NFI-A [55].

Interestingly miR-223 has its target in N protein of SARS-CoV too and it is one of the miRNAs which probably helps the virus to escape the RNAi mediated silencing at the onset of infection. Or in another way we may say that the virus encodes the target of miR-223 to escape the RNAi mediated repression and at the same time takes advantage of the effect of induction of the inflammatory chemokine responses to accelerate lung fibrosis. But the controlled expression of target CCR1 in an uninfected cell inspite of the significantly low expression of its corresponding miR-223 in BASCs suggest that there must be a cascaded fine tuned interaction masterminded by the virus which is going on within the infected cell. It is known that N protein activates NF-KB in Vero E6 cells [38]. CCR1 is one of the target genes of NF-KB which is activated at a later stage [56]. Hence, we propose that within BASCs, reduced expression of miR-223 helps N protein to enter the host cell. On successful entry, N protein might be using miR-223 once again to activate CCR1 via NF-KB at its replicative stage to enhance lung fibrosis and at the same time gets empowered to manipulate the entire gene expression of the host cell by controlling the expression of Oct-4. The exact mechanism of such miRNA-mediated events needs further investigation. Further, checking the expression levels of N-protein, CCR1 as well as miR-223 in an infected BASC will provide a more compact conclusion.

Conclusion

Based on these findings, we propose that BASCs (Sca-1+ CD34+ CD45− Pecam-) are a subset of Oct-4+ ACE2+ cells and are the chief targets of SARS-CoV infection. SARS-CoV enters the host; infects the BASCs of the lungs and controls its developmental stages via the molecular switches-miRNAs. Co-option of host miRNAs by viruses reveal their intelligent plan to control their replication in order to evade immune elimination until they undergoes successful transmission and establish a strong infection. Thereafter, they undergo rapid mutation to maximize the target-miRNA mismatches and enhance their replication before the cell reaches at a fully differentiated state. On successful replication, the virus infects resident, infiltrating, and circulating immune cells. The circulating immune cells carry the virus to other organs and causes damages to the immune cells of spleen, peripheral and central lymph nodes and other lymphoid tissues. The immune defense being weakened significantly, leads to rapid deterioration of the pneumonia. It is worth to test these cells for the expression of Oct-4 and ACE2 for further confirmatory conclusion.

The proposed host miRNA-dependent mechanism probably acts in concert with the host and viral factors for establishing a strong infection. Our proposal shifts away from the simplistic notion that host miRNAs with specific sequence complementarity to viruses are indicative of a bonafide innate antiviral immune mechanism to an established viral pathogen.

Although the immediate threat of SARS is over, we intend to continue with its therapeutic measures partly to guard for its possible return. Moreover, this will help us to transform our knowledge from SARS to other emerging viral diseases specially causing lung injuries which most certainly will be a deadly threat in future. Lot of attempts have been made and is still going on to design effective antiviral agents against this deadly disease [57]–[60]. While some of these approaches showed partial efficiency to combat SARS infection, the others need further clinical investigations to be used as a proper therapeutic reagent.

This era has witnessed an increase in the development of exogenous siRNA therapeutics against viruses. The potential therapeutic effect of modulating host-miRNA levels is also worth considering. This emerging picture of miRNA regulation in SARS-CoV infected individual will make the host-virus interaction far richer and more complex than the crisp linear pathways of the previous decade, with miRNAs participating in executive decisions at the interesting and vulnerable nodes in regulatory networks. By boosting the level of the host virus-specific miRNAs, it might be possible to turn what was a desirable outcome for the virus into a desirable outcome for the patient. Under the current hypothesis, we would predict that decreasing the level of host virus-specific miRNAs would result in promotion of viral replication, leading to immunological recognition and clearance of the virus by the host immune system. Further, increasing the level of miRNAs exclusive for the N and S protein might also block viral entry. However, manipulating the level of host miRNAs could have unintended consequences because the physiological functions of the miRNAs might be altered or viral pathology might be enhanced. Further silencing of miRNAs with antagomirs specifically within the infected cell could become another therapeutic strategy where miRNAs participate in disease aetiology. Nevertheless, these potential interventions merit further evaluation.

Materials and Methods

To identify the miRNAs involved in the cellular tropism of SARS-CoV in lungs, we undertook an intensive search for potential miRNA candidates in the BASCGAP (Broncho-alveolar stem cell genome anatomy project) sequence data. Our miRNA sequences were derived from the Sanger Institute miRBase release 12.0 (http://microrna.sanger.ac.uk/sequences). The global miRNA expression profile of our BASCGAP source sequences were derived from mouse BASCs [61].

miRNA—mouse-human (mh–miRNA)

Our search strategy was motivated by the need to apply pure-homology based selection of miRNAs so that the results represent a closer approach to mimic the molecular mechanisms that underlie disease pathway behaviors and response in a human host. Further the overall miRNA expression profile of mouse and human lungs are similar [62]. From the BASCGAP sequence data of 116 miRNAs, we screened out the BASC-miRNAs lying at the intersection of mouse-human whole genomes. This filter is based on a homology search (SSEARCH), performed with miRBase (release 12.0) and is applied on the mature candidates from the BASCGAP sequence data. E-value cutoff considered for SSEARCH was <10.0. Mature miRNA sequence homology was one of the stringent parameters considered in our search. Table 2 represents the screened output set of 95 mature candidates termed as mh-miRNAs.

mRNA targets

We have considered the developmental stage specific marker proteins of BASCs, which includes the transcription factor Pou5f1; ACE2-the most important receptor of SARS-CoV (Table 3). In addition to these, we have chosen the virulent protein coding genes, viz. Nucleocapsid (N), Spike(S), Envelope (E), Membrane (M) and Orf1a from 30 strains of SARS-CoV (Table 1 and Table 4) for our investigation. Annotated UTR sequences of the target genes were retrieved from UTRdb [63].

Target prediction

We adopted RNAhybrid [64], that predicts the most favorable hybridization sites between miRNA and UTR regions and generates energy minimized duplex structures using the Dynamic Programming technique. The resultant miRNA-mRNA duplexes are consistent with the following structural and energy constraints at the first step of screening the targets: (i) Binding energy cutoff < = −20 kcal/mol [64]; (ii) complementarity of minimum seven bases to the 5′ end of the miRNA; (iii) minimization of GU base pair within the miRNA-target duplex, a maximum of only one GU base pair between 2 and 7 position of the miRNA is allowed; (iv) minimization of loop size throughout the duplex.

Although, currently available target prediction methods are diverse, both in approach and performance and all have room for further improvement, perfect seed pairing centered on nucleotides 2–7 is considered as the most important criteria for prediction reliability. Furthermore, in order to screen out the most favorable target in the case of miRNAs having multiple targets we incorporated the constraint of site efficacy (hierarchy of site efficacy follows- 8 mer>7 mer-m8>7 mer-A1>6 mer) following the work of Grimson et al. and Friedman et al. [65], [66]. 3′ compensatory pairing of miRNA to its target have been considered to compensate for imperfect/weak seed matches [67]. We have further checked the effect of cooperativity as well as the UTR binding location.

miRNA–mRNA enrichment analysis

For bantam miRNAs, good minimum free energies (MFEs) can occur frequently by chance. The longer a putative target sequence, the better such random energies will be. Hence, statistical significance of predicted targets is assessed with an extreme value statistics of length normalized minimum free energies and a Poisson approximation of multiple binding sites.

Microarray data analysis in BASCs

miRNA expression profile of 2×106 mouse BASCs (CD45−CD31−CD34+ Sca-1+) and 4×106 control cells(CD45−CD31−CD34−Sca-1−) are determined using miRNA microarray. The data considered for our analysis is obtained by first subtracting the background and then normalizing the signals using a LOWESS filter (locally weighted regression). Background is determined using a regression-based background mapping method [68]. The miRNA dataset was filtered according to the standard procedure to exclude spots with minimum intensity and size. Further for the dual-sample experiments, the ratio of the two sets of detected signals (log2 transformed, balanced) are calculated. Thereafter, t-values are calculated for each miRNA between the test and the control groups, and P-values are computed from the theoretical t-distribution. miRNAs with P-values below a critical P-value (typically 0.01) are the differentially detected signals and were selected as our global miRNA expression profile data set (of mouse BASCs) in BASCGAP. The flowchart of methodology is provided in Figure 5.

thumbnail
Download:
Figure 5. miRNA–mRNA prediction pipelines.

This flow chart summarizes the steps followed and yields from BASCGAP consortium sequences to predict the miRNA-mRNA pairs involved in SARS pathogenesis.

https://doi.org/10.1371/journal.pone.0007837.g005

Acknowledgments

We thank the reviewers for their useful comments and suggestions to improve the manuscript. We also thank our colleagues for critical reading of the manuscript.

Author Contributions

Conceived and designed the experiments: BM ZG. Performed the experiments: BM ZG. Analyzed the data: BM ZG. Contributed reagents/materials/analysis tools: JC. Wrote the paper: BM ZG.

References

  1. 1. Donnelly CA, Ghani AC, Leung GM, Hedley AJ, Fraser C, et al. (2003) Epidemiological determinants of spread of causal agent of severe acute respiratory syndrome in Hong Kong. Lancet 361: 1761–1766.
  2. 2. Ksiazek TG, Erdman D, Goldsmith CS, Zaki SR, Peret T, et al. (2003) A novel coronavirus associated with severe acute respiratory syndrome. N Engl J Med 348(20): 1953–1966.
  3. 3. Poutanen SM, Low DE, Henry B, Finkelstein S, Rose D, et al. (2003) Identification of severe acute respiratory syndrome in Canada. N Engl J Med 348(20): 1995–2005.
  4. 4. Marra MA, Jones SJ, Astell CR, Holt RA, Brooks-Wilson A, et al. (2003) The genome sequence of the SARS-associated coronavirus. Science 300: 1399–1404.
  5. 5. Rota PA, Oberste MS, Monroe SS, Nix WA, Campagnoli R, et al. (2003) Characterization of a novel coronavirus associated with severe acute respiratory syndrome. Science 300: 1394–1399.
  6. 6. Gu J, Gong E, Zhang B, Zheng J, Gao Z, et al. (2005) Multiple organ infection and the pathogenesis of SARS. J Exp Med 202(3): 415–424.
  7. 7. Kuiken T, Fouchier RA, Schutten M, Rimmelzwaan GF, van Amerongen G, et al. (2003) Newly discovered coronavirus as the primary cause of severe acute respiratory syndrome. Lancet 362(9380): 263–270.
  8. 8. Ling TY, Kuo MD, Li CL, Yu AL, Huang YH, et al. (2006) Identification of pulmonary Oct-4+ stem/progenitor cells and demonstration of their susceptibility to SARS coronavirus (SARS-CoV) infection in vitro. Proc Natl Acad Sci USA 103(25): 9530–9535.
  9. 9. Hong KU, Reynolds SD, Watkins S, Fuchs E, Stripp BR (2004) Basal cells are a multipotent progenitor capable of renewing the bronchial epithelium. Am J Pathol 164(2): 577–588.
  10. 10. Evans MJ, Cabral-Anderson LJ, Freeman G (1978) Role of the Clara cell in renewal of the bronchiolar epithelium. Lab Invest 38: 648–653.
  11. 11. Fehrenbach H (2001) Alveolar epithelial type II cell: defender of the alveolus revisited. Respir Res 2: 33–46.
  12. 12. Kim CF, Jackson EL, Woolfenden AE, Lawrence S, Babar I, et al. (2005) Identification of bronchioalveolar stem cells in normal lung and lung cancer. Cell 121(6): 823–835.
  13. 13. Lu M, Zhang Q, Deng M, Miao J, Guo Y, et al. (2008) An analysis of human microRNA and disease associations. PLoS ONE 3(10): e3420.
  14. 14. Ghosh Z, Mallick B, Chakrabarti J (2009) Cellular versus viral microRNAs in host-virus interaction. Nucleic Acids Res 37(4): 1035–1048.
  15. 15. Bartel DP (2004) MicroRNAs: genomics, biogenesis, mechanism, and function. Cell 116(2): 281–297.
  16. 16. Ghosh Z, Chakrabarti J, Mallick B (2007) miRNomics-The bioinformatics of microRNA genes. Biochem Biophys Res Commun 37(4): 1035–1048.
  17. 17. Li W, Moore MJ, Vasilieva N, Sui J, Wong SK, et al. (2003) Angiotensin-converting enzyme 2 is a functional receptor for the SARS coronavirus. Nature 426(6965): 450–454.
  18. 18. Li X, Molina-Molina M, Abdul-Hafez A, Uhal V, Xaubet A, et al. (2008) Angiotensin converting enzyme-2 is protective but downregulated in human and experimental lung fibrosis. Am J Physiol Lung Cell Mol Physiol 295(1): L178–185.
  19. 19. Zhao X, Nicholls JM, Chen Y-G (2008) Severe Acute Respiratory Syndrome-associated coronavirus nucleocapsid protein interacts with Smad3 and modulates transforming growth factor-{beta} signaling. J Biol Chem 283: 3272–3280.
  20. 20. Yu ITS, Li Y, Wong TW, Tam W, Chan AT, et al. (2004) Evidence of airborne transmission of the Severe Acute Respiratory Syndrome virus. New Eng J Med 350: 1731–1739.
  21. 21. Jia HP, Look DC, Shi L, Hickey M, Pewe L, et al. (2005) ACE2 receptor expression and severe acute respiratory syndrome coronavirus infection depend on differentiation of human airway epithelia. J Virol 79(23): 14614–14621.
  22. 22. Campbell PA, Perez-Iratxeta C, Andrade-Navarro MA, Rudnicki MA (2007) Oct4 targets regulatory nodes to modulate stem cell function. PLoS ONE 2(6): e553.
  23. 23. Chen Y, Chan VS, Zheng B, Chan KY, Xu X, et al. (2007) A novel subset of putative stem/progenitor CD34+Oct-4+ cells is the major target for SARS coronavirus in human lung. J Exp Med 204(11): 2529–2536.
  24. 24. Yilla M, Harcourt BH, Hickman CJ, McGrew M, Tamin A, Goldsmith CS, Bellini WJ, Anderson LJ (2005) SARS-coronavirus replication in human peripheral monocytes/macrophages. Virus Research 107: 93–101.
  25. 25. Shcherbata HR, Hatfield S, Ward EJ, Reynolds S, Fischer KA, et al. (2006) The MicroRNA pathway plays a regulatory role in stem cell division. Cell Cycle 5: 172–175.
  26. 26. Zhang B, Pan X, Anderson TA (2006) MicroRNA: A new player in stem cells. J Cell Physiol 209: 266–269.
  27. 27. McQualter JL, Brouard N, Williams B, Baird BN, Sims-Lucas S, et al. (2009) Endogenous fibroblastic progenitor cells in the adult mouse lung are highly enriched in the sca-1 positive cell fraction. Stem Cells 27(3): 623–633.
  28. 28. Martin U (2008) Methods for studying stem cells: Adult stem cells for lung repair. Methods 45(2): 121–132.
  29. 29. Eramo A, Lotti F, Sette G, Pilozzi E, Biffoni M, et al. (2008) Identification and expansion of the tumorigenic lung cancer stem cell population. Cell Death Differ 15: 504–514.
  30. 30. Kim CF (2007) Paving the road for lung stem cell biology: bronchioalveolar stem cells and other putative distal lung stem cells. Am J Physiol Lung Cell Mol Physiol 293: 1092–1098.
  31. 31. Nielsen S, King LS, Christensen BM, Agre P (1997) Aquaporins in complex tissues. II. Subcellular distribution in respiratory and glandular tissues of rat. Am J Physiol 273: C1549–1561.
  32. 32. Oudit GY, Imai Y, Kuba K, Scholey JW, Penninger JM (2009) The role of ACE2 in pulmonary diseases–relevance for the nephrologist. Nephrol Dial Transplant 24(5): 1362–1365.
  33. 33. Glass WG, Subbarao K, Murphy B, Murphy PM (2004) Mechanisms of host defense following severe acute respiratory syndrome-coronavirus (SARS-CoV) pulmonary infection of mice. J Immunol 173(6): 4030–4039.
  34. 34. Tokuda A, Itakura M, Onai N, Kimura H, Kuriyama T, et al. (2000) Pivotal role of CCR1-positive Llukocytes in bleomycin- induced lung fibrosis in mice. J Immunol 164: 2745–2751.
  35. 35. Gallagher TM, Buchmeier MJ (2001) Coronavirus spike proteins in viral entry and pathogenesis. Virology 279: 371–374.
  36. 36. Wilson L, Mckinlay C, Gage P, Ewart G (2004) SARS coronavirus E protein forms cation-selective ion channels. Virology 330(1): 322–331.
  37. 37. Chana C-M, Maa C-W, Chan WY, Chan HYE (2007) The SARS-Coronavirus Membrane protein induces apoptosis through modulating the Akt survival pathway. Archiv Biochem Biophys 459(2): 197–207.
  38. 38. Milan S, Lal SK (2008) The SARS-CoV nucleocapsid protein: a protein with multifarious activities. Infection, genetics and evolution 8(4): 397–405.
  39. 39. Yan X, Hao Q, Mu Y, Timani KA, Ye L, et al. (2006) Nucleocapsid protein of SARS-CoV activates the expression of cyclooxygenase-2 by binding directly to regulatory elements for nuclear factor-kappa B and CCAAT/enhancer binding protein. Int J Biochem Cell Biol 38(8): 1417–1428.
  40. 40. Snijder EJ, Bredenbeek PJ, Dobbe JC, Thiel V, Ziebuhr J, et al. (2003) Unique and conserved features of genome and proteome of SARS-coronavirus, an early split-off from the coronavirus group 2 lineage. J Mol Biol 331(5): 991–1004.
  41. 41. Prentice E, McAuliffe J, Lu X, Subbarao K, Denison MR (2004) Identification and characterization of severe acute respiratory syndrome coronavirus replicase proteins. J Virol 78(18): 9977–9986.
  42. 42. Watanabe Y, Kishi A, Yachie N, Kanai A, Tomita M (2007) Computational analysis of microRNA-mediated antiviral defense in humans. FEBS Lett 581(24): 4603–4610.
  43. 43. Ye R, Montalto-Morrison C, Masters PS (2004) Genetic analysis of determinants for spike glycoprotein assembly into murine coronavirus virions: distinct roles for charge-rich and cysteine-rich regions of the endodomain. J Virol 78(18): 9904–9917.
  44. 44. Nicholls JM, Poon LLM, Lee KC, Ng WF, Lai SK, et al. (2003) Lung pathology of fatal severe acute respiratory syndrome. Lancet 361: 1773–1778.
  45. 45. Xiao X, Chakraborti S, Dimitrov AS, Gramatikoff K, Dimitrov DS (2003) The SARS-CoV S glycoprotein: Expression and functional characterization. Biochem Biophys Res Commun 312: 1159–1164.
  46. 46. Kuba K, Imai Y, Rao S, Gao H, Guo F, et al. (2005) A crucial role of angiotensin converting enzyme 2 (ACE2) in SARS coronavirus-induced lung injury. Nat Med 11(8): 875–879.
  47. 47. Imai Y, Kuba K, Penninger JM (2008) The discovery of angiotensin-converting enzyme 2 and its role in acute lung injury in mice. Exp Physiol 93(5): 543–548.
  48. 48. Jin DY, Zheng BJ (2009) Roles of spike protein in the pathogenesis of SARS coronavirus. Hong Kong Med J 15: Suppl 2S37–40.
  49. 49. Yan X, Hao Q, Mu Y, Timani KA, Ye L, et al. (2006) Nucleocapsid protein of SARS-CoV activates the expression of cyclooxygenase-2 by binding directly to regulatory elements for nuclear factor-kappa B and CCAAT/enhancer binding protein. Int J Biochem Cell Biol 38(8): 1417–1428.
  50. 50. Bonville CA, Easton AJ, Rosenberg HF, Domachowske JB (2003) Altered pathogenesis of severe pneumovirus infection in response to combined anti-viral and specific immunomodulatory agents. J Virol 77: 1237–1244.
  51. 51. Lee B-N, Ordonez N, Popek EJ, Lu J-G, Helfgott A, et al. (1997) Inflammatory cytokine expression is correlated with the level of human immunodeficiency virus (HIV) transcripts in HIV-infected placental trophoblastic cells. J Virol 1: 3628–3635.
  52. 52. Chang J, Renne R, Dittmer D, Ganem D (2000) Inflammatory cytokines and the reactivation of Kaposi's sarcoma-associated herpesvirus lytic replication. Virol 266(1): 17–25.
  53. 53. Asirvatham AJ, Gregorie CJ, Hu Z, Magner WJ, Tomasi TB (2008) MicroRNA targets in immune genes and the Dicer/argonaute and ARE machinery components. Mol Immunol 45: 1995–2006.
  54. 54. Baltimore D, Boldin MP, O'Connell RM, Rao DS, Taganov KD (2008) MicroRNAs: new regulators of immune cell development and function. Nat Immunol 9: 839–845.
  55. 55. Fazi F, Rosa A, Fatica A, Gelmetti V, De Marchis ML, et al. (2005) A minicircuitry comprised of microRNA-223 and transcription factors NFI-A and C/EBPalpha regulates human granulopoiesis. Cell 123: 819–831.
  56. 56. Law HK, Cheung CY, Sia SF, Chan YO, Peiris JM, et al. (2009) Toll-like receptors, chemokine receptors and death receptor ligands responses in SARS coronavirus infected human monocyte derived dendritic cells. BMC Immunol 10(1): 35.
  57. 57. Cinati J, Morgenstern B, Bauer G, Chandra P, Doerr HW, et al. (2003) Treatment of SARS with human interferons. Lancet 362: 293–294.
  58. 58. Wu CJ, Jan JT, Chen CM, Hsieh HP, Hwang DR, et al. (2004) Inhibition of Severe Acute Respiratory Syndrome Coronavirus Replication by Niclosamide. AAC 48(7): 2693–2696.
  59. 59. Wu CJ, Huang HW, Liu CY, Hong CF, Chan YL (2005) Inhibition of SARS-CoV replication by siRNA. Antiviral Research 65: 45–48.
  60. 60. Li BJ, Tang Q, Cheng Du, Qin C, Xie FY, et al. (2005) Using siRNA in prophylactic and therapeutic regimens against SARS coronavirus in Rhesus macaque. Nat Med 11(9): 944–951.
  61. 61. Qian S, Ding J, Xie R, An JH, Ao XJ, et al. (2008) MicroRNA expression profile of bronchioalveolar stem cells from mouse lung. Biochem Biophys Res Commun 377: 668–673.
  62. 62. Williams AE, Moschos SA, Perry MM, Barnes PJ, Lindsay MA (2007) Maternally imprinted microRNAs are differentially expressed during mouse and human lung development. Dev Dyn 236(2): 572–580.
  63. 63. Mignone F, Grillo G, Licciulli F, Iacono M, Liuni S, et al. (2005) UTRdb and UTRsite: a collection of sequences and regulatory motifs of the untranslated regions of eukaryotic mRNAs. Nucleic Acids Res 33: D141–146.
  64. 64. Rehmsmeier M, Steffen P, Höchsmann M, Giegerich R (2004) Fast and effective prediction of microRNA/target duplexes. RNA 10: 1507–1517.
  65. 65. Grimson A, Farh KK, Johnston WK, Garrett-Engele P, Lim LP, et al. (2007) MicroRNA targeting specificity in mammals: determinants beyond seed pairing. Mol Cell 27: 91–105.
  66. 66. Friedman RC, Farh KK, Burge CB, Bartel DP (2009) Most mammalian mRNAs are conserved targets of microRNAs. Genome Res 19: 92–105.
  67. 67. Bartel DP (2009) MicroRNAs: Target recognition and regulatory functions. Cell 136: 215–233.
  68. 68. Bolstad BM, Irizarry RA, Astrand M, Speed TP (2003) A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics 19: 185–193.