A score system for quality evaluation of RNA sequence tags: an improvement for gene expression profiling

Background: High-throughput molecular approaches for gene expression profiling, such as Serial Analysis of Gene Expression (SAGE), Massively Parallel Signature Sequencing (MPSS) or Sequencing-by-Synthesis (SBS) represent powerful techniques that provide global transcription profiles of different cell types through sequencing of short fragments of transcripts, denominated sequence tags. These techniques have improved our understanding about the relationships between these expression profiles and cellular phenotypes. Despite this, more reliable datasets are still necessary. In this work, we present a web-based tool named S3T: Score System for Sequence Tags, to index sequenced tags in accordance with their reliability. This is made through a series of evaluations based on a defined rule set. S3T allows the identification/selection of tags, considered more reliable for further gene expression analysis. Results: This methodology was applied to a public SAGE dataset. In order to compare data before and after filtering, a hierarchical clustering analysis was performed in samples from the same type of tissue, in distinct biological conditions, using these two datasets. Our results provide evidences suggesting that it is possible to find more congruous clusters after using S3T scoring system. Conclusion: These results substantiate the proposed application to generate more reliable data. This is a significant contribution for determination of global gene expression profiles. The library analysis with S3T is freely available at http://gdm.fmrp.usp.br/s3t/.S3T source code and datasets can also be downloaded from the aforementioned website.

Alternative splicing enriched cDNA libraries identify breast cancer-associated transcripts

Background: Alternative splicing (AS) is a central mechanism in the generation of genomic complexity and is a major contributor to transcriptome and proteome diversity. Alterations of the splicing process can lead to deregulation of crucial cellular processes and have been associated with a large spectrum of human diseases. Cancer-associated transcripts are potential molecular markers and may contribute to the development of more accurate diagnostic and prognostic methods and also serve as therapeutic targets. Alternative splicing-enriched cDNA libraries have been used to explore the variability generated by alternative splicing. In this study, by combining the use of trapping heteroduplexes and RNA amplification, we developed a powerful approach that enables transcriptome-wide exploration of the AS repertoire for identifying AS variants associated with breast tumor cells modulated by ERBB2 (HER-2/neu) oncogene expression. Results: The human breast cell line (C5.2) and a pool of 5 ERBB2 over-expressing breast tumor samples were used independently for the construction of two AS-enriched libraries. In total, 2,048 partial cDNA sequences were obtained, revealing 214 alternative splicing sequence-enriched tags (ASSETs). A subset with 79 multiple exon ASSETs was compared to public databases and reported 138 different AS events. A high success rate of RT-PCR validation (94.5%) was obtained, and 2 novel AS events were identified. The influence of ERBB2-mediated expression on AS regulation was evaluated by capillary electrophoresis and probe-ligation approaches in two mammary cell lines (Hb4a and C5.2) expressing different levels of ERBB2. The relative expression balance between AS variants from 3 genes was differentially modulated by ERBB2 in this model system. Conclusions: In this study, we presented a method for exploring AS from any RNA source in a transcriptome-wide format, which can be directly easily adapted to next generation sequencers. We identified AS transcripts that were differently modulated by ERBB2-mediated expression and that can be tested as molecular markers for breast cancer. Such a methodology will be useful for completely deciphering the cancer cell transcriptome diversity resulting from AS and for finding more precise molecular markers.

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

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.

Noncoding RNA's in mammals

Transcripts that lack any protein-coding potential represent at least half of the identified elements transcriptome. We review the evidence for the existence of such transcripts in the mammalian transcriptome, and argue that there may be many more noncoding RNAs (ncRNAs) still to be discovered. Relatively few ncRNA “genes” have been ascribed a function based upon mutation analysis. The review discusses possible roles of ncRNAs as cis-acting and trans-acting elements in epigenetic transcriptional control, including monoallelic gene silencing and imprinting. We also consider the evidence that the production of ncRNAs is a common feature of transcriptional enhancers.

Quantitation of alternatively spliced NMDA receptor NR1 isoform mRNA transcripts in human brain by competitive RT-PCR

The N-methyl-D-aspartate (NMDA)-selective subtype of ionotropic glutamate receptor is of importance in neuronal differentiation and synapse consolidation, activity-dependent forms of synaptic plasticity, and excitatory amino acid-mediated neuronal toxicity [Neurosci. Res. Program, Bull. 19 (1981) 1; Lab. Invest. 68 (1993) 372]. NMDA receptors exist in vivo as tetrameric or pentameric complexes comprising proteins from two families of homologous subunits, designated NR1 and NR2(A-D) [Biochem. Biophys. Res. Commun. 185 (1992) 826]. The gene coding for the human NR1 subunit (hNR1) is composed of 21 exons, three of which (4, 20 and 21) can be differentially spliced to generate a total of eight distinct subunit variants. We detail here a competitive RT-PCR (cRT-PCR) protocol to quantify endogenous levels of hNR1 splice variants in autopsied human brain. Quantitation of each hNR1 splice variant is performed using standard curve methodology in which a known amount of synthetic ribonucleic acid competitor (internal standard) is co-amplified against total RNA. This method can be used for the quantitation of hNR1 mRNA levels in response to acute or chronic disease states, in particular in the glutamatergic-associated neuronal loss observed in Alzheimer's disease [J. Neurochem. 78 (2001) 175]. Furthermore, alterations in hNR1 mRNA expression may be reflected at the translational level, resulting in functional changes in the NMDA receptor. (C) 2003 Elsevier Science B.V. All rights reserved.

The expression of tubulin cofactor A (TBCA) is regulated by a noncoding antisense Tbca RNA during testis maturation

Abstract - Recently, long noncoding RNAs have emerged as pivotal molecules for the regulation of coding genes' expression. These molecules might result from antisense transcription of functional genes originating natural antisense transcripts (NATs) or from transcriptional active pseudogenes. TBCA interacts with β-tubulin and is involved in the folding and dimerization of new tubulin heterodimers, the building blocks of microtubules. Methodology/Principal findings: We found that the mouse genome contains two structurally distinct Tbca genes located in chromosomes 13 (Tbca13) and 16 (Tbca16). Interestingly, the two Tbca genes albeit ubiquitously expressed, present differential expression during mouse testis maturation. In fact, as testis maturation progresses Tbca13 mRNA levels increase progressively, while Tbca16 mRNA levels decrease. This suggests a regulatory mechanism between the two genes and prompted us to investigate the presence of the two proteins. However, using tandem mass spectrometry we were unable to identify the TBCA16 protein in testis extracts even in those corresponding to the maturation step with the highest levels of Tbca16 transcripts. These puzzling results led us to re-analyze the expression of Tbca16. We then detected that Tbca16 transcription produces sense and natural antisense transcripts. Strikingly, the specific depletion by RNAi of these transcripts leads to an increase of Tbca13 transcript levels in a mouse spermatocyte cell line. Conclusions/Significance: Our results demonstrate that Tbca13 mRNA levels are post-transcriptionally regulated by the sense and natural antisense Tbca16 mRNA levels. We propose that this regulatory mechanism operates during spermatogenesis, a process that involves microtubule rearrangements, the assembly of specific microtubule structures and requires critical TBCA levels.

Results of An International Study of Quantitative BCR/ABL Assays Using Defined RNA Samples

Molecular monitoring of BCR/ABL transcripts by real time quantitative reverse transcription PCR (qRT-PCR) is an essential technique for clinical management of patients with BCR/ABL-positive CML and ALL. Though quantitative BCR/ABL assays are performed in hundreds of laboratories worldwide, results among these laboratories cannot be reliably compared due to heterogeneity in test methods, data analysis, reporting, and lack of quantitative standards. Recent efforts towards standardization have been limited in scope. Aliquots of RNA were sent to clinical test centers worldwide in order to evaluate methods and reporting for e1a2, b2a2, and b3a2 transcript levels using their own qRT-PCR assays. Total RNA was isolated from tissue culture cells that expressed each of the different BCR/ABL transcripts. Serial log dilutions were prepared, ranging from 100 to 10-5, in RNA isolated from HL60 cells. Laboratories performed 5 independent qRT-PCR reactions for each sample type at each dilution. In addition, 15 qRT-PCR reactions of the 10-3 b3a2 RNA dilution were run to assess reproducibility within and between laboratories. Participants were asked to run the samples following their standard protocols and to report cycle threshold (Ct), quantitative values for BCR/ABL and housekeeping genes, and ratios of BCR/ABL to housekeeping genes for each sample RNA. Thirty-seven (n=37) participants have submitted qRT-PCR results for analysis (36, 37, and 34 labs generated data for b2a2, b3a2, and e1a2, respectively). The limit of detection for this study was defined as the lowest dilution that a Ct value could be detected for all 5 replicates. For b2a2, 15, 16, 4, and 1 lab(s) showed a limit of detection at the 10-5, 10-4, 10-3, and 10-2 dilutions, respectively. For b3a2, 20, 13, and 4 labs showed a limit of detection at the 10-5, 10-4, and 10-3 dilutions, respectively. For e1a2, 10, 21, 2, and 1 lab(s) showed a limit of detection at the 10-5, 10-4, 10-3, and 10-2 dilutions, respectively. Log %BCR/ABL ratio values provided a method for comparing results between the different laboratories for each BCR/ABL dilution series. Linear regression analysis revealed concordance among the majority of participant data over the 10-1 to 10-4 dilutions. The overall slope values showed comparable results among the majority of b2a2 (mean=0.939; median=0.9627; range (0.399 - 1.1872)), b3a2 (mean=0.925; median=0.922; range (0.625 - 1.140)), and e1a2 (mean=0.897; median=0.909; range (0.5174 - 1.138)) laboratory results (Fig. 1-3)). Thirty-four (n=34) out of the 37 laboratories reported Ct values for all 15 replicates and only those with a complete data set were included in the inter-lab calculations. Eleven laboratories either did not report their copy number data or used other reporting units such as nanograms or cell numbers; therefore, only 26 laboratories were included in the overall analysis of copy numbers. The median copy number was 348.4, with a range from 15.6 to 547,000 copies (approximately a 4.5 log difference); the median intra-lab %CV was 19.2% with a range from 4.2% to 82.6%. While our international performance evaluation using serially diluted RNA samples has reinforced the fact that heterogeneity exists among clinical laboratories, it has also demonstrated that performance within a laboratory is overall very consistent. Accordingly, the availability of defined BCR/ABL RNAs may facilitate the validation of all phases of quantitative BCR/ABL analysis and may be extremely useful as a tool for monitoring assay performance. Ongoing analyses of these materials, along with the development of additional control materials, may solidify consensus around their application in routine laboratory testing and possible integration in worldwide efforts to standardize quantitative BCR/ABL testing.

RNA-driven cyclin-dependent kinase regulation: when CDK9/cyclin T subunits of P-TEFb meet their ribonucleoprotein partners.

The positive transcription elongation factor (P-TEFb) consists of CDK9, a cyclin-dependent kinase and its cyclin T partner. It is required for transcription of most class II genes. Its activity is regulated by non-coding RNAs. The 7SK cellular RNA turns the HEXIM cellular protein into a P-TEFb inhibitor that binds its cyclin T subunit. Thus, P-TEFb activity responds to variations in global cellular transcriptional activity and to physiological conditions linked to cell differentiation, proliferation or cardiac hypertrophy. In contrast, the Tat activation region RNA plays an activating role. This feature at the 5' end of the human immunodeficiency (HIV) viral transcript associates with the viral protein Tat that in turn binds cyclin T1 and recruits active P-TEFb to the HIV promoter. This results in enhanced P-TEFb activity, which is critical for an efficient production of viral transcripts. Although discovered recently, the regulation of P-TEFb becomes a paradigm for non-coding RNAs that regulate transcription factors. It is also a unique example of RNA-driven regulation of a cyclindependent kinase.

Widespread occurrence of non-canonical transcription termination by human RNA polymerase III.

Human RNA polymerase (Pol) III-transcribed genes are thought to share a simple termination signal constituted by four or more consecutive thymidine residues in the coding DNA strand, just downstream of the RNA 3'-end sequence. We found that a large set of human tRNA genes (tDNAs) do not display any T(≥4) stretch within 50 bp of 3'-flanking region. In vitro analysis of tDNAs with a distanced T(≥4) revealed the existence of non-canonical terminators resembling degenerate T(≥5) elements, which ensure significant termination but at the same time allow for the production of Pol III read-through pre-tRNAs with unusually long 3' trailers. A panel of such non-canonical signals was found to direct transcription termination of unusual Pol III-synthesized viral pre-miRNA transcripts in gammaherpesvirus 68-infected cells. Genome-wide location analysis revealed that human Pol III tends to trespass into the 3'-flanking regions of tDNAs, as expected from extensive terminator read-through. The widespread occurrence of partial termination suggests that the Pol III primary transcriptome in mammals is unexpectedly enriched in 3'-trailer sequences with the potential to contribute novel functional ncRNAs.

Production of low phytic acid rice by hairpin RNA- and artificial microRNA-mediated silencing of OsMIK in seeds

To produce agronomically competitive rice with nutritionally superior, environmentally safe phytic acid (PA) levels, hairpin RNA (hpRNA)- and artificial microRNA (amiRNA)-mediated gene silencing approaches were explored to reduce both myo-inositol kinase gene (OsMIK) expression and PA accumulation in rice seeds. hpRNA and amiRNA sequences targeted to OsMIK (hpMIK and amiMIK), under the control of a rice Ole18 promoter, were transformed into the rice cultivar Nippon-bare. Fourteen and 21 independent transgenic events were identified containing the hpMIK and amiMIK constructs, respectively, from which five stable homozygous transgenic lines of each were developed together with their null siblings. Southern blotting demonstrated transgene integration into the genome and quantitative real-time PCR showed that gene silencing was restricted to seeds. OsMIK transcripts were significantly reduced in both transgenic amiMIK and hpMIK seeds, which had PA levels reduced by 14.9-50.2 and 38.1-50.7 %, respectively, compared with their respective null siblings. There were no systematic significant differences in agronomic traits between the transgenic lines and their non-transgenic siblings, and no correlation between seed PA contents and decreased rates of seed germination and seedling emergence. The results of the present study suggest that Ole 18-driven OsMIK silencing via hpRNA and amiRNA could be an effective way to develop agronomically competitive low phytic acid rice.

Differential expression of notch signaling-related transcripts accompanies Pro-thymocyte proliferation and phenotype transition induced by epidermal growth factor plus insulin in fetal thymus organ cultures

Thymus regression upon stressing stimuli, such as infectious diseases, is followed by organ reconstitution, paralleling its development in ontogeny. A narrow window of thymus development was here studied, encompassing the pro-T lymphoid precursor expansion during specification stages, by the use of epidermal growth factor plus insulin (INS) in murine fetal thymus organ cultures. Aiming to disclose signaling pathways related to these stages, cultured thymus lobes had their RNA extracted, for the search of transcripts differentially expressed using RNAse protection assays and reverse transcriptase-polymerase chain reactions. We found no difference that could explain INS-driven thymocyte growth, in the pattern of transcripts for death/proliferation mediators, or for a series of growth factor receptors and transcriptional regulators known as essential for thymus development. Thymocyte suspensions from cultured lobes, stained for phenotype analysis by fluorescence activated cell sorting, showed a decreased staining for Notch1 protein at cell surfaces upon INS addition. We analyzed the expression of Notch-related elements, and observed the recruitment of a specific set of transcripts simultaneous and compatible with INS-driven thymocyte growth, namely, transcripts for Notch3, for its ligand Jagged2, and for Deltex1, a mediator of a poorly characterized alternative pathway downstream of the Notch receptor.

Trypanosoma cruzi: establishment of permeable cells for RNA processing analysis with drugs

Pre-mRNA maturation in trypanosomatids occurs through a process called trans-splicing which involves excision of introns and union of exons in two independent transcripts. For the first time, we present the standardization of Trypanosoma cruzi permeable cells (Y strain) as a model for trans-splicing study of mRNAs in trypanosomes, following by RNase protection reaction, which localizes the SL exon and intron. This trans-splicing reaction in vitro was also used to analyze the influence of NFOH-121, a nitrofurazone-derivative, on this mechanism. The results suggested that the prodrug affects the RNA processing in these parasites, but the trans-splicing reaction still occurred.

RT-PCR diagnosis of patients with acute nonlymphocytic leukemia and inv(16)(p13q22) and identification of new alternative splicing in CBFB-MYH11 transcripts.

As acute nonlymphocytic leukemia (ANLL) with inv(16) (p13q22) or t(16;16)(p13;q22) has been shown to result from the fusion of transcription factor subunit core binding factor (CBFB) to a myosin heavy chain (MYH11), we sought to design methods to detect this rearrangement using reverse transcriptase-polymerase chain reaction (RT-PCR). In all of 27 inv(16)(p13q22) and four t(16;16)(p13;q22) cases tested, a chimeric CBFB-MYH11 transcript coding for an in-frame fusion protein was detected. In a more extensive RT-PCR analysis with different primer pairs, we detected a second new chimeric CBFB-MYH11 transcript in 10 of 11 patients tested. The CBFB-MYH11 reading frame of the second transcript was maintained in one patient but not in the others. We show that the different CBFB-MYH11 transcripts in one patient arise from alternative splicing. Translation of the transcript in which the CBFB-MYH11 reading frame is not maintained leads to a slightly truncated CBFB protein.

A rice cis-natural antisense RNA acts as a translational enhancer for its cognate mRNA and contributes to phosphate homeostasis and plant fitness.

cis-natural antisense transcripts (cis-NATs) are widespread in plants and are often associated with downregulation of their associated sense genes. We found that a cis-NAT positively regulates the level of a protein critical for phosphate homeostasis in rice (Oryza sativa). PHOSPHATE1;2 (PHO1;2), a gene involved in phosphate loading into the xylem in rice, and its associated cis-NATPHO1;2 are both controlled by promoters active in the vascular cylinder of roots and leaves. While the PHO1;2 promoter is unresponsive to the plant phosphate status, the cis-NATPHO1;2 promoter is strongly upregulated under phosphate deficiency. Expression of both cis-NATPHO1;2 and the PHO1;2 protein increased in phosphate-deficient plants, while the PHO1;2 mRNA level remained stable. Downregulation of cis-NATPHO1;2 expression by RNA interference resulted in a decrease in PHO1;2 protein, impaired the transfer of phosphate from root to shoot, and decreased seed yield. Constitutive overexpression of NATPHO1;2 in trans led to a strong increase of PHO1;2, even under phosphate-sufficient conditions. Under all conditions, no changes occurred in the level of expression, sequence, or nuclear export of PHO1;2 mRNA. However, expression of cis-NATPHO1;2 was associated with a shift of both PHO1;2 and cis-NATPHO1;2 toward the polysomes. These findings reveal an unexpected role for cis-NATPHO1;2 in promoting PHO1;2 translation and affecting phosphate homeostasis and plant fitness.

A directed approach for the identification of transcripts harbouring the spliced leader sequence and the effect of trans-splicing knockdown in Schistosoma mansoni

Schistosomiasis is a major neglected tropical disease caused by trematodes from the genus Schistosoma. Because schistosomes exhibit a complex life cycle and numerous mechanisms for regulating gene expression, it is believed that spliced leader (SL) trans-splicing could play an important role in the biology of these parasites. The purpose of this study was to investigate the function of trans-splicing in Schistosoma mansoni through analysis of genes that may be regulated by this mechanism and via silencing SL-containing transcripts through RNA interference. Here, we report our analysis of SL transcript-enriched cDNA libraries from different S. mansoni life stages. Our results show that the trans-splicing mechanism is apparently not associated with specific genes, subcellular localisations or life stages. In cross-species comparisons, even though the sets of genes that are subject to SL trans-splicing regulation appear to differ between organisms, several commonly shared orthologues were observed. Knockdown of trans-spliced transcripts in sporocysts resulted in a systemic reduction of the expression levels of all tested trans-spliced transcripts; however, the only phenotypic effect observed was diminished larval size. Further studies involving the findings from this work will provide new insights into the role of trans-splicing in the biology of S. mansoni and other organisms. All Expressed Sequence Tags generated in this study were submitted to dbEST as five different libraries. The accessions for each library and for the individual sequences are as follows: (i) adult worms of mixed sexes (LIBEST_027999: JZ139310 - JZ139779), (ii) female adult worms (LIBEST_028000: JZ139780 - JZ140379), (iii) male adult worms (LIBEST_028001: JZ140380 - JZ141002), (iv) eggs (LIBEST_028002: JZ141003 - JZ141497) and (v) schistosomula (LIBEST_028003: JZ141498 - JZ141974).