Integrated miRNA, mRNA and protein expression analysis reveals the role of post-transcriptional regulation in controlling CHO cell growth rate
- Colin Clarke†1Email author,
- Michael Henry†1,
- Padraig Doolan1,
- Shane Kelly1,
- Sinead Aherne1, 2,
- Noelia Sanchez1,
- Paul Kelly1,
- Paula Kinsella1,
- Laura Breen1,
- Stephen F Madden1, 2,
- Lin Zhang3,
- Mark Leonard3,
- Martin Clynes1,
- Paula Meleady†1Email author and
- Niall Barron†1Email author
© Clarke et al.; licensee BioMed Central Ltd. 2012
Received: 3 August 2012
Accepted: 9 November 2012
Published: 21 November 2012
To study the role of microRNA (miRNA) in the regulation of Chinese hamster ovary (CHO) cell growth, qPCR, microarray and quantitative LC-MS/MS analysis were utilised for simultaneous expression profiling of miRNA, mRNA and protein. The sample set under investigation consisted of clones with variable cellular growth rates derived from the same population. In addition to providing a systems level perspective on cell growth, the integration of multiple profiling datasets can facilitate the identification of non-seed miRNA targets, complement computational prediction tools and reduce false positive and false negative rates.
51 miRNAs were associated with increased growth rate (35 miRNAs upregulated and 16 miRNAs downregulated). Gene ontology (GO) analysis of genes (n=432) and proteins (n=285) found to be differentially expressed (DE) identified biological processes driving proliferation including mRNA processing and translation. To investigate the influence of miRNA on these processes we combined the proteomic and transcriptomic data into two groups. The first set contained candidates where evidence of translational repression was observed (n=158). The second group was a mixture of proteins and mRNAs where evidence of translational repression was less clear (n=515). The TargetScan algorithm was utilised to predict potential targets within these two groups for anti-correlated DE miRNAs.
The evidence presented in this study indicates that biological processes such as mRNA processing and protein synthesis are correlated with growth rate in CHO cells. Through the integration of expression data from multiple levels of the biological system a number of proteins central to these processes including several hnRNPs and components of the ribosome were found to be post-transcriptionally regulated. We utilised the expression data in conjunction with in-silico tools to identify potential miRNA-mediated regulation of mRNA/proteins involved in CHO cell growth rate. These data have allowed us to prioritise candidates for cell engineering and/or biomarkers relevant to industrial cell culture. We also expect the knowledge gained from this study to be applicable to other fields investigating the role of miRNAs in mammalian cell growth.
KeywordsChinese hamster ovary CHO cells Growth rate MicroRNA mRNA Microarray Proteomics
Our understanding of the role microRNAs play in fundamental biological processes in both plants and animals has increased dramatically over the last decade . Since the discovery of miRNA in C.elegans, the miRBase data repository has expanded to hold sequence data from over 21,000 mature miRNAs across 168 species . These short, highly conserved RNA molecules (~22 nucleotides) form a layer of post-transcriptional control of gene expression, generally repressing translation  (via translational inhibition, transcriptional degradation and in some instances mRNA deadenylation ) or in rare cases actually enhancing translation . Thus far, miRNAs have been implicated in a broad range of processes from cell cycle control  to apoptosis . In addition, the effect of miRNAs on diseases such as cancer  and diabetes  has been intensively studied.
The complexity of miRNA target recognition remains a significant challenge to researchers. For instance, a single miRNA is estimated to target an average of 100–200 mRNAs  while a single mRNA transcript can be targeted by hundreds of miRNAs  and multiple miRNAs can cooperatively repress a range of targets . The laborious experimental techniques required to confirm interaction between a miRNA and mRNA have necessitated the use of in silico target prediction to prioritise targets for wet-lab confirmation and to determine the potential function of a miRNA. The most widely applied algorithms including miRanda , Pictar  and TargetScan  use a combination of sequence complementarity of a transcript to a conserved region at the 5’ end of the miRNA spanning position 2 to position 7 (known as the “seed” region), thermodynamic feasibility of hybridisation and evolutionary conservation . Until recently, active animal miRNA recognition sites were thought to be present solely within the 3’ UTR of an mRNA. However, recent evidence has suggested that sites also exist in the 5’ UTR , within the coding sequence , and in some cases these sites may be present in multiple locations within a transcript. Target prediction algorithms are undoubtedly valuable tools to researchers in the field providing a rapid appreciation of the potential processes impacted by a particular miRNA and prioritising potential direct targets for validation assays. Quantitative evaluation of algorithm performance has thus far proved difficult due to the limited number of experimentally confirmed targets. Previous studies have yielded less than encouraging false positive and false negative rates resulting from the use of algorithms [17, 18] prompting an increasing focus on combining multiple expression profiling datasets with in silico target prediction .
The Chinese hamster ovary (CHO) cell has been utilised extensively for the last 20 years in the biopharmaceutical industry and has become the cell type of choice for production of recombinant proteins for medical applications due to safety of use, rapid growth characteristics and the ability to secrete large quantities of correctly folded proteins. The majority of industrial advances have thus far stemmed from improvements in cell culture media, vectors and the design of bioreactors . Several research groups have in recent years been interested in manipulating CHO at the molecular level to improve protein production efficiency and create diagnostic tools to monitor manufacturing processes. The recent publication of the CHO genome  along with sequence from similar initiatives [22, 23] promises to facilitate analyses at all levels of the CHO biological system not only through direct analysis of genomic sequence but also through the improvement of analytical platforms such as microarrays and mass spectrometry.
Since the first report of altered miRNA expression in CHO as a result of modifying bioprocess conditions , the number of publications in this area has increased steadily. The attraction of miRNA based cell engineering arises from the potential of miRNAs to alter an entire pathway or indeed pathways, to enhance industrially beneficial phenotypes. Various studies have focussed on miRNA sequence analysis, determining homology to other species and location of genomic loci via next generation sequencing technologies [25–27]. To date, miRNAs have been associated with several important bioproduction phenotypes including growth rate , productivity  and apoptosis . In this study, we elucidate new and expand existing knowledge on the contribution of miRNA-mediated regulation to CHO cellular growth rate. Our previous work has shown that building CHO cell density to a high level in the bioreactor is intimately linked to the final volumetric titre of product and can in some cases be more important than the intrinsic productivity rate .
The experimental design and stringently controlled panel of samples used here ensures that the data presented in this study are applicable to the investigation of mammalian cell growth beyond the bioprocessing field. Firstly, the cells under investigation were selected from the same clonal population following a process of repeated passaging and monitoring of growth rate in order to minimise non-growth related variation. A set of these sister clones spanning a continuous range of low to high growth rates was selected for analysis. The second crucial aspect of our approach was to measure expression of the transcriptome (mRNA & miRNA) and proteome in parallel. While the number of studies combining miRNA, protein and mRNA expression data are limited in comparison to those comprising of two data types (i.e. miRNA & mRNA or miRNA & protein), it is evident that the use of all of 3 of these data types can facilitate the identification of direct miRNA targets  and enhance our understanding of the biological role of miRNAs . Similar to those studies, we analyse the expression data in conjunction with miRNA target prediction algorithms in order to decrease false positives, reduce the effects of experimental noise and ultimately to increase the likelihood of finding direct miRNA targets.
It is intended that this work will contribute to the understanding of those biological processes driving mammalian cell growth as well as supporting other researchers in the selection of candidates for miRNA-target confirmation assays.
The growth rates of the final sample set spanned a range of 0.011 to 0.044 hr-1 with mean productivity = 24 (± 3) pg protein/cell/day. Consistent behaviour in the samples subjected to expression profiling in terms of growth rate and productivity was ensured by monitoring over 40 passages. By choosing sister clones derived from the same transfection pool with similar recombinant protein production rates and differing only in growth rate we sought to eliminate noise and expose those variations related to the proliferation phenotype. To prioritise miRNAs, mRNAs and proteins associated with cell growth, we separated the dataset into 15 “fast” (≥ 0.025 h-1) and 15 “slow” (≤ 0.023 h-1) samples. Note: 3 outlying “slow” samples (biological replicates) within the LC-MS/MS dataset were removed following quality control using principal components analysis (data not shown). To equalise the sample numbers on both sides of the proteomic differential expression analysis a single fast growing clone (3 biological replicates) was selected at random and removed leaving 12 fast versus 12 slow samples (Additional file 1).
Prior to target prediction against the DE miRNA list, we separated protein and mRNA targets into two groups. “Group A” contains those targets where a degree of post-transcriptional regulation was observed (possibly via miRNA mediated translational repression). The proteins in Group A were DE between the fast and slow clones, their respective mRNAs were expressed above the microarray detection threshold but no change in mRNA expression was observed. The second group of candidates, referred to as “Group B” were comprised of (1) DE proteins where the probeset was under the detection threshold or not present on the chip, (2) DE mRNAs where the corresponding protein was not identified within the fraction analysed by LC-MS/MS and (3) candidates where differential expression was observed at both the mRNA and protein level (Figure 1). It is therefore likely that Group B contains a mixture of targets that could be undergoing miRNA translational repression, transcript destabilisation or indeed alteration due to non-miRNA-mediated processes. Group B candidates were considered of lower priority than Group A not because we expected a lower proportion of predicted miRNA targets, but because of incomplete evidence at the protein and mRNA levels or both the protein and mRNA were DE. The final stage in analysis involved in silico target prediction with TargetScan 6.1  for Groups A and B. Each protein/mRNA from the two groupings was compared with the predicted targets of anti-correlated DE miRNAs.
miRNA expression levels are associated with variations in the rate of CHO cell growth
Of those miRNAs found to be DE in this study, a number have been associated with cellular growth previously [34–40]. Several of these upregulated miRNAs form part of the miR-17-92 cluster (e.g. miR17, miR20a, miR20b, miR18a, miR18b and miR106a), a well studied group of miRNAs linked to cancer . It has been shown that c-myc  and N-myc  directly activate the mir-17-92 cluster, while in contrast repression of the cluster by p53 has been demonstrated . CHO is a highly proliferative cell line in general having been cultured in vitro for over 50 years  and previous studies have shown that the miR-17-92 cluster is highly expressed . The data presented in this study indicates a relationship between increased miR-17-92 cluster member expression levels and rapidly proliferating CHO cells. Several miRNAs, including miR-204, miR-338-3p, miR-497, miR-30e and miR-206 are downregulated as growth rate increases in agreement with independent studies in other systems [38, 45–47]. miR-451 (downregulated at higher growth rates in this study) has been shown to inhibit growth and induce apoptosis [37, 48]. Godlewski et al. demonstrated that in glioma cells miR-451 expression is correlated to glucose levels and upon glucose depletion a downregulation of miR-451 is observed along with slower proliferation and increased survival . In contrast, downregulation of miR-451 in CHO is observed here at higher growth rates, exemplifying the previously recognised cell-specific nature of miRNA expression.
Proteomic and transcriptomic profiling reveal several overrepresented biological processes related to mRNA processing and translation
285 non-redundant proteins (180 upregulated and 105 downregulated) were found to be DE between fast and slow growing CHO cell clones (Additional file 3). Following differential expression analysis of the mRNA data we identified 432 DE non-redundant annotated transcripts (229 upregulated and 203 downregulated) (Additional file 4). 44 proteins were also dysregulated at the mRNA level (29 upregulated and 15 down regulated, corresponding to 50 DE probesets (Figure 2C)). Correlation between both datasets was observed for 11 of the 44 proteins in this group (9 upregulated and 2 downregulated).
Protein enrichment analysis
generation of precursor metabolites and energy
coenzyme metabolic process
cell redox homeostasis
mRNA enrichment analysis
mitotic cell cycle
mRNA metabolic process
cellular macromolecular complex subunit organization
mRNA processing and splicing processes were also found overrepresented within both the protein and transcript differential expression lists. Alternative splicing is understood to represent an important stage at which regulation of translation can be mediated. Furthermore, RNA splicing is now recognised as a central step in gene expression whereby virtually all precursor mRNAs (pre-mRNAs) undergo alternative splicing, resulting in a complex level of expression regulation . Specific splicing factors are known to be important for cell cycle control, for instance multiple splice variants of the p53 tumour suppressor are DE in cancer . We observed the post-transcriptional upregulation of YBX1/YB1 in rapidly growing CHO cells; this transcription factor is associated with cancer susceptibility  and known to play a role in splice site selection . Previous studies reporting an increase or decrease in tumour cell growth upon over or under-expression of this protein respectively  are in-agreement with the evidence reported here.
In recent years, various studies have demonstrated the role of RNA-binding proteins (RBPs) in the regulation of groups of transcripts by shuttling them efficiently through processes such as mRNA splicing, transport and ultimately translation . This post-transcriptional regulation of multiple mRNAs, termed “RNA regulons”, is thought to allow the cell to respond rapidly to environmental changes. One class of RBPs, known as hnRNPs, interact with large numbers of pre-mRNAs to form hnRNP-RNA complexes containing combinations of at least 20 hnRNPs. Further examination of those RNA processing and splicing GO categories found to be enriched in this study revealed the presence of multiple hnRNPs DE at both the protein (e.g. hnRNPM, hnRNPC, hnRNPAB, and hnRNPK) and mRNA (e.g. hnRNPD and hnRNPR) level (Additional files 3 and 4).
Separation of DE mRNA and DE proteins into target Groups A and B
To compare fluctuations in protein abundance with transcript expression we identified probesets on the array representing the DE proteins identified by LC-MS/MS. In the case of multiple probesets targeting an individual protein, the probeset with the highest expression was selected to represent that protein. 554 probesets corresponding to 260 proteins were identified (probesets for 25 proteins were absent from the microarray) (Figure 2B). Group A contained 196 probesets corresponding to 158 unique DE proteins (Figure 2C) where a degree of post-transcriptional regulation was observed. Group B was comprised of 515 candidates including DE mRNAs, DE proteins where no probeset was available or the mRNA was undetected and candidates DE at both the mRNA and protein level. For the 44 Group B targets where both the mRNA and protein was DE only the protein data (i.e. direction) was utilised for miRNA target overlap.
In-silico analysis identifies a number of transcriptomic and proteomic targets of priority miRNAs
We would not expect all candidates to be predicted targets of miRNAs for a number of reasons. It is likely differential expression of mRNA and proteins is occurring due to a variety of non-miRNA related processes. Furthermore the dysregulation of proteins such as hnRNPs suggests the likelihood of non-miRNA post-transcriptional regulation. The design of the TargetScan algorithm, considered to be amongst the most stringent methods, limits the numbers of in-silico predicted targets. The 15 “star” miRNAs found to be dysregulated and correlated to growth rate in this study could not be predicted as the TargetScan database does not contain entries for these miRNAs. In addition the algorithm only considers target sites within the 3’ UTR and therefore excludes potential target sites within the coding region or at the 5’ end of a transcript .
To investigate the impact of miRNA, mRNA and protein expression on cellular growth rate we have employed an integrative methodology to combine data from several global profiling technologies with bioinformatics analysis. We sought to minimise variation through the use of carefully selected Mab-secreting CHO cell lines that, in spite of spanning a wide range of growth rates, were derived from the same transfection pool. Furthermore the abundance of miRNA, mRNA and protein was determined in parallel on identical samples to further reduce biological noise. The experimental approach may be applicable to the study of cell growth in other eukaryotic systems and prove useful in elucidating mechanisms of cellular proliferation in general.
The most crucial aspect of the experimental design is the combination of data from multiple expression profiling methods, genomic sequence and in silico prediction to study miRNA function. Both proteomic mass spectrometry and mRNA microarrays have been used previously to study miRNA function, however both methods when used in isolation, suffer from several disadvantages as noted by previous researchers in the area [32, 69]. For instance quantitative mass spectrometry based proteomics yields in the order of hundreds of DE proteins and depending on the prefractionation method may not detect many of the low abundant proteins or integral membrane proteins. On the other hand gene expression analysis using microarrays provides a wide coverage of mRNAs but post-transcriptional processes may be missed. In terms of the study of the role of miRNAs both methods when used in isolation rely heavily on computational methods to predict miRNA interaction and prioritise potential direct targets. The availability of data on multiple levels of the biological system allows us to identify targets that would not have been identified by a single dataset. Moreover, prioritisation of potential targets undergoing classical miRNA-mediated translation repression can only be achieved through the integration of both the mRNA and protein datasets.
The availability of a combined profiling approach could reduce the false negative or false positive rates associated with in silico prediction as well as enriching priority candidate cohorts for functional validation. For example, a recent study demonstrated the enhancement of ribosomal subunit translation and ribosome biogenesis upon miR-10a binding via a “non-seed” site at the transcript 5’ UTR  and confirmed that ribosome formation can be modulated by miRNA to some extent. In this study miR-10a was both upregulated (FC = 2.12) and positively correlated (PCC = 0.53) with cellular growth rate, and several ribosomal proteins previously identified as mir-10a targets showed some evidence of post-transcriptional regulation. We found no predictive evidence of miR-10a interaction with those RPs using TargetScan as the algorithm searches for the presence of sites conforming to classical miRNA seed based rules within 3’ UTR. In addition, RP 3’ UTRs tend to be relatively short.
Rab14 is an experimentally validated direct target of miR-451 and has been shown to activate tumour suppression . miR451 was downregulated and the Rab14 protein upregulated (while mRNA expression remained constant) in fast growing CHO cells. The Rab14-miR451 3’ UTR binding site is poorly conserved and therefore not ranked highly by TargetScan; in cases such as this the availability of evidence from multiple profiling datasets could be used to prioritise poorly conserved miRNA interactions.
In summary, we have analysed a range of production CHO cell clones to investigate the role of miRNAs in growth rate variation. The antibody producing clones under analysis were chosen to control for several possible confounding factors. We employed global miRNA, mRNA and proteomic expression profiling in parallel in order to integrate three levels of the biological system. Analysis of the datasets in isolation provided information on the processes involved in cell growth which included ribosome biogenesis, translation and mRNA processing, which were to some degree regulated post-transcriptionally. The ability to compare mRNA and proteins in identical samples with respect to miRNA allowed us to identify potential miRNA targets, and highlight translational repression targets which could not have been identified using a single dataset. Moreover, the use of multiple profiling datasets could permit the identification of non-seed miRNA targets complementing computational prediction tools and reducing the false positive and false negative rates. While this study is particularly relevant for the bio-pharmaceutical industry by prioritising a number of potential miRNA cell engineering candidates, the experimental design ensures that the data generated and knowledge gained on biological processes driving cellular growth are applicable to other mammalian systems.
In total, 30 clones were selected for expression profiling ranging from 0.011 to 0.044 hr-1. 15 clones were designated as “fast” (≥ 0.025 h-1) and 15 clones designated as “slow” (≤ 0.023 h-1) (Additional file 1).
qPCR screening of miRNA expression
Total RNA was extracted from 106 cells using TrireagentTM (Sigma-Aldrich), resuspended in nuclease free water, quantified on a NanodropTM spectrophotometer and checked for integrity on a Eukaryote Total RNA Nano Bioanalyzer chip (Agilent). TaqMan Low-Density Array cards (TLDAs) (Human MicroRNA A&B Cards V2.0) were run as per the manufacturer’s guidelines (Applied BioSystems). Each card consists of 384-wells containing primers designed against individual miRNAs. 100 ng of total RNA was reverse transcribed in 2 individual multiplex reactions. These cDNA mixes were subjected to 12 cycles of pre-amplification with pre-amp primer pools and then used to load the TLDA card. PCR was performed on an AB7900 real time instrument for 10 min at 95°C followed by 40 cycles of 30 sec at 97°C and 1 min at 60°C.
The R statistical software environment (http://www.r-project.org) and the HTqPCR bioconductor package were utilised for data analysis . Mamm-U6 expression was used to normalize across the samples and differential expression was calculated using the 2-ΔΔCt method. miRNAs with a fold change ≥ 1.2 in either direction between the fast and slow groups with a Benjamini-Hochberg (BH) adjusted p-value of < 0.05 were considered significant. The Pearson correlation coefficient (PCC) between miRNA ΔCT and sample growth rate (Additional file 1) was also utilised as an additional filter.
Gene expression analysis was carried out on the proprietary CHO-specific WyeHamster3a oligonucleotide microarray. The array contains 19,809 probesets corresponding to 132 control sequences, 11,857 probesets annotated to mouse, rat and/or human Unigene IDs (9,098 non-redundant genes) and 7,820 unannotated probesets. The methodology and criteria used for total RNA purification, cRNA sample processing and hybridisation to hamster microarrays have been previously described . All microarray data were pre-processed as described previously . Prior to data-analysis probesets that did not reach the detection threshold (fluorescence level ≥ log2 (100) for at least 1 sample) were identified and designated undetected. The remaining probesets were considered differentially expressed between the fast and slow groups if a fold change ≥ 1.2 in either direction along with a BH adjusted p-value < 0.05 was observed. The microarray data used in this study have been deposited in the NCBI GEO database (GSE37251).
Proteomics sample preparation, LC-MS/MS and data analysis
Sample preparation and label-free LC-MS was carried out as previously described . Data analysis was performed using Progenesis label-free LC-MS software version 3.1 (NonLinear Dynamics LTD, Newcastle upon Tyne, UK) as recommended by the manufacturer (see http://www.nonlinear.com for further background information regarding alignment, normalisation, calculation of peptide abundance, etc.). Briefly, the raw MS data is processed as follows; a run is selected that is representative of the data, to which the LC retention times of all the other samples within the experiment are aligned. The Progenesis peptide quantification algorithm calculates normalised peptide abundances as the sum of the peak areas within each peptide isotope boundary. Protein abundance is calculated from the sum of all unique peptide abundances for an individual protein on each run.
A number of criteria were used to filter the data before exporting the MS/MS output files to MASCOT (http://www.matriscience.com) for protein identification; only peptide features with a p-value < 0.05 (determined via an ANOVA) between experimental groups, mass peaks (features) with charge states from +1 to +3, and greater than 3 isotopes per peptide were retained. All MS/MS spectra were exported from Progenesis software as a MASCOT generic file (.mgf) for peptide identification with MASCOT (version 2.2) and searched against the BB-CHO specific database . The search parameters used were as follows: peptide mass tolerance set to 20 ppm, MS/MS mass tolerance set at 0.5 Da; up to two missed cleavages were allowed, carbamidomethylation set as a fixed modification and methionine oxidation set as a variable modification. Peptides with ion scores of 30 and above were re-imported into the Progenesis LC–MS software for further analysis. Only proteins with ≥ 2 peptides matched, a ≥ 1.2 fold difference in abundance in both directions and a p-value < 0.05 were considered to be DE.
Coinertia analysis is a multivariate statistical method utilised to compare datasets with different measurement sources on the same objects/samples. There have been several examples of CIA applied to “omics” data to date including a comparison of transcriptomic and proteomic data  and prediction of miRNA interactions from gene expression analysis . Briefly, CIA attempts to locate the axes of maximal co-variance between the proteomics and transcriptomic data from parallel samples. CIA was carried out using the MADE4 R package . Gene expression and proteomic data were log2 scaled and mean centred prior to analysis. In this study we employ CIA to visualise the disparity between transcript and protein abundance across the dataset in an unsupervised manner negating the requirement for arbitrary thresholds. The input for the CIA analysis was two matrices equal to n × p, where n = sample number and p = number proteins/mRNA potentially targeted by the miRNAs of interest. Note: Only 24 matched mRNA protein samples were used for CIA due to the presence of outlying samples in the LC-MS/MS dataset. Following CIA the normalised scores were plotted, each target is represented by an arrow with the circular base corresponding to the mRNA and the arrowhead corresponding to the protein. The length of the arrow relates to the difference between mRNA and protein expression across the dataset.
miRNA target prediction
Prediction of miRNA and oppositely correlated protein/mRNA interactions was performed using TargetScan 6.1 (http://www.targetscan.org/vert_61/) . Only conserved targets were utilised and each predicted target assigned a rank according to the TargetScan algorithm quality measure known as the total context+ score.
GO biological process enrichment analysis was carried out for the DE protein and mRNA lists via the DAVID interface (david.abcc.ncifcrf.gov).
Chinese hamster ovary
Liquid chromatography mass spectrometry
Pearson correlation coefficient
RNA binding protein
gene expression omnibus
Analysis of variance.
This work was supported by funding from Science Foundation Ireland (SFI) [grant number 07/IN.1/B1323].
- Bartel DP: MicroRNAs: target recognition and regulatory functions. Cell. 2009, 136 (2): 215-233. 10.1016/j.cell.2009.01.002.PubMed CentralView ArticlePubMedGoogle Scholar
- Lee RC, Feinbaum RL, Ambros V: The C. elegans heterochronic gene lin-4 encodes small RNAs with antisense complementarity to lin-14. Cell. 1993, 75 (5): 843-854. 10.1016/0092-8674(93)90529-Y.View ArticlePubMedGoogle Scholar
- Kozomara A, Griffiths-Jones S: miRBase: integrating microRNA annotation and deep-sequencing data. Nucleic Acids Res. 2011, 39 (Database issue): D152-D157.PubMed CentralView ArticlePubMedGoogle Scholar
- Wu L, Fan J, Belasco JG: MicroRNAs direct rapid deadenylation of mRNA. Proc Natl Acad Sci U S A. 2006, 103 (11): 4034-4039. 10.1073/pnas.0510928103.PubMed CentralView ArticlePubMedGoogle Scholar
- Orom UA, Nielsen FC, Lund AH: MicroRNA-10a binds the 5'UTR of ribosomal protein mRNAs and enhances their translation. Mol Cell. 2008, 30 (4): 460-471. 10.1016/j.molcel.2008.05.001.View ArticlePubMedGoogle Scholar
- Chivukula RR, Mendell JT: Circular reasoning: microRNAs and cell-cycle control. Trends Biochem Sci. 2008, 33 (10): 474-481. 10.1016/j.tibs.2008.06.008.PubMed CentralView ArticlePubMedGoogle Scholar
- Jovanovic M, Hengartner MO: miRNAs and apoptosis: RNAs to die for. Oncogene. 2006, 25 (46): 6176-6187. 10.1038/sj.onc.1209912.View ArticlePubMedGoogle Scholar
- Iorio MV, Croce CM: MicroRNA dysregulation in cancer: diagnostics, monitoring and therapeutics. A comprehensive review. EMBO Mol Med. 2012, 4 (3): 143-159. 10.1002/emmm.201100209.PubMed CentralView ArticlePubMedGoogle Scholar
- Kantharidis P, Wang B, Carew RM, Lan HY: Diabetes complications: the microRNA perspective. Diabetes. 2011, 60 (7): 1832-1837. 10.2337/db11-0082.PubMed CentralView ArticlePubMedGoogle Scholar
- Krek A, Grun D, Poy MN, Wolf R, Rosenberg L, Epstein EJ, MacMenamin P, da Piedade I, Gunsalus KC, Stoffel M, et al: Combinatorial microRNA target predictions. Nat Genet. 2005, 37 (5): 495-500. 10.1038/ng1536.View ArticlePubMedGoogle Scholar
- Friedman RC, Farh KK, Burge CB, Bartel DP: Most mammalian mRNAs are conserved targets of microRNAs. Genome Res. 2009, 19 (1): 92-105.PubMed CentralView ArticlePubMedGoogle Scholar
- Wu S, Huang S, Ding J, Zhao Y, Liang L, Liu T, Zhan R, He X: Multiple microRNAs modulate p21Cip1/Waf1 expression by directly targeting its 3' untranslated region. Oncogene. 2010, 29 (15): 2302-2308. 10.1038/onc.2010.34.View ArticlePubMedGoogle Scholar
- Enright AJ, John B, Gaul U, Tuschl T, Sander C, Marks DS: MicroRNA targets in Drosophila. Genome Biol. 2003, 5 (1): R1-10.1186/gb-2003-5-1-r1.PubMed CentralView ArticlePubMedGoogle Scholar
- Lewis BP, Burge CB, Bartel DP: Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell. 2005, 120 (1): 15-20. 10.1016/j.cell.2004.12.035.View ArticlePubMedGoogle Scholar
- Lytle JR, Yario TA, Steitz JA: Target mRNAs are repressed as efficiently by microRNA-binding sites in the 5' UTR as in the 3' UTR. Proc Natl Acad Sci U S A. 2007, 104 (23): 9667-9672. 10.1073/pnas.0703820104.PubMed CentralView ArticlePubMedGoogle Scholar
- Nelson PT, Wang WX, Mao G, Wilfred BR, Xie K, Jennings MH, Gao Z, Wang X: Specific sequence determinants of miR-15/107 microRNA gene group targets. Nucleic Acids Res. 2011, 39 (18): 8163-8172. 10.1093/nar/gkr532.PubMed CentralView ArticlePubMedGoogle Scholar
- Alexiou P, Maragkakis M, Papadopoulos GL, Reczko M, Hatzigeorgiou AG: Lost in translation: an assessment and perspective for computational microRNA target identification. Bioinformatics. 2009, 25 (23): 3049-3055. 10.1093/bioinformatics/btp565.View ArticlePubMedGoogle Scholar
- Sethupathy P, Megraw M, Hatzigeorgiou AG: A guide through present computational approaches for the identification of mammalian microRNA targets. Nat Methods. 2006, 3 (11): 881-886. 10.1038/nmeth954.View ArticlePubMedGoogle Scholar
- Muniategui A, Pey J, Planes F, Rubio A: Joint analysis of miRNA and mRNA expression data. Brief Bioinform. 2012, 10.1093/bib/bbs028.Google Scholar
- Wurm FM: Production of recombinant protein therapeutics in cultivated mammalian cells. Nat Biotechnol. 2004, 22 (11): 1393-1398. 10.1038/nbt1026.View ArticlePubMedGoogle Scholar
- Xu X, Nagarajan H, Lewis NE, Pan S, Cai Z, Liu X, Chen W, Xie M, Wang W, Hammond S, et al: The genomic sequence of the Chinese hamster ovary (CHO)-K1 cell line. Nat Biotechnol. 2011, 29 (8): 735-741. 10.1038/nbt.1932.PubMed CentralView ArticlePubMedGoogle Scholar
- Birzele F, Schaub J, Rust W, Clemens C, Baum P, Kaufmann H, Weith A, Schulz TW, Hildebrandt T: Into the unknown: expression profiling without genome sequence information in CHO by next generation sequencing. Nucleic Acids Res. 2010, 38 (12): 3999-4010. 10.1093/nar/gkq116.PubMed CentralView ArticlePubMedGoogle Scholar
- Hammond S, Swanberg JC, Kaplarevic M, Lee KH: Genomic sequencing and analysis of a Chinese hamster ovary cell line using Illumina sequencing technology. BMC Genomics. 2011, 12: 67-10.1186/1471-2164-12-67.PubMed CentralView ArticlePubMedGoogle Scholar
- Gammell P, Barron N, Kumar N, Clynes M: Initial identification of low temperature and culture stage induction of miRNA expression in suspension CHO-K1 cells. J Biotechnol. 2007, 130 (3): 213-218. 10.1016/j.jbiotec.2007.04.020.View ArticlePubMedGoogle Scholar
- Johnson KC, Jacob NM, Nissom PM, Hackl M, Lee LH, Yap M, Hu WS: Conserved microRNAs in Chinese hamster ovary cell lines. Biotechnol Bioeng. 2010, 108 (2): 475-480.View ArticleGoogle Scholar
- Hackl M, Jakobi T, Blom J, Doppmeier D, Brinkrolf K, Szczepanowski R, Bernhart SH, Siederdissen CH, Bort JA, Wieser M, et al: Next-generation sequencing of the Chinese hamster ovary microRNA transcriptome: Identification, annotation and profiling of microRNAs as targets for cellular engineering. J Biotechnol. 2011, 153 (1–2): 62-75.PubMed CentralView ArticlePubMedGoogle Scholar
- Hammond S, Swanberg JC, Polson SW, Lee KH: Profiling conserved microRNA expression in recombinant CHO cell lines using Illumina sequencing. Biotechnol Bioeng. 2012, 109 (6): 1371-1375. 10.1002/bit.24415.PubMed CentralView ArticlePubMedGoogle Scholar
- Hernandez Bort JA, Hackl M, Hoflmayer H, Jadhav V, Harreither E, Kumar N, Ernst W, Grillari J, Borth N: Dynamic mRNA and miRNA profiling of CHO-K1 suspension cell cultures. Biotechnol J. 2012, 7 (4): 500-515. 10.1002/biot.201100143.View ArticlePubMedGoogle Scholar
- Barron N, Kumar N, Sanchez N, Doolan P, Clarke C, Meleady P, O'Sullivan F, Clynes M: Engineering CHO cell growth and recombinant protein productivity by overexpression of miR-7. J Biotechnol. 2011, 151 (2): 204-211. 10.1016/j.jbiotec.2010.12.005.View ArticlePubMedGoogle Scholar
- Druz A, Chu C, Majors B, Santuary R, Betenbaugh M, Shiloach J: A novel microRNA mmu-miR-466h affects apoptosis regulation in mammalian cells. Biotechnol Bioeng. 2011, 108 (7): 1651-1661. 10.1002/bit.23092.PubMed CentralView ArticlePubMedGoogle Scholar
- Clarke C, Doolan P, Barron N, Meleady P, O'Sullivan F, Gammell P, Melville M, Leonard M, Clynes M: Large scale microarray profiling and coexpression network analysis of CHO cells identifies transcriptional modules associated with growth and productivity. J Biotechnol. 2011, 155 (3): 350-359. 10.1016/j.jbiotec.2011.07.011.View ArticlePubMedGoogle Scholar
- Elkan-Miller T, Ulitsky I, Hertzano R, Rudnicki A, Dror AA, Lenz DR, Elkon R, Irmler M, Beckers J, Shamir R, et al: Integration of transcriptomics, proteomics, and microRNA analyses reveals novel microRNA regulation of targets in the mammalian inner ear. PLoS One. 2011, 6 (4): e18195-10.1371/journal.pone.0018195.PubMed CentralView ArticlePubMedGoogle Scholar
- Somel M, Guo S, Fu N, Yan Z, Hu HY, Xu Y, Yuan Y, Ning Z, Hu Y, Menzel C, et al: MicroRNA, mRNA, and protein expression link development and aging in human and macaque brain. Genome Res. 2010, 20 (9): 1207-1218. 10.1101/gr.106849.110.PubMed CentralView ArticlePubMedGoogle Scholar
- Almeida MI, Nicoloso MS, Zeng L, Ivan C, Spizzo R, Gafa R, Xiao L, Zhang X, Vannini I, Fanini F, et al: Strand-Specific miR-28-5p and miR-28-3p Have Distinct Effects in Colorectal Cancer Cells. Gastroenterology. 2012, 142 (4): 886-896. 10.1053/j.gastro.2011.12.047. e889PubMed CentralView ArticlePubMedGoogle Scholar
- Godlewski J, Nowicki MO, Bronisz A, Nuovo G, Palatini J, De Lay M, Van Brocklyn J, Ostrowski MC, Chiocca EA, Lawler SE: MicroRNA-451 regulates LKB1/AMPK signaling and allows adaptation to metabolic stress in glioma cells. Mol Cell. 2010, 37 (5): 620-632. 10.1016/j.molcel.2010.02.018.PubMed CentralView ArticlePubMedGoogle Scholar
- Laurila EM, Sandstrom S, Rantanen LM, Autio R, Kallioniemi A: Both inhibition and enhanced expression of miR-31 lead to reduced migration and invasion of pancreatic cancer cells. Genes Chromosomes Cancer. 2012, 51 (6): 557-568. 10.1002/gcc.21941.View ArticlePubMedGoogle Scholar
- Wang R, Wang ZX, Yang JS, Pan X, De W, Chen LB: MicroRNA-451 functions as a tumor suppressor in human non-small cell lung cancer by targeting ras-related protein 14 (RAB14). Oncogene. 2011, 30 (23): 2644-2658. 10.1038/onc.2010.642.View ArticlePubMedGoogle Scholar
- Yan D, Dong Xda E, Chen X, Wang L, Lu C, Wang J, Qu J, Tu L: MicroRNA-1/206 targets c-Met and inhibits rhabdomyosarcoma development. J Biol Chem. 2009, 284 (43): 29596-29604. 10.1074/jbc.M109.020511.PubMed CentralView ArticlePubMedGoogle Scholar
- Yan HL, Xue G, Mei Q, Wang YZ, Ding FX, Liu MF, Lu MH, Tang Y, Yu HY, Sun SH: Repression of the miR-17-92 cluster by p53 has an important function in hypoxia-induced apoptosis. EMBO J. 2009, 28 (18): 2719-2732. 10.1038/emboj.2009.214.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang T, Liu M, Wang C, Lin C, Sun Y, Jin D: Down-regulation of MiR-206 promotes proliferation and invasion of laryngeal cancer by regulating VEGF expression. Anticancer Res. 2011, 31 (11): 3859-3863.PubMedGoogle Scholar
- Olive V, Jiang I, He L: mir-17-92, a cluster of miRNAs in the midst of the cancer network. Int J Biochem Cell Biol. 2010, 42 (8): 1348-1354. 10.1016/j.biocel.2010.03.004.PubMed CentralView ArticlePubMedGoogle Scholar
- O'Donnell KA, Wentzel EA, Zeller KI, Dang CV, Mendell JT: c-Myc-regulated microRNAs modulate E2F1 expression. Nature. 2005, 435 (7043): 839-843. 10.1038/nature03677.View ArticlePubMedGoogle Scholar
- Schulte JH, Horn S, Otto T, Samans B, Heukamp LC, Eilers UC, Krause M, Astrahantseff K, Klein-Hitpass L, Buettner R, et al: MYCN regulates oncogenic MicroRNAs in neuroblastoma. Int J Cancer. 2008, 122 (3): 699-704. 10.1002/ijc.23153.View ArticlePubMedGoogle Scholar
- Puck TT, Cieciura SJ, Robinson A: Genetics of somatic mammalian cells. III. Long-term cultivation of euploid cells from human and animal subjects. J Exp Med. 1958, 108 (6): 945-956.PubMedGoogle Scholar
- Lam EK, Wang X, Shin VY, Zhang S, Morrison H, Sun J, Ng EK, Yu J, Jin H: A microRNA contribution to aberrant Ras activation in gastric cancer. Am J Transl Res. 2011, 3 (2): 209-218.PubMed CentralPubMedGoogle Scholar
- Li D, Zhao Y, Liu C, Chen X, Qi Y, Jiang Y, Zou C, Zhang X, Liu S, Wang X, et al: Analysis of MiR-195 and MiR-497 expression, regulation and role in breast cancer. Clin Cancer Res. 2011, 17 (7): 1722-1730. 10.1158/1078-0432.CCR-10-1800.View ArticlePubMedGoogle Scholar
- Wu F, Zhu S, Ding Y, Beck WT, Mo YY: MicroRNA-mediated regulation of Ubc9 expression in cancer cells. Clin Cancer Res. 2009, 15 (5): 1550-1557. 10.1158/1078-0432.CCR-08-0820.PubMed CentralView ArticlePubMedGoogle Scholar
- Bian HB, Pan X, Yang JS, Wang ZX, De W: Upregulation of microRNA-451 increases cisplatin sensitivity of non-small cell lung cancer cell line (A549). J Exp Clin Cancer Res. 2011, 30: 20-10.1186/1756-9966-30-20.PubMed CentralView ArticlePubMedGoogle Scholar
- Gandin V, Miluzio A, Barbieri AM, Beugnet A, Kiyokawa H, Marchisio PC, Biffo S: Eukaryotic initiation factor 6 is rate-limiting in translation, growth and transformation. Nature. 2008, 455 (7213): 684-688. 10.1038/nature07267.PubMed CentralView ArticlePubMedGoogle Scholar
- Fingar DC, Richardson CJ, Tee AR, Cheatham L, Tsou C, Blenis J: mTOR controls cell cycle progression through its cell growth effectors S6K1 and 4E-BP1/eukaryotic translation initiation factor 4E. Mol Cell Biol. 2004, 24 (1): 200-216. 10.1128/MCB.24.1.200-216.2004.PubMed CentralView ArticlePubMedGoogle Scholar
- Warner JR, McIntosh KB: How common are extraribosomal functions of ribosomal proteins?. Mol Cell. 2009, 34 (1): 3-11. 10.1016/j.molcel.2009.03.006.PubMed CentralView ArticlePubMedGoogle Scholar
- Moss T, Langlois F, Gagnon-Kugler T, Stefanovsky V: A housekeeper with power of attorney: the rRNA genes in ribosome biogenesis. Cell Mol Life Sci. 2007, 64 (1): 29-49. 10.1007/s00018-006-6278-1.View ArticlePubMedGoogle Scholar
- Arabi A, Wu S, Ridderstrale K, Bierhoff H, Shiue C, Fatyol K, Fahlen S, Hydbring P, Soderberg O, Grummt I, et al: c-Myc associates with ribosomal DNA and activates RNA polymerase I transcription. Nat Cell Biol. 2005, 7 (3): 303-310. 10.1038/ncb1225.View ArticlePubMedGoogle Scholar
- Budde A, Grummt I: p53 represses ribosomal gene transcription. Oncogene. 1999, 18 (4): 1119-1124. 10.1038/sj.onc.1202402.View ArticlePubMedGoogle Scholar
- Henras AK, Soudet J, Gerus M, Lebaron S, Caizergues-Ferrer M, Mougin A, Henry Y: The post-transcriptional steps of eukaryotic ribosome biogenesis. Cell Mol Life Sci. 2008, 65 (15): 2334-2359. 10.1007/s00018-008-8027-0.View ArticlePubMedGoogle Scholar
- Wang ET, Sandberg R, Luo S, Khrebtukova I, Zhang L, Mayr C, Kingsmore SF, Schroth GP, Burge CB: Alternative isoform regulation in human tissue transcriptomes. Nature. 2008, 456 (7221): 470-476. 10.1038/nature07509.PubMed CentralView ArticlePubMedGoogle Scholar
- Bourdon JC: p53 and its isoforms in cancer. Br J Cancer. 2007, 97 (3): 277-282. 10.1038/sj.bjc.6603886.PubMed CentralView ArticlePubMedGoogle Scholar
- Davies AH, Barrett I, Pambid MR, Hu K, Stratford AL, Freeman S, Berquin IM, Pelech S, Hieter P, Maxwell C, et al: YB-1 evokes susceptibility to cancer through cytokinesis failure, mitotic dysfunction and HER2 amplification. Oncogene. 2011, 30 (34): 3649-3660. 10.1038/onc.2011.82.PubMed CentralView ArticlePubMedGoogle Scholar
- Raffetseder U, Frye B, Rauen T, Jurchott K, Royer HD, Jansen PL, Mertens PR: Splicing factor SRp30c interaction with Y-box protein-1 confers nuclear YB-1 shuttling and alternative splice site selection. J Biol Chem. 2003, 278 (20): 18241-18248. 10.1074/jbc.M212518200.View ArticlePubMedGoogle Scholar
- Wu J, Lee C, Yokom D, Jiang H, Cheang MC, Yorida E, Turbin D, Berquin IM, Mertens PR, Iftner T, et al: Disruption of the Y-box binding protein-1 results in suppression of the epidermal growth factor receptor and HER-2. Cancer Res. 2006, 66 (9): 4872-4879. 10.1158/0008-5472.CAN-05-3561.View ArticlePubMedGoogle Scholar
- Keene JD: RNA regulons: coordination of post-transcriptional events. Nat Rev Genet. 2007, 8 (7): 533-543. 10.1038/nrg2111.View ArticlePubMedGoogle Scholar
- Chen X, Yan Q, Li S, Zhou L, Yang H, Yang Y, Liu X, Wan X: Expression of the tumor suppressor miR-206 is associated with cellular proliferative inhibition and impairs invasion in ERalpha-positive endometrioid adenocarcinoma. Cancer Lett. 2012, 314 (1): 41-53. 10.1016/j.canlet.2011.09.014.View ArticlePubMedGoogle Scholar
- Li C, Nie H, Wang M, Su L, Li J, Yu B, Wei M, Ju J, Yu Y, Yan M, et al: MicroRNA-409-3p regulates cell proliferation and apoptosis by targeting PHF10 in gastric cancer. Cancer Lett. 2012, 320 (2): 189-197. 10.1016/j.canlet.2012.02.030.View ArticlePubMedGoogle Scholar
- Notari M, Neviani P, Santhanam R, Blaser BW, Chang JS, Galietta A, Willis AE, Roy DC, Caligiuri MA, Marcucci G, et al: A MAPK/HNRPK pathway controls BCR/ABL oncogenic potential by regulating MYC mRNA translation. Blood. 2006, 107 (6): 2507-2516. 10.1182/blood-2005-09-3732.PubMed CentralView ArticlePubMedGoogle Scholar
- Guil S, Caceres JF: The multifunctional RNA-binding protein hnRNP A1 is required for processing of miR-18a. Nat Struct Mol Biol. 2007, 14 (7): 591-596. 10.1038/nsmb1250.View ArticlePubMedGoogle Scholar
- Michlewski G, Caceres JF: Antagonistic role of hnRNP A1 and KSRP in the regulation of let-7a biogenesis. Nat Struct Mol Biol. 2010, 17 (8): 1011-1018. 10.1038/nsmb.1874.PubMed CentralView ArticlePubMedGoogle Scholar
- Eiring AM, Harb JG, Neviani P, Garton C, Oaks JJ, Spizzo R, Liu S, Schwind S, Santhanam R, Hickey CJ, et al: miR-328 functions as an RNA decoy to modulate hnRNP E2 regulation of mRNA translation in leukemic blasts. Cell. 2010, 140 (5): 652-665. 10.1016/j.cell.2010.01.007.PubMed CentralView ArticlePubMedGoogle Scholar
- Subramani D, Alahari SK: Integrin-mediated function of Rab GTPases in cancer progression. Mol Cancer. 2010, 9: 312-10.1186/1476-4598-9-312.PubMed CentralView ArticlePubMedGoogle Scholar
- Fasanaro P, Greco S, Lorenzi M, Pescatori M, Brioschi M, Kulshreshtha R, Banfi C, Stubbs A, Calin GA, Ivan M, et al: An integrated approach for experimental target identification of hypoxia-induced miR-210. J Biol Chem. 2009, 284 (50): 35134-35143. 10.1074/jbc.M109.052779.PubMed CentralView ArticlePubMedGoogle Scholar
- Dvinge H, Bertone P: HTqPCR: high-throughput analysis and visualization of quantitative real-time PCR data in R. Bioinformatics. 2009, 25 (24): 3325-3326. 10.1093/bioinformatics/btp578.PubMed CentralView ArticlePubMedGoogle Scholar
- Doolan P, Melville M, Gammell P, Sinacore M, Meleady P, McCarthy K, Francullo L, Leonard M, Charlebois T, Clynes M: Transcriptional profiling of gene expression changes in a PACE-transfected CHO DUKX cell line secreting high levels of rhBMP-2. Mol Biotechnol. 2008, 39 (3): 187-199. 10.1007/s12033-008-9039-6.View ArticlePubMedGoogle Scholar
- Meleady P, Gallagher M, Clarke C, Henry M, Sanchez N, Barron N, Clynes M: Impact of miR-7 over-expression on the proteome of Chinese hamster ovary cells. J Biotechnol. 2012, 160 (3–4): 251-262.View ArticlePubMedGoogle Scholar
- Meleady P, Hoffrogge R, Henry M, Rupp O, Bort JH, Clarke C, Brinkrolf K, Kelly S, Muller B, Doolan P, et al: Utilization and evaluation of CHO-specific sequence databases for mass spectrometry based proteomics. Biotechnol Bioeng. 2012, 109 (6): 1386-1394. 10.1002/bit.24476.View ArticlePubMedGoogle Scholar
- Fagan A, Culhane AC, Higgins DG: A multivariate analysis approach to the integration of proteomic and gene expression data. Proteomics. 2007, 7 (13): 2162-2171. 10.1002/pmic.200600898.View ArticlePubMedGoogle Scholar
- Madden SF, Carpenter SB, Jeffery IB, Bjorkbacka H, Fitzgerald KA, O'Neill LA, Higgins DG: Detecting microRNA activity from gene expression data. BMC Bioinformatics. 2010, 11: 257-10.1186/1471-2105-11-257.PubMed CentralView ArticlePubMedGoogle Scholar
- Culhane AC, Thioulouse J, Perriere G, Higgins DG: MADE4: an R package for multivariate analysis of gene expression data. Bioinformatics. 2005, 21 (11): 2789-2790. 10.1093/bioinformatics/bti394.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.