A microRNA network regulates proliferative timing and extracellular matrix synthesis during cellular quiescence in fibroblasts
Genome Biology volume 13, Article number: R121 (2012)
Although quiescence (reversible cell cycle arrest) is a key part in the life history and fate of many mammalian cell types, the mechanisms of gene regulation in quiescent cells are poorly understood. We sought to clarify the role of microRNAs as regulators of the cellular functions of quiescent human fibroblasts.
Using microarrays, we discovered that the expression of the majority of profiled microRNAs differed between proliferating and quiescent fibroblasts. Fibroblasts induced into quiescence by contact inhibition or serum starvation had similar microRNA profiles, indicating common changes induced by distinct quiescence signals. By analyzing the gene expression patterns of microRNA target genes with quiescence, we discovered a strong regulatory function for miR-29, which is downregulated with quiescence. Using microarrays and immunoblotting, we confirmed that miR-29 targets genes encoding collagen and other extracellular matrix proteins and that those target genes are induced in quiescence. In addition, overexpression of miR-29 resulted in more rapid cell cycle re-entry from quiescence. We also found that let-7 and miR-125 were upregulated in quiescent cells. Overexpression of either one alone resulted in slower cell cycle re-entry from quiescence, while the combination of both together slowed cell cycle re-entry even further.
microRNAs regulate key aspects of fibroblast quiescence including the proliferative state of the cells as well as their gene expression profiles, in particular, the induction of extracellular matrix proteins in quiescent fibroblasts.
When mammalian cells are in an environment unfavorable for continued proliferation, they can exit the cell cycle in early to mid-G1 phase at the 'restriction point'  and enter a reversible, out-of-cell cycle state denoted 'quiescence'. Many cells in the human body are quiescent, and the ability of cells to exit the cell cycle but retain their capacity to re-enter the cell cycle as needed (for instance, when required to replenish a cell lineage, mount an immune response, or heal a wound) is central to normal physiology. Failures in this process may underlie a wide range of pathologies including excessive scarring, fibrotic disease, chronic wounding, and cancer, yet we have a poor understanding of the changes that occur when cells become quiescent or the molecular basis for these changes.
Widespread gene expression changes occur when cells enter quiescence, including both repression and activation of genes [2–9]. These changes can vary among cell types and in response to different antiproliferative signals, but there are also commonalities in different types of quiescence and in different quiescent cell types [2, 7–9]. Several important regulators of the gene expression changes that occur with quiescence have been described, including the MYC and E2F family transcription factors that coordinate cell cycle re-entry and repress cell cycle genes during quiescence [5, 10–13], and the HES1 transcriptional repressor that preserves the reversibility of quiescence . There are also hundreds of genes that are upregulated when cells become quiescent, whose possible regulators include forkhead transcription factors [15, 16], ELK1, NF-κB, MEF2, IRF, AP-1, SALL2, and MXI1 . Despite these proposed factors, however, the drivers and mechanisms of many of the gene expression changes in quiescence are still not known.
In addition to regulation of quiescence by transcription factors, there is likely also regulation of quiescence gene expression changes at the post-transcriptional level. microRNAs are 20 to 23 nucleotide non-coding RNAs that regulate a wide variety of transcripts post-transcriptionally by inducing transcript degradation or inhibiting protein translation [17–19]. microRNAs have been implicated in a wide range of biological processes related to quiescence, including cell proliferation control, stem cell renewal, developmental timing, and cancer . Medina and colleagues, for example, discovered that four microRNAs were upregulated and over 100 microRNAs were downregulated as T98G glioblastoma cells progress from quiescence into the proliferative cell cycle . They and others demonstrated that miR-221 and miR-222 target the cyclin-dependent kinase inhibitors p27Kip1 and p57Kip2, such that overexpression of miR-221 and miR-222 during growth factor deprivation induces S-phase entry and triggers cell death [21–25].
Another example is the let-7 family of microRNAs, members of which are important regulators of cellular differentiation [26–34] and proliferation [29, 35–37] in mammals, C. elegans, and Drosophila melanogaster. let-7 family members can behave as tumor suppressors and antagonize oncogenes such as MYC and RAS [28, 35, 38–45].
As a final example, the miR-17-92 cluster of six microRNAs, which is induced by the MYC oncogene , can itself act as an oncogene. Enforced expression of the miR-17-92 cluster, in concert with MYC expression, can accelerate tumor development in a mouse B-cell lymphoma model . While MYC can induce transcription of E2F transcription factors, two of the members of the microRNA cluster, miR-17-5p and miR-20a, negatively regulate levels of E2F1, demonstrating a complex network of interactions that may affect the cell's commitment to proliferation or apoptosis [46–51].
We investigated the role of microRNAs in a fibroblast model of quiescence and discovered that microRNA expression is broadly and similarly altered by two different quiescence signals: contact inhibition and serum withdrawal. We further found that microRNAs regulate some of the changes in gene expression and cellular function associated with quiescence, as well as the transition between proliferation and quiescence.
microRNAs exhibit a strong quiescence signature
We have developed a model system of quiescence in primary human fibroblasts in which quiescence can be induced by either serum starvation or contact inhibition. Either condition results in an accumulation of quiescent cells, as indicated by cell cycle markers and RNA content . Using one-color microRNA microarrays, we monitored microRNA expression levels in proliferating, serum-starved, and contact-inhibited primary human dermal fibroblasts (Figure 1A). Among the 209 microRNAs detected above background, 142 (68%) were expressed at different levels in proliferating compared with either serum-starved or contact-inhibited fibroblasts at a false-discovery rate of 1% (Figure 1B). microRNA expression patterns for contact inhibition and serum starvation were extremely similar, with a 95% confidence interval (CI) Pearson's correlation of 0.952 to 0.975, much more so than the mRNA expression patterns for the same conditions (Pearson's correlation of 0.319 to 0.341, 95% CI) (Additional File 1, Figure S1). This large difference in the amount of correlation between quiescence states may be due to experimental design or microarray platform differences, but an alternative explanation is that microRNAs exhibit more of a common quiescence signature than protein-coding transcripts . microRNAs downregulated in quiescent cells included miR-18, miR-20, miR-29, and miR-7, and microRNAs upregulated with quiescence included let-7b, miR-125a, miR-30, miR-181, miR-26, and miR-199. With a stringent cutoff of greater than two-fold expression change due to quiescence, eight microRNAs were expressed at higher levels in proliferating cells and eight were expressed at higher levels in quiescent cells (Additional File 1, Table S1).
We sought to validate the changes in microRNA levels with an independent method. In collaboration with Rosetta Inpharmatics, we used massively parallel, multiplexed qRT-PCR  to monitor the abundance of 219 microRNAs in fibroblasts collected during proliferation or after 4 days of serum starvation. There was strong agreement between the fold-change values obtained via the microarray and the multiplex qRT-PCR (Pearson's correlation 0.504 to 0.751, 95% CI) (Additional File 1, Figure S2).
Targets of microRNAs change with quiescence
In order to identify microRNAs with a functional, regulatory role in quiescence, we analyzed the gene expression patterns of microRNA target genes in two whole-genome mRNA microarray timecourses comparing proliferating cells to cells induced into quiescence by contact inhibition or serum starvation (Figure 2A). In one timecourse, fibroblasts were made quiescent by serum withdrawal for 4 days and then re-stimulated with serum for 48 h . In another, fibroblasts were sampled after 7 or 14 days of contact inhibition . Using singular value decomposition of the combined timecourses, we found that the strongest orthonormal gene expression pattern ('eigengene') correlated with the proliferative state of the cell (Additional File 1, Figure S3B). This eigengene explained approximately 40% of the gene expression variation (Additional File 1, Figure S3A). The linear projection of each gene to that eigengene gave a 'proliferation index' for each gene that summarized its association with proliferation or quiescence. For each microRNA, we averaged the proliferation indexes of its predicted target genes as provided by the TargetScan algorithm [55, 56] and assigned a P value to that mean using bootstrap resampling (Figure 2B). The miR-29 family's targets had the most statistically extreme mean proliferation index, with a P value <10-4 (the lowest P value possible based on the 104 bootstrap resamplings taken). miR-29 expression is strongly associated with proliferation (Additional File 1, Figure S4), and its predicted targets are upregulated by both methods of quiescence induction.
Besides miR-29, however, there were few microRNAs with strongly anti-correlated target genes. There are multiple possible explanations. First, expression levels and activity need not be completely correlated, as microRNA activity can be affected by the cooperation or antagonism of RNA-binding proteins [57–60] as well as changing mRNA abundance, dynamics, and primary and secondary structure [61–66]. Second, the microRNAs may be affecting translation rate but not transcript abundance, in which case their effects would not be detectable by microarray analysis. Finally, many of the microRNAs investigated likely regulate too few genes to be considered significant by this whole-genome target analysis, since a small list of targets can lead to artificially low statistical significance by bootstrap analysis. Indeed, some microRNAs might regulate a small number of critical genes and thereby produce an important functional effect even without a statistically significant change in the average proliferation index for all of its targets. For these reasons, we chose to investigate further miR-29 and other candidates identified based on their previously reported associations with proliferation and cell cycle regulation: let-7 [35, 36] and miR-125 [32, 33].
miR-29 regulates collagen and collagen-chaperone genes
Gene ontology analysis of predicted, evolutionarily conserved miR-29 targets revealed an enrichment for multiple categories including collagen fibril organization and extracellular matrix formation (Additional File 1, Table S3), indicating that miR-29 most likely regulates extracellular matrix (ECM) biosynthesis in fibroblasts, consistent with previous reports on miR-29 in fibroblasts and other cell types [67–72]. We identified miR-29 targets in dermal fibroblasts by overexpressing miR-29 in asynchronously proliferating fibroblasts and analyzing the ensuing changes in gene expression by microarray analysis. As expected, genes predicted to be miR-29 targets by TargetScan were more likely to be repressed by miR-29 overexpression than genes not predicted to be miR-29 targets (Figure 3B). We identified genes that both changed significantly in the microarray analysis and contained predicted miR-29 binding sites. Of the 15 genes that met these criteria, nine are involved in extracellular matrix formation (Figure 3A and Table 1). When we plotted the behavior of these same genes in the serum starvation and contact inhibition microarray timecourse data, we discovered that these genes display a quiescence-associated gene expression pattern. The genes encoding miR-29 targets followed a general pattern of increasing expression as fibroblasts are serum-starved, decreasing expression as they are restimulated, and highest expression in cells that were contact-inhibited for 7 or 14 days (Figure 3C). These genes were therefore highly anti-correlated with the pattern of expression for miR-29 itself (Additional File 1, Figure S4). These results suggest that the downregulation of miR-29 expression levels in quiescent fibroblasts is an important contributor to the induction of extracellular matrix genes with quiescence.
We sought to confirm whether miR-29 regulates not just transcript abundance, but also protein levels of extracellular matrix components in quiescent cells. We investigated three proteins encoded by miR-29 targets (collagen I, collagen III, and collagen VI) by immunoblot analysis of protein lysates isolated from proliferating cells and cells made quiescent by mitogen (PDGF) withdrawal or contact inhibition. As anticipated, all three proteins were upregulated in both quiescence conditions compared with proliferating cells. These three miR-29 targets were also strongly repressed at the protein level by transfection of miR-29 as compared to transfection of a negative control, non-targeting microRNA, while protein levels of GAPDH and α-tubulin (two proteins from genes not targeted by miR-29) were unaffected (Figure 3D).
Autocrine TGF-ß is unlikely to mediate miR-29 expression changes in quiescence
TGF-ß signaling leads to an increase in collagen synthesis  and can repress miR-29 [69, 74, 75]. We confirmed that exogenous addition of TGF-ß repressed miR-29 expression, as measured by qRT-PCR (Additional File 1, Figure S5A), in our dermal fibroblast model. Although exogenous TGF-ß can downregulate miR-29, immunoblots for Smad3 phosphorylation levels showed no significant difference in autocrine TGF-ß signaling between proliferating and quiescent fibroblasts (Additional File 1, Figure S5B), indicating that the TGF-ß signaling pathway is unlikely to be responsible for the reduction in miR-29 expression in quiescent fibroblasts. In addition, although TGF-ß can regulate collagen expression independently of miR-29 [76, 77], the similar phospho-Smad3 levels in proliferating and quiescent fibroblasts implies that changes in TGF-ß activity are unlikely to significantly regulate collagen biosynthesis in quiescence, further emphasizing the importance of miR-29 as a regulator of quiescence-associated changes in ECM expression.
miR-29 hastens cell cycle re-entry from quiescence
We also tested whether miR-29 has a role in the cell cycle transition between proliferation and quiescence by simultaneously restimulating serum-starved fibroblasts to proliferate with full serum medium and transfecting them with miR-29. Over the next 36 h, we quantified by flow cytometry the rate of EdU nucleotide analogue incorporation by the cells and their overall DNA content, which allowed us to assign cells to G0/G1, S, and G2/M phases of the cell cycle . When compared to cells transfected with a control non-targeting microRNA, cells transfected with miR-29 contained fewer cells in G0/G1 and more cells in S phase at 20 and 24 h post transfection (Figure 4A, for 20 and 24 h timepoints, respectively). At 28 and 32 h after transfection, cells transfected with miR-29 contained fewer cells in S phase and more cells in G2/M phase than those transfected with the control ( for 28 h timepoint). miR-29 overexpression thus hastens re-entry into the cell cycle from a quiescent state.
To further explore the effects of miR-29 expression on the cell cycle, we transfected miR-29 or a negative control microRNA into asynchronously cycling fibroblasts. Forty-eight hours post transfection, miR-29 transfection led to more cells in G2/M (Figure 4B). As expected considering that cells in the G2/M phase tend to be larger than cells in other phases of the cell cycle, miR-29 transfection also led to larger cells (Figure 4D). Further experimentation revealed that miR-29 transfection resulted in fewer cells than the negative control transfection (Figure 4C, ). Thus, miR-29 transfection in proliferating cells led to G2/M arrest rather than increased mitosis. This may reflect the activity of a miR-29 target gene; indeed, one target, RCC2 (TD-60), is repressed about 57% upon miR-29 transfection (Figure 3A and Table 1), and it plays an essential role in progression through metaphase .
let-7 and miR-125 non-redundantly delay cell cycle entry from quiescence
let-7 plays roles in differentiation, cancer, and the cell cycle, as discussed above. In C. elegens, the lin-4 microRNA (miR-125 in mammals) acts in the same heterochronic pathway of temporal differentiation as let-7 . The two microRNAs are also frequently located together in microRNA clusters across many phylogenetic lineages . In multiple species, they are co-regulated and share partly overlapping roles during development [81–86]. Because both let-7 and miR-125 are upregulated in quiescence, we investigated whether let-7 and miR-125 have complementary roles in cell cycle regulation.
We monitored the functional roles of let-7 and miR-125 on cell cycle re-entry from quiescence using the same method we used for miR-29 as described above. Compared with control-transfected cells, cells transfected with let-7 contained an elevated fraction of cells in the G0/G1 phase at 20 and 24 h post transfection and fewer cells in S phase at 20 h post transfection (Figure 4A, for 20 and 24 h timepoints, respectively), indicating that cell cycle re-entry is delayed by let-7 overexpression. By 32 h post transfection, the let-7-overexpressing population contained more cells in the G2/M phase than control cells , as we have previously reported . We observed an even stronger effect on cell cycle re-entry with miR-125 than for let-7. At 20 and 24 h after transfection, cells transfected with miR-125 contained more cells in G0/G1 and fewer cells in S phase than controls (Figure 4A, for 20 and 24 h timepoints, respectively). To assess whether let-7 and miR-125 have complementary effects on cell cycle progression, we overexpressed a combination of the two microRNAs. Overexpression of let-7 and miR-125 together resulted in a further accumulation of cells in G0/G1 and even slower S phase entry than either individually (Figure 4A, compared to let-7 and miR-125, respectively, at the 20 h timepoint), implying that their cell cycle effects are non-redundant and complementary.
A microRNA quiescence program
While the predominant view of quiescent cells is that they are inactive or 'shut down', our data from several different lines of experimentation indicate that the transition into quiescence in fibroblasts is a highly regulated and active process [2, 14, 52]. We previously reported , and we again confirmed by our SVD analysis of quiescence gene expression timecourse data (Figure 2), that entry into quiescence in fibroblasts is associated with large-scale remodeling of gene expression patterns affecting a significant fraction of all genes within the genome, with comparable numbers of genes both increasing and decreasing in expression. We show here that entry into quiescence is also associated with widespread changes in the abundance of a significant number of microRNAs. microRNAs both increase and decrease in abundance upon entry into quiescence, similar to the effects on mRNA expression.
One clear distinction between microRNAs and mRNAs was noticed: while gene expression patterns have both a common component and a signal-specific component  (Figure 2A), microRNA patterns with quiescence were very similar for samples made quiescent by two distinct quiescence signals (contact inhibition and serum starvation). This finding is in accord with previous studies that indicated that microRNA profiles are extremely informative about a human cancer's developmental lineage and differentiation state, and that microRNAs are particularly valuable for classifying poorly differentiated tumors [87, 88]. Indeed, our data suggest that there may be a quiescence microRNA program that is stronger and more consistent than a quiescence gene expression program. Such a signature may facilitate the identification of universal quiescence-related pathways.
The complementarity of let-7 and miR-125
In many organisms, lin-4 (miR-125) and let-7 are both important for developmental programs involving differentiation or cell cycle arrest [26, 31]. Low levels of let-7, for example, are associated with pluripotency and proliferation, while higher let-7 levels are associated with cell cycle exit and differentiation [28, 34, 89]. In vertebrates, mature let-7 and miR-125 are largely absent from early embryos and are induced upon differentiation [84–86]. We previously reported that let-7 targets the E2 ubiquitin ligase CDC34 and that let-7 overexpression in fibroblasts results in a G2/M arrest . Here we show that, when overexpressed, both miR-125 and let-7 specifically affect the ability of quiescent fibroblasts to re-enter the proliferative cell cycle from quiescence induced by serum starvation.
Our data and the literature, taken together, support a model in which miR-125 and let-7 family members are induced upon the commitment to a cell state lineage or reversible cell cycle exit. During differentiation or quiescence, let-7 and miR-125 may actively suppress the expression of cell cycle-associated transcripts through a post-transcriptional mechanism that reinforces the out-of-cycle state established by transcriptional mechanisms. Possible candidates for these transcripts include previously reported cell cycle targets of let-7 such as RAS , CCND1 , CDC25 , and CDC34 , and miR-125 targets such as BCL3  and ETS1 . Our results indicate that in reversibly arrested cells, miR-125 and let-7 downregulate cell proliferation-promoting genes. Upon restimulation, these genes are released from let-7 and miR-125-mediated repression and are required for normal cell cycle re-entry.
Although miR-125 and let-7 are co-conserved and co-regulated in many organisms, the two microRNAs also share some overlapping target genes [33, 93, 94], which suggests the possibility that some of the functional effects on the cell cycle exerted by each microRNA are redundant. Our results demonstrate that introduction of both microRNAs together had a stronger effect on cell cycle re-entry than introduction of either one alone, suggesting that they cooperate and play non-redundant roles in suppressing the expression of proliferation-associated genes in quiescent cells. This finding helps to explain the strong evolutionary selection to retain both microRNAs. Exogenous delivery of the let-7 microRNA has been shown to cause regression of murine lung tumors through an effect on cell proliferation distinct from apoptosis . Our data indicate that administration of miR-125 or a combination of let-7 and miR-125 might have even greater effects.
miR-29's role in quiescence
One of the functional changes that we previously observed in quiescent fibroblasts is an overall induction of extracellular matrix proteins . We report here that downregulation of the microRNA miR-29 is likely regulating the induction of extracellular matrix protein expression with quiescence: as miR-29 levels decline with quiescence, levels of miR-29 targets increase, and miR-29 overexpression represses the levels of these targets. Reporter assays by multiple independent groups have found in several different cell types that miR-29 directly targets collagens COL1A1, COL3A1, and COL4A2 in a seed sequence-dependent manner [95–97]. Based on those studies and our microarray and immunoblot results, miR-29 likely also represses collagens directly in proliferating fibroblasts. The findings place miR-29 among the very few molecules discovered, along with FoxO [98–100], and FoxP [101, 102] transcription factors, and the regulators of miR-29 itself, to regulate the induction (as opposed to the repression) of genes in quiescent cells. Because our data indicate that the activity of the TGF-ß signaling pathway is similar in proliferating and quiescent fibroblasts, it is not likely that TGF-ß is regulating the changes in miR-29 expression between these states. Other possible candidates for miR-29 transcriptional regulation include NF-κB and sonic hedgehog [70, 103]. Further study is necessary to elucidate which factors are responsible in quiescence.
Repression of RCC2 could explain the G2/M arrest phenotype seen with miR-29 transfection. Targets identified in other model systems could also be relevant. miR-29 targeting of DNA methyltransferases 3A and 3B, for example, can inhibit lung cancer cell tumorigenicity . miR-29 can also induce apoptosis in cholangiocarcioma cells via the miR-29 target MCL-1 , and induce replicative senescence in HeLa cells by targeting B-MYB .
We suggest that the role of miR-29 in hastening cell cycle re-entry, however, may reflect its effects not on validated cell cycle regulators, but instead on extracellular matrix proteins. Quiescent cells, in general, are relieved of the biosynthetic requirement of synthesizing the constituents of new cells, but in our fibroblast model system they also retain a comparable rate of metabolic activity as proliferating fibroblasts . Indeed, we discovered that fibroblasts express increased levels of several extracellular matrix proteins during quiescence compared with proliferation  (Figures 3C and 3D). From this perspective, it is particularly interesting that miR-29 overexpression results in more rapid cell cycle entry. Although miR-29 has been reported to be an oncogene (transgenic mice overexpressing miR-29 in their B cells develop B-cell chronic lymphocytic leukemia ) our microarray data revealed no clear candidate cell cycle genes that would explain the early re-entry phenotype we observed in our model system.
We suggest an alternative possibility: relieved of the commitment to translate and fold extracellular matrix proteins like collagen, miR-29-overexpressing cells may be able to commit more rapidly to the cell cycle. If a competition exists for translational resources between the synthesis of proteins required for cell duplication and the synthesis of proteins targeted for secretory pathways, then miR-29 may be able to direct resources between those two processes depending on the proliferative state of the cell. Further studies, especially on fibroblast cell lines derived from patients with idiopathic pulmonary fibrosis, which are characterized by excessive secretion of extracellular matrix proteins [108, 109], will be able to elucidate whether miR-29 is an important regulator of a tradeoff between proliferative and secretory modes.
Our data indicate that quiescence is associated with widespread, consistent changes in microRNA abundance. The regulated microRNAs contribute to gene expression programs that form the characteristic attributes of quiescent cells by reinforcing the non-proliferative nature of the cells and also regulating their cell-type specific roles. As such, further investigation into microRNAs should lead to a greater understanding of both universal aspects of quiescence programs as well as the regulation of processes specific to a quiescent cell's in vivo roles. Our results support some of the ongoing efforts to administer microRNAs to patients of cancer and fibrotic disease and suggest some new strategies.
Materials and methods
We isolated primary fibroblasts from neonatal human foreskin tissue samples provided by the National Disease Research Interchange (NDRI) as described in the supplementary methods for Legesse-Miller et al.  We routinely cultured the fibroblasts aseptically at 37°C with 5% CO2 in high-glucose DMEM with 4.5 mM glutamine (Life Technologies) supplemented with 10% (v/v) fetal bovine serum (FBS) (Hyclone) and 100 μg/mL penicillin and streptomycin (Life Technologies). Cells were serum-starved by reducing the serum concentration to 0.1% (v/v). To generate contact-inhibited samples, we plated fibroblasts and changed their culture medium regularly (every 2 or 3 days) without passaging them.
Three isolates of dermal fibroblasts were harvested in proliferative conditions, that is, sparsely subcultured 2 days before harvest, after 4 days of serum starvation, or after 7 days of contact inhibition. Cells were harvested by trypsinization, centrifuged at 160 × g, and snap-frozen in liquid nitrogen. Total RNA was isolated from the frozen cells using the mirVana miRNA isolation kit (Life Technologies). RNA quality was confirmed using a Bioanalyzer 2100 (Agilent Technology) and the concentration was determined with a NanoDrop spectrophotometer (NanoDrop Technologies). 100 ng of each sample was 3′-labeled with Cy3-pCp in two separate reactions and hybridized to microarray slides using the Agilent microRNA microarray kit (Agilent, G4470A). Microarray features were extracted with Feature Extractor 22.214.171.124. We normalized arrays for total intensity and then regressed each gene's expression using the model
where i denotes the index for a microRNA, Q, S, C 1 , and C 2 are annotations for quiescence, serum starvation, and the different fibroblast cell isolates, respectively, and SVA denotes the one significant surrogate variable we found as described below. Y i is the measured log2 expression for microRNA i and m i is its baseline expression. The x variables are the given experimental variables (indexed by subscripts) with values 0 or 1, the B coefficients are the gene-specific responses to a particular x variable, and E is the error term. Surrogate variable analysis (SVA) was performed with the R package from Leek et al. , giving the one significant surrogate variable we included in the multiple regression analysis. Differential expression due to quiescence was determined with an F-test for the significance of the microRNA's response to variable x Q , with a false-discovery rate of 1% deemed statistically significant. microRNAs without statistically significant gene expression change from quiescence were not shown in Figure 1A and 1B.
We denoted the overall biological response to serum starvation and contact inhibition (plotted log2 transformed in both Figure 1A as the heat-map intensities and Figure 1B along the x-axis) as the sum of the responses and the residuals . The Pearson correlation coefficient was calculated comparing these values in the serum starvation and contact inhibition conditions.
Multiplexed real-time PCR for microRNA expression levels
We collected primary human fibroblasts over a timecourse during serum starvation. Copy number of each microRNA per 10 pg of total RNA was determined using the protocol described in . In summary, RNA was extracted using the mirVana microRNA isolation kit as described above, and a tailed, gene-specific primer was used to convert the RNA template into cDNA with a universal PCR binding site at one end. The resulting primer-extended, full-length cDNA was amplified in a highly multiplexed manner for 219 individual microRNAs. Real-time PCR was performed with a combination of an LNA-containing microRNA/siRNA-specific 'reverse' primer and a generic universal primer complementary to the universal binding site introduced during reverse transcription. Amplification was monitored with SYBR green fluorescence. The cycle number at which the signal exceeded the background was used to determine the absolute abundance of the monitored microRNA in the sample. The Pearson's correlation between the real-time PCR data and the microRNA microarray data was determined between the 4-day serum-starved data point for the qRT-PCR and the mean of the serum starvation responses from the multiple regression for the microarray.
Gene expression microarrays for quiescence and mir-29 targets
Contact-inhibited fibroblast gene expression microarrays and serum starvation and restimulation arrays have been previous described [52, 54]. To summarize briefly, total RNA was isolated from proliferating, serum-starved, and serum-stimulated fibroblasts as described above for the microRNA microarray. Total RNA from each sample, 325 ng each, was amplified and labeled using the Low RNA Input Fluorescent Labeling Kit (Agilent Technologies) to incorporate Cyanine 3-CTP (Cy-3) or Cyanine 5-CTP (Cy-5). Cy-3-labeled time zero samples were used as a reference for serum withdrawal samples, which were labeled with Cyanine 5-CTP. For serum stimulation, 4-day serum-starved fibroblasts were labeled with Cy-3 and stimulated samples were labeled with Cy-5. Labeled cRNA was mixed and co-hybridized to whole Human Genome Oligo Microarray slides (Agilent Technologies) at 60°C for 17 h and subsequently washed with the Agilent Oligo Microarray Hybridization Kit. Slides were scanned with a dual laser scanner (Agilent Technologies). The Agilent feature extraction software, in conjunction with the Princeton University Microarray database, was used to compute the log ratio of the difference between the two samples for each gene after background subtraction and dye normalization. Of the approximately 44,000 probes on the microarray, probes that generated signal in at least 80% of arrays were identified. Fluorescence data for each probe were mapped to genes based on UniGene Clusters. If multiple probes mapped to a single gene, the values were averaged. The Pearson correlation coefficient was computed between the 96h serum-starved sample and the mean of the 7d contact-inhibited samples. Hierarchical clustering was performed on centered genes via centroid linkage, and four clusters were chosen based on the resulting dendrogram. Gene ontology (GO) term enrichment was determined using the Generic Gene Ontology Term Finder . Qualitatively non-informative or redundant GO terms (for example, 'Biological process' or 'cell cycle' vs. 'cell cycle process') were removed to give a selected subset.
For miR-29 overexpression microarrays, fibroblasts were transfected as described below with Pre-miR miR-29b or Negative Control #2 oligonucleotide duplexes (Life Technologies). Forty-eight hours after transfection, total RNA from the cells was harvested and hybridized to microarrays as above. The experiments were repeated on three different dermal fibroblast isolates. Target genes annotated by TargetScan 5.1 [55, 56, 112] were considered well-conserved miR-29 targets if PCT >0.5. A one-sided t-test was used to calculate the significance of the log2 fold change between the miR-29b transfection and the control, and a gene was declared 'changing' if it was repressed greater than two-fold at 5% FDR.
Singular value decomposition to identify microRNAs with significantly changing predicted targets
The matrix of gene expression arrays was filtered to exclude genes with missing values in any array, and this matrix was decomposed by singular value decomposition (SVD) to obtain 16 eigengenes. Each gene's expression profile was then linearly projected onto the first eigengene to obtain one summarizing number, dubbed the 'proliferation index', as genes with a strong positive projection tend to be associated with proliferation and genes with a strong negative projection tend to be associated with quiescence. Sets of computationally-predicted target genes were obtained from TargetScan by excluding all predictions with context scores >-0.5 (negative numbers indicate more confident predictions). The mean projection of each of these target gene sets and its additive inverse were used as two-tailed test statistics on a null hypothesis distribution of 10,000 mean projections of randomly sampled gene sets. Each sample gene set was the same size as the original target gene set for which the linear projection was calculated.
Overexpression of microRNA mimics
Proliferating or 4-day serum-starved primary fibroblasts were reverse-transfected using Oligofectamine (Life Technologies) with a 50 nM final concentration of Pre-miR microRNA duplexes let-7b, miR-125a, miR-29a, a 1:1 combination of let-7b and miR-125a, or the Negative Control #2 non-targeting control (Life Technologies). The microRNA duplexes and Oligofectamine were diluted in OptiMEM I (Life Technologies) and incubated at room temperature for 15 min. Human fibroblasts were trypsinized, washed, and then re-suspended in OptiMEM I at a concentration of 375,000 cells/mL. One milliliter of the transfection mixture was added to 4 mL of the cell suspension and plated on a 10 cm plate. The cells were incubated for 4 h and then supplemented with 5 mL of DMEM with 20% FBS. Twenty-four hours post transfection the medium was changed to DMEM containing 10% FBS.
For the serum-restimulation timecourses, we measured the duration of serum restimulation from the moment at which DMEM with 20% FBS was added. These experiments were done in triplicate on two different days (six timecourses in total). Standard error was calculated for both G0/G1 and S phase percentages at each timepoint as the square root of the total sum of square residuals from the mean percentage on each day. Proliferating cells were harvested 48 h after transfection for the assays described below.
Cell cycle progression assay
We determined cell cycle phases using Click-iT EdU Alexa Fluor 488 according to the protocol in . Briefly, we added 10 μL of a 10 mM EdU solution (Life Technologies) in phosphate-buffered saline (PBS) (Life Technologies) directly to 10 mL of culture medium on fibroblasts for a final concentration of 10 μM. We incubated the cells for 2 h with the EdU, and then trypsinized and re-suspended them to 1 × 107 cells/mL in PBS containing 1% bovine serum albumin (BSA) (Amresco). A total of 100 μL of this cell suspension was added to 100 μL of freshly prepared 4% formaldehyde in PBS (Thermo Scientific) and incubated in the dark at room temperature for 15 min. Three milliliters of PBS with 1% BSA was added to quench the fixation. The cells were then resuspended in 100 μL of PBS containing 1% BSA and added to 100 μL of 0.2% Triton X-100 in PBS. We added to each sample 500 μL of Click-iT reaction cocktail: 100 mM Tris-Cl, pH 8.5, 2 mM CuSO4, 10 μM Alexa Fluor 488 azide (Life Technologies), and 100 mM ascorbic acid. The mixture was incubated in the dark at room temperature for 30 min. Two milliliters of wash buffer (1% BSA, 0.2% Triton X-100 in PBS) was added, the cells were pelleted at 200 × g for 5 min, and the supernatant was discarded. We then resuspended the labeled cells in 500 μL of DAPI solution containing 1 μg/mL of DAPI in 0.1% Triton X-100 in PBS and analyzed them by flow cytometry on an LSR II flow cytometer (BD Biosciences, San Jose, CA, USA). DAPI was excited at 345 nm and its emission was detected at 458 nm. Alexa Fluor 488 was excited at 494 nm and its emission was detected at 519 nm.
Statistical significance of the changes was determined using a Dirichlet likelihood ratio test that is similar to a test of independence. To summarize, the cell cycle phase proportions at each timepoint for each microRNA transfection were fit to a maximum likelihood Dirichlet distribution by an iterated, alternating mean/precision estimation method . The distributions and their log likelihoods were calculated for the null hypothesis of identical Dirichlet distributions and the alternative hypothesis of two different Dirichlet distributions for the negative control transfection and the microRNA transfection of interest. The log likelihoods of the two hypotheses were compared using the test statistic
where C is the flow cytometry data at the particular timepoint. D was then evaluated on the distribution for three degrees of freedom to calculate a P value (the alternative hypothesis calculates two three-parameter Dirichlet distributions instead of one).
Cell size and number analysis
Dermal fibroblasts were trypsinized and resuspended in PBS, and cell size was measured in triplicate for each sample using the Beckman Coulter counter. Cell numbers were determined using the Countess automated cell counter (Invitrogen). For miR-29 transfection, cell numbers were evaluated using a one-sided t test.
Immunoblotting for miR-29 targets
Fibroblasts were reverse transfected with miR-29b or a negative control microRNA as above, but cells were plated at either 7,500 cells/cm2 (proliferating and mitogen-starved conditions) or 750,000 cells/cm2 (contact-inhibited condition). Twenty-four hours post transfection, cells were washed with warm PBS and then switched to low-serum conditions for collecting extracellular matrix proteins: FBM (Lonza), insulin (Lonza), and 0.1% FBS (v/v). Proliferating and contact-inhibited conditions were additionally supplemented with 30 ng/mL recombinant human PDGF-BB.
After culturing for 4 days in low-serum medium, intracellular proteins were collected by washing cells in ice-cold PBS followed by scraping cells into a solution of 4% SDS, 100 mM Tris-HCl pH 7.5, 1 mM DTT, and an EDTA-free protease inhibitor cocktail (Roche). Lysates were vortexed, heated briefly (5-10 min at 95°C), sonicated until the solutions became non-viscous, and then centrifuged at 10,000 × g for 10 min. The soluble lysates were transferred to new tubes and insoluble pellets were discarded. Determination of soluble lysate concentration and immunoblotting conditions were otherwise identical to those previously described in Lemons et al. and Pollina et al. [52, 54]. Antibodies and specific blotting conditions used are described below.
TGF-ß treatment and signaling analysis
Cells were treated with recombinant human TGF-ß1 (Life Technologies) for 48 h in low-serum medium: high glucose DMEM with 4.5 mM glutamine, serum replacement (Sigma-Aldrich, S2640), and 30 ng/mL PDGF-BB. Cells were lysed in TRIzol, and RNA was harvested with the Direct-zol RNA Miniprep Kit (Zymo Research). The expression of miR-29b was measured in cell samples on a small scale normalized to miR-100 using the miRCURY LNA Universal RT microRNA PCR kit (Exiqon) with miR-29b and miR-100 primer sets and an ABI 7900 real-time PCR system. Total RNA samples were extracted using TRIzol reagent (Life Technologies) along with the Direct-zol RNA Miniprep kit (Zymo Research). Relative expression changes were quantified in triplicate using the method on SYBR green fluorescence. Cell lysates from proliferating, 4 days serum-starved, and 7 days contact-inhibited fibroblasts were harvested according to the procedures above.
The following primary antibodies were used for immunoblotting: rabbit polyclonal IgG against collagen I (Calbiochem, 234167), rabbit polyclonal IgG against COL3A1 (Santa Cruz Biotechnology, sc-28888), biotinylated rabbit polyclonal IgG against Collagen VI (Acris Antibodies, R1043B), rabbit monoclonal IgG against Phospho-Smad3 Ser423/425 (Cell Signaling Technology, 9520), rabbit monoclonal IgG against α-Tubulin (Cell Signaling Technology, 2125), and rabbit polyclonal IgG against GAPDH (Abcam, ab9485). Each antibody was diluted in Tris-buffered saline containing 0.1% Tween-20 and 5% BSA and incubated with immunoblot membranes overnight at 4°C.
The microarray data generated for this study (the microRNA microarrays and the miR-29 overexpression microarrays) have been deposited in the NCBI Gene Expression Omnibus (GEO)  as one SuperSeries under the accession number GSE42614. Serum starvation/restimulation timecourse microarrays  and contact inhibition microarrays  were published in prior studies and are available in GEO with accessions GSE42681 and GSE42612, respectively.
confidence interval or contact inhibition
false discovery rate
quantitative reverse-transcription polymerase chain reaction
Pardee AB: A restriction point for control of normal animal cell proliferation. Proc Natl Acad Sci USA. 1974, 71: 1286-1290. 10.1073/pnas.71.4.1286.
Coller HA, Sang L, Roberts JM: A new description of cellular quiescence. PLoS Biology. 2006, 4: e83-10.1371/journal.pbio.0040083.
Coppock DL, Kopman C, Scandalis S, Gilleran S: Preferential gene expression in quiescent human lung fibroblasts. Cell Growth Differ. 1993, 4: 483-493.
Iyer VR, Eisen MB, Ross DT, Schuler G, Moore T, Lee JC, Trent JM, Staudt LM, Hudson J, Boguski MS, Lashkari D, Shalon D, Botstein D, Brown PO: The transcriptional program in the response of human fibroblasts to serum. Science. 1999, 283: 83-87. 10.1126/science.283.5398.83.
Liu H, Adler AS, Segal E, Chang HY: A transcriptional program mediating entry into cellular quiescence. PLoS Genet. 2007, 3: e91-10.1371/journal.pgen.0030091.
Schneider C, King RM, Philipson L: Genes specifically expressed at growth arrest of mammalian cells. Cell. 1988, 54: 787-793. 10.1016/S0092-8674(88)91065-3.
Teague TK, Hildeman D, Kedl RM, Mitchell T, Rees W, Schaefer BC, Bender J, Kappler J, Marrack P: Activation changes the spectrum but not the diversity of genes expressed by T cells. Proc Natl Acad Sci USA. 1999, 96: 12691-12696. 10.1073/pnas.96.22.12691.
Glynne R, Ghandour G, Rayner J, Mack DH, Goodnow CC: B-lymphocyte quiescence, tolerance and activation as viewed by global gene expression profiling on microarrays. Immunol Rev. 2000, 176: 216-246.
Venezia TA, Merchant AA, Ramos CA, Whitehouse NL, Young AS, Shaw CA, Goodell MA: Molecular signatures of proliferation and quiescence in hematopoietic stem cells. PLoS Biol. 2004, 2: e301-10.1371/journal.pbio.0020301.
Litovchick L, Florens LA, Swanson SK, Washburn MP, DeCaprio JA: DYRK1A protein kinase promotes quiescence and senescence through DREAM complex assembly. Genes Dev. 2011, 25: 801-813. 10.1101/gad.2034211.
Litovchick L, Sadasivam S, Florens L, Zhu X, Swanson SK, Velmurugan S, Chen R, Washburn MP, Liu XS, DeCaprio JA: Evolutionarily conserved multisubunit RBL2/p130 and E2F4 protein complex represses human cell cycle-dependent genes in quiescence. Mol Cell. 2007, 26: 539-551. 10.1016/j.molcel.2007.04.015.
Perna D, Faga G, Verrecchia A, Gorski MM, Barozzi I, Narang V, Khng J, Lim KC, Sung WK, Sanges R, Stupka E, Oskarsson T, Trumpp A, Wei CL, Muller H, Amati B: Genome-wide mapping of Myc binding and gene regulation in serum-stimulated fibroblasts. Oncogene. 2012, 31: 1695-1709. 10.1038/onc.2011.359.
Grandori C, Cowley SM, James LP, Eisenman RN: The Myc/Max/Mad network and the transcriptional control of cell behavior. Annu Rev Cell Dev Biol. 2000, 16: 653-699. 10.1146/annurev.cellbio.16.1.653.
Sang L, Coller HA, Roberts JM: Control of the reversibility of cellular quiescence by the transcriptional repressor HES1. Science. 2008, 321: 1095-1100. 10.1126/science.1155998.
Essers MA, Weijzen S, de Vries-Smits AM, Saarloos I, de Ruiter ND, Bos JL, Burgering BM: FOXO transcription factor activation by oxidative stress mediated by the small GTPase Ral and JNK. EMBO J. 2004, 23: 4802-4812. 10.1038/sj.emboj.7600476.
Skon CN, Jameson SC: Fox factors fight over T cell quiescence. Nat Immunol. 2011, 12: 522-524.
Filipowicz W, Bhattacharyya SN, Sonenberg N: Mechanisms of post-transcriptional regulation by microRNAs: are the answers in sight?. Nat Rev Genet. 2008, 9: 102-114.
Wu L, Fan J, Belasco JG: MicroRNAs direct rapid deadenylation of mRNA. Proc Natl Acad Sci USA. 2006, 103: 4034-4039. 10.1073/pnas.0510928103.
Flynt AS, Lai EC: Biological principles of microRNA-mediated regulation: shared themes amid diversity. Nat Rev Genet. 2008, 9: 831-842.
Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004, 116: 281-297. 10.1016/S0092-8674(04)00045-5.
Medina R, Zaidi SK, Liu CG, Stein JL, van Wijnen AJ, Croce CM, Stein GS: MicroRNAs 221 and 222 bypass quiescence and compromise cell survival. Cancer Res. 2008, 68: 2773-2780. 10.1158/0008-5472.CAN-07-6754.
Gillies JK, Lorimer IA: Regulation of p27Kip1 by miRNA 221/222 in glioblastoma. Cell Cycle. 2007, 6: 2005-2009. 10.4161/cc.6.16.4526.
le Sage C, Nagel R, Egan DA, Schrier M, Mesman E, Mangiola A, Anile C, Maira G, Mercatelli N, Ciafre SA, Farace MG, Agami R: Regulation of the p27(Kip1) tumor suppressor by miR-221 and miR-222 promotes cancer cell proliferation. EMBO J. 2007, 26: 3699-3708. 10.1038/sj.emboj.7601790.
Galardi S, Mercatelli N, Giorda E, Massalini S, Frajese GV, Ciafre SA, Farace MG: miR-221 and miR-222 expression affects the proliferation potential of human prostate carcinoma cell lines by targeting p27Kip1. J Biol Chem. 2007, 282: 23716-23724. 10.1074/jbc.M701805200.
Visone R, Russo L, Pallante P, De Martino I, Ferraro A, Leone V, Borbone E, Petrocca F, Alder H, Croce CM, Fusco A: MicroRNAs (miR)-221 and miR-222, both overexpressed in human thyroid papillary carcinomas, regulate p27Kip1 protein levels and cell cycle. Endocr Relat Cancer. 2007, 14: 791-798. 10.1677/ERC-07-0129.
Nimmo RA, Slack FJ: An elegant miRror: microRNAs in stem cells, developmental timing and cancer. Chromosoma. 2009, 118: 405-418. 10.1007/s00412-009-0210-z.
Moss EG: Heterochronic genes and the nature of developmental time. Curr Biol. 2007, 17: R425-434. 10.1016/j.cub.2007.03.043.
Yu F, Yao H, Zhu P, Zhang X, Pan Q, Gong C, Huang Y, Hu X, Su F, Lieberman J, Song E: let-7 regulates self renewal and tumorigenicity of breast cancer cells. Cell. 2007, 131: 1109-1123. 10.1016/j.cell.2007.10.054.
Caygill EE, Johnston LA: Temporal regulation of metamorphic processes in Drosophila by the let-7 and miR-125 heterochronic microRNAs. Curr Biol. 2008, 18: 943-950. 10.1016/j.cub.2008.06.020.
Reinhart BJ, Slack FJ, Basson M, Pasquinelli AE, Bettinger JC, Rougvie AE, Horvitz HR, Ruvkun G: The 21-nucleotide let-7 RNA regulates developmental timing in Caenorhabditis elegans. Nature. 2000, 403: 901-906. 10.1038/35002607.
Ambros V: MicroRNAs and developmental timing. Curr Opin Genet Dev. 2011, 21: 511-517. 10.1016/j.gde.2011.04.003.
Lee YS, Kim HK, Chung S, Kim KS, Dutta A: Depletion of human micro-RNA miR-125b reveals that it is critical for the proliferation of differentiated cells but not for the down-regulation of putative targets during differentiation. J Biol Chem. 2005, 280: 16635-16641. 10.1074/jbc.M412247200.
Schulman BR, Esquela-Kerscher A, Slack FJ: Reciprocal expression of lin-41 and the microRNAs let-7 and mir-125 during mouse embryogenesis. Dev Dyn. 2005, 234: 1046-1054. 10.1002/dvdy.20599.
Boyerinas B, Park SM, Hau A, Murmann AE, Peter ME: The role of let-7 in cell differentiation and cancer. Endocr Relat Cancer. 2010, 17: F19-36. 10.1677/ERC-09-0184.
Johnson CD, Esquela-Kerscher A, Stefani G, Byrom M, Kelnar K, Ovcharenko D, Wilson M, Wang X, Shelton J, Shingara J, Chin L, Brown D, Slack FJ: The let-7 microRNA represses cell proliferation pathways in human cells. Cancer Res. 2007, 67: 7713-7722. 10.1158/0008-5472.CAN-07-1083.
Legesse-Miller A, Elemento O, Pfau S, Forman J, Tavazoie S, Coller H: let-7 overexpression leads to an increased fraction of cells in G2/M, direct down-regulation of Cdc34 and stabilization of Wee1 kinase in primary fibroblasts. J Biol Chem. 2009, 284: 6605-6609.
Sokol NS, Xu P, Jan Y-N, Ambros V: Drosophila let-7 microRNA is required for remodeling of the neuromusculature during metamorphosis. Genes Dev. 2008, 22: 1591-1596. 10.1101/gad.1671708.
Esquela-Kerscher A, Trang P, Wiggins JF, Patrawala L, Cheng A, Ford L, Weidhaas JB, Brown D, Bader AG, Slack FJ: The let-7 microRNA reduces tumor growth in mouse models of lung cancer. Cell Cycle. 2008, 7: 759-764. 10.4161/cc.7.6.5834.
Johnson SM, Grosshans H, Shingara J, Byrom M, Jarvis R, Cheng A, Labourier E, Reinert KL, Brown D, Slack FJ: RAS is regulated by the let-7 microRNA family. Cell. 2005, 120: 635-647. 10.1016/j.cell.2005.01.014.
Kumar MS, Erkeland SJ, Pester RE, Chen CY, Ebert MS, Sharp PA, Jacks T: Suppression of non-small cell lung tumor development by the let-7 microRNA family. Proc Natl Acad Sci USA. 2008, 105: 3903-3908. 10.1073/pnas.0712321105.
Lee YS, Dutta A: The tumor suppressor microRNA let-7 represses the HMGA2 oncogene. Genes Dev. 2007, 21: 1025-1030. 10.1101/gad.1540407.
Takamizawa J, Konishi H, Yanagisawa K, Tomida S, Osada H, Endoh H, Harano T, Yatabe Y, Nagino M, Nimura Y, Mitsudomi T, Takahashi T: Reduced expression of the let-7 microRNAs in human lung cancers in association with shortened postoperative survival. Cancer Res. 2004, 64: 3753-3756. 10.1158/0008-5472.CAN-04-0637.
Sampson VB, Rong NH, Han J, Yang Q, Aris V, Soteropoulos P, Petrelli NJ, Dunn SP, Krueger LJ: MicroRNA let-7a down-regulates MYC and reverts MYC-induced growth in Burkitt lymphoma cells. Cancer Res. 2007, 67: 9762-9770. 10.1158/0008-5472.CAN-07-2462.
Trang P, Wiggins JF, Daige CL, Cho C, Omotola M, Brown D, Weidhaas JB, Bader AG, Slack FJ: Systemic delivery of tumor suppressor microRNA mimics using a neutral lipid emulsion inhibits lung tumors in mice. Mol Ther. 2011, 19: 1116-1122. 10.1038/mt.2011.48.
Trang P, Medina PP, Wiggins JF, Ruffino L, Kelnar K, Omotola M, Homer R, Brown D, Bader AG, Weidhaas JB, Slack FJ: Regression of murine lung tumors by the let-7 microRNA. Oncogene. 2010, 29: 1580-1587. 10.1038/onc.2009.445.
O'Donnell KA, Wentzel EA, Zeller KI, Dang CV, Mendell JT: c-Myc-regulated microRNAs modulate E2F1 expression. Nature. 2005, 435: 839-843. 10.1038/nature03677.
He L, Thomson JM, Hemann MT, Hernando-Monge E, Mu D, Goodson S, Powers S, Cordon-Cardo C, Lowe SW, Hannon GJ, Hammond SM: A microRNA polycistron as a potential human oncogene. Nature. 2005, 435: 828-833. 10.1038/nature03552.
Coller HA, Forman JJ, Legesse-Miller A: "Myc'ed messages": myc induces transcription of E2F1 while inhibiting its translation via a microRNA polycistron. PLoS Genet. 2007, 3: e146-10.1371/journal.pgen.0030146.
Sylvestre Y, De Guire V, Querido E, Mukhopadhyay UK, Bourdeau V, Major F, Ferbeyre G, Chartrand P: An E2F/miR-20a autoregulatory feedback loop. J Biol Chem. 2007, 282: 2135-2143.
Li Y, Zhang H, Chen Y: MicroRNA-mediated positive feedback loop and optimized bistable switch in a cancer network involving miR-17-92. PloS One. 2011, 6: e26302-10.1371/journal.pone.0026302.
Aguda BD, Kim Y, Piper-Hunter MG, Friedman A, Marsh CB: MicroRNA regulation of a cancer network: consequences of the feedback loops involving miR-17-92, E2F, and Myc. Proc Natl Acad Sci USA. 2008, 105: 19678-19683. 10.1073/pnas.0811166106.
Lemons JMS, Feng X-J, Bennett BD, Legesse-Miller A, Johnson EL, Raitman I, Pollina EA, Rabitz HA, Rabinowitz JD, Coller HA: Quiescent fibroblasts exhibit high metabolic activity. PLoS Biology. 2010, 8: e1000514-10.1371/journal.pbio.1000514.
Raymond CK, Roberts BS, Garrett-Engele P, Lim LP, Johnson JM: Simple, quantitative primer-extension PCR assay for direct monitoring of microRNAs and short-interfering RNAs. RNA. 2005, 11: 1737-1744. 10.1261/rna.2148705.
Pollina EA, Legesse-Miller A, Haley EM, Goodpaster T, Randolph-Habecker J, Coller HA: Regulating the angiogenic balance in tissues: a potential role for the proliferative state of fibroblasts. Cell Cycle. 2008, 7: 2056-2070. 10.4161/cc.7.13.6240.
Friedman RC, Farh KK-H, Burge CB, Bartel DP: Most mammalian mRNAs are conserved targets of microRNAs. Genome Res. 2009, 19: 92-105.
Grimson A, Farh KK-H, Johnston WK, Garrett-Engele P, Lim LP, Bartel DP: MicroRNA targeting specificity in mammals: determinants beyond seed pairing. Mol Cell. 2007, 27: 91-105. 10.1016/j.molcel.2007.06.017.
Jacobsen A, Wen J, Marks DS, Krogh A: Signatures of RNA binding proteins globally coupled to effective microRNA target sites. Genome Res. 2010, 20: 1010-1019. 10.1101/gr.103259.109.
Didiano D, Hobert O: Molecular architecture of a miRNA-regulated 3'-UTR. RNA. 2008, 14: 1297-1317. 10.1261/rna.1082708.
Nolde MJ, Saka N, Reinert KL, Slack FJ: The Caenorhabditis elegans pumilio homolog, puf-9, is required for the 3'UTR-mediated repression of the let-7 microRNA target gene, hbl-1. Dev Biol. 2007, 305: 551-563. 10.1016/j.ydbio.2007.02.040.
Engels B, Jannot G, Remenyi J, Simard MJ, Hutvagner G: Polypyrimidine tract binding protein (hnRNP I) is possibly a conserved modulator of miRNA-mediated gene regulation. PloS One. 2012, 7: e33144-10.1371/journal.pone.0033144.
Sandberg R, Neilson JR, Sarma A, Sharp PA, Burge CB: Proliferating cells express mRNAs with shortened 3' untranslated regions and fewer microRNA target sites. Science. 2008, 320: 1643-1647. 10.1126/science.1155390.
Mayr C, Bartel DP: Widespread shortening of 3'UTRs by alternative cleavage and polyadenylation activates oncogenes in cancer cells. Cell. 2009, 138: 673-684. 10.1016/j.cell.2009.06.016.
Ji Z, Lee JY, Pan Z, Jiang B, Tian B: Progressive lengthening of 3' untranslated regions of mRNAs by alternative polyadenylation during mouse embryonic development. Proc Natl Acad Sci USA. 2009, 106: 7028-7033. 10.1073/pnas.0900028106.
Arvey A, Larsson E, Sander C, Leslie CS, Marks DS: Target mRNA abundance dilutes microRNA and siRNA activity. Mol Syst Biol. 2010, 6: 363-
Smibert P, Miura P, Westholm JO, Shenker S, May G, Duff MO, Zhang D, Eads BD, Carlson J, Brown JB, Eisman RC, Andrews J, Kaufman T, Cherbas P, Celniker SE, Graveley BR, Lai EC: Global patterns of tissue-specific alternative polyadenylation in Drosophila. Cell Rep. 2012, 1: 277-289. 10.1016/j.celrep.2012.01.001.
Larsson E, Sander C, Marks D: mRNA turnover rate limits siRNA and microRNA efficacy. Mol Syst Biol. 2010, 6: 433-
Luna C, Li G, Qiu J, Epstein DL, Gonzalez P: Role of miR-29b on the regulation of the extracellular matrix in human trabecular meshwork cells under chronic oxidative stress. Mol Vis. 2009, 15: 2488-2497.
Maurer B, Stanczyk J, Jüngel A, Akhmetshina A, Trenkmann M, Brock M, Kowal-Bielecka O, Gay RE, Michel BA, Distler JH, Gay S, Distler O: miR-29 is a key regulator of collagen expression in systemic sclerosis. Arthritis Rheum. 2010, 62: 1733-1743.
Villarreal G, Oh DJ, Kang MH, Rhee DJ: Coordinated regulation of extracellular matrix synthesis by the microRNA-29 family in the trabecular meshwork. Invest Ophthalmol Vis Sci. 2011, 52: 3391-3397. 10.1167/iovs.10-6165.
Kriegel AJ, Liu Y, Fang Y, Ding X, Liang M: The miR-29 family: genomics, cell biology, and relevance to renal and cardiovascular injury. Physiol Genomics. 2012, 44: 237-244. 10.1152/physiolgenomics.00141.2011.
Merk DR, Chin JT, Dake BA, Maegdefessel L, Miller MO, Kimura N, Tsao PS, Iosef C, Berry GJ, Mohr FW, Spin JM, Alvira CM, Robbins RC, Fischbein MP: miR-29b participates in early aneurysm development in Marfan syndrome. Circ Res. 2012, 110: 312-324. 10.1161/CIRCRESAHA.111.253740.
Xiao J, Meng XM, Huang XR, Chung AC, Feng YL, Hui DS, Yu CM, Sung JJ, Lan HY: miR-29 inhibits bleomycin-induced pulmonary fibrosis in mice. Mol Ther. 2012, 20: 1251-1260. 10.1038/mt.2012.36.
Sporn MB, Roberts AB, Shull JH, Smith JM, Ward JM, Sodek J: Polypeptide transforming growth factors isolated from bovine sources and used for wound healing in vivo. Science. 1983, 219: 1329-1331. 10.1126/science.6572416.
Bandyopadhyay S, Friedman RC, Marquez RT, Keck K, Kong B, Icardi MS, Brown KE, Burge CB, Schmidt WN, Wang Y, McCaffrey AP: Hepatitis C virus infection and hepatic stellate cell activation downregulate miR-29: miR-29 overexpression reduces hepatitis C viral abundance in culture. J Infect Dis. 2011, 203: 1753-1762. 10.1093/infdis/jir186.
Luna C, Li G, Qiu J, Epstein DL, Gonzalez P: Cross-talk between miR-29 and transforming growth factor-betas in trabecular meshwork cells. Invest Ophthalmol Vis Sci. 2011, 52: 3567-3572. 10.1167/iovs.10-6448.
Ghosh AK, Yuan W, Mori Y, Varga J: Smad-dependent stimulation of type I collagen gene expression in human skin fibroblasts by TGF-beta involves functional cooperation with p300/CBP transcriptional coactivators. Oncogene. 2000, 19: 3546-3555. 10.1038/sj.onc.1203693.
Cutroneo KR: How is Type I procollagen synthesis regulated at the gene level during tissue fibrosis. J Cell Biochem. 2003, 90: 1-5. 10.1002/jcb.10599.
Salic A, Mitchison TJ: A chemical method for fast and sensitive detection of DNA synthesis in vivo. Proc Natl Acad Sci USA. 2008, 105: 2415-2420. 10.1073/pnas.0712168105.
Mollinari C, Reynaud C, Martineau-Thuillier S, Monier S, Kieffer S, Garin J, Andreassen PR, Boulet A, Goud B, Kleman JP, Margolis RL: The mammalian passenger protein TD-60 is an RCC1 family member with an essential role in prometaphase to metaphase progression. Dev Cell. 2003, 5: 295-307. 10.1016/S1534-5807(03)00205-3.
Hertel J, Bartschat S, Wintsche A, Otto C, Of The Bioinformatics Computer Lab TS, Stadler PF: Evolution of the let-7 microRNA Family. RNA Biol. 2012, 9:
Pasquinelli AE, Reinhart BJ, Slack F, Martindale MQ, Kuroda MI, Maller B, Hayward DC, Ball EE, Degnan B, Muller P, Spring J, Srinivasan A, Fishman M, Finnerty J, Corbo J, Levine M, Leahy P, Davidson E, Ruvkun G: Conservation of the sequence and temporal expression of let-7 heterochronic regulatory RNA. Nature. 2000, 408: 86-89. 10.1038/35040556.
Sempere LF, Dubrovsky EB, Dubrovskaya VA, Berger EM, Ambros V: The expression of the let-7 small regulatory RNA is controlled by ecdysone during metamorphosis in Drosophila melanogaster. Dev Biol. 2002, 244: 170-179. 10.1006/dbio.2002.0594.
Bashirullah A, Pasquinelli AE, Kiger AA, Perrimon N, Ruvkun G, Thummel CS: Coordinate regulation of small temporal RNAs at the onset of Drosophila metamorphosis. Dev Biol. 2003, 259: 1-8. 10.1016/S0012-1606(03)00063-0.
Hutvagner G, McLachlan J, Pasquinelli AE, Balint E, Tuschl T, Zamore PD: A cellular function for the RNA-interference enzyme Dicer in the maturation of the let-7 small temporal RNA. Science. 2001, 293: 834-838. 10.1126/science.1062961.
Lagos-Quintana M, Rauhut R, Yalcin A, Meyer J, Lendeckel W, Tuschl T: Identification of tissue-specific microRNAs from mouse. Curr Biol. 2002, 12: 735-739. 10.1016/S0960-9822(02)00809-6.
Sempere LF, Freemantle S, Pitha-Rowe I, Moss E, Dmitrovsky E, Ambros V: Expression profiling of mammalian microRNAs uncovers a subset of brain-expressed microRNAs with possible roles in murine and human neuronal differentiation. Genome Biol. 2004, 5: R13-10.1186/gb-2004-5-3-r13.
Lu J, Getz G, Miska EA, Alvarez-Saavedra E, Lamb J, Peck D, Sweet-Cordero A, Ebert BL, Mak RH, Ferrando AA, Downing JR, Jacks T, Horvitz HR, Golub TR: MicroRNA expression profiles classify human cancers. Nature. 2005, 435: 834-838. 10.1038/nature03702.
Raponi M, Dossey L, Jatkoe T, Wu X, Chen G, Fan H, Beer DG: MicroRNA classifiers for predicting prognosis of squamous cell lung cancer. Cancer Res. 2009, 69: 5776-5783. 10.1158/0008-5472.CAN-09-0587.
Peter ME: Let-7 and miR-200 microRNAs: guardians against pluripotency and cancer progression. Cell Cycle. 2009, 8: 843-852. 10.4161/cc.8.6.7907.
Mitra D, Das PM, Huynh FC, Jones FE: Jumonji/ARID1 B (JARID1B) protein promotes breast tumor cell cycle progression through epigenetic repression of microRNA let-7e. J Biol Chem. 2011, 286: 40531-40535. 10.1074/jbc.M111.304865.
Guan Y, Yao H, Zheng Z, Qiu G, Sun K: MiR-125b targets BCL3 and suppresses ovarian cancer proliferation. Int J Cancer. 2011, 128: 2274-2283. 10.1002/ijc.25575.
Zhang Y, Yan LX, Wu QN, Du ZM, Chen J, Liao DZ, Huang MY, Hou JH, Wu QL, Zeng MS, Huang WL, Zeng YX, Shao JY: miR-125b is methylated and functions as a tumor suppressor by regulating the ETS1 proto-oncogene in human invasive breast cancer. Cancer Res. 2011, 71: 3552-3562. 10.1158/0008-5472.CAN-10-2435.
Morita K, Han M: Multiple mechanisms are involved in regulating the expression of the developmental timing regulator lin-28 in Caenorhabditis elegans. EMBO J. 2006, 25: 5794-5804. 10.1038/sj.emboj.7601451.
Wu L, Belasco JG: Micro-RNA regulation of the mammalian lin-28 gene during neuronal differentiation of embryonal carcinoma cells. Mol Cell Biol. 2005, 25: 9198-9208. 10.1128/MCB.25.21.9198-9208.2005.
Li Z, Hassan MQ, Jafferji M, Aqeilan RI, Garzon R, Croce CM, van Wijnen AJ, Stein JL, Stein GS, Lian JB: Biological functions of miR-29b contribute to positive regulation of osteoblast differentiation. J Biol Chem. 2009, 284: 15676-15684. 10.1074/jbc.M809787200.
Sengupta S, den Boon JA, Chen IH, Newton MA, Stanhope SA, Cheng YJ, Chen CJ, Hildesheim A, Sugden B, Ahlquist P: MicroRNA 29c is down-regulated in nasopharyngeal carcinomas, up-regulating mRNAs encoding extracellular matrix proteins. Proc Natl Acad Sci USA. 2008, 105: 5874-5878. 10.1073/pnas.0801130105.
Steele R, Mott JL, Ray RB: MBP-1 upregulates miR-29b that represses Mcl-1, collagens, and matrix-metalloproteinase-2 in prostate cancer cells. Genes Cancer. 2010, 1: 381-387. 10.1177/1947601910371978.
Renault VM, Rafalski VA, Morgan AA, Salih DA, Brett JO, Webb AE, Villeda SA, Thekkat PU, Guillerey C, Denko NC, Palmer TD, Butte AJ, Brunet A: FoxO3 regulates neural stem cell homeostasis. Cell Stem Cell. 2009, 5: 527-539. 10.1016/j.stem.2009.09.014.
Miyamoto K, Araki KY, Naka K, Arai F, Takubo K, Yamazaki S, Matsuoka S, Miyamoto T, Ito K, Ohmura M, Chen C, Hosokawa K, Nakauchi H, Nakayama K, Nakayama KI, Harada M, Motoyama N, Suda T, Hirao A: Foxo3a is essential for maintenance of the hematopoietic stem cell pool. Cell Stem Cell. 2007, 1: 101-112. 10.1016/j.stem.2007.02.001.
Tothova Z, Kollipara R, Huntly BJ, Lee BH, Castrillon DH, Cullen DE, McDowell EP, Lazo-Kallanian S, Williams IR, Sears C, Armstrong SA, Passegue E, DePinho RA, Gilliland DG: FoxOs are critical mediators of hematopoietic stem cell resistance to physiologic oxidative stress. Cell. 2007, 128: 325-339. 10.1016/j.cell.2007.01.003.
Feng X, Ippolito GC, Tian L, Wiehagen K, Oh S, Sambandam A, Willen J, Bunte RM, Maika SD, Harriss JV, Caton AJ, Bhandoola A, Tucker PW, Hu H: Foxp1 is an essential transcriptional regulator for the generation of quiescent naive T cells during thymocyte development. Blood. 2010, 115: 510-518. 10.1182/blood-2009-07-232694.
Feng X, Wang H, Takata H, Day TJ, Willen J, Hu H: Transcription factor Foxp1 exerts essential cell-intrinsic regulation of the quiescence of naive T cells. Nat Immunol. 2011, 12: 544-550. 10.1038/ni.2034.
Mott JL, Kurita S, Cazanave SC, Bronk SF, Werneburg NW, Fernandez-Zapico ME: Transcriptional suppression of mir-29b-1/mir-29a promoter by c-Myc, hedgehog, and NF-kappaB. J Cell Biochem. 2010, 110: 1155-1164. 10.1002/jcb.22630.
Fabbri M, Garzon R, Cimmino A, Liu Z, Zanesi N, Callegari E, Liu S, Alder H, Costinean S, Fernandez-Cymering C, Volinia S, Guler G, Morrison CD, Chan KK, Marcucci G, Calin GA, Huebner K, Croce CM: MicroRNA-29 family reverts aberrant methylation in lung cancer by targeting DNA methyltransferases 3A and 3B. Proc Natl Acad Sci USA. 2007, 104: 15805-15810. 10.1073/pnas.0707628104.
Mott JL, Kobayashi S, Bronk SF, Gores GJ: mir-29 regulates Mcl-1 protein expression and apoptosis. Oncogene. 2007, 26: 6133-6140. 10.1038/sj.onc.1210436.
Martinez I, Cazalla D, Almstead LL, Steitz JA, DiMaio D: miR-29 and miR-30 regulate B-Myb expression during cellular senescence. Proc Natl Acad Sci USA. 2011, 108: 522-527. 10.1073/pnas.1017346108.
Santanam U, Zanesi N, Efanov A, Costinean S, Palamarchuk A, Hagan JP, Volinia S, Alder H, Rassenti L, Kipps T, Croce CM, Pekarsky Y: Chronic lymphocytic leukemia modeled in mouse by targeted miR-29 expression. Proc Natl Acad Sci USA. 2010, 107: 12210-12215. 10.1073/pnas.1007186107.
King TE, Pardo A, Selman M: Idiopathic pulmonary fibrosis. Lancet. 2011, 378: 1949-1961. 10.1016/S0140-6736(11)60052-4.
Gharaee-Kermani M, Hu B, Phan SH, Gyetko MR: Recent advances in molecular targets and treatment of idiopathic pulmonary fibrosis: focus on TGFbeta signaling and the myofibroblast. Curr Med Chem. 2009, 16: 1400-1417. 10.2174/092986709787846497.
Leek JT, Storey JD: Capturing heterogeneity in gene expression studies by surrogate variable analysis. PLoS Genetics. 2007, 3: 1724-1735.
Generic Gene Ontology Term Finder. [http://go.princeton.edu/cgi-bin/GOTermFinder]
TargetScan 5.1. [http://www.targetscan.org]
Minka TP: Estimating a Dirichlet distribution. 2000, Microsoft Research
Edgar R, Domrachev M, Lash AE: Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002, 30: 207-210. 10.1093/nar/30.1.207.
HC is the Milton E. Cassel scholar of the Rita Allen Foundation. EJS and EJ are supported in part by a National Science Foundation Graduate Research Fellowship DGE-0646086. HC and ALM are supported by the NIGMS Center of Excellence grant P50 GM071508. ALM acknowledges support from National Cancer Institute K01CA128887. EJ acknowledges support from NIH Training Grant 2T32 CA009528. JF acknowledges support from NCI training grant 5T32 CA009528. TC and MK acknowledge support from the Howard Hughes Medical Institute/Princeton Summer Undergraduate Research Program. JL acknowledges support from NIH training grant T32 HG003284. This work was funded by PhRMA Foundation grant 2007RSGl9572, NIH/NIGMS 1R01 GM081686, and NIH/NIGMS 1R01 GM086465.
We wish to acknowledge Sarah Pfau (MIT), Cheng Shi (Princeton University), Liling Wang (Princeton University), Christina deCoste (Princeton University), Nithya Krishnan (Princeton University), Irene Raitman (Princeton University), Rosetta Inpharmatics, and all of the members of the Coller lab for helpful discussions. PUMAdb is funded in part by the National Institute of General Medical Sciences (NIGMS) (NIH grant P50 GM071508).
The authors declare that they have no competing interests.
EJS, ALM, and JL conducted microarray experiments. EJS performed the statistical analyses and biochemical studies. EJ, MK, ALM, EJS, and MR conducted the molecular biology assays. ESS, TC, EJS, and MR conducted the cell cycle assays. EJS, ALM, JF, and HC conceived of the study, participated in its design and coordination, and helped to draft the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
About this article
Cite this article
Suh, E.J., Remillard, M.Y., Legesse-Miller, A. et al. A microRNA network regulates proliferative timing and extracellular matrix synthesis during cellular quiescence in fibroblasts. Genome Biol 13, R121 (2012). https://doi.org/10.1186/gb-2012-13-12-r121
- Cell cycle
- Extracellular matrix