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.

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.

Thousands of putative enhancers are characterized in the human genome, yet few have been shown to have a functional role in cancer progression. Inhibiting oncokinases, such as EGFR, ALK, ERBB2, and BRAF, is a mainstay of current cancer therapy but is hindered by innate drug resistance mediated by up-regulation of the HGF receptor, MET. The mechanisms mediating such genomic responses to targeted therapy are unknown. Here, we identify lineage-specific enhancers at the MET locus for multiple common tumor types, including a melanoma lineage-specific enhancer 63 kb downstream from the MET TSS. This enhancer displays inducible chromatin looping with the MET promoter to up-regulate MET expression upon BRAF inhibition. Epigenomic analysis demonstrated that the melanocyte-specific transcription factor, MITF, mediates this enhancer function. Targeted genomic deletion (<7 bp) of the MITF motif within the MET enhancer suppressed inducible chromatin looping and innate drug resistance, while maintaining MITF-dependent, inhibitor-induced melanoma cell differentiation. Epigenomic analysis can thus guide functional disruption of regulatory DNA to decouple pro- and anti-oncogenic functions of a dominant transcription factor and block innate resistance to oncokinase therapy.

To elucidate mechanisms of cancer progression, we generated inducible human neoplasia in three-dimensionally intact epithelial tissue. Gene expression profiling of both epithelia and stroma at specific time points during tumor progression revealed sequential enrichment of genes mediating discrete biologic functions in each tissue compartment. A core cancer progression signature was distilled using the increased signaling specificity of downstream oncogene effectors and subjected to network modeling. Network topology predicted that tumor development depends on specific extracellular matrix-interacting network hubs. Blockade of one such hub, the beta1 integrin subunit, disrupted network gene expression and attenuated tumorigenesis in vivo. Thus, integrating network modeling and temporal gene expression analysis of inducible human neoplasia provides an approach to prioritize and characterize genes functioning in cancer progression.

Clear cell renal cell carcinomas (ccRCC) show a broad range of clinical behavior, and prognostic biomarkers are needed to stratify patients for appropriate management. We sought to determine whether long intergenic non-coding RNAs (lincRNAs) might predict patient survival. Candidate prognostic lincRNAs were identified by mining The Cancer Genome Atlas (TCGA) transcriptome (RNA-seq) data on 466 ccRCC cases (randomized into discovery and validation sets) annotated for ~21,000 lncRNAs. A previously uncharacterized lincRNA, SLINKY (Survival-predictive LINcRNA in KidneY cancer), was the top-ranked prognostic lincRNA, and validated in an independent University of Tokyo cohort (P=0.004). In multivariable analysis, SLINKY expression predicted overall survival independent of tumor stage and grade [TCGA HR=3.5 (CI, 2.2-5.7), P < 0.001; Tokyo HR=8.4 (CI, 1.8-40.2), P = 0.007], and by decision tree, ROC and decision curve analysis, added independent prognostic value. In ccRCC cell lines, SLINKY knockdown reduced cancer cell proliferation (with cell-cycle G1 arrest) and induced transcriptome changes enriched for cell proliferation and survival processes. Notably, the genes affected by SLINKY knockdown in cell lines were themselves prognostic and correlated with SLINKY expression in the ccRCC patient samples. From a screen for binding partners, we identified direct binding of SLINKY to Heterogeneous Nuclear Ribonucleoprotein K (HNRNPK), whose knockdown recapitulated SLINKY knockdown phenotypes. Thus, SLINKY is a robust prognostic biomarker in ccRCC, where it functions possibly together with HNRNPK in cancer cell proliferation.

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.

Gene expression is controlled by transcription factors (TFs) that bind cognate DNA motif sequences in cis-regulatory elements (CREs). The combinations of DNA motifs acting within homeostasis and disease, however, are unclear. Gene expression, chromatin accessibility, TF footprinting, and H3K27ac-dependent DNA looping data were generated and a random-forest-based model was applied to identify 7,531 cell-type-specific cis-regulatory modules (CRMs) across 15 diploid human cell types. A co-enrichment framework within CRMs nominated 838 cell-type-specific, recurrent heterotypic DNA motif combinations (DMCs), which were functionally validated using massively parallel reporter assays. Cancer cells engaged DMCs linked to neoplasia-enabling processes operative in normal cells while also activating new DMCs only seen in the neoplastic state. This integrative approach identifies cell-type-specific cis-regulatory combinatorial DNA motifs in diverse normal and diseased human cells and represents a general framework for deciphering cis-regulatory sequence logic in gene regulation.

An approach combining fluorescence-activated cell sorting and high-throughput DNA sequencing (FACS-seq) was employed to determine the efficiency of start codon recognition for all possible translation initiation sites (TIS) utilizing AUG start codons. Using FACS-seq, we measured translation from a genetic reporter library representing all 65,536 possible TIS sequences spanning the -6 to +5 positions. We found that the motif RYMRMVAUGGC enhanced start codon recognition and translation efficiency. However, dinucleotide interactions, which cannot be conveyed by a single motif, were also important for modeling TIS efficiency. Our dataset combined with modeling allowed us to predict genome-wide translation initiation efficiency for all mRNA transcripts. Additionally, we screened somatic TIS mutations associated with tumorigenesis to identify candidate driver mutations consistent with known tumor expression patterns. Finally, we implemented a quantitative leaky scanning model to predict alternative initiation sites that produce truncated protein isoforms and compared predictions with ribosome footprint profiling data. The comprehensive analysis of the TIS sequence space enables quantitative predictions of translation initiation based on genome sequence.

T cells create vast amounts of diversity in the genes that encode their T cell receptors (TCRs), which enables individual clones to recognize specific peptide-major histocompatibility complex (MHC) ligands. Here we combined sequencing of the TCR-encoding genes with assay for transposase-accessible chromatin with sequencing (ATAC-seq) analysis at the single-cell level to provide information on the TCR specificity and epigenomic state of individual T cells. By using this approach, termed transcript-indexed ATAC-seq (T-ATAC-seq), we identified epigenomic signatures in immortalized leukemic T cells, primary human T cells from healthy volunteers and primary leukemic T cells from patient samples. In peripheral blood CD4+ T cells from healthy individuals, we identified cis and trans regulators of naive and memory T cell states and found substantial heterogeneity in surface-marker-defined T cell populations. In patients with a leukemic form of cutaneous T cell lymphoma, T-ATAC-seq enabled identification of leukemic and nonleukemic regulatory pathways in T cells from the same individual by allowing separation of the signals that arose from the malignant clone from the background T cell noise. Thus, T-ATAC-seq is a new tool that enables analysis of epigenomic landscapes in clonal T cells and should be valuable for studies of T cell malignancy, immunity and immunotherapy.

Aberrations of protein-coding genes are a focus of cancer genomics; however, the impact of oncogenes on expression of the ~50% of transcripts without protein-coding potential, including long noncoding RNAs (lncRNAs), has been largely uncharacterized. Activating mutations in the BRAF oncogene are present in >70% of melanomas, 90% of which produce active mutant BRAF(V600E) protein. To define the impacts of oncogenic BRAF on the melanocyte transcriptome, massively parallel cDNA sequencing (RNA-seq) was performed on genetically matched normal human melanocytes with and without BRAF(V600E) expression. To enhance potential disease relevance by verifying expression of altered genes in BRAF-driven cancer tissue, parallel RNA-seq was also undertaken of two BRAF(V600E)-mutant human melanomas. BRAF(V600E) regulated expression of 1027 protein-coding transcripts and 39 annotated lncRNAs, as well as 70 unannotated, potentially novel, intergenic transcripts. These transcripts display both tissue-specific and multi-tissue expression profiles and harbor distinctive regulatory chromatin marks and transcription factor binding sites indicative of active transcription. Coding potential analysis of the 70 unannotated transcripts suggested that most may represent newly identified lncRNAs. BRAF-regulated lncRNA 1 (BANCR) was identified as a recurrently overexpressed, previously unannotated 693-bp transcript on chromosome 9 with a potential functional role in melanoma cell migration. BANCR knockdown reduced melanoma cell migration, and this could be rescued by the chemokine CXCL11. Combining RNA-seq of oncogene-expressing normal cells with RNA-seq of their corresponding human cancers may represent a useful approach to discover new oncogene-regulated RNA transcripts of potential clinical relevance in cancer.

Erk1/2 mitogen-activated protein kinases (MAPKs) are often hyperactivated in human cancers, where they affect multiple processes, including proliferation. However, the effects of Erk1/2 loss in normal epithelial tissue, the setting of most extracellular signal-regulated kinase (Erk)-associated neoplasms, are unknown. In epidermis, loss of Erk1 or Erk2 individually has no effect, whereas simultaneous Erk1/2 depletion inhibits cell division, demonstrating that these MAPKs are necessary for normal tissue self-renewal. Growth inhibition caused by Erk1/2 loss is rescued by reintroducing Erk2, but not by activating Erk effectors that promote G1 cell cycle progression. Unlike fibroblasts, in which Erk1/2 loss decreases cyclin D1 expression and induces G1/S arrest, Erk1/2 loss in epithelial cells reduces cyclin B1 and c-Fos expression and induces G2/M arrest while disrupting a gene regulatory network centered on cyclin B1-Cdc2. Thus, the cell cycle stages at which Erk1/2 activity is required vary by cell type, with Erk1/2 functioning in epithelial cells to enable progression through G2/M.

Progenitor cells maintain self-renewing tissues throughout life by sustaining their capacity for proliferation while suppressing cell cycle exit and terminal differentiation. DNA methylation provides a potential epigenetic mechanism for the cellular memory needed to preserve the somatic progenitor state through repeated cell divisions. DNA methyltransferase 1 (DNMT1) maintains DNA methylation patterns after cellular replication. Although dispensable for embryonic stem cell maintenance, the role for DNMT1 in maintaining the progenitor state in constantly replenished somatic tissues, such as mammalian epidermis, is unclear. Here we show that DNMT1 is essential for epidermal progenitor cell function. DNMT1 protein was found enriched in undifferentiated cells, where it was required to retain proliferative stamina and suppress differentiation. In tissue, DNMT1 depletion led to exit from the progenitor cell compartment, premature differentiation and eventual tissue loss. Genome-wide analysis showed that a significant portion of epidermal differentiation gene promoters were methylated in self-renewing conditions but were subsequently demethylated during differentiation. Furthermore, UHRF1 (refs 9, 10), a component of the DNA methylation machinery that targets DNMT1 to hemi-methylated DNA, is also necessary to suppress premature differentiation and sustain proliferation. In contrast, Gadd45A and B, which promote active DNA demethylation, are required for full epidermal differentiation gene induction. These data demonstrate that proteins involved in the dynamic regulation of DNA methylation patterns are required for progenitor maintenance and self-renewal in mammalian somatic 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.

Open chromatin regions are correlated with active regulatory elements in development and are dysregulated in diseases. The BAF (SWI/SNF) complex is essential for development, and has been demonstrated to remodel reconstituted chromatin in vitro and to control the accessibility of a few individual regions in vivo. However, it remains unclear where and how BAF controls the open chromatin landscape to regulate developmental processes, such as human epidermal differentiation.

Genome conformation is central to gene control but challenging to interrogate. Here we present HiChIP, a protein-centric chromatin conformation method. HiChIP improves the yield of conformation-informative reads by over 10-fold and lowers the input requirement over 100-fold relative to that of ChIA-PET. HiChIP of cohesin reveals multiscale genome architecture with greater signal-to-background ratios than those of in situ Hi-C.

SARS-CoV2 being highly infectious has been particularly effective in causing widespread infection globally and more variants of SARS-CoV2 are constantly being reported with increased genomic surveillance. In particular, the focus is on mutations of Spike protein, which binds human ACE2 protein enabling SARS-CoV2 entry and infection. Here we present a rapid experimental method leveraging the speed and flexibility of Mircoscale Thermopheresis (MST) to characterize the interaction between Spike Receptor Binding Domain (RBD) and human ACE2 protein. The B.1.351 variant harboring three mutations, (E484K, N501Y, and K417N) binds the ACE2 at nearly five-fold greater affinity than the original SARS-COV-2 RBD. We also find that the B.1.1.7 variant, binds two-fold more tightly to ACE2 than the SARS-COV-2 RBD.

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.

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.

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.

Somatic progenitors sustain tissue self-renewal while suppressing premature differentiation. Protein arginine methyltransferases (PRMTs) affect many processes; however, their role in progenitor function is incompletely understood. PRMT1 was found to be the most highly expressed PRMT in epidermal progenitors and the most downregulated PRMT during differentiation. In targeted mouse knockouts and in long-term regenerated human mosaic epidermis in vivo, epidermal PRMT1 loss abolished progenitor self-renewal and led to premature differentiation. Mass spectrometry of the PRMT1 protein interactome identified the CSNK1a1 kinase, which also proved essential for progenitor maintenance. CSNK1a1 directly bound and phosphorylated PRMT1 to control its genomic targeting to PRMT1-sustained proliferation genes as well as PRMT1-suppressed differentiation genes. Among the latter were GRHL3, whose derepression was required for the premature differentiation seen with PRMT1 and CSNK1a1 loss. Maintenance of the progenitors thus requires cooperation by PRMT1 and CSNK1a1 to sustain proliferation gene expression and suppress premature differentiation driven by GRHL3.