CCST: Cell clustering for spatial transcriptomics data with graph neural network

Author(s):  
Jiachen Li ◽  
Siheng Chen ◽  
Xiaoyong Pan ◽  
Ye Yuan ◽  
Hong-bin Shen

Abstract Spatial transcriptomics data can provide high-throughput gene expression profiling and spatial structure of tissues simultaneously. An essential question of its initial analysis is cell clustering. However, most existing studies rely on only gene expression information and cannot utilize spatial information efficiently. Taking advantages of two recent technical development, spatial transcriptomics and graph neural network, we thus introduce CCST, Cell Clustering for Spatial Transcriptomics data with graph neural network, an unsupervised cell clustering method based on graph convolutional network to improve ab initio cell clustering and discovering of novel sub cell types based on curated cell category annotation. CCST is a general framework for dealing with various kinds of spatially resolved transcriptomics. With application to five in vitro and in vivo spatial datasets, we show that CCST outperforms other spatial cluster approaches on spatial transcriptomics datasets, and can clearly identify all four cell cycle phases from MERFISH data of cultured cells, and find novel functional sub cell types with different micro-environments from seqFISH+ data of brain, which are all validated experimentally, inspiring novel biological hypotheses about the underlying interactions among cell state, cell type and micro-environment.

1992 ◽  
Vol 12 (3) ◽  
pp. 1202-1208
Author(s):  
R A Graves ◽  
P Tontonoz ◽  
B M Spiegelman

The molecular basis of adipocyte-specific gene expression is not well understood. We have previously identified a 518-bp enhancer from the adipocyte P2 gene that stimulates adipose-specific gene expression in both cultured cells and transgenic mice. In this analysis of the enhancer, we have defined and characterized a 122-bp DNA fragment that directs differentiation-dependent gene expression in cultured preadipocytes and adipocytes. Several cis-acting elements have been identified and shown by mutational analysis to be important for full enhancer activity. One pair of sequences, ARE2 and ARE4, binds a nuclear factor (ARF2) present in extracts derived from many cell types. Multiple copies of these elements stimulate gene expression from a minimal promoter in preadipocytes, adipocytes, and several other cultured cell lines. A second pair of elements, ARE6 and ARE7, binds a separate factor (ARF6) that is detected only in nuclear extracts derived from adipocytes. The ability of multimers of ARE6 or ARE7 to stimulate promoter activity is strictly adipocyte specific. Mutations in the ARE6 sequence greatly reduce the activity of the 518-bp enhancer. These data demonstrate that several cis- and trans-acting components contribute to the activity of the adipocyte P2 enhancer and suggest that ARF6, a novel differentiation-dependent factor, may be a key regulator of adipogenic gene expression.


1988 ◽  
Vol 90 (1) ◽  
pp. 73-77
Author(s):  
A. Harris ◽  
L. Coleman

The establishment of a tissue-culture system for epithelial cells derived from human foetal pancreas has recently been reported. Further analyses have now been made on these cells in vitro, together with parallel investigation of the distribution of different cell types within the intact foetal pancreas. Results support the view that the cultured cells are ductal in origin and nature. Pancreatic epithelial cell cultures have also been established from foetuses with cystic fibrosis.


2019 ◽  
Author(s):  
Ugur M. Ayturk ◽  
Joseph P. Scollan ◽  
Alexander Vesprey ◽  
Christina M. Jacobsen ◽  
Paola Divieti Pajevic ◽  
...  

ABSTRACTSingle cell RNA-seq (scRNA-seq) is emerging as a powerful technology to examine transcriptomes of individual cells. We determined whether scRNA-seq could be used to detect the effect of environmental and pharmacologic perturbations on osteoblasts. We began with a commonly used in vitro system in which freshly isolated neonatal mouse calvarial cells are expanded and induced to produce a mineralized matrix. We used scRNA-seq to compare the relative cell type abundances and the transcriptomes of freshly isolated cells to those that had been cultured for 12 days in vitro. We observed that the percentage of macrophage-like cells increased from 6% in freshly isolated calvarial cells to 34% in cultured cells. We also found that Bglap transcripts were abundant in freshly isolated osteoblasts but nearly undetectable in the cultured calvarial cells. Thus, scRNA-seq revealed significant differences between heterogeneity of cells in vivo and in vitro. We next performed scRNA-seq on freshly recovered long bone endocortical cells from mice that received either vehicle or Sclerostin-neutralizing antibody for 1 week. Bone anabolism-associated transcripts were also not significantly increased in immature and mature osteoblasts recovered from Sclerostin-neutralizing antibody treated mice; this is likely a consequence of being underpowered to detect modest changes in gene expression, since only 7% of the sequenced endocortical cells were osteoblasts, and a limited portion of their transcriptomes were sampled. We conclude that scRNA-seq can detect changes in cell abundance, identity, and gene expression in skeletally derived cells. In order to detect modest changes in osteoblast gene expression at the single cell level in the appendicular skeleton, larger numbers of osteoblasts from endocortical bone are required.


2019 ◽  
Author(s):  
Marcus Alvarez ◽  
Elior Rahmani ◽  
Brandon Jew ◽  
Kristina M. Garske ◽  
Zong Miao ◽  
...  

AbstractSingle-nucleus RNA sequencing (snRNA-seq) measures gene expression in individual nuclei instead of cells, allowing for unbiased cell type characterization in solid tissues. Contrary to single-cell RNA seq (scRNA-seq), we observe that snRNA-seq is commonly subject to contamination by high amounts of extranuclear background RNA, which can lead to identification of spurious cell types in downstream clustering analyses if overlooked. We present a novel approach to remove debris-contaminated droplets in snRNA-seq experiments, called Debris Identification using Expectation Maximization (DIEM). Our likelihood-based approach models the gene expression distribution of debris and cell types, which are estimated using EM. We evaluated DIEM using three snRNA-seq data sets: 1) human differentiating preadipocytes in vitro, 2) fresh mouse brain tissue, and 3) human frozen adipose tissue (AT) from six individuals. All three data sets showed various degrees of extranuclear RNA contamination. We observed that existing methods fail to account for contaminated droplets and led to spurious cell types. When compared to filtering using these state of the art methods, DIEM better removed droplets containing high levels of extranuclear RNA and led to higher quality clusters. Although DIEM was designed for snRNA-seq data, we also successfully applied DIEM to single-cell data. To conclude, our novel method DIEM removes debris-contaminated droplets from single-cell-based data fast and effectively, leading to cleaner downstream analysis. Our code is freely available for use at https://github.com/marcalva/diem.


2015 ◽  
Vol 137 (5) ◽  
Author(s):  
Jason Lee ◽  
Aaron B. Baker

In vitro systems for applying mechanical strain to cultured cells are commonly used to investigate cellular mechanotransduction pathways in a variety of cell types. These systems often apply mechanical forces to a flexible membrane on which cells are cultured. A consequence of the motion of the membrane in these systems is the generation of flow and the unintended application of shear stress to the cells. We recently described a flexible system for applying mechanical strain to cultured cells, which uses a linear motor to drive a piston array to create biaxial strain within multiwell culture plates. To better understand the fluidic stresses generated by this system and other systems of this type, we created a computational fluid dynamics model to simulate the flow during the mechanical loading cycle. Alterations in the frequency or maximal strain magnitude led to a linear increase in the average fluid velocity within the well and a nonlinear increase in the shear stress at the culture surface over the ranges tested (0.5–2.0 Hz and 1–10% maximal strain). For all cases, the applied shear stresses were relatively low and on the order of millipascal with a dynamic waveform having a primary and secondary peak in the shear stress over a single mechanical strain cycle. These findings should be considered when interpreting experimental results using these devices, particularly in the case when the cell type used is sensitive to low magnitude, oscillatory shear stresses.


2020 ◽  
Vol 35 (Supplement_3) ◽  
Author(s):  
Julie Williams ◽  
Sanlin Robinson ◽  
Babak Alaei ◽  
Kimberly Homan ◽  
Maryam Clausen ◽  
...  

Abstract Background and Aims Questions abound regarding the translation of in vitro 2D cell culture systems to the human setting. This is especially true of the kidney in which there is a complex hierarchical structure and a multitude of cell types. While it is well accepted that extracellular matrix plays a large part in directing cellular physiology emerging research has highlighted the importance of shear stresses and flow rates too. To fully recapitulate the normal gene expression and function of a particular renal cell type how important is it to completely reconstitute their in vivo surroundings? Method To answer this question, we have cultured proximal tubular (PT) epithelial cells in a 3-dimensional channel embedded within an engineered extracellular matrix (ECM) under physiological flow that is colocalised with an adjacent channel lined with renal microvascular endothelial cells that mimic a peritubular capillary. Modifications to the system were made to allow up to 12 chips to be run in parallel in an easily handleable form. After a period of maturation under continuous flow, both cell types were harvested for RNAseq analyses. RNA expression data was compared with cells cultured under static 2-dimensional conditions on plastic or the engineered ECM. Additionally, the perfusion of glucose through this 3D vascularised PT model has been investigated in the presence and absence of known diabetes modulating agents. Results PCA of RNAseq data showed that a) static non-coated, b) static matrix-coated and c) flow matrix-coated conditions separated into 3 distinct groups, while cell co-culture had less impact. Analysis of transcriptomic signatures showed that many genes were modulated by the matrix with additional genes influenced under flow conditions. Several of these genes, classified as transporters, are of particular importance when using this model to assess drug uptake and safety implications. Co-culture regulated some interesting genes, but fewer than anticipated. Preliminary experiments are underway to monitor glucose uptake and transport between tubules under different conditions. Conclusion We have developed a medium throughput system in which matrix and flow modulate gene expression. This system can be used to study the physiology of molecular cross-talk between cells. Ongoing analysis will further consider relevance to human physiology.


2007 ◽  
Vol 88 (1) ◽  
pp. 105-113 ◽  
Author(s):  
Daniel Doucet ◽  
Anic Levasseur ◽  
Catherine Béliveau ◽  
Renée Lapointe ◽  
Don Stoltz ◽  
...  

Polydnaviruses (PDVs) are dsDNA viruses transmitted by ichneumonid and braconid endoparasitoids to their lepidopteran hosts during oviposition. Wasp carriers are asymptomatic and transmit the virus to their progeny through the germ line; replication is confined to the calyx region of the wasp ovary, where the virus accumulates in the fluid bathing the eggs. In the lepidopteran host, however, no virus replication takes place, but PDV gene expression is essential for successful parasitism. Sustained gene expression in the absence of virus replication thus requires that the circular PDV genome segments persist for days within host cells. Available evidence suggests that most genome segments persist as episomes, but recent studies have indicated that some genome segments may undergo integration within lepidopteran genomic DNA, at least in vitro. In the present study, an integrated form of a Tranosema rostrale ichnovirus (TrIV) genome segment was cloned from genomic DNA extracted from infected Choristoneura fumiferana CF-124T cells and junction regions on either side of the viral DNA sequence were sequenced. This is the first proven example of integration of an ichnovirus genome segment in infected lepidopteran cells. Interestingly, circular forms of this genome segment do not appear to persist in these cells; none the less, a gene (TrFrep1) carried by this genome segment displays long-term transcription in infected cultured cells.


F1000Research ◽  
2017 ◽  
Vol 6 ◽  
pp. 1208 ◽  
Author(s):  
Michael J. Garabedian ◽  
Charles A. Harris ◽  
Freddy Jeanneteau

Glucocorticoids via the glucocorticoid receptor (GR) have effects on a variety of cell types, eliciting important physiological responses via changes in gene expression and signaling. Although decades of research have illuminated the mechanism of how this important steroid receptor controls gene expression using in vitro and cell culture–based approaches, how GR responds to changes in external signals in vivo under normal and pathological conditions remains elusive. The goal of this review is to highlight recent work on GR action in fat cells and liver to affect metabolism in vivo and the role GR ligands and receptor phosphorylation play in calibrating signaling outputs by GR in the brain in health and disease. We also suggest that both the brain and fat tissue communicate to affect physiology and behavior and that understanding this “brain-fat axis” will enable a more complete understanding of metabolic diseases and inform new ways to target them.


2016 ◽  
Vol 69 (3) ◽  
pp. 349 ◽  
Author(s):  
Yucheng Liu ◽  
Shufeng Li ◽  
Liandong Feng ◽  
Hao Yu ◽  
Xiaoliang Qi ◽  
...  

Poly(β-amino ester)s (PBAEs) have been proved to effectively transfer DNA to various cell types. However, PBAEs with high molecular weights also show considerable toxicities, partly resulting from inadequate degradation of their polyester backbone. In this study, we created novel poly(β-amino ester)s (SF-1, 2, 3, and 4; notation SFs refers to all the four polymers) which were characterised by the cleavable disulfide bonds. Moreover, a new technique, termed magnetofection that uses magnetic nanoparticles to enhance gene expression, has recently been well developed. The negatively charged magnetic nanoparticles (MNPs) with good biocompatibility in vitro were prepared here to subsequently combine with SFs and DNA via electrostatic interaction, leading to the formation of the magnetic gene complexes MNP/SFs/DNA. 3-(4,5-Dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide assays and transfection experiments were performed in A549 cells to investigate all the resulting complexes. Studies indicated that the synthesised PBAEs exhibited good biodegradation and regulated release of DNA as a result of the reductive cleavage of the disulfide bonds, giving higher transfection efficiency along with much lower cytotoxicity compared with commercially available transfection agent polyethylenimine (Mw 25 kDa). Furthermore, when MNP was involved at a MNP/DNA weight ratio of 0.5, the magnetic gene complexes MNP/SFs/DNA showed enhanced levels of gene expression while maintaining low cytotoxicity.


Reproduction ◽  
2010 ◽  
Vol 139 (4) ◽  
pp. 749-757 ◽  
Author(s):  
Gillian Cowan ◽  
Andrew J Childs ◽  
Richard A Anderson ◽  
Philippa T K Saunders

The somatic (Sertoli cell (SC), Leydig cell (LC), and peritubular myoid (PTM) cell) cells play key roles in development of the fetal testis. We established monolayer cultures from second trimester human testes and investigated the pattern of expression of cell-lineage characteristic mRNAs. Expression of some SC-associated genes (SRY, SOX9, WT1, GATA4, and SF1) was detectable up to and including passage 3 (P3), while others (anti-Müllerian hormone; desert hedgehog) present prior to dissociation were not expressed in the cultured cells. Transcripts encoding the androgen receptor were expressed but addition of dihydrotestosterone (DHT) had no impact on expression of mRNAs expressed in SC or LC. Total concentrations of mRNAs encoding smooth muscle actin (ACTA2) and desmin increased from P1 to P3; an increasing proportion of the cells in the cultures were immunopositive for ACTA2 consistent with proliferation/differentiation of PTM cells. In conclusion, somatic cell monolayer cultures were established from human fetal testes; these cultures could form the basis for future studies based on isolation of purified populations of somatic cells and manipulation of gene expression that is difficult to achieve with organ culture systems. Our results suggest that fetal SC do not maintain a fully differentiated phenotype in vitro, yet PTM (ACTA2 positive) cells readily adapt to monolayer culture conditions in the presence of DHT. This culture system provides an opportunity to study the impact of regulatory factors on gene expression in PTM cells, a population thought to play a key role in mediating androgen action within the developing testis.


Sign in / Sign up

Export Citation Format

Share Document