Sepal: Identifying Transcript Profiles with Spatial Patterns by Diffusion-based Modeling.
Ontology highlight
ABSTRACT: Motivation
Collection of spatial signals in large numbers has become a routine task in multiple omicsfields, but parsing of these rich data sets still pose certain challenges. In whole or near-full transcriptome spatial techniques, spurious expression profiles are intermixed with those exhibiting an organized structure. To distinguish profiles with spatial patterns from the background noise, a metric that enables quantification of spatial structure is desirable. Current methods designed for similar purposes tend to be built around a framework of statistical hypothesis testing, hence we were compelled to explore a fundamentally different strategy.Results
We propose an unexplored approach to analyze spatial transcriptomics data, simulating diffusion of individual transcripts to extract genes with spatial patterns. The method performed as expected when presented with synthetic data. When applied to real data, it identified genes with distinct spatial profiles, involved in key biological processes or characteristic for certain cell types. Compared to existing methods, ours seemed to be less informed by the genes' expression levels and showed better time performance when run with multiple cores.Availability
Open-source Python package with a command line interface (CLI), freely available at https://github.com/almaan/sepal under a MIT licence. A mirror of the GitHub repository can be found at Zenodo, doi: 10.5281/zenodo.4573237.Supplementary information
Supplementary data are available at Bioinformatics online.
SUBMITTER: Anderson A
PROVIDER: S-EPMC8428601 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:The food-based empirical dietary index for hyperinsulinemia (EDIH) score assesses the insulinemic potential of diet. This cross-sectional study evaluated associations between EDIH scores from food frequency questionnaires with c-peptide concentrations and with 448 metabolites, from fasting plasma samples, in multivariable linear regression analyses. Metabolites were measured with liquid chromatography tandem mass spectroscopy. Using a robust two-stage study design, discovery of metabolite associations was conducted among 1109 Women's Health Initiative (WHI) Hormone Therapy (HT) trial participants and results replicated in an independent dataset of 810 WHI Observational Study (OS) participants. In both discovery and replication datasets, statistical significance was based on the false-discovery rate adjusted P < 0.05. In the multivariable-adjusted analyses, EDIH was significantly associated with c-peptide concentrations among 919 women (HT & OS) with c-peptide data. On average, c-peptide concentrations were 18% higher (95% CI, 6%, 32%; P-trend < 0.0001) in EDIH quintile 5 compared to quintile 1. Twenty-six metabolites were significantly associated with EDIH in the discovery dataset, and 19 of these were replicated in the validation dataset. Nine metabolites were found to decrease in abundance with increasing EDIH scores and included: C14:0 CE, C16:1 CE, C18:1 CE, C18:3 CE, C20:3 CE, C20:5 CE, C36:1 PS plasmalogen, trigonelline, and eicosapentanoate, whereas the 10 metabolites observed to increase with increasing EDIH scores were: C18:2 SM, C36:3 DAG, C36:4 DAG-A, C51:3 TAG, C52:3 TAG, C52:4, TAG, C54:3 TAG, C54:4 TAG, C54:6 TAG, and C10:2 carnitine. Cholesteryl esters, phospholipids, acylglycerols, and acylcarnitines may constitute circulating metabolites that are associated with insulinemic dietary patterns.
| S-EPMC6630814 | biostudies-other
Project description:BackgroundA key issue in the analysis of many spatial processes is the choice of an appropriate scale for the analysis. Smaller geographical units are generally preferable for the study of human phenomena because they are less likely to cause heterogeneous groups to be conflated. However, it can be harder to obtain data for small units and small-number problems can frustrate quantitative analysis. This research presents a new approach that can be used to estimate the most appropriate scale at which to aggregate point data to areas.Data and methodsThe proposed method works by creating a number of regular grids with iteratively smaller cell sizes (increasing grid resolution) and estimating the similarity between two realisations of the point pattern at each resolution. The method is applied first to simulated point patterns and then to real publicly available crime data from the city of Vancouver, Canada. The crime types tested are residential burglary, commercial burglary, theft from vehicle and theft of bike.FindingsThe results provide evidence for the size of spatial unit that is the most appropriate for the different types of crime studied. Importantly, the results are dependent on both the number of events in the data and the degree of spatial clustering, so a single 'appropriate' scale is not identified. The method is nevertheless useful as a means of better estimating what spatial scale might be appropriate for a particular piece of analysis.
| S-EPMC6594600 | biostudies-literature
Project description:ObjectivesDiffusion tensor imaging (DTI) tractography has shown tract-specific pathology in temporal lobe epilepsy (TLE). This technique normally yields a single value per diffusion parameter per tract, potentially reducing the sensitivity for the detection of focal changes. Our goal was to spatially characterize diffusion abnormalities of fasciculi carrying temporal lobe connections.MethodsWe studied 30 patients with drug-resistant TLE and 21 healthy control subjects. Twenty-four patients underwent DTI toward the end of video-EEG telemetry, with an average of 50 ± 54 hours between the last seizure and DTI examination. After manual dissection of the uncinate and inferior longitudinal and arcuate bundle, they were spatially matched based on their distance to the temporal lobe, providing between-subject correspondence of tract segments. We evaluated point-wise differences in diffusion parameters along each tract at group and subject levels.ResultsOur approach localized increased mean diffusivity restricted to or more prominent within the ipsilateral temporal lobe. These abnormalities tapered off as tracts exited the temporal lobe. We observed that the shorter the interval between the last seizure and DTI, the higher the mean diffusivity (MD) of the ipsilateral tracts. Linear discriminant analysis of tract segments correctly lateralized 87% of patients.ConclusionsThe centrifugal pattern of white matter diffusion abnormalities probably reflects astrogliosis and microstructure derangement related to seizure activity in the vicinity of the focus. The negative correlation between the interval from last seizure and MD suggests a role for postictal vasogenic edema. The ability to assess tracts segmentally may contribute to a better understanding of the extent of white matter pathology in epilepsy and assist in the presurgical evaluation of patients with TLE, particularly those with unremarkable conventional imaging results.
| S-EPMC3405250 | biostudies-literature
Project description:Codon usage bias (CUB) refers to the phenomena that synonymous codons are used in different frequencies in most genes and organisms. The general assumption is that codon biases reflect a balance between mutational biases and natural selection. Today we understand that the codon content is related and can affect all gene expression steps. Starting from the 1980s, codon-based indices have been used for answering different questions in all biomedical fields, including systems biology, agriculture, medicine, and biotechnology. In general, codon usage bias indices weigh each codon or a small set of codons to estimate the fitting of a certain coding sequence to a certain phenomenon (e.g., bias in codons, adaptation to the tRNA pool, frequencies of certain codons, transcription elongation speed, etc.) and are usually easy to implement. Today there are dozens of such indices; thus, this paper aims to review and compare the different codon usage bias indices, their applications, and advantages. In addition, we perform analysis that demonstrates that most indices tend to correlate even though they aim to capture different aspects. Due to the centrality of codon usage bias on different gene expression steps, it is important to keep developing new indices that can capture additional aspects that are not modeled with the current indices.
| S-EPMC8122159 | biostudies-literature
Project description:Nitric oxide (NO) is a gaseous signaling molecule that plays an important role in neurovascular coupling. NO produced by neurons diffuses into the smooth muscle surrounding cerebral arterioles, driving vasodilation. However, the rate of NO degradation in hemoglobin is orders of magnitude higher than in brain tissue, though how this might impact NO signaling dynamics is not completely understood. We used simulations to investigate how the spatial and temporal patterns of NO generation and degradation impacted dilation of a penetrating arteriole in cortex. We found that the spatial location of NO production and the size of the vessel both played an important role in determining its responsiveness to NO. The much higher rate of NO degradation and scavenging of NO in the blood relative to the tissue drove emergent vascular dynamics. Large vasodilation events could be followed by post-stimulus constrictions driven by the increased degradation of NO by the blood, and vasomotion-like 0.1-0.3 Hz oscillations could also be generated. We found that these dynamics could be enhanced by elevation of free hemoglobin in the plasma, which occurs in diseases such as malaria and sickle cell anemia, or following blood transfusions. Finally, we show that changes in blood flow during hypoxia or hyperoxia could be explained by altered NO degradation in the parenchyma. Our simulations suggest that many common vascular dynamics may be emergent phenomena generated by NO degradation by the blood or parenchyma.
| S-EPMC7410342 | biostudies-literature
Project description:BACKGROUND: The Saccharopolyspora erythraea genome sequence was released in 2007. In order to look at the gene regulations at whole transcriptome level, an expression microarray was specifically designed on the S. erythraea strain NRRL 2338 genome sequence. Based on these data, we set out to investigate the potential transcriptional regulatory networks and their organization. METHODOLOGY/PRINCIPAL FINDINGS: In view of the hierarchical structure of bacterial transcriptional regulation, we constructed a hierarchical coexpression network at whole transcriptome level. A total of 27 modules were identified from 1255 differentially expressed transcript units (TUs) across time course, which were further classified in to four groups. Functional enrichment analysis indicated the biological significance of our hierarchical network. It was indicated that primary metabolism is activated in the first rapid growth phase (phase A), and secondary metabolism is induced when the growth is slowed down (phase B). Among the 27 modules, two are highly correlated to erythromycin production. One contains all genes in the erythromycin-biosynthetic (ery) gene cluster and the other seems to be associated with erythromycin production by sharing common intermediate metabolites. Non-concomitant correlation between production and expression regulation was observed. Especially, by calculating the partial correlation coefficients and building the network based on Gaussian graphical model, intrinsic associations between modules were found, and the association between those two erythromycin production-correlated modules was included as expected. CONCLUSIONS: This work created a hierarchical model clustering transcriptome data into coordinated modules, and modules into groups across the time course, giving insight into the concerted transcriptional regulations especially the regulation corresponding to erythromycin production of S. erythraea. This strategy may be extendable to studies on other prokaryotic microorganisms.
| S-EPMC2920828 | biostudies-literature
Project description:The spatial organization in the cell nucleus is tightly linked to genome functions such as gene regulation. Similarly, specific spatial arrangements of biological components such as macromolecular complexes, organelles and cells are involved in many biological functions. Spatial interactions among elementary components of biological systems define their relative positioning and are key determinants of spatial patterns. However, biological variability and the lack of appropriate spatial statistical methods and models limit our current ability to analyze these interactions. Here, we developed a framework to dissect spatial interactions and organization principles by combining unbiased statistical tests, multiple spatial descriptors and new spatial models. We used plant constitutive heterochromatin as a model system to demonstrate the potential of our framework. Our results challenge the common view of a peripheral organization of chromocenters, showing that chromocenters are arranged along both radial and lateral directions in the nuclear space and obey a multiscale organization with scale-dependent antagonistic effects. The proposed generic framework will be useful to identify determinants of spatial organizations and to question their interplay with biological functions.
| S-EPMC7801681 | biostudies-literature
Project description:Alternative mRNA processing mechanisms lead to multiple transcripts (i.e. splice isoforms) of a given gene which may have distinct biological functions. Microarrays like Affymetrix GeneChips measure mRNA expression of genes using sets of nucleotide probes. Until recently probe sets were not designed for transcript specificity. Nevertheless, the re-analysis of established microarray data using newly defined transcript-specific probe sets may provide information about expression levels of specific transcripts.In the present study alignment of probe sequences of the Affymetrix microarray HG-U133A with Ensembl transcript sequences was performed to define transcript-specific probe sets. Out of a total of 247,965 perfect match probes, 95,008 were designated "transcript-specific", i.e. showing complete sequence alignment, no cross-hybridization, and transcript-, not only gene-specificity. These probes were grouped into 7,941 transcript-specific probe sets and 15,619 gene-specific probe sets, respectively. The former were used to differentiate 445 alternative transcripts of 215 genes. For selected transcripts, predicted by this analysis to be differentially expressed in the human kidney, confirmatory real-time RT-PCR experiments were performed. First, the expression of two specific transcripts of the genes PPM1A (PP2CA_HUMAN and P35813) and PLG (PLMN_HUMAN and Q5TEH5) in human kidneys was determined by the transcript-specific array analysis and confirmed by real-time RT-PCR. Secondly, disease-specific differential expression of single transcripts of PLG and ABCA1 (ABCA1_HUMAN and Q5VYS0_HUMAN) was computed from the available array data sets and confirmed by transcript-specific real-time RT-PCR.Transcript-specific analysis of microarray experiments can be employed to study gene-regulation on the transcript level using conventional microarray data. In this study, predictions based on sufficient probe set size and fold-change are confirmed by independent means.
| S-EPMC2650090 | biostudies-literature
Project description:Relatively little is known about reliability of longitudinal diffusion-tensor imaging (DTI) measurements despite growing interest in using DTI to track change in white matter structure. The purpose of this study is to quantify within- and between session scan-rescan reliability of DTI-derived measures that are commonly used to describe the characteristics of neural white matter in the context of neural plasticity research. DTI data were acquired from 16 cognitively healthy older adults (mean age 68.4). We used the Tract-Based Spatial Statistics (TBSS) approach implemented in FSL, evaluating how different DTI preprocessing choices affect reliability indices. Test-Retest reliability, quantified as ICC averaged across the voxels of the TBSS skeleton, ranged from 0.524 to 0.798 depending on the specific DTI-derived measure and the applied preprocessing steps. The two main preprocessing steps that we found to improve TBSS reliability were (a) the use of a common individual template and (b) smoothing DTI data using a 1-voxel median filter. Overall our data indicate that small choices in the preprocessing pipeline have a significant effect on test-retest reliability, therefore influencing the power to detect change within a longitudinal study. Furthermore, differences in the data processing pipeline limit the comparability of results across studies.
| S-EPMC6869410 | biostudies-literature
Project description:BackgroundFusion transcripts are involved in tumourigenesis and play a crucial role in tumour heterogeneity, tumour evolution and cancer treatment resistance. However, fusion transcripts have not been studied at high spatial resolution in tissue sections due to the lack of full-length transcripts with spatial information. New high-throughput technologies like spatial transcriptomics measure the transcriptome of tissue sections on almost single-cell level. While this technique does not allow for direct detection of fusion transcripts, we show that they can be inferred using the relative poly(A) tail abundance of the involved parental genes.MethodWe present a new method STfusion, which uses spatial transcriptomics to infer the presence and absence of poly(A) tails. A fusion transcript lacks a poly(A) tail for the 5' gene and has an elevated number of poly(A) tails for the 3' gene. Its expression level is defined by the upstream promoter of the 5' gene. STfusion measures the difference between the observed and expected number of poly(A) tails with a novel C-score.ResultsWe verified the STfusion ability to predict fusion transcripts on HeLa cells with known fusions. STfusion and C-score applied to clinical prostate cancer data revealed the spatial distribution of the cis-SAGe SLC45A3-ELK4 in 12 tissue sections with almost single-cell resolution. The cis-SAGe occurred in disease areas, e.g. inflamed, prostatic intraepithelial neoplastic, or cancerous areas, and occasionally in normal glands.ConclusionsSTfusion detects fusion transcripts in cancer cell line and clinical tissue data, and distinguishes chimeric transcripts from chimeras caused by trans-splicing events. With STfusion and the use of C-scores, fusion transcripts can be spatially localised in clinical tissue sections on almost single cell level.
| S-EPMC7437936 | biostudies-literature