Identification of Required Host Factors for SARS-CoV-2 Infection in Human Cells.

To better understand host-virus genetic dependencies and find potential therapeutic targets for COVID-19, we performed a genome-scale CRISPR loss-of-function screen to identify host factors required for SARS-CoV-2 viral infection of human alveolar epithelial cells. Top-ranked genes cluster into distinct pathways, including the vacuolar ATPase proton pump, Retromer, and Commander complexes. We validate these gene targets using several orthogonal methods such as CRISPR knockout, RNA interference knockdown, and small-molecule inhibitors. Using single-cell RNA-sequencing, we identify shared transcriptional changes in cholesterol biosynthesis upon loss of top-ranked genes. In addition, given the key role of the ACE2 receptor in the early stages of viral entry, we show that loss of RAB7A reduces viral entry by sequestering the ACE2 receptor inside cells. Overall, this work provides a genome-scale, quantitative resource of the impact of the loss of each host gene on fitness/response to viral infection.

Integrated analysis of multimodal single-cell data.

The simultaneous measurement of multiple modalities represents an exciting frontier for single-cell genomics and necessitates computational methods that can define cellular states based on multimodal data. Here, we introduce "weighted-nearest neighbor" analysis, an unsupervised framework to learn the relative utility of each data type in each cell, enabling an integrative analysis of multiple modalities. We apply our procedure to a CITE-seq dataset of 211,000 human peripheral blood mononuclear cells (PBMCs) with panels extending to 228 antibodies to construct a multimodal reference atlas of the circulating immune system. Multimodal analysis substantially improves our ability to resolve cell states, allowing us to identify and validate previously unreported lymphoid subpopulations. Moreover, we demonstrate how to leverage this reference to rapidly map new datasets and to interpret immune responses to vaccination and coronavirus disease 2019 (COVID-19). Our approach represents a broadly applicable strategy to analyze single-cell multimodal datasets and to look beyond the transcriptome toward a unified and multimodal definition of cellular identity.

Comprehensive Integration of Single-Cell Data.

Single-cell transcriptomics has transformed our ability to characterize cell states, but deep biological understanding requires more than a taxonomic listing of clusters. As new methods arise to measure distinct cellular modalities, a key analytical challenge is to integrate these datasets to better understand cellular identity and function. Here, we develop a strategy to "anchor" diverse datasets together, enabling us to integrate single-cell measurements not only across scRNA-seq technologies, but also across different modalities. After demonstrating improvement over existing methods for integrating scRNA-seq data, we anchor scRNA-seq experiments with scATAC-seq to explore chromatin differences in closely related interneuron subsets and project protein expression measurements onto a bone marrow atlas to characterize lymphocyte populations. Lastly, we harmonize in situ gene expression and scRNA-seq datasets, allowing transcriptome-wide imputation of spatial gene expression patterns. Our work presents a strategy for the assembly of harmonized references and transfer of information across datasets.

Clonal lineage tracing reveals shared origin of conventional and plasmacytoid dendritic cells.

Developmental origins of dendritic cells (DCs) including conventional DCs (cDCs, comprising cDC1 and cDC2 subsets) and plasmacytoid DCs (pDCs) remain unclear. We studied DC development in unmanipulated adult mice using inducible lineage tracing combined with clonal DNA "barcoding" and single-cell transcriptome and phenotype analysis (CITE-seq). Inducible tracing of Cx3cr1+ hematopoietic progenitors in the bone marrow showed that they simultaneously produce all DC subsets including pDCs, cDC1s, and cDC2s. Clonal tracing of hematopoietic stem cells (HSCs) and of Cx3cr1+ progenitors revealed clone sharing between cDC1s and pDCs, but not between the two cDC subsets or between pDCs and B cells. Accordingly, CITE-seq analyses of differentiating HSCs and Cx3cr1+ progenitors identified progressive stages of pDC development including Cx3cr1+ Ly-6D+ pro-pDCs that were distinct from lymphoid progenitors. These results reveal the shared origin of pDCs and cDCs and suggest a revised scheme of DC development whereby pDCs share clonal relationship with cDC1s.

Mapping genotypes to chromatin accessibility profiles in single cells.

In somatic tissue differentiation, chromatin accessibility changes govern priming and precursor commitment towards cellular fates1-3. Therefore, somatic mutations are likely to alter chromatin accessibility patterns, as they disrupt differentiation topologies leading to abnormal clonal outgrowth. However, defining the impact of somatic mutations on the epigenome in human samples is challenging due to admixed mutated and wild-type cells. Here, to chart how somatic mutations disrupt epigenetic landscapes in human clonal outgrowths, we developed genotyping of targeted loci with single-cell chromatin accessibility (GoT-ChA). This high-throughput platform links genotypes to chromatin accessibility at single-cell resolution across thousands of cells within a single assay. We applied GoT-ChA to CD34+ cells from patients with myeloproliferative neoplasms with JAK2V617F-mutated haematopoiesis. Differential accessibility analysis between wild-type and JAK2V617F-mutant progenitors revealed both cell-intrinsic and cell-state-specific shifts within mutant haematopoietic precursors, including cell-intrinsic pro-inflammatory signatures in haematopoietic stem cells, and a distinct profibrotic inflammatory chromatin landscape in megakaryocytic progenitors. Integration of mitochondrial genome profiling and cell-surface protein expression measurement allowed expansion of genotyping onto DOGMA-seq through imputation, enabling single-cell capture of genotypes, chromatin accessibility, RNA expression and cell-surface protein expression. Collectively, we show that the JAK2V617F mutation leads to epigenetic rewiring in a cell-intrinsic and cell type-specific manner, influencing inflammation states and differentiation trajectories. We envision that GoT-ChA will empower broad future investigations of the critical link between somatic mutations and epigenetic alterations across clonal populations in malignant and non-malignant contexts.

A genome-scale screen for synthetic drivers of T cell proliferation.

The engineering of autologous patient T cells for adoptive cell therapies has revolutionized the treatment of several types of cancer1. However, further improvements are needed to increase response and cure rates. CRISPR-based loss-of-function screens have been limited to negative regulators of T cell functions2-4 and raise safety concerns owing to the permanent modification of the genome. Here we identify positive regulators of T cell functions through overexpression of around 12,000 barcoded human open reading frames (ORFs). The top-ranked genes increased the proliferation and activation of primary human CD4+ and CD8+ T cells and their secretion of key cytokines such as interleukin-2 and interferon-γ. In addition, we developed the single-cell genomics method OverCITE-seq for high-throughput quantification of the transcriptome and surface antigens in ORF-engineered T cells. The top-ranked ORF-lymphotoxin-ß receptor (LTBR)-is typically expressed in myeloid cells but absent in lymphocytes. When overexpressed in T cells, LTBR induced profound transcriptional and epigenomic remodelling, leading to increased T cell effector functions and resistance to exhaustion in chronic stimulation settings through constitutive activation of the canonical NF-κB pathway. LTBR and other highly ranked genes improved the antigen-specific responses of chimeric antigen receptor T cells and γδ T cells, highlighting their potential for future cancer-agnostic therapies5. Our results provide several strategies for improving next-generation T cell therapies by the induction of synthetic cell programmes.

Efficient combinatorial targeting of RNA transcripts in single cells with Cas13 RNA Perturb-seq.

Pooled CRISPR screens coupled with single-cell RNA-sequencing have enabled systematic interrogation of gene function and regulatory networks. Here, we introduce Cas13 RNA Perturb-seq (CaRPool-seq), which leverages the RNA-targeting CRISPR-Cas13d system and enables efficient combinatorial perturbations alongside multimodal single-cell profiling. CaRPool-seq encodes multiple perturbations on a cleavable CRISPR array that is associated with a detectable barcode sequence, allowing for the simultaneous targeting of multiple genes. We compared CaRPool-seq to existing Cas9-based methods, highlighting its unique strength to efficiently profile combinatorially perturbed cells. Finally, we apply CaRPool-seq to perform multiplexed combinatorial perturbations of myeloid differentiation regulators in an acute myeloid leukemia (AML) model system and identify extensive interactions between different chromatin regulators that can enhance or suppress AML differentiation phenotypes.

Somatic mutations and cell identity linked by Genotyping of Transcriptomes.

Defining the transcriptomic identity of malignant cells is challenging in the absence of surface markers that distinguish cancer clones from one another, or from admixed non-neoplastic cells. To address this challenge, here we developed Genotyping of Transcriptomes (GoT), a method to integrate genotyping with high-throughput droplet-based single-cell RNA sequencing. We apply GoT to profile 38,290 CD34+ cells from patients with CALR-mutated myeloproliferative neoplasms to study how somatic mutations corrupt the complex process of human haematopoiesis. High-resolution mapping of malignant versus normal haematopoietic progenitors revealed an increasing fitness advantage with myeloid differentiation of cells with mutated CALR. We identified the unfolded protein response as a predominant outcome of CALR mutations, with a considerable dependency on cell identity, as well as upregulation of the NF-κB pathway specifically in uncommitted stem cells. We further extended the GoT toolkit to genotype multiple targets and loci that are distant from transcript ends. Together, these findings reveal that the transcriptional output of somatic mutations in myeloproliferative neoplasms is dependent on the native cell identity.

Multimodal single-cell analysis of cutaneous T-cell lymphoma reveals distinct subclonal tissue-dependent signatures.

Cutaneous T-cell lymphoma (CTCL) is a heterogeneous group of mature T-cell neoplasms characterized by the accumulation of clonal malignant CD4+ T cells in the skin. The most common variant of CTCL, mycosis fungoides (MF ), is confined to the skin in early stages but can be accompanied by extracutaneous dissemination of malignant T cells to the blood and lymph nodes in advanced stages of disease. Sézary syndrome (SS), a leukemic form of disease, is characterized by significant blood involvement. Little is known about the transcriptional and genomic relationship between skin- and blood-residing malignant T cells in CTCL. To identify and interrogate malignant clones in matched skin and blood from patients with leukemic MF and SS, we combine T-cell receptor clonotyping with quantification of gene expression and cell surface markers at the single cell level. Our data reveal clonal evolution at a transcriptional and genetic level within the malignant populations of individual patients. We highlight highly consistent transcriptional signatures delineating skin- and blood-derived malignant T cells. Analysis of these 2 populations suggests that environmental cues, along with genetic aberrations, contribute to transcriptional profiles of malignant T cells. Our findings indicate that the skin microenvironment in CTCL promotes a transcriptional response supporting rapid malignant expansion, as opposed to the quiescent state observed in the blood, potentially influencing efficacy of therapies. These results provide insight into tissue-specific characteristics of cancerous cells and underscore the need to address the patients' individual malignant profiles at the time of therapy to eliminate all subclones.

Multiplexed detection of proteins, transcriptomes, clonotypes and CRISPR perturbations in single cells.

Multimodal single-cell assays provide high-resolution snapshots of complex cell populations, but are mostly limited to transcriptome plus an additional modality. Here, we describe expanded CRISPR-compatible cellular indexing of transcriptomes and epitopes by sequencing (ECCITE-seq) for the high-throughput characterization of at least five modalities of information from each single cell. We demonstrate application of ECCITE-seq to multimodal CRISPR screens with robust direct single-guide RNA capture and to clonotype-aware multimodal phenotyping of cancer samples.

YES1 amplification is a mechanism of acquired resistance to EGFR inhibitors identified by transposon mutagenesis and clinical genomics.

In ∼30% of patients with EGFR-mutant lung adenocarcinomas whose disease progresses on EGFR inhibitors, the basis for acquired resistance remains unclear. We have integrated transposon mutagenesis screening in an EGFR-mutant cell line and clinical genomic sequencing in cases of acquired resistance to identify mechanisms of resistance to EGFR inhibitors. The most prominent candidate genes identified by insertions in or near the genes during the screen were MET, a gene whose amplification is known to mediate resistance to EGFR inhibitors, and the gene encoding the Src family kinase YES1. Cell clones with transposon insertions that activated expression of YES1 exhibited resistance to all three generations of EGFR inhibitors and sensitivity to pharmacologic and siRNA-mediated inhibition of YES1 Analysis of clinical genomic sequencing data from cases of acquired resistance to EGFR inhibitors revealed amplification of YES1 in five cases, four of which lacked any other known mechanisms of resistance. Preinhibitor samples, available for two of the five patients, lacked YES1 amplification. None of 136 postinhibitor samples had detectable amplification of other Src family kinases (SRC and FYN). YES1 amplification was also found in 2 of 17 samples from ALK fusion-positive lung cancer patients who had progressed on ALK TKIs. Taken together, our findings identify acquired amplification of YES1 as a recurrent and targetable mechanism of resistance to EGFR inhibition in EGFR-mutant lung cancers and demonstrate the utility of transposon mutagenesis in discovering clinically relevant mechanisms of drug resistance.

Simultaneous epitope and transcriptome measurement in single cells.

High-throughput single-cell RNA sequencing has transformed our understanding of complex cell populations, but it does not provide phenotypic information such as cell-surface protein levels. Here, we describe cellular indexing of transcriptomes and epitopes by sequencing (CITE-seq), a method in which oligonucleotide-labeled antibodies are used to integrate cellular protein and transcriptome measurements into an efficient, single-cell readout. CITE-seq is compatible with existing single-cell sequencing approaches and scales readily with throughput increases.

Single-cell stabilization method identifies gonadotrope transcriptional dynamics and pituitary cell type heterogeneity.

Immediate-early response genes (IEGs) are rapidly and transiently induced following an extracellular signal. Elucidating the IEG response patterns in single cells (SCs) requires assaying large numbers of timed samples at high accuracy while minimizing handling effects. To achieve this, we developed and validated RNA stabilization Buffer for Examination of Single-cell Transcriptomes (RNA-Best), a versatile single-step cell and tissue preservation protocol that stabilizes RNA in intact SCs without perturbing transcription patterns. We characterize for the first time SC heterogeneity in IEG responses to pulsatile gonadotropin-releasing hormone (GnRH) stimuli in pituitary gonadotrope cells. Our study identifies a gene-specific hierarchical pattern of all-or-none transcript induction elicited by increasing concentrations of GnRH. This quantal pattern of gene activation raises the possibility that IEG activation, when accurately resolved at the SC level, may be mediated by gene bits that behave as pure binary switches.

Neural specificity of the RNA-binding protein Elav is achieved by post-transcriptional repression in non-neural tissues.

Drosophila Elav is the founding member of the conserved family of Hu RNA-binding proteins (RBPs), which play crucial and diverse roles in post-transcriptional regulation. Elav has long served as the canonical neuronal marker. Surprisingly, although Elav has a well-characterized neural cis-regulatory module, we find endogenous Elav is also ubiquitously transcribed and post-transcriptionally repressed in non-neural settings. Mutant clones of multiple miRNA pathway components derepress ubiquitous Elav protein. Our re-annotation of the elav transcription unit shows not only that it generates extended 3' UTR isoforms, but also that its universal 3' UTR isoform is much longer than previously believed. This longer common 3' UTR includes multiple conserved, high-affinity sites for the miR-279/996 family. Of several miRNA mutants tested, endogenous Elav and a transgenic elav 3' UTR sensor are derepressed in mutant clones of mir-279/996 We also observe cross-repression of Elav by Mei-P26, another RBP derepressed in non-neural miRNA pathway clones. Ubiquitous Elav has regulatory capacity, since derepressed Elav can stabilize an Elav-responsive sensor. Repression of Elav in non-neural territories is crucial as misexpression here has profoundly adverse consequences. Altogether, we define unexpected post-transcriptional mechanisms that direct appropriate cell type-specific expression of a conserved neural RBP.

An extensive allelic series of Drosophila kae1 mutants reveals diverse and tissue-specific requirements for t6A biogenesis.

N(6)-threonylcarbamoyl-adenosine (t6A) is one of the few RNA modifications that is universally present in life. This modification occurs at high frequency at position 37 of most tRNAs that decode ANN codons, and stabilizes cognate anticodon-codon interactions. Nearly all genetic studies of the t6A pathway have focused on single-celled organisms. In this study, we report the isolation of an extensive allelic series in the Drosophila ortholog of the core t6A biosynthesis factor Kae1. kae1 hemizygous larvae exhibit decreases in t6A that correlate with allele strength; however, we still detect substantial t6A-modified tRNAs even during the extended larval phase of null alleles. Nevertheless, complementation of Drosophila Kae1 and other t6A factors in corresponding yeast null mutants demonstrates that these metazoan genes execute t6A synthesis. Turning to the biological consequences of t6A loss, we characterize prominent kae1 melanotic masses and show that they are associated with lymph gland overgrowth and ectopic generation of lamellocytes. On the other hand, kae1 mutants exhibit other phenotypes that reflect insufficient tissue growth. Interestingly, whole-tissue and clonal analyses show that strongly mitotic tissues such as imaginal discs are exquisitely sensitive to loss of kae1, whereas nonproliferating tissues are less affected. Indeed, despite overt requirements of t6A for growth of many tissues, certain strong kae1 alleles achieve and sustain enlarged body size during their extended larval phase. Our studies highlight tissue-specific requirements of the t6A pathway in a metazoan context and provide insights into the diverse biological roles of this fundamental RNA modification during animal development and disease.

The miR-310/13 cluster antagonizes ß-catenin function in the regulation of germ and somatic cell differentiation in the Drosophila testis.

MicroRNAs (miRNAs) are regulators of global gene expression and function in a broad range of biological processes. Recent studies have suggested that miRNAs can function as tumor suppressors or oncogenes by modulating the activities of evolutionarily conserved signaling pathways that are commonly dysregulated in cancer. We report the identification of the miR-310 to miR-313 (miR-310/13) cluster as a novel antagonist of Wingless (Drosophila Wnt) pathway activity in a functional screen for Drosophila miRNAs. We demonstrate that miR-310/13 can modulate Armadillo (Arm; Drosophila ß-catenin) expression and activity by directly targeting the 3'-UTRs of arm and pangolin (Drosophila TCF) in vivo. Notably, the miR-310/13-deficient flies exhibit abnormal germ and somatic cell differentiation in the male gonad, which can be rescued by reducing Arm protein levels or activity. Our results implicate a previously unrecognized function for miR-310/13 in dampening the activity of Arm in early somatic and germline progenitor cells, whereby inappropriate/sustained activation of Arm-mediated signaling or cell adhesion may impact normal differentiation in the Drosophila male gonad.

Intertwined pathways for Argonaute-mediated microRNA biogenesis in Drosophila.

Although Dicer is essential for general microRNA (miRNA) biogenesis, vertebrate mir-451 is Dicer independent. Instead, its short pre-miRNA hairpin is 'sliced' by Ago2, then 3'-resected into mature miRNAs. Here, we show that Drosophila cells and animals generate functional small RNAs from mir-451-type precursors. However, their bulk maturation arrests as Ago-cleaved pre-miRNAs, which mostly associate with the RNAi effector AGO2. Routing of pre-mir-451 hairpins to the miRNA effector AGO1 was inhibited by Dicer-1 and its partner Loqs. Loss of these miRNA factors promoted association of pre-mir-451 with AGO1, which sliced them and permitted maturation into ∼ 23-26 nt products. The difference was due to the 3' modification of single-stranded species in AGO2 by Hen1 methyltransferase, whose depletion permitted 3' trimming of Ago-cleaved pre-miRNAs in AGO2. Surprisingly, Nibbler, a 3'-5' exoribonuclease that trims 'long' mature miRNAs in AGO1, antagonized miR-451 processing. We used an in vitro reconstitution assay to identify a soluble, EDTA-sensitive activity that resects sliced pre-miRNAs in AGO1 complexes. Finally, we use deep sequencing to show that depletion of dicer-1 increases the diversity of small RNAs in AGO1, including some candidate mir-451-like loci. Altogether, we document unexpected aspects of miRNA biogenesis and Ago sorting, and provide insights into maturation of Argonaute-cleaved miRNA substrates.

miR-33a/b contribute to the regulation of fatty acid metabolism and insulin signaling.

Cellular imbalances of cholesterol and fatty acid metabolism result in pathological processes, including atherosclerosis and metabolic syndrome. Recent work from our group and others has shown that the intronic microRNAs hsa-miR-33a and hsa-miR-33b are located within the sterol regulatory element-binding protein-2 and -1 genes, respectively, and regulate cholesterol homeostasis in concert with their host genes. Here, we show that miR-33a and -b also regulate genes involved in fatty acid metabolism and insulin signaling. miR-33a and -b target key enzymes involved in the regulation of fatty acid oxidation, including carnitine O-octaniltransferase, carnitine palmitoyltransferase 1A, hydroxyacyl-CoA-dehydrogenase, Sirtuin 6 (SIRT6), and AMP kinase subunit-α. Moreover, miR-33a and -b also target the insulin receptor substrate 2, an essential component of the insulin-signaling pathway in the liver. Overexpression of miR-33a and -b reduces both fatty acid oxidation and insulin signaling in hepatic cell lines, whereas inhibition of endogenous miR-33a and -b increases these two metabolic pathways. Together, these data establish that miR-33a and -b regulate pathways controlling three of the risk factors of metabolic syndrome, namely levels of HDL, triglycerides, and insulin signaling, and suggest that inhibitors of miR-33a and -b may be useful in the treatment of this growing health concern.

Drosophila Argonaute 1 and its miRNA biogenesis partners are required for oocyte formation and germline cell division.

Argonaute 1 (Ago1) is a member of the Argonaute/PIWI protein family involved in small RNA-mediated gene regulation. In Drosophila, Ago1 plays a specific role in microRNA (miRNA) biogenesis and function. Previous studies have demonstrated that Ago1 regulates the fate of germline stem cells. However, the function of Ago1 in other aspects of oogenesis is still elusive. Here we report the function of Ago1 in developing egg chambers. We find that Ago1 protein is enriched in the oocytes and is also highly expressed in the cytoplasm of follicle cells. Clonal analysis of multiple ago1 mutant alleles shows that many mutant egg chambers contain only 8 nurse cells without an oocyte which is phenocopied in dicer-1, pasha and drosha mutants. Our results suggest that Ago1 and its miRNA biogenesis partners play a role in oocyte determination and germline cell division in Drosophila.