Search

Oikopleura dioica is a planktonic tunicate (Appendicularia class) found extensively across the marine waters of the globe. The genome of a single male individual collected from Okinawa, Japan was sequenced using the single-molecule PacBio Hi-Fi method and assembled with NOVOLoci. The mitogenome is 39,268 bp long, featuring a large control region of around 22,000 bp. We annotated the proteins atp6, cob, cox1, cox2, cox3, nad1, nad4, and nad5, and found one more open reading frame that did not match any known gene. This study marks the first complete mitogenome assembly for an appendicularian, and reveals that A and T homopolymers cumulatively account for nearly half of its length. This reference sequence will be an asset for environmental DNA and phylogenetic studies., (© The Author(s) 2024. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Signaling pathways drive cell fate transitions largely by changing gene expression. However, the mechanisms for rapid and selective transcriptome rewiring in response to signaling cues remain elusive. Here we use deep learning to deconvolve both the sequence determinants and the trans-acting regulators that trigger extracellular signal-regulated kinase (ERK)-mitogen-activated protein kinase kinase (MEK)-induced decay of the naive pluripotency mRNAs. Timing of decay is coupled to embryo implantation through ERK-MEK phosphorylation of LIN28A, which repositions pLIN28A to the highly A+U-rich 3' untranslated region (3'UTR) termini of naive pluripotency mRNAs. Interestingly, these A+U-rich 3'UTR termini serve as poly(A)-binding protein (PABP)-binding hubs, poised for signal-induced convergence with LIN28A. The multivalency of AUU motifs determines the efficacy of pLIN28A-PABP convergence, which enhances PABP 3'UTR binding, decreases the protection of poly(A) tails and activates mRNA decay to enable progression toward primed pluripotency. Thus, the signal-induced convergence of LIN28A with PABP-RNA hubs drives the rapid selection of naive mRNAs for decay, enabling the transcriptome remodeling that ensures swift developmental progression., (© 2024. The Author(s).)

Genome structural variations within species are rare. How selective constraints preserve gene order and chromosome structure is a central question in evolutionary biology that remains unsolved. Our sequencing of several genomes of the appendicularian tunicate Oikopleura dioica around the globe reveals extreme genome scrambling caused by thousands of chromosomal rearrangements, although showing no obvious morphological differences between these animals. The breakpoint accumulation rate is an order of magnitude higher than in ascidian tunicates, nematodes, Drosophila, or mammals. Chromosome arms and sex-specific regions appear to be the primary unit of macrosynteny conservation. At the microsyntenic level, scrambling did not preserve operon structures, suggesting an absence of selective pressure to maintain them. The uncoupling of the genome scrambling with morphological conservation in O. dioica suggests the presence of previously unnoticed cryptic species and provides a new biological system that challenges our previous vision of speciation in which similar animals always share similar genome structures., (© 2024 Plessy et al.; Published by Cold Spring Harbor Laboratory Press.)

Background: CD4 + Th1 cells producing IFN-γ are required to eradicate intracellular pathogens, however if uncontrolled these cells can cause immunopathology. The cytokine IL-10 is produced by multiple immune cells including Th1 cells during infection and regulates the immune response to minimise collateral host damage. In this study we aimed to elucidate the transcriptional network of genes controlling the expression of Il10 and proinflammatory cytokines, including Ifng in Th1 cells differentiated from mouse naive CD4 + T cells., Methods: We applied computational analysis of gene regulation derived from temporal profiling of gene expression clusters obtained from bulk RNA sequencing (RNA-seq) of flow cytometry sorted naïve CD4 + T cells from mouse spleens differentiated in vitro into Th1 effector cells with IL-12 and IL-27 to produce Ifng and Il10, compared to IL-27 alone which express Il10 only , or IL-12 alone which express Ifng and no Il10, or medium control driven-CD4 + T cells which do not express effector cytokines . Data were integrated with analysis of active genomic regions from these T cells using an assay for transposase-accessible chromatin with sequencing (ATAC)-seq, integrated with literature derived-Chromatin-immunoprecipitation (ChIP)-seq data and the RNA-seq data, to elucidate the transcriptional network of genes controlling expression of Il10 and pro-inflammatory effector genes in Th1 cells. The co-dominant role for the transcription factors, Prdm1 (encoding Blimp-1) and Maf (encoding c-Maf) , in cytokine gene regulation in Th1 cells, was confirmed using T cells obtained from mice with T-cell specific deletion of these transcription factors., Results: We show that the transcription factors Blimp-1 and c-Maf each have unique and common effects on cytokine gene regulation and not only co-operate to induce Il10 gene expression in IL-12 plus IL-27 differentiated mouse Th1 cells, but additionally directly negatively regulate key proinflammatory cytokines including Ifng , thus providing mechanisms for reinforcement of regulated Th1 cell responses., Conclusions: These data show that Blimp-1 and c-Maf positively and negatively regulate a network of both unique and common anti-inflammatory and pro-inflammatory genes to reinforce a Th1 response in mice that will eradicate pathogens with minimum immunopathology., Competing Interests: No competing interests were disclosed., (Copyright: © 2023 Cox LS et al.)

The vertebrate 'neural plate border' is a transient territory located at the edge of the neural plate containing precursors for all ectodermal derivatives: the neural plate, neural crest, placodes and epidermis. Elegant functional experiments in a range of vertebrate models have provided an in-depth understanding of gene regulatory interactions within the ectoderm. However, these experiments conducted at tissue level raise seemingly contradictory models for fate allocation of individual cells. Here, we carry out single cell RNA sequencing of chick ectoderm from primitive streak to neurulation stage, to explore cell state diversity and heterogeneity. We characterise the dynamics of gene modules, allowing us to model the order of molecular events which take place as ectodermal fates segregate. Furthermore, we find that genes previously classified as neural plate border 'specifiers' typically exhibit dynamic expression patterns and are enriched in either neural, neural crest or placodal fates, revealing that the neural plate border should be seen as a heterogeneous ectodermal territory and not a discrete transitional transcriptional state. Analysis of neural, neural crest and placodal markers reveals that individual NPB cells co-express competing transcriptional programmes suggesting that their ultimate identify is not yet fixed. This population of 'border located undecided progenitors' (BLUPs) gradually diminishes as cell fate decisions take place. Considering our findings, we propose a probabilistic model for cell fate choice at the neural plate border. Our data suggest that the probability of a progenitor's daughters to contribute to a given ectodermal derivative is related to the balance of competing transcriptional programmes, which in turn are regulated by the spatiotemporal position of a progenitor., Competing Interests: AT, AB, EH, CC, NL, JB, AS No competing interests declared, (© 2023, Thiery et al.)

Crosslinking and immunoprecipitation (CLIP) technologies have become a central component of the molecular biologists' toolkit to study protein-RNA interactions and thus to uncover core principles of RNA biology. There has been a proliferation of CLIP-based experimental protocols, as well as computational tools, especially for peak-calling. Consequently, there is an urgent need for a well-documented bioinformatic pipeline that enshrines the principles of robustness, reproducibility, scalability, portability and flexibility while embracing the diversity of experimental and computational CLIP tools. To address this, we present nf-core/clipseq - a robust Nextflow pipeline for quality control and analysis of CLIP sequencing data. It is part of the international nf-core community effort to develop and curate a best-practice, gold-standard set of pipelines for data analysis. The standards enabled by Nextflow and nf-core, including workflow management, version control, continuous integration and containerisation ensure that these key needs are met. Furthermore, multiple tools are implemented ( e.g. for peak-calling), alongside visualisation of quality control metrics to empower the user to make their own informed decisions based on their data. nf-core/clipseq remains under active development, with plans to incorporate newly released tools to ensure that pipeline remains up-to-date and relevant for the community. Engagement with users and developers is encouraged through the nf-core GitHub repository and Slack channel to promote collaboration. It is available at https://nf-co.re/clipseq., Competing Interests: Competing interests: CEMW, CCh, NML, AMC and JU have no conflicts of interest to declare. CCa is employed part-time by the software company Goodwright (https://goodwright.com) to consult on their bioinformatics workflow web platform Flow (flow.bio)., (Copyright: © 2023 West C et al.)

CLIP technologies are now widely used to study RNA-protein interactions and many data sets 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 visualize a 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 normalization 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., (© 2023 Chakrabarti et al.; Published by Cold Spring Harbor Laboratory Press for the RNA Society.)

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., (© The Author(s) 2023. Published by Oxford University Press on behalf of the Guarantors of Brain. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

The structure of mRNA molecules plays an important role in its interactions with trans-acting factors, notably RNA binding proteins (RBPs), thus contributing to the functional consequences of this interplay. However, current transcriptome-wide experimental methods to chart these interactions are limited by their poor sensitivity. Here we extend the hiCLIP atlas of duplexes bound by Staufen1 (STAU1) ∼10-fold, through careful consideration of experimental assumptions, and the development of bespoke computational methods which we apply to existing data. We present Tosca, a Nextflow computational pipeline for the processing, analysis and visualisation of proximity ligation sequencing data generally. We use our extended duplex atlas to discover insights into the RNA selectivity of STAU1, revealing the importance of structural symmetry and duplex-span-dependent nucleotide composition. Furthermore, we identify heterogeneity in the relationship between transcripts with STAU1-bound 3' UTR duplexes and metabolism of the associated RNAs that we relate to RNA structure: transcripts with short-range proximal 3' UTR duplexes have high degradation rates, but those with long-range duplexes have low rates. Overall, our work enables the integrative analysis of proximity ligation data delivering insights into specific features and effects of RBP-RNA structure interactions., (© The Author(s) 2023. Published by Oxford University Press on behalf of Nucleic Acids Research.)

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., Competing Interests: KT, HL, YY, AT, AS, CA, BP, Nd, Id, MK, NM, NL, LD, AS, CS No competing interests declared, (© 2023, Trevers, Lu et al.)

Indirect development with an intermediate larva exists in all major animal lineages 1 , which makes larvae central to most scenarios of animal evolution 2-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 head 12 . Notably, this also occurs in the so-called 'head larvae' of other bilaterians 13-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-option 9,10 or innovation 11 of gene regulatory programmes to explain larva and adult origins., (© 2023. The Author(s).)

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., Competing Interests: RV, TS, CC, ZM, EC, AR, RE, IE, VC, HP, JK, PF, AS, MV, MR, JU, NL, Fv No competing interests declared, (© 2022, Varier, Sideri, Capitanchik et al.)

Background : Automation has increasingly become more commonplace in the research laboratory workspace. The introduction of articulated robotic arms allows the researcher more flexibility in the tasks a single piece of automated machinery can perform. We set out to incorporate automation in processing of genomic DNA organic extractions to increase throughput and limit researchers to the exposure of organic solvents. Methods : In order to automate the genome sequencing pipeline in our laboratory, we programmed a dual-arm anthropomorphic robot, the Robotic Biology Institute's Maholo LabDroid, to perform organic solvent-based genomic DNA extraction from cell lysates. To the best of our knowledge, this is the first time that automation of phenol-chloroform extraction has been reported. Results: We achieved routine extraction of high molecular weight genomic DNA (>100 kb) from diverse biological samples including algae cultured in sea water, bacteria, whole insects, and human cell lines. The results of pulse-field electrophoresis size analysis and the N50 sequencing metrics of reads obtained from Nanopore MinION runs verified the presence of intact DNA suitable for direct sequencing. Conclusions : We present the workflow that can be used to program similar robots and discuss the problems and solutions we encountered in developing the workflow. The protocol can be adapted to analogous methods such as RNA extraction, and there is ongoing work to incorporate further post-extraction steps such as library construction. This work shows the potential for automated robotic workflows to free molecular biological researchers from manual interventions in routine experimental work. A time-lapse movie of the entire automated run is included in this report., Competing Interests: No competing interests were disclosed., (Copyright: © 2022 Liu AW et al.)

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 SOD1 G93A 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., (© 2022 Ziff et al.; Published by Cold Spring Harbor Laboratory Press.)

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., (© The Author(s) 2021. Published by Oxford University Press on behalf of Nucleic Acids Research.)

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 m 6 A 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., (© 2021. The Author(s).)

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.

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., Competing Interests: Declaration of interests J.S. is a consultant for Dewpoint Therapeutics, Maze Therapeutics, and Vivid Sciences. B.P. is an employee of Dewpoint Therapeutics., (Copyright © 2021 The Authors. Published by Elsevier Inc. All rights reserved.)

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., (© The Author(s) (2021). Published by Oxford University Press on behalf of the Guarantors of Brain.)

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., Competing Interests: No competing interests were disclosed., (Copyright: © 2021 Wilkins OG et al.)

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.

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., (© The Author(s) 2021. Published by Oxford University Press on behalf of Nucleic Acids Research.)

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.

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.

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., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2021 The Author(s). Published by Elsevier Inc. All rights reserved.)

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., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2020 The Author(s). Published by Elsevier Inc. All rights reserved.)

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.

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 as BCL9, MCL1, ARNT (also known as HIF1B), TERT and MYC) 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 encompasses BCL9, MCL1 and ARNT), 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 (encompassing MYC) in clear cell renal cell carcinoma and chromosome 11q13.3 (encompassing CCND1) 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.

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., (© 2020 The Authors. British Journal of Dermatology published by John Wiley & Sons Ltd on behalf of British Association of Dermatologists.)

The larvacean Oikopleura dioica is a planktonic chordate and is a tunicate that belongs to the closest relatives to vertebrates. Its simple and transparent body, invariant embryonic cell lineages, and short life cycle of 5 days make it a promising model organism for the study of developmental biology. The genome browser OikoBase was established in 2013 using Norwegian O. dioica. However, genome information for other populations is not available, even though many researchers have studied local populations. In the present study, we sequenced using Illumina and PacBio RSII technologies the genome of O. dioica from a southwestern Japanese population that was cultured in our laboratory for 3 years. The genome of Japanese O. dioica was assembled into 576 scaffold sequences with a total length and N50 length of 56.6 and 1.5 Mb, respectively. A total of 18,743 gene models (transcript models) were predicted in the genome assembly, named OSKA2016. In addition, 19,277 non-redundant transcripts were assembled using RNA-seq data. The OSKA2016 has global sequence similarity of only 86.5% when compared with the OikoBase, highlighting the sequence difference between the two far distant O. dioica populations on the globe. The genome assembly, transcript assembly, and transcript models were incorporated into ANISEED (https://www.aniseed.cnrs.fr/) for genome browsing and BLAST searches. Mapping of reads obtained from male- or female-specific genome libraries yielded male-specific scaffolds in the OSKA2016 and revealed that over 2.6 Mb of sequence were included in the male-specific Y-region. The genome and transcriptome resources from two distinct populations will be useful datasets for developmental biology, evolutionary biology, and molecular ecology using this model organism., (© 2020 Japanese Society of Developmental Biologists.)

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., Competing Interests: No competing interests were disclosed., (Copyright: © 2020 Liu AW et al.)

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., Competing Interests: No competing interests were disclosed., (Copyright: © 2021 Liu AW et al.)

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.

A flurry of methods has been developed in recent years to identify N6-methyladenosine (m 6 A) sites across transcriptomes at high resolution. This raises the need to understand both the common features and those that are unique to each method. Here, we complement the analyses presented in the original papers by reviewing their various technical aspects and comparing the overlap between m 6 A-methylated messenger RNAs (mRNAs) identified by each. Specifically, we examine eight different methods that identify m 6 A sites in human cells with high resolution: two antibody-based crosslinking and immunoprecipitation (CLIP) approaches, two using endoribonuclease MazF, one based on deamination, two using Nanopore direct RNA sequencing, and finally, one based on computational predictions. We contrast the respective datasets and discuss the challenges in interpreting the overlap between them, including a prominent expression bias in detected genes. This overview will help guide researchers in making informed choices about using the available data and assist with the design of future experiments to expand our understanding of m 6 A and its regulation., (Copyright © 2020 Capitanchik, Toolan-Kerr, Luscombe and Ule.)

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.

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.