- Open Access
Effect of BRAFV600E mutation on transcription and post-transcriptional regulation in a papillary thyroid carcinoma model
Molecular Cancer volume 6, Article number: 21 (2007)
microRNAs (miRNAs) are a group of non-coding single stranded RNAs measuring approximately 22 nucleotides in length that have been found to control cell growth, differentiation and apoptosis. They negatively regulate target genes and have recently been implicated in tumourigenesis. Furthermore, miRNA expression profiling correlates with various cancers, with these genes thought to act as both tumour suppressors and oncogenes. Recently, a point mutation in the BRAF gene leading to a V600E substitution has been identified as the most common genetic change in papillary thyroid carcinoma (PTC) occurring in 29–69% of cases. This mutation leads to aberrant MAPK activation that is implicated in tumourigenesis.
The aim of this study was to identify the effect that BRAF oncogene has on post-transcriptional regulation in PTC by using microRNA analysis.
A unique miRNA expression signature differentiated between PTC cell lines with BRAF mutations and a normal thyroid cell line. 15 miRNAs were found to be upregulated and 23 miRNAs were downregulated. Several of these up/down regulated miRNAs may be involved in PTC pathogenesis. miRNA profiling will assist in the elucidation of disease pathogenesis and identification biomarkers and targets.
Papillary thyroid carcinoma (PTC) is the most commonly occurring thyroid malignancy accounting for approximately 80% of cases. Genetically PTC is defined by several alterations which cause abnormal activation of the mitogen-activated protein kinase (MAPK) pathway, the most prevalent being point mutations in the intracellular signalling kinase BRAF [1, 2]. An activating mutation in BRAF (V600E) leads to constitutive activation and hence de-regulation of MAPK pathway. Several groups have reported that PTCs harbouring mutations in BRAF have more aggressive properties, present more often with extra-thyroidal invasion and also at a more advanced stage. Histological features such as the classic PTC architecture and tall cell features are also associated with the BRAF oncogene [3, 4].
Despite this range of knowledge thyroid carcinoma still poses a significant diagnostic challenge and understanding of the complete range of genes and pathways involved in its pathogenesis is incomplete. Until recently biomarker discovery has concentrated largely on protein coding RNAs in a quest to identify cancer associated genes involved in thyroid carcinogenesis [5–7]. The recent discovery however that a large group of non-coding RNAs called microRNAs (miRNAs) may potentially play a role in cancer has led to a proliferation of miRNA related research of late including the publication of several key reviews [8–12]. These small non-coding RNAs constitute a novel class of gene regulators that function by negatively regulating gene expression by targeting mRNAs for cleavage or translational repression. Several recent studies on the regulatory functions of miRNAs suggest that they play critical roles in central processes such as development, cell proliferation and differentiation, stress resistance, metabolism and apoptosis all of which are involved in tumourigenesis [13, 15]. Not surprisingly aberrant miRNA expression has been implicated in cancer development and a hypothesis has emerged suggesting that tumours may each have a discrete "miRNA signature" [16–20]. Differential expression of miRNAs between malignant tissue and normal tissue and between different types of tumour, has been shown in several key studies, indicating that miRNAs are determinants of clinical diagnostic and prognostic significance [21, 22].
Recently several groups have utilised a range of techniques most notably northern blot analysis, cloning and oligonucleotide microarrays, in order to analyse miRNA expression profiles in different cancer types [23–26]. In this study a novel stem-looped TaqMan® RT-PCR method was used to measure the expression levels of a panel of 160 mature miRNAs in thyroid cell lines containing BRAF mutation.
Differential miRNA expression using TaqMan™ microRNA assays
Stem-looped TaqMan® RT-PCR was used to measure the expression levels of 160 miRNAs in a collection of cell lines including a normal thyroid cell lines N-thy-ori 3-1 and two cell lines containing V600E BRAF point mutation [Nthy-oriBRAF and KAT-10]. A comparison was then made between Nthy-ori and BRAF mutated cell lines to determine miRNAs whose expression was differentially expressed between the two groups. ΔΔCT method was used to calculate fold change of miRNA expression between groups. 15 miRNAs were found to be upregulated with a fold change of > 2 in BRAF mutated cell lines when compared to normal.
Three of these upregulated miRNAs showed significantly higher fold change than the other upregulated miRNAs, these are mir-200a, mir-200b and mir-141.
23 miRNAs showed underexpression of > 2 fold when compared to the normal thyroid cell line. Moreover, three miRNAs showed significantly more underexpression compared to the other downregulated miRNAs. These miRNAs are as follows mir-127, mir-130a and mir-144. (Table 1 and Figures 1 and 2).
mRNA (transcription) and miRNA (regulation) expression correlation
miRNAs interact and regulate their target genes by cleavage or translational inhibition . A current challenge is to identify biologically relevant targets that are regulated by individual miRNAs. This task is further complicated by the fact that each miRNA can potentially bind and regulate many mRNA targets and each mRNA can be bound and regulated by several miRNAs. Several publicly available databases exist for prediction of miRNA targets. MiRBase, PICTAR and TARGETSCAN were used in this study [26–28].
These three databases were utilized to yield lists of putative targets for each miRNA which are illustrated in table 2. From these lists there are several genes that are biologically relevant in the context of PTC. The miRNAs which showed the highest increase in fold change are predicted to bind genes involved in thyroid function, MAPKKK cascade, retinoic acid receptors, cell adhesion and cell structure associated genes and genes involved in cell signalling. miRNAs which showed the largest decrease in fold change are predicted to regulate and target genes involved in G-protein mediated signalling, oncogenesis and cell cycle control. One might predict that a miRNA which is overexpressed would be associated with downstream underexpression of its target. This regulation can occur either at the transcript level or the protein level.
microRNAs (miRNAs) are now recognized as an additional layer of post-transcriptional control that must be considered in order to understand the complexity of gene expression in humans. miRNAs are now known as a novel class of gene regulators that are involved in cancer related processes and deregulation or aberrant expression of miRNAs may contribute to human diseases, including cancer. Calin et al showed down-regulation of miRNA-15 and miRNA-16 in a majority of chronic lymphatic leukaemia (CLL). Altered miRNA expression has also been reported in breast cancer, glioblastoma, lung cancer and colorectal cancer [29–35].
This study utilized a novel stem-looped TaqMan RT-PCR method to measure the expression levels of 160 mature miRNAs in a collection of thyroid cell lines. This method offers several distinct advantages over conventional miRNA detection methods including a higher sensitivity and specificity and a fast and simple methodology. The aim of this study was to correlate BRAF mutant status with an associated miRNA profile, with a view to elaborating the downstream regulatory effect of BRAF mutation.
The data showed a set of differentially expressed miRNAs in cell lines with BRAF V600E mutation compared to a normal thyroid cell line (Table 1). There were 15 overexpressed miRNAs and 23 underexpressed miRNAs in these cell lines compared to the normal cell line. Of these miRNAs listed in table 1 only a few have been previously associated with other cancers. One of which is miR-127, that showed down-regulation in BRAF V600E mutated cell lines. Human miR-127 is embedded in a CpG island and is silenced in several human cancers including bladder, prostate and colon cancer and has been suggested to play a role as a possible tumour suppressor gene . miR-127 also an example of the potential utility of miRNA targeted therapy in cancer. Saito et al demonstrated that induction of miR-127 by with 5-aza-20-deoxycytidine and 4-phenylbutyric acid, which inhibits DNA methylation and histone deacetylase respectively, reduces expression of the oncogene BCL6 in bladder cancers . Its utility as a potential cancer target has yet to be explored in other cancers however this suggests the possibility of direct targeting of miRNAs that are amplified or upregulated in patient tumours. Two other miRNAs listed in table 2 have also been described in other cancers. These include miR-200b and miR-141 which have been shown to be highly overexpressed in malignant cholangiocytes and in colon carcinoma. . This suggests the possibility of miR-200b and miR-141 down-regulating possible tumour suppressor genes. Of particular interest in this study is the downregulation of miR-323 and miR-302b in BRAF mutated cell line. miRBASE target database predicted that the BRAF transcript has binding sites for miR-323 and miR-302b to potentially bind and down- regulate BRAF expression.
A number of miRNAs exhibited a difference in expression level between KAT10 cell line and normal but did not show a significant change in the BRAF transfected cell lines and normal. miRNAs specific to the PTC cell line with native mutant BRAF included, mir-10a [overexpressed with a fold change of 95 in KAT10] mir-23b [downregulated with a fold change of 2949] and mir-154 [downregulated with a fold change of 142]. It is possible these miRNAs may play a role in a later stage in tumour formation and progression; given the transfected cell line had only been exposed to mutant BRAF for a short period (3 passages) unlike KAT10 which represents a genuine PTC which has evolved with BRAT mutant as an intrinsic component. Alternatively, the difference may reflect the accumulation of several genetic insults in the process of tumourigenesis.
Recent studies in miRNA deregulation PTC found an aberrant miRNA expression profile in PTCs compared to normal thyroid tissues [36, 37]. In particular, a significant increase in mir-222, mir-221, mir-146. Our study did not find any of these miRNAs to be significantly overexpressed in the cell lines examined. miR-222 was actually found to be downregulated in the BRAF mutated cell lines. Another miRNA these two groups found to be upregulated albeit less significantly was miR-213, our study also found miR-213 to be upregulated in BRAF mutated cell lines [25 fold]. These differences may be due to the different experimental approaches; Huiling et al and Pallante et al used miRNA microarrays. They also used different chemistry and fresh tissue samples; also the RET and BRAF status of each PTC case was not disclosed.
Prediction of miRNA targets is an immediate challenge in miRNA research. Several computational approaches have been implemented for miRNA gene prediction using methods based on sequence conservation and/or structural similarity. To overcome the limitations of target prediction three of these databases MiRBase, TARGETSCAN and PICTAR were used in this study to predict possible targets of the differentially expressed miRNAs associated with BRAF mutated (Table 2). Examples of other programs developed for target identification include RNAhybrid, RNA calibrate and RNAeffective. Only a few target mRNAs have been experimentally proven and studied in vitro and not all the prediction algorithms yield overlapping lists therefore the lists may contain several false positives that may fit the criteria but do not interact with the miRNA in vitro. To overcome this, results were intersected from these databases to identify the genes commonly predicted by at least two of the algorithms. The results of the potential targets of several differentially expressed miRNAs are displayed in table 2.
It is reasonable to expect that targets of downregulated miRNAs include oncogenes or genes encoding proteins with potential oncogenic function. Indeed among putative targets, several genes with potential oncogenic functions were found include RAB35, ERBB4, MAP3K3, VASP6, KIT, FOSB, NOTCH1. Upregulated miRNAs are purported to interact with tumour suppressor genes. Those identified from this study include MAD 4 and ST6. Several extracellular membrane molecules were also identified including; myotubulin, MAD1A, Cadherins. In general overexpression of an miRNA would be associated with downregulation of its target miRNA.
Interestingly among these potential gene targets are a number of genes that have implications in thyroid carcinoma progression such as genes involved in thyroid function, cell cycle control, cell signalling, MAPKKK cascade, retinoic acid receptors, cell adhesion and cell structure associated genes. Therefore, it seems plausible that this subset of miRNAs may be important in the progression of thyroid carcinoma by working together or independently to regulate genes involved in thyroid tumourigenesis
Our group has also previously examined gene expression levels in these same cell lines. This analysis was done by evaluating genome-wide gene expression levels in BRAF mutated cell lines and matched normal thyroid cell lines using the Applied Biosystems Gene Expression Arrays. Treatment comparisons were preformed resulting in lists of differentially expressed genes (Figure 3). These gene lists were uploaded into PANTHER (Protein ANalysis THrough Evolutionary Relationships) classification system and the software statistically compared them to a reference list to look for under and over-represented molecular functions and biological process. Analysis revealed a set of genes associated with mutant BRAF expression which are involved in processes such as apoptosis, chromosomal instability, invasiveness and metastasis and proliferation. Other genes identified were shown to be involved in the insulin-like growth factor pathway (including (insulin-like growth factor binding protein) IGFBP4, (insulin receptor substrate) IRSI, GRB14, PLAG1 and IGF2) and the Wnt signaling pathway suggesting their importance in BRAF signalling.
Up and down regulation of cell adhesion and chromosome remodelling molecules were also shown. BRAF mutation has previously been shown to induce chromosome instability and decrease cell- cell interactions  and our data further corroborates the involvement of these processes in papillary thyroid carcinoma. Not surprisingly genes involved in the MAPK pathway, of which BRAF is an integral part, showed an increase in expression. Among these genes were DDIT, DUSP6, TFPI-2, Seladin-1, MAP3K12, and GRB14. These results provided an insight into the different pathways through which BRAF exerts its oncogenic effect and emphasised the importance of the MAPK and Wnt signalling pathways in thyroid tumour formation.
The list of putative miRNA targets (as defined above in table table 1) correlated with the list of differentially expressed genes on analysis. This is illustrated in table 3. It is interesting to speculate that these genes are regulated by their corresponding miRNA thereby contributing to tumour development in thyroid. These observations should be placed in context, in so far as predictive databases are incomplete and target selection is essentially speculative and involves a degree of subjectivity as these target mRNAs have not been experimentally proven and studied in vitro.
In summary this study revealed a set of miRNAs associated with expression of mutant BRAFV600E oncogene in PTC cell lines when compared to a normal thyroid cell line. The BRAF oncogene may influence miRNA expression directly or indirectly by altering the expression patterns of genes, such as transcription factors, involved in tumour progression or indeed the genes involved in the biogenesis and synthesis of miRNAs such as those from the RISC complex and miRNA machinery. These miRNAs associated with BRAF oncogene as listed in Table 1 may contribute to tumour formation by regulating genes involved in the MAPK pathway in which BRAF plays a role.
This data suggests the potential utility of miRNA expression patterns as a molecular screen for PTC to use as a diagnostic adjuvant to already established morphological criteria. A recent study by this group examined the effect the ret/PTC 1 oncogene has on miRNA expression in PTC cell lines. ret/PTC 1 rearrangements are generally associated with classic papillary architecture and indolent behaviour in thyroid carcinoma. This study showed a distinct pattern of miRNAs associated with the ret/PTC 1 oncogene when compared to a normal thyroid cell lines .
It is noteworthy that different miRNA expression profiles have been shown contingent upon ret/PTC 1 activated and BRAF mutated cell lines. This suggests that miRNA expression analysis may also be of prognostic significance. Further miRNA analysis in a larger cohort of samples will help validate and consolidate the results discussed above. Table 4 illustrates a compilation of miRNAs that are differentially expressed in PTC cell lines (either ret/PTC-1 or BRAF mutant), suggesting an important role for these miRNAs in thyroid neoplasia progression independent of molecular trigger. They may be involved in regulating the central MAPK pathway involved in thyroid tumourigenesis in which RET and BRAF oncogenes play a part.
Although the exact function of many miRNAs remains to be elucidated it is clear that miRNAs potentially have a broad influence over several diverse genetic pathways and that their deregulation is likely to contribute to disease including cancer. An exciting future prospect is that the miRNA patterns associated with a particular tumour may ultimately be of diagnostic and prognostic significance and contribute to the understanding of the molecular pathogenesis and gene regulatory processes of cancer.
Materials and methods
Cell culture and transfection
Nthy-ori 3-1 (ECACC, Wiltshire, UK) is a thyroid follicular epithelial cell line derived from normal thyroid tissue of an adult that has been transfected with a plasmid encoding for the SV40 large T gene. KAT10 is derived from a papillary thyroid carcinoma cell line with heterozygous BRAF V600E mutation. Both cell lines were grown to confluence in a humidified atmosphere containing 5% CO2 at 37°C in the following plating medium: RPMI 1640 with 2 mM L-glutamine, 10% Foetal calf serum (FCS), Penicillin (100 U/ml) and Streptomycin (100 μg/ml). A BRAF V600E expressing plasmid, pMCEF-V600EB-RAF, was transfected into Nthy-ori 3-1. Transfections were carried out using Genejuice™ transfection agent (Novagen, Germany) using the recommended protocol and transfected cells were grown in the presence of Gentamicin (Sigma Aldrich) antibiotic to yield pure cultures.
Nucleic acid extraction
Following trypsinisation of cultured cells, cell pellets (3–4 × 106 cells) were collected and total RNA was extracted using RNeasy® mini kit (Qiagen Ltd., West Sussex, UK). RNA quantity and quality are assessed using a NanoDrop® ND-1000 Spectrophotometer (Wilmington, USA) and RNA 6000 Nano LabChip® Kit in conjunction with the Agilent 2100 Bioanalyser (Agilent technologies, Waldbronn, Germany) as illustrated in figure 4. Taqman® SNP detection was used for BRAF V600E mutation detection as previously described by Smyth et al . Primers and probes used in this experiment were designed and used to the manufacturer's recommendations. The primers/probes used were as follows: 5' CAT GAA GAC CTC ACA GTA AAA ATA GGT GAT 3' [BRAF-F], 5' GGA TCC AGA CAA CTG TTC AAA CTG A 3' [BRAF-R], VIC-5' CCA TCG AGA TTT CAC TGT AG 3' [BRAF-PWT], and FAM-5' CCA TCG AGA TTT CTC TGT AG 3' [BRAF-PMUT]. Amplification and analysis was performed using an ABI Prism 7000 Sequence Detection System (Applied Biosystems, CA, USA) for 40 cycles (92°C for 15 sec, 60°C for 1 min).
Western blot analysis was carried out on the transfected cell lines to confirm BRAF protein expression.
The Applied Biosystems 1700 Expression Array System is based on a microarray design that represents the whole human genome. The V.2 array has 32,878 60-mer oligonucleotide probes for the interrogation of 29,098 individual human gene and more than 1,000 control probes. V2 arrays were used to analyse the transcriptional profiles of the cell line RNA samples in this study. Digoxigenin-UTP labeled cRNA was generated and linearly amplified from 2 μg of total RNA using Applied Biosystems Chemiluminescent RT-IVT Labelling Kit v 2.0 following the manufacturer's protocol. Array hybridization, chemiluminescence detection, image acquisition and analysis were performed using Applied Biosystems Chemiluminescence Detection Kit and 1700 Chemiluminescent Mircoarray Analyzer following manufacturer's guidelines. Each microarray was initially pre-hybridised at 55°C for 1 hr in hybridization buffer with blocking reagent. 10 μg of labeled cRNA targets were fragemented by incubating with fragmentation buffer at 60°C for 30 min, mixed with internal control target (ICT, 24-mer oligo labeled with LIZ fluorescent dye) and hybridized to each pre-hybed microarray in a 1.5 ml volume at 55°C for 16 hr. After hybridization, the arrays were washed with hybridization wash buffer and chemiluminscence rinse buffer. Enhanced chemiluminescenct signals were generated by incubating arrays with anti-digoxigenin alkaline phosphatase, enhanced with Chemiluminescence Enhancing Solution and finally by adding Chemiluminescence Substrate. Images were collected for each microarray using the 1700 analyser. Images were autogridded and the chemiluminescent signals were quantified, corrected for background and spot and spatially normalized. Replicates were performed. Microarrays were analysed using Spotfire Decision Site™ for functional genomics (Spotfire AB, Göteborg, Sweden) and R version 1.9.1 [a free language and environment for statistical analysis and graphics] (R Development Core Team, 2004).
Applied Biosystems TaqMan® microRNA (miRNA) assays are designed to detect and quantify mature miRNAs using a looped-primer real time PCR. The human early access panel used in this study contained 160 individual assays for identified human miRNAs. The assay involved 2 steps: Step one; a Stem-looped RT, and Step two; a Real Time PCR. Briefly, single stranded cDNA was generated from total RNA sample by reverse transcription using the Applied Biosystems High-Capacity cDNA Archive Kit (Applied Biosystems, CA, USA) following manufacturer's protocol. RT reactions contained 10 ng of total RNA, 50 nM stem-looped RT primer, 1 × RT buffer, 0.25 mM each of dNTPs, 3.33 U/μl Multiscribe reverse transcriptase and 0.25 U/μl RNase Inhibitor. PCR amplification was carried out using sequence specific primers on the Applied Biosystems 7900 HT Fast Real-Time PCR system. The reactions were incubated in a 96-well optical plate at 95°C for 10 min, following by 40 cycles of 95°C for 15 s and 60°C for 10 min. Analysis of relative miRNA expression data was performed using ΔΔCT method with hsa-let-7a as an endogenous control . Two negative controls were also used, ath-mir159a and cel-lin-4 as they show no expression in human tissue.
Kondo T, Ezzat S, Asa SL: Pathogenetic mechanisms in thyroid follicular-cell neoplasia. Nat Rev Cancer. 2006, 6 (4): 292-306. 10.1038/nrc1836
Sheils O: Molecular classification and biomarker discovery in papillary thyroid carcinoma. Expert Rev. Mol. Diagn. 2005, 5 (6): 927-946. 10.1586/14737220.127.116.117
Nikiforova MN, Kimura ET, Gandhi M: BRAF mutations in thyroid tumors are restricted to papillary carcinomas and anaplastic or poorly differentiated carcinomas arising from papillary carcinomas. J Clin Endocrinol Metab. 2003, 88: 5399-5404. 10.1210/jc.2003-030838
Xing M: BRAF mutation in thyroid cancer. Endocrin Rel Cancer. 2005, 12: 245-262. 10.1677/erc.1.0978.
Finley DJ, Arora N, Zhu B: Molecular profiling distinguishes papillary carcinoma from benign thyroid nodules. J Clin Endocrinol Metab. 2004, 89: 3214-3223. 10.1210/jc.2003-031811
Huang Y, Prasad M, Lemon WJ: Gene expression in papillary thyroid carcinoma reveals highly consistent profiles. Proc Natl Acad Sci USA. 2001, 98 (26): 15044-15049. 10.1073/pnas.251547398
Jarzab B, Wiench M, Fujarewicz K, Simek K, Jarzab M, Oczko-Wojciechowskam , Wloch J, Czarniecka A, Cmielik E, Lange D, Daulaczek A, Szpak S, Gubala E: Gene Expression Profile of Papillary Thyroid Cancer: Sources of Variability and Diagnostic Implications. Cancer Res. 2005, 65: 1587-1597. 10.1158/0008-5472.CAN-04-3078
Lee RC, Feinbaum RL, Ambros V: The C-elegans heterochronic gene lin-4 encodes small RNAs with antisense complementarity to lin-14. Cell. 1993, 75: 843-854. 10.1016/0092-8674(93)90529-Y
Feinbaum R, Ambros V: The Timing of lin-4 RNA accumulation controls the timing of postembryonic developmental events in Caenorhabditis elegans. Dev Biol. 1999, 210: 87-95. 10.1006/dbio.1999.9272
Zhang B, Pan X, Cobb GP, Anderson TA: microRNAs as oncogenes and tumour suppressors. Dev Biol. 2007, 302 (1): 1-12. 10.1016/j.ydbio.2006.08.028
Zhang B, Pan X, Wang Q, Cobb GP, Anderson TA: Computational identification of microRNAs and their targets. Comput Biol Chem. 2006, 30 (6): 395-407. 10.1016/j.compbiolchem.2006.08.006
Esquela-Kerscher A, Slack FJ: Oncomirs-MicroRNAs with a role in cancer. Nat Rev Cancer. 2006, 6 (4): 259-269. 10.1038/nrc1840
Calin GA, Sevignani C, Dumitru CD, Hyslop T, Noch E, Yendamuri S, Shimizu M, Rattan S, Bullrich F, Negrini M, Croce CM: Human microRNA genes are frequently located at fragile sites and genomic regions involved in cancers. Proc Natl Acad Sci USA. 2004, 101 (9): 2999-3004. 10.1073/pnas.0307323101
Hwang HW, Mendell JT: MicroRNAs in cell proliferation, cell death and tumourigensis. B J Cancer. 2006, 94: 776-780. 10.1038/sj.bjc.6603023.
Miska EA: How MicroRNAs control cell division, differentiation and death. Curr Opin Genet Dev. 2005, 15: 563-568. 10.1016/j.gde.2005.08.005
Baskerville S, Bartel DP: Microarray profiling of microRNAs reveals frequent coexpression with neighbouring miRNAs and host genes. RNA. 2005, 3: 241-247. 10.1261/rna.7240905.
Calin GA, Liu CG, Sevignani C, Ferracin M, Felli N, Dumitru CD, Shimizu M, Cimmino A, Zupo S, Dono M, Dell'Aquila ML, Alder H, Rassenti L, Kipps TJ, Bullrich F, Negrini M, Croce CM: MicroRNA profiling reveals distinct signature in B cell chronic lymphocytic leukemias. Proc Natl Acad Sci USA. 2004, 11755-11760.
Ciafre SA, Galardi S, Mangiola A, Ferracin M, Liu CG, Sabatino G, Maira G, Croce CM, Farace MG: Extensive modulation of a set of micrornas in primary glioblastoma. Biochem Biophys Res Commun. 2005, 334: 1351-8. 10.1016/j.bbrc.2005.07.030
Lu J, Getz G, Miska EA, Alvarez-Saavedra E, Lamb J, Peck D, Sweet-Cordero A, Ebert BL, Mak RH, Ferrando AA, Downing JR, Jacks T, Horvitz R, Golub TR: MicroRNA expression profiles classify human cancers. Nature. 2005, 435 (7045): 834-838. 10.1038/nature03702
Murakami Y, Yasuda T, Saigo T, Toyoda H, Okanoue T, Shimotohno K: Comprehansive analysis of microRNA expression patterns in heptaceulluar carcinoma and non tumourous tissue. Oncogene. 2006, 25: 2537-2545. 10.1038/sj.onc.1209283
Volinia S, Calin GA, Liu CG, Ambs S, Cimmino A, Petrocca F, Visone R, Iorio M, Roldo C, Ferracin M, Prueitt RL, Yanaihara N, Lanza G, Scarpa A, Vecchione A, Negrini M, Harris CC, Croce CM: A microRNA expression signature of human solid tumors defines cancer gene targets. Proc Natl Acad Sci USA. 2006, 103 (7): 2257-61. 10.1073/pnas.0510565103
Roldo C, Missiaglia E, Hagan JP, Falconi M, Capelli P, Bersani S, Calin GA, Volinia S, Liu CG, Scarpa A, Croce CM: MicroRNA expression abnormalities in pancreatic endocrine and acinar tumors are associated with distinctive pathologic features and clinical behavior. J Clin Oncol. 2006, 4677-84. 29,
Liu CG, Calin GA, Meloon B, Gamliel N, Sevignani C, Ferracin M, Dumitru C, Shimizu M, Zupo S, Dono M, Alder H, Bullrich F, Negrini M, Croce CM: An oligonucleotide microchip for genome wide microRNA profiling in human and mouse tissues. Proc Natl Acad Sci USA. 2004, 101: 9740-9744. 10.1073/pnas.0403293101
Iorio MV, Ferracin M, Liu CG, Veronese A, Spizzo R, Sabbioni S, Magri E, Pedriali M, Fabbri M, Campiglio M, Menard S, Palazzo JP, Rosenberg A, Musiani P, Volinia S, Nenci I, Calin GA, Querzoli P, Negrini M, Croce CM: MicroRNA gene expression deregulation in human breast cancer. Cancer Res. 2005, 65: 7065-7070. 10.1158/0008-5472.CAN-05-1783
Michael MZ, O'Connor Sm, Holst Pellekaan NG, Young GP, James RJ: Reduced accumulation of specific microRNAs in colorectal neoplasia. Mol Cancer Res. 2003, 1 (12): 882-91.
Krek A, Grun D, Poy M, Wolf R, Rosenburg L, Epstein EJ, MacMenamin P, da Piedade I, Grunsalus K, Stoffel M, Rajewsky N: Combinatorial MicroRNA Target Predictions. Nat Genet. 2005, 37: 495-500. 10.1038/ng1536
Lewis BP, Shih I, Jones-Rhoades MW, Bartel DP, Burge CB: Prediction of Mammalian MicroRNA Targets. Cell. 2005, 115: 787-798. 10.1016/S0092-8674(03)01018-3.
Griffiths-Jones S, Grocock RJ, Van Dongen S, Bateman A, Enright AJ: miRBase: microRNA sequences, targets and gene nomenclature. Nucleic Acid Res. 2006, 34: 140-144. 10.1093/nar/gkj112
Bandres E, Cubedo E, Agirre X, Malumbres R, Zarate R, Ramirez N, Abajo A, Navarro A, Moreno I, Monzo M, Garcia-Foncillas J: Identification by real time PCR of 13 mature mirnas differentially expressed in colorectal cancer and nontumoural tissues. Mol Cancer. 2006, 5: 29- 10.1186/1476-4598-5-29
Cimmino A, Calin GA, Fabbri M, Lorio MV, Ferracin M, Shimizu M, Wojcik SE, Aqeilan RI, Zupo S, Dono M, Rassenti L, alder H, Volinia S, Liu CG, Kipps TJ, Negrini M, Croce CM: mir-15 and miR-16 induce apoptosis by targeting BCL2. Proc Natl Acad Sci USA. 2005, 102: 13944-13949. 10.1073/pnas.0506654102
Cummins JM, Yiping H, Pagliarini R, Diaz L, Sjobiom T, Barad O, Bentawich Z, Szafranaska AE, labourier E, Raymond CK, Roberts BS, Juhl H, Kinzier KW, Vogelstein B, Velculescu VE: The colorectal microRNAome. Proc Natl Acad Sci USA. 2006, 103 (10): 3687-3692. 10.1073/pnas.0511155103
Johnson ST, Grosshans H, Shingara J, Byrom M, Jarvis R, Cheng A, Labouier E, Reinert KL, Brown D, Slack FJ: Ras is regulated by the let-7 microRNA family. Cell. 2005, 120: 635-647. 10.1016/j.cell.2005.01.014
Saito Y, Liang G, Egger G, Friedman JM, Chuang JC, Coetzee GA, Jones PA: Specific activation of microRNA-127 with downregulation of the protooncogene BCL6 by chromatin-modifying drugs in human cancer cells. Cancer Cell. 2006, 9 (6): 435-43. 10.1016/j.ccr.2006.04.020
Feili N, Fontana L, Pelosi E, Botta R, Bonci D, Facchiano F, Liuzzi F, Lulli V, Morsilli O, Santoro S, Valtieri M, Calin GA, Liu CG, Sorrentino A, Croce CM, Peschie C: MicroRNAs 221 and 222 inhibit normal erythropoiesis and erythroleukemic cell growth via kit receptor down-modulation. Proc Natl Acad Sci USA. 2005, 102: 18081-18066. 10.1073/pnas.0506216102
Yanaihara N, Caplen N, Bowman E, Seike M, Kumamoto K, Yi M, Stephens RM, Okamoto A, Yokota Y, Tanaka T, Calin GA, Liu CG, Croce CM, Harris CC: Unique microRNA molecular profiles in lung cancer diagnosis and prognosis. Cancer Cell. 2006, 9;: 189-198. 10.1016/j.ccr.2006.01.025
Huiling H, Jazdzewski K, Li W, Liyanarachchi S, Nagy R, Volinia S, Calin GA, Liu CG, Franssila K, Suster S, Kloos RT, Croce CM, De la Chapelle A: The role of microRNA genes in papillary thyroid carcinoma. Proc Natl Acad Sci USA. 102: 19075-19080.
Pallante P, Viscone R, Ferracin M, Ferraro A, Berlingieri MT, Troncone G, Chiappette G, Liu CG, Santoro M, Negrini M, Croce CM, Fusco A: MicroRNA deregulation in human papillary thyroid carcinomas. Endocr Relat Cancer. 2006, 13: 497-508. 10.1677/erc.1.01209
Mitsutake N, Knauf JA, Mitsutake Shin, Mesa C, Zhang L: Conditional BRAF V600E expression induces DNA synthesis, apoptosis, dedifferentiation and chromosomal instability in thyroid PCCL3 cells. Cancer Res. 2005, 65: 2465-2473. 10.1158/0008-5472.CAN-04-3314
Cahill S, Smyth P, Finn SP, Denning K, Flavin R, O'Regan EM, Li JH, Potratz A, Guenther SM, Henfrey R, O'Leary JJ, Sheils O: Effect of ret/PTC 1 rearrangement on transcription and post-transcriptional regulation in a papillary thyroid carcinoma model. Mol Cancer. 2006, 5: 70- 10.1186/1476-4598-5-70
Smyth P, Finn S, Cahill S, O'Regan E, Flavin R, O'Leary JJ, Sheils O: ret/PTC and BRAF act as distinct molecular, time dependent triggers in a sporadic Irish cohort of papillary thyroid carcinoma. Int J Surg Pathol. 2005, 13: 1-8.
Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)). Methods. 2001, 25 (4): 402-8. 10.1006/meth.2001.1262