Gene Expression Analysis of Forskolin Treated Basilar Papillae Identifies MicroRNA181a as a Mediator of Proliferation

Abstract

Background

Auditory hair cells spontaneously regenerate following injury in birds but not mammals. A better understanding of the molecular events underlying hair cell regeneration in birds may allow for identification and eventually manipulation of relevant pathways in mammals to stimulate regeneration and restore hearing in deaf patients.

Methodology/Principal Findings

Gene expression was profiled in forskolin treated (i.e., proliferating) and quiescent control auditory epithelia of post-hatch chicks using an Affymetrix whole-genome chicken array after 24 (n = 6), 48 (n = 6), and 72 (n = 12) hours in culture. In the forskolin-treated epithelia there was significant (p<0.05; >two-fold change) upregulation of many genes thought to be relevant to cell cycle control and inner ear development. Gene set enrichment analysis was performed on the data and identified myriad microRNAs that are likely to be upregulated in the regenerating tissue, including microRNA181a (miR181a), which is known to mediate proliferation in other systems. Functional experiments showed that miR181a overexpression is sufficient to stimulate proliferation within the basilar papilla, as assayed by BrdU incorporation. Further, some of the newly produced cells express the early hair cell marker myosin VI, suggesting that miR181a transfection can result in the production of new hair cells.

Conclusions/Significance

These studies have identified a single microRNA, miR181a, that can cause proliferation in the chicken auditory epithelium with production of new hair cells.

About the Authors:

Corey S. Frucht

Affiliations Medical Scientist Training Program, Yale School of Medicine, New Haven, Connecticut, United States of America, Interdepartmental Neuroscience Program, Yale University, New Haven, Connecticut, United States of America

Mohamed Uduman

Affiliations Interdepartmental Program in Computation Biology and Bioinformatics, Yale University, New Haven, Connecticut, United States of America, Department of Pathology, Yale School of Medicine, New Haven, Connecticut, United States of America

Jamie L. Duke

Affiliations Interdepartmental Program in Computation Biology and Bioinformatics, Yale University, New Haven, Connecticut, United States of America, Department of Pathology, Yale School of Medicine, New Haven, Connecticut, United States of America

Steven H. Kleinstein

Affiliations Interdepartmental Program in Computation Biology and Bioinformatics, Yale University, New Haven, Connecticut, United States of America, Department of Pathology, Yale School of Medicine, New Haven, Connecticut, United States of America

Joseph Santos-Sacchi

Affiliations Division of Otolaryngology, Department of Surgery, Yale School of Medicine, New Haven, Connecticut, United States of America, Department of Cellular and Molecular Physiology, Yale School of Medicine, New Haven, Connecticut, United States of America, Department of Neurobiology, Yale School of Medicine, New Haven, Connecticut, United States of America

Dhasakumar S. Navaratnam

* E-mail: [email protected]

Affiliations Department of Neurobiology, Yale School of Medicine, New Haven, Connecticut, United States of America, Department of Neurology, Yale School of Medicine, New Haven, Connecticut, United States of America

Introduction

Sensorineural hearing loss represents a major public health concern. Approximately 300 million people worldwide have moderate to profound hearing loss in both ears [1]. Loss of inner ear hair cells, which serve to transduce sound into neural impulses, is responsible for the majority of hearing loss. In humans, and other mammals, loss of hair cells is permanent since these organisms have no capacity for hair cell regeneration. In contrast, other non-mammalian vertebrates such as birds, reptiles, amphibian and fish are able to replace lost hair cells.

The basilar papilla, the avian auditory epithelium, is able to regenerate hair cells in response to hair cell loss (reviewed in [2]). However, the auditory epithelium shows no mitotic activity normally, a feature that is reminiscent of the mammalian auditory epithelium, and contrasts with the chick vestibular epithelium which shows continuous mitotic activity [3]. Thus, the avian auditory epithelium can be viewed as an intermediary in the evolution from the fish to the mammal, and we reason that the study of this epithelium will provide insight into why the mammalian auditory epithelium shows mitotic quiescence at rest (similar to the avian auditory epithelium), but is unable to proliferate in response to damage (in contrast to the avian auditory epithelium).

The basilar papilla is comprised of both sensory transducing hair cells and supporting cells. Following injury, it is the supporting cells which give rise to new hair cells [4], [5]. For example, exposure of birds to intense noise causes some supporting cells to leave growth-arrest, re-enter the cell cycle and ultimately differentiate into hair cells [6], [7], [8], [9]. New hair cells are first seen 4–5 days after the onset of exposure to an intense sound [6], [10] and undergo maturation so that by 20–28 days after stimulus onset they are virtually indistinguishable from unaffected cells [11]. In addition, some new hair cells arise from direct differentiation of supporting cells without an intervening mitotic step [12], [13], [14], [15], [16], [17], [18], [19]. After acoustic or ototoxic insult, birds initially have increased hearing thresholds, which eventually return nearly to baseline confirming that newly produced hair cells are functional [20]. It is believed that this recovery of function results from both regeneration of new hair cells as well as repair of those that have survived [21].

Though the intracellular pathways required for hair cell regeneration have not yet been fully elucidated, various pathways and signaling cascades have been implicated in this process. For example, it has been shown that treatment of the chick basilar papilla with forskolin, a potent adenylate cyclase activator that increases intracellular cAMP levels, causes a robust and widespread proliferation of supporting cells, leading to the production of new hair cells [22]. This effect is first seen after 72 hours in culture, occurs without upregulation of markers of apoptosis, and is significantly blocked by protein kinase A inhibitors. It therefore appears as though activation of this pathway can stimulate growth of new hair cells with limited cell injury. Forskolin also appears to have a mitogenic effect in the mammalian vestibular system, in which brief treatment with this compound results in an increase an supporting cell S-phase entry [23]. Moreover, the same study found that this effect is blocked by brief treatment with forskolin in the presence of monensin or bafilmycin, which inhibit recycling of membrane receptors. Increased cAMP levels may therefore lead to proliferation by causing upregulation of tyrosine kinase growth factor receptors. With age-related hearing loss in humans there is a substantial loss of hair cells within the organ of Corti [24], so it is likely that a robust proliferative process that produces new hair cells is what will ultimately lead to maximal recovery of function in hearing impaired patients. In the present study, microarrays were used to analyze gene expression in forskolin treated basilar papillae given the robustness of this effect.

In other systems, microarray analysis of differing phenotypes has revealed complex differences in mRNA expression. These complex differences are likely brought about by the combined activity of transcription factors and microRNAs (miRNAs) that regulate expression of multiple genes and specific genetic programs. In fact, a study of transcription factor expression in regenerating versus non-regenerating chicken basilar papillae revealed a number of transcription factors that were differentially expressed [25].

MiRNAs are a species of small regulatory RNA that bind to complementary sequences in mRNA and decrease gene expression by increasing degradation and blocking translation (reviewed in [26]). Though each miRNA binds a specific sequence, these molecules are promiscuous in that their complementary sequences are found in a number of different mRNAs [27]. As a result, miRNA expression may represent a way in which cells can rapidly and in a coordinated fashion up- or downregulate a large number of genes. Some miRNAs are known to be important for inner ear development [28], [29], [30], [31] and may play a role in hair cell regeneration, where there is presumably some reiteration of developmental events. However there have been few published studies on the role of miRNAs in auditory hair cell regeneration, specifically. One study revealed that the let-7 family of miRNA are downregulated in the newt inner ear after hair cell injury [32]. This particular study underscores the importance of examining more closely the role of miRNA in hair cell regeneration.

In order to gain insight into the molecular events which underlie auditory hair cell regeneration in chicken, we compared gene expression in forskolin treated (i.e., proliferating) and control (i.e., quiescent) auditory epithelia using Affymetrix whole genome chicken microarrays. The microarray data were validated by quantitative PCR. Gene set enrichment analyses were performed to identify miRNAs and transcription factors that may be important for hair cell regeneration. Many sets of genes representing predicted targets of specific miRNAs were found to be significantly enriched among genes downregulated in the forskolin-treated sensory epithelia after 72 hours, suggesting that these miRNAs were likely upregulated in the proliferating tissue. One of these microRNAs, miR181a, was selected for functional experiments given its role in proliferation and suppression of the cell cycle inhibitor p27 in human myeloid leukemia cells [33]. P27 is thought to represent a key barrier to hair cell regeneration in the mammalian inner ear [34], [35]. Basilar papillae transfected with a miR181a precursor had increased numbers of BrdU positive cells than epithelia that were transfected with a non-targeting miRNA. Some cells labeling for BrdU also expressed the early hair cell marker myosin VI, indicating that overexpression of miR181a is capable of producing proliferation with production of new hair cells.

Materials and Methods

Basilar papilla explant cultures

Animals were treated in accordance with policies set forth by the Yale Institutional Animal Care and Use Committee (protocol number 2007-10439). Under sterile conditions, cochlear ducts containing the basilar papillae were carefully dissected out of 0-day-old chicks and then individually cultured in DMEM with 10% fetal bovine serum with or without forskolin (final concentration 100 µM, delivered in 1% DMSO) for either 24, 48, or 72 hours at 37°C with 5% CO2. Control samples received DMSO at 1% as a vehicle control. At the end of 72 hours, the tegmentum vasculosum was dissected off to expose the auditory epithelium, which was delicately freed from the underlying cartilaginous plates. All explants were kept in culture for 72 hours because new hair cells are first seen in basilar papillae treated with forskolin after ∼72 hours of exposure to the small molecule [22]. Therefore, at even earlier time points in culture the molecular events that underlie hair cell proliferation are well underway. Each sample was comprised of three auditory epithelia from three different chicks and were put in 100 µL of DMEM and then immediately frozen at −80°C until RNA isolation could be performed. There were a total of 24 samples in this experiment: three 24-hour forskolin, three 24-hour control, three 48-hour forskolin, three 48-hour control, six 72-hour forskolin and six 72-hour control.

RNA isolation

Immediately upon removal of the frozen epithelia samples from the −80°C freezer, the lysis buffer from the RNAqueous® Kit (Ambion, Austin, TX) was added. Given the small amount of tissue and relatively large surface area to size ratio, no measures to actively disrupt or homogenize the tissue in each sample were necessary. Total RNA isolation then proceeded as directed by the RNAqueous® Kit manual. Sample quality was confirmed by gel analysis.

Microarray experiments

Double-stranded cDNA and biotin-labeled cRNA were synthesized from 1 to 5 µg of total RNA using a two-cycle target labeling kit [36]. Biotin-labeled cRNA was purified using the GeneChip Cleanup Module prior to fragmenting to a size of 35–200 bases by incubating at 94°C for 35 minutes in fragmentation buffer (40 mM Tris-acetate, pH 8.1, 100 mM potassium acetate, 30 mM magnesium acetate). Hybridization of the samples with Affymetrix chicken arrays was then performed at the Yale University Keck Facility according to the manufacturer's protocol [37]. All microarray data is MIAME compliant and has been uploaded to MIAMExpress (ArrayExpress accession: E-MEXP-2642, username = Reviewer_E-MEXP-2642, password = 1269882335242).

Differential gene expression analysis

The microarray gene expression analysis was performed using GeneSpring GX 9.0 software (Agilent Technologies, Santa Clara, CA). The data were GC-RMA normalized without baseline transformation, and then filtered by expression (20%–100% in at least one of the six or twelve samples at each time point) to exclude genes expressed only at very low levels in all samples. Median expression values were used for redundant probe sets. Genes differentially expressed were identified by using a fold-change cutoff of 2 and then performing an unpaired t-test after Benjamini-Hotchberg correction for multiple comparisons (p<0.05). Heat maps were generated using R version 2.10.1 (www.r-project.org) using the HeatPlus add-on in the Bioconductor package.

Quantitative PCR validation of microarray gene expression data

In order to validate the microarray data, twelve genes were selected to be verified by qPCR. Seven of the genes selected were upregulated and five were downregulated in the forskolin condition. The upregulated genes validated by qPCR were highly expressed in at least one condition and chosen to span a range in fold-change differences between the two experimental conditions at 72 hours, from 2.0 to 64.3. The 12 genes that were validated are BECN1, CCNI, CDKN1B (p27), CDKN2B, DPM1, EME1, FBXO8, ITGA4, NLGN1, OCM, PIGW, SNX1. Additionally, expression of the 18S ribosomal subunit was also assessed to allow for normalization of total RNA levels between samples using the 2−ΔΔCT method [38]. Primers were designed using MacVector (MacVector, Cary, NC) software and Primer 3 [39] and were validated by melting curve analysis. 5 µg of RNA was isolated from each sample as described above and was used to create first-strand cDNA using random hexamers and reverse transcriptase (Superscript II; Invitrogen, Carslbad, CA), according to the manufacturer's instructions. Amplification was performed using the SYBR Green Supermix kit (Bio-Rad, Hercules, CA) on an iCycler system (Bio-Rad, Hercules, CA). Three replicates were performed and the data averaged for each cDNA sample and primer pair combination. Those fold-changes with a 95% confidence interval excluding one were considered significant.

Gene sets used for enrichment analyses

Gene set enrichment analysis (GSEA) was performed to analyze the pattern of differential gene expression between forskolin treated and control basilar papillae using GSEA version 2 software [40]. Three different gene set packages were used, two of which were downloaded directly from the Broad institute website (www.broad.mit.edu/gsea).

To identify transcription factors which may be important for hair cell regeneration, a GSEA was performed using a gene set package called “c3.all.v2.5.symbols,” which was obtained from the BROAD Institute website. This package contains a series of gene sets defined by the presence of transcription factor motifs and predicted miRNA binding sites. The genes in each set of this package therefore share a cis-regulatory motif that is conserved across human, mouse, rat, and dog genomes and will be referred to hereafter as the “mammalian gene sets.” The motifs used come from Xie et al. (2005) and the TRANSFAC database and include sets of genes sharing a 3′-UTR miRNA binding motif [41]. Even though these gene sets are not based on the chicken genome, we reasoned that these results would be of interest given the known conservation of miRNA targets among mammals [42]. We further surmised that if the conservation of miRNA targets is poorly conserved between mammals and non-mammalian vertebrates then an enrichment analysis might produce no statistically significant results. Any significant findings, however, would warrant further examination.

In order to validate the GSEA results from the analyses performed using the mammalian gene sets, GSEA was also run using a gene set package that was generated by scanning the portions of the chicken genome that are conserved with frog (Xenopus tropicalis) to identify gene sets sharing particular transcription factor motifs defined in the TRANSFAC database. These gene sets will hereafter be referred to as the “chicken/frog conservation gene sets.” Only conserved regions were used to reduce the number of false positives, as there is evidence to suggest that transcription factor recognition sites tend to be enriched in conserved portions of the genome [43]. However, recent genome-wide studies have demonstrated that many transcription factor binding sites are not conserved [44], [45] and are also found outside of the proximal promoter region [44], [46]. Thus, our approach should be considered to be conservative in the sense that some true binding sites may be overlooked, but those that are identified have a higher likelihood of being functional. The frog genome in particular was used in light of this organism's close evolutionary relationship to chickens. The use of conservation between more distantly related organisms, such as chicken and mammals, would be very strict and would result in the omission of many potentially significant positive results. To create these gene sets we downloaded the transcription start site for all chicken RefSeq genes, defined by the May 2008 RefGene table [47] using the UCSC Genome Bioinformatics site (genome.ucsc.edu). The region 2 kilobases around each transcription start site was identified within the May 2006 genome-wide multiple alignment of 6 vertebrate species to the chicken genome [48], also available through UCSC Genome Bioinformatics site. In order to identify putative transcription factor binding sites, the chicken sequences, along with aligned regions from frog, were analyzed using the TRANSFAC MATCH algorithm with a cutoff chosen to minimize the sum of false positives and false negatives [49]. The analysis was performed for all vertebrate transcription factor matrices in the 2009.1 release of TRANSFAC [50], and putative binding sites were considered to be evolutionarily conserved if matches were also found at the aligned positions in both chicken and frog sequences and had no gaps present in the multiple alignment of these species.

To assess the functional relevance of the pattern of differentially expressed genes across a variety of organisms and systems, GSEA was performed using a gene set package called “c2.all.v2.5.symbols” which was also downloaded from the BROAD Institute website. This package defines curated sets of genes based on specific experimental findings from experiments on human and model animal tissue, which includes canonical pathways and chemical/genetic perturbations based gene sets. These gene sets will hereafter be referred to as the “curated gene sets.” The relevant reference is given on the BROAD website for each curated gene set for assessment of relevance.

Gene set enrichment analysis

The GeneSpring processed data from the 38,535 original probes was collapsed into 13,159 genes based on gene symbols. Some genes were assigned to one or multiple gene sets downloaded with the GSEA package. There are a total of 837 gene sets in the mammalian gene sets package, 50 of which were excluded by gene set size criteria (15–50), leaving 783 to be included in our analysis. 359 of the 566 chicken/frog conservation gene sets met this size threshold criterion, as did 1,197 out of the 1,892 curated gene sets (Table 1).

[Figure omitted. See PDF.]

Table 1. Gene set enrichment analysis results.

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

Genes were first ranked by their signal-to-noise ratio between the forskolin and control groups. In other words, the genes most consistently overexpressed in the forskolin samples were ranked at one end of the list and those most underexpressed in this group were ranked at the other end. GSEA asks whether the genes that make up a particular gene set are randomly distributed throughout the ranked list of genes, or whether they have a tendency to cluster at the top or bottom of the list. This question was answered by calculating for every gene set a running sum statistic, the maximum of which is referred to as the enrichment score (ES). Each ES is then normalized to the size of its corresponding gene set to produce a normalized enrichment score (NES). Statistical significance of each NES is determined by comparing that NES to the distribution of ESs generated by randomly permutating the genotype class labels of the data set. Therefore, the null hypothesis of a GSEA is that the distribution of a gene set throughout the ranked list of genes is random with regards to the two treatment conditions being compared. The alternative hypothesis is that location of genes in a set is associated with the treatment conditions. Both the p-value and false discovery rate (FDR; q-value) were calculated for each set. Gene sets with a FDR less than 0.25 and a p-value less than 0.05 were considered differentially expressed. The GSEA parameters used were as follows: metric = signal to noise; permutation number = 1,000; gene size minimum = 15; gene size maximum = 500; enrichment statistic = classic; permutation type = phenotype. One GSEA was performed using the 24-hour control data and 72-hour control data to negatively control for gene expression changes occurring with culture duration, as well as false positives not reflective of the treatment condition.

Pre-microRNA181a and anti-microRNA181a transfection

For functional miR181a experiments, basilar papilla explants from 0-day-old chicks were cultured with or without forskolin for 72 hours as described above. Basilar papillae from both conditions were also transfected with either a miR181a precursor (Ambion, Austin, TX) or a non-targeting miRNA which served as a negative control (Ambion, Austin, TX) at a final concentration of 100 nM. Transfections were performed for six hours in serum-free medium which was eventually replaced with regular medium. BrdU was present at a concentration of 0.01% for the entire culture duration, except for the six-hour transfection. Some basilar papillae were also cultured for 24 hours with forskolin as described above and then transfected with anti-miR181a (Ambion, Austin, TX) to suppress endogenous miR181a. Transfections were performed using the lipid-based X-tremeGENE SiRNA Transfection Reagent (Roche, Indianapolis, IN) per the manufacturer's instructions.

Antibody labeling

After 72 hours in culture basilar papillae were fixed in 4% PFA in PBS for 30 minutes. All fixative was removed by three 5-minute washes in PBS. The tissue was blocked and permeabilized using a solution of PBS with FBS (10%) and Triton-X (0.1%) for one hour. Each basilar papilla was then incubated with mouse anti-BrdU antibodies (1∶40, BD, Franklin Lakes, NJ) in PBS for one hour. After washing with PBS, alkaline phosphatase-conjugated goat anti-mouse IgG (1∶400, Santa Cruz, Santa Cruz, CA) in PBS with Triton-X (0.1%) was added for one hour. After washing in PBS, alkaline phosphatase substrate was produced using the NBT/BCIP Reagent Kit (Invitrogen, Carlsbad, CA) per the manufacturer's instructions and then added to the tissue for approximately five minutes. All substrate was then thoroughly washed with PBS.

Statistics

Basilar papillae immunohistochemically labeled for BrdU were viewed using brightfield microscopy and were imaged using a digital camera attached to a Zeiss Stemi SV 11 stereoscope. The borders of the sensory epithelium were easily viewed both by direct microscopy and in the captured images. The intraepithelial location of BrdU positive nuclei was confirmed through focal plane adjustments. Only those nuclei which were clearly within the epithelium were counted. Unpaired Student's t-tests were performed to determine the statistical significance of the comparisons discussed in the text using a cut-off of p<0.05.

Co-labeling experiments

To determine whether any newly produced BrdU positive cells express the early hair cell marker myosin VI, some basilar papillae were transfected with pre-miR181a as described above and cultured for 72 hours. At the end of the culture duration, the tegmentum vasculosum and tectorial membranes were removed prior to fixation in 4% PFA for 30 minutes. Basilar papillae were cryoprotected overnight in 30% sucrose and then embedded in PBS with 1% low melting point agarose and 18% sucrose. The tissue was cryosectioned at 10 µm and then blocked and permeabilized as described above. Sections were labeled for both myosin VI and BrdU as previously described [12]. Sections were then incubated with rabbit anti-myosin VI (1∶350, Proteus, Ramona, CA) in blocking solution for one hour. After washing, sections were then incubated in Alexa Fluor 546 conjugated anti-rabbit secondary antibodies (1∶1,000, Invitrogen, Carlsbad, CA) for one hour. After washing, basilar papilla sections were post-fixed in 4% PFA for one hour. The cryosections were then washed in PBS before 20 minute 2M HCl treatment to expose BrdU. After thorough washing, sections were incubated with mouse anti-BrdU antibodies (1∶40, BD, Franklin Lakes, NJ) for one hour. After 3 additional PBS washes, sections were incubated with Alexa Fluor 647 conjugated anti-mouse secondary antibodies (Cell Signaling Technology, Danvers, MA) for one hour. After thorough washing, sections were mounted on glass slides in Vectashield fluorescent mounting medium (Vector Laboratories, Burlingame, CA). Images were captured using a Zeiss LSM 510 confocal microscope.

Results and Discussion

Differential expression of genes in the forskolin and control groups increases as a function of time

We report here a systematic comparison of genome-wide expression between proliferating and non-proliferating chick auditory epithelia using Affymetrix whole genome chicken microarrays. The forskolin-induced proliferation paradigm was chosen for this study given the robustness of the effect, which we reasoned would facilitate identification of the molecular events resulting in the production of new hair cells. RNA was isolated from control and forskolin-treated samples cultured for 24, 48 and 72 hours. BrdU incorporation, a marker for S-phase entry, is seen at 72 hours in the forskolin paradigm. Since S-phase entry is the criterion that distinguishes the avian epithelium from that in mammals, we have concentrated our analysis on this time point.

The isolated RNA from proliferating and non-proliferating auditory epithelia was amplified, biotinylated, hybridized to Affymetrix chicken arrays and scanned as described. The normalized data were then filtered by expression to exclude all probes that were expressed below 20% of the maximum expression level in at least one of the samples at each time point. This measure was taken to ensure that differential expression of genes expressed at very low levels were not assigned a disproportionate amount of significance. 32,509 out of 38,535 probes in the microarray met this initial criterion in the 72-hour samples (Table 2). Of the 32,509 filtered probes, 5,032 had a corrected p-value less than 0.05 at 72 hours. Of these 5,032 genes, 2,824 had at least a two-fold change in expression, which was used as a criterion for functional significance. The same statistics from the 24- and 48-hour samples are also shown in Table 2. A proportional Venn diagram clearly shows differential expression of far more genes at 72 hours than at 24 or 48 hours, but expression of most of the genes differentially expressed at the two earlier time points are also affected at 72 hours (Figure 1). It is not surprising that far fewer genes are differentially expressed at 24 and 48 hours than at 72 hours, given that there must be some latency following exposure to forskolin before gene expression is affected. A heat map revealed that directionality of fold-changes for specific genes remained fairly consistent at 24, 48, and 72 hours (Figure 2). There also appeared to be a general trend toward larger fold-change differences at the two later time points.

[Figure omitted. See PDF.]

Figure 1. Far more genes were differentially expressed at 72 hours than at 24 or 48 hours of exposure to forskolin.

Shown is a proportional Venn diagram indicating the number of genes significantly (p<0.05 and fold change >2) differentially expressed at 24, 48, and 72 hours of exposure to forskolin. The number of genes in each category is indicated in parentheses. Far more genes are affected at 72 hours than at either of the two earlier time points. A little over half of the genes differentially expressed at 24 or 48 hours are also differentially expressed at 72 hours. Only one gene, AMPH, was affected at all three time points.

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

[Figure omitted. See PDF.]

Figure 2. Fold-change differences of selected genes of interest remain stable at 24, 48, and 72 hours.

Shown is a heat map depicting fold-changes between the forskolin and control conditions for a set of genes of interest the three time points indicated. The genes shown were selected because of their known expression in the inner ear (i.e., RARB, CALB1, OCM), ion channel identity (e.g., KCNAB1, KCNJ15), or relevance to the present study (i.e., E2F1). The genes span a range of fold-changes and are both up and downregulated with forskolin treatment. Fold-changes are given as log2(expression in forskolin/expression in control). The directionality of fold-changes was fairly consistent across time points, but there was a slight general trend toward increasing fold-change values from 24 to 72 hours.

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

[Figure omitted. See PDF.]

Table 2. Microarray analysis statistics.

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

Many genes thought to be related to cellular proliferation and inner ear development were differentially expressed between the two conditions at multiple time points (Table 3). The complete lists of genes differentially expressed between the control and forskolin conditions at 24, 48, and 72 hours are included as Supporting Information (Tables S1, S2, and S3, respectively). Given that cellular proliferation is first observed in the basilar papilla at 72 hours in culture with forskolin, the following discussion will focus on data from this particular time point. One gene that is differentially expressed after 72 hours of exposure to forskolin is cyclin I (CCNI), a cell cycle regulator, which is over three-fold underexpressed in the forskolin samples at 72 hours and is downregulated in proliferating murine cardiomyocytes (Table 3) [51]. Wnt5a is underexpressed in regenerating tissue, a noteworthy observation as inactivation of this gene has been linked to cellular proliferation in a particular form of leukemia [52]. Wnt proteins are also thought to play a role in the development of the lateral line in zebrafish (reviewed in [53]) and chicken inner ears [54], [55], [56]. Frizzled homolog 10, a gene likely to play a role in cellular proliferation, is over 15-fold overexpressed in proliferating versus non-proliferating tissue after 72 hours in culture [57].

[Figure omitted. See PDF.]

Table 3. QPCR validation of gene expression data.

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

Differential gene expression is confirmed by qPCR

Quantitative PCR (qPCR) was performed on select genes in order to validate the microarray data. Primers were designed for twelve genes that were found by microarray to be differentially expressed between the two conditions at 72 hours; seven genes were upregulated and five were downregulated in the forskolin condition. QPCR confirmed the directionality of the microarray findings for nine out of the twelve genes tested (75%). The fold-changes of these nine genes were all significantly different from one (95% confidence interval, data not shown), confirming their differential expression as detected by microarray. These findings validate the initial microarray results, so further analyses were performed on the gene expression data set.

Curated gene set enrichment results

Our initial analysis of expression of individual genes that were differentially expressed revealed that many are known to be involved in control of the cell cycle and hair cell differentiation. To more systematically identify genes that were correlated with specific phenotypes, we used gene set enrichment analysis (GSEA), a computational method to determine the expression of genes associated with specific empirically determined phenotypes. Gene sets packages were downloaded from the Broad Institute website and were used for GSEAs on the Affymetrix expression data from 24, 48, and 72 hours (Table 1). An additional GSEA was performed comparing the microarray data from the 24-hour control samples to that in the 72-hour control samples to serve as a negative control for changes in gene expression resulting from the culture conditions rather than exposure to forskolin. This GSEA produced no significantly (FDR <0.25 and p-value <0.05) enriched gene sets for any of the three gene set packages used (data not shown).

To assess the functional relevance of the genes that were differentially expressed, a GSEA was performed using manually curated gene sets downloaded from the Broad Institute website. This package contains sets of genes found by individual experiments to be associated with specific phenotypes. The complete GSEA results for all utilized gene sets in this and all subsequent GSEAs are included as Supporting Information (Tables S4, S5, S6, S7, S8, S9, S10, S11, S12). Consistent with our observation of individual gene expression, the GSEA performed using curated gene sets showed that in forskolin treated basilar papillae many gene sets associated with cell cycle control are significantly enriched at various time points. For example, at 24, 48, and 72 hours there is enrichment of the following gene sets: SERUM_FIBROBLAST_CELLCYCLE [58], HSA04110_CELL_CYCLE [59], [60], [61], CELL_CYCLE [62], CELL_CYCLE_KEGG [59], [60], [61], GOLDRATH_CELLCYCLE [63], CELL_CYCLE_ CHECKPOINT, G1_TO_S_CELL_CYCLE_REACTOME [64], and BRENTANI_CELL_CYCLE [65] in explants exposed to forskolin (Table 4). The enrichment statistics for all gene sets used in this GSEA are included as Supporting Information (Tables S4, S5, S6). These findings are expected and instill confidence about the rest of the enrichment analysis results, given that forskolin treatment is known to stimulate S-phase entry in normally quiescent cells [22], [23].

[Figure omitted. See PDF.]

Table 4. Curated gene set enrichment analysis results.

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

The GSEA based on the curated gene sets also revealed that targets of E2F1, a transcription factor known to be important for cell cycle control, in human primary fibroblast cells identified by ChIP analysis (REN_E2F1_TARGETS [66]) are enriched in the proliferating auditory epithelia at both 48 and 72 hours. Further, a set of genes in the retinoblastoma pathway identified by microarray analysis in human cell lines (VERNELL_PRB_CLSTR1 [67]) was found to be significantly enriched after 24, 48, and 72 hours in culture with forskolin. These findings further underscore the potentially important role of the E2F1 transcription factor and the retinoblastoma protein pathway in auditory hair cell regeneration.

A particularly interesting result of the GSEA using the curated gene sets is the observation that a set of genes found to be upregulated in mouse embryonic, neural and hematopoietic stem cells (HSC) relative to differentiated brain and bone marrow cells (STEMCELL_COMMON_UP [68]) is also enriched in the basilar papilla after culture with forskolin for 24, 48, or 72 hours. Further, a set of genes found to be downregulated in stem cells (STEMCELL_COMMON_DN [68]) in the same study were also enriched in the quiescent tissue relative to proliferating tissue after 48 or 72 hours in culture, suggesting that genes downregulated in stem cells are also downregulated in the proliferating basilar papilla.

Though not widely accepted, there is evidence for the existence of stem cells within the postnatal organ of Corti [69]. Additionally, approximately 4% of supporting cells show stem cell-like behavior in that they undergo multiple rounds of division [70]. Recent work has shown that murine embryonic and induced pluripotent stem cells retain the capacity to differentiate into mechanosensitive cells which resemble immature hair cells [71]. It is interesting that genes found to be upregulated in HSCs specifically are also upregulated in proliferating basilar papillae given that HSCs are thought to remain in a state of quiescence until stimulated to proliferate by external signals [72], [73]. The chick basilar papilla is also quiescent until injury stimulates regeneration of hair cells by signals that have not yet been identified. It is therefore very interesting that proliferating basilar papillae overexpress genes that are also overexpressed in HSCs [68]. If the basilar papilla does in fact possess stem cells that share some properties with HSCs then it may be possible to stimulate the production of new hair cells using signals known to induce normally quiescent HSCs to proliferate. These results are consistent with the hypothesis that in forskolin treated basilar papillae some supporting cells are stimulated to dedifferentiate and express stem cell markers, prior to division and production of new hair cells.

It is worth noting that at 24 and 48 hours after exposure to forskolin, there is enrichment of the following curated gene sets, all of which are associated with cell death and/or apoptosis: BRENTANI_DEATH [65], HSA04210_APOPTOSIS [59], [60], [61], APOPTOSIS [64], APOPTOSIS_GENMAPP [64]. There is also enrichment of the gene set PASSERINI_APOPTOSIS [74] after 48 hours of exposure to forskolin. It is believed that forskolin induces proliferation in the inner ear without injury and subsequent regeneration, but this has never been definitively proven. However, the proliferation that is seen in the basilar papilla following exposure to forskolin is more robust than that seen following exposure to gentamicin, which causes hair cell death and regeneration [22]. Unpublished experiments performed for the same paper also did not show a qualitative loss of hair cells at 24, 48, or 72 hours, whereas hair cell ejection from the basilar papilla is seen as early as 30 hours following subcutaneous administration of gentamicin [75]. Markers of cell death in avian hair cells can appear even earlier following injury and begin with TIAR translocation as early as 12 hours after injury, before caspase activation which can occur as early as 30 hours after injury [76], [2]. Our gene expression data do not show differential expression of key markers of hair cell death such as caspase-3 and caspase-9 [76]. Regardless, the GSEA results suggest a possible contribution of apoptosis in forskolin induced proliferation in the avian inner ear. As the aim of the present study was to identify novel genes and pathways important for the production of new hair cells in deafened animals, the reported findings are of considerable interest regardless of a possible contribution of injury with subsequent regeneration in the forskolin-induced proliferation model of hair cell regeneration. Additional work will be needed to determine just how much, if any, of this proliferation is in fact the result of injury to the auditory epithelium.

Mammalian gene set enrichment results

To determine which miRNAs were likely upregulated to give rise to the observed pattern of differential gene expression detected by microarray, an additional GSEA was performed using the mammalian gene sets which are based on genomic conservation between the human, rat, mouse, and dog genomes. These gene sets were used on our chicken gene expression data because they included not only transcription factor binding site, but also miRNA binding sites, which are fairly well conserved between species [42]. Using these mammalian gene sets and a false discovery rate cut-off of 0.25 as suggested by [40], 210 out of 783 gene sets were significantly enriched in the control condition and 37 in the forskolin condition at 72 hours (Table 1). As a negative control, enrichment of predicted targets of the transcription factor PAX6 was examined, given its known role in retinal development [77], [78] and absence of expression in the developing chicken inner ear [79]. That a pathway thought not to be important in the inner ear was not activated by forskolin exposure underscores that significantly enriched gene sets reflect real changes in gene expression (Supporting Information).

The difference in the number of gene sets associated with the control and forskolin conditions is striking, and it is worth pointing out that many of the gene sets associated with the control condition are miRNA gene sets, each of which represents a group of genes that are predicted targets of a specific miRNA. The association of these sets with the control condition implies that the miRNA target genes are down-regulated following forskolin treatment. Thus, we hypothesize that these miRNAs are upregulated in the forskolin group. We put forward this entire set of miRNAs as candidates whose role in hair cell regeneration and cellular proliferation warrants further investigation. Specific miRNAs and transcription factors of interest based on our analysis using the mammalian gene sets are presented in Table 5. A complete list of all significantly enriched gene sets from this analysis are included in the Supporting Information section (Tables S7, S8, S9).

[Figure omitted. See PDF.]

Table 5. Mammalian gene set enrichment analysis results.

https://doi.org/10.1371/journal.pone.0011502.t005

Some of the miRNAs identified by enrichment analysis are interesting candidates for further study given what is already known about their function. For example, miR141 and miR200a are expressed in the chicken inner ear epithelium, and are important regulators of epithelial-mesenchymal transition [30], [80]. It is also interesting that the miR27 family showed up as one of the most significant gene sets in this particular analysis given the known role of this miRNA in cellular proliferation. Specifically, downregulation of this miRNA seems to decrease proliferation in the context of hepatic stellate cell activation, an early event in liver fibrosis [81]. MiR27a is overexpressed in human gastric adenocarcinoma, and inhibition of this miRNA limits proliferation [82]. Another study has also shown that the oncogenic activity of miR27a in human breast cancer cells is likely due to downregulation of genes controlling the G2-M cell-cycle checkpoint [83]. MiR181a is known to be expressed in the mammalian inner ear [28], and has a pro-proliferative role in human leukemia cells, which appears to be mediated by suppression of the cell cycle inhibitor p27 [33]. Given the abundance and nature of evidence for a role of miR181a, miR27a, and miR141/miR200a in cellular proliferation or inner ear development, we put these forward as particularly intriguing candidates for further study within the context of auditory hair cell regeneration.

Chicken/frog gene set enrichment analysis

While the primary goal of the GSEA with the mammalian gene sets was to identify miRNAs that may be important for hair cell regeneration, interestingly some transcription factor gene sets were found to be up- and downregulated in the forskolin condition. This was surprising in light of the fact that these sets were based on portions conserved across only selected mammalian genomes, which were expected to differ significantly from gene sets based on the chicken genome. To verify the up- and downregulation of these gene sets in the chicken inner ear an additional GSEA was performed using gene sets based on scanning the portions of the chicken genome that are conserved with frog (Xenopus tropicalis) to look for enrichment of sets of genes sharing particular transcription factor recognition sites obtained from the TRANSFAC database. Using a q-value cutoff of 0.25 there were 41 gene sets associate with forskolin treatment, and 331 gene sets associated with forskolin treatment (see Table 1). The enrichment statistics for all gene sets included in this analysis are presented as Supporting Information (Tables S10, S11, S12). Of the six gene sets that were significantly enriched in the forskolin condition, four are sets of genes that share known E2F1 binding sites (Table 6). Additionally, some E2F1 gene sets are significantly enriched in forskolin treated tissue at 24 and 48 hours relative to untreated controls (Table 6), suggesting that upregulation of predicted E2F1 targets occurs early after forskolin treatment and is sustained until 72 hours in culture, when new cells are first seen [22]. Many E2F/E2F1 gene sets are also enriched in the forskolin treated tissue when the mammalian gene sets (Table 5) are used to compare the expression data from forskolin samples fixed after 24 and 72 hours in culture, suggesting agreement between the GSEAs performed using the mammalian gene sets and chicken/frog conservation gene sets. The identification of E2F/E2F1 by both analyses strongly suggests a role for this transcription factor in forskolin induced proliferation in the chicken inner ear.

[Figure omitted. See PDF.]

Table 6. Conservation gene set enrichment analysis results.

https://doi.org/10.1371/journal.pone.0011502.t006

E2F1 is an intriguing candidate for further study given its known role in cell cycle control (reviewed in [84]) and potential role in inner ear development (reviewed in [85]). E2Fs bind hypophosphorylated pRBs which then sequester these factors. PRBs can be phosphorylated by cyclin-dependent kinases (CDKs), rendering these proteins unable to bind E2Fs thereby releasing these factors which then become transcriptionally active [86], [87], [88]. Interestingly, overexpression of E2F1 has been shown to induce quiescent REF-52 cells to enter S-phase [89]. E2F1 forms part of an ‘activating’ E2F complex that can interact with pocket proteins (pRBs) to control the cell cycle [90]. E2F1 exclusively binds the pRB RB1, and peaks in expression during the G1-S checkpoint [90]. Given the central role of E2Fs in cell cycle control, our finding that genes targeted by E2F1 are enriched among the genes that are overexpressed in the proliferating chick inner ear suggests a role for this transcription factor in hair cell regeneration.

Functional microRNA181a studies

Of all the miRNAs whose predicted targets were significantly enriched in control relative to forskolin treated basilar papillae, miR181a was specifically selected for functional studies in light of its high level of statistical significance, known expression in the developing mammalian inner ear [28], and pro-proliferative role in human leukemia cells, which appears to be mediated by suppression of the cell cycle inhibitor p27 [33]. To test whether miR181a overexpression is sufficient to stimulate proliferation, chick basilar papilla explants were transfected with either pre-miR181a, which is converted to mature miR181a within cells, or a non-targeting negative control miRNA. Incorporation of the thymidine analog BrdU was used to assay proliferation. Basilar papillae transfected with pre-miRNA had on average 68.75 (+/− SEM = 13.54) BrdU positive cells, significantly more than organs transfected with a non-targeting miRNA which averaged only 7.00 (+/− SEM = 10.23) BrdU positive cells (p<0.01, Figure 3). MiR181a overexpression is therefore sufficient to stimulate proliferation in the normally quiescent chick auditory epithelium.

[Figure omitted. See PDF.]

Figure 3. MicroRNA181a overexpression produces proliferation in the chicken basilar papilla.

Basilar papillae from 0-day-old chicks were cultured for 72 hours in the presence of BrdU, the nuclear incorporation of which was detected by antibody labeling. In A–E sensory epithelia are outlined with dashed white lines. The neural edge of each epithelium is at the bottom, the abneural edge is at the top, the distal end is to the right and the proximal end is to the left of the image. BrdU-positive (i.e., recently divided) cells are immunohistochemically labeled and appear as purple dots. Individual cochleas were treated with 100 nM pre-miR181a (A), 100 µM forskolin (FSK) (B), FSK + 100 nM pre-miR181a (C), 100 µM FSK+ 100 nM anti-miR181a (D), or DMSO as a vector control (E). A pro-proliferative effect of miR181a can be appreciated. Further, knocking down endogenous miR181a appears to have a large suppressive effect on forskolin induced proliferation. Transfection was achieved in the first 24 hours of culture using X-tremeGENE siRNA Transfection Reagent in accordance with the manufacturer's instructions (Roche, Indianapolis, IN). The miRNA containing media was removed after 24 hours and replaced with normal medium. The summary of the average number of BrdU positive cells for the control (n = 7), miR181a (n = 8), forskolin (n = 8), forskolin + miR181a (n = 8), and forskolin + anti-miR181a (n = 3) conditions. The following comparisons were all statistically significant: control versus miR181a (p = 0.001), control versus forskolin (p<0.001), miR181a versus forskolin (p<0.001), forskolin versus forskolin plus anti-miR181a (p = 0.008) and forskolin versus forskolin plus miR181a (p = 0.005). Scale bars in panels A–E = 0.2 mm.

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

Basilar papillae cultured with forskolin for 72 hours and transfected with a non-targeting miRNA had on average 178.25 (+/− SEM = 19.02) BrdU positive cells, significantly more than is seen with miR181a overexpression (p<0.05) which suggests that miR181a upregulation may account for only a portion of forskolin's mitogenic effect (Figure 3). To determine whether miR181a expression is necessary for forskolin induced proliferation, some basilar papillae were cultured with forskolin for 24 hours, then transfected with anti-miR181a to suppress endogenous miR181a, then cultured with forskolin for an additional 24 hours. Basilar papillae transfected with only anti-miR181a without forskolin treatment had an average of only 14 (+/− SEM = 3.06) BrdU positive cells per epithelium, not significantly different from controls. Basilar papillae exposed to forskolin and transfected with anti-miR181a had an average of 67.33 (+/− SEM = 13.38) BrdU positive cells, significantly less than the basilar papillae that were cultured with forskolin and transfected with a negative control miRNA (p<0.01), suggesting that endogenous miR181a contributes to the pro-proliferative effect of forskolin in the chick inner ear.

To see whether miR181a could enhance the already robust pro-proliferative effect of forskolin, basilar papillae were transfected with pre-miR181a and cultured in the presence of forskolin for 72 hours. These basilar papillae averaged 289.63 (+/− SEM = 28.04) BrdU positive cells per basilar papilla, which is significantly higher than the forskolin effect alone (p<0.01). This observation reaffirms the finding that endogenous miR181a is an important mediator of forskolin induced proliferation in the chick inner ear and is limited in part by the ability of forskolin to stimulate upregulation of miR181a.

Given the observation that miR181a overexpression is sufficient to stimulate proliferation in the chick auditory epithelium, basilar papillae were transfected with pre-miR181a or nontargeting negative control miRNA and then cryosectioned and labeled for BrdU and the early myosin VI to determine whether any of the newly produced BrdU positive cells express the early hair cell marker. As is seen in Figure 4, some cells labeling for BrdU 72 hours after transfection with pre-miR181a are located at the most apical portion of the basilar papilla, an area occupied by only hair cells and thin supporting cell processes. Further, some cells with nuclear BrdU labeling also have cytoplasmic myosin VI labeling. No such co-labeling cells were observed in basilar papillae transfected with a non-targeting miRNA. These results suggest that miR181a overexpression results in the production of new hair cells.

[Figure omitted. See PDF.]

Figure 4. Some new cells produced by miR181a transfection express the early hair cell marker myosin VI.

Shown is a confocal image of a basilar papillae that was transfected with pre-miR181a, cultured for 72 hours, then cryosectioned and labeled for BrdU (blue) and the early hair cell marker myosin VI (red) (A). The apical surface of hair cells is to the left of the image and the neural edge of the basilar papillae is toward the bottom. The tectorial membrane has been removed to allow for antibody labeling. The inset of A is also shown in panels B (myosin VI), C (BrdU) and D (merge). The arrow in D shows a cell that expresss myosin VI and has a BrdU positive nucleus, whereas the arrowhead shows a cell that is labeling for myosin VI but not BrdU. The double-labeled cell is one that has been stimulated to divide and subsequently begin differentiating toward a hair cell phenotype. The scale bar in B also applies to C and D.

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

It is not surprising that a single miRNA is able to stimulate S-phase entry in the avian auditory epithelium given the mounting evidence that miRNAs can affect cell cycle regulation in a variety of systems. MiR181a specifically has been shown to have a pro-proliferative effect in human myeloid leukemia cells [33]. Interestingly, this effect appears to be mediated in part by producing downregulation of p27, which may represent a barrier to hair cell regeneration in the mammalian cochlea [34], [35], [91], [92]. There is also evidence that many different miRNAs, including miR181a, are expressed in the developing mouse inner ear [28]. Despite these findings, studies on the role of miRNA in hair cell regeneration are limited in number and scope. One study revealed that the let-7 family of miRNA is downregulated in the regenerating newt auditory epithelium [32]. An interesting translational correlate to this finding is that let-7 reduces tumor growth in a rodent model of lung cancer, providing additional evidence for the antiproliferative effects of this particular family of miRNA [93]. Studies such as these underscore the importance of examining more closely the role of miRNA in hair cell regeneration, as induced underexpression of antiproliferative genes may ultimately allow regeneration of hair cells in mammals.

In summary, we have identified a possible role for miR181a in avian auditory hair cell regeneration. We have further identified additional miRNAs and transcription factors as candidates that warrant further investigation as potential targets for therapeutics aimed at replacing lost hair cells. Presented additionally are results suggestive of the presence of stem cells with properties similar to HSCs in the regenerating chick inner ear. Of all the candidates identified miR181a is specifically put forward for further study in in-vivo mammalian models of induced hearing loss in light of its demonstrated proliferative effect in the chicken inner ear as well as the fact that a miRNA based therapeutic could potentially be locally delivered to the middle ear and passively absorbed through the round window. This scenario seems plausible given that small molecules such as gentamicin can be delivered to the inner ear of patients in this fashion [94]. Computational analyses of our gene expression data suggest that the transcription factor E2F1 may also play a role in hair cell regeneration when many normally quiescent supporting cells are stimulated to re-enter the cell cycle and divide prior to differentiating into new hair cells. Additional functional studies will be necessary to ultimately determine what ability, if any, these players have to replace lost hair cells and result in recovery of compromised hearing thresholds in mammals.

Supporting Information

[Figure omitted. See PDF.]

Table S1.

Genes differentially expressed (p<0.05, fold-change >2) between control and forskolin treated basilar papillae after 24 hours in culture.

https://doi.org/10.1371/journal.pone.0011502.s001

(0.02 MB XLS)

Table S2.

Genes differentially expressed (p<0.05, fold-change >2) between control and forskolin treated basilar papillae after 48 hours in culture.

https://doi.org/10.1371/journal.pone.0011502.s002

(0.03 MB XLS)

Table S3.

Genes differentially expressed (p<0.05, fold-change >2) between control and forskolin treated basilar papillae after 72 hours in culture.

https://doi.org/10.1371/journal.pone.0011502.s003

(0.52 MB XLS)

Table S4.

Gene set enrichment analysis statistics. Curated gene sets enriched in the control and forskolin samples after 24 hours.

https://doi.org/10.1371/journal.pone.0011502.s004

(0.25 MB XLS)

Table S5.

Gene set enrichment analysis statistics. Curated gene sets enriched in the control and forskolin samples after 48 hours.

https://doi.org/10.1371/journal.pone.0011502.s005

(0.25 MB XLS)

Table S6.

Gene set enrichment analysis statistics. Curated gene sets enriched in the control and forskolin samples after 72 hours.

https://doi.org/10.1371/journal.pone.0011502.s006

(0.25 MB XLS)

Table S7.

Gene set enrichment analysis statistics. Mammalian gene sets enriched in the control and forskolin samples after 24 hours.

https://doi.org/10.1371/journal.pone.0011502.s007

(0.19 MB XLS)

Table S8.

Gene set enrichment analysis statistics. Mammalian gene sets enriched in the control and forskolin samples after 48 hours.

https://doi.org/10.1371/journal.pone.0011502.s008

(0.19 MB XLS)

Table S9.

Gene set enrichment analysis statistics. Mammalian gene sets enriched in the control and forskolin samples after 72 hours.

https://doi.org/10.1371/journal.pone.0011502.s009

(0.20 MB XLS)

Table S10.

Gene set enrichment analysis statistics. Chicken/frog conservation gene sets enriched in the control and forskolin samples after 24 hours.

https://doi.org/10.1371/journal.pone.0011502.s010

(0.15 MB XLS)

Table S11.

Gene set enrichment analysis statistics. Chicken/frog conservation gene sets enriched in the control and forskolin samples after 48 hours.

https://doi.org/10.1371/journal.pone.0011502.s011

(0.15 MB XLS)

Table S12.

Gene set enrichment analysis statistics. Chicken/frog conservation gene sets enriched in the control and forskolin samples after 72 hours.

https://doi.org/10.1371/journal.pone.0011502.s012

(0.15 MB XLS)

Acknowledgments

We thank Jennifer Stone, PhD, University of Washington, for reviewing the manuscript and providing helpful suggestions.

We thank Adam Thornberg and Shrikant Mane, PhD, Yale Keck Facility, for assistance with the RNA amplification and microarray hybridizations.

We thank Guy Richardson, PhD, University of Sussex, for helpful technical advice.

Author Contributions

Conceived and designed the experiments: CSF SHK JSS DSN. Performed the experiments: CSF DSN. Analyzed the data: CSF MU SHK JSS DSN. Contributed reagents/materials/analysis tools: MU JLD. Wrote the paper: CSF SHK JSS DSN.

Citation: Frucht CS, Uduman M, Duke JL, Kleinstein SH, Santos-Sacchi J, Navaratnam DS (2010) Gene Expression Analysis of Forskolin Treated Basilar Papillae Identifies MicroRNA181a as a Mediator of Proliferation. PLoS ONE5(7): e11502. https://doi.org/10.1371/journal.pone.0011502

References

1. WHO (2005) WHO | Deafness and hearing impairment. World Health Organization. WHO2005WHO | Deafness and hearing impairmentWorld Health Organization

2. Stone J, Cotanche D (2007) Hair cell regeneration in the avian auditory epithelium. Int J Dev Biol 51: 633–647.J. StoneD. Cotanche2007Hair cell regeneration in the avian auditory epithelium.Int J Dev Biol51633647

3. Jorgensen JM, Mathiesen C (1988) The avian inner ear. Continuous production of hair cells in vestibular sensory organs, but not in the auditory papilla. Naturwissenschaften 75: 319–320.JM JorgensenC. Mathiesen1988The avian inner ear. Continuous production of hair cells in vestibular sensory organs, but not in the auditory papilla.Naturwissenschaften75319320

4. Warchol ME, Corwin JT (1996) Regenerative proliferation in organ cultures of the avian cochlea: identification of the initial progenitors and determination of the latency of the proliferative response. J Neurosci 16: 5466–5477.ME WarcholJT Corwin1996Regenerative proliferation in organ cultures of the avian cochlea: identification of the initial progenitors and determination of the latency of the proliferative response.J Neurosci1654665477

5. Balak KJ, Corwin JT, Jones JE (1990) Regenerated hair cells can originate from supporting cell progeny: evidence from phototoxicity and laser ablation experiments in the lateral line system. J Neurosci 10: 2502–2512.KJ BalakJT CorwinJE Jones1990Regenerated hair cells can originate from supporting cell progeny: evidence from phototoxicity and laser ablation experiments in the lateral line system.J Neurosci1025022512

6. Stone JS, Cotanche DA (1994) Identification of the timing of S phase and the patterns of cell proliferation during hair cell regeneration in the chick cochlea. J Comp Neurol 341: 50–67.JS StoneDA Cotanche1994Identification of the timing of S phase and the patterns of cell proliferation during hair cell regeneration in the chick cochlea.J Comp Neurol3415067

7. Raphael Y (1992) Evidence for supporting cell mitosis in response to acoustic trauma in the avian inner ear. J Neurocytol 21: 663–671.Y. Raphael1992Evidence for supporting cell mitosis in response to acoustic trauma in the avian inner ear.J Neurocytol21663671

8. Corwin JT, Cotanche DA (1988) Regeneration of sensory hair cells after acoustic trauma. Science 240: 1772–1774.JT CorwinDA Cotanche1988Regeneration of sensory hair cells after acoustic trauma.Science24017721774

9. Ryals BM, Rubel EW (1988) Hair cell regeneration after acoustic trauma in adult Coturnix quail. Science 240: 1774–1776.BM RyalsEW Rubel1988Hair cell regeneration after acoustic trauma in adult Coturnix quail.Science24017741776

10. Stone JS, Cotanche DA (1992) Synchronization of hair cell regeneration in the chick cochlea following noise damage. J Cell Sci 102 ( Pt 4): 671–680.JS StoneDA Cotanche1992Synchronization of hair cell regeneration in the chick cochlea following noise damage.J Cell Sci102 ( Pt 4)671680

11. Lee KH, Cotanche DA (1996) Localization of the hair-cell-specific protein fimbrin during regeneration in the chicken cochlea. Audiol Neurootol 1: 41–53.KH LeeDA Cotanche1996Localization of the hair-cell-specific protein fimbrin during regeneration in the chicken cochlea.Audiol Neurootol14153

12. Shang J, Cafaro J, Nehmer R, Stone J (2010) Supporting Cell Division Is Not Required for Regeneration of Auditory Hair Cells After Ototoxic Injury In Vitro. J Assoc Res Otolaryngol. J. ShangJ. CafaroR. NehmerJ. Stone2010Supporting Cell Division Is Not Required for Regeneration of Auditory Hair Cells After Ototoxic Injury In Vitro.J Assoc Res Otolaryngol

13. Adler HJ, Raphael Y (1996) New hair cells arise from supporting cell conversion in the acoustically damaged chick inner ear. Neurosci Lett 205: 17–20.HJ AdlerY. Raphael1996New hair cells arise from supporting cell conversion in the acoustically damaged chick inner ear.Neurosci Lett2051720

14. Baird RA, Burton MD, Lysakowski A, Fashena DS, Naeger RA (2000) Hair cell recovery in mitotically blocked cultures of the bullfrog saccule. Proc Natl Acad Sci USA 97: 11722–11729.RA BairdMD BurtonA. LysakowskiDS FashenaRA Naeger2000Hair cell recovery in mitotically blocked cultures of the bullfrog saccule.Proc Natl Acad Sci USA971172211729

15. Baird RA, Steyger PS, Schuff NR (1996) Mitotic and nonmitotic hair cell regeneration in the bullfrog vestibular otolith organs. Annals of the New York Academy of Sciences 781: 59–70.RA BairdPS SteygerNR Schuff1996Mitotic and nonmitotic hair cell regeneration in the bullfrog vestibular otolith organs.Annals of the New York Academy of Sciences7815970

16. Taylor RR, Forge A (2005) Hair cell regeneration in sensory epithelia from the inner ear of a urodele amphibian. J Comp Neurol 484: 105–120.RR TaylorA. Forge2005Hair cell regeneration in sensory epithelia from the inner ear of a urodele amphibian.J Comp Neurol484105120

17. Duncan LJ, Mangiardi DA, Matsui JI, Anderson JK, McLaughlin-Williamson K, et al. (2006) Differential expression of unconventional myosins in apoptotic and regenerating chick hair cells confirms two regeneration mechanisms. J Comp Neurol 499: 691–701.LJ DuncanDA MangiardiJI MatsuiJK AndersonK. McLaughlin-Williamson2006Differential expression of unconventional myosins in apoptotic and regenerating chick hair cells confirms two regeneration mechanisms.J Comp Neurol499691701

18. Roberson DW, Alosi JA, Cotanche DA (2004) Direct transdifferentiation gives rise to the earliest new hair cells in regenerating avian auditory epithelium. J Neurosci Res 78: 461–471.DW RobersonJA AlosiDA Cotanche2004Direct transdifferentiation gives rise to the earliest new hair cells in regenerating avian auditory epithelium.J Neurosci Res78461471

19. Roberson DW, Kreig CS, Rubel EW (1996) Light microscopic evidence that direct transdifferentiation gives rise to new hair cells in regenerating avian auditory epithelium. Aud Neurosci 2: 195–205.DW RobersonCS KreigEW Rubel1996Light microscopic evidence that direct transdifferentiation gives rise to new hair cells in regenerating avian auditory epithelium.Aud Neurosci2195205

20. Saunders SS, Salvi RJ (1995) Pure tone masking patterns in adult chickens before and after recovery from acoustic trauma. J Acoust Soc Am 98: 1365–1371.SS SaundersRJ Salvi1995Pure tone masking patterns in adult chickens before and after recovery from acoustic trauma.J Acoust Soc Am9813651371

21. Gale JE, Meyers JR, Periasamy A, Corwin JT (2002) Survival of bundleless hair cells and subsequent bundle replacement in the bullfrog's saccule. J Neurobiol 50: 81–92.JE GaleJR MeyersA. PeriasamyJT Corwin2002Survival of bundleless hair cells and subsequent bundle replacement in the bullfrog's saccule.J Neurobiol508192

22. Navaratnam DS, Su HS, Scott SP, Oberholtzer JC (1996) Proliferation in the auditory receptor epithelium mediated by a cyclic AMP-dependent signaling pathway. Nat Med 2: 1136–1139.DS NavaratnamHS SuSP ScottJC Oberholtzer1996Proliferation in the auditory receptor epithelium mediated by a cyclic AMP-dependent signaling pathway.Nat Med211361139

23. Montcouquiol M, Corwin JT (2001) Brief treatments with forskolin enhance s-phase entry in balance epithelia from the ears of rats. J Neurosci 21: 974–982.M. MontcouquiolJT Corwin2001Brief treatments with forskolin enhance s-phase entry in balance epithelia from the ears of rats.J Neurosci21974982

24. Schuknecht HF (1964) Further Observations on the Pathology of Presbycusis. Arch Otolaryngol 80: 369–382.HF Schuknecht1964Further Observations on the Pathology of Presbycusis.Arch Otolaryngol80369382

25. Hawkins RD, Bashiardes S, Powder KE, Sajan SA, Bhonagiri V, et al. (2007) Large scale gene expression profiles of regenerating inner ear sensory epithelia. PLoS ONE 2: e525.RD HawkinsS. BashiardesKE PowderSA SajanV. Bhonagiri2007Large scale gene expression profiles of regenerating inner ear sensory epithelia.PLoS ONE2e525

26. Jackson RJ, Standart N (2007) How do microRNAs regulate gene expression? Sci STKE 2007: re1.RJ JacksonN. Standart2007How do microRNAs regulate gene expression?Sci STKE2007re1

27. Dalmay T (2008) Identification of genes targeted by microRNAs. Biochem Soc Trans 36: 1194–1196.T. Dalmay2008Identification of genes targeted by microRNAs.Biochem Soc Trans3611941196

28. Weston MD, Pierce ML, Rocha-Sanchez S, Beisel KW, Soukup GA (2006) MicroRNA gene expression in the mouse inner ear. Brain Res 1111: 95–104.MD WestonML PierceS. Rocha-SanchezKW BeiselGA Soukup2006MicroRNA gene expression in the mouse inner ear.Brain Res111195104

29. Sacheli R, Nguyen L, Borgs L, Vandenbosch R, Bodson M, et al. (2009) Expression patterns of miR-96, miR-182 and miR-183 in the developing inner ear. Gene Expression Patterns 9: 364–370.R. SacheliL. NguyenL. BorgsR. VandenboschM. Bodson2009Expression patterns of miR-96, miR-182 and miR-183 in the developing inner ear.Gene Expression Patterns9364370

30. Soukup GA (2009) Little but loud: Small RNAs have a resounding affect on ear development. Brain Research 1277: 104–114.GA Soukup2009Little but loud: Small RNAs have a resounding affect on ear development.Brain Research1277104114

31. Friedman LM, Dror AA, Mor E, Tenne T, Toren G, et al. (2009) MicroRNAs are essential for development and function of inner ear hair cells in vertebrates. Proc Natl Acad Sci USA 106: 7915–7920.LM FriedmanAA DrorE. MorT. TenneG. Toren2009MicroRNAs are essential for development and function of inner ear hair cells in vertebrates.Proc Natl Acad Sci USA10679157920

32. Tsonis PA, Call MK, Grogg MW, Sartor MA, Taylor RR, et al. (2007) MicroRNAs and regeneration: Let-7 members as potential regulators of dedifferentiation in lens and inner ear hair cell regeneration of the adult newt. Biochem Biophys Res Commun 362: 940–945.PA TsonisMK CallMW GroggMA SartorRR Taylor2007MicroRNAs and regeneration: Let-7 members as potential regulators of dedifferentiation in lens and inner ear hair cell regeneration of the adult newt.Biochem Biophys Res Commun362940945

33. Wang X, Gocek E, Liu CG, Studzinski GP (2009) MicroRNAs181 regulate the expression of p27Kip1 in human myeloid leukemia cells induced to differentiate by 1,25-dihydroxyvitamin D3. Cell Cycle 8: 736–741.X. WangE. GocekCG LiuGP Studzinski2009MicroRNAs181 regulate the expression of p27Kip1 in human myeloid leukemia cells induced to differentiate by 1,25-dihydroxyvitamin D3.Cell Cycle8736741

34. Lee YS, Liu F, Segil N (2006) A morphogenetic wave of p27Kip1 transcription directs cell cycle exit during organ of Corti development. Development 133: 2817–2826.YS LeeF. LiuN. Segil2006A morphogenetic wave of p27Kip1 transcription directs cell cycle exit during organ of Corti development.Development13328172826

35. Lowenheim H, Furness DN, Kil J, Zinn C, Gultig K, et al. (1999) Gene disruption of p27(Kip1) allows cell proliferation in the postnatal and adult organ of corti. Proc Natl Acad Sci U S A 96: 4084–4088.H. LowenheimDN FurnessJ. KilC. ZinnK. Gultig1999Gene disruption of p27(Kip1) allows cell proliferation in the postnatal and adult organ of corti.Proc Natl Acad Sci U S A9640844088

36. Affymetrix I (2004) Eukaryotic sample and array processing. GeneChip Expression Analysis Technical Manual. Santa Clara, CA: Affymetrix, Inc. pp. 2.1.3–2.3.18.I. Affymetrix2004Eukaryotic sample and array processing. GeneChip Expression Analysis Technical ManualSanta Clara, CAAffymetrix, Inc2.1.32.3.18

37. Affymetrix I (2009) GeneChip. Expression Analysis Technical Manual. I. Affymetrix2009GeneChipExpression Analysis Technical Manual

38. Schmittgen TD, Livak KJ (2008) Analyzing real-time PCR data by the comparative C(T) method. Nat Protoc 3: 1101–1108.TD SchmittgenKJ Livak2008Analyzing real-time PCR data by the comparative C(T) method.Nat Protoc311011108

39. Rozen S, Skaletsky H (2000) Primer3 on the WWW for general users and for biologist programmers. Methods Mol Biol 132: 365–386.S. RozenH. Skaletsky2000Primer3 on the WWW for general users and for biologist programmers.Methods Mol Biol132365386

40. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, et al. (2005) Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A 102: 15545–15550.A. SubramanianP. TamayoVK MoothaS. MukherjeeBL Ebert2005Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles.Proc Natl Acad Sci U S A1021554515550

41. Xie X, Lu J, Kulbokas EJ, Golub TR, Mootha V, et al. (2005) Systematic discovery of regulatory motifs in human promoters and 3′ UTRs by comparison of several mammals. Nature 434: 338–345.X. XieJ. LuEJ KulbokasTR GolubV. Mootha2005Systematic discovery of regulatory motifs in human promoters and 3′ UTRs by comparison of several mammals.Nature434338345

42. Friedman RC, Farh KK, Burge CB, Bartel DP (2009) Most mammalian mRNAs are conserved targets of microRNAs. Genome Res 19: 92–105.RC FriedmanKK FarhCB BurgeDP Bartel2009Most mammalian mRNAs are conserved targets of microRNAs.Genome Res1992105

43. Wasserman WW, Palumbo M, Thompson W, Fickett JW, Lawrence CE (2000) Human-mouse genome comparisons to locate regulatory sites. Nat Genet 26: 225–228.WW WassermanM. PalumboW. ThompsonJW FickettCE Lawrence2000Human-mouse genome comparisons to locate regulatory sites.Nat Genet26225228

44. Schmidt D, Wilson MD, Ballester B, Schwalie PC, Brown GD, et al. Five-Vertebrate ChIP-seq Reveals the Evolutionary Dynamics of Transcription Factor Binding. Science. D. SchmidtMD WilsonB. BallesterPC SchwalieGD BrownFive-Vertebrate ChIP-seq Reveals the Evolutionary Dynamics of Transcription Factor Binding.Science

45. Borneman AR, Gianoulis TA, Zhang ZD, Yu H, Rozowsky J, et al. (2007) Divergence of transcription factor binding sites across related yeast species. Science 317: 815–819.AR BornemanTA GianoulisZD ZhangH. YuJ. Rozowsky2007Divergence of transcription factor binding sites across related yeast species.Science317815819

46. Hua S, Kallen CB, Dhar R, Baquero MT, Mason CE, et al. (2008) Genomic analysis of estrogen cascade reveals histone variant H2A.Z associated with breast cancer progression. Mol Syst Biol 4: 188.S. HuaCB KallenR. DharMT BaqueroCE Mason2008Genomic analysis of estrogen cascade reveals histone variant H2A.Z associated with breast cancer progression.Mol Syst Biol4188

47. Pruitt KD, Tatusova T, Maglott DR (2005) NCBI Reference Sequence (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Res 33: D501–504.KD PruittT. TatusovaDR Maglott2005NCBI Reference Sequence (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins.Nucleic Acids Res33D501504

48. Blanchette M, Kent WJ, Riemer C, Elnitski L, Smit AF, et al. (2004) Aligning multiple genomic sequences with the threaded blockset aligner. Genome Res 14: 708–715.M. BlanchetteWJ KentC. RiemerL. ElnitskiAF Smit2004Aligning multiple genomic sequences with the threaded blockset aligner.Genome Res14708715

49. Kel AE, Gossling E, Reuter I, Cheremushkin E, Kel-Margoulis OV, et al. (2003) MATCH: A tool for searching transcription factor binding sites in DNA sequences. Nucleic Acids Res 31: 3576–3579.AE KelE. GosslingI. ReuterE. CheremushkinOV Kel-Margoulis2003MATCH: A tool for searching transcription factor binding sites in DNA sequences.Nucleic Acids Res3135763579

50. Matys V, Fricke E, Geffers R, Gossling E, Haubrock M, et al. (2003) TRANSFAC: transcriptional regulation, from patterns to profiles. Nucleic Acids Res 31: 374–378.V. MatysE. FrickeR. GeffersE. GosslingM. Haubrock2003TRANSFAC: transcriptional regulation, from patterns to profiles.Nucleic Acids Res31374378

51. Liu Y, Tang MK, Cai DQ, Li M, Wong WM, et al. (2007) Cyclin I and p53 are differentially expressed during the terminal differentiation of the postnatal mouse heart. Proteomics 7: 23–32.Y. LiuMK TangDQ CaiM. LiWM Wong2007Cyclin I and p53 are differentially expressed during the terminal differentiation of the postnatal mouse heart.Proteomics72332

52. Roman-Gomez J, Jimenez-Velasco A, Cordeu L, Vilas-Zornoza A, San Jose-Eneriz E, et al. (2007) WNT5A, a putative tumour suppressor of lymphoid malignancies, is inactivated by aberrant methylation in acute lymphoblastic leukaemia. Eur J Cancer 43: 2736–2746.J. Roman-GomezA. Jimenez-VelascoL. CordeuA. Vilas-ZornozaE. San Jose-Eneriz2007WNT5A, a putative tumour suppressor of lymphoid malignancies, is inactivated by aberrant methylation in acute lymphoblastic leukaemia.Eur J Cancer4327362746

53. Ma EY, Raible DW (2009) Signaling Pathways Regulating Zebrafish Lateral Line Development. Current Biology 19: R381–R386.EY MaDW Raible2009Signaling Pathways Regulating Zebrafish Lateral Line Development.Current Biology19R381R386

54. Sienknecht UJ, Fekete DM (2008) Comprehensive Wnt-related gene expression during cochlear duct development in chicken. J Comp Neurol 510: 378–395.UJ SienknechtDM Fekete2008Comprehensive Wnt-related gene expression during cochlear duct development in chicken.J Comp Neurol510378395

55. Sienknecht UJ, Fekete DM (2009) Mapping of Wnt, frizzled, and Wnt inhibitor gene expression domains in the avian otic primordium. J Comp Neurol 517: 751–764.UJ SienknechtDM Fekete2009Mapping of Wnt, frizzled, and Wnt inhibitor gene expression domains in the avian otic primordium.J Comp Neurol517751764

56. Stevens CB, Davies AL, Battista S, Lewis JH, Fekete DM (2003) Forced activation of Wnt signaling alters morphogenesis and sensory organ identity in the chicken inner ear. Developmental Biology 261: 149–164.CB StevensAL DaviesS. BattistaJH LewisDM Fekete2003Forced activation of Wnt signaling alters morphogenesis and sensory organ identity in the chicken inner ear.Developmental Biology261149164

57. Yang-Snyder J, Miller JR, Brown JD, Lai CJ, Moon RT (1996) A frizzled homolog functions in a vertebrate Wnt signaling pathway. Curr Biol 6: 1302–1306.J. Yang-SnyderJR MillerJD BrownCJ LaiRT Moon1996A frizzled homolog functions in a vertebrate Wnt signaling pathway.Curr Biol613021306

58. Chang HY, Sneddon JB, Alizadeh AA, Sood R, West RB, et al. (2004) Gene expression signature of fibroblast serum response predicts human cancer progression: similarities between tumors and wounds. PLoS Biol 2: E7.HY ChangJB SneddonAA AlizadehR. SoodRB West2004Gene expression signature of fibroblast serum response predicts human cancer progression: similarities between tumors and wounds.PLoS Biol2E7

59. Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, et al. (2008) KEGG for linking genomes to life and the environment. Nucleic Acids Res 36: D480–484.M. KanehisaM. ArakiS. GotoM. HattoriM. Hirakawa2008KEGG for linking genomes to life and the environment.Nucleic Acids Res36D480484

60. Kanehisa M, Goto S (2000) KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Research 28: 27–30.M. KanehisaS. Goto2000KEGG: kyoto encyclopedia of genes and genomes.Nucleic Acids Research282730

61. Kanehisa M, Goto S, Hattori M, Aoki-Kinoshita KF, Itoh M, et al. (2006) From genomics to chemical genomics: new developments in KEGG. Nucleic Acids Research 34: D354–357.M. KanehisaS. GotoM. HattoriKF Aoki-KinoshitaM. Itoh2006From genomics to chemical genomics: new developments in KEGG.Nucleic Acids Research34D354357

62. Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, et al. (2000) Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet 25: 25–29.M. AshburnerCA BallJA BlakeD. BotsteinH. Butler2000Gene ontology: tool for the unification of biology. The Gene Ontology Consortium.Nat Genet252529

63. Goldrath AW, Luckey CJ, Park R, Benoist C, Mathis D (2004) The molecular program induced in T cells undergoing homeostatic proliferation. Proc Natl Acad Sci USA 101: 16885–16890.AW GoldrathCJ LuckeyR. ParkC. BenoistD. Mathis2004The molecular program induced in T cells undergoing homeostatic proliferation.Proc Natl Acad Sci USA1011688516890

64. Doniger SW, Salomonis N, Dahlquist KD, Vranizan K, Lawlor SC, et al. (2003) MAPPFinder: using Gene Ontology and GenMAPP to create a global gene-expression profile from microarray data. Genome Biol 4: R7.SW DonigerN. SalomonisKD DahlquistK. VranizanSC Lawlor2003MAPPFinder: using Gene Ontology and GenMAPP to create a global gene-expression profile from microarray data.Genome Biol4R7

65. Brentani H, Caballero OL, Camargo AA, da Silva AM, da Silva WA, et al. (2003) The generation and utilization of a cancer-oriented representation of the human transcriptome by using expressed sequence tags. Proc Natl Acad Sci USA 100: 13418–13423.H. BrentaniOL CaballeroAA CamargoAM da SilvaWA da Silva2003The generation and utilization of a cancer-oriented representation of the human transcriptome by using expressed sequence tags.Proc Natl Acad Sci USA1001341813423

66. Ren B, Cam H, Takahashi Y, Volkert T, Terragni J, et al. (2002) E2F integrates cell cycle progression with DNA repair, replication, and G(2)/M checkpoints. Genes Dev 16: 245–256.B. RenH. CamY. TakahashiT. VolkertJ. Terragni2002E2F integrates cell cycle progression with DNA repair, replication, and G(2)/M checkpoints.Genes Dev16245256

67. Vernell R, Helin K, Müller H (2003) Identification of target genes of the p16INK4A-pRB-E2F pathway. J Biol Chem 278: 46124–46137.R. VernellK. HelinH. Müller2003Identification of target genes of the p16INK4A-pRB-E2F pathway.J Biol Chem2784612446137

68. Ramalho-Santos M, Yoon S, Matsuzaki Y, Mulligan RC, Melton DA (2002) “Stemness”: transcriptional profiling of embryonic and adult stem cells. Science 298: 597–600.M. Ramalho-SantosS. YoonY. MatsuzakiRC MulliganDA Melton2002“Stemness”: transcriptional profiling of embryonic and adult stem cells.Science298597600

69. Oshima K, Grimm CM, Corrales CE, Senn P, Martinez Monedero R, et al. (2007) Differential distribution of stem cells in the auditory and vestibular organs of the inner ear. J Assoc Res Otolaryngol 8: 18–31.K. OshimaCM GrimmCE CorralesP. SennR. Martinez Monedero2007Differential distribution of stem cells in the auditory and vestibular organs of the inner ear.J Assoc Res Otolaryngol81831

70. Stone JS, Rubel EW (1999) Delta1 expression during avian hair cell regeneration. Development 126: 961–973.JS StoneEW Rubel1999Delta1 expression during avian hair cell regeneration.Development126961973

71. Oshima K, Shin K, Diensthuber M, Peng AW, Ricci AJ, et al. (2010) Mechanosensitive Hair Cell-like Cells from Embryonic and Induced Pluripotent Stem Cells. Cell 141: 704–716.K. OshimaK. ShinM. DiensthuberAW PengAJ Ricci2010Mechanosensitive Hair Cell-like Cells from Embryonic and Induced Pluripotent Stem Cells.Cell141704716

72. Jude CD, Gaudet JJ, Speck NA, Ernst P (2008) Leukemia and hematopoietic stem cells: balancing proliferation and quiescence. Cell Cycle 7: 586–591.CD JudeJJ GaudetNA SpeckP. Ernst2008Leukemia and hematopoietic stem cells: balancing proliferation and quiescence.Cell Cycle7586591

73. Narbonne P, Roy R (2008) Genes that affect both cell growth and polarity mediate stem cell quiescence. Front Biosci 13: 995–1002.P. NarbonneR. Roy2008Genes that affect both cell growth and polarity mediate stem cell quiescence.Front Biosci139951002

74. Passerini AG, Polacek DC, Shi C, Francesco NM, Manduchi E, et al. (2004) Coexisting proinflammatory and antioxidative endothelial transcription profiles in a disturbed flow region of the adult porcine aorta. Proc Natl Acad Sci USA 101: 2482–2487.AG PasseriniDC PolacekC. ShiNM FrancescoE. Manduchi2004Coexisting proinflammatory and antioxidative endothelial transcription profiles in a disturbed flow region of the adult porcine aorta.Proc Natl Acad Sci USA10124822487

75. Mangiardi DA, McLaughlin-Williamson K, May KE, Messana EP, Mountain DC, et al. (2004) Progression of hair cell ejection and molecular markers of apoptosis in the avian cochlea following gentamicin treatment. J Comp Neurol 475: 1–18.DA MangiardiK. McLaughlin-WilliamsonKE MayEP MessanaDC Mountain2004Progression of hair cell ejection and molecular markers of apoptosis in the avian cochlea following gentamicin treatment.J Comp Neurol475118

76. Mangiardi DA, McLaughlin-Williamson K, May KE, Messana EP, Mountain DC, et al. (2004) Progression of hair cell ejection and molecular markers of apoptosis in the avian cochlea following gentamicin treatment. J Comp Neurol 475: 1–18.DA MangiardiK. McLaughlin-WilliamsonKE MayEP MessanaDC Mountain2004Progression of hair cell ejection and molecular markers of apoptosis in the avian cochlea following gentamicin treatment.J Comp Neurol475118

77. Ashery-Padan R, Gruss P (2001) Pax6 lights-up the way for eye development. Curr Opin Cell Biol 13: 706–714.R. Ashery-PadanP. Gruss2001Pax6 lights-up the way for eye development.Curr Opin Cell Biol13706714

78. Kumar JP (2001) Signalling pathways in Drosophila and vertebrate retinal development. Nat Rev Genet 2: 846–857.JP Kumar2001Signalling pathways in Drosophila and vertebrate retinal development.Nat Rev Genet2846857

79. Heanue TA, Davis RJ, Rowitch DH, Kispert A, McMahon AP, et al. (2002) Dach1, a vertebrate homologue of Drosophila dachshund, is expressed in the developing eye and ear of both chick and mouse and is regulated independently of Pax and Eya genes. Mech Dev 111: 75–87.TA HeanueRJ DavisDH RowitchA. KispertAP McMahon2002Dach1, a vertebrate homologue of Drosophila dachshund, is expressed in the developing eye and ear of both chick and mouse and is regulated independently of Pax and Eya genes.Mech Dev1117587

80. Gregory PA, Bracken CP, Bert AG, Goodall GJ (2008) MicroRNAs as regulators of epithelial-mesenchymal transition. Cell Cycle 7: 3112–3118.PA GregoryCP BrackenAG BertGJ Goodall2008MicroRNAs as regulators of epithelial-mesenchymal transition.Cell Cycle731123118

81. Ji J, Zhang J, Huang G, Qian J, Wang X, et al. (2009) Over-expressed microRNA-27a and 27b influence fat accumulation and cell proliferation during rat hepatic stellate cell activation. FEBS Lett 583: 759–766.J. JiJ. ZhangG. HuangJ. QianX. Wang2009Over-expressed microRNA-27a and 27b influence fat accumulation and cell proliferation during rat hepatic stellate cell activation.FEBS Lett583759766

82. Liu T, Tang H, Lang Y, Liu M, Li X (2009) MicroRNA-27a functions as an oncogene in gastric adenocarcinoma by targeting prohibitin. Cancer Lett 273: 233–242.T. LiuH. TangY. LangM. LiuX. Li2009MicroRNA-27a functions as an oncogene in gastric adenocarcinoma by targeting prohibitin.Cancer Lett273233242

83. Mertens-Talcott S, Chintharlapalli S, Li X, Safe S (2007) The Oncogenic microRNA-27a Targets Genes That Regulate Specificity Protein Transcription Factors and the G2-M Checkpoint in MDA-MB-231 Breast Cancer Cells. Cancer Res 67: 11001–11011.S. Mertens-TalcottS. ChintharlapalliX. LiS. Safe2007The Oncogenic microRNA-27a Targets Genes That Regulate Specificity Protein Transcription Factors and the G2-M Checkpoint in MDA-MB-231 Breast Cancer Cells.Cancer Res671100111011

84. Johnson DG, Degregori J (2006) Putting the Oncogenic and Tumor Suppressive Activities of E2F into Context. Curr Mol Med 6: 731–738.DG JohnsonJ. Degregori2006Putting the Oncogenic and Tumor Suppressive Activities of E2F into Context.Curr Mol Med6731738

85. Rocha-Sanchez SM, Beisel KW (2007) Pocket proteins and cell cycle regulation in inner ear development. Int J Dev Biol 51: 585–595.SM Rocha-SanchezKW Beisel2007Pocket proteins and cell cycle regulation in inner ear development.Int J Dev Biol51585595

86. Cobrinik D (2005) Pocket proteins and cell cycle control. Oncogene 24: 2796–2809.D. Cobrinik2005Pocket proteins and cell cycle control.Oncogene2427962809

87. Dannenberg JH, Schuijff L, Dekker M, van der Valk M, te Riele H (2004) Tissue-specific tumor suppressor activity of retinoblastoma gene homologs p107 and p130. Genes Dev 18: 2952–2962.JH DannenbergL. SchuijffM. DekkerM. van der ValkH. te Riele2004Tissue-specific tumor suppressor activity of retinoblastoma gene homologs p107 and p130.Genes Dev1829522962

88. Li JM, Hu PP, Shen X, Yu Y, Wang XF (1997) E2F4-RB and E2F4-p107 complexes suppress gene expression by transforming growth factor beta through E2F binding sites. Proc Natl Acad Sci U S A 94: 4948–4953.JM LiPP HuX. ShenY. YuXF Wang1997E2F4-RB and E2F4-p107 complexes suppress gene expression by transforming growth factor beta through E2F binding sites.Proc Natl Acad Sci U S A9449484953

89. Johnson DG, Schwarz JK, Cress WD, Nevins JR (1993) Expression of transcription factor E2F1 induces quiescent cells to enter S phase. Nature 365: 349–352.DG JohnsonJK SchwarzWD CressJR Nevins1993Expression of transcription factor E2F1 induces quiescent cells to enter S phase.Nature365349352

90. Dyson N (1998) The regulation of E2F by pRB-family proteins. Genes Dev 12: 2245–2262.N. Dyson1998The regulation of E2F by pRB-family proteins.Genes Dev1222452262

91. Chen P, Segil N (1999) p27(Kip1) links cell proliferation to morphogenesis in the developing organ of Corti. Development 126: 1581–1590.P. ChenN. Segil1999p27(Kip1) links cell proliferation to morphogenesis in the developing organ of Corti.Development12615811590

92. White PM, Doetzlhofer A, Lee YS, Groves AK, Segil N (2006) Mammalian cochlear supporting cells can divide and trans-differentiate into hair cells. Nature 441: 984–987.PM WhiteA. DoetzlhoferYS LeeAK GrovesN. Segil2006Mammalian cochlear supporting cells can divide and trans-differentiate into hair cells.Nature441984987

93. Esquela-Kerscher A, Trang P, Wiggins JF, Patrawala L, Cheng A, et al. (2008) The let-7 microRNA reduces tumor growth in mouse models of lung cancer. Cell Cycle 7: 759–764.A. Esquela-KerscherP. TrangJF WigginsL. PatrawalaA. Cheng2008The let-7 microRNA reduces tumor growth in mouse models of lung cancer.Cell Cycle7759764

94. Nedzelski JM, Chiong CM, Fradet G, Schessel DA, Bryce GE, et al. (1993) Intratympanic gentamicin instillation as treatment of unilateral Menière's disease: update of an ongoing study. The American journal of otology 14: 278–282.JM NedzelskiCM ChiongG. FradetDA SchesselGE Bryce1993Intratympanic gentamicin instillation as treatment of unilateral Menière's disease: update of an ongoing study.The American journal of otology14278282