Genome-wide kinetic profiling of pre-mRNA 3' end cleavage.

Cleavage and polyadenylation is necessary for the formation of mature mRNA molecules. The rate at which this process occurs can determine the temporal availability of mRNA for subsequent function throughout the cell and is likely tightly regulated. Despite advances in high-throughput approaches for global kinetic profiling of RNA maturation, genome-wide 3' end cleavage rates have never been measured. Here, we describe a novel approach to estimate the rates of cleavage, using metabolic labeling of nascent RNA, high-throughput sequencing, and mathematical modeling. Using in silico simulations of nascent RNA-seq data, we show that our approach can accurately and precisely estimate cleavage half-lives for both constitutive and alternative sites. We find that 3' end cleavage is fast on average, with half-lives under a minute, but highly variable across individual sites. Rapid cleavage is promoted by the presence of canonical sequence elements and an increased density of polyadenylation signals near a cleavage site. Finally, we find that cleavage rates are associated with the localization of RNA polymerase II at the end of a gene, and faster cleavage leads to quicker degradation of downstream readthrough RNA. Our findings shed light on the features important for efficient 3' end cleavage and the regulation of transcription termination.

mRNA initiation and termination are spatially coordinated.

The expression of a precise mRNA transcriptome is crucial for establishing cell identity and function, with dozens of alternative isoforms produced for a single gene sequence. The regulation of mRNA isoform usage occurs by the coordination of co-transcriptional mRNA processing mechanisms across a gene. Decisions involved in mRNA initiation and termination underlie the largest extent of mRNA isoform diversity, but little is known about any relationships between decisions at both ends of mRNA molecules. Here, we systematically profile the joint usage of mRNA transcription start sites (TSSs) and polyadenylation sites (PASs) across tissues and species. Using both short and long read RNA-seq data, we observe that mRNAs preferentially using upstream TSSs also tend to use upstream PASs, and congruently, the usage of downstream sites is similarly paired. This observation suggests that mRNA 5' end choice may directly influence mRNA 3' ends. Our results suggest a novel "Positional Initiation-Termination Axis" (PITA), in which the usage of alternative terminal sites are coupled based on the order in which they appear in the genome. PITA isoforms are more likely to encode alternative protein domains and use conserved sites. PITA is strongly associated with the length of genomic features, such that PITA is enriched in longer genes with more area devoted to regions that regulate alternative 5' or 3' ends. Strikingly, we found that PITA genes are more likely than non-PITA genes to have multiple, overlapping chromatin structural domains related to pairing of ordinally coupled start and end sites. In turn, PITA coupling is also associated with fast RNA Polymerase II (RNAPII) trafficking across these long gene regions. Our findings indicate that a combination of spatial and kinetic mechanisms couple transcription initiation and mRNA 3' end decisions based on ordinal position to define the expression mRNA isoforms.

Challenges in identifying mRNA transcript starts and ends from long-read sequencing data.

Long-read sequencing (LRS) technologies have the potential to revolutionize scientific discoveries in RNA biology, especially by enabling the comprehensive identification and quantification of full length mRNA isoforms. However, inherently high error rates make the analysis of long-read sequencing data challenging. While these error rates have been characterized for sequence and splice site identification, it is still unclear how accurately LRS reads represent transcript start and end sites. Here, we systematically assess the variability and accuracy of mRNA terminal ends identified by LRS reads across multiple sequencing platforms. We find substantial inconsistencies in both the start and end coordinates of LRS reads spanning a gene, such that LRS reads often fail to accurately recapitulate annotated or empirically derived terminal ends of mRNA molecules. To address this challenge, we introduce an approach to condition reads based on empirically derived terminal ends and identified a subset of reads that are more likely to represent full-length transcripts. Our approach can improve transcriptome analyses by enhancing the fidelity of transcript terminal end identification, but may result in lower power to quantify genes or discover novel isoforms. Thus, it is necessary to be cautious when selecting sequencing approaches and/or interpreting data from long-read RNA sequencing.

G-rich motifs within phosphorothioate-based antisense oligonucleotides (ASOs) drive activation of FXN expression through indirect effects.

Friedreich's ataxia is an incurable disease caused by frataxin (FXN) protein deficiency, which is mostly induced by GAA repeat expansion in intron 1 of the FXN gene. Here, we identified antisense oligonucleotides (ASOs), complementary to two regions within the first intron of FXN pre-mRNA, which could increase FXN mRNA by ∼2-fold in patient fibroblasts. The increase in FXN mRNA was confirmed by the identification of multiple overlapping FXN-activating ASOs at each region, two independent RNA quantification assays, and normalization by multiple housekeeping genes. Experiments on cells with the ASO-binding sites deleted indicate that the ASO-induced FXN activation was driven by indirect effects. RNA sequencing analyses showed that the two ASOs induced similar transcriptome-wide changes, which did not resemble the transcriptome of wild-type cells. This RNA-seq analysis did not identify directly base-paired off-target genes shared across ASOs. Mismatch studies identified two guanosine-rich motifs (CCGG and G4) within the ASOs that were required for FXN activation. The phosphorodiamidate morpholino oligomer analogs of our ASOs did not activate FXN, pointing to a PS-backbone-mediated effect. Our study demonstrates the importance of multiple, detailed control experiments and target validation in oligonucleotide studies employing novel mechanisms such as gene activation.

High-throughput analysis of ANRIL circRNA isoforms in human pancreatic islets.

The antisense non-coding RNA in the INK locus (ANRIL) is a hotspot for genetic variants associated with cardiometabolic disease. We recently found increased ANRIL abundance in human pancreatic islets from donors with certain Type II Diabetes (T2D) risk-SNPs, including a T2D risk-SNP located within ANRIL exon 2 associated with beta cell proliferation. Recent studies have found that expression of circular species of ANRIL is linked to the regulation of cardiovascular phenotypes. Less is known about how the abundance of circular ANRIL may influence T2D phenotypes. Herein, we sequence circular RNA in pancreatic islets to characterize circular isoforms of ANRIL. We identify several consistently expressed circular ANRIL isoforms whose expression is correlated across dozens of individuals and characterize ANRIL splice sites that are commonly involved in back-splicing. We find that samples with the T2D risk allele in ANRIL exon 2 had higher ratios of circular to linear ANRIL compared to protective-allele carriers, and that higher circular:linear ANRIL was associated with decreased beta cell proliferation. Our study points to a combined involvement of both linear and circular ANRIL species in T2D phenotypes and opens the door for future studies of the molecular mechanisms by which ANRIL impacts cellular function in pancreatic islets.

Widespread occurrence of hybrid internal-terminal exons in human transcriptomes.

Messenger RNA isoform differences are predominantly driven by alternative first, internal, and last exons. Despite the importance of classifying exons to understand isoform structure, few tools examine isoform-specific exon usage. We recently observed that alternative transcription start sites often arise near internal exons, often creating "hybrid" first/internal exons. To systematically detect hybrid exons, we built the hybrid-internal-terminal (HIT) pipeline to classify exons depending on their isoform-specific usage. On the basis of splice junction reads in RNA sequencing data and probabilistic modeling, the HIT index identified thousands of previously misclassified hybrid first-internal and internal-last exons. Hybrid exons are enriched in long genes and genes involved in RNA splicing and have longer flanking introns and strong splice sites. Their usage varies considerably across human tissues. By developing the first method to classify exons according to isoform contexts, our findings document the occurrence of hybrid exons, a common quirk of the human transcriptome.

HERC1 Regulates Breast Cancer Cells Migration and Invasion.

Tumor cell migration and invasion into adjacent tissues is one of the hallmarks of cancer and the first step towards secondary tumors formation, which represents the leading cause of cancer-related deaths. This process is considered an unmet clinical need in the treatment of this disease, particularly in breast cancers characterized by high aggressiveness and metastatic potential. To identify and characterize genes with novel functions as regulators of tumor cell migration and invasion, we performed a genetic loss-of-function screen using a shRNA library directed against the Ubiquitin Proteasome System (UPS) in a highly invasive breast cancer derived cell line. Among the candidates, we validated HERC1 as a gene regulating cell migration and invasion. Furthermore, using animal models, our results indicate that HERC1 silencing affects primary tumor growth and lung colonization. Finally, we conducted an in silico analysis using publicly available protein expression data and observed an inverse correlation between HERC1 expression levels and breast cancer patients' overall survival. Altogether, our findings demonstrate that HERC1 might represent a novel therapeutic target for the development or improvement of breast cancer treatment.

USP19 modulates cancer cell migration and invasion and acts as a novel prognostic marker in patients with early breast cancer.

Tumor cell dissemination in cancer patients is associated with a significant reduction in their survival and quality of life. The ubiquitination pathway plays a fundamental role in the maintenance of protein homeostasis both in normal and stressed conditions and its dysregulation has been associated with malignant transformation and invasive potential of tumor cells, thus highlighting its value as a potential therapeutic target. In order to identify novel molecular targets of tumor cell migration and invasion we performed a genetic screen with an shRNA library against ubiquitination pathway-related genes. To this end, we set up a protocol to specifically enrich positive migration regulator candidates. We identified the deubiquitinase USP19 and demonstrated that its silencing reduces the migratory and invasive potential of highly invasive breast cancer cell lines. We extended our investigation in vivo and confirmed that mice injected with USP19 depleted cells display increased tumor-free survival, as well as a delay in the onset of the tumor formation and a significant reduction in the appearance of metastatic foci, indicating that tumor cell invasion and dissemination is impaired. In contrast, overexpression of USP19 increased cell invasiveness both in vitro and in vivo, further validating our findings. More importantly, we demonstrated that USP19 catalytic activity is important for the control of tumor cell migration and invasion, and that its molecular mechanism of action involves LRP6, a Wnt co-receptor. Finally, we showed that USP19 overexpression is a surrogate prognostic marker of distant relapse in patients with early breast cancer. Altogether, these findings demonstrate that USP19 might represent a novel therapeutic target in breast cancer.