Effect of ret/PTC 1 rearrangement on transcription and post-transcriptional regulation in a papillary thyroid carcinoma model

Background microRNAs (miRNAs) are a group of non-coding single stranded RNAs measuring approximately 22 nt in length that have been found to control cell growth, differentiation and apoptosis. miRNAs negatively regulate their target genes and recently have been implicated in tumourigenesis. Furthermore, miRNA expression profiling correlates with various cancers, with these genes thought to act as both tumour suppressors and oncogenes. ret/PTC 1 is an oncogene with constitutive kinase activity implicated in the development of papillary thyroid carcinoma (PTC). This rearrangement leads to aberrant MAPK activation that is implicated in PTC tumourigenesis. Aim The aim of this study was to identify the effect that ret/PTC 1 has on transcription and post-transcriptional regulation in PTC by using DNA microarray and microRNA analysis. Results DNA microarray analysis revealed a group of genes differentially expressed between normal thyroid cell lines and those harbouring a ret/PTC 1 rearrangement. Furthermore, a unique miRNA expression signature differentiated between PTC cell lines with ret/PTC 1 and a normal thyroid cell line. 21 miRNAs showed significant overexpression and 14 miRNAs showed underexpression in these cell lines when compared to normal thyroid. Several of these up/down regulated miRNAs may be involved in PTC pathogenesis.


Background
In 1985, RET was first identified as a novel oncogene. RET encodes a receptor tyrosine kinase, the ligands of which belong to the glial cell line-derived neurotropic factor (GDNF) family [1]. Chimeric RET oncogenes, designated ret/PTC, display a ligand independent constitutive tyro-sine kinase activity and are implicated in the development of papillary thyroid carcinoma (PTC). They are formed from the juxtaposition of the genomic region coding for the tyrosine kinase domain of RET with the 5'-promoter regions of a variety of unrelated genes. At least 15 chimeric mRNAs involving ten distinct donor genes have been described. Of these, the most commonly occurring in PTC are ret/PTC 1 [H4-(CCDC6)-RET] [2] and ret/PTC 3 [ELE1-RET] [3].
Papillary thyroid carcinoma is the most frequently occurring type of thyroid malignancy accounting for approximately 80%-90% of cases. A unique feature of papillary thyroid cancers is that oncogenes such as ret/PTC, activate effectors that signal along the mitogen activated protein kinase (MAPK) pathway playing important roles in the pathogenesis of thyroid tumourogenesis. Genetic studies suggest that ret/PTC activation is one of the key first steps in thyroid cancer pathogenesis. The reported prevalence of ret/PTC rearrangements in papillary carcinoma varies from 3% to as high as 85% with the highest incidence in radiation associated PTC. This variation reflects the different sensitivities of the techniques used as well as geographic variations and the influence of environmental factors such as ionizing radiation exposure [4,5]. Several reports have shown that ret/PTC rearrangements are in general associated with classic papillary architecture in thyroid carcinoma and RET has been suggested to cause structural and nuclear peculiarities of PTC. PTCs expressing ret/PTC tend to be small and do not exhibit a tendency to progress to an aggressive phenotype [6][7][8]. However, the exact role of ret/PTC in the context of papillary thyroid carcinoma remains unclear and the complete repertoire of genes and signalling pathways involved in pathogenesis of thyroid disease remains poorly defined. A further shortcoming remains in the inadequacy of diagnostic and prognostic markers with frequent reliance on conventional morphological assessment.
Recently it has been proposed that microRNAs (miRNAs), small non coding RNAs, may play a role in tumour formation including PTC development. Ambrose et al described Lin-4 over a decade ago and there has been an explosion of interest in the area since [9,10]. microRNAs are cleaved from hairpin shaped precursors to yield mature miRNAs of approximately 22 nucleotides in length. miRNAs negatively regulate gene expression and function by targeting mRNAs for cleavage or translational repression depending on the degree of complementarity between the miRNA and the target mRNA [11,12]. Because of their relatively recent discovery, little is known about the exact processes in which miRNAs function or the exact nature of their regulation. Studies on the regulatory roles of miRNAs in a variety of organisms suggest that they have critical roles in central biological processes and are regulators of in devel-opment, cell proliferation, cell differentiation, stress resistance, metabolism and apoptosis, all of which are involved in tumourigenesis [13,14]. Not surprisingly as a consequence aberrant regulation of miRNA expression has been implicated in human disease including cancer. Furthermore, a recent study showed that greater than 50% of annotated human miRNAs are located in areas of the genome known as fragile sites that are associated with cancer [15]. This further supports a role for miRNAs in cancer progression. Differential expression of miRNAs has been shown between malignant and normal tissues of the same tissue type and between different types of tumour. This has lead to the hypothesis that tumours may each have a discrete miRNA signature [16][17][18]. It has also been reported that numerous oncogenes and tumour suppressor genes are potentially regulated by miRNAs, or indeed miRNAs may act as oncogenes or tumour suppressor genes themselves. miRNA species may therefore be potentially used as biomarkers for tumour diagnostics and treatment.
Several groups have used a variety of techniques in order to quantify miRNAs the most common including oligonucleotide microarray technology, bead based flow cytometric analysis, a modified invader assay and a new single molecule technique. Each of these methods has distinct advantages and disadvantages. Real-time PCR however, has unparalleled sensitivity and specificity. It is also more quantative than DNA microarray technology which is the preferred method in recent publications. The TaqMan ® microRNA RT-PCR method used in this study is a new real time quantification method to accurately detect mature miRNAs. The stem-loop structure of the primers is specific to the 3' end of the mature miRNA and creates a steric hindrance to prevent priming of the precursor miRNA. The aim of this study was to identify the effect the ret/PTC 1 oncogene has on transcription and post-transcriptional regulation in tumourigenesis by using a combination of DNA microarray and microRNA analysis.

Differential gene expression using DNA Microarrays Upregulated genes
In order to investigate differential gene expression in papillary thyroid carcinoma cell lines, a comparison was made between the two cell lines harbouring ret/PTC 1 rearrangement [Nthy-ori transfected with ret/PTC 1 and TPC-1] and normal thyroid cell lines. To identify potential markers of malignancy t-test with FDR correction (<0.1) were used to compare the two groups [19]. A p-value cut off of < 0.01 and a fold change difference of >2 resulted in a list of 47 genes (Table 1).

Downregulated genes
Similar cut off criteria were applied in analysing downregulated mRNA species. 29 genes were found to be statistically significant, see Table 2. Hierarchical clustering was performed based on statistically different genes (Fig. 1).

Differential miRNA expression using TaqMan ® microRNA assays
Upregulated miRNAs were considered to be significantly upregulated when there was a fold change of > 2 fold in PTC cell lines  Fig 3)

mRNA (transcription) and miRNA (regulation) expression correlation
miRNAs interact and regulate their target genes by cleavage or translational inhibition.
The current challenge is to identify biologically relevant targets that are regulated by individual miRNAs. Further complicating this task is the fact that every miRNA can potentially bind and regulate many mRNA targets and each mRNA can be bound and regulated by several miR-NAs. miRNAs usually bind to their targets with incomplete complimentarity within the 3' untranslated region (UTR) of the mRNA target. Because of this the identification of targets is impossible with a simple BLAST search. Most bioinformatic methods use the first 2-8 bases of the mature miRNA sequence to search for complementarity to the 3' UTR of all expressed genes. This will yield a list of genes with diverse functions that a single miRNA can bind to. There are three publicly available databases for prediction of miRNA targets. These are miRBase [20], PICTAR [21], and TARGETSCAN [22]. For example, all animal miRNA sequences from the miRBase Sequence database are scanned against 3'-UTRs predicted from all available species in Ensembl along with Caenorhabditis briggsae and Drosophila pseudoobscura. P-values are assigned to individual miRNA-target binding sites, multiple sites in a single UTR, and sites that appear to be conserved in multiple species. The interface connects each miRNA to a list of pre-  Table 3. This is illustrated in Table 5. Thus, there was correlation between the data sets obtained from the functional expression (microarray) and the regulatory (miRNA) data sets.

Discussion
The primary aim of this study was to try to elucidate the effects of the RET proto-oncogene on gene expression and transcriptional regulation in a matrix of cells lines harbouring ret/PTC 1 rearrangement. This was performed using Applied Biosystems whole genome expression array technology and TaqMan ® microRNA assays. The sequencing of the human genome and advances in microarray technology has significantly increased understanding of the molecular pathogenesis of thyroid tumour formation in recent years. On order to gain a better understanding of the genes and pathways involved in PTC and unveil potential biomarkers, several groups have applied DNA microarray technology to show the genetic profiles in PTC [23][24][25][26]. These profiles uncovered genes that were sug-gested to be involved in processes such cell adhesion and cell structure, inhibition and induction of apoptosis, cell proliferation and differentiation.
This study showed a group of differentially expressed genes in two cell lines contingent on a rearrangement of the RET oncogene when compared to normal thyroid. Among the genes found to be over-expressed, those involved in biological processes such as cell differentiation and proliferation were over-represented. Such genes included CEBPB, CCNG1, IFITM3 and HTRA1. CEBPB is a bZip transcription factor that cooperates with Stat 6 to activate transcription [27]. It has been shown to have a critical role in Ras mediated tumourigenesis and cell survival and has been implicated as a target for tumour inhibition. Studies have shown that it is expressed in the cytoplasm of PTC and is implicated to play a role in the regulation of thyroid specific genes and processes and has altered function in thyroid carcinoma [28]. It is involved in the positive regulation of genes involved in the immune and inflammatory response, a group which has been shown to be associated with ret/PTC 1 and tumour formation in this study and in previous thyroid literature [29]. CEBPB also has been shown to be a regulator of Cox-2; Cox-2 has recently been implicated as a potential biomarker in thyroid neoplasm [30].
HTRA1 [31] is a member of the trypsin family of serine proteases and is a regulator of cell growth. It has been implicated in the pathology of several diseases. Several papers have implicated matrix metalloproteases in PTC tumour progression; during tumour invasion they degrade protein components of the extracellular matrix aiding angiogenesis [32]. HTRA has been proven to cause the induction of matrix metalloproteases 1 and 3.
Cell signalling also showed an increase in ret/PTC 1 cell lines, SEMA3F is a protein involved in signalling. Endogenous SEMA3F causes cells to round and loose extracellular contacts a process which may contribute to the classic appearance of PTC. Membralin (C19orf6) is a novel gene originally cloned from a human ovarian cancer cell line. It is suggested to be a tumour associated marker in ovarian and colorectal cancer [33,34]. It was also shown to be over-expressed in ret/PTC 1 cell lines when compared to normal. It is interesting to hypothesize that there may be a link between membralin overexpression and the papillary architecture seen in PTC and in ovarian cancers.
Furthermore, genes involved in p53 pathway were overrepresented. P300/CBP-associated factor (PCAF) was one such gene and is a co-activator of the tumour suppressor p53. PCAF was found to be induced by p53 in breast can-cer cell lines and indirectly regulated ERK1 effecting cellular proliferation [35].
Among the genes shown to be down-regulated in ret/PTC 1 cell lines when compared to normal, those involved in cell structure and motility were over-represented. Examples of such genes include DCTN5, TPM1, TPM3, CRP1 and Keratin type 1. Tumour cells are held together by direct cell -cell contact and by adhesion to the extracellular matrix. Loss of adhesion is thought to promote tumour invasiveness and increase the metastatic potential of carcinoma. Loss of these cell structure genes might also be involved in the papillary pattern of PTC. Tropomysoin (TM) has a role in regulating cellular functions associated with cytoskeletal remodelling. Suppression of TMs is a prominent feature of many transformed cells. TM1 expression is abolished in human breast cancer cell lines and also in the ret/PTC 1 cell lines in this study. This downregulation may destabilize microfilament architecture and facilitate survival of neoplastic cells and promote malignant growth [36]. PSMD2 (TRAP2) participates in TNF signalling pathway and was found to be downregulated. This suggested a role for PSMD2 in cell survival. Rab32 inactivation may represent a component of the oncogenic pathway of microsatellite-unstable gastrointestinal adenocarcinoma. It encodes an A-kinase anchoring protein and is often hypermethylated in gastric and colon cancers and was found to also be downregulated in PTC cell lines. Rab32 has also been shown to induce apoptosis [37,38].
Of particular interest in this study was the discovery of the underexpression of DROSHA in PTC cell lines. This is the core nuclease that executes the initiation step of miRNA processing in the nucleus. DROSHA collaborates with Dicer in stepwise processing of miRNAs and has a key role in miRNA mediated gene regulation processes such as development and differentiation [39]. Components of the miRNA-machinery have been implicated in tumourigenesis and reduced levels have been shown to be associated with cancer [40]. Its differential expression may effect the expression of certain miRNAs such as those in table 3 and consequently play a role in gene expression regulation.
This data corroborates the importance of these pathways in PTC. The genes involved may play a role in tumour initiation and in the morphological characteristic of classic PTC due to the ret/PTC 1 oncogene. Several reports have shown that ret/PTC 1 rearrangements are in general associated with classic papillary architecture in thyroid carcinoma. The profiles examined contained genes that were suggested to be involved in processes such as cell adhesion and cell structure, such as tropomysoins, and inhibition and induction of apoptosis. This aberrant expression has an impact on the molecular structure and the extracellular matrix which may facilitate destabilization of the microfilament architecture. While also adding to tumour spread this phenomenon may contribute to the classic papillary appearance associated with this rearrangement. Gene expression has recently been shown to be regulated by small noncoding RNAs called microRNAs (miRNAs). It is believed that miRNAs are responsible for fine tuning gene expression and are thought to regulate approximately 30% of the human genome, making them one of the largest classes of gene regulator. Several clues have pointed to the potential role of miRNAs in cancer and suggest that aberrations in miRNA expression may be important in tumour progression by regulating genes and pathways involved in cancer associated processes. Studies have reported that reduced levels of the Let-7 family of miRNAs have an association with lung cancer tumourigenesis [41,42]. The miRNA cluster of miR-17-miR-92 is highly expressed in many B-cell chronic lymphomas [43]. Furthermore, Mir 15a-Mir 16 expression is upregulated in B-cell chronic lymphoma and leukaemia and falls within a region thought to harbour a tumour suppressor gene. Glioblastomas and breast tumours have also been reported to show over expression of miR-21 [44,45]. This study examined, by real-time PCR, the expression of miRNA in papillary thyroid carcinoma cell lines with ret/ PTC1 rearrangement. Tables 3 shows a set of differentially expressed miRNAs in 2 cell lines with ret/PTC1 rearrangement when compared to normal thyroid cell line. 21 miRNAs were found to be overexpressed and 14 miRNAs were underexpressed in both cell lines when compared to Nthy-ori 3-1. These miRNAs are effected by direct ret/PTC 1 expression or by the expression of other genes that are induced by ret/PTC 1 such as transcription factors. Little is known about the function of these miRNAs and only a subset has been found to be differentially expressed in other cancers. Of the miRNAs in table 3, miR-181a expression was downregulated and expression of miR-200b and miR-141 was up-regulated. Ciafre et al found miR-181a expression to be downregulated in glioblastoma [45]. Studies have reported that miR-181 upregulated during differentiation and suggest that miR-181 downregulates the homeobox protein Hox-A11 (a repressor of the differentiation process), thus establishing a functional link between miR-181 and the complex process of differentiation. In fully differ-   LogNthy LogNthyRet

Fold Change of miRNAs in ret/PTC harbouring cell lines Vs Normal
LogTPC-1 entiated tissue miR-181 was found to be downregulated, PTCs are well-differentiated carcinomas therefore this may explain downregulation of miR-181a in this study.
miR-200b and miR-141 have also been shown to be highly overexpressed in malignant cholangiocytes and in colon carcinoma [46,47]. Studies showed that their inhi-bition increased sensitivity to the anticancer drug gemcitabine and also decreased cell growth suggesting a role for their involvement in tumour formation.
An immediate obstacle facing research groups at present is the lack of confirmed miRNA gene targets. Although multiple approaches have been suggested from bioinformat- Heat map of miRNA expression Figure 3 Heat map of miRNA expression. miRNA expression. Each miRNA listed was detected as significantly differentially expressed between ret/PTC 1 cell lines and N-thy-ori. The delta CT values for each miRNA were used to create the heat map. ics including many prediction algorithms accurate target prediction and validation are still major obstacles facing miRNA researchers. Each miRNA found to be significantly up/down regulated from table 3 was imputed into the algorithms, miRBase, PICTAR and TARGETSCAN to search for potential miRNA:miRNA target pairings. The potential binding partners of 6 of these miRNAs are listed in table 4. It is interesting to note that these differentially expressed miRNAs potentially bind and regulate a number of genes that have implications in thyroid carcinoma progression such as genes involved in thyroid function, MAPKKK cascade, retinoic acid receptors, cell adhesion and cell structure associated genes, genes involved in cell signalling, G-protein mediated signalling, oncogenesis and cell cycle control. Examples of these genes include thyroid hormone interacting protein, protocadherin, MAP3K10, KIT, PDGFRA and exportin 4. Therefore, it seems plausible that this subset of miRNAs may be important in the progression of thyroid carcinoma by regulating genes involved in thyroid neoplasia such as those found in this study discussed above.
Analysis of DNA microarray data showed a group of genes that are differentially expressed in ret/PTC1 harbouring cell lines (Table 1 and 2). Comparing this list of genes with the list of putative miRNA targets yielded from analysis using PICTAR, TARGETSCAN and miRBase showed that 15 of the differentially expressed genes are targeted by several miRNAs that are significantly up/downregulated in cell lines with ret/PTC 1 rearrangement. This data is shown in Table 5. Again it is interesting to speculate that these genes are regulated by their corresponding miRNA thereby contributing to tumour development in thyroid. It remains difficult to identify the precise target mRNA of each miRNA species. Therefore without further experimental investigation any mRNA:miRNA partners are essentially speculative and selection involves a degree of subjectivity.
Recent studies on miRNA deregulation in PTC found an aberrant miRNA expression profile in PTCs compared to normal thyroid tissues. In particular, a significant increase in miR-222, miR-221, miR-146 and miR-181b was seen [48,49]. In this study miR-146 was upregulated 500 fold in TPC-1 cell line only. It is possible this miRNA may play a role in a late stages of tumour formation and progression; given the transfected cell line had only been exposed to ret/PTC1 for a short period (3 passages). Alternatively, the difference may reflect the accumulation of several genetic insults. The other miRNAs [miR-222, miR-221, miR-181b], however, were not found to be upregulated significantly. This may be due to the differences in experimental approaches; de la Chapelle et al and Fusco A et al used miRNA microarrays. They also used different chemistry and fresh tissue samples; also the RET status of each PTC case was not disclosed.
This field of miRNA research is young and many important questions still remain. It is clear that miRNAs exert a regulatory role on protein coding genes in a variety of species and their deregulation may contribute to the development of human disease. However, the precise biological function of miRNAs has yet to be fully elucidated.
Even without more in-depth knowledge of their exact function or role in PTC, it seems reasonable that a distinct signature of miRNAs or miRNA expression profiles assigned for thyroid tumours might be used as an adjunct in diagnosis and prognostication of thyroid neoplasia. It has been suggested that miRNA signatures may be more effective than profiles of protein coding genes in distinguishing tumours from benign lesions. Fortunately miR- NAs, unlike mRNA targets remain intact in routinely collected paraffin embedded samples, a fact which increases their potential use as clinical biomarkers. miR-NAs also have the potential as possible treatment targets. One example of potential miRNA targeted therapy has been shown with miR-127 in bladder cancer. This suggests the possibility of direct targeting of miRNAs that are amplified or upregulated in patient tumours.

Conclusion
This study is an in vitro model for the biological and regulatory processes that occur in human thyroid diseases due to ret/PTC1 rearrangement. As miRNAs are stable, abundant and easily detectable they represent ideal candidates for effective diagnostic biomarkers. Future exploration on a larger cohort of samples will hopefully determine the correlation between these miRNAs and their host genes in PTC and help aid in the identification of miRNA biomarkers.
The cell line is derived from normal thyroid tissue of an adult that has been transfected with a plasmid encoding for the SV40 large T gene.
TPC-1 is a cell line derived from a papillary thyroid carcinoma. This cell line expresses the ret/PTC-1 oncogene.
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). ret/PTC 1 transcript was cloned into the multiple cloning site of pcDNA4/TO using standard molecular techniques involving restriction and ligation. This is part of the TREx tetracycline-regulated mammalian expression system (Invitrogen, UK). Expression of the gene of interest from pcDNA4/TO is controlled by the CMV promoter (into which 2 copies of the tet operator 2 (TetO2) sequence have been inserted in tandem.) Each TetO2 sequence serves as a binding site for 2 molecules of the Tet repressor. The second major component of the T-REx system is the pcDNA6/TR © regulatory vector which expresses high levels of the TetR gene under the control of the human CMV promoter. In the absence of tetracycline, the Tet repressor forms a homodimer that binds with high affinity to each TetO2 sequence in the promoter of the inducible expression vector. Upon addition, tetracycline binds with high affinity to each Tet repressor homodimer and causes a conformational change in the repressor that renders it unable to bind to the Tet operator. The Tet repressor:tetracycline complex then dissociates from the Tet operator and allows induction of transcription from the gene of interest. pcDNA6/ TR and pcDNA4/TO-ret/ptc1 were tranfected into N-thyori cells using Genejuice™ transfection agent (Novagen, Germany) using the recommended protocol and grown in the presence selection agents, Blasticidin (Sigma Aldrich) and Zeocin (Sigma Aldrich) to yield pure cultures.

Microarray analysis
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 genes 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.

Statistical analysis
Microarrays were analysed using R version 1.9.1 [a free language and environment for statistical analysis and graphics] (R Development Core Team, 2004). Arrays were normalised using Quantile normalisation. The primary means of identifying genes that were differentially expressed between groups was based on fold change and statistical significance. An ANOVA test was used to generate p-values for statistical difference between groups. Pvalues were then adjusted for multiple comparisons using a technique described by Benjamini and Hochberg with a cut of FDR <0.1 [19]. Genes were deemed statistically different between groups fold change >2 fold. Hierarchical clustering was performed based on statistically different genes (Fig. 1). Gene ontology analysis was performed using an online database known as PANTHER (panther.appliedbiosystems.com). Its functions include the ability to merge and differentiate genes within lists and the ability to check for over-representation of genes with a particular biological function, or involved in a specific biological pathway.

miRNA analysis
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 covering many of the 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/Ol Multiscribe reverse transcriptase and 0.25 U/Ol 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 15s 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.