SpatialDE: identification of spatially variable genes.
Ontology highlight
ABSTRACT:
SUBMITTER: Svensson V
PROVIDER: S-EPMC6350895 | biostudies-literature | 2018 May
REPOSITORIES: biostudies-literature
Publications
SpatialDE: identification of spatially variable genes.
Svensson Valentine V Teichmann Sarah A SA Stegle Oliver O
Nature methods 20180319 5
Technological advances have made it possible to measure spatially resolved gene expression at high throughput. However, methods to analyze these data are not established. Here we describe SpatialDE, a statistical test to identify genes with spatial patterns of expression variation from multiplexed imaging or spatial RNA-sequencing data. SpatialDE also implements 'automatic expression histology', a spatial gene-clustering approach that enables expression-based tissue histology. ...[more]
Similar Datasets
Project description:Analytical solutions of partial differential equation (PDE) models describing reactive transport phenomena in saturated porous media are often used as screening tools to provide insight into contaminant fate and transport processes. While many practical modelling scenarios involve spatially variable coefficients, such as spatially variable flow velocity, v(x), or spatially variable decay rate, k(x), most analytical models deal with constant coefficients. Here we present a framework for constructing exact solutions of PDE models of reactive transport. Our approach is relevant for advection-dominant problems, and is based on a regular perturbation technique. We present a description of the solution technique for a range of one-dimensional scenarios involving constant and variable coefficients, and we show that the solutions compare well with numerical approximations. Our general approach applies to a range of initial conditions and various forms of v(x) and k(x). Instead of simply documenting specific solutions for particular cases, we present a symbolic worksheet, as supplementary material, which enables the solution to be evaluated for different choices of the initial condition, v(x) and k(x). We also discuss how the technique generalizes to apply to models of coupled multispecies reactive transport as well as higher dimensional problems.
| S-EPMC4453259 | biostudies-literature
Project description:In the studies of genomics, it is essential to select a small number of genes that are more significant than the others for research ranging from candidate gene studies to genome-wide association studies. In this study, we proposed a Bayesian method for identifying the promising candidate genes that are significantly more influential than the others. We employed the framework of variable selection and a Gibbs sampling based technique to identify significant genes. The proposed approach was applied to a genomics study for persons with chronic fatigue syndrome. Our studies show that the proposed Bayesian methodology is effective for deriving models for genomic studies and for providing information on significant genes.
| S-EPMC3169938 | biostudies-literature
Project description:Incorporating the values of the services that ecosystems provide into decision making is becoming increasingly common in nature conservation and resource management policies, both locally and globally. Yet with limited funds for conservation of threatened species and ecosystems there is a desire to identify priority areas where investment efficiently conserves multiple ecosystem services. We mapped four mangrove ecosystems services (coastal protection, fisheries, biodiversity, and carbon storage) across Fiji. Using a cost-effectiveness analysis, we prioritised mangrove areas for each service, where the effectiveness was a function of the benefits provided to the local communities, and the costs were associated with restricting specific uses of mangroves. We demonstrate that, although priority mangrove areas (top 20%) for each service can be managed at relatively low opportunity costs (ranging from 4.5 to 11.3% of overall opportunity costs), prioritising for a single service yields relatively low co-benefits due to limited geographical overlap with priority areas for other services. None-the-less, prioritisation of mangrove areas provides greater overlap of benefits than if sites were selected randomly for most ecosystem services. We discuss deficiencies in the mapping of ecosystems services in data poor regions and how this may impact upon the equity of managing mangroves for particular services across the urban-rural divide in developing countries. Finally we discuss how our maps may aid decision-makers to direct funding for mangrove management from various sources to localities that best meet funding objectives, as well as how this knowledge can aid in creating a national mangrove zoning scheme.
| S-EPMC4805192 | biostudies-other
Project description:To develop efficient algorithms for fast voxel-by-voxel quantification of tissue longitudinal relaxation time (T(1)) from variable flip angles magnetic resonance images (MRI) to reduce voxel-level noise without blurring tissue edges.T(1) estimations regularized by total variation (TV) and quadratic penalty are developed to measure T(1) from fast variable flip angles MRI and to reduce voxel-level noise without decreasing the accuracy of the estimates. First, a quadratic surrogate for a log likelihood cost function of T(1) estimation is derived based upon the majorization principle, and then the TV-regularized surrogate function is optimized by the fast iterative shrinkage thresholding algorithm. A fast optimization algorithm for the quadratically regularized T(1) estimation is also presented. The proposed methods are evaluated by the simulated and experimental MR data.The means of the T(1) values in the simulated brain data estimated by the conventional, TV-regularized, and quadratically regularized methods have less than 3% error from the true T(1) in both GM and WM tissues with image noise up to 9%. The relative standard deviations (SDs) of the T(1) values estimated by the conventional method are more than 12% and 15% when the images have 7% and 9% noise, respectively. In comparison, the TV-regularized and quadratically regularized methods are able to suppress the relative SDs of the estimated T(1) to be less than 2% and 3%, respectively, regardless of the image noise level. However, the quadratically regularized method tends to overblur the edges compared to the TV-regularized method.The spatially regularized methods improve quality of T(1) estimation from multiflip angles MRI. Quantification of dynamic contrast-enhanced MRI can benefit from the high quality measurement of native T(1).
| S-EPMC3390049 | biostudies-other
Project description:Streams are highly heterogeneous ecosystems, in terms of both geomorphology and hydrodynamics. While flow is recognized to shape the physical architecture of benthic biofilms, we do not yet understand what drives community assembly and biodiversity of benthic biofilms in the heterogeneous flow landscapes of streams. Within a metacommunity ecology framework, we experimented with streambed landscapes constructed from bedforms in large-scale flumes to illuminate the role of spatial flow heterogeneity in biofilm community composition and biodiversity in streams. Our results show that the spatial variation of hydrodynamics explained a remarkable percentage (up to 47%) of the variation in community composition along bedforms. This suggests species sorting as a model of metacommunity dynamics in stream biofilms, though natural biofilm communities will clearly not conform to a single model offered by metacommunity ecology. The spatial variation induced by the hydrodynamics along the bedforms resulted in a gradient of bacterial beta diversity, measured by a range of diversity and similarity indices, that increased with bedform height and hence with spatial flow heterogeneity at the flume level. Our results underscore the necessity to maintain small-scale physical heterogeneity for community composition and biodiversity of biofilms in stream ecosystems.
| S-EPMC2786531 | biostudies-literature
Project description:The purpose of this study is to define new genes for family risks of developing colon cancer.
| 2048176 | ecrin-mdr-crc
Project description:We consider a Bayesian hierarchical model for the integration of gene expression levels with comparative genomic hybridization (CGH) array measurements collected on the same subjects. The approach defines a measurement error model that relates the gene expression levels to latent copy number states. In turn, the latent states are related to the observed surrogate CGH measurements via a hidden Markov model. The model further incorporates variable selection with a spatial prior based on a probit link that exploits dependencies across adjacent DNA segments. Posterior inference is carried out via Markov chain Monte Carlo stochastic search techniques. We study the performance of the model in simulations and show better results than those achieved with recently proposed alternative priors. We also show an application to data from a genomic study on lung squamous cell carcinoma, where we identify potential candidates of associations between copy number variants and the transcriptional activity of target genes. Gene ontology (GO) analyses of our findings reveal enrichments in genes that code for proteins involved in cancer. Our model also identifies a number of potential candidate biomarkers for further experimental validation.
| S-EPMC4179607 | biostudies-literature
Project description:The basolateral amygdala complex (BLA), extensively connected with both local amygdalar nuclei as well as long-range circuits, is involved in a diverse array of functional roles. Understanding the mechanisms of such functional diversity will be greatly informed by understanding the cell-type-specific landscape of the BLA. Here, beginning with single-cell RNA sequencing, we identified both discrete and graded continuous gene-expression differences within the mouse BLA. Via in situ hybridization, we next mapped this discrete transcriptomic heterogeneity onto a sharp spatial border between the basal and lateral amygdala nuclei, and identified continuous spatial gene-expression gradients within each of these regions. These discrete and continuous spatial transformations of transcriptomic cell-type identity were recapitulated by local morphology as well as long-range connectivity. Thus, BLA excitatory neurons are a highly heterogenous collection of neurons that spatially covary in molecular, cellular, and circuit properties. This heterogeneity likely drives pronounced spatial variation in BLA computation and function.
| S-EPMC7486123 | biostudies-literature
Project description:Chromatin structure at the length scale encompassing local nucleosome-nucleosome interactions is thought to play a crucial role in regulating transcription and access to DNA. However, this secondary structure of chromatin remains poorly understood compared with the primary structure of single nucleosomes or the tertiary structure of long-range looping interactions. Here we report the first genome-wide map of chromatin conformation in human cells at the 1-3 nucleosome (50-500?bp) scale, obtained using ionizing radiation-induced spatially correlated cleavage of DNA with sequencing (RICC-seq) to identify DNA-DNA contacts that are spatially proximal. Unbiased analysis of RICC-seq signal reveals regional enrichment of DNA fragments characteristic of alternating rather than adjacent nucleosome interactions in tri-nucleosome units, particularly in H3K9me3-marked heterochromatin. We infer differences in the likelihood of nucleosome-nucleosome contacts among open chromatin, H3K27me3-marked, and H3K9me3-marked repressed chromatin regions. After calibrating RICC-seq signal to three-dimensional distances, we show that compact two-start helical fibre structures with stacked alternating nucleosomes are consistent with RICC-seq fragmentation patterns from H3K9me3-marked chromatin, while non-compact structures and solenoid structures are consistent with open chromatin. Our data support a model of chromatin architecture in intact interphase nuclei consistent with variable longitudinal compaction of two-start helical fibres.
| S-EPMC5526328 | biostudies-literature
Project description:The placenta is the principal organ nurturing the fetus during pregnancy and was traditionally considered to be sterile. Recent work has suggested that the placenta harbours microbial communities, however the location and possible function of these microbes remain to be confirmed and elucidated. Here, we employed genomic DNA sequencing of multiple variable (V) regions of the bacterial 16S ribosomal gene, to interrogate microbial profiles in term pregnancies, from the basal plate, which is in direct contact with maternal uterine, endothelial, and immune cells; placental villi, which are bathed in maternal blood, and fetal membranes, which encapsulate the amniotic cavity. QIIME, R package "Phyloseq" analysis was used to assess alpha and beta diversity and absolute abundance of the 16S rRNA gene per location. We demonstrate that (1) microbiota exhibit spatially distinct profiles depending on the location within the placenta and (2) "semi-composite" 16S profiles using multiple V regions validated by quantitative PCR analysis confirmed that distinct bacterial taxa dominate in different placental niches. Finally, profiles are not altered by mode of delivery. Together these findings suggest that there is niche-specificity to the placental microbiota and placental microbiome studies should consider regional differences, which may affect maternal, fetal, and/or neonatal health and physiology.
| S-EPMC5593928 | biostudies-literature