Search

Valosin-containing protein (VCP) is a hexameric ATPase associated with diverse cellular activities. Genetic mutations in VCP are associated with several forms of muscular and neuronal degeneration, including amyotrophic lateral sclerosis (ALS). Moreover, VCP mediates UV-induced proteolysis of RNA polymerase II (RNAPII), but little is known about the effects of VCP mutations on the transcriptional machinery. Here, we used silica particle-assisted chromatin enrichment and mass spectrometry to study proteins co-localized with RNAPII in precursor neurons differentiated from VCP-mutant or control induced pluripotent stem cells. Remarkably, we observed diminished RNAPII binding of proteins involved in transcription elongation and mRNA splicing in mutant cells. One of these is SART3, a recycling factor of the splicing machinery, whose knockdown leads to perturbed intron retention in several ALS-associated genes. Additional reduced proteins are RBM45, EIF5A and RNF220, mutations in which are associated with various neurodegenerative disorders and are linked to TDP-43 aggregation. Conversely, we observed increased RNAPII binding of heat shock proteins such as HSPB1. Together, these findings shed light on how transcription and splicing machinery are impaired by VCP mutations, which might contribute to aberrant alternative splicing and proteinopathy in neurodegeneration.

SUMMARYGenes are not randomly distributed throughout chromosomes. How gene order evolves and how selective constraints act to preserve or vary gene order, both at the macrosyntenic level of whole chromosomes or microsyntenic level of gene blocks, are central questions of evolutionary biology and genomics that remain largely unsolved. Here, after sequencing several genomes of the appendicularian tunicateOikopleura dioicafrom different locations around the globe, we show an unprecedented amount of genome scrambling in animals with no obvious morphological differences, consistent with cryptic speciation. Our assemblies suggest that all members of this clade possess a common 3-chromosome karyotype, and that different species largely preserve gene content, despite the presence of thousands of rearrangements in gene order. The movements of genes are largely restricted to chromosome arms and sex-specific regions, which appear to be the primary unit of macrosynteny conservation, and examples of these within-arm movements can be seen in theHoxandFgfgene families. Our approach employing whole-genome alignments demonstrates that segments containing protein-coding elements tend to be preserved at the microsyntenic scale, consistent with strong purifying selection, with appreciably less preservation of non-coding elements. Unexpectedly, scrambling did not preserve operon structure across species, suggesting an absence of selective pressure to maintain operon structure. As well, genome scrambling does not occur uniformly across all chromosomes, as short chromosome arms possess shorter genes, smaller operons, more breakpoints, and elevated dN/dS values compared to long chromosome arms. Estimation of divergence times among the crypticO. dioicalineages yielded an estimated breakpoint accumulation rate of 6 to 25 breakpoints per megabase per million years, which is an order of magnitude higher than the rates for other ascidian tunicates orDrosophilaspecies. Therefore,O. dioicaappears to be an attractive animal system to unravel the mechanisms that underlie gene order and synteny conservation, as well as exploring the limits of genome scrambling without an apparent impact on phenotypic evolution.

During early vertebrate development, signals from a special region of the embryo, the organizer, can redirect the fate of non-neural ectoderm cells to form a complete, patterned nervous system. This is called neural induction and has generally been imagined as a single signalling event, causing a switch of fate. Here, we undertake a comprehensive analysis, in very fine time course, of the events following exposure of competent ectoderm of the chick to the organizer (the tip of the primitive streak, Hensen’s node). Using transcriptomics and epigenomics we generate a gene regulatory network comprising 175 transcriptional regulators and 5614 predicted interactions between them, with fine temporal dynamics from initial exposure to the signals to expression of mature neural plate markers. Using in situ hybridization, single-cell RNA-sequencing, and reporter assays, we show that the gene regulatory hierarchy of responses to a grafted organizer closely resembles the events of normal neural plate development. The study is accompanied by an extensive resource, including information about conservation of the predicted enhancers in other vertebrates.

Integrator is a multi-subunit complex that directly interacts with the C-terminal domain (CTD) of RNA polymerase II (RNAPII). Through its RNA endonuclease activity, Integrator is required for 3′-end processing of both non-coding and coding transcripts. Here we demonstrate that depleting Integrator subunit 11 (INTS11), the main catalytic subunit of the Integrator complex, leads to a global elongation defect as a result of decreased polymerase processivity. We observe this defect in the region approximately 12 to 35 kb downstream of the transcription start site (TSS), where RNAPII normally transitions to its maximum processivity. We also identify an important role for INTS11, possibly in association with RNAPII CTD phospho-Tyr1, in repressing antisense transcription upstream of active promoters, as well as repressing transcription of genic regions near AsiSI-induced double-strand breaks.Altogether, this study points toward a novel function of Integrator in promoting termination of incompetent RNAPII molecules while facilitating the transition to fully processive polymerase in order to enable efficient elongation.

RNA-binding proteins (RBPs) play diverse roles in regulating co-transcriptional RNA-processing and chromatin functions, but our knowledge of the repertoire of chromatin-associated RBPs (caRBPs) and their interactions with chromatin remains limited. Here, we developed SPACE (Silica Particle Assisted Chromatin Enrichment) to isolate global and regional chromatin components with high specificity and sensitivity, and SPACEmap to identify the chromatin-contact regions in proteins. Applied to mouse embryonic stem cells, SPACE identified 1459 chromatin-associated proteins, ∼48% of which are annotated as RBPs, indicating their dual roles in chromatin and RNA-binding. Additionally, SPACEmap stringently verified chromatin-binding of 403 RBPs and identified their chromatin-contact regions. Notably, SPACEmap showed that about 40% of the caRBPs bind chromatin by intrinsically disordered regions (IDRs). Studying SPACE and total proteome dynamics from mES cells grown in 2iL and serum medium indicates significant correlation (R = 0.62). One of the most dynamic caRBPs is Dazl, which we find co-localized with PRC2 at transcription start sites of genes that are distinct from Dazl mRNA binding. Dazl and other PRC2-colocalised caRBPs are rich in intrinsically disordered regions (IDRs), which could contribute to the formation and regulation of phase-separated PRC condensates. Together, our approach provides an unprecedented insight into IDR-mediated interactions and caRBPs with moonlighting functions in native chromatin., Graphical Abstract Graphical AbstractSPACE, SPACEmap and ChIP-SPACE are highly sensitive and stringent approaches for identification of chromatin-associated proteins and their chromatin-contact regions.

We recently described aberrantly increased cytoplasmic SFPQ intron-retaining transcripts (IRTs) and concurrent SFPQ protein mislocalization as new hallmarks of amyotrophic lateral sclerosis (ALS). However, the generalizability and potential roles of cytoplasmic IRTs in health and disease remain unclear. Here, using time-resolved deep sequencing of nuclear and cytoplasmic fractions of human induced pluripotent stem cells undergoing motor neurogenesis, we reveal that ALS-causing VCP gene mutations lead to compartment-specific aberrant accumulation of IRTs. Specifically, we identify >100 IRTs with increased cytoplasmic abundance in ALS samples. Furthermore, these aberrant cytoplasmic IRTs possess sequence-specific attributes and differential predicted binding affinity to RNA binding proteins. Remarkably, TDP-43, SFPQ and FUS—RNA binding proteins known for nuclear-to-cytoplasmic mislocalization in ALS—abundantly and specifically bind to this aberrant cytoplasmic pool of IRTs. Our data are therefore consistent with a novel role for cytoplasmic IRTs in regulating compartment-specific protein abundance. This study provides new molecular insight into potential pathomechanisms underlying ALS and highlights aberrant cytoplasmic IRTs as potential therapeutic targets., Aberrant cytoplasmic SFPQ intron-retaining transcripts (IRTs) have recently been described in ALS, but their roles remained unclear. Tyzack, Neeves et al. now identify >100 cytoplasmic IRTs and show that they are ‘blueprints’ for ALS-related RNA binding protein mislocalization, as well as potential therapeutic targets.

Background The larvacean Oikopleura dioica is an abundant tunicate plankton with the smallest (65–70 Mbp) non-parasitic, non-extremophile animal genome identified to date. Currently, there are two genomes available for the Bergen (OdB3) and Osaka (OSKA2016) O. dioica laboratory strains. Both assemblies have full genome coverage and high sequence accuracy. However, a chromosome-scale assembly has not yet been achieved. Results Here, we present a chromosome-scale genome assembly (OKI2018_I69) of the Okinawan O. dioica produced using long-read Nanopore and short-read Illumina sequencing data from a single male, combined with Hi-C chromosomal conformation capture data for scaffolding. The OKI2018_I69 assembly has a total length of 64.3 Mbp distributed among 19 scaffolds. 99% of the assembly is contained within five megabase-scale scaffolds. We found telomeres on both ends of the two largest scaffolds, which represent assemblies of two fully contiguous autosomal chromosomes. Each of the other three large scaffolds have telomeres at one end only and we propose that they correspond to sex chromosomes split into a pseudo-autosomal region and X-specific or Y-specific regions. Indeed, these five scaffolds mostly correspond to equivalent linkage groups in OdB3, suggesting overall agreement in chromosomal organization between the two populations. At a more detailed level, the OKI2018_I69 assembly possesses similar genomic features in gene content and repetitive elements reported for OdB3. The Hi-C map suggests few reciprocal interactions between chromosome arms. At the sequence level, multiple genomic features such as GC content and repetitive elements are distributed differently along the short and long arms of the same chromosome. Conclusions We show that a hybrid approach of integrating multiple sequencing technologies with chromosome conformation information results in an accurate de novo chromosome-scale assembly of O. dioica’s highly polymorphic genome. This genome assembly opens up the possibility of cross-genome comparison between O. dioica populations, as well as of studies of chromosomal evolution in this lineage.

Indirect development with an intermediate larva exists in all major animal lineages1, which makes larvae central to most scenarios of animal evolution2–11. Yet how larvae evolved remains disputed. Here we show that temporal shifts (that is, heterochronies) in trunk formation underpin the diversification of larvae and bilaterian life cycles. We performed chromosome-scale genome sequencing in the annelid Owenia fusiformis with transcriptomic and epigenomic profiling during the life cycles of this and two other annelids. We found that trunk development is deferred to pre-metamorphic stages in the feeding larva of O. fusiformis but starts after gastrulation in the non-feeding larva with gradual metamorphosis of Capitella teleta and the direct developing embryo of Dimorphilus gyrociliatus. Accordingly, the embryos of O. fusiformis develop first into an enlarged anterior domain that forms larval tissues and the adult head12. Notably, this also occurs in the so-called ‘head larvae’ of other bilaterians13–17, with which the O. fusiformis larva shows extensive transcriptomic similarities. Together, our findings suggest that the temporal decoupling of head and trunk formation, as maximally observed in head larvae, facilitated larval evolution in Bilateria. This diverges from prevailing scenarios that propose either co-option9,10 or innovation11 of gene regulatory programmes to explain larva and adult origins.

Nucleosome organization has been suggested to affect local mutation rates in the genome. However, the lack of de novo mutation and high-resolution nucleosome data has limited the investigation of this hypothesis. Additionally, analyses using indirect mutation rate measurements have yielded contradictory and potentially confounding results. Here, we combine data on >300,000 human de novo mutations with high-resolution nucleosome maps and find substantially elevated mutation rates around translationally stable (‘strong’) nucleosomes. We show that the mutational mechanisms affected by strong nucleosomes are low-fidelity replication, insufficient mismatch repair and increased double-strand breaks. Strong nucleosomes preferentially locate within young SINE/LINE transposons, suggesting that when subject to increased mutation rates, transposons are then more rapidly inactivated. Depletion of strong nucleosomes in older transposons suggests frequent positioning changes during evolution. The findings have important implications for human genetics and genome evolution., Nucleosome organization has been suggested to affect local mutation rates in the genome. Here, the authors analyse data on >300,000 human de novo mutations and high-resolution nucleosome maps and provide evidence that nucleosome positioning stability modulates germline mutation rate variation across the human genome.

Summary In the past, the skills required to make an accurate dermatological diagnosis have required exposure to thousands of patients over many years. However, in recent years, artificial intelligence (AI) has made enormous advances, particularly in the area of image classification. This has led computer scientists to apply these techniques to develop algorithms that are able to recognize skin lesions, particularly melanoma. Since 2017, there have been numerous studies assessing the accuracy of algorithms, with some reporting that the accuracy matches or surpasses that of a dermatologist. While the principles underlying these methods are relatively straightforward, it can be challenging for the practising dermatologist to make sense of a plethora of unfamiliar terms in this domain. Here we explain the concepts of AI, machine learning, neural networks and deep learning, and explore the principles of how these tasks are accomplished. We critically evaluate the studies that have assessed the efficacy of these methods and discuss limitations and potential ethical issues. The burden of skin cancer is growing within the Western world, with major implications for both population skin health and the provision of dermatology services. AI has the potential to assist in the diagnosis of skin lesions and may have particular value at the interface between primary and secondary care. The emerging technology represents an exciting opportunity for dermatologists, who are the individuals best informed to explore the utility of this powerful novel diagnostic tool, and facilitate its safe and ethical implementation within healthcare systems., What is already known about this topic? There is considerable interest in the application of artificial intelligence to medicine.Several publications in recent years have described computer algorithms that can diagnose melanoma or skin lesions.Multiple groups have independently evaluated algorithms for the diagnosis of melanoma and skin lesions. What does this study add? We combine an introduction to the field with a summary of studies comparing dermatologists against artificial intelligence algorithms with the aim of providing a comprehensive resource for clinicians.This review will equip clinicians with the relevant knowledge to critically appraise future studies, and also assess the clinical utility of this technology.A better informed and engaged cohort of clinicians will ensure that the technology is applied effectively and ethically. Plain language summary available online

N6-methyladenosine (m6A) RNA modification impacts mRNA fate primarily via reader proteins, which dictate processes in development, stress, and disease. Yet little is known about m6A function in Saccharomyces cerevisiae, which occurs solely during early meiosis. Here we perform a multifaceted analysis of the m6A reader protein Pho92/Mrb1. Cross-linking immunoprecipitation analysis reveals that Pho92 associates with the 3’end of meiotic mRNAs in both an m6A-dependent and independent manner. Within cells, Pho92 transitions from the nucleus to the cytoplasm, and associates with translating ribosomes. In the nucleus Pho92 associates with target loci through its interaction with transcriptional elongator Paf1C. Functionally, we show that Pho92 promotes and links protein synthesis to mRNA decay. As such, the Pho92-mediated m6A-mRNA decay is contingent on active translation and the CCR4-NOT complex. We propose that the m6A reader Pho92 is loaded co-transcriptionally to facilitate protein synthesis and subsequent decay of m6A modified transcripts, and thereby promotes meiosis.

Appendicularian tunicates are some of the most abundant mesozooplankton organisms with key roles in marine trophic webs and global carbon flux. Like most appendicularians with cosmopolitan distributions, Oikopleura dioica Fol, 1872 is considered a single species worldwide based on morphological features that distinguish them from other appendicularians. Despite their abundance however, there are still only ∼70 described appendicularian species, compared with over 2,800 ascidian tunicates. Here we perform a molecular phylogenetic, morphological, and reproductive assessment of O. dioica specimens collected from the Ryukyu Archipelago, mainland Japan, and Europe. The specimens are morphologically very similar, with only detailed examination of the oikoplastic epithelium and quantitative measurements revealing minor distinguishing characteristics. Phylogenetic analyses of the ribosomal gene loci and mitochondrial cytochrome oxidase I (COI) gene strongly indicate that they form three separate genetic clades despite their morphological similarities. Finally, in vitro crosses between the Ryukyu and mainland Japanese specimens show total prezygotic reproductive isolation. Our results reveal that the current taxonomic O. dioica classification likely hides multiple cryptic species, highlighting the genetic diversity and complexity of their population structures. Cryptic organisms are often hidden under a single species name because their morphological similarities make them difficult to disinguish and their correct identification is fundamental to understanding Earth’s biodiversity. O. dioica is an attractive model to understand how morphological conservation can be maintained despite pronounced genetic divergence.

N6- methyladenosine (m6A) RNA modification impacts mRNA fate primarily via reader proteins, which dictate processes in development, stress, and disease. Yet little is known about m6A function in Saccharomyces cerevisiae, which occurs solely during early meiosis. Here, we perform a multifaceted analysis of the m6A reader protein Pho92/Mrb1. Cross-linking immunoprecipitation analysis reveals that Pho92 associates with the 3’end of meiotic mRNAs in both an m6A-dependent and independent manner. Within cells, Pho92 transitions from the nucleus to the cytoplasm, and associates with translating ribosomes. In the nucleus Pho92 associates with target loci through its interaction with transcriptional elongator Paf1C. Functionally, we show that Pho92 promotes and links protein synthesis to mRNA decay. As such, the Pho92-mediated m6A-mRNA decay is contingent on active translation and the CCR4-NOT complex. We propose that the m6A reader Pho92 is loaded co-transcriptionally to facilitate protein synthesis and subsequent decay of m6A modified transcripts, and thereby promotes meiosis.

Summary Mutations causing amyotrophic lateral sclerosis (ALS) often affect the condensation properties of RNA-binding proteins (RBPs). However, the role of RBP condensation in the specificity and function of protein-RNA complexes remains unclear. We created a series of TDP-43 C-terminal domain (CTD) variants that exhibited a gradient of low to high condensation propensity, as observed in vitro and by nuclear mobility and foci formation. Notably, a capacity for condensation was required for efficient TDP-43 assembly on subsets of RNA-binding regions, which contain unusually long clusters of motifs of characteristic types and density. These “binding-region condensates” are promoted by homomeric CTD-driven interactions and required for efficient regulation of a subset of bound transcripts, including autoregulation of TDP-43 mRNA. We establish that RBP condensation can occur in a binding-region-specific manner to selectively modulate transcriptome-wide RNA regulation, which has implications for remodeling RNA networks in the context of signaling, disease, and evolution., Graphical abstract, Highlights • TDP-43 mutants affect condensation properties to a similar extent at multiple scales • Binding-region condensates form on long RNA regions with dispersed UG-rich motifs • RBPchimera-CLIP indicates homomeric interactions promote molecular-scale condensates • Condensation selectively tunes the regulatory capacity of TDP-43; e.g., autoregulation, The condensation propensity of an RNA-binding protein tunes its binding to specific RNA regions across the transcriptome and affects its RNA processing functions. Formation of these “binding-region condensates,” promoted by specific motif types that are dispersed across long RNA regions, expands the ways in which RNA binding can be selectively controlled beyond canonical RNA-binding domains.

RNA molecules undergo a vast array of chemical post-transcriptional modifications (PTMs) that can affect their structure and interaction properties. In recent years, a growing number of PTMs have been successfully mapped to the transcriptome using experimental approaches relying on high-throughput sequencing. Oxford Nanopore direct-RNA sequencing has been shown to be sensitive to RNA modifications. We developed and validated Nanocompore, a robust analytical framework that identifies modifications from these data. Our strategy compares an RNA sample of interest against a non-modified control sample, not requiring a training set and allowing the use of replicates. We show that Nanocompore can detect different RNA modifications with position accuracy in vitro, and we apply it to profile m6A in vivo in yeast and human RNAs, as well as in targeted non-coding RNAs. We confirm our results with orthogonal methods and provide novel insights on the co-occurrence of multiple modified residues on individual RNA molecules., Nanopore direct RNA Sequencing data contain information about the presence of RNA modifications, but their detection poses substantial challenges. Here the authors introduce Nanocompore, a new methodology for modification detection from Nanopore data.

CLIP technologies are now widely used to study RNA-protein interactions and many datasets are now publicly available. An important first step in CLIP data exploration is the visual inspection and assessment of processed genomic data on selected genes or regions and performing comparisons: either across conditions within a particular project, or incorporating publicly available data. However, the output files produced by data processing pipelines or preprocessed files available to download from data repositories are often not suitable for direct comparison and usually need further processing. Furthermore, to derive biological insight it is usually necessary to visualise CLIP signal alongside other data such as annotations, or orthogonal functional genomic data (e.g. RNA-seq). We have developed a simple, but powerful, command-line tool: clipplotr, which facilitates these visual comparative and integrative analyses with normalisation and smoothing options for CLIP data and the ability to show these alongside reference annotation tracks and functional genomic data. These data can be supplied as input to clipplotr in a range of file formats, which will output a publication quality figure. It is written in R and can both run on a laptop computer independently, or be integrated into computational workflows on a high-performance cluster. Releases, source code and documentation are freely available at: https://github.com/ulelab/clipplotr.

Astrocytes contribute to motor neuron death in amyotrophic lateral sclerosis (ALS), but whether they adopt deleterious features consistent with inflammatory reactive states remains incompletely resolved. To identify inflammatory reactive features in ALS human induced pluripotent stem cell (hiPSC)–derived astrocytes, we examined transcriptomics, proteomics, and glutamate uptake in VCP-mutant astrocytes. We complemented this by examining other ALS mutations and models using a systematic meta-analysis of all publicly-available ALS astrocyte sequencing data, which included hiPSC-derived astrocytes carrying SOD1, C9orf72, and FUS gene mutations as well as mouse ALS astrocyte models with SOD1G93A mutation, Tardbp deletion, and Tmem259 (also known as membralin) deletion. ALS astrocytes were characterized by up-regulation of genes involved in the extracellular matrix, endoplasmic reticulum stress, and the immune response and down-regulation of synaptic integrity, glutamate uptake, and other neuronal support processes. We identify activation of the TGFB, Wnt, and hypoxia signaling pathways in both hiPSC and mouse ALS astrocytes. ALS changes positively correlate with TNF, IL1A, and complement pathway component C1q-treated inflammatory reactive astrocytes, with significant overlap of differentially expressed genes. By contrasting ALS changes with models of protective reactive astrocytes, including middle cerebral artery occlusion and spinal cord injury, we uncover a cluster of genes changing in opposing directions, which may represent down-regulated homeostatic genes and up-regulated deleterious genes in ALS astrocytes. These observations indicate that ALS astrocytes augment inflammatory processes while concomitantly suppressing neuronal supporting mechanisms, thus resembling inflammatory reactive states and offering potential therapeutic targets.

Background: The first step of virtually all next generation sequencing analysis involves the splitting of the raw sequencing data into separate files using sample-specific barcodes, a process known as “demultiplexing”. However, we found that existing software for this purpose was either too inflexible or too computationally intensive for fast, streamlined processing of raw, single end fastq files containing combinatorial barcodes. Results: Here, we introduce a fast and uniquely flexible demultiplexer, named Ultraplex, which splits a raw FASTQ file containing barcodes either at a single end or at both 5’ and 3’ ends of reads, trims the sequencing adaptors and low-quality bases, and moves unique molecular identifiers (UMIs) into the read header, allowing subsequent removal of PCR duplicates. Ultraplex is able to perform such single or combinatorial demultiplexing on both single- and paired-end sequencing data, and can process an entire Illumina HiSeq lane, consisting of nearly 500 million reads, in less than 20 minutes. Conclusions: Ultraplex greatly reduces computational burden and pipeline complexity for the demultiplexing of complex sequencing libraries, such as those produced by various CLIP and ribosome profiling protocols, and is also very user friendly, enabling streamlined, robust data processing. Ultraplex is available on PyPi and Conda and via Github.

Background: The first step of virtually all next generation sequencing analysis involves the splitting of the raw sequencing data into separate files using sample-specific barcodes, a process known as “demultiplexing”. However, we found that existing software for this purpose was either too inflexible or too computationally intensive for fast, streamlined processing of raw, single end fastq files containing combinatorial barcodes. Results: Here, we introduce a fast and uniquely flexible demultiplexer, named Ultraplex, which splits a raw FASTQ file containing barcodes either at a single end or at both 5’ and 3’ ends of reads, trims the sequencing adaptors and low-quality bases, and moves unique molecular identifiers (UMIs) into the read header, allowing subsequent removal of PCR duplicates. Ultraplex is able to perform such single or combinatorial demultiplexing on both single- and paired-end sequencing data, and can process an entire Illumina HiSeq lane, consisting of nearly 500 million reads, in less than 20 minutes. Conclusions: Ultraplex greatly reduces computational burden and pipeline complexity for the demultiplexing of complex sequencing libraries, such as those produced by various CLIP and ribosome profiling protocols, and is also very user friendly, enabling streamlined, robust data processing. Ultraplex is available on PyPi and Conda and via Github.

Background Eukaryotic genomes undergo pervasive transcription, leading to the production of many types of stable and unstable RNAs. Transcription is not restricted to regions with annotated gene features but includes almost any genomic context. Currently, the source and function of most RNAs originating from intergenic regions in the human genome remain unclear. Results We hypothesize that many intergenic RNAs can be ascribed to the presence of as-yet unannotated genes or the “fuzzy” transcription of known genes that extends beyond the annotated boundaries. To elucidate the contributions of these two sources, we assemble a dataset of more than 2.5 billion publicly available RNA-seq reads across 5 human cell lines and multiple cellular compartments to annotate transcriptional units in the human genome. About 80% of transcripts from unannotated intergenic regions can be attributed to the fuzzy transcription of existing genes; the remaining transcripts originate mainly from putative long non-coding RNA loci that are rarely spliced. We validate the transcriptional activity of these intergenic RNAs using independent measurements, including transcriptional start sites, chromatin signatures, and genomic occupancies of RNA polymerase II in various phosphorylation states. We also analyze the nuclear localization and sensitivities of intergenic transcripts to nucleases to illustrate that they tend to be rapidly degraded either on-chromatin by XRN2 or off-chromatin by the exosome. Conclusions We provide a curated atlas of intergenic RNAs that distinguishes between alternative processing of well-annotated genes from independent transcriptional units based on the combined analysis of chromatin signatures, nuclear RNA localization, and degradation pathways.

See Vidal and Atkin (doi:10.1093/brain/awz256) for a scientific commentary on this article. ALS remains incurable due to an incomplete understanding of its molecular pathogenesis. Tyzack, Luisier et al. report that FUS protein, previously thought to mislocalise only in patients with FUS mutations, is also mislocalised in many cases of sporadic ALS. They further propose an underlying molecular mechanism for the mislocalisation., Mutations causing amyotrophic lateral sclerosis (ALS) clearly implicate ubiquitously expressed and predominantly nuclear RNA binding proteins, which form pathological cytoplasmic inclusions in this context. However, the possibility that wild-type RNA binding proteins mislocalize without necessarily becoming constituents of cytoplasmic inclusions themselves remains relatively unexplored. We hypothesized that nuclear-to-cytoplasmic mislocalization of the RNA binding protein fused in sarcoma (FUS), in an unaggregated state, may occur more widely in ALS than previously recognized. To address this hypothesis, we analysed motor neurons from a human ALS induced-pluripotent stem cell model caused by the VCP mutation. Additionally, we examined mouse transgenic models and post-mortem tissue from human sporadic ALS cases. We report nuclear-to-cytoplasmic mislocalization of FUS in both VCP-mutation related ALS and, crucially, in sporadic ALS spinal cord tissue from multiple cases. Furthermore, we provide evidence that FUS protein binds to an aberrantly retained intron within the SFPQ transcript, which is exported from the nucleus into the cytoplasm. Collectively, these data support a model for ALS pathogenesis whereby aberrant intron retention in SFPQ transcripts contributes to FUS mislocalization through their direct interaction and nuclear export. In summary, we report widespread mislocalization of the FUS protein in ALS and propose a putative underlying mechanism for this process.

Chromatin transactions are typically studied in vivo, or in vitro using artificial chromatin lacking the epigenetic complexity of the natural material. Attempting to bridge the gap between these approaches, we established a system for isolating the yeast genome as a library of mononucleosomes harboring the natural epigenetic signature, suitable for biochemical manipulation. Combined with deep sequencing, this library was used to investigate the stability of individual nucleosomes and, as proof of principle, the nucleosome preference of the chromatin remodeling complex, RSC. This approach uncovered a distinct preference of RSC for nucleosomes derived from regions with a high density of histone variant H2AZ, and this preference is indeed markedly diminished using nucleosomes from cells lacking H2AZ. The preference for H2AZ remodeling/nucleosome ejection can also be reconstituted with recombinant nucleosome arrays. Together, our data indicate that, despite being separated from their genomic context, individual nucleosomes can retain their original identity as promoter- or transcription start site (TSS)-nucleosomes. Besides shedding new light on substrate preference of the chromatin remodeler RSC, the simple experimental system outlined here should be generally applicable to the study of chromatin transactions.