Deep RNA sequencing at single base-pair resolution reveals high complexity of the rice transcriptome

Understanding the dynamics of eukaryotic transcriptome is essential for studying the complexity of transcriptional regulation and its impact on phenotype. However, comprehensive studies of transcriptomes at single base resolution are rare, even for modern organisms, and lacking for rice. Here, we present the first transcriptome atlas for eight organs of cultivated rice. Using high-throughput paired-end RNA-seq, we unambiguously detected transcripts expressing at an extremely low level, as well as a substantial number of novel transcripts, exons, and untranslated regions. An analysis of alternative splicing in the rice transcriptome revealed that alternative cis-splicing occurred in similar to 33% of all rice genes. This is far more than previously reported. In addition, we also identified 234 putative chimeric transcripts that seem to be produced by trans-splicing, indicating that transcript fusion events are more common than expected. In-depth analysis revealed a multitude of fusion transcripts that might be by-products of alternative splicing. Validation and chimeric transcript structural analysis provided evidence that some of these transcripts are likely to be functional in the cell. Taken together, our data provide extensive evidence that transcriptional regulation in rice is vastly more complex than previously believed.

Knock down of gfp and no tail expression in zebrafish embryo by in vivo-transcribed short hairpin RNA with T7 plasmid system

A short-hairpin RNA (shRNA) expression system, based on T7 RNA polymerase (T7RP) directed transcription machinery, has been developed and used to generate a knock down effect in zebrafish embryos by targeting green fluorescent protein (gfp) and no tail (ntl) mRNA. The vector pCMVT7R harboring T7RP driven by CMV promoter was introduced into zebrafish embryos and the germline transmitted transgenic individuals were screened out for subsequent RNAi application. The shRNA transcription vectors pT7shRNA were constructed and validated by in vivo transcription assay. When pT7shGFP vector was injected into the transgenic embryos stably expressing T7RP, gfp relative expression level showed a decrease of 68% by analysis of fluorescence real time RT-PCR. As a control, injection of chemical synthesized siRNA resulted in expression level of 40% lower than the control when the injection dose was as high as 2 mu g/mu l. More importantly, injection of pT7shNTL vector in zebrafish embryos expressing T7RP led to partial absence of endogenous ntl transcripts in 30% of the injected embryos when detected by whole mount in situ hybridization. Herein, the T7 transcription system could be used to drive the expression of shRNA in zebrafish embryos and result in gene knock down effect, suggesting a potential role for its application in RNAi studies in zebrafish embryos.

Molecular characterization and effect of RNA interference of retinoid X receptor (RXR) on E75 and chitinase gene expression in Chinese shrimp Fenneropenaeus chinensis

Retinoid X receptor (RXR)/ultraspiracle (USP) is the heterodimeric partner of ecdysteroid receptor and is required for the molting process of arthropods. To better understand the molecular aspects governing the process of molting in shrimp, the full-length cDNA of two RXRs, named as FcRXR-1 and FcRXR-2 were obtained from Chinese shrimp Fenneropenaeus chinensis which were of 1715 and 1700 bp long, revealed a 1315 and 1300 bp open reading frame (ORF) respectively. Quantitative Real time PCR analysis showed a marked tissue-specific difference in the expression of FcRXR transcript, which revealed that the expression of FcRXR Could be regulated in a tissue-specific manner. Moreover, high expression of FcRXR mRNAs was observed in late pre-molt period (D3) and post molt stages (A-B) of shrimp. Among the two isoforms, FcRXR-2 appeared in a considerably high level in all the stages compared to the FcRXR-1. In addition, we examined the temporal expression of two chitinase genes: FcChitinase (FcChi) and FcChitinase-1 (FcChi-1) during the molt cycle of F chinensis. Both the FcChi and FcChi-1 transcripts were detected in all stages of molting, although considerable fluctuations observed through the molt cycle. Injection of double stranded RXR (dsRXR) into juvenile shrimp resulted in a maximum silencing effect at 48 h post injection. We analyzed the expression levels of FcChi, FcChi-1 and the ecdysone inducible gene E75 (FcE75) in samples of dsRXR injected shrimp. Significant reduction in levels of both FcE75, FcChi and FcChi-1 transcripts Occurred in the silenced shrimp. This correlation suggested that RXR might involve in the downstream regulation of E75 and chitinase gene transcription in the ecdysone signaling pathway of decapod crustaceans. (C) 2009 Published by Elsevier Inc.

The prevalence and regulation of antisense transcripts in Schizosaccharomyces pombe.

A strand-specific transcriptome sequencing strategy, directional ligation sequencing or DeLi-seq, was employed to profile antisense transcriptome of Schizosaccharomyces pombe. Under both normal and heat shock conditions, we found that polyadenylated antisense transcripts are broadly expressed while distinct expression patterns were observed for protein-coding and non-coding loci. Dominant antisense expression is enriched in protein-coding genes involved in meiosis or stress response pathways. Detailed analyses further suggest that antisense transcripts are independently regulated with respect to their sense transcripts, and diverse mechanisms might be potentially involved in the biogenesis and degradation of antisense RNAs. Taken together, antisense transcription may have profound impacts on global gene regulation in S. pombe.

Short interfering RNA-mediated gene silencing in Globodera pallida and Meloidogyne incognita infective stage juveniles

The analysis of gene function through RNA interference (RNAi)-based reverse genetics in plant parasitic nematodes (PPNs) remains inexplicably reliant on the use of long double-stranded RNA (dsRNA) silencing triggers; a practice inherently disadvantageous due to the introduction of superfluous dsRNA sequence. increasing chances of aberrant or off-target gene silencing through interactions between nascent short interfering RNAs (siRNAs) and non-cognate mRNA targets. Recently, we have shown that non-nematode, long dsRNAs have a propensity to elicit profound impacts on the phenotype and migrational abilities of both root knot and cyst nematodes. This study presents, to our knowledge for the first time, gene-specific knockdown of FMRFamide-like peptide (flp) transcripts, using discrete 21 bp siRNAs in potato cyst nematode Globodera pallida, and root knot nematode Meloidogyne incognita infective (J2) stage juveniles. Both knockdown at the transcript level through quantitative (q)PCR analysis and functional data derived from migration assay, indicate that siRNAs targeting certain areas of the FMRFamide-like peptide (FLP) transcripts are potent and specific in the silencing of gene function. In addition, we present a method of manipulating siRNA activity through the management of strand thermodynamics. Initial evaluation of strand thermodynamics as a determinant of RNA-induced Silencing Complex (RISC) strand selection (inferred from knockdown efficacy) in the siRNAs presented here suggested that the purported influence of 5' stand stability on guide incorporation may be somewhat promiscuous. However, we have found that on strategically incorporating base mismatches in the sense strand of a G. pallida-specific siRNA we could specifically increase or decrease the knockdown of its target (specific to the antisense strand), presumably through creating more favourable thermodynamic profiles for incorporation of either the sense (non-target-specific) or antisense (target-specific) strand into a cleavage-competent RISC. Whilst the efficacy of similar approaches to siRNA modification has been demonstrated in the context of Drosophila whole-cell lysate preparations and in mammalian cell cultures, it remained to be seen how these sense strand mismatches may impact on gene silencing in vivo, in relation to different targets and in different sequence contexts. This work presents the first application of such an approach in a whole organism; initial results show promise. (C) 2009 Australian Society for Parasitology Inc. Published by Elsevier Ltd. All rights reserved.

Short interfering RNA-mediated knockdown of drosha and pasha in undifferentiated Meloidogyne incognita eggs leads to irregular growth and embryonic lethality

Micro-(mi)RNAs play a pivotal role in the developmental regulation of plants and animals. We reasoned that disruption of normal heterochronic activity in differentiating Meloidogyne incognita eggs may lead to irregular development, lethality and by extension, represent a novel target for parasite control On silencing the nuclear RNase III enzyme drosha, a critical effector of miRNA maturation in animals, we found a significant inhibition of normal development and hatching in short interfering (sORNA-soaked M incognita eggs Developing juveniles presented with highly irregular tissue patterning within the egg, and we found that unlike our previous gene silencing efforts focused on FMRFamide (Phe-Met-Arg-Phe-NH2)-like peptides (FLPs), there was no observable phenotypic recovery following removal of the environmental siRNA. Aberrant phenotypes were exacerbated over time, and drosha knockdown proved embryonically lethal Subsequently, we identified and silenced the drosha cofactor pasha, revealing a comparable inhibition of normal embryonic development within the eggs to that of drosha-silenced eggs, eventually leading to embryonic lethality To further probe the link between normal embryonic development and the M. incognita RNA interference (RNAi) pathway, we attempted to examine the impact of silencing the cytosolic RNase III enzyme dicer. Unexpectedly, we found a substantial up-regulation of dicer transcript abundance, which did not impact on egg differentiation or hatching rates. Silencing of the individual transcripts in hatched J2s was significantly less successful and resulted in temporary phenotypic aberration of the J2s. which recovered within 24 h to normal movement and posture on washing out the siRNA. Soaking the J2s in dicer siRNA resulted in a modest decrease in dicer transcript abundance which had no observable impact on phenotype or behaviour within 48 h of initial exposure to siRNA. We propose that drosha, pasha and their ancillary factors may represent excellent targets for novel nematicides and/or in planta controls aimed at M incognita, and potentially other parasitic nematodes, through disruption of miRNA-directed developmental pathways. In addition, we have identified a putative Mi-en-I transcript which encodes an RNAi-inhibiting siRNA exonuclease We observe a marked up-regulation of MI-en-I transcript abundance in response to exogenously introduced siRNA, and reason that this may impact on the interpretation of RN/NI-based reverse genetic screens in plant parasitic nematodes. (C) 2010 Australian Society for Parasitology Inc. Published by Elsevier Ltd. All rights reserved.

Targeted 3' processing of antisense transcripts triggers Arabidopsis FLC chromatin silencing.

Noncoding RNA is emerging as an important regulator of gene expression in many organisms. We are characterizing RNA-mediated chromatin silencing of the Arabidopsis major floral repressor gene, FLC. Through suppressor mutagenesis, we identify a requirement for CstF64 and CstF77, two conserved RNA 3'-end-processing factors, in FLC silencing. However, FLC sense transcript 3' processing is not affected in the mutants. Instead, CstF64 and CstF77 are required for 3' processing of FLC antisense transcripts. A specific RNA-binding protein directs their activity to a proximal antisense polyadenylation site. This targeted processing triggers localized histone demethylase activity and results in reduced FLC sense transcription. Targeted 3' processing of antisense transcripts may be a common mechanism triggering transcriptional silencing of the corresponding sense gene.

Cold-induced silencing by long antisense transcripts of an Arabidopsis Polycomb target.

Transcription in eukaryotic genomes generates an extensive array of non-protein-coding RNA, the functional significance of which is mostly unknown. We are investigating the link between non-coding RNA and chromatin regulation through analysis of FLC - a regulator of flowering time in Arabidopsis and a target of several chromatin pathways. Here we use an unbiased strategy to characterize non-coding transcripts of FLC and show that sense/antisense transcript levels correlate in a range of mutants and treatments, but change independently in cold-treated plants. Prolonged cold epigenetically silences FLC in a Polycomb-mediated process called vernalization. Our data indicate that upregulation of long non-coding antisense transcripts covering the entire FLC locus may be part of the cold-sensing mechanism. Induction of these antisense transcripts occurs earlier than, and is independent of, other vernalization markers and coincides with a reduction in sense transcription. We show that addition of the FLC antisense promoter sequences to a reporter gene is sufficient to confer cold-induced silencing of the reporter. Our data indicate that cold-induced FLC antisense transcripts have an early role in the epigenetic silencing of FLC, acting to silence FLC transcription transiently. Recruitment of the Polycomb machinery then confers the epigenetic memory. Antisense transcription events originating from 3' ends of genes might be a general mechanism to regulate the corresponding sense transcription in a condition/stage-dependent manner.

Guidelines for the measurement of BCR-ABL1 transcripts in chronic myeloid leukaemia

Molecular testing for the BCR-ABL1 fusion gene by real time quantitative polymerase chain reaction (RT-qPCR) is the most sensitive routine approach for monitoring the response to therapy of patients with chronic myeloid leukaemia. In the context of tyrosine kinase inhibitor (TKI) therapy, the technique is most appropriate for patients who have achieved complete cytogenetic remission and can be used to define specific therapeutic milestones. To achieve this effectively, standardization of the laboratory procedures and the interpretation of results are essential. We present here consensus best practice guidelines for RT-qPCR testing, data interpretation and reporting that have been drawn up and agreed by a consortium of 21 testing laboratories in the United Kingdom and Ireland in accordance with the procedures of the UK Clinical Molecular Genetics Society.

Resolving candidate genes of mouse skeletal muscle QTL via RNA-Seq and expression network analyses

BACKGROUND:

We have recently identified a number of Quantitative Trait Loci (QTL) contributing to the 2-fold muscle weight difference between the LG/J and SM/J mouse strains and refined their confidence intervals. To facilitate nomination of the candidate genes responsible for these differences we examined the transcriptome of the tibialis anterior (TA) muscle of each strain by RNA-Seq.

13,726 genes were expressed in mouse skeletal muscle. Intersection of a set of 1061 differentially expressed transcripts with a mouse muscle Bayesian Network identified a coherent set of differentially expressed genes that we term the LG/J and SM/J Regulatory Network (LSRN). The integration of the QTL, transcriptome and the network analyses identified eight key drivers of the LSRN (Kdr, Plbd1, Mgp, Fah, Prss23, 2310014F06Rik, Grtp1, Stk10) residing within five QTL regions, which were either polymorphic or differentially expressed between the two strains and are strong candidates for quantitative trait genes (QTGs) underlying muscle mass. The insight gained from network analysis including the ability to make testable predictions is illustrated by annotating the LSRN with knowledge-based signatures and showing that the SM/J state of the network corresponds to a more oxidative state. We validated this prediction by NADH tetrazolium reductase staining in the TA muscle revealing higher oxidative potential of the SM/J compared to the LG/J strain (p<0.03).

Thus, integration of fine resolution QTL mapping, RNA-Seq transcriptome information and mouse muscle Bayesian Network analysis provides a novel and unbiased strategy for nomination of muscle QTGs.

The Arabidopsis RNA-binding protein FCA requires a lysine-specific demethylase 1 homolog to downregulate FLC.

A repressor of the transition to flowering in Arabidopsis is the MADS box protein FLOWERING LOCUS C (FLC). FCA, an RNA-binding protein, and FY, a homolog of the yeast RNA 3' processing factor Pfs2p, downregulate FLC expression and therefore promote flowering. FCA/FY physically interact and alter polyadenylation/3' processing to negatively autoregulate FCA. Here, we show that FCA requires FLOWERING LOCUS D (FLD), a homolog of the human lysine-specific demethylase 1 (LSD1) for FLC downregulation. FCA also partially depends on DICER-LIKE 3, involved in chromatin silencing. fca mutations increased levels of unspliced sense FLC transcript, altered processing of antisense FLC transcripts, and increased H3K4 dimethylation in the central region of FLC. These data support a close association of FCA and FLD in mediating H3K4 demethylation and thus transcriptional silencing of FLC and reveal roles for antisense RNA processing and DCL3 function in this regulation.

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.

Exploring TERRA (TElomeric Repeat-containing RNA) Expression and Regulation During Cell Growth in Saccharomyces cerevisiae

Please find the referenced videos attached

Temporal regulation of expression of immediate early and second phase transcripts by endothelin-1 in cardiomyocytes

Background: Endothelin-1 stimulates Gq protein-coupled receptors to promote proliferation in dividing cells or hypertrophy in terminally differentiated cardiomyocytes. In cardiomyocytes, endothelin-1 rapidly (within minutes) stimulates protein kinase signaling, including extracellular-signal regulated kinases 1/2 (ERK1/2; though not ERK5), with phenotypic/physiological changes developing from approximately 12 h. Hypertrophy is associated with changes in mRNA/protein expression, presumably consequent to protein kinase signaling, but the connections between early, transient signaling events and developed hypertrophy are unknown. Results: Using microarrays, we defined the early transcriptional responses of neonatal rat cardiomyocytes to endothelin-1 over 4 h, differentiating between immediate early gene (IEG) and second phase RNAs with cycloheximide. IEGs exhibited differential temporal and transient regulation, with expression of second phase RNAs within 1 h. Of transcripts upregulated at 30 minutes encoding established proteins, 28 were inhibited >50% by U0126 (which inhibits ERK1/2/5 signaling), with 9 inhibited 25-50%. Expression of only four transcripts was not inhibited. At 1 h, most RNAs (approximately 67%) were equally changed in total and polysomal RNA with approximately 17% of transcripts increased to a greater extent in polysomes. Thus, changes in expression of most protein-coding RNAs should be reflected in protein synthesis. However, approximately 16% of transcripts were essentially excluded from the polysomes, including some protein-coding mRNAs, presumably inefficiently translated. Conclusion: The phasic, temporal regulation of early transcriptional responses induced by endothelin-1 in cardiomyocytes indicates that, even in terminally differentiated cells, signals are propagated beyond the primary signaling pathways through transcriptional networks leading to phenotypic changes (that is, hypertrophy). Furthermore, ERK1/2 signaling plays a major role in this response.

Cell walls of developing wheat starchy endosperm: comparison of composition and RNA-Seq transcriptome

The transcriptome of the developing starchy endosperm of hexaploid wheat (Triticum aestivum) was determined using RNA-Seq isolated at five stages during grain fill. This resource represents an excellent way to identify candidate genes responsible for the starchy endosperm cell wall, which is dominated by arabinoxylan (AX), accounting for 70% of the cell wall polysaccharides, with 20% (1,3; 1,4)-beta-D-glucan, 7% glucomannan, and 4% cellulose. A complete inventory of transcripts of 124 glycosyltransferase (GT) and 72 glycosylhydrolase (GH) genes associated with cell walls is presented. The most highly expressed GT transcript (excluding those known to be involved in starch synthesis) was a GT47 family transcript similar to Arabidopsis (Arabidopsis thaliana) IRX10 involved in xylan extension, and the second most abundant was a GT61. Profiles for GT43 IRX9 and IRX14 putative orthologs were consistent with roles in AX synthesis. Low abundances were found for transcripts from genes in the acyl-coA transferase BAHD family, for which a role in AX feruloylation has been postulated. The relative expression of these was much greater in whole grain compared with starchy endosperm, correlating with the levels of bound ferulate. Transcripts associated with callose (GSL), cellulose (CESA), pectin (GAUT), and glucomannan (CSLA) synthesis were also abundant in starchy endosperm, while the corresponding cell wall polysaccharides were confirmed as low abundance (glucomannan and callose) or undetectable (pectin) in these samples. Abundant transcripts from GH families associated with the hydrolysis of these polysaccharides were also present, suggesting that they may be rapidly turned over. Abundant transcripts in the GT31 family may be responsible for the addition of Gal residues to arabinogalactan peptide.