RESUMO
Jointly profiling the transcriptome, chromatin accessibility and other molecular properties of single cells offers a powerful way to study cellular diversity. Here we present MultiVI, a probabilistic model to analyze such multiomic data and leverage it to enhance single-modality datasets. MultiVI creates a joint representation that allows an analysis of all modalities included in the multiomic input data, even for cells for which one or more modalities are missing. It is available at scvi-tools.org .
Assuntos
Modelos Estatísticos , TranscriptomaRESUMO
MOTIVATION: Gene regulatory networks define regulatory relationships between transcription factors and target genes within a biological system, and reconstructing them is essential for understanding cellular growth and function. Methods for inferring and reconstructing networks from genomics data have evolved rapidly over the last decade in response to advances in sequencing technology and machine learning. The scale of data collection has increased dramatically; the largest genome-wide gene expression datasets have grown from thousands of measurements to millions of single cells, and new technologies are on the horizon to increase to tens of millions of cells and above. RESULTS: In this work, we present the Inferelator 3.0, which has been significantly updated to integrate data from distinct cell types to learn context-specific regulatory networks and aggregate them into a shared regulatory network, while retaining the functionality of the previous versions. The Inferelator is able to integrate the largest single-cell datasets and learn cell-type-specific gene regulatory networks. Compared to other network inference methods, the Inferelator learns new and informative Saccharomyces cerevisiae networks from single-cell gene expression data, measured by recovery of a known gold standard. We demonstrate its scaling capabilities by learning networks for multiple distinct neuronal and glial cell types in the developing Mus musculus brain at E18 from a large (1.3 million) single-cell gene expression dataset with paired single-cell chromatin accessibility data. AVAILABILITY AND IMPLEMENTATION: The inferelator software is available on GitHub (https://github.com/flatironinstitute/inferelator) under the MIT license and has been released as python packages with associated documentation (https://inferelator.readthedocs.io/). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.
Assuntos
Redes Reguladoras de Genes , Software , Animais , Camundongos , Genômica , Genoma , CromatinaRESUMO
The analysis of single-cell genomics data presents several statistical challenges, and extensive efforts have been made to produce methods for the analysis of this data that impute missing values, address sampling issues and quantify and correct for noise. In spite of such efforts, no consensus on best practices has been established and all current approaches vary substantially based on the available data and empirical tests. The k-Nearest Neighbor Graph (kNN-G) is often used to infer the identities of, and relationships between, cells and is the basis of many widely used dimensionality-reduction and projection methods. The kNN-G has also been the basis for imputation methods using, e.g., neighbor averaging and graph diffusion. However, due to the lack of an agreed-upon optimal objective function for choosing hyperparameters, these methods tend to oversmooth data, thereby resulting in a loss of information with regard to cell identity and the specific gene-to-gene patterns underlying regulatory mechanisms. In this paper, we investigate the tuning of kNN- and diffusion-based denoising methods with a novel non-stochastic method for optimally preserving biologically relevant informative variance in single-cell data. The framework, Denoising Expression data with a Weighted Affinity Kernel and Self-Supervision (DEWÄKSS), uses a self-supervised technique to tune its parameters. We demonstrate that denoising with optimal parameters selected by our objective function (i) is robust to preprocessing methods using data from established benchmarks, (ii) disentangles cellular identity and maintains robust clusters over dimension-reduction methods, (iii) maintains variance along several expression dimensions, unlike previous heuristic-based methods that tend to oversmooth data variance, and (iv) rarely involves diffusion but rather uses a fixed weighted kNN graph for denoising. Together, these findings provide a new understanding of kNN- and diffusion-based denoising methods. Code and example data for DEWÄKSS is available at https://gitlab.com/Xparx/dewakss/-/tree/Tjarnberg2020branch.
Assuntos
Algoritmos , Genômica/métodos , Análise de Célula Única/métodos , Aprendizado de Máquina Supervisionado , Animais , Linhagem Celular , Bases de Dados Genéticas , Humanos , Camundongos , RNA-Seq , Saccharomyces cerevisiaeRESUMO
The cardinal classes are a useful simplification of cortical interneuron diversity, but such broad subgroupings gloss over the molecular, morphological, and circuit specificity of interneuron subtypes, most notably among the somatostatin interneuron class. Although there is evidence that this diversity is functionally relevant, the circuit implications of this diversity are unknown. To address this knowledge gap, we designed a series of genetic strategies to target the breadth of somatostatin interneuron subtypes and found that each subtype possesses a unique laminar organization and stereotyped axonal projection pattern. Using these strategies, we examined the afferent and efferent connectivity of three subtypes (two Martinotti and one non-Martinotti) and demonstrated that they possess selective connectivity with intratelecephalic or pyramidal tract neurons. Even when two subtypes targeted the same pyramidal cell type, their synaptic targeting proved selective for particular dendritic compartments. We thus provide evidence that subtypes of somatostatin interneurons form cell-type-specific cortical circuits.
Assuntos
Interneurônios , Neurônios , Interneurônios/fisiologia , Neurônios/fisiologia , Células Piramidais/fisiologia , Axônios/metabolismo , Somatostatina/metabolismo , Parvalbuminas/metabolismoRESUMO
Understanding how gene expression programs are controlled requires identifying regulatory relationships between transcription factors and target genes. Gene regulatory networks are typically constructed from gene expression data acquired following genetic perturbation or environmental stimulus. Single-cell RNA sequencing (scRNAseq) captures the gene expression state of thousands of individual cells in a single experiment, offering advantages in combinatorial experimental design, large numbers of independent measurements, and accessing the interaction between the cell cycle and environmental responses that is hidden by population-level analysis of gene expression. To leverage these advantages, we developed a method for scRNAseq in budding yeast (Saccharomyces cerevisiae). We pooled diverse transcriptionally barcoded gene deletion mutants in 11 different environmental conditions and determined their expression state by sequencing 38,285 individual cells. We benchmarked a framework for learning gene regulatory networks from scRNAseq data that incorporates multitask learning and constructed a global gene regulatory network comprising 12,228 interactions.
Organisms switch their genes on and off to adapt to changing environments. This takes place thanks to complex networks of regulators that control which genes are actively 'read' by the cell to create the RNA molecules that are needed at the time. Piecing together these networks is key to fully understand the inner workings of living organisms, and how to potentially modify or artificially create them. Single-cell RNA sequencing is a powerful new tool that can measure which genes are turned on (or 'expressed') in an individual cell. Datasets with millions of gene expression profiles for individual cells now exist for organisms such as mice or humans. Yet, it is difficult to use these data to reconstruct networks of regulators; this is partly because scientists are not sure if the computational methods normally used to build these networks also work for single-cell RNA sequencing data. One way to check if this is the case is to use the methods on single-cell datasets from organisms where the networks of regulators are already known, and check whether the computational tools help to reach the same conclusion. Unfortunately, the regulatory networks in the organisms for which scientists have a lot of single-cell RNA sequencing data are still poorly known. There are living beings in which the networks are well characterised such as yeast but it has been difficult to do single-cell sequencing in them at the scale seen in other organisms. Jackson, Castro et al. first adapted a system for single-cell sequencing so that it would work in yeast. This generated a gene expression dataset of over 40,000 yeast cells. They then used a computational method (called the Inferelator) on these data to construct networks of regulators, and the results showed that the method performed well. This allowed Jackson, Castro et al. to start mapping how different networks connect, for example those that control the response to the environment and cell division. This is one of the benefits of single-cell RNA methods: cell division for example is not a process that can be examined at the level of a population, since the cells may all be at different life stages. In the future, the dataset will also be useful to scientists to benchmark a variety of single cell computational tools.
Assuntos
Código de Barras de DNA Taxonômico , Redes Reguladoras de Genes , Genótipo , Saccharomyces cerevisiae/genética , Análise de Sequência de RNA/métodos , Análise de Célula Única/métodos , Deleção de Genes , Regulação Fúngica da Expressão Gênica , Genes Fúngicos , Saccharomyces cerevisiae/metabolismo , Fatores de Transcrição/metabolismoRESUMO
A fundamental impediment to understanding the brain is the availability of inexpensive and robust methods for targeting and manipulating specific neuronal populations. The need to overcome this barrier is pressing because there are considerable anatomical, physiological, cognitive and behavioral differences between mice and higher mammalian species in which it is difficult to specifically target and manipulate genetically defined functional cell types. In particular, it is unclear the degree to which insights from mouse models can shed light on the neural mechanisms that mediate cognitive functions in higher species, including humans. Here we describe a novel recombinant adeno-associated virus that restricts gene expression to GABAergic interneurons within the telencephalon. We demonstrate that the viral expression is specific and robust, allowing for morphological visualization, activity monitoring and functional manipulation of interneurons in both mice and non-genetically tractable species, thus opening the possibility to study GABAergic function in virtually any vertebrate species.