Nir Hacohen

Full-length transcriptome assembly from RNA-Seq data without a reference genome

Massively parallel sequencing of cDNA has enabled deep and efficient probing of transcriptomes. Current approaches for transcript reconstruction from such data often rely on aligning reads to a reference genome, and are thus unsuitable for samples with a partial or missing reference genome. Here we present the Trinity method for de novo assembly of full-length transcripts and evaluate it on samples from fission yeast, mouse and whitefly, whose reference genome is not yet available. By efficiently constructing and analyzing sets of de Bruijn graphs, Trinity fully reconstructs a large fraction of transcripts, including alternatively spliced isoforms and transcripts from recently duplicated genes. Compared with other de novo transcriptome assemblers, Trinity recovers more full-length transcripts across a broad range of expression levels, with a sensitivity similar to methods that rely on genome alignments. Our approach provides a unified solution for transcriptome reconstruction in any sample, especially in the absence of a reference genome.Massively parallel sequencing of cDNA has enabled deep and efficient probing of transcriptomes. Current approaches for transcript reconstruction from such data often rely on aligning reads to a reference genome, and are thus unsuitable for samples with a partial or missing reference genome. Here we present the Trinity method for de novo assembly of full-length transcripts and evaluate it on samples from fission yeast, mouse and whitefly, whose reference genome is not yet available. By efficiently constructing and analyzing sets of de Bruijn graphs, Trinity fully reconstructs a large fraction of transcripts, including alternatively spliced isoforms and transcripts from recently duplicated genes. Compared with other de novo transcriptome assemblers, Trinity recovers more full-length transcripts across a broad range of expression levels, with a sensitivity similar to methods that rely on genome alignments. Our approach provides a unified solution for transcriptome reconstruction in any sample, especially in the absence of a reference genome.

Chromatin signature reveals over a thousand highly conserved large non-coding RNAs in mammals

There is growing recognition that mammalian cells produce many thousands of large intergenic transcripts. However, the functional significance of these transcripts has been particularly controversial. Although there are some well-characterized examples, most (>95%) show little evidence of evolutionary conservation and have been suggested to represent transcriptional noise. Here we report a new approach to identifying large non-coding RNAs using chromatin-state maps to discover discrete transcriptional units intervening known protein-coding loci. Our approach identified ∼1,600 large multi-exonic RNAs across four mouse cell types. In sharp contrast to previous collections, these large intervening non-coding RNAs (lincRNAs) show strong purifying selection in their genomic loci, exonic sequences and promoter regions, with greater than 95% showing clear evolutionary conservation. We also developed a functional genomics approach that assigns putative functions to each lincRNA, demonstrating a diverse range of roles for lincRNAs in processes from embryonic stem cell pluripotency to cell proliferation. We obtained independent functional validation for the predictions for over 100 lincRNAs, using cell-based assays. In particular, we demonstrate that specific lincRNAs are transcriptionally regulated by key transcription factors in these processes such as p53, NFκB, Sox2, Oct4 (also known as Pou5f1) and Nanog. Together, these results define a unique collection of functional lincRNAs that are highly conserved and implicated in diverse biological processes.

A Lentiviral RNAi Library for Human and Mouse Genes Applied to an Arrayed Viral High-Content Screen

To enable arrayed or pooled loss-of-function screens in a wide range of mammalian cell types, including primary and nondividing cells, we are developing lentiviral short hairpin RNA (shRNA) libraries targeting the human and murine genomes. The libraries currently contain 104,000 vectors, targeting each of 22,000 human and mouse genes with multiple sequence-verified constructs. To test the utility of the library for arrayed screens, we developed a screen based on high-content imaging to identify genes required for mitotic progression in human cancer cells and applied it to an arrayed set of 5,000 unique shRNA-expressing lentiviruses that target 1,028 human genes. The screen identified several known and approximately 100 candidate regulators of mitotic progression and proliferation; the availability of multiple shRNAs targeting the same gene facilitated functional validation of putative hits. This work provides a widely applicable resource for loss-of-function screens, as well as a roadmap for its application to biological discovery.

Rab27a and Rab27b control different steps of the exosome secretion pathway

Exosomes are secreted membrane vesicles that share structural and biochemical characteristics with intraluminal vesicles of multivesicular endosomes (MVEs). Exosomes could be involved in intercellular communication and in the pathogenesis of infectious and degenerative diseases. The molecular mechanisms of exosome biogenesis and secretion are, however, poorly understood. Using an RNA interference (RNAi) screen, we identified five Rab GTPases that promote exosome secretion in HeLa cells. Among these, Rab27a and Rab27b were found to function in MVE docking at the plasma membrane. The size of MVEs was strongly increased by Rab27a silencing, whereas MVEs were redistributed towards the perinuclear region upon Rab27b silencing. Thus, the two Rab27 isoforms have different roles in the exosomal pathway. In addition, silencing two known Rab27 effectors, Slp4 (also known as SYTL4, synaptotagmin-like 4) and Slac2b (also known as EXPH5, exophilin 5), inhibited exosome secretion and phenocopied silencing of Rab27a and Rab27b, respectively. Our results therefore strengthen the link between MVEs and exosomes, and introduce ways of manipulating exosome secretion in vivo.

Evolution and Impact of Subclonal Mutations in Chronic Lymphocytic Leukemia

Clonal evolution is a key feature of cancer progression and relapse. We studied intratumoral heterogeneity in 149 chronic lymphocytic leukemia (CLL) cases by integrating whole-exome sequence and copy number to measure the fraction of cancer cells harboring each somatic mutation. We identified driver mutations as predominantly clonal (e.g., MYD88, trisomy 12, and del(13q)) or subclonal (e.g., SF3B1 and TP53), corresponding to earlier and later events in CLL evolution. We sampled leukemia cells from 18 patients at two time points. Ten of twelve CLL cases treated with chemotherapy (but only one of six without treatment) underwent clonal evolution, predominantly involving subclones with driver mutations (e.g., SF3B1 and TP53) that expanded over time. Furthermore, presence of a subclonal driver mutation was an independent risk factor for rapid disease progression. Our study thus uncovers patterns of clonal evolution in CLL, providing insights into its stepwise transformation, and links the presence of subclones with adverse clinical outcomes.

SF3B1 and Other Novel Cancer Genes in Chronic Lymphocytic Leukemia

BACKGROUND The somatic genetic basis of chronic lymphocytic leukemia, a common and clinically heterogeneous leukemia occurring in adults, remains poorly understood. METHODS We obtained DNA samples from leukemia cells in 91 patients with chronic lymphocytic leukemia and performed massively parallel sequencing of 88 whole exomes and whole genomes, together with sequencing of matched germline DNA, to characterize the spectrum of somatic mutations in this disease. RESULTS Nine genes that are mutated at significant frequencies were identified, including four with established roles in chronic lymphocytic leukemia (TP53 in 15% of patients, ATM in 9%, MYD88 in 10%, and NOTCH1 in 4%) and five with unestablished roles (SF3B1, ZMYM3, MAPK1, FBXW7, and DDX3X). SF3B1, which functions at the catalytic core of the spliceosome, was the second most frequently mutated gene (with mutations occurring in 15% of patients). SF3B1 mutations occurred primarily in tumors with deletions in chromosome 11q, which are associated with a poor prognosis in patients with chronic lymphocytic leukemia. We further discovered that tumor samples with mutations in SF3B1 had alterations in pre-messenger RNA (mRNA) splicing. CONCLUSIONS Our study defines the landscape of somatic mutations in chronic lymphocytic leukemia and highlights pre-mRNA splicing as a critical cellular process contributing to chronic lymphocytic leukemia.

sprouty encodes a novel antagonist of FGF signaling that patterns apical branching of the Drosophila airways.

Antagonists of several growth factor signaling pathways play important roles in developmental patterning by limiting the range of the cognate inducer. Here, we describe an antagonist of FGF signaling that patterns apical branching of the Drosophila airways. In wild-type embryos, the Branchless FGF induces secondary branching by activating the Breathless FGF receptor near the tips of growing primary branches. In sprouty mutants, the FGF pathway is overactive and ectopic branches are induced on the stalks of primary branches. We show that FGF signaling induces sprouty expression in the nearby tip cells, and sprouty acts nonautonomously and in a competitive fashion to block signaling to the more distant stalk cells. sprouty encodes a novel cysteine-rich protein that defines a new family of putative signaling molecules that may similarly function as FGF antagonists in vertebrate development.

CRISPR-Cas9 Knockin Mice for Genome Editing and Cancer Modeling

CRISPR-Cas9 is a versatile genome editing technology for studying the functions of genetic elements. To broadly enable the application of Cas9 in vivo, we established a Cre-dependent Cas9 knockin mouse. We demonstrated in vivo as well as ex vivo genome editing using adeno-associated virus (AAV)-, lentivirus-, or particle-mediated delivery of guide RNA in neurons, immune cells, and endothelial cells. Using these mice, we simultaneously modeled the dynamics of KRAS, p53, and LKB1, the top three significantly mutated genes in lung adenocarcinoma. Delivery of a single AAV vector in the lung generated loss-of-function mutations in p53 and Lkb1, as well as homology-directed repair-mediated Kras(G12D) mutations, leading to macroscopic tumors of adenocarcinoma pathology. Together, these results suggest that Cas9 mice empower a wide range of biological and disease modeling applications.

Single-cell transcriptomics reveals bimodality in expression and splicing in immune cells

Recent molecular studies have shown that, even when derived from a seemingly homogenous population, individual cells can exhibit substantial differences in gene expression, protein levels and phenotypic output, with important functional consequences. Existing studies of cellular heterogeneity, however, have typically measured only a few pre-selected RNAs or proteins simultaneously, because genomic profiling methods could not be applied to single cells until very recently. Here we use single-cell RNA sequencing to investigate heterogeneity in the response of mouse bone-marrow-derived dendritic cells (BMDCs) to lipopolysaccharide. We find extensive, and previously unobserved, bimodal variation in messenger RNA abundance and splicing patterns, which we validate by RNA-fluorescence in situ hybridization for select transcripts. In particular, hundreds of key immune genes are bimodally expressed across cells, surprisingly even for genes that are very highly expressed at the population average. Moreover, splicing patterns demonstrate previously unobserved levels of heterogeneity between cells. Some of the observed bimodality can be attributed to closely related, yet distinct, known maturity states of BMDCs; other portions reflect differences in the usage of key regulatory circuits. For example, we identify a module of 137 highly variable, yet co-regulated, antiviral response genes. Using cells from knockout mice, we show that variability in this module may be propagated through an interferon feedback circuit, involving the transcriptional regulators Stat2 and Irf7. Our study demonstrates the power and promise of single-cell genomics in uncovering functional diversity between cells and in deciphering cell states and circuits.

Peroxisomes are signaling platforms for antiviral innate immunity.

Peroxisomes have long been established to play a central role in regulating various metabolic activities in mammalian cells. These organelles act in concert with mitochondria to control the metabolism of lipids and reactive oxygen species. However, while mitochondria have emerged as an important site of antiviral signal transduction, a role for peroxisomes in immune defense is unknown. Here, we report that the RIG-I-like receptor (RLR) adaptor protein MAVS is located on peroxisomes and mitochondria. We find that peroxisomal and mitochondrial MAVS act sequentially to create an antiviral cellular state. Upon viral infection, peroxisomal MAVS induces the rapid interferon-independent expression of defense factors that provide short-term protection, whereas mitochondrial MAVS activates an interferon-dependent signaling pathway with delayed kinetics, which amplifies and stabilizes the antiviral response. The interferon regulatory factor IRF1 plays a crucial role in regulating MAVS-dependent signaling from peroxisomes. These results establish that peroxisomes are an important site of antiviral signal transduction.