PreTSA: computationally efficient modeling of temporal and spatial gene expression patterns.
Ontology highlight
ABSTRACT:
SUBMITTER: Zhuang H
PROVIDER: S-EPMC10996487 | biostudies-literature | 2024 Mar
REPOSITORIES: biostudies-literature
Publications
PreTSA: computationally efficient modeling of temporal and spatial gene expression patterns.
Zhuang Haotian H Ji Zhicheng Z
bioRxiv : the preprint server for biology 20240325
Modeling temporal and spatial gene expression patterns in large-scale single-cell and spatial transcriptomics data is a computationally intensive task. We present PreTSA, a method that offers computational efficiency in modeling these patterns and is applicable to single-cell and spatial transcriptomics data comprising millions of cells. PreTSA consistently matches the results of state-of-the-art methods while significantly reducing computational time. PreTSA provides a unique solution for study ...[more]
Similar Datasets
Project description:The ongoing global change and the increased interest in macroecological processes call for the analysis of spatially extensive data on species communities to understand and forecast distributional changes of biodiversity. Recently developed joint species distribution models can deal with numerous species efficiently, while explicitly accounting for spatial structure in the data. However, their applicability is generally limited to relatively small spatial data sets because of their severe computational scaling as the number of spatial locations increases. In this work, we propose a practical alleviation of this scalability constraint for joint species modeling by exploiting two spatial-statistics techniques that facilitate the analysis of large spatial data sets: Gaussian predictive process and nearest-neighbor Gaussian process. We devised an efficient Gibbs posterior sampling algorithm for Bayesian model fitting that allows us to analyze community data sets consisting of hundreds of species sampled from up to hundreds of thousands of spatial units. The performance of these methods is demonstrated using an extensive plant data set of 30,955 spatial units as a case study. We provide an implementation of the presented methods as an extension to the hierarchical modeling of species communities framework.
| S-EPMC7027487 | biostudies-literature
Project description:Thousands of Pseudomonas aeruginosa RNA sequencing (RNA-seq) gene expression profiles are publicly available via the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA). In this work, the transcriptional profiles from hundreds of studies performed by over 75 research groups were reanalyzed in aggregate to create a powerful tool for hypothesis generation and testing. Raw sequence data were uniformly processed using the Salmon pseudoaligner, and this read mapping method was validated by comparison to a direct alignment method. We developed filtering criteria to exclude samples with aberrant levels of housekeeping gene expression or an unexpected number of genes with no reported values and normalized the filtered compendia using the ratio-of-medians method. The filtering and normalization steps greatly improved gene expression correlations for genes within the same operon or regulon across the 2,333 samples. Since the RNA-seq data were generated using diverse strains, we report the effects of mapping samples to noncognate reference genomes by separately analyzing all samples mapped to cDNA reference genomes for strains PAO1 and PA14, two divergent strains that were used to generate most of the samples. Finally, we developed an algorithm to incorporate new data as they are deposited into the SRA. Our processing and quality control methods provide a scalable framework for taking advantage of the troves of biological information hibernating in the depths of microbial gene expression data and yield useful tools for P. aeruginosa RNA-seq data to be leveraged for diverse research goals. IMPORTANCE Pseudomonas aeruginosa is a causative agent of a wide range of infections, including chronic infections associated with cystic fibrosis. These P. aeruginosa infections are difficult to treat and often have negative outcomes. To aid in the study of this problematic pathogen, we mapped, filtered for quality, and normalized thousands of P. aeruginosa RNA-seq gene expression profiles that were publicly available via the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA). The resulting compendia facilitate analyses across experiments, strains, and conditions. Ultimately, the workflow that we present could be applied to analyses of other microbial species.
| S-EPMC9948711 | biostudies-literature
Project description:Allen Brain Atlas (ABA) provides a valuable resource of spatial/temporal gene expressions in mammalian brains. Despite rich information extracted from this database, current analyses suffer from several limitations. First, most studies are either gene-centric or region-centric, thus are inadequate to capture the superposition of multiple spatial-temporal patterns. Second, standard tools of expression analysis such as matrix factorization can capture those patterns but do not explicitly incorporate spatial dependency. To overcome those limitations, we proposed a computational method to detect recurrent patterns in the spatial-temporal gene expression data of developing mouse brains. We demonstrated that regional distinction in brain development could be revealed by localized gene expression patterns. The patterns expressed in the forebrain, medullary and pontomedullary, and basal ganglia are enriched with genes involved in forebrain development, locomotory behavior, and dopamine metabolism respectively. In addition, the timing of global gene expression patterns reflects the general trends of molecular events in mouse brain development. Furthermore, we validated functional implications of the inferred patterns by showing genes sharing similar spatial-temporal expression patterns with Lhx2 exhibited differential expression in the embryonic forebrains of Lhx2 mutant mice. These analysis outcomes confirm the utility of recurrent expression patterns in studying brain development.
| S-EPMC4726224 | biostudies-literature
Project description:BackgroundThe sea urchin is a basal deuterostome that is more closely related to vertebrates than many organisms traditionally used to study neurogenesis. This phylogenetic position means that the sea urchin can provide insights into the evolution of the nervous system by helping resolve which developmental processes are deuterostome innovations, which are innovations in other clades, and which are ancestral. However, the nervous system of echinoderms is one of the least understood of all major metazoan phyla. To gain insights into echinoderm neurogenesis, spatial and temporal gene expression data are essential. Then, functional data will enable the building of a detailed gene regulatory network for neurogenesis in the sea urchin that can be compared across metazoans to resolve questions about how nervous systems evolved.ResultsHere, we analyze spatiotemporal gene expression during sea urchin neurogenesis for genes that have been shown to be neurogenic in one or more species. We report the expression of 21 genes expressed in areas of neurogenesis in the sea urchin embryo from blastula stage (just before neural progenitors begin their specification sequence) through pluteus larval stage (when much of the nervous system has been patterned). Among those 21 gene expression patterns, we report expression of 11 transcription factors and 2 axon guidance genes, each expressed in discrete domains in the neuroectoderm or in the endoderm. Most of these genes are expressed in and around the ciliary band. Some including the transcription factors Lv-mbx, Lv-dmrt, Lv-islet, and Lv-atbf1, the nuclear protein Lv-prohibitin, and the guidance molecule Lv-semaa are expressed in the endoderm where they are presumably involved in neurogenesis in the gut.ConclusionsThis study builds a foundation to study how neurons are specified and evolved by analyzing spatial and temporal gene expression during neurogenesis in a basal deuterostome. With these expression patterns, we will be able to understand what genes are required for neural development in the sea urchin. These data can be used as a starting point to (1) build a spatial gene regulatory network for sea urchin neurogenesis, (2) identify how subtypes of neurons are specified, (3) perform comparative studies with the sea urchin, protostome, and vertebrate organisms.
| S-EPMC6371548 | biostudies-literature
Project description:MotivationIn the pursuits of mechanistic understanding of cell differentiation, it is often necessary to compare multiple differentiation processes triggered by different external stimuli and internal perturbations. Available methods for comparing temporal gene expression patterns are limited to a gene-by-gene approach, which ignores co-expression information and thus is sensitive to measurement noise.MethodsWe present a method for co-expression network based comparison of temporal expression patterns (NACEP). NACEP compares the temporal patterns of a gene between two experimental conditions, taking into consideration all of the possible co-expression modules that this gene may participate in. The NACEP program is available at http://biocomp.bioen.uiuc.edu/nacep.ResultsWe applied NACEP to analyze retinoid acid (RA)-induced differentiation of embryonic stem (ES) cells. The analysis suggests that RA may facilitate neural differentiation by inducing the shh and insulin receptor pathways. NACEP was also applied to compare the temporal responses of seven RNA inhibition (RNAi) experiments. As proof of concept, we demonstrate that the difference in the temporal responses to RNAi treatments can be used to derive interaction relationships of transcription factors (TFs), and therefore infer regulatory modules within a transcription network. In particular, the analysis suggested a novel regulatory relationship between two pluripotency regulators, Esrrb and Tbx3, which was supported by in vivo binding of Esrrb to the promoter of Tbx3.AvailabilityThe NACEP program and the supplementary documents are available at http://biocomp.bioen.uiuc.edu/nacep.Contactszhong@illinois.eduSupplementary informationSupplementary data are available at Bioinformatics online.
| S-EPMC2982154 | biostudies-literature
Project description:The aim of the present study was to generate hypotheses on the involvement of uncharacterized genes in biological processes. To this end, supervised learning was used to analyze microarray-derived time-series gene expression data. Our method was objectively evaluated on known genes using cross-validation and provided high-precision Gene Ontology biological process classifications for 211 of the 213 uncharacterized genes in the data set used. In addition, new roles in biological process were hypothesized for known genes. Our method uses biological knowledge expressed by Gene Ontology and generates a rule model associating this knowledge with minimal characteristic features of temporal gene expression profiles. This model allows learning and classification of multiple biological process roles for each gene and can predict participation of genes in a biological process even though the genes of this class exhibit a wide variety of gene expression profiles including inverse coregulation. A considerable number of the hypothesized new roles for known genes were confirmed by literature search. In addition, many biological process roles hypothesized for uncharacterized genes were found to agree with assumptions based on homology information. To our knowledge, a gene classifier of similar scope and functionality has not been reported earlier.
| S-EPMC430886 | biostudies-literature
Project description:The transcriptional mechanism through which chondrocytes control the spatial and temporal composition of the cartilage tissue has remained largely elusive. The central aim of this study was to identify whether transcriptional enhancers played a role in the organisation of the chondrocytes in cartilaginous tissue. We focused on the Aggrecan gene (Acan) as it is essential for the normal structure and function of cartilage and it is expressed developmentally in different stages of chondrocyte maturation. Using transgenic reporter studies in mice we identified four elements, two of which showed individual chondrocyte developmental stage specificity. In particular, one enhancer (-80) distinguishes itself from the others by being predominantly active in adult cartilage. Furthermore, the -62 element uniquely drove reporter activity in early chondrocytes. The remaining chondrocyte specific enhancers, +28 and -30, showed no preference to chondrocyte type. The transcription factor SOX9 interacted with all the enhancers in vitro and mutation of SOX9 binding sites in one of the enhancers (-30) resulted in a loss of its chondrocyte specificity and ectopic enhancer reporter activity. Thus, the Acan enhancers orchestrate the precise spatiotemporal expression of this gene in cartilage types at different stages of development and adulthood.
| S-EPMC5772622 | biostudies-literature
Project description:Gene expression in bacteria is a remarkably controlled and intricate process impacted by many factors. One such factor is the genomic position of a gene within a bacterial genome. Genes located near the origin of replication generally have a higher expression level, increased dosage, and are often more conserved than genes located farther from the origin of replication. The majority of the studies involved with these findings have only noted this phenomenon in a single gene or cluster of genes that was re-located to pre-determined positions within a bacterial genome. In this work, we look at the overall expression levels from eleven bacterial data sets from Escherichia coli, Bacillus subtilis, Streptomyces, and Sinorhizobium meliloti. We have confirmed that gene expression tends to decrease when moving away from the origin of replication in majority of the replicons analysed in this study. This study sheds light on the impact of genomic location on molecular trends such as gene expression and highlights the importance of accounting for spatial trends in bacterial molecular analysis.
| S-EPMC7324424 | biostudies-literature
Project description:Spatial transcriptomics (ST) provides critical insights into the complex spatial organization of gene expression in tissues, enabling researchers to unravel the intricate relationship between cellular environments and biological function. Identifying spatial domains within tissues is essential for understanding tissue architecture and the mechanisms underlying various biological processes, including development and disease progression. Here, we present Randomized Spatial PCA (RASP), a novel spatially aware dimensionality reduction method for spatial transcriptomics (ST) data. RASP is designed to be orders-of-magnitude faster than existing techniques, scale to ST data with hundreds-of-thousands of locations, support the flexible integration of non-transcriptomic covariates, and enable the reconstruction of de-noised and spatially smoothed expression values for individual genes. To achieve these goals, RASP uses a randomized two-stage principal component analysis (PCA) framework that leverages sparse matrix operations and configurable spatial smoothing. We compared the performance of RASP against five alternative methods (BASS, GraphST, SEDR, spatialPCA, and STAGATE) on four publicly available ST datasets generated using diverse techniques and resolutions (10x Visium, Stereo-Seq, MERFISH, and 10x Xenium) on human and mouse tissues. Our results demonstrate that RASP achieves tissue domain detection performance comparable or superior to existing methods with a several orders-of-magnitude improvement in computational speed. The efficiency of RASP enhances the analysis of complex ST data by facilitating the exploration of increasingly high-resolution subcellular ST datasets that are being generated.
| S-EPMC11703158 | biostudies-literature
Project description:The C. elegans genome has been completely sequenced, and the developmental anatomy of this model organism is described at single-cell resolution. Here we utilize strategies that exploit this precisely defined architecture to link gene expression to cell type. We obtained RNAs from specific cells and from each developmental stage using tissue-specific promoters to mark cells for isolation by FACS or for mRNA extraction by the mRNA-tagging method. We then generated gene expression profiles of more than 30 different cells and developmental stages using tiling arrays. Machine-learning-based analysis detected transcripts corresponding to established gene models and revealed novel transcriptionally active regions (TARs) in noncoding domains that comprise at least 10% of the total C. elegans genome. Our results show that about 75% of transcripts with detectable expression are differentially expressed among developmental stages and across cell types. Examination of known tissue- and cell-specific transcripts validates these data sets and suggests that newly identified TARs may exercise cell-specific functions. Additionally, we used self-organizing maps to define groups of coregulated transcripts and applied regulatory element analysis to identify known transcription factor- and miRNA-binding sites, as well as novel motifs that likely function to control subsets of these genes. By using cell-specific, whole-genome profiling strategies, we have detected a large number of novel transcripts and produced high-resolution gene expression maps that provide a basis for establishing the roles of individual genes in cellular differentiation.
| S-EPMC3032935 | biostudies-literature