Intron Retention and TE Exonization Events in ZRANB2

Transcription start site signal profiling improves transposable element RNA expression analysis at locus-level

The transcriptional activity of Transposable Elements (TEs) has been involved in numerous pathological processes, including neurodegenerative diseases such as amyotrophic lateral sclerosis and frontotemporal lobar degeneration. The TE expression analysis from short-read sequencing technologies is, however, challenging due to the multitude of similar sequences derived from singular TEs subfamilies and the exaptation of TEs within longer coding or non-coding RNAs. Specialised tools have been developed to quantify the expression of TEs that either relies on probabilistic re-distribution of multimapper count fractions or allow for discarding multimappers altogether. Until now, the benchmarking across those tools was largely limited to aggregated expression estimates over whole TEs subfamilies. Here, we compared the performance of recently published tools (SQuIRE, TElocal, SalmonTE) with simplistic quantification strategies (featureCounts in unique, fraction and random modes) at the individual loci level. Using simulated datasets, we examined the false discovery rate and the primary driver of those false positive hits in the optimal quantification strategy. Our findings suggest a high false discovery number that exceeds the total number of correctly recovered active loci for all the quantification strategies, including the best performing tool TElocal. As a remedy, filtering based on the minimum number of read counts or baseMean expression improves the F1 score and decreases the number of false positives. Finally, we demonstrate that additional profiling of Transcription Start Site mapping statistics (using a k-means clustering approach) significantly improves the performance of TElocal while reporting a reliable set of detected and differentially expressed TEs in human simulated RNA-seq data.

The Dual Role of miR-186 in Cancers: Oncomir Battling With Tumor Suppressor miRNA

MicroRNAs (miRNAs) are a class of small non-coding RNAs which regulate gene expression at post-transcriptional level. Alterations of miR-186 expression were demonstrated in numerous cancers, shown to play a vital role in oncogenesis, invasion, metastasis, apoptosis, and drug resistance. MiR-186 was documented as a tumor suppressor miRNA in the majority of studies, while conflicting reports verified miR-186 as an oncomir. The contradictory role in cancers may impede the application of miR-186, as well as other dual-functional miRNAs, as a diagnostic and therapeutic target. This review emphasizes the alterations and functions of miR-186 in cancers and discusses the mechanisms behind the contradictory findings. Among these, target abundance and dose-dependent effects of miR-186 are highlighted. The paper aims to review the challenges involved in developing diagnostic and therapeutic strategies for cancer treatment based on dual-functional miRNAs.

Gain of a New Exon by a Lineage-Specific Alu Element-Integration Event in the BCS1L Gene during Primate Evolution

BCS1L gene encodes mitochondrial protein and is a member of conserved AAA protein family. This gene is involved in the incorporation of Rieske FeS and Qcr10p into complex III of respiratory chain. In our previous study, AluYRa2-derived alternative transcript in rhesus monkey genome was identified. However, this transcript has not been reported in human genome. In present study, we conducted evolutionary analysis of AluYRa2-exonized transcript with various primate genomic DNAs and cDNAs from humans, rhesus monkeys, and crab-eating monkeys. Remarkably, our results show that AluYRa2 element has only been integrated into genomes of Macaca species. This Macaca lineage-specific integration of AluYRa2 element led to exonization event in the first intron region of BCS1L gene by producing a conserved 3' splice site. Intriguingly, in rhesus and crab-eating monkeys, more diverse transcript variants by alternative splicing (AS) events, including exon skipping and different 5' splice sites from humans, were identified. Alignment of amino acid sequences revealed that AluYRa2-exonized transcript has short N-terminal peptides. Therefore, AS events play a major role in the generation of various transcripts and proteins during primate evolution. In particular, lineage-specific integration of Alu elements and species-specific Alu-derived exonization events could be important sources of gene diversification in primates.

Alternative splicing acting as a bridge in evolution

BACKGROUND

Alternative splicing (AS) regulates diverse cellular and developmental functions through alternative protein structures of different isoforms. Alternative exons dominate AS in vertebrates; however, very little is known about the extent and function of AS in lower eukaryotes. To understand the role of introns in gene evolution, we examined AS from a green algal and five fungal genomes using a novel EST-based gene-modeling algorithm (COMBEST).

METHODS

AS from each genome was classified with COMBEST that maps EST sequences to genomes to build gene models. Various aspects of AS were analyzed through statistical methods. The interplay of intron 3n length, phase, coding property, and intron retention (RI) were examined with Chi-square testing.

RESULTS

With 3 to 834 times EST coverage, we identified up to 73% of AS in intron-containing genes and found preponderance of RI among 11 types of AS. The number of exons, expression level, and maximum intron length correlated with number of AS per gene (NAG), and intron-rich genes suppressed AS. Genes with AS were more ancient, and AS was conserved among fungal genomes. Among stopless introns, non-retained introns (NRI) avoided, but major RI preferred 3n length. In contrast, stop-containing introns showed uniform distribution among 3n, 3n+1, and 3n+2 lengths. We found a clue to the intron phase enigma: it was the coding function of introns involved in AS that dictates the intron phase bias.

CONCLUSIONS

Majority of AS is non-functional, and the extent of AS is suppressed for intron-rich genes. RI through 3n length, stop codon, and phase bias bridges the transition from functionless to functional alternative isoforms.

miR-186 inhibits muscle cell differentiation through myogenin regulation

The complex process of skeletal muscle differentiation is organized by the myogenic regulatory factors (MRFs), Myf5, MyoD, Myf6, and myogenin, where myogenin plays a critical role in the regulation of the final stage of muscle differentiation. In an effort to investigate the role microRNAs (miRNAs) play in regulating myogenin, a bioinformatics approach was used and six miRNAs (miR-182, miR-186, miR-135, miR-491, miR-329, and miR-96) were predicted to bind the myogenin 3'-untranslated region (UTR). However, luciferase assays showed only miR-186 inhibited translation and 3'-UTR mutagenesis analysis confirmed this interaction was specific. Interestingly, the expression of miR-186 mirrored that of its host gene, ZRANB2, during development. Functional studies demonstrated that miR-186 overexpression inhibited the differentiation of C2C12 and primary muscle cells. Our findings therefore identify miR-186 as a novel regulator of myogenic differentiation.

RNA sequencing of the exercise transcriptome in equine athletes

The horse is an optimal model organism for studying the genomic response to exercise-induced stress, due to its natural aptitude for athletic performance and the relative homogeneity of its genetic and environmental backgrounds. Here, we applied RNA-sequencing analysis through the use of SOLiD technology in an experimental framework centered on exercise-induced stress during endurance races in equine athletes. We monitored the transcriptional landscape by comparing gene expression levels between animals at rest and after competition. Overall, we observed a shift from coding to non-coding regions, suggesting that the stress response involves the differential expression of not annotated regions. Notably, we observed significant post-race increases of reads that correspond to repeats, especially the intergenic and intronic L1 and L2 transposable elements. We also observed increased expression of the antisense strands compared to the sense strands in intronic and regulatory regions (1 kb up- and downstream) of the genes, suggesting that antisense transcription could be one of the main mechanisms for transposon regulation in the horse under stress conditions. We identified a large number of transcripts corresponding to intergenic and intronic regions putatively associated with new transcriptional elements. Gene expression and pathway analysis allowed us to identify several biological processes and molecular functions that may be involved with exercise-induced stress. Ontology clustering reflected mechanisms that are already known to be stress activated (e.g., chemokine-type cytokines, Toll-like receptors, and kinases), as well as "nucleic acid binding" and "signal transduction activity" functions. There was also a general and transient decrease in the global rates of protein synthesis, which would be expected after strenuous global stress. In sum, our network analysis points toward the involvement of specific gene clusters in equine exercise-induced stress, including those involved in inflammation, cell signaling, and immune interactions.