scholarly journals Identification of spatial co-expression patterns and intra-tissue heterogeneity in spatially resolved transcriptomics by region-specific denoising

Author(s):  
Linhua Wang ◽  
Zhandong Liu

Abstract We are pleased to introduce a first-of-its-kind tool that combines in-silico region detection and missing value estimation for spatially resolved transcriptomics. Spatial transcriptomics by 10X Visium (ST) is a new technology used to dissect gene and cell spatial organization. Analyzing this new type of data has two main challenges: automatically annotating the major tissue regions and excessive zero values of gene expression due to high dropout rates. We developed a computational tool—MIST—that addresses both challenges by automatically identifying tissue regions and estimating missing gene-expression values for each detected region. We validated MIST detected regions across multiple datasets using manual annotation on the histological staining images as references. We also demonstrated that MIST can accurately recover ST’s missing values through hold-out experiments. Furthermore, we showed that MIST could identify intra-tissue heterogeneity and recover spatial gene-gene co-expression signals. We therefore strongly encourage using MIST before downstream ST analysis because it provides unbiased region annotations and enables accurately denoised spatial gene-expression profiles.

2021 ◽  
Author(s):  
Linhua Wang ◽  
Zhandong Liu

We are pleased to introduce a first–of–its–kind algorithm that combines in–silico region detection and spatial gene expression imputation. Spatial transcriptomics by 10X Visium (ST) is a new technology used to dissect gene and cell spatial organization. Analyzing this new type of data has two main challenges: automatically annotating the major tissue regions and excessive zero values of gene expression due to high dropout rates. We developed a computational tool—MIST—that addresses both challenges by automatically identifying tissue regions and estimating missing gene expression values for individual tissue regions. We validated MIST detected regions across multiple datasets using manual annotation on the histological staining images as references. We also demonstrated that MIST can accurately recover ST's missing values through hold–out experiments. Furthermore, we showed that MIST could identify subtle intra–tissue heterogeneity and recover spatial gene–gene interaction signals. We therefore strongly encourage using MIST prior to downstream ST analysis because it provides unbiased region annotations and enables accurately de–noised spatial gene expression profiles.


2021 ◽  
Author(s):  
Kangning Dong ◽  
Shihua Zhang

Recent advances in spatially resolved transcriptomics have enabled comprehensive measurements of gene expression patterns while retaining spatial context of tissue microenvironment. Deciphering the spatial context of spots in a tissue needs to use their spatial information carefully. To this end, we developed a graph attention auto- encoder framework STGATE to accurately identify spatial domains by learning low-dimensional latent embeddings via integrating spatial information and gene expression profiles. To better characterize the spatial similarity at the boundary of spatial domains, STGATE adopts an attention mechanism to adaptively learn the similarity of neighboring spots, and an optional cell type-aware module through integrating the pre-clustering of gene expressions. We validated STGATE on diverse spatial transcriptomics datasets generated by different platforms with different spatial resolutions. STGATE could substantially improve the identification accuracy of spatial domains, and denoise the data while preserving spatial expression patterns. Importantly, STGATE could be extended to multiple consecutive sections for reducing batch effects between sections and extracting 3D expression domains from the reconstructed 3D tissue effectively.


2008 ◽  
Vol 5 (2) ◽  
Author(s):  
Li Teng ◽  
Laiwan Chan

SummaryTraditional analysis of gene expression profiles use clustering to find groups of coexpressed genes which have similar expression patterns. However clustering is time consuming and could be diffcult for very large scale dataset. We proposed the idea of Discovering Distinct Patterns (DDP) in gene expression profiles. Since patterns showing by the gene expressions reveal their regulate mechanisms. It is significant to find all different patterns existing in the dataset when there is little prior knowledge. It is also a helpful start before taking on further analysis. We propose an algorithm for DDP by iteratively picking out pairs of gene expression patterns which have the largest dissimilarities. This method can also be used as preprocessing to initialize centers for clustering methods, like K-means. Experiments on both synthetic dataset and real gene expression datasets show our method is very effective in finding distinct patterns which have gene functional significance and is also effcient.


2005 ◽  
Vol 289 (4) ◽  
pp. L545-L553 ◽  
Author(s):  
Joseph Zabner ◽  
Todd E. Scheetz ◽  
Hakeem G. Almabrazi ◽  
Thomas L. Casavant ◽  
Jian Huang ◽  
...  

Cystic fibrosis (CF) is caused by mutations in the cystic fibrosis transmembrane conductance regulator (CFTR), an epithelial chloride channel regulated by phosphorylation. Most of the disease-associated morbidity is the consequence of chronic lung infection with progressive tissue destruction. As an approach to investigate the cellular effects of CFTR mutations, we used large-scale microarray hybridization to contrast the gene expression profiles of well-differentiated primary cultures of human CF and non-CF airway epithelia grown under resting culture conditions. We surveyed the expression profiles for 10 non-CF and 10 ΔF508 homozygote samples. Of the 22,283 genes represented on the Affymetrix U133A GeneChip, we found evidence of significant changes in expression in 24 genes by two-sample t-test ( P < 0.00001). A second, three-filter method of comparative analysis found no significant differences between the groups. The levels of CFTR mRNA were comparable in both groups. There were no significant differences in the gene expression patterns between male and female CF specimens. There were 18 genes with significant increases and 6 genes with decreases in CF relative to non-CF samples. Although the function of many of the differentially expressed genes is unknown, one transcript that was elevated in CF, the KCl cotransporter (KCC4), is a candidate for further study. Overall, the results indicate that CFTR dysfunction has little direct impact on airway epithelial gene expression in samples grown under these conditions.


2020 ◽  
Author(s):  
Alexander Calderwood ◽  
Jo Hepworth ◽  
Shannon Woodhouse ◽  
Lorelei Bilham ◽  
D. Marc Jones ◽  
...  

AbstractThe timing of the floral transition affects reproduction and yield, however its regulation in crops remains poorly understood. Here, we use RNA-Seq to determine and compare gene expression dynamics through the floral transition in the model species Arabidopsis thaliana and the closely related crop Brassica rapa. A direct comparison of gene expression over time between species shows little similarity, which could lead to the inference that different gene regulatory networks are at play. However, these differences can be largely resolved by synchronisation, through curve registration, of gene expression profiles. We find that different registration functions are required for different genes, indicating that there is no common ‘developmental time’ to which Arabidopsis and B. rapa can be mapped through gene expression. Instead, the expression patterns of different genes progress at different rates. We find that co-regulated genes show similar changes in synchronisation between species, suggesting that similar gene regulatory sub-network structures may be active with different wiring between them. A detailed comparison of the regulation of the floral transition between Arabidopsis and B. rapa, and between two B. rapa accessions reveals different modes of regulation of the key floral integrator SOC1, and that the floral transition in the B. rapa accessions is triggered by different pathways, even when grown under the same environmental conditions. Our study adds to the mechanistic understanding of the regulatory network of flowering time in rapid cycling B. rapa under long days and highlights the importance of registration methods for the comparison of developmental gene expression data.


2021 ◽  
Author(s):  
Anna E Backhaus ◽  
Ashleigh Lister ◽  
Melissa Tomkins ◽  
Nikolai M. Adamski ◽  
James Simmonds ◽  
...  

Spikelets are the fundamental building blocks of Poaceae inflorescences and their development and branching patterns determine the various inflorescence architectures and grain yield of grasses. In wheat, the central spikelets produce the most and largest grains, while spikelet size gradually decreases acro- and basipetally, giving rise to the characteristic lanceolate shape of wheat spikes. The acropetal gradient correlates with the developmental age of spikelets, however the basal spikelets are developed first and the cause of their small size and rudimentary development is unclear. Here, we adapted G&T-seq, a low-input transcriptomics approach, to characterise gene expression profiles within spatial sections of individual spikes before and after the establishment of the lanceolate shape. We observed larger differences in gene expression profiles between the apical, central and basal sections of a single spike than between any section belonging to consecutive developmental timepoints. We found that SVP MADS-box transcription factors, including VRT-A2, are expressed highest in the basal section of the wheat spike and display the opposite expression gradient to flowering E-class SEP1 genes. Based on multi-year field trials and transgenic lines we show that higher expression of VRT-A2 in the basal sections of the spike is associated with increased numbers of rudimentary basal spikelets. Our results, supported by computational modelling, suggest that the delayed transition of basal spikelets from vegetative to floral developmental programmes results in the lanceolate shape of wheat spikes. This study highlights the value of spatially resolved transcriptomics to gain new insights into developmental genetics pathways of grass inflorescences.


Author(s):  
Ana M Mesa ◽  
Jiude Mao ◽  
Theresa I Medrano ◽  
Nathan J Bivens ◽  
Alexander Jurkevich ◽  
...  

Abstract Histone proteins undergo various modifications that alter chromatin structure, including addition of methyl groups. Enhancer of homolog 2 (EZH2), is a histone methyltransferase that methylates lysine residue 27, and thereby, suppresses gene expression. EZH2 plays integral role in the uterus and other reproductive organs. We have previously shown that conditional deletion of uterine EZH2 results in increased proliferation of luminal and glandular epithelial cells, and RNAseq analyses reveal several uterine transcriptomic changes in Ezh2 conditional (c) knockout (KO) mice that can affect estrogen signaling pathways. To pinpoint the origin of such gene expression changes, we used the recently developed spatial transcriptomics (ST) method with the hypotheses that Ezh2cKO mice would predominantly demonstrate changes in epithelial cells and/or ablation of this gene would disrupt normal epithelial/stromal gene expression patterns. Uteri were collected from ovariectomized adult WT and Ezh2cKO mice and analyzed by ST. Asb4, Cxcl14, Dio2, and Igfbp5 were increased, Sult1d1, Mt3, and Lcn2 were reduced in Ezh2cKO uterine epithelium vs. WT epithelium. For Ezh2cKO uterine stroma, differentially expressed key hub genes included Cald1, Fbln1, Myh11, Acta2, and Tagln. Conditional loss of uterine Ezh2 also appears to shift the balance of gene expression profiles in epithelial vs. stromal tissue toward uterine epithelial cell and gland development and proliferation, consistent with uterine gland hyperplasia in these mice. Current findings provide further insight into how EZH2 may selectively affect uterine epithelial and stromal compartments. Additionally, these transcriptome data might provide the mechanistic understanding and valuable biomarkers for human endometrial disorders with epigenetic underpinnings.


Author(s):  
Crescenzio Gallo

The possible applications of modeling and simulation in the field of bioinformatics are very extensive, ranging from understanding basic metabolic paths to exploring genetic variability. Experimental results carried out with DNA microarrays allow researchers to measure expression levels for thousands of genes simultaneously, across different conditions and over time. A key step in the analysis of gene expression data is the detection of groups of genes that manifest similar expression patterns. In this chapter, the authors examine various methods for analyzing gene expression data, addressing the important topics of (1) selecting the most differentially expressed genes, (2) grouping them by means of their relationships, and (3) classifying samples based on gene expressions.


2020 ◽  
Vol 21 (22) ◽  
pp. 8520
Author(s):  
Ling Qin ◽  
Erying Chen ◽  
Feifei Li ◽  
Xiao Yu ◽  
Zhenyu Liu ◽  
...  

Foxtail millet (Setaria italica (L.) P. Beauv) is an important food and forage crop because of its health benefits and adaptation to drought stress; however, reports of transcriptomic analysis of genes responding to re-watering after drought stress in foxtail millet are rare. The present study evaluated physiological parameters, such as proline content, p5cs enzyme activity, anti-oxidation enzyme activities, and investigated gene expression patterns using RNA sequencing of the drought-tolerant foxtail millet variety (Jigu 16) treated with drought stress and rehydration. The results indicated that drought stress-responsive genes were related to many multiple metabolic processes, such as photosynthesis, signal transduction, phenylpropanoid biosynthesis, starch and sucrose metabolism, and osmotic adjustment. Furthermore, the Δ1-pyrroline-5-carboxylate synthetase genes, SiP5CS1 and SiP5CS2, were remarkably upregulated in foxtail millet under drought stress conditions. Foxtail millet can also recover well on rehydration after drought stress through gene regulation. Our data demonstrate that recovery on rehydration primarily involves proline metabolism, sugar metabolism, hormone signal transduction, water transport, and detoxification, plus reversal of the expression direction of most drought-responsive genes. Our results provided a detailed description of the comparative transcriptome response of foxtail millet variety Jigu 16 under drought and rehydration environments. Furthermore, we identify SiP5CS2 as an important gene likely involved in the drought tolerance of foxtail millet.


2019 ◽  
Vol 3 (Supplement_1) ◽  
Author(s):  
Zhi Chai ◽  
Yafei Lyu ◽  
Qiuyan Chen ◽  
Cheng-Hsin Wei ◽  
Lindsay Snyder ◽  
...  

Abstract Objectives To characterize and compare the impact of vitamin A (VA) deficiency on gene expression patterns in the small intestine (SI) and the colon, and to discover novel target genes in VA-related biological pathways. Methods vitamin A deficient (VAD) mice were generated by feeding VAD diet to pregnant C57/BL6 dams and their post-weaning offspring. Total mRNA extracted from SI and colon were sequenced using Illumina HiSeq 2500 platform. Differentially Expressed Gene (DEG), Gene Ontology (GO) enrichment, and Weighted Gene Co-expression Network Analysis (WGCNA) were performed to characterize expression patterns and co-expression patterns. Results The comparison between vitamin A sufficient (VAS) and VAD groups detected 49 and 94 DEGs in SI and colon, respectively. According to GO information, DEGs in the SI demonstrated significant enrichment in categories relevant to retinoid metabolic process, molecule binding, and immune function. Immunity related pathways, such as “humoral immune response” and “complement activation,” were positively associated with VA in SI. On the contrary, in colon, “cell division” was the only enriched category and was negatively associated with VA. WGCNA identified modules significantly correlated with VA status in SI and in colon. One of those modules contained five known retinoic acid targets. Therefore we have prioritized the other module members (e.g., Mbl2, Mmp9, Mmp13, Cxcl14 and Pkd1l2) to be investigated as candidate genes regulated by VA. Comparison of co-expression modules between SI and colon indicated distinct VA effects on these two organs. Conclusions The results show that VA deficiency alters the gene expression profiles in SI and colon quite differently. Some immune-related genes (Mbl2, Mmp9, Mmp13, Cxcl14 and Pkd1l2) may be novel targets under the control of VA in SI. Funding Sources NIH training grant and NIH research grant. Supporting Tables, Images and/or Graphs


Sign in / Sign up

Export Citation Format

Share Document