Skip to main content
  • Methodology article
  • Open access
  • Published:

Comprehensive annotation of microRNA expression profiles

Abstract

Background

MicroRNAs (miRNAs) regulate many biological processes by post-translational gene silencing. Analysis of miRNA expression profiles is a reliable method for investigating particular biological processes due to the stability of miRNA and the development of advanced sequencing methods. However, this approach is limited by the broad specificity of miRNAs, which may target several mRNAs.

Result

In this study, we developed a method for comprehensive annotation of miRNA array or deep sequencing data for investigation of cellular biological effects. Using this method, the specific pathways and biological processes involved in Alzheimer’s disease were predicted with high correlation in four independent samples. Furthermore, this method was validated for evaluation of cadmium telluride (CdTe) nanomaterial cytotoxicity. As a result, apoptosis pathways were selected as the top pathways associated with CdTe nanoparticle exposure, which is consistent with previous studies.

Conclusions

Our findings contribute to the validation of miRNA microarray or deep sequencing results for early diagnosis of disease and evaluation of the biological safety of new materials and drugs.

Background

MicroRNAs (miRNAs) are short ribonucleic acid (RNA) molecules with an average length of 22 nucleotides (nt), which exhibit higher stability than messenger RNAs (mRNAs) [1, 2]. They are post-transcriptional regulators that bind to complementary sequences on target mRNA transcripts, usually resulting in translational repression or target degradation and gene silencing [35]. MiRNAs regulate numerous biological processes, including cell viability, proliferation, development and differentiation [1, 2]. Similar to mRNA microarray techniques, methods for studying miRNA expression profiles have been developed including deep sequencing techniques [6, 7]. In addition to the evaluation of the stability of miRNA during sample processing, the assessment of differential miRNA expression profiles has been identified as a reliable method for the investigation of mRNAs, proteins and mechanistic pathways involved in particular biological processes, such as differentiation, carcinogenesis and cytotoxicity. Indeed, methods have been developed to validate miRNA microarray or deep sequencing data for biological research purposes, such as the prediction of miRNA function and activity, interaction of miRNA and mRNAs and the investigation of miRNA regulatory influences on sub-pathways [810].

Generally, miRNAs shown to be expressed at significantly different levels by miRNAs array/sequencing are selected out through fold-change analysis. The target genes of selected miRNA are predicted by tools such as PicTar, myMIR, TargetScan and miRanda [1115]. Subsequently, these predicted genes are enriched in KEGG pathway or Gene Ontology (GO) analyses [1618]. The KEGG pathway database records networks of molecular interactions in cells and the GO analysis provides the ontology of defined terms that represent gene product properties. Three domains are covered by GO: biological processes, molecular functions and sets of molecular events with a defined beginning and end (Figure 1A).

Figure 1
figure 1

The global repression effect of a miRNA expression profile on a specific miRNA. Generally, for miRNA microarray or deep sequencing data analysis, miRNAs (circle) with significantly changed expression profiles are selected through fold-change analysis. The target mRNAs (triangle) of the selected miRNA are predicted by miRNA target prediction tools such as miRanda (A). Theoretically, a mRNA may be targeted by several miRNAs with different context scores, which are regarded as the repression ratio score of a miRNA for its target (B). Thus, the total repression ratio score of the total miRNA expression profile for a specific mRNA can be calculated based on the combined repression effects of the miRNAs on their target mRNA (C, see formula (1)).

However, miRNAs exhibit broad specificity and may target several mRNAs in a given cell. The methods for miRNA target prediction had been intensively studied [1922] and the tools available generate several scores for different potential targets of a particular miRNA. This reflects the complex regulatory networks of miRNAs and mRNAs. Thus, the repression of a particular mRNA by a given set of miRNAs is mediated by the combined effects of each individual miRNA on this common target mRNA [6]. Due to this combinatorial repression effect, the abundance of each miRNA will also contribute to the repression capacity of a given set of miRNAs (Figure 1A).

Based on this study, a mathematical model has been proposed that can comprehensively predict genes that are effectively regulated by a given miRNA. This is achieved by integrating the effect of each miRNA on their target mRNA based on the combined repressive effects of the relevant miRNAs. This model can be used to elucidate the combined effects of a miRNA profile.

Methods

Data source

MiRNA microarray data from four Alzheimer’s disease subjects were obtained from the NCBI Gene Expression Omnibus (GEO), as reported previously [23]. In this study, postmortem human brain samples were obtained and RNA was extracted from parietal lobes of postmortem brains of Alzheimer’s disease patients and controls. The mRNA array measurements were performed at the UCLA microarray core using Affymetrix HG-U133 Plus 2.0 arrays. MicroRNAs were assayed by LC Sciences using a custom -Paraflo array containing probes for 470 miRNAs from Sanger miRBase and 419 miRNAs predicted by miRNAMap.

MiRNA deep sequencing data were obtained in a previous study that evaluated CdTe nanoparticle cytotoxicity [6]. In this study, NIH/3 T3 cells were exposed to cadmium telluride quantum dots (CdTe QD) to extract the small RNAs. After exposure to nanomaterials for 24 h, the cells were harvested with the trypsin to extract the miRNAs. Then, the small RNAs in a sample were converted into a double-stranded cDNA library. The results of SOLiD sequencing were in the form of nucleotide sequences and their coverage. The registered miRNAs were screened out by comparing them in GenBank (http://www.ncbi.nlm.gov/genbank/) and miRbase (http://www.mirbase.org/).

Methods

A mathematical model to evaluate the comprehensive repression rate of specific mRNAs using total miRNA expression profiles

The repression (Re) of mRNA is directly proportional to the inhibitory effect (IE) of miRNA. For instance, for a specific mRNA, mRNAα, regulated by specific miRNAs, the formula is as follows:

Re mRNA α i = a n IE i α
(1)

Where IE i-α indicate the inhibitory effect of a specific miRNAi (i = a, b, c,……n) to mRNAα.

For a particular miRNA, miRNAa, IE a is directly proportional to the distribution of miRNAa (DU a ) on its target mRNAα (represented as DU a-α in the following formula) and the repression score (RS) of miRNAa for its target mRNAα, which can be obtained from miRanda or other tools:

IE a α DU a α × RS a α
(2)

DU a-α is directly proportional to the abundance of miRNAa (AU a ) in the given miRNA array or sequencing database and the ratio of miRNAa combined with mRNAα (CO a,α ):

DU a α AU a × CO a α
(3)

A single miRNAa may target several mRNAs. By induction of Pa − α as the proportion of miRNAa combined with mRNA α , COa − α can be calculated as:

CO a α = P a α i = α n P a i
(4)

Scoring in prediction tools is designed to reveal the ability of miRNA to bind complementary regions of mRNAs [24] and previous work has shown that contexts of 7-nt or 8-nt matches appear sufficient for miRNA-like regulation [25, 26]. Studies have concluded that additional recognition features, such as pairing with the remainder of the miRNA, accessible mRNA structures and protein-binding sites are usually dispensable or occur so frequently that they impart little overall specificity [27]. Hence,

P a α RS a α
(5)

In formula (3), AU a indicates the abundance of miRNAa in the given miRNA array or sequencing database. We designated C a as the counts of miRNAa in the sequencing data (or the signal value in microarray data), and TC miRNAs as the total counts (or the signal value for microarray analysis) of miRNA provided by the sequencing data. Thus, AU a can be calculated as:

AU a = Ca TCmiRNAS
(6)

By combining formulae (3), (5) and (6), DU miRNAa can be calculated as:

DU a α Ca TCmiRNAs × RS a α i = α n RS a i
(7)

Hence, according to formula (2), IE miRNAa can be calculated as:

IE a α Ca × RS a α 2 TCmiRNAs × i = α n RS a i
(8)

By combining formulae (8) and (1), and induction of K as a coefficient factor, the total regression of a given miRNA sequencing data to mRNAα can be calculated as:

Re mRNA α = K × i = a n Ca × R S i α 2 TCmiRNAs × i = α n RS a i
(9)

Statistical analysis

Z-tests were performed to investigate differences in the repression of target mRNAs by miRNAs. p0 was calculated using formula (10), which was used to investigate the null hypothesis.

p 0 α = Re mRNA α Control + Re mRNA α Test TC miRN A Control + TC miRN A Test
(10)

Here, RemRNAα Control and RemRNAα Test indicate the repression of mRNAα in the control and test groups, whereas TCmiRNA Control and TCmiRNA Test indicate the total counts of miRNA in the control and test groups.

Then, the Z-test was performed according to the following formula:

Z a = Re mRN A Test Re mRNA α Control p 0 α 1 p 0 α 1 TC miRN A Control + 1 TC miRN A Test
(11)

The null hypothesis was rejected at Z-values >2.58 or < −2.58, which indicated significant differences between the repression of mRNAs in the test and control groups.

Results and discussion

The global repression effect of a miRNA expression profile on a specific mRNA

It has been reported that an individual miRNA s may target different mRNAs. Indeed, context scores for ranking the predicted targets of each miRNA have been provided by previous studies [10, 28, 29]. In contrast, one mRNA may be targeted by several miRNAs with different context scores, which are represented by the repression ratio score of an individual miRNA for its target (Figure 1B) [30]. To validate miRNA microarray or deep sequencing results for the prediction of changes in proteins or mechanistic pathways, the total repression ratio score of the total miRNA expression profile for a specific mRNA was calculated based on the combined miRNA repression effects on their target mRNA (Figure 1C, see formula (1)).

The repression capacity of an individual miRNA on its target mRNA is affected by miRNA concentration and the repression score identified using programs such as PicTar. MiRNAs bind targets with different efficiencies and therefore, the distribution of an individual miRNA on the specific target should be considered. As scoring in prediction tools is designed to reveal the capacity of miRNA to bind complementary regions of mRNAs, the repression ability of a miRNA is directly proportional to the square of the repression score as shown in formula (9). Formula (9) shows that, for a given miRNA, the efficiency of binding to a given miRNA affects the repression capacity in an exponential manner. Based on the total repression effect of a given miRNA expression profile on a specific mRNA, the significantly regulated mRNA can be selected out by Z-tests for further analysis.

Prediction of biological pathway and process regulation using miRNA expression profiles based on array data

In order to validate our method for identification of specific pathways or biological processes regulated by a given miRNA profile, the significantly regulated mRNA selected by our algorithm was enriched in KEGG pathways or GO terms.

Data from four Alzheimer’s disease subjects and one control subject were used in this study. In order to apply KEGG pathway or GO terms analyses for the precise and comprehensive elucidation of the effects of miRNAs on biological processes in Alzheimer’s disease, the repression effect of Alzheimer’s disease miRNA profiles on a specific mRNA were calculated according to the repression value obtained using miRNA target prediction tools and the interaction properties of miRNA and mRNA. Formula (9) gives a total repression score of a specific mRNA, which is directly proportional to the abundance of the related miRNA and the square of the repression value obtained using miRNA target prediction tools. For example, in sample S3, beta-catenin mRNA (ENTREZ_GENE_ID:1499; NM_001904) is regulated by several miRNAs with different repression scores given by miRanda, such as hsa-miR-139, hsa-miR-200a and hsa-miR-320. The total repression rate calculated using formula (9) was 10,296 in S3 compared with 5,096 in the control group. Based on the repression rate, a Z-test was performed for selection of mRNAs for further analysis.

The selected mRNAs were enriched into KEGG by using the web-based GO analysis tool, DAVID and it was shown that the four Alzheimer’s disease samples shared most of the KEGG pathways that were found to be significantly altered under miRNA regulation (Figure 2A). Among these significantly regulated KEGG pathways, four pathways (WNT signaling pathway, MAPK signaling pathway, axon guidance and pathways involved in cancer) were highly regulated by miRNA in all of the four Alzheimer’s disease samples (Table 1). Other pathways, such as endocytosis, focal adhesion, neurotrophin signaling pathway and regulation of the actin cytoskeleton also showed a significant difference between the four Alzheimer’s disease samples and the control samples (P < 10-5). Studies have illustrated the close correlation of these pathways with Alzheimer’s disease [3138]. Several other selected pathways (P < 0.01) are also shown in Figure 3. Most of these, such as gap junctions, have been reported in the studies on Alzheimer’s disease [39].

Figure 2
figure 2

KEGG pathways and GO terms significantly regulated by miRNAs in Alzheimer’s disease. KEGG pathways significantly regulated by miRNAs are summarized in (A). Most KEGG pathways were shared by all four Alzheimer’s disease samples (64%, 32/50) (B). GO terms significantly regulated by miRNAs are summarized in (C). Four Alzheimer’s disease samples shared most of these biological processes (72%, 42/58) (D). (+) indicates the positive regulation while (−) indicates negative regulation.

Table 1 The top 5 pathways that were most affected by the miRNAs in the four Alzheimer’s disease samples
Figure 3
figure 3

KEGG pathways and GO terms significantly regulated by miRNAs after CdTe QD treatment. SOLiD sequencing based miRNA expression profile data of the CdTe QD treated cells were applied to our algorithm. Selected mRNAs indicated to be under significant repression by miRNA profiling were enriched using KEGG pathways (A) and GO terms (B). (+) indicates the positive regulation while (−) indicates negative regulation.

The GO provides the ontology of defined terms representing gene product properties. Among the three domains covered by the GO, biological processes, operations or sets of molecular events with either a defined beginning or end, can relevantly describe the functioning of integrated genes. The biological processes involved in Alzheimer’s disease were assessed based on the significantly regulated genes, using DAVID. Cellular processes were widely regulated in Alzheimer’s disease compared with the control group. Among these cellular processes, the nervous system development process was markedly altered in all four Alzheimer’s samples (P < 10-25). In this process, genes selected using our algorithm included presenilin (NM_000021), superoxide dismutase (NM_000384) and the oxytocin receptor (NM_000794), which have been reported to be involved in Alzheimer’s disease [4042]. Other processes implicated in Alzheimer’s disease are also presented in Figure 2A. Some processes were significantly regulated in only one sample. For example, dysregulation of cell death and programmed cell death were detected only in S3.

Most of the selected pathways and GO terms predicted by our algorithm were shared by all four Alzheimer’s disease samples (Figure 2B and D), indicating that this algorithm can be used to predict biological processes and pathways in cells or tissues based on their miRNA expression profiles.

Evaluation of nanomaterial cytotoxicity using miRNA expression profiles based on deep sequencing data

Due to their unique properties and diverse application in the life sciences, nanomaterials have attracted considerable interest recently [4345]. However, knowledge of the cellular effects of nanomaterials, such as cytotoxicity, is limited compared with the rapid increase in biological and medical applications [6]. The lack of reliable methods to assess the overall cellular effects of nanomaterials as opposed to consideration of conventional toxicity assays remains a crucial challenge. MiRNAs have been shown to repress gene expression at the post-transcriptional level and to participate in a wide range of cellular processes. Moreover, it is conceivable that miRNAs participate in the cytotoxic activity of nanomaterials, such as apoptosis-like cell death [46]. Combined with the higher stability of miRNA relative to that of mRNA, investigation of miRNA expression profiles represents a beneficial technique for elucidation of the biological effects and the biocompatibility of nanomaterials.

Previously we have reported that miRNAs may participate in the cytotoxicity of cadmium telluride quantum dots (CdTe QD) [6, 16]. The expression patterns of miRNAs were extensively affected after CdTe QD treatment, resulting in apoptosis-like cell death. SOLiD sequencing based miRNA expression profile data were applied to our algorithm. Selected mRNAs indicated to be under significant repression by miRNA profiling were enriched using KEGG pathways and GO terms. KEGG pathways significantly regulated after CdTe QD exposure are shown in Figure 3A. The top pathway on this list was “apoptosis” (P < 0.001). Figure 4A summarizes the factors involved in apoptosis pathways. According to this diagram, CdTe QD treatment induced cell apoptosis via the caspase-3 pathway. Other factors including Fas, IL1 or calcium-related pathways may participate in regulating capsase-3. Furthermore, this diagram also displays a cell self-protection mechanism against apoptosis mediated via the inhibitors of apoptosis (IAP) pathway (Figure 4A).

Figure 4
figure 4

Biological processes in GO terms significantly regulated by miRNAs after CdTe QD treatment.

The functions of the mRNAs selected as significantly regulated genes by miRNA profiling with GO analysis were then annotated. As shown in Figure 3B, at the biological process level 4, the top five processes induced by CdTe QDs were: negative regulation of growth, regulation of developmental growth, negative regulation of developmental growth, regulation of cell motion and negative regulation of developmental growth. Other processes included the regulation of cell death, regulation of programmed cell death, negative regulation of cell migration, apoptosis and endocytosis (Figure 3B).

In KEGG pathways and GO terms analyses, one of the main effects induced by CdTe QDs was apoptosis-like cell death and the apoptosis-related proteins were shown to be modulated. According to KEGG pathway analysis, apoptosis proteins were significantly regulated and GO terms analysis annotated CdTe QDs exposure may arrest cell growth and induce apoptosis, which is in accordance with other reports that CdTe QDs induce cell death in variety of cell types [45, 4749].

MiRNAs play a major role in cellular biological processes, including viability, proliferation, development and differentiation. Due to the stability of miRNAs during analysis and the development of sequencing methods, miRNA expression profiling identified as a reliable method for investigation of mRNAs, proteins and pathways involved in particular biological processes. In this study, we proposed a new method for validation of miRNA microarray or deep sequencing results for prediction of proteins and pathways under regulation. Our findings may contribute to early diagnosis of disease and assessment of the biological safety of new materials and drugs.

Conclusion

MiRNAs play a major role in cellular biological processes, including viability, proliferation, development and differentiation. Due to the stability of miRNAs during analysis and the development of sequencing methods, miRNA expression profiling identified as a reliable method for investigation of mRNAs, proteins and pathways involved in particular biological processes. In this study, we proposed a new method for validation of miRNA microarray or deep sequencing results for prediction of proteins and pathways under regulation. Using our method, the pathways and biological processes involved in Alzheimer’s disease were predicted with high correlation in four independent samples. Moreover, this method was successfully used for annotation of miRNA expression profiles from deep sequencing data for evaluation of CdTe nanomaterial cytotoxicity. As a result, apoptosis pathways were selected as the top pathways involved in CdTe nanoparticle treatment. Our findings may contribute to early diagnosis of disease and assessment of the biological safety of new materials and drugs.

Abbreviations

GO:

Gene ontology

3′:

UTR 3′untranslated region

CdTe:

QD cadmium telluride quantum dots

IAP:

Inhibitors of apoptosis

GEO:

Gene expression omnibus.

References

  1. Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004, 116 (2): 281-297. 10.1016/S0092-8674(04)00045-5.

    Article  CAS  PubMed  Google Scholar 

  2. Ambros V: The functions of animal microRNAs. Nature. 2004, 431 (7006): 350-355. 10.1038/nature02871.

    Article  CAS  PubMed  Google Scholar 

  3. Huntzinger E, Izaurralde E: Gene silencing by microRNAs: contributions of translational repression and mRNA decay. Nat Rev Genet. 2011, 12 (2): 99-110. 10.1038/nrg2936.

    Article  CAS  PubMed  Google Scholar 

  4. Hu W, Coller J: What comes first: translational repression or mRNA degradation? The deepening mystery of microRNA function. Cell Res. 2012, 22 (9): 1322-1324. 10.1038/cr.2012.80.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  5. Morita T, Mochizuki Y, Aiba H: Translational repression is sufficient for gene silencing by bacterial small noncoding RNAs in the absence of mRNA destruction. Proc Natl Acad Sci USA. 2006, 103 (13): 4858-4863. 10.1073/pnas.0509638103.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  6. Li S, Wang H, Qi Y, Tu J, Bai Y, Tian T, Huang N, Wang Y, Xiong F, Lu Z, et al: Assessment of nanomaterial cytotoxicity with SOLiD sequencing-based microRNA expression profiling. Biomaterials. 2011, 32 (34): 9021-9030. 10.1016/j.biomaterials.2011.08.033.

    Article  CAS  PubMed  Google Scholar 

  7. Git A, Dvinge H, Salmon-Divon M, Osborne M, Kutter C, Hadfield J, Bertone P, Caldas C: Systematic comparison of microarray profiling, real-time PCR, and next-generation sequencing technologies for measuring differential microRNA expression. RNA. 2010, 16 (5): 991-1006. 10.1261/rna.1947110.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  8. Li X, Jiang W, Li W, Lian B, Wang S, Liao M, Chen X, Wang Y, Lv Y, Wang S, et al: Dissection of human MiRNA regulatory influence to subpathway. Brief Bioinform. 2012, 13 (2): 175-186. 10.1093/bib/bbr043.

    Article  PubMed  Google Scholar 

  9. Wang Y-P, Li K-B: Correlation of expression profiles between microRNAs and mRNA targets using NCI-60 data. BMC Genomics. 2009, 10 (1): 218-10.1186/1471-2164-10-218.

    Article  PubMed Central  PubMed  Google Scholar 

  10. Ulitsky I, Laurent LC, Shamir R: Towards computational prediction of microRNA function and activity. Nucleic Acids Res. 2010, 38 (15): e160-10.1093/nar/gkq570.

    Article  PubMed Central  PubMed  Google Scholar 

  11. Corrada D, Viti F, Merelli I, Battaglia C, Milanesi L: myMIR: a genome-wide microRNA targets identification and annotation tool. Brief Bioinform. 2011, 12 (6): 588-600. 10.1093/bib/bbr062.

    Article  CAS  PubMed  Google Scholar 

  12. Sales G, Coppe A, Bicciato S, Bortoluzzi S, Romualdi C: Impact of probe annotation on the integration of miRNA–mRNA expression profiles for miRNA target detection. Nucleic Acids Res. 2010, 38 (7): e97-10.1093/nar/gkp1239.

    Article  PubMed Central  PubMed  Google Scholar 

  13. 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.

    Article  PubMed Central  PubMed  Google Scholar 

  14. Lewis BP, Shih I, Jones-Rhoades MW, Bartel DP, Burge CB: Prediction of mammalian MicroRNA targets. Cell. 2003, 115 (7): 787-798. 10.1016/S0092-8674(03)01018-3.

    Article  CAS  PubMed  Google Scholar 

  15. Rajewsky N: microRNA target predictions in animals. Nat Genet. 2006

    Google Scholar 

  16. Li S, Wang Y, Wang H, Bai Y, Liang G, Wang Y, Huang N, Xiao Z: MicroRNAs as participants in cytotoxicity of CdTe quantum dots in NIH/3 T3 cells. Biomaterials. 2011, 32 (15): 3807-3814. 10.1016/j.biomaterials.2011.01.074.

    Article  CAS  PubMed  Google Scholar 

  17. Gusev Y: Computational methods for analysis of cellular functions and pathways collectively targeted by differentially expressed microRNA. Methods. 2008, 44 (1): 61-72. 10.1016/j.ymeth.2007.10.005.

    Article  CAS  PubMed  Google Scholar 

  18. Hua Y-J, Tang Z-Y, Tu K, Zhu L, Li Y-X, Xie L, Xiao H-S: Identification and target prediction of miRNAs specifically expressed in rat neural tissue. BMC Genomics. 2009, 10 (1): 214-10.1186/1471-2164-10-214.

    Article  PubMed Central  PubMed  Google Scholar 

  19. Yan X, Chao T, Tu K, Zhang Y, Xie L, Gong Y, Yuan J, Qiang B, Peng X: Improving the prediction of human microRNA target genes by using ensemble algorithm. FEBS Lett. 2007, 581 (8): 1587-1593. 10.1016/j.febslet.2007.03.022.

    Article  CAS  PubMed  Google Scholar 

  20. Betel D, Koppal A, Agius P, Sander C, Leslie C: Comprehensive modeling of microRNA targets predicts functional non-conserved and non-canonical sites. Genome Biol. 2010, 11 (8): R90-10.1186/gb-2010-11-8-r90.

    Article  PubMed Central  PubMed  Google Scholar 

  21. Cao S, Chen S-J: Predicting kissing interactions in microRNA–target complex and assessment of microRNA activity. Nucleic Acids Res. 2012, 40 (10): 4681-4690. 10.1093/nar/gks052.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  22. Chi SW, Hannon GJ, Darnell RB: An alternative mode of microRNA target recognition. Nat Struct Mol Biol. 2012, 19 (3): 321-327. 10.1038/nsmb.2230.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  23. Nunez-Iglesias J, Liu C-C, Morgan TE, Finch CE, Zhou XJ: Joint genome-wide profiling of miRNA and mRNA expression in Alzheimer’s disease cortex reveals altered miRNA regulation. PloS one. 2010, 5 (2): e8898-10.1371/journal.pone.0008898.

    Article  PubMed Central  PubMed  Google Scholar 

  24. Krek A, Grun D, Poy MN, Wolf R, Rosenberg L, Epstein EJ, MacMenamin P, da Piedade I, Gunsalus KC, Stoffel M, et al: Combinatorial microRNA target predictions. Nat Genet. 2005, 37 (5): 495-500. 10.1038/ng1536.

    Article  CAS  PubMed  Google Scholar 

  25. Brennecke J, Stark A, Russell RB, Cohen SM: Principles of MicroRNA–target recognition. PLoS Biol. 2005, 3 (3): e85-10.1371/journal.pbio.0030085.

    Article  PubMed Central  PubMed  Google Scholar 

  26. Lai EC, Tam B, Rubin GM: Pervasive regulation of Drosophila Notch target genes by GY-box-, Brd-box-, and K-box-class microRNAs. Genes Dev. 2005, 19 (9): 1067-1080. 10.1101/gad.1291905.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  27. Farh KK-H, Grimson A, Jan C, Lewis BP, Johnston WK, Lim LP, Burge CB, Bartel DP: The widespread impact of mammalian MicroRNAs on mRNA repression and evolution. Science. 2005, 310 (5755): 1817-1821. 10.1126/science.1121158.

    Article  CAS  PubMed  Google Scholar 

  28. Sahoo S, Albrecht AA: Ranking of microRNA target prediction scores by Pareto front analysis. Comput Biol Chem. 2010, 34 (5–6): 284-292.

    Article  CAS  PubMed  Google Scholar 

  29. Friedman RC, Farh KK-H, Burge CB, Bartel DP: Most mammalian mRNAs are conserved targets of microRNAs. Genome Res. 2009, 19 (1): 92-105.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  30. Garcia DM, Baek D, Shin C, Bell GW, Grimson A, Bartel DP: Weak seed-pairing stability and high target-site abundance decrease the proficiency of lsy-6 and other microRNAs. Nat Struct Mol Biol. 2011, 18 (10): 1139-1146. 10.1038/nsmb.2115.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  31. De Ferrari GV, Inestrosa NC: Wnt signaling function in Alzheimer’s disease. Brain Res Rev. 2000, 33 (1): 1-12. 10.1016/S0165-0173(00)00021-7.

    Article  CAS  PubMed  Google Scholar 

  32. Luo J, Chen J, Deng Z-L, Luo X, Song W-X, Sharff KA, Tang N, Haydon RC, Luu HH, He T-C: Wnt signaling and human diseases: what are the therapeutic implications?. Lab Invest. 2007, 87 (2): 97-103. 10.1038/labinvest.3700509.

    Article  CAS  PubMed  Google Scholar 

  33. Kim EK, Choi E-J: Pathological roles of MAPK signaling pathways in human diseases. Biochimica et Biophysica Acta. 2010, 1802 (4): 396-405. 10.1016/j.bbadis.2009.12.009.

    Article  CAS  PubMed  Google Scholar 

  34. Dineley KT, Westerman M, Bui D, Bell K, Ashe KH, Sweatt JD: β-amyloid activates the mitogen-activated protein kinase cascade via hippocampal α7 nicotinic acetylcholine receptors:in vitro and in vivo mechanisms related to Alzheimer’s disease. J Neurosci. 2001, 21 (12): 4125-4133.

    CAS  PubMed  Google Scholar 

  35. Uchida Y, Ohshima T, Sasaki Y, Suzuki H, Yanai S, Yamashita N, Nakamura F, Takei K, Ihara Y, Mikoshiba K, et al: Semaphorin3A signalling is mediated via sequential Cdk5 and GSK3β phosphorylation of CRMP2: implication of common phosphorylating mechanism underlying axon guidance and Alzheimer’s disease. Genes Cells. 2005, 10 (2): 165-179. 10.1111/j.1365-2443.2005.00827.x.

    Article  CAS  PubMed  Google Scholar 

  36. Lesnick TG, Papapetropoulos S, Mash DC, Ffrench-Mullen J, Shehadeh L, de Andrade M, Henley JR, Rocca WA, Ahlskog JE, Maraganore DM: A genomic pathway approach to a complex disease: axon guidance and parkinson disease. PLoS Genet. 2007, 3 (6): e98-10.1371/journal.pgen.0030098.

    Article  PubMed Central  PubMed  Google Scholar 

  37. Keating DJ, Chen C, Pritchard MA: Alzheimer’s disease and endocytic dysfunction: clues from the down syndrome-related proteins, DSCR1 and ITSN1. Ageing Res Rev. 2006, 5 (4): 388-401. 10.1016/j.arr.2005.11.001.

    Article  CAS  PubMed  Google Scholar 

  38. Zhang C, Lambert MP, Bunch C, Barber K, Wade WS, Krafft GA, Klein WL: Focal adhesion kinase expressed by nerve cell lines shows increased tyrosine phosphorylation in response to Alzheimer’s a beta peptide. J Biol Chem. 1994, 269 (41): 25247-25250.

    CAS  PubMed  Google Scholar 

  39. Nakase T, Naus CCG: Gap junctions and neurological disorders of the central nervous system. Biochimica et Biophysica Acta. 2004, 1662 (1–2): 149-158.

    Article  CAS  PubMed  Google Scholar 

  40. Scheuner D, Eckman C, Jensen M, Song X, Citron M, Suzuki N, Bird TD, Hardy J, Hutton M, Kukull W, et al: Secreted amyloid [beta]-protein similar to that in the senile plaques of Alzheimer’s disease is increased in vivo by the presenilin 1 and 2 and APP mutations linked to familial Alzheimer’s disease. Nat Med. 1996, 2 (8): 864-870. 10.1038/nm0896-864.

    Article  CAS  PubMed  Google Scholar 

  41. Curtain CC, Ali F, Volitakis I, Cherny RA, Norton RS, Beyreuther K, Barrow CJ, Masters CL, Bush AI, Barnham KJ: Alzheimer’s disease amyloid-β binds copper and zinc to generate an allosterically ordered membrane-penetrating structure containing superoxide dismutase-like subunits. J Biol Chem. 2001, 276 (23): 20466-20473. 10.1074/jbc.M100175200.

    Article  CAS  PubMed  Google Scholar 

  42. Mazurek MF, Beal MF, Bird ED, Martin JB: Oxytocin in Alzheimer’s disease. Neurology. 1987, 37 (6): 1001-10.1212/WNL.37.6.1001.

    Article  CAS  PubMed  Google Scholar 

  43. Kim J-W, Galanzha EI, Shashkov EV, Moon H-M, Zharov VP: Golden carbon nanotubes as multimodal photoacoustic and photothermal high-contrast molecular agents. Nat Nano. 2009, 4 (10): 688-694. 10.1038/nnano.2009.231.

    Article  CAS  Google Scholar 

  44. Galanzha EI, Shashkov EV, Kelly T, Kim J-W, Yang L, Zharov VP: In vivo magnetic enrichment and multiplex photoacoustic detection of circulating tumour cells. Nat Nano. 2009, 4 (12): 855-860. 10.1038/nnano.2009.333.

    Article  CAS  Google Scholar 

  45. Medintz IL, Uyeda HT, Goldman ER, Mattoussi H: Quantum dot bioconjugates for imaging, labelling and sensing. Nat Mater. 2005, 4 (6): 435-446. 10.1038/nmat1390.

    Article  CAS  PubMed  Google Scholar 

  46. Lovrić J, Cho SJ, Winnik FM, Maysinger D: Unmodified cadmium telluride quantum dots induce reactive oxygen species formation leading to multiple organelle damage and cell death. Chem Biol. 2005, 12 (11): 1227-1234. 10.1016/j.chembiol.2005.09.008.

    Article  PubMed  Google Scholar 

  47. Cho SJ, Maysinger D, Jain M, Röder B, Hackbarth S, Winnik FM: Long-term exposure to CdTe quantum dots causes functional impairments in live cells. Langmuir. 2007, 23 (4): 1974-1980. 10.1021/la060093j.

    Article  CAS  PubMed  Google Scholar 

  48. Lovrić J, Bazzi HS, Cuie Y, Fortin GRA, Winnik FM, Maysinger D: Differences in subcellular distribution and toxicity of green and red emitting CdTe quantum dots. J Mol Med. 2005, 83 (5): 377-385. 10.1007/s00109-004-0629-x.

    Article  PubMed  Google Scholar 

  49. Khatchadourian A, Maysinger D: Lipid droplets: their role in nanoparticle-induced oxidative stress. Mol Pharm. 2009, 6 (4): 1125-1137. 10.1021/mp900098p.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgments

This work was supported by grants from the Major State Basic Research Development Program of China (973 Program) (No. 2013CB932902), the National Natural Science Foundation of China (NSFC) (No. 61071047), and the Natural Science Foundation of Jiangsu Province (No. BK2012122).

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Bo Sun or Zhong-Dang Xiao.

Additional information

Competing interests

We declare that there are no competing interests.

Authors’ contributions

BS, NPH, ZDX conceived the project. BS, ZDX designed and implement the algorithm and performed the analysis. FY, FHH contributed to the computational analysis. BS, ZDX wrote the paper. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Sun, B., Yang, F., Hu, FH. et al. Comprehensive annotation of microRNA expression profiles. BMC Genet 14, 120 (2013). https://doi.org/10.1186/1471-2156-14-120

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2156-14-120

Keywords