- Open Access
Profiles of the auditory epithelia related microRNA expression in neonatal and adult rats
European Journal of Medical Researchvolume 19, Article number: 48 (2014)
The impact of miRNA differential expression on the auditory epithelium stem cell development in postnatal rats is not clear. The present study was designed to analyze miRNA expression in the organ of Corti of neonatal and adult rats.
The cochleae of newborn (P0) and adult (P30) Sprague-Dawley rats were dissected in cold PBS to collect the sensory epithelia. Small RNAs were extracted using the mirVana RNA Isolation kit. Then, miRNA expression profiling was performed with RNAs from three newborns and three adult rats utilizing the TaqMan Array Rodent MicroRNA Panel.
Eighteen miRNAs were found be differentially expressed, 16 were unregulated in mature cochleae with the fold changes ranging from 17 to 600 folds. The expression levels of two miRNAs were reduced in the mature rat cochleae. GO analysis and signaling pathway analysis revealed the potential involvement of the miRNAs in the regulation of Wnt and TGF-β signaling pathways in hair cell development.
Our results provided novel insights into the functional significance of miRNAs in the basilar membrane cells development, and revealed the potential importance of miRNAs in the hair cell by regulation of Wnt and TGF-β signaling.
The inner ear, located in the temporal bone, contains the cochlea responsible for hearing and the vestibule responsible for balance. In the cochlea, the organ of Corti is a specialized structure that responds to fluid-borne vibrations. The organ of Corti is a complex organ that comprises a highly-ordered cellular mosaic of sensory hair cells (HCs) and non-sensory supporting cells (SCs). The generation of new HCs occurs throughout life in the auditory and vestibular sensory receptor of fishes and amphibians ,. In mammals, embryonic HC and SC proliferation within the sensory epithelium culminates between embryonic day 13 (E13) and E15, but stops after birth . It has been shown that acutely dissociated cells from the newborn rat or young rat organ of Corti can develop into otospheres consisting of 98% nestin cells when plated on a non-adherent substratum in the presence of either epidermal growth factors or fibroblast growth factors ,. Li and colleagues have shown that the adult utricular sensory epithelium contains cells that display the biological features of stem cells including self-renewal, sphere formation and capability of differentiating into hair cell-like cells . However, the replacement of lost hair cells does not occur spontaneously in pathological conditions such as age-related cochlear degeneration. A number of genes that affect various aspects of inner ear development have been identified and include transcription factors, morphogens, growth factors, receptors, and so on.
microRNAs (miRNAs) were discovered by Lee and colleague in Caenorhabditis elegans in 1993 as novel molecules that play an important role in gene expression regulation . miRNAs are small noncoding RNA molecules (approximately 22 nucleotides) that regulate posttranscriptional gene expression by relatively nonspecific binding to the 3’-untranslated region of mRNA . A single miRNA may regulate several genes because of the sequence similarity. It has been proposed that over one third of all protein-encoding genes are under translational control by miRNA . miRNAs are involved in a variety of cellular processes, including cellular differentiation, proliferation and apoptosis . miRNAs play an essential role in inner ear development . A recent study using conditionally knocked out Dicer only in the inner ear, SE hair and SCs after their normal differentiation from progenitor cells revealed the importance of miRNAs in inner ear development and function in vertebrates . Using an in silico prediction model that integrates miRNAs, mRNA and protein expression, Elkan-Miller and co-workers discovered the expression of 157 miRNAs in the inner ear sensory epithelium, with 53 miRNAs differentially expressed between the cochlea and the vestibule. Six miRNA families appear to be functionally important in the inner ear .
Zhang and colleagues  identified the miRNAs involved in degeneration of the organ of Corti during age-related hearing loss. They showed that 111 and 71 miRNAs exhibited differential expression in the C57 and CBA mice aged from postnatal day 21 to 16 months, respectively, and that downregulated miRNAs substantially outnumbered upregulated miRNAs during aging. However, comparisons of miRNA differential expression in the organ of Corti between newborn and adult rats, representing the early development of the inner ear sensory epithelium, have not yet been investigated. Therefore, in this study, we characterized the miRNAs expression profile in the auditory epithelia in both newborn and adult rats in order to examine the patterns and potential roles of miRNA differential expression in the early development of the inner ear sensory epithelium. The results showed that 16 differentially expressed miRNAs were identified. GO (Gene Ontology) term analysis revealed the importance of Wnt and transforming growth factor (TGF)-β signaling in the hair cell development. Understanding the miRNA and gene interaction network shed light on their roles on the development of normal and impaired hearing, and the mechanisms leading towards deafness.
Neonatal (P0) and adult (P30) Sprague-Dawley (SD) rats were approved by the Institutional Animal Care and Use Committees of Chinese PLA General Hospital.
The cochleae of new born (P0) and adult (P30) SD rats were dissected in cold PBS (10 mM Na2HPO4,,1.7 mM KH2PO4,137 mM NaCl, 2.7 mM KCl, pH 7.4) to collect the sensory epithelia. The collected tissues were stored in RNAlater (Ambion, Austin, TX, USA) until use. Small RNAs (<200 nucleotides) were extracted using the mirVana RNA Isolation kit (Ambion, Austin, TX, USA) according to the manufacture’s instruction. The quality and quantity of the RNA preparations were determined using a 2100 Agilent BioAnalyzer and a NanoDrop ND-1000 spectrophotometer (Thermo Scientific, Wilmington, DE, USA).
miRNA expression profiling was performed with RNAs from three newborns and three adult rats utilizing the TaqMan Array Rodent MicroRNA Panel (Applied Biosystems, Foster City, CA, USA) using 50 ng of RNA per port for a total of 400 ng. This array contains 365 miRNA targets as well as endogenous controls. Normalization was performed with the small nuclear RNAs (snRNAs) U44 and U48. These snRNAs are stably-expressed reference genes suitable for normalization of miRNAs. The qRT-PCR for the assessment of gene expression levels was performed using an ABI Prism 7900HT Sequence detection system (Applied Biosystems, Foster City, CA, USA).
Two class differential
We applied the random variance model (RVM) t-test to filter the differentially expressed miRNAs for the newborn and adult groups because the RVM t-test can raise degrees of freedom effectively in the cases of small samples. After the significant analysis and FDR analysis, we selected the differentially expressed genes according to the P-value threshold .
GO analysis was applied to analyze the main function of the differential expression genes according to the Gene Ontology which is the key functional classification of National Center for Biotechnology Information (NCBI) . Generally, Fisher’s exact test and χ2 test were used to classify the GO category, and the false discovery rate (FDR)  was calculated to correct the P-value, the smaller the FDR, the smaller the error in judging the P-value. The FDR was defined as:
where N k refers to the number of Fisher’s test P-values less than χ2 test P- values. We computed the P-values for the GOs of all the differentially expressed genes. This enrichment analysis provides a measure of the significance of the function: as the enrichment increases, the corresponding function is more specific, which enabled us to identify those GOs with more concrete function description in the experiment. Within the significant category, the enrichment Re was given by:
where n f is the number of differential genes within the particular category, n is the total number of genes within the same category, N f is the number of differential genes in the entire microarray, and N is the total number of genes in the microarray .
Pathway analysis was used to identify the significant pathway of the differentially expressed genes according to the Kyoto Encyclopedia of Genes and Genomes (KEGG), Biocarta and Reatome. Again, we used the Fisher’s exact test and χ2 test to select significant pathways. The threshold of significance was defined by P-value and FDR. The enrichment Re was calculated using the equation described above -. The relationship of the miRNAs and genes was counted by their differential expression values, and according to the interactions of miRNA and genes in Sanger MicroRNA database to build the MicroRNA-Gene-Network. The adjacency matrix of MicroRNA and genes A = [ai,j] is made by the attribute relationships among genes and MicroRNA, and [ai,j] represents the relation between the weight of gene i and MicroRNA j. In the MicroRNA-Gene-Network, circles represent genes and squares represent miRNAs, and their relationship is represented by one edge. The center of the network is represented by degree. Degree means the contribution made by one miRNA to the genes around or the contribution made by one gene to the miRNAs around. The key miRNAs and genes in the network always have the largest degrees ,.
The miRNA-GO-network is built according to the relationship between significant GOs and genes and the relationships among miRNAs and genes. The adjacency matrix of MicroRNA and genes A = [ai, j] is made by the attribute relationships among GOs and miRNAs, and [ai,j] represents the relation between the weight of GO i and MicroRNA j. In the MicroRNA-Gene-Network, circles represent genes and squares represent miRNAs, and their relationship is represented by one edge. The center of the network is represented by degree. Degree means the contribution one miRNA to the GOs around or the contribution one GO to the MicroRNAs around. The key miRNA and gene in the network always have the largest degrees.
miRNA expression profile analysis
To gain insights into the role of miRNAs that may be associated with the proliferative ability of cochlear cells during maturation of the cochlea, we examined the global expression pattern of mature miRNA using TaqMan® Rodent MicroRNA Arrays V2.0 (Applied Biosystems, Foster City, CA, USA). A total of 18 miRNAs exhibited expression changes in adult cochleae (Table 1) when compared with the newborn cochleae. Among these miRNAs, 16 were unregulated in mature cochleae with the fold changes ranging from 17 to 600 folds. The expression levels of two miRNAs (rno-miR-29c and rno-miR-29a) were reduced in the mature rat cochleae. This observation suggests miRNAs are involved in development process of cochleae.
Microarray-based GO analysis
According to the threshold of GOs significantly regulated by miRNAs, the P-value and FDR was < 0.001 and < 0.05, respectively. The high-enrichment GOs targeted by over-expressed miRNAs included negative regulation of epithelial cell differentiation, common-partner SMAD protein phosphorylation, mesenchymal-epithelial cell signaling, regulation of TGF-β2 (Figure 1). In contrast, significant GOs corresponding to downregulated miRNAs included protein heterotrimerization, negative regulation of phosphatidylinositol biosynthetic process and regulation of mitosis. Among these cellular processes, the maximum-enriched-GO relating to TGF-β2 and SMAD signaling suggested that they might have an important role in the proliferation potient of SE. Additionally, the miRNA-mRNA network analysis integrated these miRNAs and GOs by outlining the interactions of miRNA and GO-related genes (Figure 2). The miRNA-mRNA regulatory networks were established (Figure 3) and distinguished the putative target mRNAs between overexpressed and under-expressed miRNAs. Seven overexpressed miRNAs (miR-20a, miR-199a-5p, miR-199, miR-323, miR-301a, miR-301b and miR-130b) showed the most target mRNAs. The miRNAs including miR-301a, miR-301b and miR-130b regulated some important genes, including TGF-βR1, TGF-βR2, Smad2, Smad4 and Smad5 and, therefore, might be of great importance to the activation of the organ of Corti.
Signaling pathways regulated by differentially expressed miRNAs
Functional analysis of miRNAs by KEGG revealed that 19 signal transduction pathways were upregulated (Figure 4A) and 14 were downregulated (Figure 4B). The upregulated signaling pathways including Wnt, TGF-β and mitogen-activated protein kinase (MAPK) have been shown to participate in the activation of stem cells. A wide variety of cellular processes, including regulation of actin cytoskeleton, MAPK and gonadotropin releasing hormone (GnRH) signaling pathway, also featured the functions of significant signaling pathways.
miRNAs have become an area of intense study because of their involvement in human diseases. The lack of inner ear hair cell proliferation contributes to hearing loss in the aging population. In the present study, we compared the miRNA expression profiles between the newborn and adult cochlear sensory epithelia. Our results revealed that several miRNAs were differentially expressed in these two cochleae age groups. The difference in miRNA expression may contribute to the loss of proliferation of sensory cells in the organ of Corti in adult cochleae. These results provide novel insights into the functional significance of miRNAs in the basilar membrane cells development.
Given that embryonic HC and SC proliferation within the sensory epithelia culminates between embryonic day 13 (E13) and E15, but stops after birth in mammalian ears, the fundamental mechanisms behind the embryonic HC and SC proliferation potential are probably lost after birth. One of the possibilities is that coordinated and tightly controlled gene expression programs orchestrate the developmental process. miRNAs as key regulators might play important roles during this phenotypic transition, further adding another layer of complexity to the regulatory network for basilar membrane proliferation. Our miRNAs microarray data suggested that the expression profile of miRNAs in the rat inner ear appears to be well established by P0, consistent with the fact that early inner ear development and cell fate specification mostly occur embryonically . Our profiling data identified two distinct expression patterns of miRNAs between newborn and adult rat basilar membrane. Such differences appeared to be associated with basilar membrane proliferation. The miRNAs expression profiling identified 18 differentially expressed miRNAs with 16 miRNAs being increased and two miRNAs decreased in the mature rat compared to that of the newborn rat. miRNAs that increased mostly in the adult basilar membrane includes miR-296, mi-130b and miR-183 and that those that decreased mostly include miR-29c and miR-29a. MiR-296 has been demonstrated to modulate the pluripotency of embryonic stem cells (ESCs) by repressing the expression of Oct4, Sox2, and Nanog . In vertebrates, the expression domain of conserved miRNA-183 (miR-183) family members appears to be restricted to ciliated neurosensory epithelial cells and certain cranial and spinal ganglia ,. In zebrafish the miRNAs are detected in the eye, nasal epithelium, and sensory hair cells of the ear and neuromasts , and injection of miR-183 and miR-200 family members into zebrafish embryos have been demonstrated to impact development and affect neuromast migration . Additionally, expression of miR-183 family members in mouse eye and aural sensory hair cells of the ear has been previously demonstrated .
Compared with the differential expression pattern of miRNAs between newborn (P0) and adult rats (P30), Zhang and colleagues  reported different or even opposite patterns when they questioned which miRNAs are involved in age-related (from P21 to 16 months) degeneration of the organ of Corti, the auditory sensory epithelium that transduces mechanical stimuli to electrical activity in the inner ear. They showed that 111 and 71 miRNAs exhibited differential expression in the C57 and CBA mice, respectively, and that downregulated miRNAs substantially outnumbered upregulated miRNAs during aging. miRNAs that had approximately 2-fold upregulation included members of miR-29 family and miR-34 family and that were downregulated by about 2-fold were members of the miR-181 family and miR-183 family. The inconsistency between Zhang’s report and our study suggested that miRNA patterns in the organ of Corti change with aging and that miRNAs such as miR-183 and miR29 play different roles in the development of organ of Corti in newborn, younger and older animals.
The present GO analysis and signaling pathway analysis showed that the high-enrichment GOs targeted by over-expressed miRNAs in young adult (P30) compared with newborn rats (P0) included negative regulation of epithelial cell differentiation, common-partner SMAD protein phosphorylation, mesenchymal-epithelial cell signaling, regulation of TGF-β2 production. In contrast, significant GOs corresponding to downregulated miRNAs included protein heterotrimerization, negative regulation of phosphatidylinositol biosynthetic process and regulation of mitosis. Among these cellular processes, the maximum-enriched-GO relating to TGF-β2 and SMAD signaling suggested that they have an important role in the proliferation potient of SE. However, Zhang and colleagues  reported that miRNAs upregulated in age-related mice (from P21 to 16 months) are known regulators of pro-apoptotic pathways. In contrast, downregulated miRNAs are known to be important for proliferation and differentiation, respectively. The author concluded that the shift of miRNA expression favoring apoptosis occurred earlier than detectable hearing threshold elevation and hair cell loss. The authors suggested that changes in miRNA expression precede morphological and functional changes, and that upregulation of pro-apoptotic miRNAs and downregulation of miRNAs promoting proliferation and differentiation are both involved in age-related degeneration of the organ of Corti. The inconsistency in functions between Zhang’s report and our study can be explained by the different roles of miRNAs in the development of the organ of Corti in newborn, younger and older animals.
Establishment of primitive streak cells upon differentiation of ESCs depends on the presence of active Wnt and TGF-β/nodal/activin signaling, which recapitulates early events that led to germ-layer induction in the mammalian embryo . In our study, we found that miRNAs that inhibit Wnt and TGF-β signaling pathway were decreased in adult rats. Recently, Oshima and colleagues generated mechanosensitive sensory hair cell-like cells from embryonic and induced pluripotent stem cells by the combination of Wnt inhibitor Dkk1, selective inhibitor of Smad3 (SIS3) which interferes with TGF-β signaling, and insulin-like growth factor 1 (IGF-1) . Consistent with previous studies, our findings revealed the potential importance of miRNAs in the hair cell by regulation of Wnt and TGF-β signaling.
Our results provided novel insights into the functional significance of miRNAs in the basilar membrane cells development, and revealed the potential importance of miRNAs in the hair cell by regulation of Wnt and TGF-β signaling.
WG and SY conceived and designed the study. WG, YY, ZH performed the experiments and analyzed the data. XW provided data on the device and reviewed the manuscript. All authors read and approved the final manuscript.
embryonic stem cells
false discovery rate
sensory hair cells
insulin-like growth factor 1
Kyoto Encyclopedia of Genes and Genomes
mitogen-activated protein kinase
gonadotropin releasing hormone
National Center for Biotechnology Information
polymerase chain reaction
non-sensory supporting cells
small nuclear RNAs
transforming growth factor
Corwin JT: Perpetual production of hair cells and maturational changes in hair cell ultrastructure accompany postembryonic growth in an amphibian ear. Proc Natl Acad Sci U S A 1985, 82: 3911–3915. 10.1073/pnas.82.11.3911
Corwin JT, Cotanche DA: Regeneration of sensory hair cells after acoustic trauma. Science 1988, 240: 1772–1774. 10.1126/science.3381100
Ruben RJ: Development of the inner ear of the mouse: a radioautographic study of terminal mitoses. Acta Otolaryngol 1967, Suppl- 220: 1–44.
Malgrange B, Belachew S, Thiry M, Nguyen L, Rogister B, Alvarez ML, Rigo JM, Van De Water TR, Moonen G, Lefebvre PP: Proliferative generation of mammalian auditory hair cells in culture. Mech Dev 2002, 112: 79–88. 10.1016/S0925-4773(01)00642-6
Lou X, Zhang Y, Yuan C: Multipotent stem cells from the young rat inner ear. Neurosci Lett 2007, 416: 28–33. 10.1016/j.neulet.2006.12.061
Li H, Liu H, Heller S: Pluripotent stem cells from the adult mouse inner ear. Nat Med 2003, 9: 1293–1299. 10.1038/nm925
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: 843–854. 10.1016/0092-8674(93)90529-Y
Scaria V, Hariharan M, Pillai B, Maiti S, Brahmachari SK: Host-virus genome interactions: macro roles for microRNAs. Cell Microbiol 2007, 9: 2784–2794. 10.1111/j.1462-5822.2007.01050.x
Friedman RC, Farh KK, Burge CB, Bartel DP: Most mammalian mRNAs are conserved targets of microRNAs. Genome Res 2009, 19: 92–105. 10.1101/gr.082701.108
Croce CM, Calin GA: miRNAs, cancer, and stem cell division. Cell 2005, 122: 6–7. 10.1016/j.cell.2005.06.036
Li H, Fekete DM: MicroRNAs in hair cell development and deafness. Curr Opin Otolaryngol Head Neck Surg 2010, 18(5):459–465. 10.1097/MOO.0b013e32833e0601
Friedman LM, Dror AA, Mor E, Tenne T, Toren G, Satoh T, Biesemeier DJ, Shomron N, Fekete DM, Hornstein E, Avraham KB: MicroRNAs are essential for development and function of inner ear hair cells in vertebrates. Proc Natl Acad Sci U S A 2009, 106: 7915–7920. 10.1073/pnas.0812446106
Elkan-Miller T, Ulitsky I, Hertzano R, Rudnicki A, Dror AA, Lenz DR, Elkon R, Irmler M, Beckers J, Shamir R, Avraham KB: Integration of transcriptomics, proteomics, and microRNA analyses reveals novel microRNA regulation of targets in the mammalian inner ear. PLoS One 2011, 6: e18195. 10.1371/journal.pone.0018195
Zhang Q, Liu H, McGee J, Walsh EJ, Soukup GA, He DZ: Identifying microRNAs involved in degeneration of the organ of Corti during age-related hearing loss. PLoS One 2013, 8: e62786. 10.1371/journal.pone.0062786
Wright GW, Simon RM: A random variance model for detection of differential gene expression in small microarray experiments. Bioinformatics 2003, 19: 2448–2455. 10.1093/bioinformatics/btg345
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G: Gene ontology: tool for the unification of biology. The gene ontology consortium. Nat Genet 2003, 25: 25–29. 10.1038/75556
Dupuy D, Bertin N, Hidalgo CA, Venkatesan K, Tu D, Lee D, Rosenberg J, Svrzikapa N, Blanc A, Carnec A, Carvunis AR, Pulak R, Shingles J, Reece-Hoyes J, Hunt-Newbury R, Viveiros R, Mohler WA, Tasan M, Roth FP, Le Peuch C, Hope IA, Johnsen R, Moerman DG, Barabási AL, Baillie D, Vidal M: Genome-scale analysis of in vivo spatiotemporal promoter activity in Caenorhabditis elegans . Nat Biotechnol 2007, 25: 663–668. 10.1038/nbt1305
Schlitt T, Palin K, Rung J, Dietmann S, Lappe M, Ukkonen E, Brazma A: From gene networks to gene function. Genome Res 2003, 13: 2568–2576. 10.1101/gr.1111403
Kanehisa M, Goto S, Sato Y, Furumichi M, Tanabe M: KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res 2012, 40: D109-D114. 10.1093/nar/gkr988
Yi M, Horton JD, Cohen JC, Hobbs HH, Stephens RM: Whole pathway scope: a comprehensive pathway-based analysis tool for high-throughput data. BMC Bioinformatics 2006, 7: 30. 10.1186/1471-2105-7-30
Draghici S, Khatri P, Tarca AL, Amin K, Done A, Voichita C, Georgescu C, Romero R: A systems biology approach for pathway level analysis. Genome Res 2007, 17: 1537–1545. 10.1101/gr.6202607
Joung JG, Hwang KB, Nam JW, Kim SJ, Zhang BT: Discovery of microRNA-mRNA modules via population-based probabilistic learning. Bioinformatics 2007, 23: 1141–1147. 10.1093/bioinformatics/btm045
Shalgi R, Lieber D, Oren M, Pilpel Y: Global and local architecture of the mammalian microRNA-transcription factor regulatory network. PLoS Comput Biol 2007, 3: e131. 10.1371/journal.pcbi.0030131
Barald KF, Kelley MW: From placode to polarization: new tunes in inner ear development. Development 2004, 131: 4119–4130. 10.1242/dev.01339
Houbaviy HB, Murray MF, Sharp PA: Embryonic stem cell-specific MicroRNAs. Dev Cell 2003, 5: 351–358. 10.1016/S1534-5807(03)00227-2
Wienholds E, Kloosterman WP, Miska E, Alvarez-Saavedra E, Berezikov E, de Bruijn E, Horvitz HR, Kauppinen S, Plasterk RH: MicroRNA expression in zebrafish embryonic development. Science 2005, 309: 310–311. 10.1126/science.1114519
Kloosterman WP, Wienholds E, De Bruijn E, Kauppinen S, Plasterk RH: In situ detection of miRNAs in animal embryos using LNA-modified oligonucleotide probes. Nat Methods 2006, 3: 27–29. 10.1038/nmeth843
Ason B, Darnell DK, Wittbrodt B, Berezikov E, Kloosterman WP, Wittbrodt J, Antin PB, Plasterk RH: Differences in vertebrate microRNA expression. Proc Natl Acad Sci U S A 2006, 103: 14385–14389. 10.1073/pnas.0603529103
Weston MD, Pierce ML, Rocha-Sanchez S, Beisel KW, Soukup GA: MicroRNA gene expression in the mouse inner ear. Brain Res 2006, 1111: 95–104. 10.1016/j.brainres.2006.07.006
Gadue P, Huber TL, Paddison PJ, Keller GM: Wnt and TGF-beta signaling are required for the induction of an in vitro model of primitive streak formation using embryonic stem cells. Proc Natl Acad Sci U S A 2006, 103: 16806–16811. 10.1073/pnas.0603916103
Oshima K, Shin K, Diensthuber M, Peng AW, Ricci AJ, Heller S: Mechanosensitive hair cell-like cells from embryonic and induced pluripotent stem cells. Cell 2010, 141: 704–716. 10.1016/j.cell.2010.03.035
This work was supported by grants from the Major State Basic Research Development Program of China (973 Program) (#2011CBA01000), the National Basic Research Program of China (973 Program) (#2012CB967900), the National Natural Science Foundation of China (NSFC) (#81000483), and the New Star Project of Science and Technology of Beijing (2010B083) to ZHH.
The authors declare that they have no competing interests.