Bioinformatics approach to spatially resolved transcriptomics

Author(s):  
Ivan Krešimir Lukić

Spatially resolved transcriptomics encompasses a growing number of methods developed to enable gene expression profiling of individual cells within a tissue. Different technologies are available and they vary with respect to: the method used to define regions of interest, the method used to assess gene expression, and resolution. Since techniques based on next-generation sequencing are the most prevalent, and provide single-cell resolution, many bioinformatics tools for spatially resolved data are shared with single-cell RNA-seq. The analysis pipelines diverge at the level of quantification matrix, downstream of which spatial techniques require specific tools to answer key biological questions. Those questions include: (i) cell type classification; (ii) detection of genes with specific spatial distribution; (iii) identification of novel tissue regions based on gene expression patterns; (iv) cell–cell interactions. On the other hand, analysis of spatially resolved data is burdened by several specific challenges. Defining regions of interest, e.g. neoplastic tissue, often calls for manual annotation of images, which then poses a bottleneck in the pipeline. Another specific issue is the third spatial dimension and the need to expand the analysis beyond a single slice. Despite the problems, it can be predicted that the popularity of spatial techniques will keep growing until they replace single-cell assays (which will remain limited to specific cases, like blood). As soon as the computational protocol reach the maturity (e.g. bulk RNA-seq), one can foresee the expansion of spatial techniques beyond basic or translational research, even into routine medical diagnostics.

2020 ◽  
Author(s):  
Timothy J. Durham ◽  
Riza M. Daza ◽  
Louis Gevirtzman ◽  
Darren A. Cusanovich ◽  
William Stafford Noble ◽  
...  

AbstractRecently developed single cell technologies allow researchers to characterize cell states at ever greater resolution and scale. C. elegans is a particularly tractable system for studying development, and recent single cell RNA-seq studies characterized the gene expression patterns for nearly every cell type in the embryo and at the second larval stage (L2). Gene expression patterns are useful for learning about gene function and give insight into the biochemical state of different cell types; however, in order to understand these cell types, we must also determine how these gene expression levels are regulated. We present the first single cell ATAC-seq study in C. elegans. We collected data in L2 larvae to match the available single cell RNA-seq data set, and we identify tissue-specific chromatin accessibility patterns that align well with existing data, including the L2 single cell RNA-seq results. Using a novel implementation of the latent Dirichlet allocation algorithm, we leverage the single-cell resolution of the sci-ATAC-seq data to identify accessible loci at the level of individual cell types, providing new maps of putative cell type-specific gene regulatory sites, with promise for better understanding of cellular differentiation and gene regulation in the worm.


2021 ◽  
Vol 5 (Supplement_1) ◽  
pp. A753-A754
Author(s):  
Jingfei Chen ◽  
Jialei Duan ◽  
Alina Montalbano ◽  
Boxun Li ◽  
Gary Hon ◽  
...  

Abstract Single-cell RNA-seq of Mouse Decidual Leukocytes Reveals Intriguing Gestational Changes in the Immune Cell Landscape and Effects of SRC-1/-2 Double-DeficiencyOur previous findings suggest that the fetus signals the mother when it is ready to be born through secretion of surfactant components/immune modulators, surfactant protein A (SP-A) and platelet-activating factor (PAF) by the fetal lung into amniotic fluid (AF). This occurs with increased proinflammatory cytokine expression by fetal AF macrophages (Mϕ), increased myometrial NF-κB activation and contractile (CAP) gene expression. Steroid receptor coactivator (Src)-1 and Src-2 are critical for developmental induction SP-A and PAF by the fetal lung. The finding that pregnant wild-type (WT) mice carrying Src-1 and Src-2 double-deficient fetuses (Src-1/-2d/d) manifested a marked delay (~38h) in parturition, further suggests that signals for parturition arise from the fetus and that Src-1 and Src-2 serve a critical role. Infiltrating leukocytes at the maternal-fetal interface (MFI)/decidua are known to play a central role in pregnancy maintenance and parturition timing. However, there is limited knowledge regarding gestational changes in immune cell composition toward term. To analyze gestational changes in the composition of immune cells within decidua and effects of Src-1/-2d/d, we conducted single-cell RNA-seq of ~17,000 decidual leukocytes (CD45+) from pregnant mice at 15.5 and 18.5 dpc carrying WT or Src-1/-2d/d fetuses. Unsupervised clustering identified 22 distinct decidual immune cell clusters, comprised of Mϕ, B cells, natural killer cells, neutrophils, dendritic cells and monocytes. Significant differences in cell type composition and transcriptional profiles were found across study groups (WT @ 15.5 dpc vs. 18.5 dpc; Src-1/-2d/dvs. WT @15.5 dpc and 18.5 dpc). Interestingly, in deciduas of pregnant mice carrying WT fetuses, Mϕ and B cell clusters markedly increased between 15.5 and 18.5 dpc, whereas, neutrophils declined; however, these gestational changes did not occur in pregnant mice carrying Src-1/-2d/d fetuses. Differential gene expression and gene ontology enrichment analyses revealed specific gene expression patterns distinguishing these immune cell subtypes to uncover their putative functions. These findings highlight the complexity and dynamics of the decidual immune cell landscape during the transition from myometrial quiescence to contractility, and the fetal-maternal interactions leading to parturition. Support: NIH grants R01-HL050022 (CRM) and P01-HD087150 (CRM), Burroughs Wellcome Preterm Birth Grant #1019823 (CRM).


2017 ◽  
Author(s):  
Garth R. Ilsley ◽  
Ritsuko Suyama ◽  
Takeshi Noda ◽  
Nori Satoh ◽  
Nicholas M. Luscombe

AbstractSingle-cell RNA-seq has been established as a reliable and accessible technique enabling new types of analyses, such as identifying cell types and studying spatial and temporal gene expression variation and change at single-cell resolution. Recently, single-cell RNA-seq has been applied to developing embryos, which offers great potential for finding and characterising genes controlling the course of development along with their expression patterns. In this study, we applied single-cell RNA-seq to the 16-cell stage of the Ciona embryo, a marine chordate and performed a computational search for cell-specific gene expression patterns. We recovered many known expression patterns from our single-cell RNA-seq data and despite extensive previous screens, we succeeded in finding new cell-specific patterns, which we validated by in situ and single-cell qPCR.


2018 ◽  
Author(s):  
Michael L. Mucenski ◽  
Robert Mahoney ◽  
Mike Adam ◽  
Andrew S. Potter ◽  
S. Steven Potter

AbstractThe uterus is a remarkable organ that must guard against infections while maintaining the ability to support growth of a fetus without rejection. The Hoxa10 and Hoxa11 genes have previously been shown to play essential roles in uterus development and function. In this report we show that the Hoxc9,10,11 genes play a redundant role in the formation of uterine glands. In addition, we use single cell RNA-seq to create a high resolution gene expression atlas of the developing wild type mouse uterus. Cell types and subtypes are defined, for example dividing endothelial cells into arterial, venous, capillary, and lymphatic, while epithelial cells separate into luminal and glandular subtypes. Further, a surprising heterogeneity of stromal and myocyte cell types are identified. Transcription factor codes and ligand/receptor interactions are characterized. We also used single cell RNA-seq to globally define the altered gene expression patterns in all developing uterus cell types for two Hox mutants, with 8 or 9 mutant Hox genes. The mutants show a striking disruption of Wnt signaling as well as the Cxcl12/Cxcr4 ligand/receptor axis.Summary statementA single cell RNA-seq study of the developing mouse uterus defines cellular heterogeneities, lineage specific gene expression programs and perturbed pathways in Hox9,10,11 mutants.


2020 ◽  
Author(s):  
Ping An ◽  
Paul G. Cantalupo ◽  
Wenshan Zheng ◽  
Maria Teresa Sáenz-Robles ◽  
Alexis M. Duray ◽  
...  

BKV is a human polyomavirus that is generally harmless but can cause devastating disease in immunosuppressed individuals. BKV infection of renal cells is a common problem for kidney transplant patients undergoing immunosuppressive therapy. In cultured primary human renal proximal tubule epithelial cells (RPTE), BKV undergoes a productive infection. The BKV-encoded large T antigen (LT) induces cell cycle entry, resulting in the upregulation of numerous genes associated with cell proliferation. Consistently, microarray and RNA-seq experiments performed in bulk infected cell populations identified several proliferation-related pathways that are upregulated by BKV. These studies revealed few genes that are downregulated. In this study, we analyzed viral and cellular transcripts in single mock or BKV-infected cells. We found that the levels of viral mRNAs vary widely among infected cells, resulting in different levels of LT and viral capsid protein expression. Cells expressing the highest levels of viral transcripts account for approximately 20% of the culture and have a gene expression pattern that is distinct from cells expressing lower levels of viral mRNAs. Surprisingly, cells expressing low levels of viral mRNA do not progress with time to high expression, suggesting that the two cellular responses are determined prior to or shortly following infection. Finally, comparison of cellular gene expression patterns of cells expressing high levels of viral mRNA with mock-infected cells, or with cells expressing low levels of viral mRNA, revealed previously unidentified pathways that are downregulated by BKV. Among these are pathways associated with drug metabolism and detoxification, TNF-signaling, energy metabolism, and translation. IMPORTANCE The outcome of viral infection is determined by the ability of the virus to redirect cellular systems towards progeny production countered by the ability of the cell to block these viral actions. Thus, an infected culture consists of thousands of cells, each fighting their own individual battle. Bulk measurements, such as PCR or RNA-seq, measure the average of these individual responses to infection. Single cell transcriptomics provides a window to the one-on-one battle between BKV and each cell. Our studies reveal that only a minority of infected cells are overwhelmed by the virus and produce large amounts of BKV mRNAs and proteins, while the infection appears to be restricted in the remaining cells. Correlation of viral transcript levels with cellular gene expression patterns reveals pathways manipulated by BKV that may play a role in limiting infection.


2020 ◽  
Vol 3 (4) ◽  
pp. 72
Author(s):  
Anupama Prakash ◽  
Antónia Monteiro

Butterflies are well known for their beautiful wings and have been great systems to understand the ecology, evolution, genetics, and development of patterning and coloration. These color patterns are mosaics on the wing created by the tiling of individual units called scales, which develop from single cells. Traditionally, bulk RNA sequencing (RNA-seq) has been used extensively to identify the loci involved in wing color development and pattern formation. RNA-seq provides an averaged gene expression landscape of the entire wing tissue or of small dissected wing regions under consideration. However, to understand the gene expression patterns of the units of color, which are the scales, and to identify different scale cell types within a wing that produce different colors and scale structures, it is necessary to study single cells. This has recently been facilitated by the advent of single-cell sequencing. Here, we provide a detailed protocol for the dissociation of cells from Bicyclus anynana pupal wings to obtain a viable single-cell suspension for downstream single-cell sequencing. We outline our experimental design and the use of fluorescence-activated cell sorting (FACS) to obtain putative scale-building and socket cells based on size. Finally, we discuss some of the current challenges of this technique in studying single-cell scale development and suggest future avenues to address these challenges.


2021 ◽  
Vol 22 (24) ◽  
pp. 13623
Author(s):  
Braulio Valdebenito-Maturana ◽  
Cristina Guatimosim ◽  
Mónica Alejandra Carrasco ◽  
Juan Carlos Tapia

Spatial transcriptomics (ST) is transforming the way we can study gene expression and its regulation through position-specific resolution within tissues. However, as in bulk RNA-Seq, transposable elements (TEs) are not being studied due to their highly repetitive nature. In recent years, TEs have been recognized as important regulators of gene expression, and thus, TE expression analysis in a spatially resolved manner could further help to understand their role in gene regulation within tissues. We present SpatialTE, a tool to analyze TE expression from ST datasets and show its application in somatic and diseased tissues. The results indicate that TEs have spatially regulated expression patterns and that their expression profiles are spatially altered in ALS disease, indicating that TEs might perform differential regulatory functions within tissue organs. We have made SpatialTE publicly available as open-source software under an MIT license.


Genes ◽  
2021 ◽  
Vol 12 (10) ◽  
pp. 1640
Author(s):  
Zhipeng Li ◽  
Xinhui Song ◽  
Shan Yin ◽  
Jiageng Yan ◽  
Peiru Lv ◽  
...  

Donkeys are an important domesticated animal, providing labor, meat, milk, and medicinal materials for humans. However, the donkey population is continuously declining and even at risk of extinction. The application of modern animal production technology, such as oocyte in vitro maturation, is a promising method to improve the donkey population. In this study, we explore the gene expression patterns of donkey germinal vesicle (GV) and in vitro matured metaphase II (MII) oocytes using single cell RNA-seq of the candidate genes along with the regulatory mechanisms that affect donkey oocyte maturation. We identified a total of 24,164 oocyte genes of which 9073 were significant differentially expressed in the GV and MII oocytes. Further Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis indicated that these genes were associated with the meiotic cell cycle, mitochondrion activity, and N-glycan biosynthesis, which might be the key genes and regulatory mechanisms affecting the maturation of donkey oocytes. Our study provides considerable understanding regarding the maturation of donkey oocytes and serves as a theoretical basis for improving the development of donkey oocytes, which could ultimately benefit the expansion of the donkey population and conservation of biodiversity and genetic resources.


2021 ◽  
Vol 22 (S3) ◽  
Author(s):  
Yuanyuan Li ◽  
Ping Luo ◽  
Yi Lu ◽  
Fang-Xiang Wu

Abstract Background With the development of the technology of single-cell sequence, revealing homogeneity and heterogeneity between cells has become a new area of computational systems biology research. However, the clustering of cell types becomes more complex with the mutual penetration between different types of cells and the instability of gene expression. One way of overcoming this problem is to group similar, related single cells together by the means of various clustering analysis methods. Although some methods such as spectral clustering can do well in the identification of cell types, they only consider the similarities between cells and ignore the influence of dissimilarities on clustering results. This methodology may limit the performance of most of the conventional clustering algorithms for the identification of clusters, it needs to develop special methods for high-dimensional sparse categorical data. Results Inspired by the phenomenon that same type cells have similar gene expression patterns, but different types of cells evoke dissimilar gene expression patterns, we improve the existing spectral clustering method for clustering single-cell data that is based on both similarities and dissimilarities between cells. The method first measures the similarity/dissimilarity among cells, then constructs the incidence matrix by fusing similarity matrix with dissimilarity matrix, and, finally, uses the eigenvalues of the incidence matrix to perform dimensionality reduction and employs the K-means algorithm in the low dimensional space to achieve clustering. The proposed improved spectral clustering method is compared with the conventional spectral clustering method in recognizing cell types on several real single-cell RNA-seq datasets. Conclusions In summary, we show that adding intercellular dissimilarity can effectively improve accuracy and achieve robustness and that improved spectral clustering method outperforms the traditional spectral clustering method in grouping cells.


Genes ◽  
2021 ◽  
Vol 12 (2) ◽  
pp. 311
Author(s):  
Zhenqiu Liu

Single-cell RNA-seq (scRNA-seq) is a powerful tool to measure the expression patterns of individual cells and discover heterogeneity and functional diversity among cell populations. Due to variability, it is challenging to analyze such data efficiently. Many clustering methods have been developed using at least one free parameter. Different choices for free parameters may lead to substantially different visualizations and clusters. Tuning free parameters is also time consuming. Thus there is need for a simple, robust, and efficient clustering method. In this paper, we propose a new regularized Gaussian graphical clustering (RGGC) method for scRNA-seq data. RGGC is based on high-order (partial) correlations and subspace learning, and is robust over a wide-range of a regularized parameter λ. Therefore, we can simply set λ=2 or λ=log(p) for AIC (Akaike information criterion) or BIC (Bayesian information criterion) without cross-validation. Cell subpopulations are discovered by the Louvain community detection algorithm that determines the number of clusters automatically. There is no free parameter to be tuned with RGGC. When evaluated with simulated and benchmark scRNA-seq data sets against widely used methods, RGGC is computationally efficient and one of the top performers. It can detect inter-sample cell heterogeneity, when applied to glioblastoma scRNA-seq data.


Sign in / Sign up

Export Citation Format

Share Document