Dataset Information

Broad distribution spectrum from Gaussian to power law appears in stochastic variations in RNA-seq data.

ABSTRACT: Gene expression levels exhibit stochastic variations among genetically identical organisms under the same environmental conditions. In many recent transcriptome analyses based on RNA sequencing (RNA-seq), variations in gene expression levels among replicates were assumed to follow a negative binomial distribution, although the physiological basis of this assumption remains unclear. In this study, RNA-seq data were obtained from Arabidopsis thaliana under eight conditions (21-27 replicates), and the characteristics of gene-dependent empirical probability density function (ePDF) profiles of gene expression levels were analyzed. For A. thaliana and Saccharomyces cerevisiae, various types of ePDF of gene expression levels were obtained that were classified as Gaussian, power law-like containing a long tail, or intermediate. These ePDF profiles were well fitted with a Gauss-power mixing distribution function derived from a simple model of a stochastic transcriptional network containing a feedback loop. The fitting function suggested that gene expression levels with long-tailed ePDFs would be strongly influenced by feedback regulation. Furthermore, the features of gene expression levels are correlated with their functions, with the levels of essential genes tending to follow a Gaussian-like ePDF while those of genes encoding nucleic acid-binding proteins and transcription factors exhibit long-tailed ePDF.

SUBMITTER: Awazu A

PROVIDER: S-EPMC5974282 | biostudies-literature | 2018 May

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Broad distribution spectrum from Gaussian to power law appears in stochastic variations in RNA-seq data.

Awazu Akinori A Tanabe Takahiro T Kamitani Mari M Tezuka Ayumi A Nagano Atsushi J AJ

Scientific reports 20180529 1

Gene expression levels exhibit stochastic variations among genetically identical organisms under the same environmental conditions. In many recent transcriptome analyses based on RNA sequencing (RNA-seq), variations in gene expression levels among replicates were assumed to follow a negative binomial distribution, although the physiological basis of this assumption remains unclear. In this study, RNA-seq data were obtained from Arabidopsis thaliana under eight conditions (21-27 replicates), and ...[more]

PMID: 29844539

Similar Datasets

Project description:To determine the x-ray photon energy dependence of the anatomic power spectrum of the breast when imaged with dedicated breast computed tomography (CT).A theoretical framework for scaling the empirically determined anatomic power spectrum at one x-ray photon energy to that at any given x-ray photon energy when imaged with dedicated breast CT was developed. Theory predicted that when the anatomic power spectrum is fitted with a power curve of the form k f(-β), where k and β are fit coefficients and f is spatial frequency, the exponent β would be independent of x-ray photon energy (E), and the amplitude k scales with the square of the difference in energy-dependent linear attenuation coefficients of fibroglandular and adipose tissues. Twenty mastectomy specimens based numerical phantoms that were previously imaged with a benchtop flat-panel cone-beam CT system were converted to 3D distribution of glandular weight fraction (f(g)) and were used to verify the theoretical findings. The 3D power spectrum was computed in terms of f(g) and after converting to linear attenuation coefficients at monoenergetic x-ray photon energies of 20-80 keV in 5 keV intervals. The 1D power spectra along the axes were extracted and fitted with a power curve of the form k f(-β). The energy dependence of k and β were analyzed.For the 20 mastectomy specimen based numerical phantoms used in the study, the exponent β was found to be in the range of 2.34-2.42, depending on the axis of measurement. Numerical simulations agreed with the theoretical predictions that for a power-law anatomic spectrum of the form k f(-β), β was independent of E and k(E) = k(1)[μ(g)(E) - μ(a)(E)](2), where k(1) is a constant, and μ(g)(E) and μ(a)(E) represent the energy-dependent linear attenuation coefficients of fibroglandular and adipose tissues, respectively.Numerical simulations confirmed the theoretical predictions that in dedicated breast CT, the spatial frequency dependence of the anatomic power spectrum will be independent of x-ray photon energy, and the amplitude of the anatomic power spectrum scales by the square of difference in linear attenuation coefficients of fibroglandular and adipose tissues.

Project description:Viscoelastic characteristics of many materials falling under the category of soft glassy substances, including biological tissue, often exhibit a mechanical complex modulus Y(?) well described by a fractional derivative model: Y(?) = E(i?/?)k, where E = a generalized viscoelastic stiffness; i = (-1)1/2; ? = angular frequency; ? = scaling factor; and k = an exponent valued between 0 and 1. The term "fractional derivative" refers to the value of k: when k = 0 the viscoelastic response is purely elastic, and when k = 1 the response is purely viscous. We provide an analytical derivation of the fractional derivative complex modulus based on the hypothesis that the viscoelastic response arises from many intermittent molecular crosslinks, whose lifetimes longer than a critical threshold lifetime, tcrit, are distributed with an inverse power law proportional to t-(k+2). We demonstrate that E is proportional to the number and stiffness of crosslinks formed at any moment; the scaling factor ? is equivalent to reciprocal of tcrit; and the relative mean lifetime of the attached crosslinks is inversely proportional to the parameter k. To test whether electrostatic molecular bonds could be responsible for the fractional derivative viscoelasticity, we used chemically skinned human skeletal muscle as a one-dimensional model of a soft glassy substance. A reduction in ionic strength from 175 to 110 mEq resulted in a larger E with no change in k, consistent with a higher probability of interfilament molecular interactions. Thick to thin filament spacing was reduced by applying 4% w/v of the osmolyte Dextran T500, which also resulted in a larger E, indicating a greater probability of crosslink formation in proportion to proximity. A 10°C increase in temperature resulted in an increase in k, which corresponded to a decrease in cross-bridge attachment lifetime expected with higher temperatures. These theoretical and experimental results suggest that the fractional derivative viscoelasticity observed in some biological tissue arises as a mechanical consequence of electrostatic interactions, whose longest lifetimes are distributed with an inverse power law.

Dataset Information

Broad distribution spectrum from Gaussian to power law appears in stochastic variations in RNA-seq data.

Publications

Broad distribution spectrum from Gaussian to power law appears in stochastic variations in RNA-seq data.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets