Liver kinase B1 (LKB1) regulates the epigenetic landscape of mouse pancreatic beta cells.

Liver kinase B1 (LKB1/STK11) is an important regulator of pancreatic ß-cell identity and function. Elimination of Lkb1 from the ß-cell results in improved glucose-stimulated insulin secretion and is accompanied by profound changes in gene expression, including the upregulation of several neuronal genes. The mechanisms through which LKB1 controls gene expression are, at present, poorly understood. Here, we explore the impact of ß cell-selective deletion of Lkb1 on chromatin accessibility in mouse pancreatic islets. To characterize the role of LKB1 in the regulation of gene expression at the transcriptional level, we combine these data with a map of islet active transcription start sites and histone marks. We demonstrate that LKB1 elimination from ß-cells results in widespread changes in chromatin accessibility, correlating with changes in transcript levels. Changes occurred in hundreds of promoter and enhancer regions, many of which were close to neuronal genes. We reveal that dysregulated enhancers are enriched in binding motifs for transcription factors (TFs) important for ß-cell identity, such as FOXA, MAFA or RFX6, and we identify microRNAs (miRNAs) that are regulated by LKB1 at the transcriptional level. Overall, our study provides important new insights into the epigenetic mechanisms by which LKB1 regulates ß-cell identity and function.

Liver kinase B1 (LKB1) regulates the epigenetic landscape of mouse pancreatic beta cells.

Liver kinase B1 (LKB1/STK11) is an important regulator of pancreatic ß-cell identity and function. Elimination of Lkb1 from the ß-cell results in improved glucose-stimulated insulin secretion and is accompanied by profound changes in gene expression, including the upregulation of several neuronal genes. The mechanisms through which LKB1 controls gene expression are, at present, poorly understood. Here, we explore the impact of ß cell- selective deletion of Lkb1 on chromatin accessibility in mouse pancreatic islets. To characterize the role of LKB1 in the regulation of gene expression at the transcriptional level, we combine these data with a map of islet active transcription start sites and histone marks. We demonstrate that LKB1 elimination from ß-cells results in widespread changes in chromatin accessibility, correlating with changes in transcript levels. Changes occurred in hundreds of promoter and enhancer regions, many of which were close to neuronal genes. We reveal that dysregulated enhancers are enriched in binding motifs for transcription factors important for ß-cell identity, such as FOXA, MAFA or RFX6 and we identify microRNAs (miRNAs) that are regulated by LKB1 at the transcriptional level. Overall, our study provides important new insights into the epigenetic mechanisms by which LKB1 regulates ß-cell identity and function.

Liver RBFOX2 regulates cholesterol homeostasis via Scarb1 alternative splicing in mice.

RNA alternative splicing (AS) expands the regulatory potential of eukaryotic genomes. The mechanisms regulating liver-specific AS profiles and their contribution to liver function are poorly understood. Here, we identify a key role for the splicing factor RNA-binding Fox protein 2 (RBFOX2) in maintaining cholesterol homeostasis in a lipogenic environment in the liver. Using enhanced individual-nucleotide-resolution ultra-violet cross-linking and immunoprecipitation, we identify physiologically relevant targets of RBFOX2 in mouse liver, including the scavenger receptor class B type I (Scarb1). RBFOX2 function is decreased in the liver in diet-induced obesity, causing a Scarb1 isoform switch and alteration of hepatocyte lipid homeostasis. Our findings demonstrate that specific AS programmes actively maintain liver physiology, and underlie the lipotoxic effects of obesogenic diets when dysregulated. Splice-switching oligonucleotides targeting this network alleviate obesity-induced inflammation in the liver and promote an anti-atherogenic lipoprotein profile in the blood, underscoring the potential of isoform-specific RNA therapeutics for treating metabolism-associated diseases.

Dysregulated RNA polyadenylation contributes to metabolic impairment in non-alcoholic fatty liver disease.

Pre-mRNA processing is an essential mechanism for the generation of mature mRNA and the regulation of gene expression in eukaryotic cells. While defects in pre-mRNA processing have been implicated in a number of diseases their involvement in metabolic pathologies is still unclear. Here, we show that both alternative splicing and alternative polyadenylation, two major steps in pre-mRNA processing, are significantly altered in non-alcoholic fatty liver disease (NAFLD). Moreover, we find that Serine and Arginine Rich Splicing Factor 10 (SRSF10) binding is enriched adjacent to consensus polyadenylation motifs and its expression is significantly decreased in NAFLD, suggesting a role mediating pre-mRNA dysregulation in this condition. Consistently, inactivation of SRSF10 in mouse and human hepatocytes in vitro, and in mouse liver in vivo, was found to dysregulate polyadenylation of key metabolic genes such as peroxisome proliferator-activated receptor alpha (PPARA) and exacerbate diet-induced metabolic dysfunction. Collectively our work implicates dysregulated pre-mRNA polyadenylation in obesity-induced liver disease and uncovers a novel role for SRSF10 in this process.

Author Correction: A systems view of spliceosomal assembly and branchpoints with iCLIP.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

A systems view of spliceosomal assembly and branchpoints with iCLIP.

Studies of spliceosomal interactions are challenging due to their dynamic nature. Here we used spliceosome iCLIP, which immunoprecipitates SmB along with small nuclear ribonucleoprotein particles and auxiliary RNA binding proteins, to map spliceosome engagement with pre-messenger RNAs in human cell lines. This revealed seven peaks of spliceosomal crosslinking around branchpoints (BPs) and splice sites. We identified RNA binding proteins that crosslink to each peak, including known and candidate splicing factors. Moreover, we detected the use of over 40,000 BPs with strong sequence consensus and structural accessibility, which align well to nearby crosslinking peaks. We show how the position and strength of BPs affect the crosslinking patterns of spliceosomal factors, which bind more efficiently upstream of strong or proximally located BPs and downstream of weak or distally located BPs. These insights exemplify spliceosome iCLIP as a broadly applicable method for transcriptomic studies of splicing mechanisms.

The Transcriptome-wide Landscape and Modalities of EJC Binding in Adult Drosophila.

Exon junction complex (EJC) assembles after splicing at specific positions upstream of exon-exon junctions in mRNAs of all higher eukaryotes, affecting major regulatory events. In mammalian cell cytoplasm, EJC is essential for efficient RNA surveillance, while in Drosophila, EJC is essential for localization of oskar mRNA. Here we developed a method for isolation of protein complexes and associated RNA targets (ipaRt) to explore the EJC RNA-binding landscape in a transcriptome-wide manner in adult Drosophila. We find the EJC at canonical positions, preferably on mRNAs from genes comprising multiple splice sites and long introns. Moreover, EJC occupancy is highest at junctions adjacent to strong splice sites, CG-rich hexamers, and RNA structures. Highly occupied mRNAs tend to be maternally localized and derive from genes involved in differentiation or development. These modalities, which have not been reported in mammals, specify EJC assembly on a biologically coherent set of transcripts in Drosophila.

Splicing of long non-coding RNAs primarily depends on polypyrimidine tract and 5' splice-site sequences due to weak interactions with SR proteins.

Many nascent long non-coding RNAs (lncRNAs) undergo the same maturation steps as pre-mRNAs of protein-coding genes (PCGs), but they are often poorly spliced. To identify the underlying mechanisms for this phenomenon, we searched for putative splicing inhibitory sequences using the ncRNA-a2 as a model. Genome-wide analyses of intergenic lncRNAs (lincRNAs) revealed that lincRNA splicing efficiency positively correlates with 5'ss strength while no such correlation was identified for PCGs. In addition, efficiently spliced lincRNAs have higher thymidine content in the polypyrimidine tract (PPT) compared to efficiently spliced PCGs. Using model lincRNAs, we provide experimental evidence that strengthening the 5'ss and increasing the T content in PPT significantly enhances lincRNA splicing. We further showed that lincRNA exons contain less putative binding sites for SR proteins. To map binding of SR proteins to lincRNAs, we performed iCLIP with SRSF2, SRSF5 and SRSF6 and analyzed eCLIP data for SRSF1, SRSF7 and SRSF9. All examined SR proteins bind lincRNA exons to a much lower extent than expression-matched PCGs. We propose that lincRNAs lack the cooperative interaction network that enhances splicing, which renders their splicing outcome more dependent on the optimality of splice sites.

Exon Junction Complex Shapes the Transcriptome by Repressing Recursive Splicing.

Recursive splicing (RS) starts by defining an "RS-exon," which is then spliced to the preceding exon, thus creating a recursive 5' splice site (RS-5ss). Previous studies focused on cryptic RS-exons, and now we find that the exon junction complex (EJC) represses RS of hundreds of annotated, mainly constitutive RS-exons. The core EJC factors, and the peripheral factors PNN and RNPS1, maintain RS-exon inclusion by repressing spliceosomal assembly on RS-5ss. The EJC also blocks 5ss located near exon-exon junctions, thus repressing inclusion of cryptic microexons. The prevalence of annotated RS-exons is high in deuterostomes, while the cryptic RS-exons are more prevalent in Drosophila, where EJC appears less capable of repressing RS. Notably, incomplete repression of RS also contributes to physiological alternative splicing of several human RS-exons. Finally, haploinsufficiency of the EJC factor Magoh in mice is associated with skipping of RS-exons in the brain, with relevance to the microcephaly phenotype and human diseases.

Heteromeric RNP Assembly at LINEs Controls Lineage-Specific RNA Processing.

Long mammalian introns make it challenging for the RNA processing machinery to identify exons accurately. We find that LINE-derived sequences (LINEs) contribute to this selection by recruiting dozens of RNA-binding proteins (RBPs) to introns. This includes MATR3, which promotes binding of PTBP1 to multivalent binding sites within LINEs. Both RBPs repress splicing and 3' end processing within and around LINEs. Notably, repressive RBPs preferentially bind to evolutionarily young LINEs, which are located far from exons. These RBPs insulate the LINEs and the surrounding intronic regions from RNA processing. Upon evolutionary divergence, changes in RNA motifs within LINEs lead to gradual loss of their insulation. Hence, older LINEs are located closer to exons, are a common source of tissue-specific exons, and increasingly bind to RBPs that enhance RNA processing. Thus, LINEs are hubs for the assembly of repressive RBPs and also contribute to the evolution of new, lineage-specific transcripts in mammals. VIDEO ABSTRACT.

'Read-through marking' reveals differential nucleotide composition of read-through and truncated cDNAs in iCLIP.

We established a modified iCLIP protocol, called 'read-through marking', which facilitates the detection of cDNAs that have not been truncated upon encountering the RNA-peptide complex during reverse transcription (read-through cDNAs). A large proportion of these cDNAs would be undesirable in an iCLIP library, as it could affect the resolution of the method. To this end, we added an oligonucleotide to the 5'-end of RNA fragments-a 5'-marker-to mark the read-through cDNAs. By applying this modified iCLIP protocol to PTBP1 and eIF4A3, we found that the start sites of read-through cDNAs are enriched in adenosines, while the remaining cDNAs have a markedly different sequence content at their starts, preferentially containing thymidines. This finding in turn indicates that most of the reads in our iCLIP libraries are a product of truncation with valuable information regarding the proteins' RNA-binding sites. Thus, cDNA start sites confidently identify a protein's RNA-crosslink sites and we can account for the impact of read-through cDNAs by commonly adding a 5'-marker.

Mitochondria maintain controlled activation state of epithelial-resident T lymphocytes.

Epithelial-resident T lymphocytes, such as intraepithelial lymphocytes (IELs) located at the intestinal barrier, can offer swift protection against invading pathogens. Lymphocyte activation is strictly regulated because of its potential harmful nature and metabolic cost, and most lymphocytes are maintained in a quiescent state. However, IELs are kept in a heightened state of activation resembling effector T cells but without cytokine production or clonal proliferation. We show that this controlled activation state correlates with alterations in the IEL mitochondrial membrane, especially the cardiolipin composition. Upon inflammation, the cardiolipin composition is altered to support IEL proliferation and effector function. Furthermore, we show that cardiolipin makeup can particularly restrict swift IEL proliferation and effector functions, reducing microbial containment capability. These findings uncover an alternative mechanism to control cellular activity, special to epithelial-resident T cells, and a novel role for mitochondria, maintaining cells in a metabolically poised state while enabling rapid progression to full functionality.

Data Science Issues in Studying Protein-RNA Interactions with CLIP Technologies.

An interplay of experimental and computational methods is required to achieve a comprehensive understanding of protein-RNA interactions. UV crosslinking and immunoprecipitation (CLIP) identifies endogenous interactions by sequencing RNA fragments that copurify with a selected RNA-binding protein under stringent conditions. Here we focus on approaches for the analysis of the resulting data and appraise the methods for peak calling, visualization, analysis, and computational modeling of protein-RNA binding sites. We advocate that the sensitivity and specificity of data be assessed in combination for computational quality control. Moreover, we demonstrate the value of analyzing sequence motif enrichment in peaks assigned from CLIP data and of visualizing RNA maps, which examine the positional distribution of peaks around regulated landmarks in transcripts. We use these to assess how variations in CLIP data quality and in different peak calling methods affect the insights into regulatory mechanisms. We conclude by discussing future opportunities for the computational analysis of protein-RNA interaction experiments.

High-Resolution RNA Maps Suggest Common Principles of Splicing and Polyadenylation Regulation by TDP-43.

Many RNA-binding proteins (RBPs) regulate both alternative exons and poly(A) site selection. To understand their regulatory principles, we developed expressRNA, a web platform encompassing computational tools for integration of iCLIP and RNA motif analyses with RNA-seq and 3' mRNA sequencing. This reveals at nucleotide resolution the "RNA maps" describing how the RNA binding positions of RBPs relate to their regulatory functions. We use this approach to examine how TDP-43, an RBP involved in several neurodegenerative diseases, binds around its regulated poly(A) sites. Binding close to the poly(A) site generally represses, whereas binding further downstream enhances use of the site, which is similar to TDP-43 binding around regulated exons. Our RNAmotifs2 software also identifies sequence motifs that cluster together with the binding motifs of TDP-43. We conclude that TDP-43 directly regulates diverse types of pre-mRNA processing according to common position-dependent principles.

Insights into the design and interpretation of iCLIP experiments.

BACKGROUND: Ultraviolet (UV) crosslinking and immunoprecipitation (CLIP) identifies the sites on RNAs that are in direct contact with RNA-binding proteins (RBPs). Several variants of CLIP exist, which require different computational approaches for analysis. This variety of approaches can create challenges for a novice user and can hamper insights from multi-study comparisons. Here, we produce data with multiple variants of CLIP and evaluate the data with various computational methods to better understand their suitability. RESULTS: We perform experiments for PTBP1 and eIF4A3 using individual-nucleotide resolution CLIP (iCLIP), employing either UV-C or photoactivatable 4-thiouridine (4SU) combined with UV-A crosslinking and compare the results with published data. As previously noted, the positions of complementary DNA (cDNA)-starts depend on cDNA length in several iCLIP experiments and we now find that this is caused by constrained cDNA-ends, which can result from the sequence and structure constraints of RNA fragmentation. These constraints are overcome when fragmentation by RNase I is efficient and when a broad cDNA size range is obtained. Our study also shows that if RNase does not efficiently cut within the binding sites, the original CLIP method is less capable of identifying the longer binding sites of RBPs. In contrast, we show that a broad size range of cDNAs in iCLIP allows the cDNA-starts to efficiently delineate the complete RNA-binding sites. CONCLUSIONS: We demonstrate the advantage of iCLIP and related methods that can amplify cDNAs that truncate at crosslink sites and we show that computational analyses based on cDNAs-starts are appropriate for such methods.

Cholinergic Surveillance over Hippocampal RNA Metabolism and Alzheimer's-Like Pathology.

The relationship between long-term cholinergic dysfunction and risk of developing dementia is poorly understood. Here we used mice with deletion of the vesicular acetylcholine transporter (VAChT) in the forebrain to model cholinergic abnormalities observed in dementia. Whole-genome RNA sequencing of hippocampal samples revealed that cholinergic failure causes changes in RNA metabolism. Remarkably, key transcripts related to Alzheimer's disease are affected. BACE1, for instance, shows abnormal splicing caused by decreased expression of the splicing regulator hnRNPA2/B1. Resulting BACE1 overexpression leads to increased APP processing and accumulation of soluble Aß1-42. This is accompanied by age-related increases in GSK3 activation, tau hyperphosphorylation, caspase-3 activation, decreased synaptic markers, increased neuronal death, and deteriorating cognition. Pharmacological inhibition of GSK3 hyperactivation reversed deficits in synaptic markers and tau hyperphosphorylation induced by cholinergic dysfunction, indicating a key role for GSK3 in some of these pathological changes. Interestingly, in human brains there was a high correlation between decreased levels of VAChT and hnRNPA2/B1 levels with increased tau hyperphosphorylation. These results suggest that changes in RNA processing caused by cholinergic loss can facilitate Alzheimer's-like pathology in mice, providing a mechanism by which decreased cholinergic tone may increase risk of dementia.

Splicing repression allows the gradual emergence of new Alu-exons in primate evolution.

Alu elements are retrotransposons that frequently form new exons during primate evolution. Here, we assess the interplay of splicing repression by hnRNPC and nonsense-mediated mRNA decay (NMD) in the quality control and evolution of new Alu-exons. We identify 3100 new Alu-exons and show that NMD more efficiently recognises transcripts with Alu-exons compared to other exons with premature termination codons. However, some Alu-exons escape NMD, especially when an adjacent intron is retained, highlighting the importance of concerted repression by splicing and NMD. We show that evolutionary progression of 3' splice sites is coupled with longer repressive uridine tracts. Once the 3' splice site at ancient Alu-exons reaches a stable phase, splicing repression by hnRNPC decreases, but the exons generally remain sensitive to NMD. We conclude that repressive motifs are strongest next to cryptic exons and that gradual weakening of these motifs contributes to the evolutionary emergence of new alternative exons.

Recursive splicing in long vertebrate genes.

It is generally believed that splicing removes introns as single units from precursor messenger RNA transcripts. However, some long Drosophila melanogaster introns contain a cryptic site, known as a recursive splice site (RS-site), that enables a multi-step process of intron removal termed recursive splicing. The extent to which recursive splicing occurs in other species and its mechanistic basis have not been examined. Here we identify highly conserved RS-sites in genes expressed in the mammalian brain that encode proteins functioning in neuronal development. Moreover, the RS-sites are found in some of the longest introns across vertebrates. We find that vertebrate recursive splicing requires initial definition of an 'RS-exon' that follows the RS-site. The RS-exon is then excluded from the dominant mRNA isoform owing to competition with a reconstituted 5' splice site formed at the RS-site after the first splicing step. Conversely, the RS-exon is included when preceded by cryptic promoters or exons that fail to reconstitute an efficient 5' splice site. Most RS-exons contain a premature stop codon such that their inclusion can decrease mRNA stability. Thus, by establishing a binary splicing switch, RS-sites demarcate different mRNA isoforms emerging from long genes by coupling cryptic elements with inclusion of RS-exons.