Efficient RNAi-induced Protein Knockdown in Somatic Cells Using Diced or Chemically Produced Small Interfering RNAs (siRNA)

Background: RNA interference (RNAi) is a post-transcriptional gene silencing process in which double-stranded RNA (dsRNA) directs the degradation of a specific corresponding target mRNA. The mediators of this process are small dsRNAs of approximately 21 to 23 bp in length, called small interfering RNAs (siRNAs), which can be prepared in vitro and used to direct the degradation of specific mRNAs inside cells. Hence, siRNAs represent a powerful tool to study and control gene and cell function. Rapid progress has been made in the use of siRNA as a means to attenuate the expression of any protein for which the cDNA sequence is known. Individual siRNAs can be chemically synthesized, in vitro-transcribed, or expressed in cells from siRNA expression vectors. However, screening for the most efficient siRNAs for post-transcriptional gene silencing in cells in culture is a laborious and expensive process. In this study, the effectiveness of two siRNA production strategies for the attenuation of abundant proteins for DNA repair were compared in human cells: (a) the in vitro production of siRNA mixtures by the Dicer enzyme (Diced siRNAs); and (b) the chemical synthesis of very specific and unique siRNA sequences (Stealth RNai (TM)). Materials, Methods & Results: For in vitro-produced siRNAs, two segments of the human Ku70 (167 bp in exon 5; and 249 bp in exon 13; NM001469) and Xrcc4 (172 bp in exon 2; and 108 bp in exon 6; NM003401) genes were chosen to generate dsRNA for subsequent "Dicing" to create mixtures of siRNAs. The Diced fragments of siRNA for each gene sequence were pooled and stored at -80 degrees C. Alternatively, chemically synthesized Stealth siRNAs were designed and generated to match two very specific gene sequence regions for each target gene of interest (Ku70 and Xrcc4). HCT116 cells were plated at 30% confluence in 24- or 6-well culture plates. The next day, cells were transfected by lipofection with either Diced or Stealth siRNAs for Ku70 or Xrcc4, in duplicate, at various doses, with blank and sham transfections used as controls. Cells were harvested at 0, 24, 48, 72 and 96 h post-transfection for protein determination. The knockdown of specific targeted gene products was quantified by Western blot using GAPDH as control. Transfection of gene-specific siRNA to either Ku70 or Xrcc4 with both Diced and Stealth siRNAs resulted in a down regulation of the targeted proteins to approximately 10 to 20% of control levels 48 h after transfection, with recovery to pre-treatment levels by 96 h. Discussion: By transfecting cells with Diced or chemically synthesized Stealth siRNAs, Ku70 and Xrcc4, two highly expressed proteins in cells, were effectively attenuated, demonstrating the great potential for the use of both siRNA production strategies as tools to perform loss of function experiments in mammalian cells. In fact, down-regulation of Ku70 and Xrcc4 has been shown to reduce the activity of the non-homologous end joining DNA pathway, a very desirable approach for the use of homologous recombination technology for gene targeting or knockout studies. Stealth RNAi (TM) was developed to achieve high specificity and greater stability when compared with mixtures of enzymatically-produced (Diced) siRNA fragments. In this study, both siRNA approaches inhibited the expression of Ku70 and Xrcc4 gene products, with no detectable toxic effects to the cells in culture. However, similar knockdown effects using Diced siRNAs were only attained at concentrations 10-fold higher than with Stealth siRNAs. The application of RNAi technology will expand and continue to provide new insights into gene regulation and as potential applications for new therapies, transgenic animal production and basic research.

In silico identification of conserved intercoding sequences in Leishmania genomes: Unraveling putative cis-regulatory elements

In silico analyses of Leishmania spp. genome data are a powerful resource to improve the understanding of these pathogens' biology. Trypanosomatids such as Leishmania spp. have their protein-coding genes grouped in long polycistronic units of functionally unrelated genes. The control of gene expression happens by a variety of posttranscriptional mechanisms. The high degree of synteny among Leishmania species is accompanied by highly conserved coding sequences (CDS) and poorly conserved intercoding untranslated sequences. To identify the elements involved in the control of gene expression, we conducted an in silico investigation to find conserved intercoding sequences (CICS) in the genomes of L major, L infantum, and L braziliensis. We used a combination of computational tools, such as Linux-Shell, PERL and R languages, BLAST, MSPcrunch, SSAKE, and Pred-A-Term algorithms to construct a pipeline which was able to: (i) search for conservation in target-regions, (ii) eliminate CICS redundancy and mask repeat elements, (iii) predict the mRNA's extremities, (iv) analyze the distribution of orthologous genes within the generated LeishCICS-clusters, (v) assign GO terms to the LeishCICS-clusters. and (vi) provide statistical support for the gene-enrichment annotation. We associated the LeishCICS-cluster data, generated at the end of the pipeline, with the expression profile oft. donovani genes during promastigote-amastigote differentiation, as previously evaluated by others (GEO accession: GSE21936). A Pearson's correlation coefficient greater than 0.5 was observed for 730 LeishCICS-clusters containing from 2 to 17 genes. The designed computational pipeline is a useful tool and its application identified potential regulatory cis elements and putative regulons in Leishmania. (C) 2012 Elsevier B.V. All rights reserved.

MicroRNA miR-146b-5p regulates signal transduction of TGF-beta by repressing SMAD4 in thyroid cancer

MicroRNAs (miRNA) are small non-coding RNAs involved in post-transcriptional gene regulation that have crucial roles in several types of tumors, including papillary thyroid carcinoma (PTC). miR-146b-5p is overexpressed in PTCs and is regarded as a relevant diagnostic marker for this type of cancer. A computational search revealed that miR-146b-5p putatively binds to the 3' untranslated region (UTR) of SMAD4, an important member of the transforming growth factor beta (TGF-beta) signaling pathway. The TGF-beta pathway is a negative regulator of thyroid follicular cell growth, and the mechanism by which thyroid cancer cells evade its inhibitory signal remains unclear. We questioned whether the modulation of the TGF-beta pathway by miR-146b-5p can contribute to thyroid tumorigenesis. Luciferase reporter assay confirmed the direct binding of miR-146b-5p on the SMAD4 3'UTR. Specific inhibition of miR-146b-5p with a locked nucleic acid-modified anti-miR-146b oligonucleotide significantly increased SMAD4 levels in the human papillary carcinoma cell lines, TPC-1 and BCPAP. Moreover, suppression of miR-146b-5p increased the cellular response to the TGF-beta anti-proliferative signal, significantly decreasing the proliferation rate. The overexpression of miR-146b-5p in normal rat follicular PCCL3 cells decreased SMAD4 levels and disrupted TGF-beta signal transduction. MiR-146b-5p overexpression in PCCL3 cells also significantly increased cell proliferation in the absence of thyroid-stimulating hormone and conferred resistance to TGF-beta-mediated cell-cycle arrest. Additionally, the activation of thyroid most common oncogenes RET/PTC3 and BRAF in PCCL3 cells upregulated miR-146b-5p expression. Our results confirm the oncogenic role of miR-146b-5p in thyroid follicular cells and contribute to knowledge regarding the modulation of TGF-beta signal transduction by miRNAs in PTCs. Oncogene (2012) 31, 1910-1922; doi:10.1038/onc.2011.381; published online 29 August 2011

T-3 Rapidly Increases SLC2A4 Gene Expression and GLUT4 Trafficking to the Plasma Membrane in Skeletal Muscle of Rat and Improves Glucose Homeostasis

Background: Glucose transporter 4 (GLUT4) is highly expressed in muscle and fat tissue, where triiodothyronine (T-3) induces solute carrier family 2 facilitated glucose transporter member 4 (SLC2A4) gene transcription. T-3 was also shown to rapidly increase glucose uptake in myocytes exposed to cycloheximide, indicating that it might act nongenomically to regulate GLUT4 availability. We tested this hypothesis by evaluating, in thyroidectomized rats (Tx rats), the acute and/or chronic T-3 effects on GLUT4 mRNA expression and polyadenylation, protein content, and trafficking to the plasma membrane (PM) in skeletal muscle, as well as on blood glucose disappearance rate (kITT) after insulin administration. Methods: Rats were surgically thyroidectomized and treated with T-3 (0.3 to 100 mu g/100 g body weight) from 10 minutes to 5 days, and killed thereafter. Sham-operated (SO) rats were used as controls. Total RNA was extracted from the skeletal muscles (soleus [SOL] and extensorum digitalis longus [EDL]) and subjected to Northern blotting analysis using rat GLUT4 cDNA probe. Total protein was extracted and subjected to specific centrifugations for subcellular fractionation, and PM as well as microsomal (M) fractions were subjected to Western blotting analysis, using anti-GLUT4 antiserum as a probe. GLUT4 mRNA polyadenylation was examined by a rapid amplification of cDNA ends-poly(A) test (RACE-PAT). Results: Thyroidectomy reduced skeletal muscle GLUT4 mRNA, mRNA poly(A) tail length, protein content, and trafficking to the PM, as well as the kITT. The acute T-3 treatment rapidly (30 minutes) increased all these parameters compared with Tx rats. The 5-day T-3 treatment increased GLUT4 mRNA and protein expression, and restored GLUT4 trafficking to the PM and kITT to SO values. Conclusions: The results presented here show for the first time that, in parallel to its transcriptional action on the SLC2A4 gene, T-3 exerts a rapid post-transcriptional effect on GLUT4 mRNA polyadenylation, which might increase transcript stability and translation efficiency, leading to the increased GLUT4 content and availability to skeletal muscle, as well as on GLUT4 translocation to the PM, improving the insulin sensitivity, as shown by the kITT.

Exome analysis of HIV patients submitted to dendritic cells therapeutic vaccine reveals an association of CNOT1 gene with response to the treatment

INTRODUCTION: With the aim of searching genetic factors associated with the response to an immune treatment based on autologous monocyte-derived dendritic cells pulsed with autologous inactivated HIV, we performed exome analysis by screening more than 240,000 putative functional exonic variants in 18 HIV-positive Brazilian patients that underwent the immune treatment. METHODS: Exome analysis has been performed using the ILLUMINA Infinium HumanExome BeadChip. zCall algorithm allowed us to recall rare variants. Quality control and SNP-centred analysis were done with GenABEL R package. An in-house implementation of the Wang method permitted gene-centred analysis. RESULTS: CCR4-NOT transcription complex, subunit 1 (CNOT1) gene (16q21), showed the strongest association with the modification of the response to the therapeutic vaccine (p=0.00075). CNOT1 SNP rs7188697 A/G was significantly associated with DC treatment response. The presence of a G allele indicated poor response to the therapeutic vaccine (p=0.0031; OR=33.00; CI=1.74-624.66), and the SNP behaved in a dominant model (A/A vs. A/G+G/G p=0.0009; OR=107.66; 95% CI=3.85-3013.31), being the A/G genotype present only in weak/transient responders, conferring susceptibility to poor response to the immune treatment. DISCUSSION: CNOT1 is known to be involved in the control of mRNA deadenylation and mRNA decay. Moreover, CNOT1 has been recently described as being involved in the regulation of inflammatory processes mediated by tristetraprolin (TTP). The TTP-CCR4-NOT complex (CNOT1 in the CCR4-NOT complex is the binding site for TTP) has been reported as interfering with HIV replication, through post-transcriptional control. Therefore, we can hypothesize that genetic variation occurring in the CNOT1 gene could impair the TTP-CCR4-NOT complex, thus interfering with HIV replication and/or host immune response. CONCLUSIONS: Being aware that our findings are exclusive to the 18 patients studied with a need for replication, and that the genetic variant of CNOT1 gene, localized at intron 3, has no known functional effect, we propose a novel potential candidate locus for the modulation of the response to the immune treatment, and open a discussion on the necessity to consider the host genome as another potential variant to be evaluated when designing an immune therapy study

Analysis of plant LTR-retrotransposons at the fine-scale family level reveals individual molecular patterns

Background: Sugarcane is an important crop worldwide for sugar production and increasingly, as a renewable energy source. Modern cultivars have polyploid, large complex genomes, with highly unequal contributions from ancestral genomes. Long Terminal Repeat retrotransposons (LTR-RTs) are the single largest components of most plant genomes and can substantially impact the genome in many ways. It is therefore crucial to understand their contribution to the genome and transcriptome, however a detailed study of LTR-RTs in sugarcane has not been previously carried out. Results: Sixty complete LTR-RT elements were classified into 35 families within four Copia and three Gypsy lineages. Structurally, within lineages elements were similar, between lineages there were large size differences. FISH analysis resulted in the expected pattern of Gypsy/heterochromatin, Copia/euchromatin, but in two lineages there was localized clustering on some chromosomes. Analysis of related ESTs and RT-PCR showed transcriptional variation between tissues and families. Four distinct patterns were observed in sRNA mapping, the most unusual of which was that of Ale1, with very large numbers of 24nt sRNAs in the coding region. The results presented support the conclusion that distinct small RNA-regulated pathways in sugarcane target the lineages of LTR-RT elements. Conclusions: Individual LTR-RT sugarcane families have distinct structures, and transcriptional and regulatory signatures. Our results indicate that in sugarcane individual LTR-RT families have distinct behaviors and can potentially impact the genome in diverse ways. For instance, these transposable elements may affect nearby genes by generating a diverse set of small RNA's that trigger gene silencing mechanisms. There is also some evidence that ancestral genomes contribute significantly different element numbers from particular LTR-RT lineages to the modern sugarcane cultivar genome.

Inactivation of the putative suppressor gene DOK1 by promoter hypermethylation in primary human cancers

The DOK1 gene is a putative tumour suppressor gene located on the human chromosome 2p13 which is frequently rearranged in leukaemia and other human tumours. We previously reported that the DOK1 gene can be mutated and its expression down-regulated in human malignancies. However, the mechanism underlying DOK1 silencing remains largely unknown. We show here that unscheduled silencing of DOK1 expression through aberrant hypermethylation is a frequent event in a variety of human malignancies. DOK1 was found to be silenced in nine head and neck cancer (HNC) cell lines studied and DOK1 CpG hypermethylation correlated with loss of gene expression in these cells. DOK1 expression could be restored via demethylating treatment using 5-aza-2'deoxycytidine. In addition, transduction of cancer cell lines with DOK1 impaired their proliferation, consistent with the critical role of epigenetic silencing of DOK1 in the development and maintenance of malignant cells. We further observed that DOK1 hypermethylation occurs frequently in a variety of primary human neoplasm including solid tumours (93% in HNC, 81% in lung cancer) and haematopoietic malignancy (64% in Burkitt's lymphoma). Control blood samples and exfoliated mouth epithelial cells from healthy individuals showed a low level of DOK1 methylation, suggesting that DOK1 hypermethylation is a tumour specific event. Finally, an inverse correlation was observed between the level of DOK1 gene methylation and its expression in tumour and adjacent non tumour tissues. Thus, hypermethylation of DOK1 is a potentially critical event in human carcinogenesis, and may be a potential cancer biomarker and an attractive target for epigenetic-based therapy.

Human leukocyte antigen-G 3 ' untranslated region polymorphisms are associated with better kidney allograft acceptance

Human leukocyte antigen-G (FILA-G) plays a well-recognized role in the modulation of the immune response, and HLA-G expression has been associated with increased graft survival and decreased rejection episodes. To investigate the role of the HLA-G 3' untranslated region (3'UTR) in renal transplantation, we evaluated several polymorphic sites (14-bp Del/Ins +3003T/C, +3010C/G, +3027C/A, +3035C/T, +3142G/C, and +3187A/G) in patients exhibiting or not exhibiting rejection episodes. A total of 104 patients (15 with acute and 48 with chronic rejection, and 41 with no rejection) and 142 healthy individuals were studied. HLA-G 3'UTR was typed by direct sequencing. The +3035C-C genotype was more frequent in patients exhibiting chronic rejection compared with healthy controls, and the +3035C-T genotype was less frequent in chronic rejection compared with patients without rejection (acute plus chronic) or compared with healthy controls. The +3187G-A genotype, in which the A allele is associated with increased mRNA degradation, showed increased frequency in the rejection group (acute plus chronic) when compared with healthy controls. The 14 base pair Deletion/Insertion genotype was marginally increased in patients with acute rejection. This is the first study to show associations among numerous polymorphic sites in the HLA-G 3'UTR in kidney allotransplantation, which may contribute to the understanding of HLA-G post-transcriptional mechanisms. (C) 2012 American Society for Histocompatibility and Immunogenetics. Published by Elsevier Inc. All rights reserved.

Trypanosoma cruzi mitochondrial maxicircles display species- and strain-specific variation and a conserved element in the non-coding region

Abstract Background The mitochondrial DNA of kinetoplastid flagellates is distinctive in the eukaryotic world due to its massive size, complex form and large sequence content. Comprised of catenated maxicircles that contain rRNA and protein-coding genes and thousands of heterogeneous minicircles encoding small guide RNAs, the kinetoplast network has evolved along with an extreme form of mRNA processing in the form of uridine insertion and deletion RNA editing. Many maxicircle-encoded mRNAs cannot be translated without this post-transcriptional sequence modification. Results We present the complete sequence and annotation of the Trypanosoma cruzi maxicircles for the CL Brener and Esmeraldo strains. Gene order is syntenic with Trypanosoma brucei and Leishmania tarentolae maxicircles. The non-coding components have strain-specific repetitive regions and a variable region that is unique for each strain with the exception of a conserved sequence element that may serve as an origin of replication, but shows no sequence identity with L. tarentolae or T. brucei. Alternative assemblies of the variable region demonstrate intra-strain heterogeneity of the maxicircle population. The extent of mRNA editing required for particular genes approximates that seen in T. brucei. Extensively edited genes were more divergent among the genera than non-edited and rRNA genes. Esmeraldo contains a unique 236-bp deletion that removes the 5'-ends of ND4 and CR4 and the intergenic region. Esmeraldo shows additional insertions and deletions outside of areas edited in other species in ND5, MURF1, and MURF2, while CL Brener has a distinct insertion in MURF2. Conclusion The CL Brener and Esmeraldo maxicircles represent two of three previously defined maxicircle clades and promise utility as taxonomic markers. Restoration of the disrupted reading frames might be accomplished by strain-specific RNA editing. Elements in the non-coding region may be important for replication, transcription, and anchoring of the maxicircle within the kinetoplast network.

A quantitative view of the transcriptome of Schistosoma mansoni adult-worms using SAGE

Abstract Background Five species of the genus Schistosoma, a parasitic trematode flatworm, are causative agents of Schistosomiasis, a disease that is endemic in a large number of developing countries, affecting millions of patients around the world. By using SAGE (Serial Analysis of Gene Expression) we describe here the first large-scale quantitative analysis of the Schistosoma mansoni transcriptome, one of the most epidemiologically relevant species of this genus. Results After extracting mRNA from pooled male and female adult-worms, a SAGE library was constructed and sequenced, generating 68,238 tags that covered more than 6,000 genes expressed in this developmental stage. An analysis of the ordered tag-list shows the genes of F10 eggshell protein, pol-polyprotein, HSP86, 14-3-3 and a transcript yet to be identified to be the five top most abundant genes in pooled adult worms. Whereas only 8% of the 100 most abundant tags found in adult worms of S. mansoni could not be assigned to transcripts of this parasite, 46.9% of the total ditags could not be mapped, demonstrating that the 3 sequence of most of the rarest transcripts are still to be identified. Mapping of our SAGE tags to S. mansoni genes suggested the occurrence of alternative-polyadenylation in at least 13 gene transcripts. Most of these events seem to shorten the 3 UTR of the mRNAs, which may have consequences over their stability and regulation. Conclusion SAGE revealed the frequency of expression of the majority of the S. mansoni genes. Transcriptome data suggests that alternative polyadenylation is likely to be used in the control of mRNA stability in this organism. When transcriptome was compared with the proteomic data available, we observed a correlation of about 50%, suggesting that both transcriptional and post-transcriptional regulation are important for determining protein abundance in S. mansoni. The generation of SAGE tags from other life-cycle stages should contribute to reveal the dynamics of gene expression in this important parasite.

Expression profiling and cross-species RNA interference (RNAi) of desiccation-induced transcripts in the anhydrobiotic nematode Aphelenchus avenae

Abstract Background Some organisms can survive extreme desiccation by entering a state of suspended animation known as anhydrobiosis. The free-living mycophagous nematode Aphelenchus avenae can be induced to enter anhydrobiosis by pre-exposure to moderate reductions in relative humidity (RH) prior to extreme desiccation. This preconditioning phase is thought to allow modification of the transcriptome by activation of genes required for desiccation tolerance. Results To identify such genes, a panel of expressed sequence tags (ESTs) enriched for sequences upregulated in A. avenae during preconditioning was created. A subset of 30 genes with significant matches in databases, together with a number of apparently novel sequences, were chosen for further study. Several of the recognisable genes are associated with water stress, encoding, for example, two new hydrophilic proteins related to the late embryogenesis abundant (LEA) protein family. Expression studies confirmed EST panel members to be upregulated by evaporative water loss, and the majority of genes was also induced by osmotic stress and cold, but rather fewer by heat. We attempted to use RNA interference (RNAi) to demonstrate the importance of this gene set for anhydrobiosis, but found A. avenae to be recalcitrant with the techniques used. Instead, therefore, we developed a cross-species RNAi procedure using A. avenae sequences in another anhydrobiotic nematode, Panagrolaimus superbus, which is amenable to gene silencing. Of 20 A. avenae ESTs screened, a significant reduction in survival of desiccation in treated P. superbus populations was observed with two sequences, one of which was novel, while the other encoded a glutathione peroxidase. To confirm a role for glutathione peroxidases in anhydrobiosis, RNAi with cognate sequences from P. superbus was performed and was also shown to reduce desiccation tolerance in this species. Conclusions This study has identified and characterised the expression profiles of members of the anhydrobiotic gene set in A. avenae. It also demonstrates the potential of RNAi for the analysis of anhydrobiosis and provides the first genetic data to underline the importance of effective antioxidant systems in metazoan desiccation tolerance.

Identification and expression analysis of microRNAs and targets in the biofuel crop sugarcane

Abstract Background MicroRNAs (miRNAs) are small regulatory RNAs, some of which are conserved in diverse plant genomes. Therefore, computational identification and further experimental validation of miRNAs from non-model organisms is both feasible and instrumental for addressing miRNA-based gene regulation and evolution. Sugarcane (Saccharum spp.) is an important biofuel crop with publicly available expressed sequence tag and genomic survey sequence databases, but little is known about miRNAs and their targets in this highly polyploid species. Results In this study, we have computationally identified 19 distinct sugarcane miRNA precursors, of which several are highly similar with their sorghum homologs at both nucleotide and secondary structure levels. The accumulation pattern of mature miRNAs varies in organs/tissues from the commercial sugarcane hybrid as well as in its corresponding founder species S. officinarum and S. spontaneum. Using sugarcane MIR827 as a query, we found a novel MIR827 precursor in the sorghum genome. Based on our computational tool, a total of 46 potential targets were identified for the 19 sugarcane miRNAs. Several targets for highly conserved miRNAs are transcription factors that play important roles in plant development. Conversely, target genes of lineage-specific miRNAs seem to play roles in diverse physiological processes, such as SsCBP1. SsCBP1 was experimentally confirmed to be a target for the monocot-specific miR528. Our findings support the notion that the regulation of SsCBP1 by miR528 is shared at least within graminaceous monocots, and this miRNA-based post-transcriptional regulation evolved exclusively within the monocots lineage after the divergence from eudicots. Conclusions Using publicly available nucleotide databases, 19 sugarcane miRNA precursors and one new sorghum miRNA precursor were identified and classified into 14 families. Comparative analyses between sugarcane and sorghum suggest that these two species retain homologous miRNAs and targets in their genomes. Such conservation may help to clarify specific aspects of miRNA regulation and evolution in the polyploid sugarcane. Finally, our dataset provides a framework for future studies on sugarcane RNAi-dependent regulatory mechanisms.

Identification of protein-coding and non-coding RNA expression profiles in CD34+ and in stromal cells in refractory anemia with ringed sideroblasts

Abstract Background Myelodysplastic syndromes (MDS) are a group of clonal hematological disorders characterized by ineffective hematopoiesis with morphological evidence of marrow cell dysplasia resulting in peripheral blood cytopenia. Microarray technology has permitted a refined high-throughput mapping of the transcriptional activity in the human genome. Non-coding RNAs (ncRNAs) transcribed from intronic regions of genes are involved in a number of processes related to post-transcriptional control of gene expression, and in the regulation of exon-skipping and intron retention. Characterization of ncRNAs in progenitor cells and stromal cells of MDS patients could be strategic for understanding gene expression regulation in this disease. Methods In this study, gene expression profiles of CD34+ cells of 4 patients with MDS of refractory anemia with ringed sideroblasts (RARS) subgroup and stromal cells of 3 patients with MDS-RARS were compared with healthy individuals using 44 k combined intron-exon oligoarrays, which included probes for exons of protein-coding genes, and for non-coding RNAs transcribed from intronic regions in either the sense or antisense strands. Real-time RT-PCR was performed to confirm the expression levels of selected transcripts. Results In CD34+ cells of MDS-RARS patients, 216 genes were significantly differentially expressed (q-value ≤ 0.01) in comparison to healthy individuals, of which 65 (30%) were non-coding transcripts. In stromal cells of MDS-RARS, 12 genes were significantly differentially expressed (q-value ≤ 0.05) in comparison to healthy individuals, of which 3 (25%) were non-coding transcripts. Conclusions These results demonstrated, for the first time, the differential ncRNA expression profile between MDS-RARS and healthy individuals, in CD34+ cells and stromal cells, suggesting that ncRNAs may play an important role during the development of myelodysplastic syndromes.

Genome mapping and expression analyses of human intronic noncoding RNAs reveal tissue-specific patterns and enrichment in genes related to regulation of transcription

Abstract Background RNAs transcribed from intronic regions of genes are involved in a number of processes related to post-transcriptional control of gene expression. However, the complement of human genes in which introns are transcribed, and the number of intronic transcriptional units and their tissue expression patterns are not known. Results A survey of mRNA and EST public databases revealed more than 55,000 totally intronic noncoding (TIN) RNAs transcribed from the introns of 74% of all unique RefSeq genes. Guided by this information, we designed an oligoarray platform containing sense and antisense probes for each of 7,135 randomly selected TIN transcripts plus the corresponding protein-coding genes. We identified exonic and intronic tissue-specific expression signatures for human liver, prostate and kidney. The most highly expressed antisense TIN RNAs were transcribed from introns of protein-coding genes significantly enriched (p = 0.002 to 0.022) in the 'Regulation of transcription' Gene Ontology category. RNA polymerase II inhibition resulted in increased expression of a fraction of intronic RNAs in cell cultures, suggesting that other RNA polymerases may be involved in their biosynthesis. Members of a subset of intronic and protein-coding signatures transcribed from the same genomic loci have correlated expression patterns, suggesting that intronic RNAs regulate the abundance or the pattern of exon usage in protein-coding messages. Conclusion We have identified diverse intronic RNA expression patterns, pointing to distinct regulatory roles. This gene-oriented approach, using a combined intron-exon oligoarray, should permit further comparative analysis of intronic transcription under various physiological and pathological conditions, thus advancing current knowledge about the biological functions of these noncoding RNAs.

High-throughput sequencing of small RNA transcriptomes reveals critical biological features targeted by microRNAs in cell models used for squamous cell cancer research

Abstract Background The implication of post-transcriptional regulation by microRNAs in molecular mechanisms underlying cancer disease is well documented. However, their interference at the cellular level is not fully explored. Functional in vitro studies are fundamental for the comprehension of their role; nevertheless results are highly dependable on the adopted cellular model. Next generation small RNA transcriptomic sequencing data of a tumor cell line and keratinocytes derived from primary culture was generated in order to characterize the microRNA content of these systems, thus helping in their understanding. Both constitute cell models for functional studies of microRNAs in head and neck squamous cell carcinoma (HNSCC), a smoking-related cancer. Known microRNAs were quantified and analyzed in the context of gene regulation. New microRNAs were investigated using similarity and structural search, ab initio classification, and prediction of the location of mature microRNAs within would-be precursor sequences. Results were compared with small RNA transcriptomic sequences from HNSCC samples in order to access the applicability of these cell models for cancer phenotype comprehension and for novel molecule discovery. Results Ten miRNAs represented over 70% of the mature molecules present in each of the cell types. The most expressed molecules were miR-21, miR-24 and miR-205, Accordingly; miR-21 and miR-205 have been previously shown to play a role in epithelial cell biology. Although miR-21 has been implicated in cancer development, and evaluated as a biomarker in HNSCC progression, no significant expression differences were seen between cell types. We demonstrate that differentially expressed mature miRNAs target cell differentiation and apoptosis related biological processes, indicating that they might represent, with acceptable accuracy, the genetic context from which they derive. Most miRNAs identified in the cancer cell line and in keratinocytes were present in tumor samples and cancer-free samples, respectively, with miR-21, miR-24 and miR-205 still among the most prevalent molecules at all instances. Thirteen miRNA-like structures, containing reads identified by the deep sequencing, were predicted from putative miRNA precursor sequences. Strong evidences suggest that one of them could be a new miRNA. This molecule was mostly expressed in the tumor cell line and HNSCC samples indicating a possible biological function in cancer. Conclusions Critical biological features of cells must be fully understood before they can be chosen as models for functional studies. Expression levels of miRNAs relate to cell type and tissue context. This study provides insights on miRNA content of two cell models used for cancer research. Pathways commonly deregulated in HNSCC might be targeted by most expressed and also by differentially expressed miRNAs. Results indicate that the use of cell models for cancer research demands careful assessment of underlying molecular characteristics for proper data interpretation. Additionally, one new miRNA-like molecule with a potential role in cancer was identified in the cell lines and clinical samples.