- Research article
- Open Access
The deregulated microRNAome contributes to the cellular response to aneuploidy
BMC Genomics volume 19, Article number: 197 (2018)
Aneuploidy, or abnormal chromosome numbers, severely alters cell physiology and is widespread in cancers and other pathologies. Using model cell lines engineered to carry one or more extra chromosomes, it has been demonstrated that aneuploidy per se impairs proliferation, leads to proteotoxic as well as replication stress and triggers conserved transcriptome and proteome changes.
In this study, we analysed for the first time miRNAs and demonstrate that their expression is altered in response to chromosome gain. The miRNA deregulation is independent of the identity of the extra chromosome and specific to individual cell lines. By cross-omics analysis we demonstrate that although the deregulated miRNAs differ among individual aneuploid cell lines, their known targets are predominantly associated with cell development, growth and proliferation, pathways known to be inhibited in response to chromosome gain. Indeed, we show that up to 72% of these targets are downregulated and the associated miRNAs are overexpressed in aneuploid cells, suggesting that the miRNA changes contribute to the global transcription changes triggered by aneuploidy. We identified hsa-miR-10a-5p to be overexpressed in majority of aneuploid cells. Hsa-miR-10a-5p enhances translation of a subset of mRNAs that contain so called 5’TOP motif and we show that its upregulation in aneuploids provides resistance to starvation-induced shut down of ribosomal protein translation.
Our work suggests that the changes of the microRNAome contribute on one hand to the adverse effects of aneuploidy on cell physiology, and on the other hand to the adaptation to aneuploidy by supporting translation under adverse conditions.
A balanced karyotype is essential for cell viability and therefore aneuploidy, characterized by unbalanced changes in chromosome numbers and sub-chromosomal structural variations, has often profound detrimental consequences for cell physiology. In humans, aneuploidy is the major cause of spontaneous abortions and the few trisomies compatible with survival result in severe developmental defects . Aneuploidy in somatic cells is frequently associated with cancer, as 70% of haematopoietic and 90% of solid cancers show an abnormal karyotype [2, 3]. Recently developed aneuploid model systems in several different species have accelerated research on the effects of aneuploidy per se. The physiological changes in response to an unbalanced karyotype are multifold, including impaired proliferation, replication stress and proteotoxic stress that is characterized by changes in protein stoichiometry, reduced protein folding capacity and by activation of autophagy [4,5,6,7,8,9,10,11,12,13]. The physiological response to aneuploidy goes hand in hand with a transcriptional response that is manifested in a conserved pathway deregulation [14, 15]. What triggers these changes in gene expression and how exactly this specific response is modulated has not been clarified so far.
Given that the triggers of the transcriptional deregulations remain unclear, we asked whether microRNA (miRNA) regulation is involved in the response to aneuploidy. MiRNAs are small non-coding RNA molecules that posttranscriptionally modulate gene expression of over 60% of protein coding genes . The regulatory function of miRNAs is mediated by their binding to the target 3’untranslated region (UTR) via partially complementary sequences, thereby inducing translational repression and/or mRNA decay. Alternatively, some miRNAs bind to sites within the coding region or the 5’UTR [17, 18]. This binding can potentially stabilize the mRNA or enhance its translation . One mRNA can be affected by multiple miRNAs that cooperatively translationally repress or degrade the target mRNA or compete for target regulations . In turn, single miRNA can repress hundreds of mRNAs and lead to large-scale transcriptome changes that for example play a key role in stem cell differentiation . Moreover, the complex miRNA-target network allows to fine-tune diverse cellular processes by modulating the amount of transcripts translated particularly upon stress conditions or changes in the environment [22, 23].
In the disease context such as in cancer, alterations of miRNA expression are common and specific miRNA deregulation profiles are sufficient to distinguish cancerous from non-cancerous tissues and to predict invasiveness and aggressiveness of various cancers [24,25,26]. The widespread miRNA deregulation in cancers has been attributed to genomic copy number changes. For instance, a gain of chromosome 1q relates to miRNA expression changes in cervical cancer  and downregulation of let-7 family was associated with copy number changes in medulloblastoma, breast and ovarian cancer . In some cancers, deregulated miRNAs promote genomic and chromosomal instability by targeting the mitotic checkpoint or DNA damage repair components [29, 30]. Yet, beyond these few examples, we lack a deeper understanding of how aneuploidy in cancer and miRNAs expression changes are linked. Moreover, the relation of miRNA and aneuploidy per se has not been studied so far.
We asked whether aneuploidy affects miRNA abundance and analysed the phenotypic consequences of miRNA changes in aneuploid cells (Fig. 1a). To this end, we used a series of human trisomic and tetrasomic cell lines that very previously constructed in our laboratory by micronuclei-mediated chromosome transfer into diploid, chromosomally stable human cell lines HCT116 and RPE1 (, for further details see Methods). By transferring the individual chromosome, we obtained homogeneous population with a defined aberrant karyotype. Our model cell lines therefore allow study the general consequences of chromosome gain by direct comparison of different trisomies and tetrasomies to their isogenic parental cells. Using this model system, we found that chromosome gains strongly deregulate the expression of more than 25% of miRNAs, although only a few individual miRNAs were commonly altered among different aneuploidies. Most of the identified miRNA deregulations negatively affect cell growth and proliferation, thus suggesting that miRNAs suppress proliferation of aneuploid cells. We identified miRNA hsa-miR-10a-5p to be overexpressed in majority of analysed cell lines. This miRNA acts through both the canonical 3’UTR targeting as well as via the 5’UTR of target mRNAs. We show that the alternative 5’UTR targeting provides the aneuploid cells with resistance to starvation by ensuring sustained translation of ribosomal mRNAs upon serum withdrawal. Our analysis of the global miRNAome changes in response to chromosome gain reveals a complex interplay of the gene expression regulatory mechanisms that shape the global transcriptome and proteome dynamics in aneuploid cells.
Deregulation of miRNAome in human aneuploid model cell lines
To determine the effects of chromosome gain on miRNA expression in human cells, we used a series of cells derived from HCT116 and RPE1 cell lines that contain one or more extra copies of different chromosomes ([5, 9, 12], Fig. 1a, see Methods for more details). Following cell lines were subjected to small RNA sequencing: four HCT116- derived cell lines that gained chromosome 3, 8, 5 or 18, respectively, the parental HCT116 cell line, three RPE1- derived cell lines trisomic for chromosome 7, 21, cell line trisomic for both chromosome 5 and 12 and the parental RPE1 cell line. The miRBase repository , where 1881 mature human miRNAs are listed to date, was used as the miRNA data source for mapping. Mapping of the raw sequences to the human genome and subsequent identification of miRNAs with mirdeep2  resulted in at least 554 and up to 719 identified mature miRNAs in the sequenced cell lines (Table 1).
While the numbers of identified miRNAs were similar for all sequenced aneuploid cell lines, the extent of miRNA deregulation differed. We analyzed the differential miRNA expression in aneuploid cell lines in comparison to their diploid counterparts by applying DESeq2 for normalization and statistical testing . Since miRNAs with very low expression levels are characterized by low read counts and an inherent large variability on the logarithmic scale, these miRNAs were pre-filtered . This filtering for miRNAs with a read count greater than 10 reduced the number of identified miRNAs to 231–293 (Table 2). Using this dataset, we then determined the significantly deregulated miRNAs as miRNAs with the aneuploidy/diploid log2 fold change +/− 0.6 at adjusted p-value < 0.05. We identified from 23 (in RPE1* 21/3) to 74 (in RPE1 5/3 12/3 and HCT116 5/4) significantly deregulated miRNAs (Table 2). The percentage of deregulated miRNAs ranged from 10% (in RPE1* 21/3) up to 31.2% (in RPE1 5/3 12/3). There was only a weak correlation between the percentages of deregulated miRNAs with the amount of extra DNA (Additional file 1: Figure S1). Thus, gain of even a single chromosome deregulates up to one third of cellular miRNAs.
Cell line specific and general miRNA deregulation in response to aneuploidy
Comparison of the significantly altered miRNAs in aneuploid cell lines revealed a pattern of miRNA expression changes largely distinct for each specific cell line (Fig. 1b). The Euclidean distance clustering was not defined by the parental cell type, but rather by distinct and cell line specific miRNA expression clusters individual for each aneuploid cell line. HCT116 5/4 showed the most prominent clusters of down- and upregulated miRNAs and appeared as an outlier in the Euclidean distance cluster. The upregulated miRNA cluster in HCT116 5/4 partially overlapped with the upregulated miRNAs in RPE1 5/3 12/3; among them most prominently the miRNA family miR-192/215, the nearby located miR-194 and miRNAs hsa-miR-29b and hsa-miR-29c, which are encoded in a cluster on chromosome 1. Only a few similarities were found among different cell lines. For instance, hsa-miR-21-5p and hsa-miR-24-2-5p showed similar upregulation in all HCT116-derived cell lines (Fig. 1c). Moreover, five out of seven sequenced aneuploid cell lines shared the upregulation of hsa-miR-139-5p and hsa-miR-10a-5p (Fig. 1d). The latter upregulation did not affect the entire miRNA family and seems to be specific for the individual miRNA: while hsa-miR-10a-5p was upregulated in five out of seven of the analyzed cell lines, the family member hsa-miR-10b-5p was upregulated only in RPE1 5/3 12/3. Only a few globally downregulated miRNAs were identified, such as hsa-miR-374b-5p that was downregulated in all cell lines except RPE1 5/3 12/3 (Fig. 1e).
The presence of extra chromosomes in cells generally leads to an elevated expression of the genes that are located on the supernumerary chromosomes; thus the abundance of transcripts scales with the gene copy numbers [6, 12,13,14,15, 34, 35]. To determine whether this holds also true for miRNA expression, we ordered the miRNAs according to their chromosome location. We then compared the distribution of the miRNA expression from the trisomic/tetrasomic chromosomes to miRNA expression from the disomes and tested whether the elevated expression occurred just by chance. The chromosome-specific miRNA expression was significantly higher than the disomic miRNA expression only for RPE1* 7/3, RPE1* 21/3 and for the miRNAs encoded on chromosome 5 in RPE1 5/3 12/3 (Fig. 2a–c). In all other aneuploid cell lines, the miRNA expression distribution did not increase significantly with the number of chromosome copies (Fig. 2c–g, Additional file 2: Table S1). In addition, the median miRNA expression differed from the expected levels even for some disomic chromosomes in several cell lines (Additional file 1: Figure S2A). For instance, the median expression of miRNAs encoded on chromosome 11 and 17 were significantly increased in HCT116 5/4 despite their disomic status that was confirmed by DNA sequencing. Taken together, gain of even a single chromosome leads to strong deregulation of up to 30% of miRNAs. The largely increased abundance of the deregulated miRNAs cannot be explained by the chromosome copy number changes.
Effects of the deregulated miRNAome on transcriptome
We asked to what extent the mRNA expression changes in response to aneuploidy can be explained by miRNA deregulation. The transcriptome of HCT116 5/4, HCT116* 8/3 c7, RPE1 5/3 12/3 and RPE1* 21/3 was sequenced. The percentage of deregulated mRNAs (log2 fold change +/− 0.6, adjusted p-value < 0.05) highly differed among the cell lines, with 3.3% in RPE1* 21/3, 16.6% in RPE1 5/3 12/3, 3.0% in HCT116* 8/3 c7 and 20.8% in HCT116 5/4 (Table 2). To model the miRNA-target network, we retrieved reported miRNA targets from miRTarBase (v6.1) for all deregulated miRNAs of each cell line . We then filtered these miRNA targets for those that were altered in response to aneuploidy. Although in principle each target can be affected by several miRNAs, we simplified the complex miRNA-mRNA relation and filtered for pairs of individual deregulated mRNA targets that were inversely expressed relative to the individual corresponding miRNAs. This analysis revealed that 38% of mRNAs that are deregulated in HCT116 5/4 are potentially targeted by the deregulated miRNAs, as judged by their inverse expression, 24% in RPE1* 21/3, 27% in RPE1 5/3 12/3 and only 4.8% HCT116* 8/3 c7 (Table 3). Thus, under these approximations, only a minor fraction of the mRNA changes in response to aneuploidy might be directly influenced by miRNA regulation.
Deregulated miRNAs negatively affect cell cycle associated pathways
To determine which cellular functions are affected by the changes in the miRNAome in aneuploid cells, we performed functional analysis based on the manually curated knowledge base of Ingenuity Pathway Analysis , where the miRNA annotations are derived from their published functional roles, independently of their validated and predicted targets. Pathway “Cell Cycle” was affected in all aneuploid cell lines and additional four molecular and cellular functions affected by the deregulated miRNAs were shared in at least five out of seven aneuploid cell lines (Fig. 3a). The miRNAome of HCT116* 3/3 showed less functional associations and a low number of involved miRNAs, although pathways “Cellular Development” and “Cell Cycle” were also affected (Fig. 3a).
We employed IPA to infer an activation state of the biological functions as an activation z-score. The IPA activation z-score is a quantitative measure of the predicted effect of the deregulated miRNAome on a biological function . It is calculated based on the direction of miRNA deregulation and the literature-derived positive or negative effects of these miRNAs on a biological function. The available miRNA dataset allowed to retrieve the z-score for the sub-functional annotation terms defined by IPA within the categories “Cell Death and Survival”, “Cellular Development, Cellular Growth and Proliferation”, “Cellular Growth and Proliferation” and “Cellular Movement” (Fig. 3b). The effect of the deregulated miRNAome was mostly negative, as the sub-functional annotation terms within the four parent categories were predicted to be inhibited or deactivated in HCT116 18/3, HCT116 5/4 and RPE1* 7/3. The activation z-score of HCT116* 3/3 miRNA deregulation predicted the changes for the category “Cell Death and Survival”, which is in concordance with the identity of the five top molecular and cellular functions (Fig. 3). Therein, “apoptosis” was predicted to be activated with a z-score of 1.8 and 2 and “cell death” as well as “necrosis” with a z-score of 1.5 and 1.4, respectively. Thus, the negative impact of miRNA deregulation might be executed in HCT116* 3/3 rather via upregulation of “Cell Death” pathways than by downregulation of “Cellular Growth and Proliferation” pathways.
To examine the effect of identified miRNAs on cell cycle, we transfected mimics and inhibitors of the only two miRNA candidates that were upregulated in at least five out of seven analysed aneuploid cell lines: miR-10a-5p and miR-139-5p (Fig. 1b-d). Whereas transfections of the hsa-miRNA-10a-5p mimics and inhibitors had only a negligible effect on proliferation and expression of cell cycle markers (data not shown), we found that inhibition of miRNA 139-5p indeed decreased the levels of negative cell cycle regulators CDKN1A (p21) and CDKN1B (p27) (Fig. 3c). MiRNA-139-5p is located on chromosome 11 at 11q13.4 and its expression is downregulated in a variety of cancers including hepatocellular carcinoma, gastric cancer, breast cancer and glioblastoma (e.g ). In summary, although the deregulated miRNAs differed in individual cells lines, we found a striking overlap in their predicted, largely negative impact on cellular growth and proliferation.
Integration of miRNA, RNA and protein expression data confirms the negative effect on cellular development, cellular growth and proliferation
To illuminate how the miRNAs execute their negative impact on cellular development, growth and proliferation in aneuploids, we analyzed their targets. For this purpose we exploited the fact that complete miRNAome, transcriptome and proteome data were available for three cell lines (HCT116 5/4, RPE1 5/3 12/3 and RPE1* 21/3, this work and ). In total, we identified 40 aneuploidy-responsive miRNAs to be associated with the cellular development, growth and proliferation category (Additional file 1: Figure S3A). Comparison with experimentally validated miRNA-target interactions revealed 604 interactions with 325 unique targets for HCT116 5/4 (Fig. 4a), 238 in RPE1 5/3 12/3 and 208 in RPE1* 21/3 (Additional file 1: Figure S3B, C). Chromosome alignment of the miRNAs revealed no bias for the superfluous chromosome in any of the cell lines.
To infer the effect of miRNA regulation on experimentally validated targets, we matched the transcriptome and proteome data to the targets of miRNAs within the cellular development, growth and proliferation category. Targets were filtered for their expression changes above or below a threshold ratio of +/− 0.6 (log2 fold change). Only experimentally validated targets with a potential miRNA relation (that means with inverse expression to the interacting miRNA) were considered for the annotation enrichment analysis. Strikingly, the protein abundance of 72% of the filtered targets was downregulated in HCT116 5/4, many of them being known cell cycle and proliferation- associated proteins (Fig. 4b). This observation holds also true for RPE1 5/3 12/3 and RPE1* 21/3 (Additional file 1: Figure S3D, E). Thus, the target gene expression analysis substantiates the predicted activation state of the molecular function and suggests a global negative effect of the miRNA changes on cell cycle and proliferation in response to aneuploidy.
Additionally, we performed functional annotation clustering of the miRNA targets with the DAVID functional annotation tool. This analysis confirmed the cell cycle functions as the most affected cluster (Additional file 3: Table S2). Among the deregulated targets in HCT116 5/4 we found key players of proliferation, such as the mitotic checkpoint serine/threonine kinases (BUB1), cell division cycling 20 (CDC20), cyclin-dependent kinase 4 (CDK4), checkpoint kinase 1 (CHEK1), high motility group A2 (HMGA2) and regulator of chromosome condensation (RCC2) (Fig. 4b). Cell cycle associated factors were also targeted in RPE1 5/3 12/3 and RPE1* 21/3, such as the minichromosome maintenance protein complex members (MCM2, 3 and 6) and Cyclin E1 (Additional file 1: Figure S3D, E). Immunoblotting of some of the targets confirmed their downregulation in response to aneuploidy (Fig. 4c). Taken together, deregulated miRNAome in aneuploid cells negatively affects expression of cell cycle factors and thereby may contribute to the proliferation defect observed in aneuploid model cell lines.
Hsa-miR-10a-5p activity is increased in human aneuploid model cell lines
While the majority of deregulated miRNAs were not shared among aneuploid cell lines, the miRNA hsa-miR-10a-5p was significantly upregulated in five out of seven aneuploid cell lines. Quantitative real-time PCR of 18 different aneuploid model cell lines (for details see Material and Methods) revealed that the mean expression of hsa-miR-10a-5p was significantly increased above the levels of the parental cell line in 7 out of 18 cell lines and elevated in additional 8 cell lines (Fig. 5a). Interestingly, the trisomies of chromosome 13, 18 and 21 that are compatible with survival in a whole organismal context (Patau, Edwards and Down syndrome, respectively) showed no or only negligible increase of hsa-miR-10a-5p expression. The increased expression of hsa-miR-10a-5p might be a secondary effect of the expression of HOXB3, in whose locus the miRNA sequence is located. Comparison of the hsa-miR-10a-5p and HOXB3 expression levels in four cell lines revealed that both miR-10a-5p and HOXB3 were significantly upregulated only in two cell lines, HCT116 5/4 and RPE1* 21/3, while their levels were independent in the other cell lines (Additional file 1: Figure S4). Thus, the upregulation of hsa-miR-10a-5p cannot be generally explained by increased expression of the host gene HOXB3.
To determine whether hsa-miR-10a-5p overexpression affects gene expression in aneuploid cells, we used a luciferase reporter assay to evaluate its translational repression efficiency (Ørom et al., 2008). The luciferase reporter consists of a miR-10a binding site in a synthetic 3’untranslated region (UTR) of the Renilla luciferase gene (Fig. 5b). The reporter additionally contains a Firefly luciferase gene that allows normalization for transfection efficiency. Indeed, the normalized luciferase signal was significantly lower in four out of five tested aneuploid cell lines (Fig. 5c). No changes in luciferase signal were detected in the HCT116* 8/3 c7 cell line; this is likely due to a relatively minor increase in the hsa-miR-10a-5p expression in this cell line (see Fig. 5a).
The target database lists 298 unique target genes for hsa-miR-10a-5p with only nine previously validated targets. Only six of the targets were downregulated on both transcriptome and proteome level in HCT116 5/4 (Additional file 1: Figure S5A). Only two of the 270 hsa-miR-10a targets were downregulated on transcriptome and proteome in both RPE1 5/3 12/3 and RPE1* 21/3 (Additional file 1: Figure S5B,C). This is in contrast with the other identified deregulated miRNAs, where majority of the targets showed inversed expression levels and suggests that the hsa-miR-10a-5p function in aneuploids might be independent of the 3’UTR target binding.
Hsa-miR-10a-5p overexpression governs resistance to translation stress
Next, we considered the fact that hsa-miR-10a-5p binds also to the 5’UTR downstream of the 5’TOP motif that is found in mRNAs of ribosomal proteins (RP) and other mRNAs of translation associated proteins . This binding results in enhanced translation of RP mRNAs and alleviates redistribution of RP mRNAs from active polysomes to inactive RNP complexes upon amino acid starvation. To analyse the translation efficiency of mRNAs with the 5’TOP motif, we employed a luciferase reporter containing the transcriptional start site of the ribosomal protein S16 (Rps16) and 29 nt of the exon 1 including the 5’TOP motif in the 5’UTR of a luciferase gene (pS16-wt-luc) (Fig. 6a, ). We first determined the steady state levels of 5’TOP motif mRNA translation with the pS16-wt-luc luciferase reporter and determined normalized Firefly luciferase signal similar to the wild type (Fig. 6b, mock control). Next, we measured the luciferase activity after overexpression of the hsa-miR-10a-5p mimic (Additional file 1: Figure S6A, B). Using the pS16-wt-luc luciferase reporter, we observed that hsa-miR-10a-5p overexpression enhanced the 5’UTR-mediated effect on the translation of 5’TOP motif mRNAs in the parental HCT116 cell line as well as in three out of four tested aneuploids (Fig. 6b). Thus, the expression of mRNAs with the 5’TOP motifs is sensitive to the changes in the hsa-miR-10a-5p abundance.
The translation of mRNAs with the 5’TOP motif is tightly regulated and cellular stresses such as starvation result in their translational repression and redistribution into inactive ribonucleoprotein complexes . We hypothesized that upregulation of hsa-miR-10a-5p in aneuploid cell lines might protect the cells from the translational repression. To test this hypothesis, we measured the 5’TOP motif mRNA translation via the luciferase reporter assay in serum-deprived cells. Serum starvation led to a significant decrease of the luciferase activity in the parental HCT116 cell line. In contrast, all aneuploid cell lines were resistant to the translational reduction of the 5’TOP motif mRNAs (Fig. 6c). Aneuploidy is known to inhibit pathways related to ribosome and translation [12,13,14,15]. We found that even a short treatment with cycloheximide, a eukaryotic protein synthesis inhibitor, induced the expression of hsa-miR-10a-5p, whereas other stress conditions did not significantly affect its expression (Fig. 6d, e). Thus, we propose that the negative effect of aneuploidy on translation leads to induction of hsa-miR-10a-5p that may in turn facilitate survival of aneuploids by its positive effect on translation.
Using an integrated approach of large-scale miRNA, RNA and protein expression data analysis, we document deregulation of the miRNAome and its consequences in human aneuploid model cell lines (Fig. 1a). For the first time, we show that addition of one or two chromosomes results in extensive genome-wide miRNA expression changes in human cells (Fig. 1b-e). The data analysis reveals that the deregulated miRNAome inflicts negative effects on development, growth and proliferation, and may thereby contribute to the cellular response to aneuploidy. Moreover, our evidence suggests that upregulation of hsa-miR-10a-5p in aneuploid cell lines might present an adaptation to aneuploidy by providing means to selectively increase translation.
Karyotype copy number alterations were previously associated with deregulated miRNA expression in some cancers and miRNA coding regions show high frequency of copy number variations in ovarian, breast, melanoma and lung cancer [41, 42]. Interestingly, some studies report no correlation between the copy number status and miRNA expression in haematopoetic cancer and acute myeloid leukemia [43, 44]. We observed that the miRNA abundance does not readily increase according to the chromosomes copy number (Fig. 2, Additional file 1: Figure S2), which is in a contrast with mRNA expression changes that largely scale up with the chromosome copy number changes [12,13,14,15]. The lack of general correlation between miRNAs abundance and chromosome copy number changes suggests that the copy number changes are counteracted by a tight regulation of miRNA expression.
Aneuploidy remodels the cellular miRNAome, as 9 to 30% of miRNAs were differentially expressed in aneuploids compared to isogenic controls (Tables 1 and 2). Although the deregulated miRNAs are mostly unique to the individual aneuploid cell lines, the affected cellular pathways appear the same: analysis of the cellular functions affected by the deregulated miRNAs identified cellular development, growth and proliferation as common targets (Fig. 3a, b). The integrated data analysis showed that the majority of targets, which show an inverse expression to the deregulated miRNAs, are indeed downregulated on mRNA or protein level (Fig. 4, Additional file 1: Figure S3). Among the downregulated targets are proteins that are crucial for the cell cycle progression such as BUB1, CDC20, CDK4 and CHEK1 in HCT116 5/4 or Cyclin E1 and RB1 in RPE1-derived aneuploid cell lines. One of the targets, HMGA2, is strongly downregulated on both mRNA and protein level in HCT116 5/4 (Fig. 4b). HMGA2 is a validated target of hsa-miR-26a-5p, let-7a-5p and hsa-miR-125b-5p, all of which were found to affect cellular growth and proliferation. Downregulation of HMGA2 by hsa-miR-26a negatively affects cell proliferation and has tumour suppressive effects in gall bladder cancer . Another strongly downregulated target is RCC2, which is targeted by hsa-miR-192-5p and hsa-miR-7-5p in HCT1165/4 and RPE1 5/3 12/3 (Fig. 4b, c). RCC2 is essential for mitotic spindle assembly and recently a novel function of RCC2 for G1-S transition has been recognized [46,47,48]. Aneuploidy impairs proliferation in most model systems [12, 13, 34] and this is reflected in the conserved transcriptional changes that include the downregulation of cell cycle associated pathways [14, 15]. Molecular mechanisms that lead to down-regulation of pro-proliferative genes and impaired proliferation in response to chromosome gain are poorly understood. Based on our results we propose that the miRNA deregulation contributes to the gene expression changes and cell cycle impairments observed in aneuploid cells.
Only a few miRNAs were commonly upregulated in response to aneuploidy, such as miR-139-5p or miR-10a-5p. Strikingly, many of these upregulated miRNAs are often downregulated in cancers. For example, hsa-miR-139-5p, a recognized tumor-suppressing miRNA, is downregulated in non-small cell lung cancer, gastric cancer, breast cancer and colorectal carcinoma and often associates with poor prognosis . We found that inhibition of miR-139-5p decreases the expression of cell cycle inhibitors CDKN1A (p21) and CDKN1B (p27) (Fig. 3c). The opposite direction of expression changes in model aneuploids was previously observed for changes of mRNA levels, where the transcriptional response to chromosome gain is inversed to the cancer transcriptome . Since gain of chromosomes has a potent tumor-suppressing effect , our findings suggest that the miRNAs deregulation contributes to this phenotype.
Among the few commonly upregulated miRNAs, we found hsa-miR-10a-5p and confirmed its increased endogenous activity using a specific luciferase-based reporter system (Fig. 6a, b). Only a few hsa-miR-10a-5p targets showed a concordant downregulation on protein or mRNA level in human aneuploid model cell lines. These few affected targets were involved in diverse cellular functions such as cytoskeleton stabilization, transcriptional regulation and metabolic processes (Additional file 1: Figure S5, Additional file 4: Table S3). This finding prompted us to investigate another possible function of hsa-miR-10a-5p. Hsa-miR-10a-5p binds upstream of the mRNA 5’TOP motif and affects the translation of mRNAs carrying this motif. The 5’TOP motif is a cis-regulatory motif present in transcripts encoding for the translational machinery, mostly ribosomal proteins [40, 51].
What is the function of the upregulated hsa-miR-10a-5p in aneuploid cells? Upon nutrition deprivation or cell cycle arrest the translation of 5’ TOP mRNAs is repressed, which is indicated by their shift to inactive ribonucleoprotein complexes . We observed that the 5’TOP motif luciferase activity is not repressed upon starvation in aneuploid cells, in contrast to the diploid parental cell lines, where the 5’TOP motif luciferase activity is sensitive to starvation (Fig. 6c-e). As aneuploidy results in a plethora of cellular stresses, such as proteotoxic stress, metabolic and replication stress , we propose that the upregulation of hsa-miR-10a-5p in aneuploid cell lines might be a protective adaption to adverse stress conditions.
Our unique data set describes the aneuploid cells from an “omics” perspective and allows a global, quantitative analysis of the consequences of aneuploidy in different cell lines. Our system demonstrates that an addition of extra chromosomes largely affects cell physiology on multiple levels, including marked deregulation of the miRNAome. The altered miRNAome landscape may explain some of the transcriptome and proteome changes that occur in response to aneuploidy and generally emphasizes the detrimental consequences of chromosome gains.
The retinal pigment epithelial cell line RPE1 hTERT and RPE1 hTERT H2B-GFP were a gift from Stefan Taylor (University of Manchester, UK). The human colorectal cell line HCT116 was obtained from American Type Culture Collection (no. CCL-247). HCT116 H2B-GFP was generated previously by lipofection (FugeneHD, Roche) transfection of pBOS-H2B-GFP (BD Pharmingen) according to manufacturer’s protocols . The tetrasomic cell line HCT116 5/4 were a kind gift of Minoru Koi (Baylor University Medical Centre, Dallas, Texas, USA). All other trisomic and tetrasomic cell lines were generated by microcell-mediated chromosome transfer as described previously [5, 9, 12]. Cell lines originating from individual clonal populations arising from a single cell after by microcell-mediated chromosome transfer are indicated with their clone number (c#). The genome sequencing data of all cell lines are available in the public repository of NCBI under the accession number GSE102855. Summary of all cell line information is in Table 4.
RNA sequencing and data processing
Following aneuploid cell lines and parental diploid cell lines were subjected to RNA sequencing: HCT116, HCT116 H2B-GFP (cell lines containing H2B-GFP will be indicated by * in the cell line name), HCT116 5/4, HCT116* 8/3 clone7, RPE1, RPE*, RPE1 5/3 12/3, RPE1* 21/3.
Total RNA was extracted using RNeasy Mini Kit (QIAGEN). TruSeq RNA library preparation and Illumina HiSeq2500 sequencing with 25 million 100 bp single reads per library were performed by the Max Planck-Genome-Center Cologne, Germany. Subsequently, sequencing adapters were removed from the raw sequences with cutadapt and sequencing reads were mapped to the human genome (hg19) using TopHat (v2.0.10) with the following parameters: “tophat2 -g1 -G”. RefSeq information in the GTF file was downloaded from the UCSC genome browser. featureCounts (v1.4.3) was used to generate the count matrix with the same GTF file as for the alignment with the following parameters: “-t exon -g gene_id”. Normalization and differential expression analysis was performed using the R/Bioconductor package DESeq2 . For differential expression analysis, trisomic and tetrasomic cell lines were compared to the parental diploid cell line. The NGS data discussed in this publication have been deposited in NCBI’s Gene Expression Omnibus and are accessible through GEO Series accession number GSE102855 (https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/geo/query/acc.cgi?acc=GSE102855).”
Small RNA sequencing and data processing
Following aneuploid cell lines and parental diploid cell lines were subjected to small RNA sequencing: HCT116, HCT116*, HCT116 5/4, HCT116 18/3 clone2, HCT116* 8/3 clone7, HCT116* 3/3, RPE1, RPE*, RPE1 5/3 12/3, RPE1* 21/3 and RPE1* 7/3.
Total RNA, including small RNAs from 18 nucleotides upwards, was extracted using miRNeasy Mini Kit (QIAGEN). TruSeq small RNA library preparation and Illumina HiSeq2500 sequencing with 25 million 100 bp or 75 bp single reads per library were performed by the Max Planck-Genome-Center Cologne, Germany. If necessary, raw sequencing read ends were trimmed for low quality base calls using the FASTQ quality trimmer with the following parameters: “-Q33 -t 20 -l 17” [28, 32]. Sequencing adapters were removed with cutadapt. Mapping to the human genome hg19 and miRNA identification was performed using mirdeep2 (v126.96.36.199)  with the following commands: “mapper.pl -e -h -i -j -l 18 -m -p hg19 -q” and “miRDeep2.pl reads.fa hg19.fa genome.arf human_mature_rmspace.fa others_mature_rmspace.fa human_hairpin_rmspace.fa -d -t Human”. The three biological replicates of aneuploid and the corresponding parental cell line were analyzed within the same analysis run. For subsequent normalization and differential expression analysis miRNA_expressed_all_samples.csv was used as an input for DESeq2. Differential expression analysis was performed between aneuploid cell lines and corresponding parental diploid cell line. To account for sequencing batch effects, paired replicate information in addition to condition information was input to the DESeq2 analysis for HCT116 5/4, RPE1 5/3 12/3 and RPE1* 21/3 and corresponding parental cell lines. The genome sequencing data of all cell lines are available in the public repository of NCBI under the accession number GSE102855.
Integrative mRNA, miRNA and target analysis
Data analysis, such as integration and data visualization was performed using the computing environment R. Mapping of identifiers as well as genome locations was performed with BioMart using biomaRt package in R [54, 55]. miRNA (miRNA) identifiers were retrieved from miRBase (v21). miRNA target information was retrieved from miRTarBase (v6.1) [36, 56].
Functional annotation analysis
miRNA differential expression datasets were analyzed using the Ingenuity Pathway Analysis Software . Core analysis was performed with the differential expressed miRNAs (log2 fold change +/− 0.6, adjusted p-value < 0.05) against the Ingenuity knowledge base considering direct experimentally observed relationships in human species. Functional annotation analysis results were exported and visualized in R. Functional annotation analysis of target genes was conducted using the Database for Annotation, Visualization and Integrated Discovery (DAVID v6.8) by applying functional annotation clustering.
Quantitative real-time PCR of miRNAs
Total RNA, including small RNAs from 18 nucleotides upwards, was extracted using miRNeasy Mini Kit (QIAGEN). Reverse transcription was performed using the universal cDNA synthesis kit miRCURY LNA™ (Exiqon) according to manufacturers protocol. UniSp6 RNA spike-in control was added to each sample in equal amounts. Quantitative PCR was conducted using the Light Cycler 480 System (Roche Diagnostics) using the ExiLENT SYBR® Green master mix miRCURY LNA™ (Exiqon) and miRNA specific LNA™ PCR primer sets (Exiqon) as well as the UniSp6 RNA spike control primer set (Exiqon). Absolute quantification with an external standard was performed and negative non-template controls were included in all experiments. The specificity of the primer product amplification was confirmed in each run by melting curve analysis. miRNA expression was normalized to the control spike RNA and corresponding diploid miRNA expression.
hsa-miR-10a-5p 3’UTR luciferase reporter (psiCheck2-10a) was constructed by cloning the complementary miRNA target sequence in the synthetic 3’UTR of psiCHECK™-2 Vector (Promega). PCR cloning was performed with the following primers:
top 5’-CACAAATTCGGATCTACAGGGTAGTTTAAACCTAGAGCGGCCGCT-‘3 and.
pS16-WT luciferase vector containing the transcriptional start site of the ribosomal protein S16 (Rps16) and 29 nt of the exon 1 including the 5’TOP motif was a kind gift of Anders H. Lund (University of Copenhagen, Denmark). pRL-TK Renilla vector was a kind gift of Reinhard Fässler (Max Planck Institute of Biochemistry, Martinsried).
Endogenous hsa-miR-10a-5p mediated 3’UTR repression assay
Cell lines were forward transfected with psiCheck2-10a using Lipofectamine 200 (Thermo Fisher Scientific) according to manufacturer’s protocol. Cells were reseeded into 96 well plates 24 h (hrs) post transfection at 20.000/ well HCT116- derived cell lines and 10.000/ well RPE1 derived cell lines. Luciferase activity was monitored 48 h post transfection using the Dual-Glo® Luciferase Assay System (Promega). Renilla luciferase activity values were normalized to the Firefly luciferase activity and subsequently to the parental cell line. Statistical testing and data plotting was performed in GraphPad Prism 6.
Hsa-miR-10a-5p mediated 5’TOP motif mRNA translation assay
Cell lines were reverse transfected with miRCURY LNA™ miRNA Mimic or miRCURY LNA™ miRNA Mimic Negative Control (Exiqon) at 50 nM using Lipofectamine 200 (Thermo Fisher Scientific) according to manufacturers protocol. Cells were forward transfected with the pS16-wt-luc Firefly luciferase reporter construct and pRL-TK Renilla luciferase control vector 24 h post mimic transfection. Forty-eight hours post mimic transfection, cells were reseeded into 96 well plates at 20.000/ well HCT116- derived cell lines and 10.000/ well RPE1- derived cell lines. Luciferase activity was monitored 72 h post transfection using Dual-Glo® Luciferase Assay System (Promega). Starvation was performed by replacing the Dulbecco’s Modified Eagle Medium (DMEM) supplemented with 10% fetal bovine serum and 1% penicillin/ streptomycin with DMEM without supplements for 3 h prior measurement of luciferase activity. Firefly luciferase activity values were normalized to Renilla luciferase activity and subsequently to the parental cell line. Statistical testing and data plotting was performed in GraphPad Prism 6.
3′ untranslated region
5′ untranslated region
Ingenuity pathway analysis
micro ribonucleic acid
Polymerase chain reaction
Colnaghi R, Carpenter G, Volker M, O’driscoll M. The consequences of structural genomic alterations in humans: genomic disorders, genomic instability and cancer. Semin Cell Dev Biol. 2011;22:875–85.
Mitelman F, Johansson B, F M: Mitelman Database of Chromosome Aberrations and Gene Fusions in Cancer In: http://cgap.nci.nih.gov/Chromosomes/Mitelman;. 2016.
Weaver BA, Cleveland DW. Does aneuploidy cause cancer? Curr Opin Cell Biol. 2006;18:658–67.
Ariyoshi K, Miura T, Kasai K, Fujishima Y, Oshimura M, Yoshida MA. Induction of genomic instability and activation of autophagy in artificial human aneuploid cells. Mutat Res. 2016;790:19–30.
Donnelly N, Passerini V, Dürrbaum M, Stingele S, Storchova Z. HSF1 deficiency and impaired HSP90-dependent protein folding are hallmarks of aneuploid human cells. EMBO J. 2014;33:2374–87.
Nawata H, Kashino G, Tano K, Daino K, Shimada Y, Kugoh H, Oshimura M, Watanabe M. Dysregulation of gene expression in the artificial human trisomy cells of chromosome 8 associated with transformed cell phenotypes. PLoS One. 2011;6:e25319.
Niwa O, Tange Y, Kurabayashi A. Growth arrest and chromosome instability in aneuploid yeast. Yeast. 2006;23:937–50.
Oromendia AB, Dodgson SE, Amon A. Aneuploidy causes proteotoxic stress in yeast. Genes Dev. 2012;26:2696–708.
Passerini V, Ozeri-Galai E, de Pagter MS, Donnelly N, Schmalbrock S, Kloosterman WP, Kerem B, Storchova Z. The presence of extra chromosomes leads to genomic instability. Nat Commun. 2016;7:10754.
Pavelka N, Rancati G, Zhu J, Bradford WD, Saraf A, Florens L, Sanderson BW, Hattem GL, Li R: Aneuploidy confers quantitative proteome changes and phenotypic variation in budding yeast. In: Nature. vol. 468: Nature Publishing Group; 2010: 321–325.
Stingele S, Stoehr G, Storchova Z. Activation of autophagy in cells with abnormal karyotype. Autophagy. 2013;9:246–8.
Stingele S, Stoehr G, Peplowska K, Jur C, Mann M, Storchova Z. Global analysis of genome, transcriptome and proteome reveals the response to aneuploidy in human cells. Mol Syst Biol. 2012;8:1–12.
Torres EM, Sokolsky T, Tucker CM, Chan LY, Boselli M, Dunham MJ, Amon A. Effects of aneuploidy on cellular physiology and cell division in haploid yeast. Science. 2007;317:916–24.
Sheltzer JM, Torres EM, Dunham MJ, Amon A. Transcriptional consequences of aneuploidy. Proc Natl Acad Sci. 2012;109:12644–9.
Dürrbaum M, Kuznetsova AY, Passerini V, Stingele S, Stoehr G, Storchova Z. Unique features of the transcriptional response to model aneuploidy in human cells. BMC Genomics. 2014;15:1–14.
Friedman RC, Farh KK-H, Burge CB, Bartel DP. Most mammalian mRNAs are conserved targets of microRNAs. Genome Res. 2009;19:92–105.
Liu C, Mallick B, Long D, Rennie WA, Wolenc A, Carmack CS, Ding Y. CLIP-based prediction of mammalian microRNA binding sites. Nucleic Acids Res. 2013;41:e138.
Schnall-Levin M, Rissland OS, Johnston WK, Perrimon N, Bartel DP, Berger B. Unusually effective microRNA targeting within repeat-rich coding regions of mammalian mRNAs. Genome Res. 2011;21:1395–403.
Ørom UA, Nielsen FC, Lund AH. MicroRNA-10a binds the 5′UTR of ribosomal protein mRNAs and enhances their translation. Mol Cell. 2008;30:460–71.
Selbach M, Schwanhäusser B, Thierfelder N, Fang Z, Khanin R, Rajewsky N. Widespread changes in protein synthesis induced by microRNAs. Nature. 2008;455:58–63.
Shenoy A, Blelloch RH: Regulation of microRNA function in somatic stem cell proliferation and differentiation. In: Nature Reviews Molecular Cell Biology. vol. 15: Nature Publishing Group; 2014: 565–576.
Hwang H-W, Mendell JT. MicroRNAs in cell proliferation, cell death, and tumorigenesis. Br J Cancer. 2006;94:776–80.
Leung AKL, Sharp PA: MicroRNA Functions in Stress Responses. In: Mol Cell. vol. 40: Elsevier Inc.; 2010: 205–215.
Kurozumi A, Goto Y, Okato A, Ichikawa T, Seki N: Aberrantly expressed microRNAs in bladder cancer and renal cell carcinoma. In.: Nature Publishing Group; 2016: 1–8.
Thomas J, Ohtsuka M, Pichler M, Ling H. MicroRNAs: clinical relevance in colorectal cancer. IJMS. 2015;16:28063–76.
Yeh C-H, Moles R, Nicot C. Clinical significance of microRNAs in chronic and acute human leukemia. Mol Cancer. 2016;15:37.
Wilting SM, Snijders PJF, Verlaat W, Jaspers A, van de Wiel MA, van Wieringen WN, Meijer GA, Kenter GG, Yi Y, le Sage C et al: Altered microRNA expression associated with chromosomal changes contributes to cervical carcinogenesis. In: Oncogene. vol. 32: Nature Publishing Group; 2012: 106–116.
Wang Y, Hu X, Greshock J, Shen L, Yang X, Shao Z, Liang S, Tanyi JL, Sood AK, Zhang L. Genomic DNA copy-number alterations of the let-7 family in human cancers. PLoS One. 2012;7:e44399.
Choi YE, Pan Y, Park E, Konstantinopoulos P, De S, D’Andrea A, Chowdhury D. MicroRNAs down-regulate homologous recombination in the G1 phase of cycling cells to maintain genomic stability. elife. 2014;3:e02445.
Hell MP, Thoma CR, Fankhauser N, Christinat Y, Weber TC, Krek W. miR-28-5p promotes chromosomal instability in VHL-associated cancers by inhibiting Mad2 translation. Cancer Res. 2014;74:2432–43.
Kozomara A, Griffiths-Jones S. miRBase: annotating high confidence microRNAs using deep sequencing data. Nucleic Acids Res. 2013;42:D68–73.
Friedländer M, Mackowiak SD, Li N, Chen W, Rajewsky N. miRDeep2. Nucleic Acids Res. 2011;40:1–16.
Love MI, Wolfgang H, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:31.
Williams BR, Prabhu VR, Hunter KE, Glazier CM, Whittaker CA, Housman DE, Amon A. Aneuploidy affects proliferation and spontaneous immortalization in mammalian cells. Science. 2008;322:703–9.
Upender MB, Habermann JK, McShane LM, Korn EL, Barrett JC, Difilippantonio MJ, Ried T. Chromosome transfer induced aneuploidy results in complex dysregulation of the cellular transcriptome in immortalized and cancer cells. Cancer Res. 2004;64:6941–9.
Vergoulis T, Vlachos IS, Alexiou P, Georgakilas G, Maragkakis M, Reczko M, Gerangelos S, Koziris N, Dalamagas T, Hatzigeorgiou AG. TarBase 6.0: capturing the exponential growth of miRNA targets with experimental support. Nucleic Acids Res. 2012;40:D222–9.
Quiagen: Ingenuity Pathway Analysis Software http://www.ingenuity.com.
Kraemer A, Green J, Pollard JJ, Tugendreich S. Causal analysis approaches in ingenuity pathway analysis. Bioinformatics. 2014;30:523–30.
Wong CC-L, Wong C-M, Tung EK-K, Au SL-K, Lee JM-F, Poon RT-P, Man K, Ng IO-L. The microRNA miR-139 suppresses metastasis and progression of hepatocellular carcinoma by down-regulating rho-kinase 2. Gastroenterology. 2011;140:322–31.
Meyuhas O, Kahan T: The race to decipher the top secrets of TOP mRNAs. In: BBA - Gene Regulatory Mechanisms. vol. 1849: Elsevier B.V.; 2015: 801–811.
Czubak K, Lewandowska MA, Klonowska K, Roszkowski K, Kowalewski J, Figlerowicz M, Kozlowski P. High copy number variation of cancer-related microRNA genes and frequent amplification of DICER1 and DROSHA in lung cancer. Oncotarget. 2015;6:23399–416.
Zhang L, Huang J, Yang N, Greshock J, Megraw MS, Giannakakis A, Liang S, Naylor TL, Barchetti A, Ward MR, et al. microRNAs exhibit high frequency genomic alterations in human cancer. Proc Natl Acad Sci U S A. 2006;103:9136–41.
Ramsingh G, Jacoby MA, Shao J, De Jesus Pizzaro RE, Shen D, Trissal M, Getz AH, Ley TJ, Walter MJ, Link DC. Acquired copy number alterations of miRNA genes in acute myeloid leukemia are uncommon. Blood. 2013;122(15):e44–51.
Veigaard C, Kjeldsen E. Exploring the genome-wide relation between copy number status and microRNA expression. Genomics. 2014;104(4):271–8.
Zhou H, Guo W, Zhao Y, Wang Y, Zha R, Ding J, Liang L, Hu J, Shen H, Chen Z, et al. MicroRNA-26a acts as a tumor suppressor inhibiting gallbladder cancer cell proliferation by directly targeting HMGA2. Int J Oncol. 2014;44:2050–8.
Mollinari C, Reynaud C, Martineau-Thuillier S, Monier S, Kieffer S, Garin J, Andreassen PR, Boulet A, Goud B, Kleman J-P, et al. 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.
Rosasco-Nitcher SE, Lan W, Khorasanizadeh S, Stukenberg PT. Centromeric aurora-B activation requires TD-60, microtubules, and substrate priming phosphorylation. Science. 2008;319:469–72.
Yenjerla M, Panopoulos A, Reynaud C, Fotedar R, Margolis RL: TD-60 is required for interphase cell cycle progression. In: cc. vol. 12; 2013: 837–841.
Sheltzer JM, Transcriptional A. Metabolic signature of primary aneuploidy is present in chromosomally unstable cancer cells and informs clinical prognosis. Cancer Res. 2013;73:6401–12.
Sheltzer JM, Ko JH, Replogle JM, Habibe Burgos NC, Chung ES, Meehl CM, Sayles NM, Passerini V, Storchova Z, Amon A. Single-chromosome gains commonly function as tumor suppressors. Cancer Cell. 2017;31(2):240–55.
Levy S, Avni D, Hariharan N, Perry RP, Meyuhas O. Oligopyrimidine tract at the 5' end of mammalian ribosomal protein mRNAs is required for their translational control. Proc Natl Acad Sci U S A. 1991;88:3319–23.
Santaguida S, Amon A. Short- and long-term effects of chromosome mis-segregation and aneuploidy. Nat Rev Mol Cell Biol. 2015;16:473–85.
Kuffer C, Kuznetsova AY, Storchova Z. Abnormal mitosis triggers p53-dependent cell cycle arrest in human tetraploid cells. Chromosoma. 2013;122:305–18.
Durinck S, Spellman PT, Birney E, Wolfgang H. Mapping identifiers for the integration of genomic datasets with the R/Bioconductor package biomaRt. Nat Protoc. 2009;4:1184–91.
Smedley D, Haider S, Durinck S, Pandini L, Provero P, Allen J, Arnaiz O, Awedh MH, Baldock R, Barbiera G, et al. The BioMart community portal: an innovative alternative to large, centralized data repositories. Nucleic Acids Res. 2015;43:W589–98.
Chou C-H, Chang N-W, Shrestha S, Hsu S-D, Lin Y-L, Lee W-H, Yang C-D, Hong H-C, Wei T-Y, Tu S-J, et al. miRTarBase 2016: updates to the experimentally validated miRNA-target interactions database. Nucleic Acids Res. 2016;44:D239–47.
We are thankful to Aline Sewo de Campos and Isabell Kirchner for their technical support, as well as the members of the Storchova and Habermann groups for scientific discussions and suggestions.
This work was supported by a DFG grant to Z. S. (STO918/2–2) and by the Cluster of Excellence – Center for Integrative Protein Science Munich (CIPSM).
Availability of data and materials
All the data and material are available in the NCBI public repository under the accession number GSE102855 (https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/geo/query/acc.cgi?acc=GSE102855) or upon request.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. Correlation of extra DNA and the percentage of deregulated microRNAs for each aneuploid cell line. Figure S2. microRNA expression aligned according to their chromosome position. Summary of all microRNA expression in analyzed aneuploid cell lines. Figure S3. microRNAs and their targets deregulated in aneuploids. Graphic presentation of the miRNAs deregulation relative to their targets in aneuploid human cells. Figure S4. Comparison of hsa-miR-10a-5p expression and HOXB3 mRNA expression. The correlation between expression of hsa-miR-10a-5p and its hosting genomic sequence. Figure S5. Target mRNA and protein expression levels. Graphic presentation of the hsa-miR-10a deregulation relative to its targets in aneuploid human cells. Figure S6. Overexpression of hsa-miR-10a-5p leads to repression of luciferase activity. Validation of the hsa-miR-10a-5p overexpression in human cells. (PDF 250 kb)
Table S1. Total number and number of deregulated microRNAs per chromosome (log2 fold change > 0.6 or <(− 0.6)). Overview of all deregulated microRNAs. (DOCX 16 kb)
Table S2. Functional annotation clustering. Complete results of the DAVID analysis of functional annotation clustering of microRNAs targets deregulation. (XLSX 114 kb)
Table S3. Functional annotation clustering of hsa-miR-10a-5p targets. Complete results of the DAVID analysis of functional annotation clustering of hsa-miR-10a-5p target deregulation. (XLSX 47 kb)
About this article
Cite this article
Dürrbaum, M., Kruse, C., Nieken, K.J. et al. The deregulated microRNAome contributes to the cellular response to aneuploidy. BMC Genomics 19, 197 (2018) doi:10.1186/s12864-018-4556-6