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.

Expression analysis and in silico characterization of intronic long noncoding RNAs in renal cell carcinoma: emerging functional associations

Abstract Background Intronic and intergenic long noncoding RNAs (lncRNAs) are emerging gene expression regulators. The molecular pathogenesis of renal cell carcinoma (RCC) is still poorly understood, and in particular, limited studies are available for intronic lncRNAs expressed in RCC Methods Microarray experiments were performed with custom-designed arrays enriched with probes for lncRNAs mapping to intronic genomic regions. Samples from 18 primary RCC tumors and 11 nontumor adjacent matched tissues were analyzed. Meta-analyses were performed with microarray expression data from three additional human tissues (normal liver, prostate tumor and kidney nontumor samples), and with large-scale public data for epigenetic regulatory marks and for evolutionarily conserved sequences. Results A signature of 29 intronic lncRNAs differentially expressed between RCC and nontumor samples was obtained (false discovery rate (FDR) <5%). A signature of 26 intronic lncRNAs significantly correlated with the RCC five-year patient survival outcome was identified (FDR <5%, p-value ≤0.01). We identified 4303 intronic antisense lncRNAs expressed in RCC, of which 22% were significantly (p <0.05) cis correlated with the expression of the mRNA in the same locus across RCC and three other human tissues. Gene Ontology (GO) analysis of those loci pointed to 'regulation of biological processes’ as the main enriched category. A module map analysis of the protein-coding genes significantly (p <0.05) trans correlated with the 20% most abundant lncRNAs, identified 51 enriched GO terms (p <0.05). We determined that 60% of the expressed lncRNAs are evolutionarily conserved. At the genomic loci containing the intronic RCC-expressed lncRNAs, a strong association (p <0.001) was found between their transcription start sites and genomic marks such as CpG islands, RNA Pol II binding and histones methylation and acetylation. Conclusion Intronic antisense lncRNAs are widely expressed in RCC tumors. Some of them are significantly altered in RCC in comparison with nontumor samples. The majority of these lncRNAs is evolutionarily conserved and possibly modulated by epigenetic modifications. Our data suggest that these RCC lncRNAs may contribute to the complex network of regulatory RNAs playing a role in renal cell malignant transformation.

Long noncoding intronic RNAs are differentially expressed in primary and metastatic pancreatic cancer

Abstract Background Pancreatic ductal adenocarcinoma (PDAC) is known by its aggressiveness and lack of effective therapeutic options. Thus, improvement in current knowledge of molecular changes associated with pancreatic cancer is urgently needed to explore novel venues of diagnostics and treatment of this dismal disease. While there is mounting evidence that long noncoding RNAs (lncRNAs) transcribed from intronic and intergenic regions of the human genome may play different roles in the regulation of gene expression in normal and cancer cells, their expression pattern and biological relevance in pancreatic cancer is currently unknown. In the present work we investigated the relative abundance of a collection of lncRNAs in patients' pancreatic tissue samples aiming at identifying gene expression profiles correlated to pancreatic cancer and metastasis. Methods Custom 3,355-element spotted cDNA microarray interrogating protein-coding genes and putative lncRNA were used to obtain expression profiles from 38 clinical samples of tumor and non-tumor pancreatic tissues. Bioinformatics analyses were performed to characterize structure and conservation of lncRNAs expressed in pancreatic tissues, as well as to identify expression signatures correlated to tissue histology. Strand-specific reverse transcription followed by PCR and qRT-PCR were employed to determine strandedness of lncRNAs and to validate microarray results, respectively. Results We show that subsets of intronic/intergenic lncRNAs are expressed across tumor and non-tumor pancreatic tissue samples. Enrichment of promoter-associated chromatin marks and over-representation of conserved DNA elements and stable secondary structure predictions suggest that these transcripts are generated from independent transcriptional units and that at least a fraction is under evolutionary selection, and thus potentially functional. Statistically significant expression signatures comprising protein-coding mRNAs and lncRNAs that correlate to PDAC or to pancreatic cancer metastasis were identified. Interestingly, loci harboring intronic lncRNAs differentially expressed in PDAC metastases were enriched in genes associated to the MAPK pathway. Orientation-specific RT-PCR documented that intronic transcripts are expressed in sense, antisense or both orientations relative to protein-coding mRNAs. Differential expression of a subset of intronic lncRNAs (PPP3CB, MAP3K14 and DAPK1 loci) in metastatic samples was confirmed by Real-Time PCR. Conclusion Our findings reveal sets of intronic lncRNAs expressed in pancreatic tissues whose abundance is correlated to PDAC or metastasis, thus pointing to the potential relevance of this class of transcripts in biological processes related to malignant transformation and metastasis in pancreatic cancer.

Identification of sense and antisense transcripts regulated by drought in sugarcane

Sugarcane is an important sugar and energy crop that can be used efficiently for biofuels production. The development of sugarcane cultivars tolerant to drought could allow for the expansion of plantations to sub-prime regions. Knowledge on the mechanisms underlying drought responses and its relationship with carbon partition would greatly help to define routes to increase yield. In this work we studied sugarcane responses to drought using a custom designed oligonucleotide array with 21,901 different probes. The oligoarrays were designed to contain probes that detect transcription in both sense and antisense orientation. We validated the results obtained using quantitative real-time PCR (qPCR). A total of 987 genes were differentially expressed in at least one sample of sugarcane plants submitted to drought for 24, 72 and 120 h. Among them, 928 were sense transcripts and 59 were antisense transcripts. Genes related to Carbohydrate Metabolism, RNA Metabolism and Signal Transduction were selected for gene expression validation by qPCR that indicated a validation percentage of 90 %. From the probes presented on the array, 75 % of the sense probes and 11.9 % of the antisense probes have signal above background and can be classified as expressed sequences. Our custom sugarcane oligonucleotide array provides sensitivity and good coverage of sugarcane transcripts for the identification of a representative proportion of natural antisense transcripts (NATs) and sense-antisense transcript pairs (SATs). The antisense transcriptome showed, in most cases, co-expression with respective sense transcripts.

The humankind genome: from genetic diversity to the origin of human diseases

Genome-wide association studies have failed to establish common variant risk for the majority of common human diseases. The underlying reasons for this failure are explained by recent studies of resequencing and comparison of over 1200 human genomes and 10 000 exomes, together with the delineation of DNA methylation patterns (epigenome) and full characterization of coding and noncoding RNAs (transcriptome) being transcribed. These studies have provided the most comprehensive catalogues of functional elements and genetic variants that are now available for global integrative analysis and experimental validation in prospective cohort studies. With these datasets, researchers will have unparalleled opportunities for the alignment, mining, and testing of hypotheses for the roles of specific genetic variants, including copy number variations, single nucleotide polymorphisms, and indels as the cause of specific phenotypes and diseases. Through the use of next-generation sequencing technologies for genotyping and standardized ontological annotation to systematically analyze the effects of genomic variation on humans and model organism phenotypes, we will be able to find candidate genes and new clues for disease’s etiology and treatment. This article describes essential concepts in genetics and genomic technologies as well as the emerging computational framework to comprehensively search websites and platforms available for the analysis and interpretation of genomic data.

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.

Expression of Mitochondrial Non-coding RNAs (ncRNAs) Is Modulated by High Risk Human Papillomavirus (HPV) Oncogenes

The study of RNA and DNA oncogenic viruses has proved invaluable in the discovery of key cellular pathways that are rendered dysfunctional during cancer progression. An example is high risk human papillomavirus (HPV), the etiological agent of cervical cancer. The role of HPV oncogenes in cellular immortalization and transformation has been extensively investigated. We reported the differential expression of a family of human mitochondrial non-coding RNAs (ncRNAs) between normal and cancer cells. Normal cells express a sense mitochondrial ncRNA (SncmtRNA) that seems to be required for cell proliferation and two antisense transcripts (ASncmtRNAs). In contrast, the ASncmtRNAs are down-regulated in cancer cells. To shed some light on the mechanisms that trigger down-regulation of the ASncmtRNAs, we studied human keratinocytes (HFK) immortalized with HPV. Here we show that immortalization of HFK with HPV-16 or 18 causes down-regulation of the ASncmtRNAs and induces the expression of a new sense transcript named SncmtRNA-2. Transduction of HFK with both E6 and E7 is sufficient to induce expression of SncmtRNA-2. Moreover, E2 oncogene is involved in down-regulation of the ASncmtRNAs. Knockdown of E2 in immortalized cells reestablishes in a reversible manner the expression of the ASncmtRNAs, suggesting that endogenous cellular factors(s) could play functions analogous to E2 during non-HPV-induced oncogenesis.

Androgen responsive intronic non-coding RNAs

Abstract Background Transcription of large numbers of non-coding RNAs originating from intronic regions of human genes has been recently reported, but mechanisms governing their biosynthesis and biological functions are largely unknown. In this work, we evaluated the existence of a common mechanism of transcription regulation shared by protein-coding mRNAs and intronic RNAs by measuring the effect of androgen on the transcriptional profile of a prostate cancer cell line. Results Using a custom-built cDNA microarray enriched in intronic transcribed sequences, we found 39 intronic non-coding RNAs for which levels were significantly regulated by androgen exposure. Orientation-specific reverse transcription-PCR indicated that 10 of the 13 were transcribed in the antisense direction. These transcripts are long (0.5–5 kb), unspliced and apparently do not code for proteins. Interestingly, we found that the relative levels of androgen-regulated intronic transcripts could be correlated with the levels of the corresponding protein-coding gene (asGAS6 and asDNAJC3) or with the alternative usage of exons (asKDELR2 and asITGA6) in the corresponding protein-coding transcripts. Binding of the androgen receptor to a putative regulatory region upstream from asMYO5A, an androgen-regulated antisense intronic transcript, was confirmed by chromatin immunoprecipitation. Conclusion Altogether, these results indicate that at least a fraction of naturally transcribed intronic non-coding RNAs may be regulated by common physiological signals such as hormones, and further corroborate the notion that the intronic complement of the transcriptome play functional roles in the human gene-expression program.

The intronic long noncoding RNA ANRASSF1 recruits PRC2 to the RASSF1A promoter, reducing the expression of RASSF1A and increasing cell proliferation

The down-regulation of the tumor-suppressor gene RASSF1A has been shown to increase cell proliferation in several tumors. RASSF1A expression is regulated through epigenetic events involving the polycomb repressive complex 2 (PRC2); however, the molecular mechanisms modulating the recruitment of this epigenetic modifier to the RASSF1 locus remain largely unknown. Here, we identify and characterize ANRASSF1, an endogenous unspliced long noncoding RNA (lncRNA) that is transcribed from the opposite strand on the RASSF1 gene locus in several cell lines and tissues and binds PRC2. ANRASSF1 is transcribed through RNA polymerase II and is 5'-capped and polyadenylated; it exhibits nuclear localization and has a shorter half-life compared with other lncRNAs that bind PRC2. ANRASSF1 endogenous expression is higher in breast and prostate tumor cell lines compared with non-tumor, and an opposite pattern is observed for RASSF1A. ANRASSF1 ectopic overexpression reduces RASSF1A abundance and increases the proliferation of HeLa cells, whereas ANRASSF1 silencing causes the opposite effects. These changes in ANRASSF1 levels do not affect the RASSF1C isoform abundance. ANRASSF1 overexpression causes a marked increase in both PRC2 occupancy and histone H3K27me3 repressive marks, specifically at the RASSF1A promoter region. No effect of ANRASSF1 overexpression was detected on PRC2 occupancy and histone H3K27me3 at the promoter regions of RASSF1C and the four other neighboring genes, including two well-characterized tumor suppressor genes. Additionally, we demonstrated that ANRASSF1 forms an RNA/DNA hybrid and recruits PRC2 to the RASSF1A promoter. Together, these results demonstrate a novel mechanism of epigenetic repression of the RASSF1A tumor suppressor gene involving antisense unspliced lncRNA, in which ANRASSF1 selectively represses the expression of the RASSF1 isoform overlapping the antisense transcript in a location-specific manner. In a broader perspective, our findings suggest that other non-characterized unspliced intronic lncRNAs transcribed in the human genome might contribute to a location-specific epigenetic modulation of genes.

Structural synaptic elements are differentially regulated in superior temporal cortex of schizophrenia patients

Inaccurate wiring and synaptic pathology appear to be major hallmarks of schizophrenia. A variety of gene products involved in synaptic neurotransmission and receptor signaling are differentially expressed in brains of schizophrenia patients. However, synaptic pathology may also develop by improper expression of intra- and extra-cellular structural elements weakening synaptic stability. Therefore, we have investigated transcription of these elements in the left superior temporal gyrus of 10 schizophrenia patients and 10 healthy controls by genome-wide microarrays (Illumina). Fourteen up-regulated and 22 downregulated genes encoding structural elements were chosen from the lists of differentially regulated genes for further qRT-PCR analysis. Almost all genes confirmed by this method were downregulated. Their gene products belonged to vesicle-associated proteins, that is, synaptotagmin 6 and syntaxin 12, to cytoskeletal proteins, like myosin 6, pleckstrin, or to proteins of the extracellular matrix, such as collagens, or laminin C3. Our results underline the pivotal roles of structural genes that control formation and stabilization of pre- and post-synaptic elements or influence axon guidance in schizophrenia. The glial origin of collagen or laminin highlights the close interrelationship between neurons and glial cells in establishment and maintenance of synaptic strength and plasticity. It is hypothesized that abnormal expression of these and related genes has a major impact on the pathophysiology of schizophrenia.

Ionotropic Glutamate Receptors in Hypothalamic Paraventricular and Supraoptic Nuclei Mediate Vasopressin and Oxytocin Release in Unanesthetized Rats

We report changes in plasma arginine vasopressin (AVP) and oxytocin (OT) concentrations evoked by the microinjection of L-glutamate (L-glu) into the hypothalamic supraoptic nucleus (SON) and paraventricular nucleus(PVN) of unanesthetized rats, as well as which local mechanisms are involved in their mediation. L-Glu microinjection (10 nmol/100 nl) into the SON increased the circulating levels of both AVP and OT. The AVP increases were blocked by local pretreatment with the selective non-N-methyl-D-aspartate (NMDA) receptor antagonist 2,3-dioxo-6-nitro-1,2,3,4-tetrahydrobenzo[f]quinoxaline-7-sulfonamide (NBQX) (2 nmol/100 nl), but it was not affected by pretreatment with the NMDA-receptor antagonist LY235959 (2 nmol/100 nl). The OT response to L-glu microinjection into the SON was blocked by local pretreatment with either NBQX or LY235959. Furthermore, the administration of either the non-NMDA receptor agonist (+/-)-alpha-amino-3-hydroxy-5-methylisoxazole-4-propionic acid hydrobromide (AMPA) (5 nmol/100 nl) or NMDA receptor agonist NMDA (5 nmol/100 nl) into the SON had no effect on OT baseline plasma levels, but when both agonists were microinjected together these levels were increased. L-Glu microinjection into the PVN did not change circulating levels of either AVP or OT. However, after local pretreatment with LY235959, the L-glu microinjection increased plasma levels of the hormones. The L-glu microinjection into the PVN after the local treatment with NBQX did not affect the circulating AVP and OT levels. Therefore, results suggest the AVP release from the SON is mediated by activation of non-NMDA glutamate receptors, whereas the OT release from this nucleus is mediated by an interaction of NMDA and non-NMDA receptors. The present study also suggests an inhibitory role for NMDA receptors in the PVN on the release of AVP and OT. (Endocrinology 153: 2323-2331, 2012)

Maternal exposure to picrotoxin modifies the response of the GABA(A) receptor during sexual behavior of adult male rat offspring

This study investigated whether perinatal exposure to picrotoxin, a GABA(A) antagonist, modifies the effect of muscimol, a GABA(A) agonist, on the sexual behavior of adult male rats. Two hours after birth and then once daily during the next 9 days of lactation, dams received picrotoxin (0.75 mg/kg subcutaneously) or saline (1 ml/kg subcutaneously). The adult male offspring from the picrotoxin and saline groups received saline (1 ml/kg intraperitoneally) or muscimol (1 mg/kg intraperitoneally), and 15 min later, their sexual behavior was assessed. Muscimol treatment in the saline-exposed group increased the mount and intromission latencies. However, these effects were absent in the picrotoxin-exposed groups. The latencies to first ejaculation, postejaculatory mount, and intromission were decreased in both picrotoxin-exposed groups relative to the saline-exposed groups. The picrotoxin + muscimol-treated rats required more intromissions to ejaculate and the picrotoxin-exposed groups made more ejaculations than the saline-exposed groups. Thus, muscimol treatment did not increase the mount and intromission latencies following picrotoxin exposure, but increased the ejaculation frequency, which did not differ between the picrotoxin + muscimol and the picrotoxin + saline groups. These data indicate that perinatal picrotoxin treatment interfered with GABA(A) receptor development Behavioural Pharmacology 23:703-709 (c) 2012 Wolters Kluwer Health vertical bar Lippincott Williams & Wilkins.

MicroRNAs miR-146-5p and let-7f as prognostic tools for aggressive papillary thyroid carcinoma: a case report

Papillary thyroid cancer (PTC) is the most incident histotype of thyroid cancer. A certain fraction of PTC cases (5%) are irresponsive to conventional treatment, and refractory to radioiodine therapy. The current prognostic factors for aggressiveness are mainly based on tumor size, the presence of lymph node metastasis, extrathyroidal invasion and, more recently, the presence of the BRAFT(1799A) mutation. MicroRNAs (miRNAs) have been described as promising molecular markers for cancer as their deregulation is observed in a wide range of tumors. Recent studies indicate that the over-expression of miR-146b-5p is associated with aggressiveness and BRAFT(1799A) mutation. Furthermore, down-regulation of let-7f is observed in several types of tumors, including PTC. In this study, we evaluated the miR146b-5p and let-7f status in a young male patient with aggressive, BRAFT(1799A)-positive papillary thyroid carcinoma, with extensive lymph node metastases and short-time recurrence. The analysis of miR-146b-5p and let-7f expression revealed a distinct pattern from a cohort of PTC patients, suggesting caution in evaluating miRNA expression data as molecular markers of PTC diagnosis and prognosis. Arq Bras Endocrinol Metab. 2012;56(8):552-7

Analysis of energetically biased transcripts of viruses and transposable elements

RNA interference (RNAi) is a natural endogenous process by which double-stranded RNA molecules trigger potent and specific gene silencing in eukaryotic cells and is characterized by target RNA cleavage. In mammals, small interfering RNAs (siRNAs) are the trigger molecules of choice and constitute a new class of RNA-based antiviral agents. In an efficient RNAi response, the antisense strand of siRNAs must enter the RNA-induced silencing complex (RISC) in a process mediated by thermodynamic features. In this report, we hypothesize that silent mutations capable of inverting thermodynamic properties can promote resistance to siRNAs. Extensive computational analyses were used to assess whether continuous selective pressure that promotes such mutations could lead to the emergence of viral strains completely resistant to RNAi (i.e., prone to transfer only the sense strands to RISC). Based on our findings, we propose that, although synonymous mutations may produce functional resistance, this strategy cannot be systematically adopted by viruses since the longest RNAi-refractory sequence is only 10 nt long. This finding also suggests that all mRNAs display fluctuating thermodynamic landscapes and that, in terms of thermodynamic features, RNAi is a very efficient antiviral system since there will always be sites susceptible to siRNAs.

Effect of human TGF-beta on the gene expression profile of Schistosoma mansoni adult worms

Schistosoma mansoni is responsible for schistosomiasis, a parasitic disease that affects 200 million people worldwide. Molecular mechanisms of host-parasite interaction are complex and involve a crosstalk between host signals and parasite receptors. TGF-beta signaling pathway has been shown to play an important role in S. mansoni development and embryogenesis. In particular human (h) TGF-beta has been shown to bind to a S. mansoni receptor, transduce a signal that regulates the expression of a schistosome target gene. Here we describe 381 parasite genes whose expression levels are affected by in vitro treatment with hTGF-beta. Among these differentially expressed genes we highlight genes related to morphology, development and cell cycle that could be players of cytokine effects on the parasite. We confirm by qPCR the expression changes detected with microarrays for 5 out of 7 selected genes. We also highlight a set of non-coding RNAs transcribed from the same loci of protein-coding genes that are differentially expressed upon hTCF-beta treatment. These datasets offer potential targets to be explored in order to understand the molecular mechanisms behind the possible role of hTGF-beta effects on parasite biology. (C) 2012 Elsevier B.V. All rights reserved.