Search

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.

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 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.

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.

During early vertebrate development, signals from a special region of the embryo, the organizer, can re-direct 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 ectoderm to the organizer. Using transcriptomics and epigenomics we generate a Gene Regulatory Network comprising 175 transcriptional regulators and 5,614 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 neural induction by a grafted organizer closely resembles normal neural plate development. The study is accompanied by a comprehensive resource including information about conservation of the predicted enhancers in different vertebrate systems.

RNA helicases remodel the spliceosome to enable pre-mRNA splicing, but their binding and mechanism of action remain poorly understood. To define helicase-RNA contacts in specific spliceosomal states, we develop purified spliceosome iCLIP (psiCLIP), which reveals dynamic helicase-RNA contacts during splicing catalysis. The helicase Prp16 binds along the entire available single-stranded RNA region between the branchpoint and 3′-splice site, while Prp22 binds diffusely downstream of the branchpoint before exon ligation, but then switches to more narrow binding in the downstream exon after exon ligation, arguing against a mechanism of processive translocation. Depletion of the exon-ligation factor Prp18 destabilizes Prp22 binding to the pre-mRNA, suggesting that proofreading by Prp22 may sense the stability of the spliceosome during exon ligation. Thus, psiCLIP complements structural studies by providing key insights into the binding and proofreading activity of spliceosomal RNA helicases., ATP-dependent helicases remodel the spliceosome and proofread splice site recognition. A new method – Purified Spliceosome iCLIP (psiCLIP) – probes protein-RNA interactions in defined spliceosome complexes to reveal how the helicases Prp16 and Prp22 promote correct mRNA synthesis through dynamic binding on their RNA substrates.

Oikopleura dioica is a ubiquitous marine zooplankton of biological interest owing to features that include dioecious reproduction, a short life cycle, conserved chordate body plan, and a compact genome. It is an important tunicate model for evolutionary and developmental research, as well as investigations into marine ecosystems. The genome of north Atlantic O. dioica comprises three chromosomes. However, comparisons with the genomes of O. dioica sampled from mainland and southern Japan revealed extensive sequence differences. Moreover, historical studies have reported widely varying chromosome counts. We recently initiated a project to study the genomes of O. dioica individuals collected from the coastline of the Ryukyu (Okinawa) Islands in southern Japan. Given the potentially large extent of genomic diversity, we employed karyological techniques to count individual animals’ chromosomes in situ using centromere-specific antibodies directed against H3S28P, a prophase-metaphase cell cycle-specific marker of histone H3. Epifluorescence and confocal images were obtained of embryos and oocytes stained with two commercial anti-H3S28P antibodies (Abcam ab10543 and Thermo Fisher 07-145). The data lead us to conclude that diploid cells from Okinawan O. dioica contain three pairs of chromosomes, in line with the north Atlantic populations. The finding facilitates the telomere-to-telomere assembly of Okinawan O. dioica genome sequences and gives insight into the genomic diversity of O. dioica from different geographical locations. The data deposited in the EBI BioImage Archive provide representative images of the antibodies’ staining properties for use in epifluorescent and confocal based fluorescent microscopy.

Reactive astrocytes are implicated in amyotrophic lateral sclerosis (ALS), although the mechanisms controlling reactive transformation are unknown. We show that decreased intron retention (IR) is common to human-induced pluripotent stem cell (hiPSC)-derived astrocytes carrying ALS-causing mutations in VCP, SOD1 and C9orf72. Notably, transcripts with decreased IR and increased expression are overrepresented in reactivity processes including cell adhesion, stress response and immune activation. This was recapitulated in public-datasets for (i) hiPSC-derived astrocytes stimulated with cytokines to undergo reactive transformation and (ii) in vivo astrocytes following selective deletion of TDP-43. We also re-examined public translatome sequencing (TRAP-seq) of astrocytes from a SOD1 mouse model, which revealed that transcripts upregulated in translation significantly overlap with transcripts exhibiting decreased IR. Using nucleocytoplasmic fractionation of VCP mutant astrocytes coupled with mRNA sequencing and proteomics, we identify that decreased IR in nuclear transcripts is associated with enhanced nonsense mediated decay and increased cytoplasmic expression of transcripts and proteins regulating reactive transformation. These findings are consistent with a molecular model for reactive transformation in astrocytes whereby poised nuclear reactivity-related IR transcripts are spliced, undergo nuclear-to-cytoplasmic translocation and translation. Our study therefore provides new insights into the molecular regulation of reactive transformation in astrocytes., Graphical Abstract Graphical abstractHealthy astrocytes (left, blue) have prevalent nuclear intron retention in reactivity transcripts, promoting their nuclear confinement (introns shown as red rectangles). Conversely, ALS astrocytes (right, red) exhibit decreased intron retention in reactivity transcripts, enabling cytoplasmic translocation and translation upon engaging ribosomes, underlying their reactive transformation. A small number of intron retaining transcripts escape to the cytoplasm where they are degraded by nonsense mediated decay (NMD), which is enhanced in ALS astrocytes.

Background The start and end sites of messenger RNAs (TSSs and TESs) are highly regulated, often in a cell-type-specific manner. Yet the contribution of transcript diversity in regulating gene expression remains largely elusive. We perform an integrative analysis of multiple highly synchronized cell-fate transitions and quantitative genomic techniques in Saccharomyces cerevisiae to identify regulatory functions associated with transcribing alternative isoforms. Results Cell-fate transitions feature widespread elevated expression of alternative TSS and, to a lesser degree, TES usage. These dynamically regulated alternative TSSs are located mostly upstream of canonical TSSs, but also within gene bodies possibly encoding for protein isoforms. Increased upstream alternative TSS usage is linked to various effects on canonical TSS levels, which range from co-activation to repression. We identified two key features linked to these outcomes: an interplay between alternative and canonical promoter strengths, and distance between alternative and canonical TSSs. These two regulatory properties give a plausible explanation of how locally transcribed alternative TSSs control gene transcription. Additionally, we find that specific chromatin modifiers Set2, Set3, and FACT play an important role in mediating gene repression via alternative TSSs, further supporting that the act of upstream transcription drives the local changes in gene transcription. Conclusions The integrative analysis of multiple cell-fate transitions suggests the presence of a regulatory control system of alternative TSSs that is important for dynamic tuning of gene expression. Our work provides a framework for understanding how TSS heterogeneity governs eukaryotic gene expression, particularly during cell-fate changes.

RNA-binding proteins (RBPs) act as posttranscriptional regulators controlling the fate of target mRNAs. Unraveling how RNAs are recognized by RBPs and in turn are assembled into neuronal RNA granules is therefore key to understanding the underlying mechanism. While RNA sequence elements have been extensively characterized, the functional impact of RNA secondary structures is only recently being explored. Here, we show that Staufen2 binds complex, long-ranged RNA hairpins in the 3′-untranslated region (UTR) of its targets. These structures are involved in the assembly of Staufen2 into RNA granules. Furthermore, we provide direct evidence that a defined Rgs4 RNA duplex regulates Staufen2-dependent RNA localization to distal dendrites. Importantly, disrupting the RNA hairpin impairs the observed effects. Finally, we show that these secondary structures differently affect protein expression in neurons. In conclusion, our data reveal the importance of RNA secondary structure in regulating RNA granule assembly, localization and eventually translation. It is therefore tempting to speculate that secondary structures represent an important code for cells to control the intracellular fate of their mRNAs.

Transcription factors (TFs) regulate gene expression by recognising and binding specific DNA sequences. At times, these regulatory elements may be occluded by nucleosomes, making them inaccessible for TF-binding. The competition for DNA occupancy between TFs and nucleosomes, and associated gene regulatory outputs, are important consequences of the cis-regulatory information encoded in the genome. However, these sequence patterns are subtle and remain difficult to interpret. Here, we introduce ChromWave, a deep-learning model that, for the first time, predicts the competing profiles for TF and nucleosomes occupancies with remarkable accuracy. Models trained using short- and long-fragment MNase-Seq data successfully learn the sequence preferences underlying TF and nucleosome occupancies across the entire yeast genome. They recapitulate nucleosome evictions from regions containing “strong” TF binding sites and knock-out simulations show nucleosomes gaining occupancy in the absence of these TFs, accompanied by lateral rearrangement of adjacent nucleosomes. At a local level, models anticipate with high accuracy the outcomes of detailed experimental analysis of partially unwrapped nucleosomes at the GAL4 UAS locus. Finally, we trained a ChromWave model that successfully predicts nucleosome positions at promoters in the human genome. We find that human promoters generally contain few sites at which simple sequence changes can alter nucleosome occupancies and that these positions align well with causal variants linked to DNase hypersensitivity. ChromWave is readily combined with diverse genomic datasets and can be trained to predict any output that is linked to the underlying genomic sequence. ChromWave’s application is limited only by the user's imagination and availability of training data.

Chromosomal instability in cancer consists of dynamic changes to the number and structure of chromosomes(1,2). The resulting diversity in somatic copy number alterations (SCNAs) may provide the variation necessary for tumour evolution(1,3,4). Here we use multi-sample phasing and SCNA analysis of 1,421 samples from 394 tumours across 22 tumour types to show that continuous chromosomal instability results in pervasive SCNA heterogeneity. Parallel evolutionary events, which cause disruption in the same genes (such asBCL9, MCL1,ARNT(also known asHIF1B),TERTandMYC) within separate subclones, were present in 37% of tumours. Most recurrent losses probably occurred before whole-genome doubling, that was found as a clonal event in 49% of tumours. However, loss of heterozygosity at the human leukocyte antigen (HLA) locus and loss of chromosome 8p to a single haploid copy recurred at substantial subclonal frequencies, even in tumours with whole-genome doubling, indicating ongoing karyotype remodelling. Focal amplifications that affected chromosomes 1q21 (which encompassesBCL9, MCL1andARNT), 5p15.33 (TERT), 11q13.3 (CCND1), 19q12 (CCNE1) and 8q24.1 (MYC) were frequently subclonal yet appeared to be clonal within single samples. Analysis of an independent series of 1,024 metastatic samples revealed that 13 focal SCNAs were enriched in metastatic samples, including gains in chromosome 8q24.1 (encompassingMYC) in clear cell renal cell carcinoma and chromosome 11q13.3 (encompassingCCND1) in HER2(+)breast cancer. Chromosomal instability may enable the continuous selection of SCNAs, which are established as ordered events that often occur in parallel, throughout tumour evolution.Chromosomal instability enables the continuous selection of somatic copy number alterations, which are established as ordered events that often occur in parallel, throughout tumour evolution and metastasis.

BackgroundThe 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.ResultsHere, 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 in 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 of 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.ConclusionsWe 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 assembly will be a useful resource for genome-wide comparative studies between O. dioica and other species, as well as studies of chromosomal evolution in this lineage.

SUMMARYWe recently described aberrant cytoplasmic SFPQ intron-retaining transcripts (IRTs) and concurrent SFPQ protein mislocalization as a new hallmark 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 hiPSCs 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 (but not nuclear) abundance in ALS samples. Furthermore, these aberrant cytoplasmic IRTs possess sequence-specific attributes and differential predicted binding affinity to RNA binding proteins (RBPs). Remarkably, TDP-43, SFPQ and FUS – RBPs known for nuclear-to-cytoplasmic mislocalization in ALS – avidly and specifically bind to this aberrant cytoplasmic pool of IRTs, as opposed to any individual IRT. 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.Abstract Figure

Oikopleura dioica is a ubiquitous marine tunicate of biological interest due to features that include dioecious reproduction, short life cycle, and vertebrate-like dorsal notochord while possessing a relatively compact genome. The use of tunicates as model organisms, particularly with these characteristics, offers the advantage of facilitating studies in evolutionary development and furthering understanding of enduring attributes found in the more complex vertebrates. At present, we are undertaking an initiative to sequence the genomes of Oikopleura individuals in populations found among the seas surrounding the Ryukyu Islands in southern Japan. To facilitate and validate genome assemblies, karyotyping was employed to count individual animals’ chromosomes in situ using centromere-specific antibodies directed against H3S28P, a prophase-metaphase cell cycle-specific marker of histone H3. New imaging data of embryos and oocytes stained with two different antibodies were obtained; interpretation of these data lead us to conclude that the Okinawan Oikopleura dioica has three pairs of chromosomes, akin to previous results from genomic assemblies in Atlantic populations. The imaging data have been deposited to the open-access EBI BioImage Archive for reuse while additionally providing representative images of two commercially available anti-H3S28P antibodies’ staining properties for use in epifluorescent and confocal based fluorescent microscopy.

Oikopleura dioica is a planktonic chordate with exceptional filter-feeding ability, rapid generation time, conserved early development, and a compact genome. For these reasons, it is considered a useful model organism for marine ecological studies, evolutionary developmental biology, and genomics. As research often requires a steady supply of animal resources, it is useful to establish a reliable, low-maintenance culture system. Here we describe a step-by-step method for establishing an O. dioica culture. We describe how to select potential sampling sites, collection methods, target animal identification, and the set-up of the culturing system. We provide troubleshooting advice based on our own experiences. We also highlight critical factors that help sustain a robust culture system. Although the culture protocol provided here is optimized for O. dioica, we hope our sampling technique and culture setup will inspire new ideas for maintaining other fragile pelagic invertebrates.

Coronavirus disease 2019 (COVID-19) is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and is associated with a wide range of systemic manifestations. Several observations support a role for vascular endothelial dysfunction in the pathogenesis including an increased incidence of thrombotic events and coagulopathy and the presence of vascular risk factors as an independent predictor of poor prognosis. It has recently been reported that endothelitis is associated with viral inclusion bodies suggesting a direct role for SARS-CoV-2 in the pathogenesis. The ACE2 receptor has been shown to mediate SARS-CoV-2 uptake and it has been proposed that CD147 (BSG) can function as an alternative cell surface receptor. To define the endothelial cell populations that are susceptible to infection with SARS-CoV-2, we investigated the expression of ACE2 as well as other genes implicated in the cellular entry of SARS-Cov-2 in the vascular endothelium through the analysis of single cell sequencing data derived from multiple human tissues (skin, liver, kidney, lung and intestine). We found that CD147 (BSG) but not ACE2 is detectable in vascular endothelial cells within single cell sequencing datasets derived from multiple tissues in healthy individuals. This implies that either ACE2 is not expressed in healthy tissue but is instead induced in response to SARS-Cov2 or that SARS-Cov2 enters endothelial cells via an alternative receptor such as CD147.

Many long noncoding RNAs (lncRNAs) act in cis through transcription-coupled chromatin alterations that drive changes in local gene expression. How some cis-acting lncRNAs promote and others repress gene expression remains poorly understood. Here we report that in S. cerevisiae transcription levels of the lncRNA IRT2, located upstream in the promoter of the inducer of meiosis gene, regulate opposing chromatin and transcription states. Low IRT2 transcription displays enhancer RNA-like features. At these levels, IRT2 promotes histone exchange delivering acetylated histone H3 lysine 56 to chromatin thereby facilitating recruitment of a transcription factor and consequently activating transcription. Conversely, increasing IRT2 transcription enhances chromatin assembly and transcriptional repression. The opposing functions of IRT2 direct a regulatory circuit, which ensures that cells expressing opposite, but not one of either, mating-type loci enter meiosis. Our data demonstrate that the transcription levels of an lncRNA are key to controlling gene expression and cell fate outcomes.

BackgroundAscidians, a tunicate class, use a mitochondrial genetic code that is distinct from vertebrates and other invertebrates. Though it has been used to translate the coding sequences from other tunicate species on a case-by-case basis, it is has not been investigated whether this can be done systematically. This is an important because a) some tunicate mitochondrial sequences are currently translated with the invertebrate code by repositories such as NCBI’s GenBank, and b) uncertainties about the genetic code to use can complicate or introduce errors in phylogenetic studies based on translated mitochondrial protein sequences.MethodsWe collected publicly available nucleotide sequences for non-ascidian tunicates including appendicularians such asOikopleura dioica, translated them using the ascidian mitochondrial code, and built multiple sequence alignments covering all tunicate classes.ResultsAll tunicates studied here appear to translate AGR codons to glycine instead of serine (invertebrates) or as a stop codon (vertebrates), as initially described in ascidians. Among Oikopleuridae, we suggest further possible changes in the use of the ATA (Ile → Met) and TGA (Trp → Arg) codons.ConclusionsWe recommend using the ascidian mitochondrial code in automatic translation pipelines of mitochondrial sequences for all tunicates. Further investigation is required for additional species-specific differences.

The structure of RNA molecules has a critical role in regulating gene expression, largely through influencing their interactions with RNA-binding proteins (RBPs). RNA hybrid and individual-nucleotide resolution UV cross-linking and immunoprecipitation (hiCLIP) is a transcriptome-wide method of monitoring these interactions by identifying RNA duplexes bound by a specific RBP. The hiCLIP protocol consists of the following steps: in vivo cross-linking of RBPs to their bound RNAs; partial RNA digestion and purification of RNA duplexes interacting with the specific RBP using immunoprecipitation; ligation of the two arms of RNA duplexes via a linker; reverse transcription; cDNA library amplification; and finally high-throughput DNA sequencing. Mapping of the sequenced arms to a reference transcriptome identifies the exact locations of duplexes. hiCLIP data can directly identify all types of RNA duplexes bound by RBPs, including those that are challenging to predict computationally, such as intermolecular and long-range intramolecular duplexes. Moreover, the use of an adaptor that links the two arms of the RNA duplex permits hiCLIP to unambiguously identify the duplexes. Here we describe in detail the procedure for a hiCLIP experiment and the subsequent streamlined data analysis with an R package, 'hiclipr' (https://github.com/luslab/hiclipr/). Preparation of the library for high-throughput DNA sequencing takes ∼7 d and the basic bioinformatic pipeline takes 1 d.

Rare amyotrophic lateral sclerosis (ALS)-causing mutations clearly implicate ubiquitously expressed and predominantly nuclear RNA binding proteins (RBPs), which form pathological cytoplasmic inclusions in this context. However, it remains possible that RBPs can mislocalize without necessarily becoming constituents of cytoplasmic ALS inclusions. We hypothesized that nuclear-to-cytoplasmic mislocalization of the RBP Fused in Sarcoma (FUS), in an unaggregated state, may occur more widely in ALS that previously recognized. To address this hypothesis, we analysed motor neurons (MNs) from human control and patient-specific VCP mutant induced-pluripotent stem cells (iPSCs), in addition to mouse transgenic models and post-mortem tissue from human sporadic ALS. We report nuclear-to-cytoplasmic mislocalization of FUS in both VCP-mutation related ALS and in sporadic ALS spinal cord tissue from multiple cases. In addition to finding evidence of associated DNA damage, we further link this mislocalization of FUS protein mechanistically to aberrant intron retention in the SFPQ transcript. Specifically, FUS protein binds extensively to the aberrantly retained intron within the SFPQ transcript, which is exported from the nucleus into the cytoplasm. Collectively, these data support a model 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.

Understanding the patterns and genesis of germlinede novomutations is important for studying genome evolution and human diseases. Nucleosome organization is suggested to be a contributing factor to mutation rate variation across the genome. However, the small number of publishedde novomutations and the low resolution of earlier nucleosome maps limited our understanding of how nucleosome organization affects germline mutation rates in the human genome. Here, we systematically investigated the relationship between nucleosome organization and fine-scale mutation rate variation by analyzing >300,000de novomutations from whole-genome trio sequencing and high-resolution nucleosome maps in human. We found thatde novomutation rates are elevated around strong, translationally stable nucleosomes, a previously under-appreciated aspect. We confirmed this observation having controlled for local sequence context and other potential confounding factors. Analysis of the underlying mutational processes suggests that the increased mutation rates around strong nucleosomes are shaped by a combination of low-fidelity replication, frequent DNA damage and insufficient/error-prone repair in these regions. Interestingly, strong nucleosomes are preferentially located in young SINE/LINE elements, implying frequent nucleosome re-positioning (i.e. shifting of dyad position) and their contribution to hypermutation at new retrotransposons during evolution. These findings provide novel insights into how chromatin organization affects germline mutation rates and have important implications in human genetics and genome evolution.

Summary The CRISPR-Cas9 system has successfully been adapted to edit the genome of various organisms. However, our ability to predict the editing outcome at specific sites is limited. Here, we examined indel profiles at over 1,000 genomic sites in human cells and uncovered general principles guiding CRISPR-mediated DNA editing. We find that precision of DNA editing (i.e., recurrence of a specific indel) varies considerably among sites, with some targets showing one highly preferred indel and others displaying numerous infrequent indels. Editing precision correlates with editing efficiency and a preference for single-nucleotide homologous insertions. Precise targets and editing outcome can be predicted based on simple rules that mainly depend on the fourth nucleotide upstream of the protospacer adjacent motif (PAM). Indel profiles are robust, but they can be influenced by chromatin features. Our findings have important implications for clinical applications of CRISPR technology and reveal general patterns of broken end joining that can provide insights into DNA repair mechanisms., Graphical Abstract, Highlights • The outcome of CRISPR-mediated editing can be predicted • Not all target sites are edited in a predictable manner • The precision of DNA editing is mainly determined by the fourth nucleotide upstream of the PAM site • Chromatin states affect editing of imprecise, but not precise, target sites, Chakrabarti, Henser-Brownhill, Monserrat et al. show that the genome-editing outcome can be predicted based on simple rules that mainly depend on the target site sequence. Since editing precision varies considerably across sites, careful selection of a predictable target is critical to induce a desired modification in a cell-type-independent manner.

Summary The 3′ untranslated regions (3′ UTRs) of messenger RNAs (mRNAs) are non-coding sequences involved in many aspects of mRNA metabolism, including intracellular localization and translation. Incorrect processing and delivery of mRNA cause severe developmental defects and have been implicated in many neurological disorders. Here, we use deep sequencing to show that in sympathetic neuron axons, the 3′ UTRs of many transcripts undergo cleavage, generating isoforms that express the coding sequence with a short 3′ UTR and stable 3′ UTR-derived fragments of unknown function. Cleavage of the long 3′ UTR of Inositol Monophosphatase 1 (IMPA1) mediated by a protein complex containing the endonuclease argonaute 2 (Ago2) generates a translatable isoform that is necessary for maintaining the integrity of sympathetic neuron axons. Thus, our study provides a mechanism of mRNA metabolism that simultaneously regulates local protein synthesis and generates an additional class of 3′ UTR-derived RNAs., Graphical Abstract, Highlights • Axons and cell bodies of sympathetic neurons express distinct 3′ UTR isoforms • Axon-specific short 3′ UTR isoforms are generated by local cleavage of longer 3′ UTRs • A protein complex containing Ago2, Upf1, HuD, and Pabpc4 mediates the 3′ UTR cleavage, Andreassi et al. show widespread differential usage of 3′ UTR in axons and cell bodies of sympathetic neurons. In axons, the cleavage of a longer 3′ UTR of Impa1 generates a shorter isoform that is stable, polyadenylated, and necessary for maintaining axon integrity.

Summary Transcription through noncoding regions of the genome is pervasive. How these transcription events regulate gene expression remains poorly understood. Here, we report that, in S. cerevisiae, the levels of transcription through a noncoding region, IRT2, located upstream in the promoter of the inducer of meiosis, IME1, regulate opposing chromatin and transcription states. At low levels, the act of IRT2 transcription promotes histone exchange, delivering acetylated histone H3 lysine 56 to chromatin locally. The subsequent open chromatin state directs transcription factor recruitment and induces downstream transcription to repress the IME1 promoter and meiotic entry. Conversely, increasing transcription turns IRT2 into a repressor by promoting transcription-coupled chromatin assembly. The two opposing functions of IRT2 transcription shape a regulatory circuit, which ensures a robust cell-type-specific control of IME1 expression and yeast meiosis. Our data illustrate how intergenic transcription levels are key to controlling local chromatin state, gene expression, and cell fate outcomes., Graphical Abstract, Highlights • Transcription levels of the lncRNA IRT2 regulate opposing local transcription states • Low levels of IRT2 promote H3K56ac, chromatin disassembly, and TF recruitment • Increasing IRT2 levels promote chromatin assembly and transcriptional repression • IRT2 transcription shapes the regulatory circuit for cell-type-specific control of yeast meiosis, Moretto et al. demonstrate that transcription levels of a lncRNA named IRT2 can have opposing effects on local gene transcription, thereby controlling entry into yeast meiosis. Low levels of IRT2 promote H3K56ac, chromatin disassembly, transcription factor recruitment, and, thus, transcriptional activation. Increasing levels promote chromatin assembly and, consequently, transcriptional repression.

An amendment to this paper has been published and can be accessed via the original article.

Backmasking is a recording technique used to hide a sound or message in a music track in reverse, meaning that it is only audible when the record is played backwards. Analogously, the compact yeast genome encodes for diverse sources of information such as overlapping coding and non-coding transcripts, and protein-binding sites on the two complementary DNA strands. Examples are the consensus binding site sequences of the RNA-binding proteins Nrd1 and Nab3 that target non-coding transcripts for degradation. Here, by examining the overlap of stable (SUTs, stable unannotated transcripts) and unstable (CUTs, cryptic unstable transcripts) transcripts with protein-coding genes, we show that the predicted Nrd1 and Nab3-binding site sequences occur at differing frequencies. They are always depleted in the sense direction of protein-coding genes, thus avoiding degradation of the transcript. However in the antisense direction, predicted binding sites occur at high frequencies in genes with overlapping unstable ncRNAs (CUTs), so limiting the availability of non-functional transcripts. In contrast they are depleted in genes with overlapping stable ncRNAs (SUTs), presumably to avoid degrading the non-coding transcript. The protein-coding genes maintain similar amino-acid contents, but they display distinct codon usages so that Nrd1 and Nab3-binding sites can arise at differing frequencies in antisense depending on the overlapping transcript type. Our study demonstrates how yeast has evolved to encode multiple layers of information-protein-coding genes in one strand and the relative chance of degrading antisense RNA in the other strand-in the same regions of a compact genome.

Oikopleura dioica is a ubiquitous marine tunicate of biological interest due to features that include dioecious reproduction, short life cycle, and vertebrate-like dorsal notochord while possessing a relatively compact genome. The use of tunicates as model organisms, particularly with these characteristics, offers the advantage of facilitating studies in evolutionary development and furthering understanding of enduring attributes found in the more complex vertebrates. At present, we are undertaking an initiative to sequence the genomes of Oikopleura individuals in populations found among the seas surrounding the Ryukyu Islands in southern Japan. To facilitate and validate genome assemblies, karyotyping was employed to count individual animals’ chromosomes in situ using centromere-specific antibodies directed against H3S28P, a prophase-metaphase cell cycle-specific marker of histone H3. New imaging data of embryos and oocytes stained with two different antibodies were obtained; interpretation of these data lead us to conclude that the Okinawan Oikopleura dioica has three pairs of chromosomes, akin to previous results from genomic assemblies in Atlantic populations. The imaging data have been deposited to the open-access EBI BioImage Archive for reuse while additionally providing representative images of two commercially available anti-H3S28P antibodies’ staining properties for use in epifluorescent and confocal based fluorescent microscopy.

Single-cell RNA-seq has been established as a reliable and accessible technique enabling new types of analyses, such as identifying cell types and studying spatial and temporal gene expression variation and change at single-cell resolution. Recently, single-cell RNA-seq has been applied to developing embryos, which offers great potential for finding and characterising genes controlling the course of development along with their expression patterns. In this study, we applied single-cell RNA-seq to the 16-cell stage of the Ciona embryo, a marine chordate and performed a computational search for cell-specific gene expression patterns. We recovered many known expression patterns from our single-cell RNA-seq data and despite extensive previous screens, we succeeded in finding new cell-specific patterns, which we validated by in situ and single-cell qPCR.

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 mono-nucleosomes harboring the natural epigenetic signature, suitable for biochemical manipulation. Combined with deep sequencing, this library was used to investigate the intrinsic stability of individual nucleosomes, and – as proof of principle - the nucleosome preference of the chromatin remodeling complex, RSC. Our data indicate that the natural stability of nucleosomes differs greatly, with nucleosomes on tRNA genes and on promoters of protein-coding genes standing out as intrinsically unstable. Interestingly, RSC shows a distinct preference for nucleosomes derived from regions with a high density of histone variant H2A.Z, and this preference is indeed markedly diminished using nucleosomes from cells lacking H2A.Z. Importantly, the preference for H2A.Z 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 TSS-nucleosomes. Besides shedding new light on nucleosome stability and the chromatin remodeler RSC, the simple experimental system outlined here should be generally applicable to the study of chromatin transactions.

The transcription factor c-Maf induces the anti-inflammatory cytokine IL-10 in CD4+ T cells in vitro. However, the global effects of c-Maf on diverse immune responses in vivo are unknown. Here we found that c-Maf regulated IL-10 production in CD4+ T cells in disease models involving the TH1 subset of helper T cells (malaria), TH2 cells (allergy) and TH17 cells (autoimmunity) in vivo. Although mice with c-Maf deficiency targeted to T cells showed greater pathology in TH1 and TH2 responses, TH17 cell–mediated pathology was reduced in this context, with an accompanying decrease in TH17 cells and increase in Foxp3+ regulatory T cells. Bivariate genomic footprinting elucidated the c-Maf transcription-factor network, including enhanced activity of NFAT; this led to the identification and validation of c-Maf as a negative regulator of IL-2. The decreased expression of the gene encoding the transcription factor RORγt (Rorc) that resulted from c-Maf deficiency was dependent on IL-2, which explained the in vivo observations. Thus, c-Maf is a positive and negative regulator of the expression of cytokine-encoding genes, with context-specific effects that allow each immune response to occur in a controlled yet effective manner. The transcription factor c-Maf controls IL-10 production in T cells. O’Garra and colleagues use systems and in vivo functional analysis of T cell subsets to reveal distinct context-specific roles for c-Maf.

It is challenging for RNA processing machineries to select exons within long intronic regions. We find that intronic LINE repeat sequences (LINEs) contribute to this selection by recruiting dozens of RNA-binding proteins (RBPs). This includes MATR3, which promotes binding of PTBP1 to multivalent binding sites in LINEs. Both RBPs repress splicing and 3’ end processing within and around LINEs, as demonstrated in cultured human cells and mouse brain. Notably, repressive RBPs preferentially bind to evolutionarily young LINEs, which are confined to deep intronic regions. These RBPs insulate both LINEs and surrounding regions from RNA processing. Upon evolutionary divergence, gradual loss of insulation diversifies the roles of LINEs. 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 assembly of repressive RBPs, and contribute to evolution of new, lineage-specific transcripts in mammals.

Large cohorts of human iPSCs from healthy donors are potentially a powerful tool for investigating the relationship between genetic variants and cellular phenotypes. Here we integrate high content imaging, gene expression and DNA sequence datasets for over 100 human iPSC lines to identify the genetic basis of inter-individual variability in cell behaviour. By applying a dimensionality reduction approach, Probabilistic Estimation of Expression Residuals (PEER), we identified genes that correlated in expression with intrinsic (genetic) and extrinsic (ECM) factors. However, variation in mRNA levels could not account for outlier cell behaviour. Instead, we identified rare, deleterious SNVs in the coding sequence of genes involved in ECM adhesion that occurred in cell lines that were outliers for one or more phenotypes such as cell spreading. These also correlated with altered germ layer differentiation on micropatterned surfaces. Our study thus establishes a strategy for integrating genetic and cell biological measurements for high-throughput analysis.

The evolutionary features of clear-cell renal cell carcinoma (ccRCC) have not been systematically studied to date. We analyzed 1,206 primary tumor regions from 101 patients recruited into the multi-center prospective study, TRACERx Renal. We observe up to 30 driver events per tumor and show that subclonal diversification is associated with known prognostic parameters. By resolving the patterns of driver event ordering, co-occurrence, and mutual exclusivity at clone level, we show the deterministic nature of clonal evolution. ccRCC can be grouped into seven evolutionary subtypes, ranging from tumors characterized by early fixation of multiple mutational and copy number drivers and rapid metastases to highly branched tumors with > 10 subclonal drivers and extensive parallel evolution associated with attenuated progression. We identify genetic diversity and chromosomal complexity as determinants of patient outcome. Our insights reconcile the variable clinical behavior of ccRCC and suggest evolutionary potential as a biomarker for both intervention and surveillance. We thank Aida Murra, Naheed Shaikh, Justine Korteweg, Jeremy Tai, Eleanor Carlyle, Leonora Conneely, KimEdmonds, Karla Lingard, Karen O'Meara, Helen Breeze, Sarah Sarker, Lesley Cooper, Linda Shephard, Susie Slater, and Catherine Rogers for study support. We thank the patients and their families. S.T. and H.X. are funded by Cancer Research UK (CRUK) (C50947/A18176). S.T., T.C., J.L., and M.G. are funded by the National Institute for Health Research (NIHR) Biomedical Research Centre (BRC) at the Royal Marsden Hospital and Institute of Cancer Research (A109). J.I.L. is funded by the Ministerio de Economia y Competitividad (MINECO, SAF2016-79847-R). M.S., A.S., J. Lynch, R.F., L.A., and L.S. are funded by the Royal Marsden Cancer Charity. K.L. is funded by UK Medical Research Council (MR/P014712/1). T.B.K.W. is funded by the European Union Seventh Framework Programme (FP7-People-2013-ITN). M.J.-H. is funded by the NIHR. N.M.L. is a Winton Group Leader in recognition of the Winton Charitable Foundation's support towards the establishment of The Francis Crick Institute. N.M.L. is additionally funded by a Wellcome Trust Joint Investigator Award (103760/Z/14/Z) and the MRC eMedLab Medical Bioinformatics Infrastructure Award (MR/L016311/1). M.K. is funded by the Danish Cancer Society grant (R90-A6213). P.C. is funded by the Wellcome Trust (WT088340MA). N.M. receives funding from CRUK, Rosetrees, and the NIHR BRC at University College London Hospitals. C.S. is a Royal Society Napier Research Professor. C.S. is funded by Cancer Research UK (TRACERx and CRUK Cancer Immunotherapy Catalyst Network), the CRUK Lung Cancer Centre of Excellence, Stand Up 2 Cancer (SU2C), the Rosetrees and Stoneygate Trusts, NovoNordisk Foundation (THESEUS), Marie Curie Network PloidyNet, the NIHR BRC at University College London Hospitals, and the CRUK University College London Experimental Cancer Medicine Centre. C.S., O.P., D.R., I.C., and Z.S. are funded by NovoNordisk Foundation (16584). O.P., D.R., and I.C. are funded by the National Research, Development and Innovation Office of Hungary (NVKP_16-1-20160004). This work was supported by CRUK (Clinical Scientist Fellowship to S.T., C50947/A18176), The Francis Crick Institute, which receives its core funding from Cancer Reseach UK (FC010110), the UK Medical Research Council (FC010110), the Wellcome Trust (FC010110), and NIHR BRC at the Royal Marsden Hospital and Institute of Cancer Research (A109). High-performance computing was supported by eMedLab. The results published here are in whole or part based upon data generated by the TCGA Research Network: http://cancergenome.nih.gov/.

SummaryNeural induction in vertebrates generates a central nervous system that extends the rostral-caudal length of the body. The prevailing view is that neural cells are initially induced with anterior (forebrain) identity, with caudalising signals then converting a proportion to posterior fates (spinal cord). To test this model, we used chromatin accessibility assays to define how cells adopt region-specific neural fates. Together with genetic and biochemical perturbations this identified a developmental time window in which genome-wide chromatin remodeling events preconfigure epiblast cells for neural induction. Contrary to the established model, this revealed that cells commit to a regional identity before acquiring neural identity. This “primary regionalization” allocates cells to anterior or posterior regions of the nervous system, explaining how cranial and spinal neurons are generated at appropriate axial positions. These findings prompt a revision to models of neural induction and support the proposed dual evolutionary origin of the vertebrate central nervous system.

Summary Neural induction in vertebrates generates a CNS that extends the rostral-caudal length of the body. The prevailing view is that neural cells are initially induced with anterior (forebrain) identity; caudalizing signals then convert a proportion to posterior fates (spinal cord). To test this model, we used chromatin accessibility to define how cells adopt region-specific neural fates. Together with genetic and biochemical perturbations, this identified a developmental time window in which genome-wide chromatin-remodeling events preconfigure epiblast cells for neural induction. Contrary to the established model, this revealed that cells commit to a regional identity before acquiring neural identity. This “primary regionalization” allocates cells to anterior or posterior regions of the nervous system, explaining how cranial and spinal neurons are generated at appropriate axial positions. These findings prompt a revision to models of neural induction and support the proposed dual evolutionary origin of the vertebrate CNS., Graphical Abstract, Highlights • Chromatin accessibility defines neural progenitor identity • A limited developmental window exists to establish spinal cord competency • Cells acquire axial identity prior to neural identity, Regional identity precedes neural identity of the anterior and posterior nervous systems.

The structure of messenger RNA is important for post-transcriptional regulation, mainly because it affects binding of trans-acting factors. However, little is known about the in vivo structure of full-length mRNAs. Here we present hiCLIP, a biochemical technique for transcriptome-wide identification of RNA secondary structures interacting with RNA-binding proteins (RBPs). Using this technique to investigate RNA structures bound by Staufen 1 (STAU1) in human cells, we uncover a dominance of intra-molecular RNA duplexes, a depletion of duplexes from coding regions of highly translated mRNAs, an unexpected prevalence of long-range duplexes in 3' untranslated regions (UTRs), and a decreased incidence of single nucleotide polymorphisms in duplex-forming regions. We also discover a duplex spanning 858 nucleotides in the 3' UTR of the X-box binding protein 1 (XBP1) mRNA that regulates its cytoplasmic splicing and stability. Our study reveals the fundamental role of mRNA secondary structures in gene expression and introduces hiCLIP as a widely applicable method for discovering new, especially long-range, RNA duplexes.

Cellular morphology is an essential determinant of cellular function in all kingdoms of life, yet little is known about how cell shape is controlled. Here we describe a molecular program that controls the early morphology of neurons through a metazoan-specific zinc finger protein, Unkempt. Depletion of Unkempt in mouse embryos disrupts the shape of migrating neurons, while ectopic expression confers neuronal-like morphology to cells of different nonneuronal lineages. We found that Unkempt is a sequence-specific RNA-binding protein and identified its precise binding sites within coding regions of mRNAs linked to protein metabolism and trafficking. RNA binding is required for Unkempt-induced remodeling of cellular shape and is directly coupled to a reduced production of the encoded proteins. These findings link post-transcriptional regulation of gene expression with cellular shape and have general implications for the development and disease of multicellular organisms.

An interplay of experimental and computational methods is required to achieve a comprehensive understanding of protein-RNA interactions. Crosslinking and immunoprecipitation (CLIP) identifies endogenous interactions by sequencing RNA fragments that co-purify with a selected RBP under stringent conditions. Here we focus on approaches for the analysis of resulting data and appraise the methods for peak calling, visualisation, analysis and computational modelling of protein-RNA binding sites. We advocate a combined assessment of cDNA complexity and specificity for data quality control. Moreover, we demonstrate the value of analysing sequence motif enrichment in peaks assigned from CLIP data, and of visualising 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.

Background Canonical Wnt/beta-catenin signalling regulates self-renewal and lineage selection within the mammalian epidermis. Although the transcriptional response of keratinocytes that receive a Wnt signal is well characterized, little is known about the mechanism by which keratinocytes in proximity to the Wnt-receiving cell are co-opted to undergo a change in cell fate. Results To address this, we perform single-cell RNA-sequencing on mouse keratinocytes co-cultured with and without beta-catenin-activated neighbouring cells. We identify five distinct cell states in cultures that had not been exposed to the beta-catenin stimulus and show that the stimulus redistributes wild-type subpopulation proportions. Using temporal single-cell analysis, we reconstruct the cell fate change induced by Wnt activation from neighbouring cells. Gene expression heterogeneity is reduced in neighbouring cells and this effect is most dramatic for protein synthesis-associated genes. Changes in gene expression are accompanied by a shift to a more proliferative stem cell state. By integrating imaging and reconstructed sequential gene expression changes during the state transition we identify transcription factors, including Smad4 and Bcl3, that are responsible for effecting the transition in a contact-dependent manner. Conclusions Our data indicate that non-cell autonomous Wnt/beta-catenin signalling decreases transcriptional heterogeneity. This furthers our understanding of how epidermal Wnt signalling orchestrates regeneration and self-renewal. Electronic supplementary material The online version of this article (doi:10.1186/s13059-017-1384-y) contains supplementary material, which is available to authorized users.

Most genomes, including yeast Saccharomyces cerevisiae, are pervasively transcribed producing numerous non-coding RNAs, many of which are unstable and eliminated by nuclear or cytoplasmic surveillance pathways. We previously showed that accumulation of PHO84 antisense RNA (asRNA), in cells lacking the nuclear exosome component Rrp6, is paralleled by repression of sense transcription in a process dependent on the Hda1 histone deacetylase (HDAC) and the H3K4 histone methyl transferase Set1. Here we investigate this process genome-wide and measure the whole transcriptome of various histone modification mutants in a Δrrp6 strain using tiling arrays. We confirm widespread occurrence of potentially antisense-dependent gene regulation and identify three functionally distinct classes of genes that accumulate asRNAs in the absence of Rrp6. These classes differ in whether the genes are silenced by the asRNA and whether the silencing is HDACs and histone methyl transferase-dependent. Among the distinguishing features of asRNAs with regulatory potential, we identify weak early termination by Nrd1/Nab3/Sen1, extension of the asRNA into the open reading frame promoter and dependence of the silencing capacity on Set1 and the HDACs Hda1 and Rpd3 particularly at promoters undergoing extensive chromatin remodelling. Finally, depending on the efficiency of Nrd1/Nab3/Sen1 early termination, asRNA levels are modulated and their capability of silencing is changed.

In the version of this article initially published, the Supplementary Data file was an incorrect version. The correct version is now provided. The error has been corrected in the HTML and PDF version of the article.

Background DNA is subject to constant chemical modification and damage, which eventually results in variable mutation rates throughout the genome. Although detailed molecular mechanisms of DNA damage and repair are well understood, damage impact and execution of repair across a genome remain poorly defined. Results To bridge the gap between our understanding of DNA repair and mutation distributions, we developed a novel method, AP-seq, capable of mapping apurinic sites and 8-oxo-7,8-dihydroguanine bases at approximately 250-bp resolution on a genome-wide scale. We directly demonstrate that the accumulation rate of apurinic sites varies widely across the genome, with hot spots acquiring many times more damage than cold spots. Unlike single nucleotide variants (SNVs) in cancers, damage burden correlates with marks for open chromatin notably H3K9ac and H3K4me2. Apurinic sites and oxidative damage are also highly enriched in transposable elements and other repetitive sequences. In contrast, we observe a reduction at chromatin loop anchors with increased damage load towards inactive compartments. Less damage is found at promoters, exons, and termination sites, but not introns, in a seemingly transcription-independent but GC content-dependent manner. Leveraging cancer genomic data, we also find locally reduced SNV rates in promoters, coding sequence, and other functional elements. Conclusions Our study reveals that oxidative DNA damage accumulation and repair differ strongly across the genome, but culminate in a previously unappreciated mechanism that safeguards the regulatory and coding regions of genes from mutations. Electronic supplementary material The online version of this article (10.1186/s13059-018-1582-2) contains supplementary material, which is available to authorized users.