Glucose dissociates DDX21 dimers to regulate mRNA splicing and tissue differentiation.

Glucose is a universal bioenergy source; however, its role in controlling protein interactions is unappreciated, as are its actions during differentiation-associated intracellular glucose elevation. Azido-glucose click chemistry identified glucose binding to a variety of RNA binding proteins (RBPs), including the DDX21 RNA helicase, which was found to be essential for epidermal differentiation. Glucose bound the ATP-binding domain of DDX21, altering protein conformation, inhibiting helicase activity, and dissociating DDX21 dimers. Glucose elevation during differentiation was associated with DDX21 re-localization from the nucleolus to the nucleoplasm where DDX21 assembled into larger protein complexes containing RNA splicing factors. DDX21 localized to specific SCUGSDGC motif in mRNA introns in a glucose-dependent manner and promoted the splicing of key pro-differentiation genes, including GRHL3, KLF4, OVOL1, and RBPJ. These findings uncover a biochemical mechanism of action for glucose in modulating the dimerization and function of an RNA helicase essential for tissue differentiation.

Multimodal Analysis of Composition and Spatial Architecture in Human Squamous Cell Carcinoma.

To define the cellular composition and architecture of cutaneous squamous cell carcinoma (cSCC), we combined single-cell RNA sequencing with spatial transcriptomics and multiplexed ion beam imaging from a series of human cSCCs and matched normal skin. cSCC exhibited four tumor subpopulations, three recapitulating normal epidermal states, and a tumor-specific keratinocyte (TSK) population unique to cancer, which localized to a fibrovascular niche. Integration of single-cell and spatial data mapped ligand-receptor networks to specific cell types, revealing TSK cells as a hub for intercellular communication. Multiple features of potential immunosuppression were observed, including T regulatory cell (Treg) co-localization with CD8 T cells in compartmentalized tumor stroma. Finally, single-cell characterization of human tumor xenografts and in vivo CRISPR screens identified essential roles for specific tumor subpopulation-enriched gene networks in tumorigenesis. These data define cSCC tumor and stromal cell subpopulations, the spatial niches where they interact, and the communicating gene networks that they engage in cancer.

Coupled Single-Cell CRISPR Screening and Epigenomic Profiling Reveals Causal Gene Regulatory Networks.

Here, we present Perturb-ATAC, a method that combines multiplexed CRISPR interference or knockout with genome-wide chromatin accessibility profiling in single cells based on the simultaneous detection of CRISPR guide RNAs and open chromatin sites by assay of transposase-accessible chromatin with sequencing (ATAC-seq). We applied Perturb-ATAC to transcription factors (TFs), chromatin-modifying factors, and noncoding RNAs (ncRNAs) in ∼4,300 single cells, encompassing more than 63 genotype-phenotype relationships. Perturb-ATAC in human B lymphocytes uncovered regulators of chromatin accessibility, TF occupancy, and nucleosome positioning and identified a hierarchy of TFs that govern B cell state, variation, and disease-associated cis-regulatory elements. Perturb-ATAC in primary human epidermal cells revealed three sequential modules of cis-elements that specify keratinocyte fate. Combinatorial deletion of all pairs of these TFs uncovered their epistatic relationships and highlighted genomic co-localization as a basis for synergistic interactions. Thus, Perturb-ATAC is a powerful strategy to dissect gene regulatory networks in development and disease.

The functions and unique features of long intergenic non-coding RNA.

Long intergenic non-coding RNA (lincRNA) genes have diverse features that distinguish them from mRNA-encoding genes and exercise functions such as remodelling chromatin and genome architecture, RNA stabilization and transcription regulation, including enhancer-associated activity. Some genes currently annotated as encoding lincRNAs include small open reading frames (smORFs) and encode functional peptides and thus may be more properly classified as coding RNAs. lincRNAs may broadly serve to fine-tune the expression of neighbouring genes with remarkable tissue specificity through a diversity of mechanisms, highlighting our rapidly evolving understanding of the non-coding genome.

Spatially resolved clonal copy number alterations in benign and malignant tissue.

Defining the transition from benign to malignant tissue is fundamental to improving early diagnosis of cancer1. Here we use a systematic approach to study spatial genome integrity in situ and describe previously unidentified clonal relationships. We used spatially resolved transcriptomics2 to infer spatial copy number variations in >120,000 regions across multiple organs, in benign and malignant tissues. We demonstrate that genome-wide copy number variation reveals distinct clonal patterns within tumours and in nearby benign tissue using an organ-wide approach focused on the prostate. Our results suggest a model for how genomic instability arises in histologically benign tissue that may represent early events in cancer evolution. We highlight the power of capturing the molecular and spatial continuums in a tissue context and challenge the rationale for treatment paradigms, including focal therapy.

Integrating single-cell and spatial transcriptomics to elucidate intercellular tissue dynamics.

Single-cell RNA sequencing (scRNA-seq) identifies cell subpopulations within tissue but does not capture their spatial distribution nor reveal local networks of intercellular communication acting in situ. A suite of recently developed techniques that localize RNA within tissue, including multiplexed in situ hybridization and in situ sequencing (here defined as high-plex RNA imaging) and spatial barcoding, can help address this issue. However, no method currently provides as complete a scope of the transcriptome as does scRNA-seq, underscoring the need for approaches to integrate single-cell and spatial data. Here, we review efforts to integrate scRNA-seq with spatial transcriptomics, including emerging integrative computational methods, and propose ways to effectively combine current methodologies.

The Functional Proximal Proteome of Oncogenic Ras Includes mTORC2.

Proximity-dependent biotin labeling (BioID) may identify new targets for cancers driven by difficult-to-drug oncogenes such as Ras. Therefore, BioID was used with wild-type (WT) and oncogenic mutant (MT) H-, K-, and N-Ras, identifying known interactors, including Raf and PI3K, as well as a common set of 130 novel proteins proximal to all Ras isoforms. A CRISPR screen of these proteins for Ras dependence identified mTOR, which was also found proximal to MT Ras in human tumors. Oncogenic Ras directly bound two mTOR complex 2 (mTORC2) components, mTOR and MAPKAP1, to promote mTORC2 kinase activity at the plasma membrane. mTORC2 enabled the Ras pro-proliferative cell cycle transcriptional program, and perturbing the Ras-mTORC2 interaction impaired Ras-dependent neoplasia in vivo. Combining proximity-dependent proteomics with CRISPR screening identified a new set of functional Ras-associated proteins, defined mTORC2 as a new direct Ras effector, and offers a strategy for finding new proteins that cooperate with dominant oncogenes.

PROBER identifies proteins associated with programmable sequence-specific DNA in living cells.

DNA-protein interactions mediate physiologic gene regulation and may be altered by DNA variants linked to polygenic disease. To enhance the speed and signal-to-noise ratio (SNR) in the identification and quantification of proteins associated with specific DNA sequences in living cells, we developed proximal biotinylation by episomal recruitment (PROBER). PROBER uses high-copy episomes to amplify SNR, and proximity proteomics (BioID) to identify the transcription factors and additional gene regulators associated with short DNA sequences of interest. PROBER quantified both constitutive and inducible association of transcription factors and corresponding chromatin regulators to target DNA sequences and binding quantitative trait loci due to single-nucleotide variants. PROBER identified alterations in regulator associations due to cancer hotspot mutations in the hTERT promoter, indicating that these mutations increase promoter association with specific gene activators. PROBER provides an approach to rapidly identify proteins associated with specific DNA sequences and their variants in living cells.

The proximal proteome of 17 SARS-CoV-2 proteins links to disrupted antiviral signaling and host translation.

Viral proteins localize within subcellular compartments to subvert host machinery and promote pathogenesis. To study SARS-CoV-2 biology, we generated an atlas of 2422 human proteins vicinal to 17 SARS-CoV-2 viral proteins using proximity proteomics. This identified viral proteins at specific intracellular locations, such as association of accessary proteins with intracellular membranes, and projected SARS-CoV-2 impacts on innate immune signaling, ER-Golgi transport, and protein translation. It identified viral protein adjacency to specific host proteins whose regulatory variants are linked to COVID-19 severity, including the TRIM4 interferon signaling regulator which was found proximal to the SARS-CoV-2 M protein. Viral NSP1 protein adjacency to the EIF3 complex was associated with inhibited host protein translation whereas ORF6 localization with MAVS was associated with inhibited RIG-I 2CARD-mediated IFNB1 promoter activation. Quantitative proteomics identified candidate host targets for the NSP5 protease, with specific functional cleavage sequences in host proteins CWC22 and FANCD2. This data resource identifies host factors proximal to viral proteins in living human cells and nominates pathogenic mechanisms employed by SARS-CoV-2.

CALML5 is a ZNF750- and TINCR-induced protein that binds stratifin to regulate epidermal differentiation.

Outward migration of epidermal progenitors occurs with induction of hundreds of differentiation genes, but the identities of all regulators required for this process are unknown. We used laser capture microdissection followed by RNA sequencing to identify calmodulin-like 5 (CALML5) as the most enriched gene in differentiating outer epidermis. CALML5 mRNA was up-regulated by the ZNF750 transcription factor and then stabilized by the long noncoding RNA TINCR. CALML5 knockout impaired differentiation, abolished keratohyalin granules, and disrupted epidermal barrier function. Mass spectrometry identified SFN (stratifin/14-3-3σ) as a CALML5-binding protein. CALML5 interacts with SFN in suprabasal epidermis, cocontrols 13% of late differentiation genes, and modulates interaction of SFN to some of its binding partners. A ZNF750-TINCR-CALML5-SFN network is thus essential for epidermal differentiation.

Targeted Proteomic Approaches for Proteome-Wide Characterizations of the AMP-Binding Capacities of Kinases.

Kinases play important roles in cell signaling, and adenosine monophosphate (AMP) is known to modulate cellular energy homeostasis through AMP-activated protein kinase (AMPK). Here, we explored novel AMP-binding kinases by employing a desthiobiotin-conjugated AMP acyl-phosphate probe to enrich efficiently AMP-binding proteins. Together with a parallel-reaction monitoring-based targeted proteomic approach, we uncovered 195 candidate AMP-binding kinases. We also enriched desthiobiotin-labeled peptides from adenine nucleotide-binding sites of kinases and analyzed them using LC-MS/MS in the multiple-reaction monitoring mode, which resulted in the identification of 44 peptides derived from 43 kinases displaying comparable or better binding affinities toward AMP relative to adenosine triphosphate (ATP). Moreover, our proteomic data revealed a potential involvement of AMP in the MAPK pathway through binding directly to the relevant kinases, especially MEK2 and MEK3. Together, we revealed the AMP-binding capacities of a large number of kinases, and our work built a strong foundation for understanding how AMP functions as a second messenger to modulate cell signaling.

Author Correction: Methods to study RNA-protein interactions.

In the version of this paper originally published, three references were accidentally omitted: Schwartz, J. C. et al. Cell Rep. 5, 918-925 (2013); Tundup, S. et al. FEBS Lett. 580, 1285-1293 (2006); and Itri, F. et al. Biochem. Biophys. Res. Commun. 492, 67-73 (2017). The PDF and HTML versions of the paper now include these as references 58, 59, and 60, respectively, and subsequent references have been renumbered accordingly.

Methods to study RNA-protein interactions.

Noncoding RNA sequences, including long noncoding RNAs, small nucleolar RNAs, and untranslated mRNA regions, accomplish many of their diverse functions through direct interactions with RNA-binding proteins (RBPs). Recent efforts have identified hundreds of new RBPs that lack known RNA-binding domains, thus underscoring the complexity and diversity of RNA-protein complexes. Recent progress has expanded the number of methods for studying RNA-protein interactions in two general categories: approaches that characterize proteins bound to an RNA of interest (RNA-centric), and those that examine RNAs bound to a protein of interest (protein-centric). Each method has unique strengths and limitations, which makes it important to select optimal approaches for the biological question being addressed. Here we review methods for the study of RNA-protein interactions, with a focus on their suitability for specific applications.

HiChIRP reveals RNA-associated chromosome conformation.

Modular domains of long non-coding RNAs can serve as scaffolds to bring distant regions of the linear genome into spatial proximity. Here, we present HiChIRP, a method leveraging bio-orthogonal chemistry and optimized chromosome conformation capture conditions, which enables interrogation of chromatin architecture focused around a specific RNA of interest down to approximately ten copies per cell. HiChIRP of three nuclear RNAs reveals insights into promoter interactions (7SK), telomere biology (telomerase RNA component) and inflammatory gene regulation (lincRNA-EPS).

ZNF750 interacts with KLF4 and RCOR1, KDM1A, and CTBP1/2 chromatin regulators to repress epidermal progenitor genes and induce differentiation genes.

ZNF750 controls epithelial homeostasis by inhibiting progenitor genes while inducing differentiation genes, a role underscored by pathogenic ZNF750 mutations in cancer and psoriasis. How ZNF750 accomplishes these dual gene regulatory impacts is unknown. Here, we characterized ZNF750 as a transcription factor that binds both the progenitor and differentiation genes that it controls at a CCNNAGGC DNA motif. ZNF750 interacts with the pluripotency transcription factor KLF4 and chromatin regulators RCOR1, KDM1A, and CTBP1/2 through conserved PLNLS sequences. ChIP-seq (chromatin immunoprecipitation [ChIP] followed by high-throughput sequencing) and gene depletion revealed that KLF4 colocalizes ∼ 10 base pairs from ZNF750 at differentiation target genes to facilitate their activation but is unnecessary for ZNF750-mediated progenitor gene repression. In contrast, KDM1A colocalizes with ZNF750 at progenitor genes and facilitates their repression but is unnecessary for ZNF750-driven differentiation. ZNF750 thus controls differentiation in concert with RCOR1 and CTBP1/2 by acting with either KDM1A to repress progenitor genes or KLF4 to induce differentiation genes.

RNA-protein interaction detection in living cells.

RNA-protein interactions play numerous roles in cellular function and disease. Here we describe RNA-protein interaction detection (RaPID), which uses proximity-dependent protein labeling, based on the BirA* biotin ligase, to rapidly identify the proteins that bind RNA sequences of interest in living cells. RaPID displays utility in multiple applications, including in evaluating protein binding to mutant RNA motifs in human genetic disorders, in uncovering potential post-transcriptional networks in breast cancer, and in discovering essential host proteins that interact with Zika virus RNA. To improve the BirA*-labeling component of RaPID, moreover, a new mutant BirA* was engineered from Bacillus subtilis, termed BASU, that enables >1,000-fold faster kinetics and >30-fold increased signal-to-noise ratio over the prior standard Escherichia coli BirA*, thereby enabling direct study of RNA-protein interactions in living cells on a timescale as short as 1 min.

Impact of a patient-derived hepatitis C viral RNA genome with a mutated microRNA binding site.

Hepatitis C virus (HCV) depends on liver-specific microRNA miR-122 for efficient viral RNA amplification in liver cells. This microRNA interacts with two different conserved sites at the very 5' end of the viral RNA, enhancing miR-122 stability and promoting replication of the viral RNA. Treatment of HCV patients with oligonucleotides that sequester miR-122 resulted in profound loss of viral RNA in phase II clinical trials. However, some patients accumulated in their sera a viral RNA genome that contained a single cytidine to uridine mutation at the third nucleotide from the 5' genomic end. It is shown here that this C3U variant indeed displayed higher rates of replication than that of wild-type HCV when miR-122 abundance is low in liver cells. However, when miR-122 abundance is high, binding of miR-122 to site 1, most proximal to the 5' end in the C3U variant RNA, is impaired without disrupting the binding of miR-122 to site 2. As a result, C3U RNA displays a much lower rate of replication than wild-type mRNA when miR-122 abundance is high in the liver. This phenotype was accompanied by binding of a different set of cellular proteins to the 5' end of the C3U RNA genome. In particular, binding of RNA helicase DDX6 was important for displaying the C3U RNA replication phenotype in liver cells. These findings suggest that sequestration of miR-122 leads to a resistance-associated mutation that has only been observed in treated patients so far, and raises the question about the function of the C3U variant in the peripheral blood.

Profiling of rotavirus 3'UTR-binding proteins reveals the ATP synthase subunit ATP5B as a host factor that supports late-stage virus replication.

Genome replication and virion assembly of segmented RNA viruses are highly coordinated events, tightly regulated by sequence and structural elements in the UTRs of viral RNA. This process is poorly defined and likely requires the participation of host proteins in concert with viral proteins. In this study, we employed a proteomics-based approach, named RNA-protein interaction detection (RaPID), to comprehensively screen for host proteins that bind to a conserved motif within the rotavirus (RV) 3' terminus. Using this assay, we identified ATP5B, a core subunit of the mitochondrial ATP synthase, as having high affinity to the RV 3'UTR consensus sequences. During RV infection, ATP5B bound to the RV 3'UTR and co-localized with viral RNA and viroplasm. Functionally, siRNA-mediated genetic depletion of ATP5B or other ATP synthase subunits such as ATP5A1 and ATP5O reduced the production of infectious viral progeny without significant alteration of intracellular viral RNA levels or RNA translation. Chemical inhibition of ATP synthase diminished RV yield in both conventional cell culture and in human intestinal enteroids, indicating that ATP5B positively regulates late-stage RV maturation in primary intestinal epithelial cells. Collectively, our results shed light on the role of host proteins in RV genome assembly and particle formation and identify ATP5B as a novel pro-RV RNA-binding protein, contributing to our understanding of how host ATP synthases may galvanize virus growth and pathogenesis.

An improved ATAC-seq protocol reduces background and enables interrogation of frozen tissues.

We present Omni-ATAC, an improved ATAC-seq protocol for chromatin accessibility profiling that works across multiple applications with substantial improvement of signal-to-background ratio and information content. The Omni-ATAC protocol generates chromatin accessibility profiles from archival frozen tissue samples and 50-µm sections, revealing the activities of disease-associated DNA elements in distinct human brain structures. The Omni-ATAC protocol enables the interrogation of personal regulomes in tissue context and translational studies.