HALO--a Java framework for precise transcript half-life determination

Transcriptomic differences between fibrotic and non-fibrotic testicular tissue reveal possible key players in Klinefelter syndrome-related testicular fibrosis

Klinefelter syndrome (KS; 47,XXY) affects 1-2 in 1000 males. Most men with KS suffer from an early germ cell loss and testicular fibrosis from puberty onwards. Mechanisms responsible for these processes remain unknown. Previous genomics studies on testis tissue from men with KS focused on germ cell loss, while a transcriptomic analysis focused on testicular fibrosis has not yet been performed. This study aimed to identify factors involved in the fibrotic remodelling of KS testes by analysing the transcriptome of fibrotic and non-fibrotic testicular tissue. RNA sequencing was performed to compare the genes expressed in testicular samples with (KS and testis atrophy) and without (Sertoli cell-only syndrome and fertile controls) fibrosis (n = 5, each). Additionally, differentially expressed genes (DEGs) between KS and testis atrophy samples were studied to reveal KS-specific fibrotic genes. DEGs were considered significant when p < 0.01 and log2FC > 2. Next, downstream analyses (GO and KEGG) were performed. Lastly, RNA in situ hybridization was performed to validate the results. The first analysis (fibrotic vs non-fibrotic) resulted in 734 significant DEGs (167 up- and 567 down-regulated). Genes involved in the extracellular structure organization (e.g. VCAM1) were found up-regulated. KEGG analysis showed an up-regulation of genes involved in the TGF-β pathway. The KS vs testis atrophy analysis resulted in 539 significant DEGs (59 up- and 480 down-regulated). Chronic inflammatory response genes were found up-regulated. The overlap of X-linked DEGs from the two analyses revealed three genes: matrix-remodelling associated 5 (MXRA5), doublecortin (DCX) and variable charge X-Linked 3B (VCX3B). RNA in situ hybridization showed an overexpression of VCAM1, MXRA5 and DCX within the fibrotic group compared with the non-fibrotic group. To summarize, this study revealed DEGs between fibrotic and non-fibrotic testis tissue, including VCAM1. In addition, X-linked fibrotic genes were revealed, e.g. MXRA5, DCX and VCX3B. Their potential role in KS-related testicular fibrosis needs further study.

Circular RNAs are long-lived and display only minimal early alterations in response to a growth factor

Circular RNAs (circRNAs) are widespread circles of non-coding RNAs with largely unknown function. Because stimulation of mammary cells with the epidermal growth factor (EGF) leads to dynamic changes in the abundance of coding and non-coding RNA molecules, and culminates in the acquisition of a robust migratory phenotype, this cellular model might disclose functions of circRNAs. Here we show that circRNAs of EGF-stimulated mammary cells are stably expressed, while mRNAs and microRNAs change within minutes. In general, the circRNAs we detected are relatively long-lived and weakly expressed. Interestingly, they are almost ubiquitously co-expressed with the corresponding linear transcripts, and the respective, shared promoter regions are more active compared to genes producing linear isoforms with no detectable circRNAs. These findings imply that altered abundance of circRNAs, unlike changes in the levels of other RNAs, might not play critical roles in signaling cascades and downstream transcriptional networks that rapidly commit cells to specific outcomes.

Metabolic labeling of newly transcribed RNA for high resolution gene expression profiling of RNA synthesis, processing and decay in cell culture

The development of whole-transcriptome microarrays and next-generation sequencing has revolutionized our understanding of the complexity of cellular gene expression. Along with a better understanding of the involved molecular mechanisms, precise measurements of the underlying kinetics have become increasingly important. Here, these powerful methodologies face major limitations due to intrinsic properties of the template samples they study, i.e. total cellular RNA. In many cases changes in total cellular RNA occur either too slowly or too quickly to represent the underlying molecular events and their kinetics with sufficient resolution. In addition, the contribution of alterations in RNA synthesis, processing, and decay are not readily differentiated. We recently developed high-resolution gene expression profiling to overcome these limitations. Our approach is based on metabolic labeling of newly transcribed RNA with 4-thiouridine (thus also referred to as 4sU-tagging) followed by rigorous purification of newly transcribed RNA using thiol-specific biotinylation and streptavidin-coated magnetic beads. It is applicable to a broad range of organisms including vertebrates, Drosophila, and yeast. We successfully applied 4sU-tagging to study real-time kinetics of transcription factor activities, provide precise measurements of RNA half-lives, and obtain novel insights into the kinetics of RNA processing. Finally, computational modeling can be employed to generate an integrated, comprehensive analysis of the underlying molecular mechanisms.

Global analysis of mRNA half-lives and de novo transcription in a dinoflagellate, Karenia brevis

Dinoflagellates possess many physiological processes that appear to be under post-transcriptional control. However, the extent to which their genes are regulated post-transcriptionally remains unresolved. To gain insight into the roles of differential mRNA stability and de novo transcription in dinoflagellates, we biosynthetically labeled RNA with 4-thiouracil to isolate newly transcribed and pre-existing RNA pools in Karenia brevis. These isolated fractions were then used for analysis of global mRNA stability and de novo transcription by hybridization to a K. brevis microarray. Global K. brevis mRNA half-lives were calculated from the ratio of newly transcribed to pre-existing RNA for 7086 array features using the online software HALO (Half-life Organizer). Overall, mRNA half-lives were substantially longer than reported in other organisms studied at the global level, ranging from 42 minutes to greater than 144 h, with a median of 33 hours. Consistent with well-documented trends observed in other organisms, housekeeping processes, including energy metabolism and transport, were significantly enriched in the most highly stable messages. Shorter-lived transcripts included a higher proportion of transcriptional regulation, stress response, and other response/regulatory processes. One such family of proteins involved in post-transcriptional regulation in chloroplasts and mitochondria, the pentatricopeptide repeat (PPR) proteins, had dramatically shorter half-lives when compared to the arrayed transcriptome. As transcript abundances for PPR proteins were previously observed to rapidly increase in response to nutrient addition, we queried the newly synthesized RNA pools at 1 and 4 h following nitrate addition to N-depleted cultures. Transcriptome-wide there was little evidence of increases in the rate of de novo transcription during the first 4 h, relative to that in N-depleted cells, and no evidence for increased PPR protein transcription. These results lend support to the growing consensus of post-transcriptional control of gene expression in dinoflagellates.

Measurement of genome-wide RNA synthesis and decay rates with Dynamic Transcriptome Analysis (DTA)

Standard transcriptomics measures total cellular RNA levels. Our understanding of gene regulation would be greatly improved if we could measure RNA synthesis and decay rates on a genome-wide level. To that end, the Dynamic Transcriptome Analysis (DTA) method has been developed. DTA combines metabolic RNA labeling with standard transcriptomics to measure RNA synthesis and decay rates in a precise and non-perturbing manner. Here, we present the open source R/Bioconductor software package DTA. It implements all required bioinformatics steps that allow the accurate absolute quantification and comparison of RNA turnover.