Combinatorial single-cell CRISPR screens by direct guide RNA capture and targeted sequencing.

Single-cell CRISPR screens enable the exploration of mammalian gene function and genetic regulatory networks. However, use of this technology has been limited by reliance on indirect indexing of single-guide RNAs (sgRNAs). Here we present direct-capture Perturb-seq, a versatile screening approach in which expressed sgRNAs are sequenced alongside single-cell transcriptomes. Direct-capture Perturb-seq enables detection of multiple distinct sgRNA sequences from individual cells and thus allows pooled single-cell CRISPR screens to be easily paired with combinatorial perturbation libraries that contain dual-guide expression vectors. We demonstrate the utility of this approach for high-throughput investigations of genetic interactions and, leveraging this ability, dissect epistatic interactions between cholesterol biogenesis and DNA repair. Using direct capture Perturb-seq, we also show that targeting individual genes with multiple sgRNAs per cell improves efficacy of CRISPR interference and activation, facilitating the use of compact, highly active CRISPR libraries for single-cell screens. Last, we show that hybridization-based target enrichment permits sensitive, specific sequencing of informative transcripts from single-cell RNA-seq experiments.

Rare Pulmonary Neuroendocrine Cells Are Stem Cells Regulated by Rb, p53, and Notch.

Pulmonary neuroendocrine (NE) cells are neurosensory cells sparsely distributed throughout the bronchial epithelium, many in innervated clusters of 20-30 cells. Following lung injury, NE cells proliferate and generate other cell types to promote epithelial repair. Here, we show that only rare NE cells, typically 2-4 per cluster, function as stem cells. These fully differentiated cells display features of classical stem cells. Most proliferate (self-renew) following injury, and some migrate into the injured area. A week later, individual cells, often just one per cluster, lose NE identity (deprogram), transit amplify, and reprogram to other fates, creating large clonal repair patches. Small cell lung cancer (SCLC) tumor suppressors regulate the stem cells: Rb and p53 suppress self-renewal, whereas Notch marks the stem cells and initiates deprogramming and transit amplification. We propose that NE stem cells give rise to SCLC, and transformation results from constitutive activation of stem cell renewal and inhibition of deprogramming.

Automated cell-type classification in intact tissues by single-cell molecular profiling.

A major challenge in biology is identifying distinct cell classes and mapping their interactions in vivo. Tissue-dissociative technologies enable deep single cell molecular profiling but do not provide spatial information. We developed a proximity ligation in situ hybridization technology (PLISH) with exceptional signal strength, specificity, and sensitivity in tissue. Multiplexed data sets can be acquired using barcoded probes and rapid label-image-erase cycles, with automated calculation of single cell profiles, enabling clustering and anatomical re-mapping of cells. We apply PLISH to expression profile ~2900 cells in intact mouse lung, which identifies and localizes known cell types, including rare ones. Unsupervised classification of the cells indicates differential expression of 'housekeeping' genes between cell types, and re-mapping of two sub-classes of Club cells highlights their segregated spatial domains in terminal airways. By enabling single cell profiling of various RNA species in situ, PLISH can impact many areas of basic and medical research.

Ultrasensitive optical imaging with lanthanide lumiphores.

In principle, the millisecond emission lifetimes of lanthanide chelates should enable their ultrasensitive detection in biological systems by time-resolved optical microscopy. In practice, however, lanthanide imaging techniques have provided no better sensitivity than conventional fluorescence microscopy. Here, we identified three fundamental problems that have impeded lanthanide microscopy: low photon flux, inefficient excitation, and optics-derived background luminescence. We overcame these limitations with a new lanthanide imaging modality, transreflected illumination with luminescence resonance energy transfer (trLRET), which increases the time-integrated signal intensities of lanthanide lumiphores by 170-fold and the signal-to-background ratios by 75-fold. We demonstrate that trLRET provides at least an order-of-magnitude increase in detection sensitivity over that of conventional epifluorescence microscopy when used to visualize endogenous protein expression in zebrafish embryos. We also show that trLRET can be used to optically detect molecular interactions in vivo. trLRET promises to unlock the full potential of lanthanide lumiphores for ultrasensitive, autofluorescence-free biological imaging.

G Protein-Coupled Receptor Endocytosis Confers Uniformity in Responses to Chemically Distinct Ligands.

The ability of chemically distinct ligands to produce different effects on the same G protein-coupled receptor (GPCR) has interesting therapeutic implications, but, if excessively propagated downstream, would introduce biologic noise compromising cognate ligand detection. We asked whether cells have the ability to limit the degree to which chemical diversity imposed at the ligand-GPCR interface is propagated to the downstream signal. We carried out an unbiased analysis of the integrated cellular response elicited by two chemically and pharmacodynamically diverse ß-adrenoceptor agonists, isoproterenol and salmeterol. We show that both ligands generate an identical integrated response, and that this stereotyped output requires endocytosis. We further demonstrate that the endosomal ß2-adrenergic receptor signal confers uniformity on the downstream response because it is highly sensitive and saturable. Based on these findings, we propose that GPCR signaling from endosomes functions as a biologic noise filter to enhance reliability of cognate ligand detection.

Pericytes are progenitors for coronary artery smooth muscle.

Epicardial cells on the heart's surface give rise to coronary artery smooth muscle cells (caSMCs) located deep in the myocardium. However, the differentiation steps between epicardial cells and caSMCs are unknown as are the final maturation signals at coronary arteries. Here, we use clonal analysis and lineage tracing to show that caSMCs derive from pericytes, mural cells associated with microvessels, and that these cells are present in adults. During development following the onset of blood flow, pericytes at arterial remodeling sites upregulate Notch3 while endothelial cells express Jagged-1. Deletion of Notch3 disrupts caSMC differentiation. Our data support a model wherein epicardial-derived pericytes populate the entire coronary microvasculature, but differentiate into caSMCs at arterial remodeling zones in response to Notch signaling. Our data are the first demonstration that pericytes are progenitors for smooth muscle, and their presence in adult hearts reveals a new potential cell type for targeting during cardiovascular disease.

Automated Analysis and Classification of Histological Tissue Features by Multi-Dimensional Microscopic Molecular Profiling.

Characterization of the molecular attributes and spatial arrangements of cells and features within complex human tissues provides a critical basis for understanding processes involved in development and disease. Moreover, the ability to automate steps in the analysis and interpretation of histological images that currently require manual inspection by pathologists could revolutionize medical diagnostics. Toward this end, we developed a new imaging approach called multidimensional microscopic molecular profiling (MMMP) that can measure several independent molecular properties in situ at subcellular resolution for the same tissue specimen. MMMP involves repeated cycles of antibody or histochemical staining, imaging, and signal removal, which ultimately can generate information analogous to a multidimensional flow cytometry analysis on intact tissue sections. We performed a MMMP analysis on a tissue microarray containing a diverse set of 102 human tissues using a panel of 15 informative antibody and 5 histochemical stains plus DAPI. Large-scale unsupervised analysis of MMMP data, and visualization of the resulting classifications, identified molecular profiles that were associated with functional tissue features. We then directly annotated H&E images from this MMMP series such that canonical histological features of interest (e.g. blood vessels, epithelium, red blood cells) were individually labeled. By integrating image annotation data, we identified molecular signatures that were associated with specific histological annotations and we developed statistical models for automatically classifying these features. The classification accuracy for automated histology labeling was objectively evaluated using a cross-validation strategy, and significant accuracy (with a median per-pixel rate of 77% per feature from 15 annotated samples) for de novo feature prediction was obtained. These results suggest that high-dimensional profiling may advance the development of computer-based systems for automatically parsing relevant histological and cellular features from molecular imaging data of arbitrary human tissue samples, and can provide a framework and resource to spur the optimization of these technologies.

Transcriptome-wide mapping of pseudouridines: pseudouridine synthases modify specific mRNAs in S. cerevisiae.

We developed a novel technique, called pseudouridine site identification sequencing (PSI-seq), for the transcriptome-wide mapping of pseudouridylation sites with single-base resolution from cellular RNAs based on the induced termination of reverse transcription specifically at pseudouridines following CMCT treatment. PSI-seq analysis of RNA samples from S. cerevisiae correctly detected all of the 43 known pseudouridines in yeast 18S and 25S ribosomal RNA with high specificity. Moreover, application of PSI-seq to the yeast transcriptome revealed the presence of site-specific pseudouridylation within dozens of mRNAs, including RPL11a, TEF1, and other genes implicated in translation. To identify the mechanisms responsible for mRNA pseudouridylation, we genetically deleted candidate pseudouridine synthase (Pus) enzymes and reconstituted their activities in vitro. These experiments demonstrated that the Pus1 enzyme was necessary and sufficient for pseudouridylation of RPL11a mRNA, whereas Pus4 modified TEF1 mRNA, and Pus6 pseudouridylated KAR2 mRNA. Finally, we determined that modification of RPL11a at Ψ -68 was observed in RNA from the related yeast S. mikitae, and Ψ -239 in TEF1 mRNA was maintained in S. mikitae as well as S. pombe, indicating that these pseudouridylations are ancient, evolutionarily conserved RNA modifications. This work establishes that site-specific pseudouridylation of eukaryotic mRNAs is a genetically programmed RNA modification that naturally occurs in multiple yeast transcripts via distinct mechanisms, suggesting that mRNA pseudouridylation may provide an important novel regulatory function. The approach and strategies that we report here should be generally applicable to the discovery of pseudouridylation, or other RNA modifications, in diverse biological contexts.

The yeast Rab GTPase Ypt1 modulates unfolded protein response dynamics by regulating the stability of HAC1 RNA.

The unfolded protein response (UPR) is a conserved mechanism that mitigates accumulation of unfolded proteins in the ER. The yeast UPR is subject to intricate post-transcriptional regulation, involving recruitment of the RNA encoding the Hac1 transcription factor to the ER and its unconventional splicing. To investigate the mechanisms underlying regulation of the UPR, we screened the yeast proteome for proteins that specifically interact with HAC1 RNA. Protein microarray experiments revealed that HAC1 interacts specifically with small ras GTPases of the Ypt family. We characterized the interaction of HAC1 RNA with one of these proteins, the yeast Rab1 homolog Ypt1. We found that Ypt1 protein specifically associated in vivo with unspliced HAC1 RNA. This association was disrupted by conditions that impaired protein folding in the ER and induced the UPR. Also, the Ypt1-HAC1 interaction depended on IRE1 and ADA5, the two genes critical for UPR activation. Decreasing expression of the Ypt1 protein resulted in a reduced rate of HAC1 RNA decay, leading to significantly increased levels of both unspliced and spliced HAC1 RNA, and delayed attenuation of the UPR, when ER stress was relieved. Our findings establish that Ypt1 contributes to regulation of UPR signaling dynamics by promoting the decay of HAC1 RNA, suggesting a potential regulatory mechanism for linking vesicle trafficking to the UPR and ER homeostasis.

Identification of RNA recognition elements in the Saccharomyces cerevisiae transcriptome.

Post-transcriptional regulation of gene expression, including mRNA localization, translation and decay, is ubiquitous yet still largely unexplored. How is the post-transcriptional regulatory program of each mRNA encoded in its sequence? Hundreds of specific RNA-binding proteins (RBPs) appear to play roles in mediating the post-transcriptional regulatory program, akin to the roles of specific DNA-binding proteins in transcription. As a step toward decoding the regulatory programs encoded in each mRNA, we focused on specific mRNA-protein interactions. We computationally analyzed the sequences of Saccharomyces cerevisiae mRNAs bound in vivo by 29 specific RBPs, identifying eight novel candidate motifs and confirming or extending six earlier reported recognition elements. Biochemical selections for RNA sequences selectively recognized by 12 yeast RBPs yielded novel motifs bound by Pin4, Nsr1, Hrb1, Gbp2, Sgn1 and Mrn1, and recovered the known recognition elements for Puf3, She2, Vts1 and Whi3. Most of the RNA elements we uncovered were associated with coherent mRNA expression changes and were significantly conserved in related yeasts, supporting their functional importance and suggesting that the corresponding RNA-protein interactions are evolutionarily conserved.

Diverse RNA-binding proteins interact with functionally related sets of RNAs, suggesting an extensive regulatory system.

RNA-binding proteins (RBPs) have roles in the regulation of many post-transcriptional steps in gene expression, but relatively few RBPs have been systematically studied. We searched for the RNA targets of 40 proteins in the yeast Saccharomyces cerevisiae: a selective sample of the approximately 600 annotated and predicted RBPs, as well as several proteins not annotated as RBPs. At least 33 of these 40 proteins, including three of the four proteins that were not previously known or predicted to be RBPs, were reproducibly associated with specific sets of a few to several hundred RNAs. Remarkably, many of the RBPs we studied bound mRNAs whose protein products share identifiable functional or cytotopic features. We identified specific sequences or predicted structures significantly enriched in target mRNAs of 16 RBPs. These potential RNA-recognition elements were diverse in sequence, structure, and location: some were found predominantly in 3'-untranslated regions, others in 5'-untranslated regions, some in coding sequences, and many in two or more of these features. Although this study only examined a small fraction of the universe of yeast RBPs, 70% of the mRNA transcriptome had significant associations with at least one of these RBPs, and on average, each distinct yeast mRNA interacted with three of the RBPs, suggesting the potential for a rich, multidimensional network of regulation. These results strongly suggest that combinatorial binding of RBPs to specific recognition elements in mRNAs is a pervasive mechanism for multi-dimensional regulation of their post-transcriptional fate.

Flexibility and constraint in the nucleosome core landscape of Caenorhabditis elegans chromatin.

Nucleosome positions within the chromatin landscape are known to serve as a major determinant of DNA accessibility to transcription factors and other interacting components. To delineate nucleosomal patterns in a model genetic organism, Caenorhabditis elegans, we have carried out a genome-wide analysis in which DNA fragments corresponding to nucleosome cores were liberated using an enzyme (micrococcal nuclease) with a strong preference for cleavage in non-nucleosomal regions. Sequence analysis of 284,091 putative nucleosome cores obtained in this manner from a mixed-stage population of C. elegans reveals a combined picture of flexibility and constraint in nucleosome positioning. As has previously been observed in studies of individual loci in diverse biological systems, we observe areas in the genome where nucleosomes can adopt a wide variety of positions in a given region, areas with little or no nucleosome coverage, and areas where nucleosomes reproducibly adopt a specific positional pattern. In addition to illuminating numerous aspects of chromatin structure for C. elegans, this analysis provides a reference from which to begin an investigation of relationships between the nucleosomal pattern, chromosomal architecture, and lineage-based gene activity on a genome-wide scale.