- Research article
- Open Access
Identification of differentially expressed miRNAs in chicken lung and trachea with avian influenza virus infection by a deep sequencing approach
BMC Genomicsvolume 10, Article number: 512 (2009)
The Erratum to this article has been published in BMC Genomics 2010 11:373
MicroRNAs (miRNAs) play critical roles in a wide spectrum of biological processes and have been shown to be important effectors in the intricate host-pathogen interaction networks. Avian influenza virus (AIV) not only causes significant economic losses in poultry production, but also is of great concern to human health. The objective of this study was to identify miRNAs associated with AIV infections in chickens.
Total RNAs were isolated from lung and trachea of low pathogenic H5N3 infected and non-infected SPF chickens at 4 days post-infection. A total of 278,398 and 340,726 reads were obtained from lung and trachea, respectively. And 377 miRNAs were detected in lungs and 149 in tracheae from a total of 474 distinct chicken miRNAs available at the miRBase, respectively. Seventy-three and thirty-six miRNAs were differentially expressed between infected and non-infected chickens in lungs and tracheae, respectively. There were more miRNAs highly expressed in non-infected tissues than in infected tissues. Interestingly, some of these differentially expressed miRNAs, including miR-146, have been previously reported to be associated with immune-related signal pathways in mammals.
To our knowledge, this is the first study on miRNA gene expression in AIV infected chickens using a deep sequencing approach. During AIV infection, many host miRNAs were differentially regulated, supporting the hypothesis that certain miRNAs might be essential in the host-pathogen interactions. Elucidation of the mechanism of these miRNAs on the regulation of host-AIV interaction will lead to the development of new control strategies to prevent or treat AIV infections in poultry.
Avian influenza virus (AIV) is a type A virus of the family Orthomyxovirida. Although wild aquatic birds such as water fowls and sea gulls are their natural reservoir , land-based birds including chickens may also be infected, which cause significant economic losses to the poultry industry, and raise a great public health threat due to potential host jump from animals to humans .
miRNAs are non-coding, single-stranded RNAs of 19~23 nucleotides which represent a novel class of gene regulators and play important roles in a variety of biological processes in both plants and animals [3–5]. miRNAs modulate gene expression largely at the post-transcriptional level by different mechanisms including direct cleavage of targeted mRNAs , inhibition of translation  or even up-regulation of translation . miRNAs are involved in different biological activities such as development, differentiation, growth and metabolism [8–11]. Recently, in mammals, miRNAs have been reported to participate in the regulation of immunity, including development and differentiation of lymphocytes, monocytes and neutrophils, and modulation of inflammation . MiR-150 expresses in mature B and T cells derived from mouse hematopoietic stem cells, and is able to block early B cell development when expressed prematurely . MiR-181a is an intrinsic modulator of T cell sensitivity and selection in mice . After exposure of THP-1 (human acute monocytic leukemia cell line) cells to lipopolysaccharides (LPS), miR-146 was identified as an inhibitor of signalling proteins of the innate immune responses by NF-kappaB . miRNAs have also been found to be critical effectors in the regulation of viral pathogenesis. Two human encoded miRNAs (miR-136 and miR-507), have been shown to have potential binding sites for the genes that code for the polymerase basic 2 (PB2) and hemmagglutinin (HA) proteins and are reported to be involved in the pathogenesis of H5N1 AIV . All of these evidences suggest that certain miRNAs might be important in the modulation of AIV infections in chickens.
In order to effectively control AIV infection in poultry, it is essential to elucidate the mechanisms of virus pathogenesis in chickens. However, how host cells interact with AIVs during infection in poultry remains poorly understood. Identification of differentially expressed miRNAs in AIV infected chickens will pave a novel avenue to understand host-virus interaction. With the development of next generation sequencing, massively parallel sequencing holds great promise for expression profiling  and it can provide a superior sensitivity at high sequencing depth to discover especially those miRNAs with low abundance and novel miRNAs that are not able to be identified using traditional cloning approaches. Deep sequencing has been previously used to profile both chicken miRNAs and Marek's disease virus miRNAs [17, 18]. In the current study, a Solexa Sequencer was used to deep sequence differentially regulated chicken miRNAs in H5N3 infected and non-infected SPF chickens. Our results will expand the list of miRNAs which might be related to the host immune responses in animals.
Virus replication was examined by real-time RT-PCR for influenza matrix gene from total RNAs of lung and tracheae at 4 dpi. The titer of infected samples was 12.29 log10 EID50/ml in lung, and 3.89 log10 EID50/ml in tracheae. Both non-infected lung and tracheae samples were negative.
Small RNA libraries
A total of 278,398 and 340,726 filtered high quality reads were obtained from chicken lungs and tracheae, respectively (Table 1). In the libraries of chicken lungs, 98,849 and 179,549 reads were obtained from infected and non-infected lungs, respectively. Out of these reads, 52,363 of these high quality reads were exact matches while another 9,357 reads were loose matches to known chicken miRNAs. All reads with a perfect match to mature miRNA sequences from chicken deposited in miRBase http://microrna.sanger.ac.uk/ with insertions or deletions of 1-4 nucleotides was considered as a loose match to represent dicer-processing products from each of the chicken miRNA precursors. An example from gga-mir-181a is shown in Figure 1. Here we saw that the sample had 77 copies of a sequence that was identical to that of the mature gga-mir-181a (denoted by a '*'). In addition we observed significant alternate processing at the 3'-end that was characteristic of miRNAs in various copy numbers. Loose matches were defined by sequence reads that aligned with chicken miRNA consensus sequence with 1-4 mismatches. These may represent sequencing errors (when occurring in low copy numbers), mutations and/or RNA editing events. In the libraries of chicken trachea, 250,460 reads were obtained from infected tracheae and 90,266 were obtained from non-infected tracheae. Out of these reads, 44,243 of these high quality reads were exact matches and another 4,178 reads were loose matches to known chicken miRNAs. The sum of exact and loose match reads was used as the total number of reads for each miRNA.
To display the distribution of miRNAs at each library, reads of individual miRNA within each group were transformed to log10. The plots of distributions of transformed reads for each miRNA from each group are shown in Figure 2. Most miRNAs had around 10-100 reads in these four groups. There was no significant difference in the medians of the four libraries (P > 0.05) indicating they had similar distributions of the miRNA reads.
Of the 474 distinct Gallus gallus (gga) miRNA entries in miRBase [19–22], 377 miRNAs were identified in chicken lungs and 149 miRNAs in chicken tracheae. In addition, we identified 87 potential novel miRNA sequences [see Additional file 1]. Seven novel miRNAs were selected and all of them were confirmed by Real-time quantification PCR in chicken lungs and tracheae, respectively, which is shown in Table 2.
miRNAs expression profiling analysis
In the present study there were four different groups: two tissues (lung and trachea) and two states per tissue (infected and non-infected states). Two comparisons were made between infected vs. non-infected within tissue and another two comparisons between lung and trachea at either infected or non-infected state. Differentially expressed miRNAs in each comparison were identified (P < 0.05, Q < 0.05 and fold change > 2). Within tissues, 73 and 36 miRNAs were differentially expressed between infected and non-infected groups in lungs and tracheae, respectively (Tables 3, 4). When between tissues were compared, 57 and 78 miRNAs were differentially expressed at infected and non-infected states, respectively (Tables 5, 6).
More miRNAs (60 out of 73 miRNAs in lungs and 27 out of 36 miRNAs in tracheae) were down-regulated than up-regulated with AIV infection in both lungs and tracheae. When infected vs. non-infected was compared, 5 miRNAs (miR-106, miR-1729, miR-1798, miR-429 and miR-1711) were only expressed in the non-infected lungs, while 2 miRNAs (miR-1576 and miR-1636) were only expressed in infected lungs (Table 3). Between infected and non-infected tracheae, two (miR-1612 and miR-1458) out of 27 down-regulated miRNAs were expressed only in non-infected tracheae (Table 4).
In the comparisons between tissues, only few miRNAs (6 out of 57 miRNAs in the infected state and 1 out of 78 miRNAs in the non-infected state) were highly expressed in tracheae compared to lungs (Tables 5, 6). Under the infected state, 28 miRNAs were specifically expressed in lungs and 23 miRNAs were expressed at higher levels in lungs than in tracheae. In the non-infected state, 11 miRNAs specifically expressed in lungs and 66 miRNAs were expressed at higher levels in lungs than tracheae. Of particular interest, miR-1a, miR-140, and miR-449, which were highly expressed in infected tracheas than the non-infected ones, and also were differentially expressed between infected tissues (higher expression levels in infected tracheae than infected lungs). In the tissue comparison under the non-infected state, miR-206 was the only miRNA that had higher expression level in tracheae than in lungs. In general, those highly abundant miRNAs were observed across all four groups examined (Table 7).
Confirmation of differentially expressed miRNAs
TaqMan miRNA assays were used to confirm the expression pattern of differentially expressed miRNAs in lungs. There were general consistency between TaqMan assay and deep sequence analysis in three miRNAs (miR-1a, miR-125b and miR-146a) in terms of directions of regulation and significance. Specifically, there was a 1.16 fold up-regulation (2.12 fold in deep sequencing analysis) in miR-1a, 2.13 fold down-regulation (8.33 fold in deep sequencing analysis) in miR-125b, and 3.03 fold down-regulation (3.45 fold in deep sequencing analysis) in miR-146a (P < 0.05).
Clustering of chicken miRNAs
Chromosomal positions of differentially expressed miRNAs revealed that some of them were very close to each other. According to a previous report , miRNAs can be grouped as one cluster if they are less than 1,000 bp apart on the same chromosome. Based on the miRBase 13.0 [19–22], there are 20 miRNA clusters in the chicken genome according to the criteria above. Eighteen of these clusters were detected in lungs and 12 clusters in tracheae, respectively (Table 8). Each cluster contained at least two miRNAs, and total of 47 miRNAs were included in these clusters. Within these clusters, the mir-92-mir-19b-mir-20a-mir-19a-mir-18a-mir-17, which is equivalent to the mammalian mir-17-92 cluster, and mir-302b-mir-302c-mir-1811-mir-302a-mir-302d-mir-367 cluster were the biggest clusters containing six miRNAs. Both of them were detected in lungs. There were only seven clusters differentially expressed (all miRNAs within the cluster differentially expressed in the comparison of infection vs. non-infection or between tissues). Clusters mir-16-1-mir-15a, let-7f-let-7a-1, mir-181a-1-mir-181b-1, let-7j-let-7k, mir-23b-mir-27b-mir-24, and mir-16-2-mir-15b were down-regulated in lungs and mir-181a-1-mir-181b-1 was also down-regulated in tracheae with AIV infection. Cluster mir-34b-mir-34c was up-regulated in tracheae.
Gene ontology analysis
Potential target genes of differentially expressed miRNAs in each comparison were predicted by miRanda . In brief, each differentially expressed miRNA was submitted to miRanda individually and all of its targets predicted in miRanda were used for the following gene ontology (GO) analysis. For each comparison, target genes of induced and repressed miRNAs were separately analyzed. All targets of induced miRNAs in each comparison were submitted to DAVID program  and so were the targets of repressed miRNAs. Functional category enrichment based on the GO terms was evaluated on the targets of these differentially expressed miRNAs. All enriched GO terms in biological process of each comparison were shown in Additional file 2. Immune related GO terms of each comparison are presented in Figures 3, 4, 5 and 6.
In the comparison of infected vs. non-infected lungs, 15 immune related GO terms in biological process were significantly enriched (P < 0.05) (Figure 3) from the targets of down-regulated miRNAs. Response to virus was the GO term with the highest fold enrichment (15 folds). Other functional terms including immune system response, lymphocytes, and lung development were also identified.
In the comparison of infected vs. non-infected tracheae, six immune-related GO terms in biological process were significantly enriched (P < 0.05) (Figure 4) from the targets of down-regulated miRNAs, in which interleukin-12 production and interleukin-12 biosynthetic process had the highest fold enrichment (9.8 folds). For the targets of up-regulated miRNAs between the infected vs. non-infected both in lung and trachea comparisons, only two immune related GO terms: immunoglobulin I set and immunoglobulin subtype 2 in INTERPRO category were significantly enriched.
For the tissue comparison under the infected state, targets of induced miRNAs in lungs were associated with ten immune related GO terms in biological process including tissue development, interleukin-12, lymphocytes cytokines and immune system development (P < 0.05) (Figure 5). Targets of repressed miRNAs in lungs were associated with immune system development and immunoglobulin subtype 2 in the INTERPRO category. Under the non-infected state, 15 immune related terms were significantly enriched in biological process in the targets of miRNAs highly expressed in lungs. These GO terms were related to lung development and host immune system (Figure 6). The GO term NF-kappaB binding was also enriched in molecular function (P < 0.05).
The impact of miRNAs expression on the understanding of molecular mechanisms in gene regulations has been remarkable. Although thousands of small RNAs have been identified over the last decade, the challenge remains to fully identify all small nuclear RNAs, especially very low abundant ones and to determine their individual functions. The majority of known miRNAs have been identified through traditional cloning method, which is both time consuming and labor intensive. The advantages of next-generation sequencing technologies have provided an innovative tool to look into the genome with unprecedented depth of coverage. Solexa deep sequencing is one of these high throughput technologies, by which miRNAs can be detected in any organism without prior sequence or secondary structure information. This technology has been used in many species including human, mice and birds [17, 18, 26–28]. Expression of miRNAs varies in different developmental stages [29–31]. Chicken miRNAs identified in the present study provided novel information in the profiling of miRNAs not only in AIV infected chickens, but also the two tissues (lung and trachea) that have not been previously examined for miRNA profiling in chickens. To our knowledge, this is the first study to profile chicken miRNAs in AIV infected chickens by deep sequencing approach. There are 474 chicken miRNAs predicted in miRBase 13.0 [19–22]. The deep sequencing results in the current study experimentally confirmed 377 miRNAs in chicken lungs and 149 miRNAs in chicken tracheae, and the approach is more powerful than other conventional technologies previously used in birds . The identification of these chicken miRNAs will be very useful in further investigating the functions and regulatory mechanisms of miRNAs in the chicken.
Growing evidence has suggested a relationship between differential miRNA expression and human diseases [32, 33]. miRNAs can regulate many aspects of the immune response, including the development and differentiation of B and T cells, proliferation of monocytes and neutrophils, antibody switching and the release of inflammatory mediators by regulating basic component of host immune system [34–36]. MiR-155 has been reported by several groups to play important roles in both innate and adaptive immune responses in mammals [11, 37, 38]. MiR-155 deficient mice lacked the capability to generate defensive immune responses and to develop lymphocytes, especially B-cell, and antigen-presenting cell functions . However, miR-155 showed very low abundance in both lungs and tracheae and no significant differential expression was observed in the present study. Over expression of miR-181a in mature mouse T cells can augment the sensitivity to peptide antigens, while suppressing miR-181a expression can reduce sensitivity and impair both positive and negative selection . Selective expressions of miR-181a in the thymus and miR-223 in the bone marrow have been shown to be involved in the differentiation of pluripotent hematopoietic stem cells into the various blood cells lineages including B and T cells [40, 41]. In the present study, miR-223 was not significantly regulated while miR-181a was down-regulated in both infected lungs and tracheae. In addition, miR-181a had a higher expression level in lungs than in tracheae under both infected and non-infected states. The expression levels of miR-181a, 181a* and 181b were investigated in LPS activated and CD40-lignad activated macrophages of chickens, respectively . Only miR-181b was expressed in the macrophage cell line HD11 as well in the spleen adherent cells and that its expression increased after activation by LPS or CD40-ligand . In the current study, miR-181b had same expression pattern with miR-181a in both lung and trachea comparisons. These results suggest that miR-181a and miR-181b may be strong miRNA candidates that regulate host response to AIV infection, and warrant further investigation of their targets and regulation mechanism in chickens.
Although the interaction between miRNA expression and virus infection remains to be elucidated, we speculated that miRNA might target immune related genes or modulate virus replication. Sequencing of chicken miRNAs in Marek's disease virus (MDV) infected and non-infected chicken embryo fibroblast (CEF) indicated that more miRNAs were up-regulated in MDV infected cells . These results differ from the current study in which most differentially expressed miRNAs (55 out of 73 in lungs and 27 out of 36 in tracheae) were down-regulated in AIV infected tissues. These results indicate that the mechanisms of miRNA regulation of the host response to different types of virus in chickens are different. Chicken miR-221 and miR-222, the most abundant miRNAs in the CEF small RNA libraries, had significantly higher reads in MDV infected than non-infected CEF . While both miR-221 and 222 had relatively lower abundance in the present study. These results demonstrate that miRNA expression can be tissue-specific with high abundance of miR-221 and 222 in the CEF libraries and low abundances in lungs and tracheae. It can also be speculated that host miRNAs expression may be suppressed by AIV replication based on the miRNA expression patterns observed in the current study.
Some miRNAs have been shown to be directly involved in virus replication. A liver specific miRNA (miR-122) was shown to be required for Hepatitis C virus (HCV) replication in humans . MiR-122 can positively affect the viral replication and has become a therapeutic target for the treatment of HCV infection . In the current study, miR-122 specifically expressed in chicken lungs compared to tracheae under both infected and non-infected states. These data suggest miR-122 might play a more important role in tissue distribution than the responses to AIV infection in chickens. Another two human miRNAs (miR-507 and miR-136) have potential target binding sites in polymerase basic 2 (PB2) and hemagglutinin (HA) genes of AIV, respectively . Unfortunately, these two miRNAs are absent in the chicken genome, which might indicate different infectivity and lethality of the virus between chickens and humans.
Although in the present study most differentially expressed miRNA were down-regulated during AIV infection, some miRNAs were also up-regulated. MiR-1a, miR-140 and miR-449 were significantly up-regulated in both tissues, while miR-455, miR-34b and miR-34c were only up-regulated with AIV infection in tracheae. This suggests different miRNA regulation mechanisms might exist on host response to virus infection. These up-regulated miRNAs might inhibit gene expression of their target genes; therefore down-regulation of these target genes might help the host to inhibit virus replication.
Different tissues serve different biological functions in animals and the expression patterns of miRNAs can vary in different tissues [23, 30]. miRNAs in bursa and spleen of developing chicken embryo have been recently identified, and diverse expression patterns of these miRNAs between different immune organs were observed, suggesting that miRNAs may function as dynamic regulators of the vertebrate immune system . Some miRNAs show tissue-specific distribution in mouse, suggesting specific functions within these tissues . In the current study, chicken lung and trachea were examined, as they are both part of the respiratory system and important sites for AIV replication. There was a significant difference in miRNA expression between lung and trachea with more miRNAs expressed in lungs (377 miRNAs identified) than tracheae (149 miRNAs identified), although only small percentage of miRNAs (19% in lung and 24% in trachea) were significantly differentially expressed in AIV infected samples.
When tissues in the state of virus infection were compared, 28 and 23 miRNAs were specifically and highly expressed in lungs, respectively, and only 6 miRNAs (miR-1a-1 and 2, miR-1b, miR-34b, 34c and miR-449) were highly expressed in tracheae (Table 5). When tissues were compared under the non-infected state, all differentially expressed miRNAs were expressed at higher levels in lungs than tracheae with the only exception of miR-206, which showed a higher expression level in non-infected trachea than lung (Table 6). More interestingly, miR-206 was up-regulated in virus infected vs. non-infected lungs and was down-regulated in infected vs. non-infected tracheae. We can conclude that miR-206 has an opposite regulatory role in lungs and tracheae or might have different targets in different tissues and therefore play different roles in host-virus interactions. MiR-1458 and miR-1612 were up-regulated in AIV infected chicken lungs, while they were specifically expressed in non-infected tracheae not the infected one (Tables 3 and 4). The different regulation of miR-1458 and miR-1612 between lung and trachea suggests they may also have different mechanisms in response to AIV infection between tissues.
We hypothesize that miR-34b, miR-34c, miR-206, miR-1458 and miR-1612 might be some of the most important miRNAs associated with AIV infection. Significantly different miRNA expression pattern between lung and trachea suggests the regulatory mechanism of miRNAs on host response to the AIV infection between lung and trachea is distinct. However, similar regulatory mechanism might also exist in these two tissues. Within the down-regulated miRNAs in infected vs. non-infected lungs and tracheae, there were 18 miRNAs which overlapped in both tissues. This suggests that these 18 miRNAs might have common modulation mechanisms with the AIV infection in chickens.
GO term enrichment analysis has been widely used in functional analysis and allows the identification of important categories associated with functions of interests. GO terms enriched by the target genes of differentially expressed miRNAs in the current study can provide useful information for the follow-up study to elucidate the regulatory mechanism of miRNAs in host immune response to AIV infection. During AIV infection, the host immune system is stimulated to develop a defensive mechanism, which might be the reason why genes involved in immune system development were enriched in all comparisons.
With virus infection, more immune related GO terms were enriched by the targets of repressed miRNAs in lungs than in tracheae (15 terms in lung comparison and 6 terms in trachea comparison) (Figures 3 and 4). Response to virus was identified as the most enriched term (15 fold enrichment) in lung comparison, confirming that genes related to virus infection were regulated by miRNAs. The hyperinduction of proinflammatory cytokines such as TNF-α and IFN-β in human macrophages and respiratory epithelial cells by the highly pathogenic AIV H5N1 was believed to contribute to its high pathogenecity . Lymphocytes were also reported to be suppressed by AIV . Enrichments of T-cell and leukocytes activation and cytokines activities terms identified in the comparison of infected vs. non-infected lungs might be an indication of host immune system response against virus infection. Meanwhile, GO terms involved in lung development and epithelium morphogenesis were enriched, suggesting the genes associated with lung epithelium development in lungs may be important for the recovery from AIV infection in chickens.
It was interesting that two GO terms, interleuklin-12 production and biosynthetic process, were enriched in the infected tissue comparison, which were not included in the non-infected comparison. These two terms were also enriched in the comparison of infected vs. non-infected trachea instead of the lung comparison. IL-12 plays a pivotal regulatory role in the anti-viral response due to its induction of IFN-γ, an anti-viral cytokine . These may suggest that a different defensive mechanism against virus infection might occur in trachea compared to lungs.
The two terms, response to virus and T-cell activation were also enriched by immune related genes differentially expressed in the early immune responses to H9N2 infection in tracheal organ cultures (TOC) . Host immune response, showed as adaptive immune responses in the current study, was enriched by differentially expressed genes in H5N1 infected chicken embryo fibroblasts (CEF) as well . Influenza virus triggered a cascade of both innate and specific immune responses. Then both immune related genes and miRNAs who might regulate these genes maybe involved in similar biological processes with the same GO terms.
Of special note, NF-KappaB binding was also enriched in the comparison between lung and trachea under non-infected state. A similar GO term, regulation of NF-KappaB, was enriched in the previous TOC model with the infection of AIV H9N2 . Activation of NF-KappaB pathway is an essential immediate early step of immune activation. Many viruses have developed strategies to manipulate NF-KappaB signalling through the use of multifunctional viral proteins that target the host innate immune response pathways . Enrichment of GO term NF-KappaB binding suggests these two tissues might utilize this signal pathway differently.
Post-transcriptional gene activity can be regulated through the interaction of regulatory RNA-binding proteins and small non-coding RNAs such as miRNAs [4, 53]. miRNAs can modulate protein activities by altering mRNA stability, translational efficiency or localization [53, 54]. The 3' untranslated regions (3' UTR) are widely accepted as important post-transcriptional regulatory regions of mRNAs, which are particularly rich in cis-acting regulatory elements [55, 56]. miRNAs can regulate their target genes through the cis-acting regulatory elements . miRNAs within the same cluster might share the same cis-regulatory elements , and therefore, might have the same regulatory mechanism for their target genes. Out of the 18 miRNA clusters identified in lungs and 12 miRNA clusters identified in tracheae, there were 7 miRNA clusters differentially expressed in different comparisons. The miRNAs from five of these clusters (mir-16-1-mir-15a, mir-16-2-mir-15b, let-7f-let-7a-1, let-7j-let-7k and mir-23b-mir-27b-mir-24) identified in both lungs and tracheae were significantly down-regulated in infected lungs compared to non-infected lungs and also had higher expression levels in non-infected lungs than non-infected tracheae. The mir-181a-1-mir181b-1 cluster was significantly down-regulated in both infected lungs and tracheae. And the mir-34b-mir-34c cluster was the only significantly up-regulated cluster in the AIV infected trachea. Different miRNA clusters had different regulation direction in AIV infected tissues in the present study. This illustrates that, during AIV infection, different modulation mechanisms among different miRNA clusters might coexist in both lungs and tracheae.
It is interesting to note that when considering the miRNA clusters that were most active in chicken lung and trachea, mir-17-92 cluster (consisting of six miRNAs) and mir-302b-mir-302c-mir-1811-mir-302a-mir-302d-mir-367 cluster are highly associated with cell proliferation and self-renewal of stem cells and cancer cells [58–61]. In addition the miRNAs clusters that were significantly down-regulated miR-15/16 and let-7 are typically down-regulated in stem cells and cancer [62–64]. These results suggest that AIV infection in chickens may instigate cell proliferation and self-renewal like behaviour in chicken lung epithelium and the newly recruited T lymphocytes.
Modulation of target genes by miRNA is one of most critical steps for gene expression regulation. The targeted genes for some differentially expressed miRNAs in the current study were predicted using miRanda [24, 65]. Interestingly, many of the target genes were involved in the host immune system. The potential target genes for miR-1a and miR-1b are the T-cell immuno-modulatory protein. MiR-34b and miR34c, whose target genes are B-cell CLL-pymphoma 2 & 11, might be involved in the B-cell differentiation. Target genes for miR-206 were associated with monocyte macrophage differentiation, suggesting they maybe associated with antigen presentation. Based on other immune related miRNA studies in mammals [11, 66], differentially expressed miRNAs of their mammalian homologs and their targets are presented in Table 9. MiR-15a, miR-21 and miR-181a have important functions in lymphocytes development and modulations while miR-122 and miR-24 are related to virus infection and miR-146a, induced by macrophages, can activate Toll like receptor (TLR) and expose antigens to interleukin-1 beta. Although the exact functions of these miRNAs in the AIV infected chickens remains to be determined, candidate miRNAs and their potential targets identified in the current study provide strong evidence of their roles and warrant further investigation. Whether these chicken miRNAs have the same function as mammals or not need to be validated in the future studies. On-going efforts in the author's laboratory focusing on gene expressions of these target genes and determination of target genes for these differentially expressed miRNAs will provide new insights of miRNA regulations on AIV infection in chickens.
miRNAs have recently been implicated in the intricate cross-talk between host and pathogen in viral infections and are critical in viral pathogenesis. In AIV infected tissues, expression patterns of some host miRNAs were significantly differentially regulated in the current study, supporting the hypothesis that certain miRNAs are essential in the host-pathogen interactions. Once the role of these miRNAs in the regulation of host-AIV interaction has been determined, it will improve the protective strategies in AIV infection in poultry.
Sample collection and RNA isolation
One week old commercial Leghorn SPF chickens were randomly divided into two groups (4 chickens per group), housed in a negative pressure Horsfall-Bauer, temperature control isolation unit, and provided with water and commercial feed ad libitum. At three weeks of age, one group was inoculated with 0.2 ml H5N3 virus containing 106.0 EID50/ml, while the other group was inoculated with PBS by the intra-choanal cleft route. Based on the pilot study at 4 dpi, depression and severely congested lungs and trachea were observed. Therefore, all chickens were euthanized at four days post-inoculation, and lung and trachea epithelial layers were collected for RNA isolation. The animal experiment was performed according to the guidelines approved by the Institutional Animal Care and Use Committee, Texas A&M University.
Two pools of total RNA samples (two random chickens per pool) of each tissue from each group were generated. Total RNAs were isolated using Trizol (Invitrogen, Carlsbad, CA) following the manufacturer's protocol. Dnase I (Ambion, Austin, TX) digestion was carried out after RNA isolation according to manufacturer's instructions. The RNA concentration and purity were determined by measuring absorbance at 260 nm and A260/A280 ratio using a NanoDrop ND-1000 spectrophotometer (Nanodrop Technologies, Wilmington, DE). RNA samples were stored at -80°C until further use.
Virus replication at 4 dpi was determined by real-time RT-PCR for influenza matrix gene using AgPath-ID™ AIV- M kit (Ambion, Austin, TX) following the manufacturer's instructions. Control RNA was extracted from serially diluted H5N3 virus (101.5-105.5 log10 EID50/ml). Standard curve was generated with control viral RNAs. The amount of RNA in the samples was converted into log10 EID50/ml by interpolation as described previously .
Small RNA sequencing and analysis
For small RNA library construction, RNA samples were prepared using the DGE-Small RNA Sample Prep Kit (Illumina, San Diego, CA). In brief, RNA was purified by polyacrylamide gel electrophoresis (PAGE), to enrich for molecules in the range of 18-30 nt, and ligated with proprietary adapters to both 5' and 3' termini of the RNA. Ligated samples were used as templates for cDNA synthesis and then amplified with 15 PCR cycles to produce sequencing libraries. A total of eight Solexa-ready small RNA templates were prepared through two gel purification steps to eliminate concatenated adaptors without inserts. Purified cDNAs were quantified using the Quant-iT PicoGreen dsDNA Kit (Invitrogen, Carlsbad, CA) and diluted to 10 nM for sequencing on an Illumina 1G Genome Analyzer at the Genome Sequencing Center of Baylor College of Medicine. Cluster generation was performed and clusters were sequenced.
For each sample, sequences were first passed through an adaptor filter that searched for sequences that were followed by at least 6 nucleotides of the 3' adaptor. Out of the total reads, any reads without a perfect 10-nt linker subsequence were directly discarded adjoining the insert, yielding of length 10 nt or longer that were subject to further processing. All full-length, exact sequence matches to E. coli (k12, o157:h7, o157:h7 edl933, cft073) were discarded to eliminate possible sequence artifacts arising from the amplification process. All unique sequence reads with a minimum read count of 10 were aligned with precursor chicken miRNA sequences from miRBase version 13.0 [19–22]. Reads of each miRNA were the sum of exact and loose matches (± 4 bp) to known miRNAs. For each sample, counts were normalized to the total number of small RNA sequences, and then for each miRNA, the normalized number of counts was compared between groups or between tissues.
Fisher's Exact test was used to identify differentially expressed miRNAs at a 5% false discovery rate. False discovery rate (FDR) (Q values) was calculated by R program according to Benjamin and Hochberg's method . Ratios were calculated as the ratio of normalized reads of infected over non-infected group or lung over trachea. Statistics related to over representation of functional categories were performed using DAVID, which is based upon a Fisher Exact statistic methodology similar to that described by Al-Shahrour et al . A P < 0.05 was considered as significant.
Novel miRNAs from both lung and trachea libraries were identified using the method by Creighton et al. 2009 . In brief, the first step is to take the sequence reads that did not map to known miRNA precursors, mapped them to the chicken genome, and got an exact sequence match along with 100 bases flanking either side. About 220-bp sequence was then tested for miRNA-like hairpin structure, and folded with the Vienna package . The miRNA hairpin structures that meet the Ambros  criteria were identified. Specifically, the putative miRNA must lie on one arm of a single-loop hairpin with minimum free energy less than -25 kcal/mol. The sequence reads that were appropriately placed in these miRNA-like hairpins were considered as 'putative mature miRNAs' (pmms). Then we examined the pmms for cross-species conservation of the hairpin structure. The sequence reads with strong conservation of the mature miRNA, significant conservation of the hairpin arm opposite the mature miRNA, and little or no conservation of the hairpin loops were considered as novel miRNAs.
Confirmation of novel miRNAs by miScript real-time quantitative RT-PCR
Novel potential miRNA expression was determined by using the same total RNA samples for small RNA library constructions. Total RNA (1 μg) was reverse-transcribed with miScript Reverse Transcription Kit from QIAGEN (Valencia, CA). The real-time quantification of a selected subset of novel miRNAs (C1N1: AAGCUGCCAGUUGAAGAACU; C1N2: AAGGUCCAACCUCACAUGUCC; C1N3: UUGGUGGUUCAGUGGUAGAA; C1N4: AGAAUUGCGUUUGGACAAUC; C1N5: CACAAGAAUUGCGUUUGGACAA; C1N6: UUGACAUCAUCAUACUUGGGAU and C1N7: UGGCAGUGCGUGUUAGCUGGCUGUU) was carried out with customer designed miScript Primer Assays and miScript SYBR Green PCR Kit from QIAGEN (Valencia, CA). Chicken small nuclear RNA U6 (GCAGGGGCCAUGCUAAUCUUCUCUGUAUCG) was used for normalization. The expression levels of novel miRNAs were measured in terms of threshold cycle value (CT) and normalized to U6 using 2-ΔΔCT .
Confirmation of differentially expressed miRNAs by TaqMan MicroRNA Assay
To determine the expression of miRNAs by quantitative RT-PCR (qRT-PCR), TaqMan microRNA assay was performed. The specific stem-loop RT primers of miR-1a, miR125b, miR-146a and U6 were obtained commercially from Applied Biosystems (Foster City, Calif., USA). In brief, cDNA was synthesized from total RNA by using the gene specific primers according to the protocol of TaqMan Micro RNA Assays (Applied Biosystems, CA, USA). Reverse transcriptase reactions contained 10 ng of RNA samples, 3 μl stem loop RT primer and reagents from a TaqMan MicroRNA Reverse Transcription Kit (Applied Biosystems, CA, USA). The 15 μl reactions were incubated for 30 min at 16°C, 30 min at 42°C and 5 min at 85°C, and then held at 4°C.
Real-time PCR was performed by using gene specific probes and a pair of primers (TaqMan MicroRNA Assays, Applied Biosystems) and reagents of TaqMan 2* Universal PCR Master Mix (No AmpErase UNG) (Applied Biosystems, CA, USA). The 20 μl PCR reactions included 1.33 μl RT-PCR product, 10 μl PCR master mix, and 1 μl 20* TaqMan MicroRNA Assay mix (Applied Biosystems, CA, USA). These reactions were incubated at 95°C for 10 min, followed by 40 cycles at 95°C for 10 s, 60°C for 40 s and 72°C for 1 s by ABI 7900 Realtime PCR instrument (Applied Biosystems, CA, USA). All reactions were run in triplicates. The threshold cycle was defined as the fractional cycle number at which the fluorescence passes the fixed threshold. The expression levels of miR-1a, miR-125b and miR-146a in each sample were measured in terms of threshold cycle value and normalized to U6 using 2-ΔΔCT . U6 was used as an internal control.
Webster RG, Bean WJ, Gorman OT, Chambers TM, Kawaoka Y: Evolution and ecology of influenza A viruses. Microbiol Rev. 1992, 56 (1): 152-179.
Webby RJ, Webster RG: Are we ready for pandemic influenza?. Science. 2003, 302 (5650): 1519-1522. 10.1126/science.1090350.
Ambros V: MicroRNA pathways in flies and worms: growth, death, fat, stress, and timing. Cell. 2003, 113 (6): 673-676. 10.1016/S0092-8674(03)00428-8.
Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004, 116 (2): 281-297. 10.1016/S0092-8674(04)00045-5.
Carrington JC, Ambros V: Role of microRNAs in plant and animal development. Science. 2003, 301 (5631): 336-338. 10.1126/science.1085242.
Zhang B, Pan X, Cobb GP, Anderson TA: Plant microRNA: a small regulatory molecule with big impact. Dev Biol. 2006, 289 (1): 3-16. 10.1016/j.ydbio.2005.10.036.
Vasudevan S, Tong Y, Steitz JA: Switching from repression to activation: microRNAs can up-regulate translation. Science. 2007, 318 (5858): 1931-1934. 10.1126/science.1149460.
Guo HS, Xie Q, Fei JF, Chua NH: MicroRNA directs mRNA cleavage of the transcription factor NAC1 to downregulate auxin signals for arabidopsis lateral root development. Plant Cell. 2005, 17 (5): 1376-1386. 10.1105/tpc.105.030841.
Hatfield SD, Shcherbata HR, Fischer KA, Nakahara K, Carthew RW, Ruohola-Baker H: Stem cell division is regulated by the microRNA pathway. Nature. 2005, 435 (7044): 974-978. 10.1038/nature03816.
Lee RC, Feinbaum RL, Ambros V: The C. elegans heterochronic gene lin-4 encodes small RNAs with antisense complementarity to lin-14. Cell. 1993, 75 (5): 843-854. 10.1016/0092-8674(93)90529-Y.
Lindsay MA: microRNAs and the immune response. Trends Immunol. 2008, 29 (7): 343-351. 10.1016/j.it.2008.04.004.
Zhou B, Wang S, Mayr C, Bartel DP, Lodish HF: miR-150, a microRNA expressed in mature B and T cells, blocks early B cell development when expressed prematurely. Proc Natl Acad Sci USA. 2007, 104 (17): 7080-7085. 10.1073/pnas.0702409104.
Li QJ, Chau J, Ebert PJ, Sylvester G, Min H, Liu G, Braich R, Manoharan M, Soutschek J, Skare P, et al: miR-181a is an intrinsic modulator of T cell sensitivity and selection. Cell. 2007, 129 (1): 147-161. 10.1016/j.cell.2007.03.008.
Taganov KD, Boldin MP, Chang KJ, Baltimore D: NF-kappaB-dependent induction of microRNA miR-146, an inhibitor targeted to signaling proteins of innate immune responses. Proc Natl Acad Sci USA. 2006, 103 (33): 12481-12486. 10.1073/pnas.0605298103.
Scaria V, Hariharan M, Maiti S, Pillai B, Brahmachari SK: Host-virus interaction: a new role for microRNAs. Retrovirology. 2006, 3: 68-10.1186/1742-4690-3-68.
Torres TT, Metta M, Ottenwalder B, Schlotterer C: Gene expression profiling by massively parallel sequencing. Genome Res. 2008, 18 (1): 172-177. 10.1101/gr.6984908.
Burnside J, Ouyang M, Anderson A, Bernberg E, Lu C, Meyers BC, Green PJ, Markis M, Isaacs G, Huang E, et al: Deep sequencing of chicken microRNAs. BMC Genomics. 2008, 9: 185-10.1186/1471-2164-9-185.
Glazov EA, Cottee PA, Barris WC, Moore RJ, Dalrymple BP, Tizard ML: A microRNA catalog of the developing chicken embryo identified by a deep sequencing approach. Genome Res. 2008, 18 (6): 957-964. 10.1101/gr.074740.107.
Griffiths-Jones S, Saini HK, van Dongen S, Enright AJ: miRBase: tools for microRNA genomics. Nucleic Acids Res. 2008, D154-158. 36 Database
Griffiths-Jones S, Grocock RJ, van Dongen S, Bateman A, Enright AJ: miRBase: microRNA sequences, targets and gene nomenclature. Nucleic Acids Res. 2006, D140-144. 10.1093/nar/gkj112. 34 Database
Griffiths-Jones S: The microRNA Registry. Nucleic Acids Res. 2004, D109-111. 10.1093/nar/gkh023. 32 Database
Ambros V, Bartel B, Bartel DP, Burge CB, Carrington JC, Chen X, Dreyfuss G, Eddy SR, Griffiths-Jones S, Marshall M, et al: A uniform system for microRNA annotation. RNA. 2003, 9 (3): 277-279. 10.1261/rna.2183803.
Gu Z, Eleswarapu S, Jiang H: Identification and characterization of microRNAs from the bovine adipose tissue and mammary gland. FEBS Lett. 2007, 581 (5): 981-988. 10.1016/j.febslet.2007.01.081.
John B, Enright AJ, Aravin A, Tuschl T, Sander C, Marks DS: Human MicroRNA targets. PLoS Biol. 2004, 2 (11): e363-10.1371/journal.pbio.0020363.
Dennis G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA: DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003, 4 (5): P3-10.1186/gb-2003-4-5-p3.
Friedlander MR, Chen W, Adamidi C, Maaskola J, Einspanier R, Knespel S, Rajewsky N: Discovering microRNAs from deep sequencing data using miRDeep. Nat Biotechnol. 2008, 26 (4): 407-415. 10.1038/nbt1394.
Hafner M, Landgraf P, Ludwig J, Rice A, Ojo T, Lin C, Holoch D, Lim C, Tuschl T: Identification of microRNAs and other small regulatory RNAs using cDNA library sequencing. Methods. 2008, 44 (1): 3-12. 10.1016/j.ymeth.2007.09.009.
t Hoen PA, Ariyurek Y, Thygesen HH, Vreugdenhil E, Vossen RH, de Menezes RX, Boer JM, van Ommen GJ, den Dunnen JT: Deep sequencing-based expression analysis shows major advances in robustness, resolution and inter-lab portability over five microarray platforms. Nucleic Acids Res. 2008, 36 (21): e141-10.1093/nar/gkn705.
Goldsmith ZG, Dhanasekaran N: The microrevolution: applications and impacts of microarray technology on molecular biology and medicine (review). Int J Mol Med. 2004, 13 (4): 483-495.
Xu H, Wang X, Du Z, Li N: Identification of microRNAs from different tissues of chicken embryo and adult chicken. FEBS Lett. 2006, 580 (15): 3610-3616. 10.1016/j.febslet.2006.05.044.
Darnell DK, Kaur S, Stanislaw S, Konieczka JH, Yatskievych TA, Antin PB: MicroRNA expression during chick embryo development. Dev Dyn. 2006, 235 (11): 3156-3165. 10.1002/dvdy.20956.
Jiang J, Gusev Y, Aderca I, Mettler TA, Nagorney DM, Brackett DJ, Roberts LR, Schmittgen TD: Association of MicroRNA expression in hepatocellular carcinomas with hepatitis infection, cirrhosis, and patient survival. Clin Cancer Res. 2008, 14 (2): 419-427. 10.1158/1078-0432.CCR-07-0523.
Jones SW, Watkins G, Le Good N, Roberts S, Murphy CL, Brockbank SM, Needham MR, Read SJ, Newham P: The identification of differentially expressed microRNA in osteoarthritic tissue that modulate the production of TNF-alpha and MMP13. Osteoarthritis Cartilage. 2009, 17 (4): 464-72. 10.1016/j.joca.2008.09.012.
Sonkoly E, Pivarcsi A: Advances in microRNAs: implications for immunity and inflammatory diseases. J Cell Mol Med. 2009, 13 (1): 24-38. 10.1111/j.1582-4934.2008.00534.x.
Sonkoly E, Stahle M, Pivarcsi A: MicroRNAs and immunity: novel players in the regulation of normal immune function and inflammation. Semin Cancer Biol. 2008, 18 (2): 131-140. 10.1016/j.semcancer.2008.01.005.
Turner M, Vigorito E: Regulation of B- and T-cell differentiation by a single microRNA. Biochem Soc Trans. 2008, 36 (Pt 3): 531-533. 10.1042/BST0360531.
Romania P, Lulli V, Pelosi E, Biffoni M, Peschle C, Marziali G: MicroRNA 155 modulates megakaryopoiesis at progenitor and precursor level by targeting Ets-1 and Meis1 transcription factors. Br J Haematol. 2008, 143 (4): 570-580.
O'Connell RMTK, Boldin MP, Cheng G, Baltimore D: MicroRNA-155 is induced during the macrophage inflammatory response. Proc Natl Acad Sci USA. 2006, 104 (5): 6-
Thai TH, Calado DP, Casola S, Ansel KM, Xiao C, Xue Y, Murphy A, Frendewey D, Valenzuela D, Kutok JL, et al: Regulation of the germinal center response by microRNA-155. Science. 2007, 316 (5824): 604-608. 10.1126/science.1141229.
Chen CZ, Li L, Lodish HF, Bartel DP: MicroRNAs modulate hematopoietic lineage differentiation. Science. 2004, 303 (5654): 83-86. 10.1126/science.1091903.
Cobb BS, Hertweck A, Smith J, O'Connor E, Graf D, Cook T, Smale ST, Sakaguchi S, Livesey FJ, Fisher AG, et al: A role for Dicer in immune regulation. J Exp Med. 2006, 203 (11): 2519-2527. 10.1084/jem.20061692.
Ahanda ML, Ruby T, Wittzell H, Bed'Hom B, Chausse AM, Morin V, Oudin A, Chevalier C, Young JR, Zoorob R: Non-coding RNAs revealed during identification of genes involved in chicken immune responses. Immunogenetics. 2009, 61 (1): 55-70. 10.1007/s00251-008-0337-8.
Jopling CL, Yi M, Lancaster AM, Lemon SM, Sarnow P: Modulation of hepatitis C virus RNA abundance by a liver-specific MicroRNA. Science. 2005, 309 (5740): 1577-1581. 10.1126/science.1113329.
Pan QW, Henry SD, Scholte BJ, Tilanus HW, Janssen HL, Laan van der LJ: New therapeutic opportunities for hepatitis C based on small RNA. World J Gastroenterol. 2007, 13 (33): 4431-4436.
Hicks JA, Tembhurne PA, Liu HC: Identification of microRNA in the developing chick immune organs. Immunogenetics. 2009, 61 (3): 231-240. 10.1007/s00251-009-0355-1.
Lagos-Quintana M, Rauhut R, Yalcin A, Meyer J, Lendeckel W, Tuschl T: Identification of tissue-specific microRNAs from mouse. Curr Biol. 2002, 12 (9): 735-739. 10.1016/S0960-9822(02)00809-6.
Hui KP, Lee SM, Cheung CY, Ng IH, Poon LL, Guan Y, Ip NY, Lau AS, Peiris JS: Induction of proinflammatory cytokines in primary human macrophages by influenza A virus (H5N1) is selectively regulated by IFN regulatory factor 3 and p38 MAPK. J Immunol. 2009, 182 (2): 1088-1098.
Wiwanitkit V: Leucocyte and lymphocyte count in cases of bird flu infection in Thailand. J Indian Med Assoc. 2008, 106 (3): 168-186
Thomas JD, Morris KR, Godfrey DI, Lowenthal JW, Bean AG: Expression, purification and characterisation of recombinant Escherichia coli derived chicken interleukin-12. Vet Immunol Immunopathol. 2008, 126 (3-4): 403-406. 10.1016/j.vetimm.2008.08.004.
Reemers SS, Groot Koerkamp MJ, Holstege FC, van Eden W, Vervelde L: Cellular host transcriptional responses to influenza A virus in chicken tracheal organ cultures differ from responses in in vivo infected trachea. Vet Immunol Immunopathol. 2009,
Sarmento L, Afonso CL, Estevez C, Wasilenko J, Pantin-Jackwood M: Differential host gene expression in cells infected with highly pathogenic H5N1 avian influenza viruses. Vet Immunol Immunopathol. 2008, 125 (3-4): 291-302. 10.1016/j.vetimm.2008.05.021.
Hiscott J, Nguyen TL, Arguello M, Nakhaei P, Paz S: Manipulation of the nuclear factor-kappaB pathway and the innate immune response by viruses. Oncogene. 2006, 25 (51): 6844-6867. 10.1038/sj.onc.1209941.
Keene JD: RNA regulons: coordination of post-transcriptional events. Nat Rev Genet. 2007, 8 (7): 533-543. 10.1038/nrg2111.
Chabanon H, Mickleburgh I, Hesketh J: Zipcodes and postage stamps: mRNA localisation signals and their trans-acting binding proteins. Brief Funct Genomic Proteomic. 2004, 3 (3): 240-256. 10.1093/bfgp/3.3.240.
Chen JM, Ferec C, Cooper DN: A systematic analysis of disease-associated variants in the 3' regulatory regions of human protein-coding genes II: the importance of mRNA secondary structure in assessing the functionality of 3' UTR variants. Hum Genet. 2006, 120 (3): 301-333. 10.1007/s00439-006-0218-x.
Chen JM, Ferec C, Cooper DN: A systematic analysis of disease-associated variants in the 3' regulatory regions of human protein-coding genes I: general principles and overview. Hum Genet. 2006, 120 (1): 1-21. 10.1007/s00439-006-0180-7.
Xie X, Lu J, Kulbokas EJ, Golub TR, Mootha V, Lindblad-Toh K, Lander ES, Kellis M: Systematic discovery of regulatory motifs in human promoters and 3' UTRs by comparison of several mammals. Nature. 2005, 434 (7031): 338-345. 10.1038/nature03441.
He L, Thomson JM, Hemann MT, Hernando-Monge E, Mu D, Goodson S, Powers S, Cordon-Cardo C, Lowe SW, Hannon GJ, et al: A microRNA polycistron as a potential human oncogene. Nature. 2005, 435 (7043): 828-833. 10.1038/nature03552.
Wang Y, Lee CG: MicroRNA and cancer--focus on apoptosis. J Cell Mol Med. 2009, 13 (1): 12-23. 10.1111/j.1582-4934.2008.00510.x.
Aguda BD, Kim Y, Piper-Hunter MG, Friedman A, Marsh CB: MicroRNA regulation of a cancer network: consequences of the feedback loops involving miR-17-92, E2F, and Myc. Proc Natl Acad Sci USA. 2008, 105 (50): 19678-19683. 10.1073/pnas.0811166106.
Mendell JT: miRiad roles for the miR-17-92 cluster in development and disease. Cell. 2008, 133 (2): 217-222. 10.1016/j.cell.2008.04.001.
Cho WC: OncomiRs: the discovery and progress of microRNAs in cancers. Mol Cancer. 2007, 6: 60-10.1186/1476-4598-6-60.
Cimmino A, Calin GA, Fabbri M, Iorio MV, Ferracin M, Shimizu M, Wojcik SE, Aqeilan RI, Zupo S, Dono M, et al: miR-15 and miR-16 induce apoptosis by targeting BCL2. Proc Natl Acad Sci USA. 2005, 102 (39): 13944-13949. 10.1073/pnas.0506654102.
Nimmo RA, Slack FJ: An elegant miRror: microRNAs in stem cells, developmental timing and cancer. Chromosoma. 2009, 118 (4): 405-18. 10.1007/s00412-009-0210-z.
Enright AJ, John B, Gaul U, Tuschl T, Sander C, Marks DS: MicroRNA targets in Drosophila. Genome Biol. 2003, 5 (1): R1-10.1186/gb-2003-5-1-r1.
Chen XM, Splinter PL, O'Hara SP, LaRusso NF: A cellular micro-RNA, let-7i, regulates Toll-like receptor 4 expression and contributes to cholangiocyte immune responses against Cryptosporidium parvum infection. J Biol Chem. 2007, 282 (39): 28929-28938. 10.1074/jbc.M702633200.
Lee CW, Suarez DL: Application of real-time RT-PCR for the quantitation and competitive replication study of H5 and H7 subtype avian influenza virus. J Virol Methods. 2004, 119 (2): 151-158. 10.1016/j.jviromet.2004.03.014.
Benjamini YHY: Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. Journal of the Royal Statistical Society Series B (Methodological). 1995, 57 (1): 12-
Al-Shahrour F, Diaz-Uriarte R, Dopazo J: FatiGO: a web tool for finding significant associations of Gene Ontology terms with groups of genes. Bioinformatics. 2004, 20 (4): 578-580. 10.1093/bioinformatics/btg455.
Creighton CJ, Nagaraja AK, Hanash SM, Matzuk MM, Gunaratne PH: A bioinformatics tool for linking gene expression profiling results with public databases of microRNA target predictions. RNA. 2008, 14 (11): 2290-2296. 10.1261/rna.1188208.
Hofacker IFW, Stadler P, Bonhoeffer L, Tacker M, Schuster P: Fast folding and comparison of RNA secondary structure. Chem Month. 1994, 125: 35-10.1007/BF00811671.
Schmittgen TD, Livak KJ: Analyzing real-time PCR data by the comparative C(T) method. Nat Protoc. 2008, 3 (6): 1101-1108. 10.1038/nprot.2008.73.
This project was supported by the Formula Animal Health Grant no. 07-024 and partially supported by the Formula Animal Health Grant no. TEX09244 from the USDA Cooperative State Research, Education, and Extension Service. This research was also partially supported by the National Natural Science Foundation of China (No. 30828024 and 30972156). We are grateful for Drs. Xianyao Li and Hsin-I Chiang for their helps in the animal trial.
YW carried out the RNA isolation, small RNA library construction preparation, analyzed data and drafted the manuscript. VB was responsible for the animal trial. BL and SR contribute to experiment design. BY, HFZ, JHK, JW, and PG contributed to the analysis of miRNA. RC ran the miRNA deep sequencing. HZ provided the concepts of the study, and revised the manuscript. All authors submitted comments, read and approved the final manuscript.
An erratum to this article is available at http://dx.doi.org/10.1186/1471-2164-11-373.