Search

Abstract RNA-binding proteins (RBPs) have pivotal functions in RNA metabolism, but current methods are limited in retrieving RBP-RNA interactions within endogenous biological contexts. Here, we develop INSCRIBE (IN situ Sensitive Capture of RNA-protein Interactions in Biological Environments), circumventing the challenges through in situ RNA labeling by precisely directing a purified APOBEC1-nanobody fusion to the RBP of interest. This method enables highly specific RNA-binding site identification across a diverse range of fixed biological samples such as HEK293T cells and mouse brain tissue and accurately identifies the canonical binding motifs of RBFOX2 (UGCAUG) and TDP-43 (UGUGUG) in native cellular environments. Applicable to any RBP with available primary antibodies, INSCRIBE enables sensitive capture of RBP-RNA interactions from ultra-low input equivalent to ~5 cells. The robust, versatile, and sensitive INSCRIBE workflow is particularly beneficial for precious tissues such as clinical samples, empowering the exploration of genuine RBP-RNA interactions in RNA-related disease contexts.

Abstract RNA-binding proteins (RBPs) regulate key aspects of RNA processing including alternative splicing, mRNA degradation and localization by physically binding RNA molecules. Current methods to map these interactions, such as CLIP, rely on purifying single proteins at a time. Our new method, ePRINT, maps RBP-RNA interaction networks on a global scale without purifying individual RBPs. ePRINT uses exoribonuclease XRN1 to precisely map the 5′ end of the RBP binding site and uncovers direct and indirect targets of an RBP of interest. Importantly, ePRINT can also uncover RBPs that are differentially activated between cell fate transitions, including neural progenitor differentiation into neurons.

Abstract The spatial organization of molecules in a cell is essential for their functions. While current methods focus on discerning tissue architecture, cell–cell interactions, and spatial expression patterns, they are limited to the multicellular scale. We present Bento, a Python toolkit that takes advantage of single-molecule information to enable spatial analysis at the subcellular scale. Bento ingests molecular coordinates and segmentation boundaries to perform three analyses: defining subcellular domains, annotating localization patterns, and quantifying gene–gene colocalization. We demonstrate MERFISH, seqFISH + , Molecular Cartography, and Xenium datasets. Bento is part of the open-source Scverse ecosystem, enabling integration with other single-cell analysis tools.

Summary: Reactive changes of glial cells during neuroinflammation impact brain disorders and disease progression. Elucidating the mechanisms that control reactive gliosis may help us to understand brain pathophysiology and improve outcomes. Here, we report that adult ablation of autism spectrum disorder (ASD)-associated CHD8 in astrocytes attenuates reactive gliosis via remodeling chromatin accessibility, changing gene expression. Conditional Chd8 deletion in astrocytes, but not microglia, suppresses reactive gliosis by impeding astrocyte proliferation and morphological elaboration. Astrocyte Chd8 ablation alleviates lipopolysaccharide-induced neuroinflammation and septic-associated hypothermia in mice. Astrocytic CHD8 plays an important role in neuroinflammation by altering the chromatin landscape, regulating metabolic and lipid-associated pathways, and astrocyte-microglia crosstalk. Moreover, we show that reactive gliosis can be directly mitigated in vivo using an adeno-associated virus (AAV)-mediated Chd8 gene editing strategy. These findings uncover a role of ASD-associated CHD8 in the adult brain, which may warrant future exploration of targeting chromatin remodelers in reactive gliosis and neuroinflammation in injury and neurological diseases.

Abstract RNA binding proteins (RBPs) are key regulators of RNA processing and cellular function. Technologies to discover RNA targets of RBPs such as TRIBE (targets of RNA binding proteins identified by editing) and STAMP (surveying targets by APOBEC1 mediated profiling) utilize fusions of RNA base-editors (rBEs) to RBPs to circumvent the limitations of immunoprecipitation (CLIP)-based methods that require enzymatic digestion and large amounts of input material. To broaden the repertoire of rBEs suitable for editing-based RBP-RNA interaction studies, we have devised experimental and computational assays in a framework called PRINTER (protein-RNA interaction-based triaging of enzymes that edit RNA) to assess over thirty A-to-I and C-to-U rBEs, allowing us to identify rBEs that expand the characterization of binding patterns for both sequence-specific and broad-binding RBPs. We also propose specific rBEs suitable for dual-RBP applications. We show that the choice between single or multiple rBEs to fuse with a given RBP or pair of RBPs hinges on the editing biases of the rBEs and the binding preferences of the RBPs themselves. We believe our study streamlines and enhances the selection of rBEs for the next generation of RBP-RNA target discovery.

Summary: Here, we present a protocol for using Skipper, a pipeline designed to process crosslinking and immunoprecipitation (CLIP) data into annotated binding sites. We describe steps for partitioning annotated transcript regions and fitting data to a beta-binomial model to call windows of enriched binding. From raw CLIP data, we detail how users can map reproducible RNA-binding sites to call enriched windows and perform downstream analysis. This protocol supports optional customizations for different use cases.For complete details on the use and execution of this protocol, please refer to Boyle et al.1 : Publisher’s note: Undertaking any experimental protocol requires adherence to local institutional guidelines for laboratory safety and ethics.

Summary: The RNA-binding protein PARP13 is a primary factor in the innate antiviral response, which suppresses translation and drives decay of bound viral and host RNA. PARP13 interacts with many proteins encoded by interferon-stimulated genes (ISG) to activate antiviral pathways including co-translational addition of ISG15, or ISGylation. We performed enhanced crosslinking immunoprecipitation (eCLIP) and RNA-seq in human cells to investigate PARP13’s role in transcriptome regulation for both basal and antiviral states. We find that the antiviral response shifts PARP13 target localization, but not its binding preferences, and that PARP13 supports the expression of ISGylation-related genes, including PARP13’s cofactor, TRIM25. PARP13 associates with TRIM25 via RNA-protein interactions, and we elucidate a transcriptome-wide periodicity of PARP13 binding around TRIM25. Taken together, our study implicates PARP13 in creating and maintaining a cellular environment poised for an antiviral response through limiting PARP13 translation, regulating access to distinct mRNA pools, and elevating ISGylation machinery expression.

Abstract Background Fusion of RNA-binding proteins (RBPs) to RNA base-editing enzymes (such as APOBEC1 or ADAR) has emerged as a powerful tool for the discovery of RBP binding sites. However, current methods that analyze sequencing data from RNA-base editing experiments are vulnerable to false positives due to off-target editing, genetic variation and sequencing errors. Results We present FLagging Areas of RNA-editing Enrichment (FLARE), a Snakemake-based pipeline that builds on the outputs of the SAILOR edit site discovery tool to identify regions statistically enriched for RNA editing. FLARE can be configured to analyze any type of RNA editing, including C to U and A to I. We applied FLARE to C-to-U editing data from a RBFOX2-APOBEC1 STAMP experiment, to show that our approach attains high specificity for detecting RBFOX2 binding sites. We also applied FLARE to detect regions of exogenously introduced as well as endogenous A-to-I editing. Conclusions FLARE is a fast and flexible workflow that identifies significantly edited regions from RNA-seq data. The FLARE codebase is available at https://github.com/YeoLab/FLARE .

Abstract Pancreatic ductal adenocarcinoma (PDAC) is a lethal malignancy in need of new therapeutic options. Using unbiased analyses of super-enhancers (SEs) as sentinels of core genes involved in cell-specific function, here we uncover a druggable SE-mediated RNA-binding protein (RBP) cascade that supports PDAC growth through enhanced mRNA translation. This cascade is driven by a SE associated with the RBP heterogeneous nuclear ribonucleoprotein F, which stabilizes protein arginine methyltransferase 1 (PRMT1) to, in turn, control the translational mediator ubiquitin-associated protein 2-like. All three of these genes and the regulatory SE are essential for PDAC growth and coordinately regulated by the Myc oncogene. In line with this, modulation of the RBP network by PRMT1 inhibition reveals a unique vulnerability in Myc-high PDAC patient organoids and markedly reduces tumor growth in male mice. Our study highlights a functional link between epigenetic regulation and mRNA translation and identifies components that comprise unexpected therapeutic targets for PDAC.

Nonalcoholic steatohepatitis (NASH) is an advanced form of fatty liver disease with complex pathogenic mechanisms. Here, the authors report that SRSF1 deficiency in mice livers provokes deleterious R-loop formation and genotoxicity, which impedes hepatocellular gene expression, metabolism, and lipid trafficking, resulting in NASH-like pathology.

MYC amplification is an independent prognostic factor for the most aggressive subgroup (Group 3) of pediatric medulloblastoma (G3 MB). Here, the authors highlight the role of the RNA-binding protein, Musashi-1 (MSI1) in G3 MB and identify MSI1-bound targets sharing MYC associated pathways.

Leukemic stem cells (LSCs) drive chemoresistance and relapse in acute myeloid leukemia. Here, the authors show that the splicing factor RBM17 supports LSCs through avoiding nonsense-mediated decay of pro-leukaemic factors such as the translation initiation factor EIF4A2.

Summary: Background: Schools are high-risk settings for SARS-CoV-2 transmission, but necessary for children's educational and social-emotional wellbeing. Previous research suggests that wastewater monitoring can detect SARS-CoV-2 infections in controlled residential settings with high levels of accuracy. However, its effective accuracy, cost, and feasibility in non-residential community settings is unknown. Methods: The objective of this study was to determine the effectiveness and accuracy of community-based passive wastewater and surface (environmental) surveillance to detect SARS-CoV-2 infection in neighborhood schools compared to weekly diagnostic (PCR) testing. We implemented an environmental surveillance system in nine elementary schools with 1700 regularly present staff and students in southern California. The system was validated from November 2020 to March 2021. Findings: In 447 data collection days across the nine sites 89 individuals tested positive for COVID-19, and SARS-CoV-2 was detected in 374 surface samples and 133 wastewater samples. Ninety-three percent of identified cases were associated with an environmental sample (95% CI: 88%–98%); 67% were associated with a positive wastewater sample (95% CI: 57%–77%), and 40% were associated with a positive surface sample (95% CI: 29%–52%). The techniques we utilized allowed for near-complete genomic sequencing of wastewater and surface samples. Interpretation: Passive environmental surveillance can detect the presence of COVID-19 cases in non-residential community school settings with a high degree of accuracy. Funding: County of San Diego, Health and Human Services Agency, National Institutes of Health, National Science Foundation, Centers for Disease Control.

During a microbial infection, responding CD8+ T cells give rise to effector cells that provide acute host defense and memory cells that provide sustained protection. An alternative outcome is exhaustion, a state of T cell dysfunction that occurs in the context of chronic infections and cancer. Although it is evident that exhausted CD8+ T (TEX) cells are phenotypically and molecularly distinct from effector and memory CD8+ T cells, the factors regulating the earliest events in the differentiation process of TEX cells remain incompletely understood. Here, we performed single-cell RNA-sequencing and single-cell ATAC-sequencing of CD8+ T cells responding to LCMV-Armstrong (LCMV-Arm) or LCMV-Clone 13 (LCMV-Cl13), which result in acute or chronic infections, respectively. Compared to CD8+ T cells that had undergone their first division in response to LCMV-Arm (Div1ARM) cells, CD8+ T cells that had undergone their first division in response to LCMV-Cl13 (Div1CL13) expressed higher levels of genes encoding transcription factors previously associated with exhaustion, along with higher levels of Ezh2, the catalytic component of the Polycomb Repressive Complex 2 (PRC2) complex, which mediates epigenetic silencing. Modulation of Ezh2 resulted in altered expression of exhaustion-associated molecules by CD8+ T cells responding to LCMV-Cl13, though the specific cellular and infectious contexts, rather than simply the level of Ezh2 expression, likely determine the eventual outcome. Taken together, these findings suggest that the differentiation paths of CD8+ T cells responding to acute versus chronic infections may diverge earlier than previously appreciated. How early can the process of T cell exhaustion begin? This study shows that CD8+ T cells that have undergone their first cellular division in response to chronic versus acute infection may have already diverged along transcriptionally and epigenetically distinct differentiation paths.

The off-target effects of CRISPR-Cas9 are thought to be mediated by its cognate guide RNA. Here the authors show that Cas9 independently interacts with the human transcriptome, correlating with elevated RNA editing even under guide RNA co-expression.

Abstract Throughout the COVID-19 pandemic, massive sequencing and data sharing efforts enabled the real-time surveillance of novel SARS-CoV-2 strains throughout the world, the results of which provided public health officials with actionable information to prevent the spread of the virus. However, with great sequencing comes great computation, and while cloud computing platforms bring high-performance computing directly into the hands of all who seek it, optimal design and configuration of a cloud compute cluster requires significant system administration expertise. We developed ViReflow, a user-friendly viral consensus sequence reconstruction pipeline enabling rapid analysis of viral sequence datasets leveraging Amazon Web Services (AWS) cloud compute resources and the Reflow system. ViReflow was developed specifically in response to the COVID-19 pandemic, but it is general to any viral pathogen. Importantly, when utilized with sufficient compute resources, ViReflow can trim, map, call variants, and call consensus sequences from amplicon sequence data from 1000 SARS-CoV-2 samples at 1000X depth in

Background & Aims: Hyperbaric oxygen therapy (HBOT) is a promising treatment for moderate-to-severe ulcerative colitis. However, our current understanding of the host and microbial response to HBOT remains unclear. This study examined the molecular mechanisms underpinning HBOT using a multi-omic strategy. Methods: Pre- and post-intervention mucosal biopsies, tissue, and fecal samples were collected from HBOT phase 2 clinical trials. Biopsies and fecal samples were subjected to shotgun metaproteomics, metabolomics, 16s rRNA sequencing, and metagenomics. Tissue was subjected to bulk RNA sequencing and digital spatial profiling (DSP) for single-cell RNA and protein analysis, and immunohistochemistry was performed. Fecal samples were also used for colonization experiments in IL10-/- germ-free UC mouse models. Results: Proteomics identified negative associations between HBOT response and neutrophil azurophilic granule abundance. DSP identified an HBOT-specific reduction of neutrophil STAT3, which was confirmed by immunohistochemistry. HBOT decreased microbial diversity with a proportional increase in Firmicutes and a secondary bile acid lithocholic acid. A major source of the reduction in diversity was the loss of mucus-adherent taxa, resulting in increased MUC2 levels post-HBOT. Targeted database searching revealed strain-level associations between Akkermansia muciniphila and HBOT response status. Colonization of IL10-/- with stool obtained from HBOT responders resulted in lower colitis activity compared with non-responders, with no differences in STAT3 expression, suggesting complementary but independent host and microbial responses. Conclusions: HBOT reduces host neutrophil STAT3 and azurophilic granule activity in UC patients and changes in microbial composition and metabolism in ways that improve colitis activity. Intestinal microbiota, especially strain level variations in A muciniphila, may contribute to HBOT non-response.

Mutations in the splicing factor RBM20 cause aggressive Dilated Cardiomyopathy. Here the authors generated RBM20 R636S mutants and knockout in human iPSC-derived cardiomyocytes. Mutant RBM20 showed different target RNA binding, altered splicing and localization to cytoplasmic processing bodies.

Synapse dysfunction and loss are central features of neurodegenerative diseases, caused in part by the accumulation of protein oligomers. Amyloid-β, tau, prion, and α-synuclein oligomers bind to the cellular prion protein (PrPC), resulting in the activation of macromolecular complexes and signaling at the post-synapse, yet the early signaling events are unclear. Here we sought to determine the early transcript and protein alterations in the hippocampus during the pre-clinical stages of prion disease. We used a transcriptomic approach focused on the early-stage, prion-infected hippocampus of male wild-type mice, and identify immediate early genes, including the synaptic activity response gene, Arc/Arg3.1, as significantly upregulated. In a longitudinal study of male, prion-infected mice, Arc/Arg-3.1 protein was increased early (40% of the incubation period), and by mid-disease (pre-clinical), phosphorylated AMPA receptors (pGluA1-S845) were increased and metabotropic glutamate receptors (mGluR5 dimers) were markedly reduced in the hippocampus. Notably, sporadic Creutzfeldt-Jakob disease (sCJD) post-mortem cortical samples also showed low levels of mGluR5 dimers. Together, these findings suggest that prions trigger an early Arc response, followed by an increase in phosphorylated GluA1 and a reduction in mGluR5 receptors.

Disordered RNA-binding proteins and repetitive RNA sequences are the main genetic causes of several neurodegenerative diseases, including amyotrophic lateral sclerosis and Huntington’s disease. Importantly, these components also seed the formation of cytoplasmic liquid-like granules, like stress granules and P bodies. Emerging evidence demonstrates that healthy granules formed via liquid-liquid phase separation can mature into solid- or gel-like inclusions that persist within the cell. These solidified inclusions are a precursor to the aggregates identified in patients, demonstrating that dysregulation of RNA granule biology is an important component of neurodegeneration. Here, we review recent literature highlighting how RNA molecules seed proteinaceous granules, the mechanisms of healthy turnover of RNA granules in cells, which biophysical properties underly a transition to solid- or gel-like material states, and why persistent granules disrupt the cellular homeostasis of neurons. We also identify various methods that will illuminate the contributions of disordered proteins and RNAs to neurodegeneration in ongoing research efforts.

Lipid transfer proteins mediate the exchange of lipids between closely apposed membranes at organelle contact sites and play key roles in lipid metabolism, membrane homeostasis, and cellular signaling. A recently discovered novel family of lipid transfer proteins, which includes the VPS13 proteins (VPS13A-D), adopt a rod-like bridge conformation with an extended hydrophobic groove that enables the bulk transfer of membrane lipids for membrane growth. Loss of function mutations in VPS13A and VPS13C cause chorea acanthocytosis and Parkinson's disease, respectively. VPS13A and VPS13C localize to multiple organelle contact sites, including endoplasmic reticulum (ER) – lipid droplet (LD) contact sites, but the functional roles of these proteins in LD regulation remains mostly unexplored. Here we employ CRISPR-Cas9 genome editing to generate VPS13A and VPS13C knockout cell lines in U-2 OS cells via deletion of exon 2 and introduction of an early frameshift. Analysis of LD content in these cell lines revealed that loss of either VPS13A or VPS13C results in reduced LD abundance under oleate-stimulated conditions. These data implicate two lipid transfer proteins, VPS13A and VPS13C, in LD regulation.

ABSTRACT A promising approach to help students safely return to in person learning is through the application of sentinel cards for accurate high resolution environmental monitoring of SARS-CoV-2 traces indoors. Because SARS-CoV-2 RNA can persist for up to a week on several indoor surface materials, there is a need for increased temporal resolution to determine whether consecutive surface positives arise from new infection events or continue to report past events. Cleaning sentinel cards after sampling would provide the needed resolution but might interfere with assay performance. We tested the effect of three cleaning solutions (BZK wipes, Wet Wipes, RNase Away) at three different viral loads: “high” (4 × 104 GE/mL), “medium” (1 × 104 GE/mL), and “low” (2.5 × 103 GE/mL). RNase Away, chosen as a positive control, was the most effective cleaning solution on all three viral loads. Wet Wipes were found to be more effective than BZK wipes in the medium viral load condition. The low viral load condition was easily reset with all three cleaning solutions. These findings will enable temporal SARS-CoV-2 monitoring in indoor environments where transmission risk of the virus is high and the need to avoid individual-level sampling for privacy or compliance reasons exists. IMPORTANCE Because SARS-CoV-2, the virus that causes COVID-19, persists on surfaces, testing swabs taken from surfaces is useful as a monitoring tool. This approach is especially valuable in school settings, where there are cost and privacy concerns that are eliminated by taking a single sample from a classroom. However, the virus persists for days to weeks on surface samples, so it is impossible to tell whether positive detection events on consecutive days are a persistent signal or new infectious cases and therefore whether the positive individuals have been successfully removed from the classroom. We compare several methods for cleaning “sentinel cards” to show that this approach can be used to identify new SARS-CoV-2 signals day to day. The results are important for determining how to monitor classrooms and other indoor environments for SARS-CoV-2 virus.

ABSTRACT Surface sampling for SARS-CoV-2 RNA detection has shown considerable promise to detect exposure of built environments to infected individuals shedding virus who would not otherwise be detected. Here, we compare two popular sampling media (VTM and SDS) and two popular workflows (Thermo and PerkinElmer) for implementation of a surface sampling program suitable for environmental monitoring in public schools. We find that the SDS/Thermo pipeline shows superior sensitivity and specificity, but that the VTM/PerkinElmer pipeline is still sufficient to support surface surveillance in any indoor setting with stable cohorts of occupants (e.g., schools, prisons, group homes, etc.) and may be used to leverage existing investments in infrastructure. IMPORTANCE The ongoing COVID-19 pandemic has claimed the lives of over 5 million people worldwide. Due to high density occupancy of indoor spaces for prolonged periods of time, schools are often of concern for transmission, leading to widespread school closings to combat pandemic spread when cases rise. Since pediatric clinical testing is expensive and difficult from a consent perspective, we have deployed surface sampling in SASEA (Safer at School Early Alert), which allows for detection of SARS-CoV-2 from surfaces within a classroom. In this previous work, we developed a high-throughput method which requires robotic automation and specific reagents that are often not available for public health laboratories such as the San Diego County Public Health Laboratory (SDPHL). Therefore, we benchmarked our method (Thermo pipeline) against SDPHL’s (PerkinElmer) more widely used method for the detection and prediction of SARS-CoV-2 exposure. While our method shows superior sensitivity (false-negative rate of 9% versus 27% for SDPHL), the SDPHL pipeline is sufficient to support surface surveillance in indoor settings. These findings are important since they show that existing investments in infrastructure can be leveraged to slow the spread of SARS-CoV-2 not in just the classroom but also in prisons, nursing homes, and other high-risk, indoor settings.

Base editors (BEs) are genome editing agents that install point mutations with high efficiency and specificity. Due to their reliance on uracil and inosine DNA damage intermediates (rather than double-strand DNA breaks, or DSBs), it has been hypothesized that BEs rely on more ubiquitous DNA repair pathways than DSB-reliant genome editing methods, which require processes that are only active during certain phases of the cell cycle. We report here the first systematic study of the cell cycle-dependence of base editing using cell synchronization experiments. We find that nickase-derived BEs (which introduce DNA backbone nicks opposite the uracil or inosine base) function independently of the cell cycle, while non-nicking BEs are highly dependent on S-phase (DNA synthesis phase). We found that synchronization in G1 (growth phase) during the process of cytosine base editing causes significant increases in C•G to A•T “byproduct” introduction rates, which can be leveraged to discover new strategies for precise C•G to A•T base editing. We observe that endogenous expression levels of DNA damage repair pathways are sufficient to process base editing intermediates into desired editing outcomes, and the process of base editing does not significantly perturb transcription levels. Overall, our study provides mechanistic data demonstrating the robustness of nickase-derived BEs for performing genome editing across the cell cycle.

ABSTRACT Monitoring severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) on surfaces is emerging as an important tool for identifying past exposure to individuals shedding viral RNA. Our past work demonstrated that SARS-CoV-2 reverse transcription-quantitative PCR (RT-qPCR) signals from surfaces can identify when infected individuals have touched surfaces and when they have been present in hospital rooms or schools. However, the sensitivity and specificity of surface sampling as a method for detecting the presence of a SARS-CoV-2 positive individual, as well as guidance about where to sample, has not been established. To address these questions and to test whether our past observations linking SARS-CoV-2 abundance to Rothia sp. in hospitals also hold in a residential setting, we performed a detailed spatial sampling of three isolation housing units, assessing each sample for SARS-CoV-2 abundance by RT-qPCR, linking the results to 16S rRNA gene amplicon sequences (to assess the bacterial community at each location), and to the Cq value of the contemporaneous clinical test. Our results showed that the highest SARS-CoV-2 load in this setting is on touched surfaces, such as light switches and faucets, but a detectable signal was present in many untouched surfaces (e.g., floors) that may be more relevant in settings, such as schools where mask-wearing is enforced. As in past studies, the bacterial community predicts which samples are positive for SARS-CoV-2, with Rothia sp. showing a positive association. IMPORTANCE Surface sampling for detecting SARS-CoV-2, the virus that causes coronavirus disease 2019 (COVID-19), is increasingly being used to locate infected individuals. We tested which indoor surfaces had high versus low viral loads by collecting 381 samples from three residential units where infected individuals resided, and interpreted the results in terms of whether SARS-CoV-2 was likely transmitted directly (e.g., touching a light switch) or indirectly (e.g., by droplets or aerosols settling). We found the highest loads where the subject touched the surface directly, although enough virus was detected on indirectly contacted surfaces to make such locations useful for sampling (e.g., in schools, where students did not touch the light switches and also wore masks such that they had no opportunity to touch their face and then the object). We also documented links between the bacteria present in a sample and the SARS-CoV-2 virus, consistent with earlier studies.

Abstract Background Circular RNAs (circRNAs) are stable, often highly expressed RNA transcripts with potential to modulate other regulatory RNAs. A few circRNAs have been shown to bind RNA-binding proteins (RBPs); however, little is known about the prevalence and distribution of these interactions in different biological contexts. Methods We conduct an extensive screen of circRNA-RBP interactions in the ENCODE cell lines HepG2 and K562. We profile circRNAs in deep-sequenced total RNA samples and analyze circRNA-RBP interactions using a large set of eCLIP data with binding sites of 150 RBPs. We validate interactions for select circRNAs and RBPs by performing RNA immunoprecipitation and functionally characterize our most interesting candidates by conducting knockdown studies followed by RNA-Seq. Results We generate a comprehensive catalog of circRNA-RBP interactions in HepG2 and K562 cells. We show that KHSRP binding sites are enriched in flanking introns of circRNAs and that KHSRP depletion affects circRNA biogenesis. We identify circRNAs that are highly covered by RBP binding sites and experimentally validate individual circRNA-RBP interactions. We show that circCDYL, a highly expressed circRNA with clinical and functional implications in bladder cancer, is almost completely covered with GRWD1 binding sites in HepG2 cells, and that circCDYL depletion counteracts the effect of GRWD1 depletion. Furthermore, we confirm interactions between circCDYL and RBPs in bladder cancer cells and demonstrate that circCDYL depletion affects hallmarks of cancer and perturbs the expression of key cancer genes, e.g., TP53. Finally, we show that elevated levels of circCDYL are associated with overall survival of bladder cancer patients. Conclusions Our study demonstrates transcriptome-wide and cell-type-specific circRNA-RBP interactions that could play important regulatory roles in tumorigenesis.

Summary: Gene expression profiling and proteome analysis of normal and malignant hematopoietic stem cells (HSCs) point to shared core stemness properties. However, discordance between mRNA and protein signatures highlights an important role for post-transcriptional regulation by microRNAs (miRNAs) in governing this critical nexus. Here, we identify miR-130a as a regulator of HSC self-renewal and differentiation. Enforced expression of miR-130a impairs B lymphoid differentiation and expands long-term HSCs. Integration of protein mass spectrometry and chimeric AGO2 crosslinking and immunoprecipitation (CLIP) identifies TBL1XR1 as a primary miR-130a target, whose loss of function phenocopies miR-130a overexpression. Moreover, we report that miR-130a is highly expressed in t(8;21) acute myeloid leukemia (AML), where it is critical for maintaining the oncogenic molecular program mediated by the AML1-ETO complex. Our study establishes that identification of the comprehensive miRNA targetome within primary cells enables discovery of genes and molecular networks underpinning stemness properties of normal and leukemic cells.

Abstract Background A critical step in uncovering rules of RNA processing is to study the in vivo regulatory networks of RNA binding proteins (RBPs). Crosslinking and immunoprecipitation (CLIP) methods enable mapping RBP targets transcriptome-wide, but methodological differences present challenges to large-scale analysis across datasets. The development of enhanced CLIP (eCLIP) enabled the mapping of targets for 150 RBPs in K562 and HepG2, creating a unique resource of RBP interactomes profiled with a standardized methodology in the same cell types. Results Our analysis of 223 eCLIP datasets reveals a range of binding modalities, including highly resolved positioning around splicing signals and mRNA untranslated regions that associate with distinct RBP functions. Quantification of enrichment for repetitive and abundant multicopy elements reveals 70% of RBPs have enrichment for non-mRNA element classes, enables identification of novel ribosomal RNA processing factors and sites, and suggests that association with retrotransposable elements reflects multiple RBP mechanisms of action. Analysis of spliceosomal RBPs indicates that eCLIP resolves AQR association after intronic lariat formation, enabling identification of branch points with single-nucleotide resolution, and provides genome-wide validation for a branch point-based scanning model for 3′ splice site recognition. Finally, we show that eCLIP peak co-occurrences across RBPs enable the discovery of novel co-interacting RBPs. Conclusions This work reveals novel insights into RNA biology by integrated analysis of eCLIP profiling of 150 RBPs with distinct functions. Further, our quantification of both mRNA and other element association will enable further research to identify novel roles of RBPs in regulating RNA processing.

The dearth of glioblastoma model systems that accurately recapitulate the disease remains a challenge. Here, the authors develop cancer avatars using genetically engineered human induced pluripotent cells.

ABSTRACT Environmental monitoring in public spaces can be used to identify surfaces contaminated by persons with coronavirus disease 2019 (COVID-19) and inform appropriate infection mitigation responses. Research groups have reported detection of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) on surfaces days or weeks after the virus has been deposited, making it difficult to estimate when an infected individual may have shed virus onto a SARS-CoV-2-positive surface, which in turn complicates the process of establishing effective quarantine measures. In this study, we determined that reverse transcription-quantitative PCR (RT-qPCR) detection of viral RNA from heat-inactivated particles experiences minimal decay over 7 days of monitoring on eight out of nine surfaces tested. The properties of the studied surfaces result in RT-qPCR signatures that can be segregated into two material categories, rough and smooth, where smooth surfaces have a lower limit of detection. RT-qPCR signal intensity (average quantification cycle [Cq]) can be correlated with surface viral load using only one linear regression model per material category. The same experiment was performed with untreated viral particles on one surface from each category, with essentially identical results. The stability of RT-qPCR viral signal demonstrates the need to clean monitored surfaces after sampling to establish temporal resolution. Additionally, these findings can be used to minimize the number of materials and time points tested and allow for the use of heat-inactivated viral particles when optimizing environmental monitoring methods. IMPORTANCE Environmental monitoring is an important tool for public health surveillance, particularly in settings with low rates of diagnostic testing. Time between sampling public environments, such as hospitals or schools, and notifying stakeholders of the results should be minimal, allowing decisions to be made toward containing outbreaks of coronavirus disease 2019 (COVID-19). The Safer At School Early Alert program (SASEA) (https://saseasystem.org/), a large-scale environmental monitoring effort in elementary school and child care settings, has processed >13,000 surface samples for SARS-CoV-2, detecting viral signals from 574 samples. However, consecutive detection events necessitated the present study to establish appropriate response practices around persistent viral signals on classroom surfaces. Other research groups and clinical labs developing environmental monitoring methods may need to establish their own correlation between RT-qPCR results and viral load, but this work provides evidence justifying simplified experimental designs, like reduced testing materials and the use of heat-inactivated viral particles.

Summary: Selective 2′-hydroxyl acylation analyzed by primer extension (SHAPE) structure probing techniques characterize the secondary structure of RNA molecules, which influence their functions and interactions. A variation of SHAPE, footprinting SHAPE (fSHAPE), probes RNA in the presence and absence of protein to identify RNA bases that hydrogen-bond with protein. SHAPE or fSHAPE coupled with enhanced crosslinking and immunoprecipitation (SHAPE-eCLIP or fSHAPE-eCLIP) pulls down RNAs bound by any protein of interest and returns their structure or protein interaction information, respectively. Here, we describe detailed protocols for SHAPE-eCLIP and fSHAPE-eCLIP and an analysis protocol for fSHAPE.For complete details on the use and execution of these protocols, please refer to Corley et al. (2020).

Summary: RNA binding protein (RBP) expression is finite. For RBPs that are vastly outnumbered by their potential target sites, a simple competition for binding can set the magnitude of post-transcriptional control. Here, we show that LIN28, best known for its direct regulation of let-7 miRNA biogenesis, is also indirectly regulated by its widespread binding of non-miRNA transcripts. Approximately 99% of LIN28 binding sites are found on non-miRNA transcripts, like protein coding and ribosomal RNAs. These sites are bound specifically and strongly, but they do not appear to mediate direct post-transcriptional regulation. Instead, non-miRNA sites act to sequester LIN28 protein and effectively change its functional availability, thus impeding the regulation of let-7 in cells. Together, these data show that the binding properties of the transcriptome broadly influence the ability of an RBP to mediate changes in RNA metabolism and gene expression.

Differential binding of RNA-binding proteins mediated by genetic variants (GVs) can influence posttranscriptional regulation. Here, the authors develop BEAPR, a computational approach to identify allele-specific binding events in eCLIP-Seq data.

ABSTRACT Although SARS-CoV-2 is primarily transmitted by respiratory droplets and aerosols, transmission by fomites remains plausible. During Halloween, a major event for children in numerous countries, SARS-CoV-2 transmission risk via candy fomites worries many parents. To address this concern, we enrolled 10 recently diagnosed asymptomatic or mildly/moderately symptomatic COVID-19 patients to handle typical Halloween candy (pieces individually wrapped) under three conditions: normal handling with unwashed hands, deliberate coughing and extensive touching, and normal handling following handwashing. We then used a factorial design to subject the candies to two posthandling treatments: no washing (untreated) and household dishwashing detergent. We measured SARS-CoV-2 load by reverse transcriptase quantitative PCR (RT-qPCR) and loop-mediated isothermal amplification (LAMP). From the candies not washed posthandling, we detected SARS-CoV-2 on 60% of candies that were deliberately coughed on, 60% of candies normally handled with unwashed hands, but only 10% of candies handled after hand washing. We found that treating candy with dishwashing detergent reduced SARS-CoV-2 load by 62.1% in comparison to untreated candy. Taken together, these results suggest that although the risk of transmission of SARS-CoV-2 by fomites is low even from known COVID-19 patients, viral RNA load can be reduced to near zero by the combination of handwashing by the infected patient and ≥1 min detergent treatment after collection. We also found that the inexpensive and fast LAMP protocol was more than 80% concordant with RT-qPCR. IMPORTANCE The COVID-19 pandemic is leading to important tradeoffs between risk of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) transmission and mental health due to deprivation from normal activities, with these impacts being especially profound in children. Due to the ongoing pandemic, Halloween activities will be curtailed as a result of the concern that candy from strangers might act as fomites. Here, we demonstrate that these risks can be mitigated by ensuring that, prior to handling candy, the candy giver washes their hands and, after receipt, by washing candy with household dishwashing detergent. Even in the most extreme case, with candy deliberately coughed on by known COVID-19 patients, viral load was reduced dramatically after washing with household detergent. We conclude that with reasonable precautions, even if followed only by either the candy giver or the candy recipient, the risk of viral transmission by this route is very low.

Summary: Site-directed RNA editing approaches offer great potential to correct genetic mutations in somatic cells while avoiding permanent off-target genomic edits. Nuclease-dead RNA-targeting CRISPR-Cas systems recruit functional effectors to RNA molecules in a programmable fashion. Here, we demonstrate a Streptococcus pyogenes Cas9-ADAR2 fusion system that uses a 3′ modified guide RNA (gRNA) to enable adenosine-to-inosine (A-to-I) editing of specific bases on reporter and endogenously expressed mRNAs. Due to the sufficient nature of the 3′ gRNA extension sequence, we observe that Cas9 gRNA spacer sequences are dispensable for directed RNA editing, revealing that Cas9 can act as an RNA-aptamer-binding protein. We demonstrate that Cas9-based A-to-I editing is comparable in on-target efficiency and off-target specificity with Cas13 RNA editing versions. This study provides a systematic benchmarking of RNA-targeting CRISPR-Cas designs for reversible nucleotide-level conversion at the transcriptome level.

Summary: Large-scale collections of induced pluripotent stem cells (iPSCs) could serve as powerful model systems for examining how genetic variation affects biology and disease. Here we describe the iPSCORE resource: a collection of systematically derived and characterized iPSC lines from 222 ethnically diverse individuals that allows for both familial and association-based genetic studies. iPSCORE lines are pluripotent with high genomic integrity (no or low numbers of somatic copy-number variants) as determined using high-throughput RNA-sequencing and genotyping arrays, respectively. Using iPSCs from a family of individuals, we show that iPSC-derived cardiomyocytes demonstrate gene expression patterns that cluster by genetic background, and can be used to examine variants associated with physiological and disease phenotypes. The iPSCORE collection contains representative individuals for risk and non-risk alleles for 95% of SNPs associated with human phenotypes through genome-wide association studies. Our study demonstrates the utility of iPSCORE for examining how genetic variants influence molecular and physiological traits in iPSCs and derived cell lines. : Working as part of the NHLBI NextGen consortium, Panopoulos and colleagues report the derivation and characterization of 222 publicly available iPSCs from ethnically diverse individuals with corresponding genomic data including SNP arrays, RNA-seq, and whole-genome sequencing. This collection provides a powerful resource to investigate the function of genetic variants. Keywords: iPSCORE, iPSC, GWAS, molecular traits, physiological traits, cardiac disease, NHLBI Next Gen, LQT2, KCNH2, iPSC-derived cardiomyocytes

Summary: Reprogramming somatic cells to induced pluripotent stem cells (iPSCs) offers the possibility of studying the molecular mechanisms underlying human diseases in cell types difficult to extract from living patients, such as neurons and cardiomyocytes. To date, studies have been published that use small panels of iPSC-derived cell lines to study monogenic diseases. However, to study complex diseases, where the genetic variation underlying the disorder is unknown, a sizable number of patient-specific iPSC lines and controls need to be generated. Currently the methods for deriving and characterizing iPSCs are time consuming, expensive, and, in some cases, descriptive but not quantitative. Here we set out to develop a set of simple methods that reduce cost and increase throughput in the characterization of iPSC lines. Specifically, we outline methods for high-throughput quantification of surface markers, gene expression analysis of in vitro differentiation potential, and evaluation of karyotype with markedly reduced cost. : Working as part of the NHLBI NextGen consortium, D'Antonio and colleagues developed three simple methods that reduce cost and increase throughput in the characterization of iPSCs. These methods include: (1) fluorescent cell barcoding flow cytometry to investigate heterogeneity; (2) gene expression analysis to examine in vitro differentiation potential; and (3) high-resolution digital karyotyping to detect chromosomal aberrations. Keywords: induced pluripotent stem cells, high-throughput methods, pluripotency characterization, differentiation potential, flow cytometry, qPCR, SNP arrays, digital karyotyping, fluorescent cell barcoding

Summary: Stress granule (SG) formation is frequently accompanied by ubiquitin proteasome system (UPS) impairment and ubiquitylated protein accumulation. SGs, ubiquitin, and UPS components co-localize, but the relationship between the ubiquitin pathway and SGs has not been systematically characterized. We utilize pharmacological inhibition of either the ubiquitin- or NEDD8-activating enzyme (UAE or NAE) to probe whether active ubiquitylation or neddylation modulate SG dynamics. We show that UAE inhibition results in rapid loss of global protein ubiquitylation using ubiquitin-specific proteomics. Critically, inhibiting neither UAE nor NAE significantly affected SG formation or disassembly, indicating that active protein ubiquitylation or neddylation is dispensable for SG dynamics. Using antibodies with varying preference for free ubiquitin or polyubiquitin and fluorescently tagged ubiquitin variants in combination with UAE inhibition, we show that SGs co-localize primarily with unconjugated ubiquitin rather than polyubiquitylated proteins. These findings clarify the role of ubiquitin in SG biology and suggest that free ubiquitin may alter SG protein interactions. : Protein ubiquitylation has been implicated in pathways by which cellular stress induces the formation of stress granules (SGs) and affects protein homeostasis through the ubiquitin proteasome system. Markmiller et al. show that ubiquitylation is dispensable for SG dynamics and that SGs co-localize primarily with free ubiquitin rather than polyubiquitylated proteins. Keywords: ubiquitin, stress granule, Nedd8, oxidative stress, sodium arsenite, protein homeostasis, centrosome, proteasome, protein aggregation, neurodegeneration

Abnormal functions of RNA-binding proteins TAF15, FUS and TDP43 are associated with amyotrophic lateral sclerosis. Here, Kapeli et al. characterize the RNA targets of TAF15 and identify points of convergence and divergence between the targets of TAF15, FUS and TDP43 in several neuronal systems.

Alternative splicing (AS) defects that adversely affect gene expression and function have been identified in diabetic hearts; however, the mechanisms responsible are largely unknown. Here, we show that the RNA-binding protein RBFOX2 contributes to transcriptome changes under diabetic conditions. RBFOX2 controls AS of genes with important roles in heart function relevant to diabetic cardiomyopathy. RBFOX2 protein levels are elevated in diabetic hearts despite low RBFOX2 AS activity. A dominant-negative (DN) isoform of RBFOX2 that blocks RBFOX2-mediated AS is generated in diabetic hearts. DN RBFOX2 interacts with wild-type (WT) RBFOX2, and ectopic expression of DN RBFOX2 inhibits AS of RBFOX2 targets. Notably, DN RBFOX2 expression is specific to diabetes and occurs at early stages before cardiomyopathy symptoms appear. Importantly, DN RBFOX2 expression impairs intracellular calcium release in cardiomyocytes. Our results demonstrate that RBFOX2 dysregulation by DN RBFOX2 is an early pathogenic event in diabetic hearts.

Human pluripotent stem cells (hPSCs) require precise control of post-transcriptional RNA networks to maintain proliferation and survival. Using enhanced UV crosslinking and immunoprecipitation (eCLIP), we identify RNA targets of the IMP/IGF2BP family of RNA-binding proteins in hPSCs. At the broad region and binding site levels, IMP1 and IMP2 show reproducible binding to a large and overlapping set of 3′ UTR-enriched targets. RNA Bind-N-seq applied to recombinant full-length IMP1 and IMP2 reveals CA-rich motifs that are enriched in eCLIP-defined binding sites. We observe that IMP1 loss in hPSCs recapitulates IMP1 phenotypes, including a reduction in cell adhesion and increase in cell death. For cell adhesion, we find IMP1 maintains levels of integrin mRNA specifically regulating RNA stability of ITGB5 in hPSCs. Additionally, we show that IMP1 can be linked to hPSC survival via direct target BCL2. Thus, transcriptome-wide binding profiles identify hPSC targets modulating well-characterized IMP1 roles.

The generation of pancreas, liver, and intestine from a common pool of progenitors in the foregut endoderm requires the establishment of organ boundaries. How dorsal foregut progenitors activate pancreatic genes and evade the intestinal lineage choice remains unclear. Here, we identify Pdx1 and Sox9 as cooperative inducers of a gene regulatory network that distinguishes the pancreatic from the intestinal lineage. Genetic studies demonstrate dual and cooperative functions for Pdx1 and Sox9 in pancreatic lineage induction and repression of the intestinal lineage choice. Pdx1 and Sox9 bind to regulatory sequences near pancreatic and intestinal differentiation genes and jointly regulate their expression, revealing direct cooperative roles for Pdx1 and Sox9 in gene activation and repression. Our study identifies Pdx1 and Sox9 as important regulators of a transcription factor network that initiates pancreatic fate and sheds light on the gene regulatory circuitry that governs the development of distinct organs from multi-lineage-competent foregut progenitors.

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

The human brain is composed of a complex assembly of about 171 billion heterogeneous cellular units (86 billion neurons and 85 billion non-neuronal glia cells). A comprehensive description of brain cells is necessary to understand the nervous system in health and disease. Recently, advances in genomics have permitted the accurate analysis of the full transcriptome of single cells (scRNA-seq). We have built upon such technical progress to combine scRNA-seq with patch-clamping electrophysiological recording and morphological analysis of single human neurons in vitro. This new powerful method, referred to as Patch-seq, enables a thorough, multimodal profiling of neurons and permits us to expose the links between functional properties, morphology, and gene expression. Here, we present a detailed Patch-seq protocol for isolating single neurons from in vitro neuronal cultures. We have validated the Patch-seq whole-transcriptome profiling method with human neurons generated from embryonic and induced pluripotent stem cells (ESCs/iPSCs) derived from healthy subjects, but the procedure may be applied to any kind of cell type in vitro. Patch-seq may be used on neurons in vitro to profile cell types and states in depth to unravel the human molecular basis of neuronal diversity and investigate the cellular mechanisms underlying brain disorders.

Inadequate adenosine-to-inosine editing of noncoding regions occurs in disease but is often uncorrelated with ADAR levels, underscoring the need to study deaminase-independent control of editing. C. elegans have two ADAR proteins, ADR-2 and the theoretically catalytically inactive ADR-1. Using high-throughput RNA sequencing of wild-type and adr mutant worms, we expand the repertoire of C. elegans edited transcripts over 5-fold and confirm that ADR-2 is the only active deaminase in vivo. Despite lacking deaminase function, ADR-1 affects editing of over 60 adenosines within the 3′ UTRs of 16 different mRNAs. Furthermore, ADR-1 interacts directly with ADR-2 substrates, even in the absence of ADR-2, and mutations within its double-stranded RNA (dsRNA) binding domains abolish both binding and editing regulation. We conclude that ADR-1 acts as a major regulator of editing by binding ADR-2 substrates in vivo. These results raise the possibility that other dsRNA binding proteins, including the inactive human ADARs, regulate RNA editing through deaminase-independent mechanisms.

Understanding how RNA binding proteins control the splicing code is fundamental to human biology and disease. Here, we present a comprehensive study to elucidate how heterogeneous nuclear ribonucleoparticle (hnRNP) proteins, among the most abundant RNA binding proteins, coordinate to regulate alternative pre-mRNA splicing (AS) in human cells. Using splicing-sensitive microarrays, crosslinking and immunoprecipitation coupled with high-throughput sequencing (CLIP-seq), and cDNA sequencing, we find that more than half of all AS events are regulated by multiple hnRNP proteins and that some combinations of hnRNP proteins exhibit significant synergy, whereas others act antagonistically. Our analyses reveal position-dependent RNA splicing maps, in vivo consensus binding sites, a surprising level of cross- and autoregulation among hnRNP proteins, and the coordinated regulation by hnRNP proteins of dozens of other RNA binding proteins and genes associated with cancer. Our findings define an unprecedented degree of complexity and compensatory relationships among hnRNP proteins and their splicing targets that likely confer robustness to cells.