Small RNAs Establish Delays and Temporal Thresholds in Gene Expression

Non-coding RNAs are crucial regulators of gene expression in prokaryotes and eukaryotes, but it remains poorly understood how they affect the dynamics of transcriptional networks. We analyzed the temporal characteristics of the cyanobacterial iron stress response by mathematical modeling and quantitative experimental analyses, and focused on the role of a recently discovered small non-coding RNA, IsrR. We found that IsrR is responsible for a pronounced delay in the accumulation of isiA mRNA encoding the late-phase stress protein, IsiA, and that it ensures a rapid decline in isiA levels once external stress triggers are removed. These kinetic properties allow the system to selectively respond to sustained (as opposed to transient) stimuli, and thus establish a temporal threshold, which prevents energetically costly IsiA accumulation under short-term stress conditions. Biological information is frequently encoded in the quantitative aspects of intracellular signals (e.g., amplitude and duration). Our simulations reveal that competitive inhibition and regulated degradation allow intracellular regulatory networks to efficiently discriminate between transient and sustained inputs

A stem cell zoo uncovers intracellular scaling of developmental tempo across mammals

Differential speeds in biochemical reactions have been proposed to be responsible for the differences in developmental tempo between mice and humans. However, the underlying mechanism controlling the species-specific kinetics remains to be determined. Using in vitro differentiation of pluripotent stem cells, we recapitulated the segmentation clocks of diverse mammalian species varying in body weight and taxa: marmoset, rabbit, cattle, and rhinoceros. Together with mousee and human, the segmentation clock periods of the six species did not scale with the animal body weight, but with the embryogenesis length. The biochemical kinetics of the core clock gene HES7 displayed clear scaling with the species-specific segmentation clock period. However, the cellular metabolic rates did not show an evident correlation. Instead, genes involving biochemical reactions showed an expression pattern that scales with the segmentation clock period. Altogether, our stem cell zoo uncovered general scaling laws governing species-specific developmental tempo

The Kunitz-Like Modulatory Protein Haemangin Is Vital for Hard Tick Blood-Feeding Success

Ticks are serious haematophagus arthropod pests and are only second to mosquitoes as vectors of diseases of humans and animals. The salivary glands of the slower feeding hard ticks such as Haemaphysalis longicornis are a rich source of bioactive molecules and are critical to their biologic success, yet distinct molecules that help prolong parasitism on robust mammalian hosts and achieve blood-meals remain unidentified. Here, we report on the molecular and biochemical features and precise functions of a novel Kunitz inhibitor from H. longicornis salivary glands, termed Haemangin, in the modulation of angiogenesis and in persistent blood-feeding. Haemangin was shown to disrupt angiogenesis and wound healing via inhibition of vascular endothelial cell proliferation and induction of apoptosis. Further, this compound potently inactivated trypsin, chymotrypsin, and plasmin, indicating its antiproteolytic potential on angiogenic cascades. Analysis of Haemangin-specific gene expression kinetics at different blood-feeding stages of adult ticks revealed a dramatic up-regulation prior to complete feeding, which appears to be functionally linked to the acquisition of blood-meals. Notably, disruption of Haemangin-specific mRNA by a reverse genetic tool significantly diminished engorgement of adult H. longicornis, while the knock-down ticks failed to impair angiogenesis in vivo. To our knowledge, we have provided the first insights into transcriptional responses of human microvascular endothelial cells to Haemangin. DNA microarray data revealed that Haemangin altered the expression of 3,267 genes, including those of angiogenic significance, further substantiating the antiangiogenic function of Haemangin. We establish the vital roles of Haemangin in the hard tick blood-feeding process. Moreover, our results provide novel insights into the blood-feeding strategies that enable hard ticks to persistently feed and ensure full blood-meals through the modulation of angiogenesis and wound healing processes

Threshold-Free Population Analysis Identifies Larger DRG Neurons to Respond Stronger to NGF Stimulation

Sensory neurons in dorsal root ganglia (DRG) are highly heterogeneous in terms of cell size, protein expression, and signaling activity. To analyze their heterogeneity, threshold-based methods are commonly used, which often yield highly variable results due to the subjectivity of the individual investigator. In this work, we introduce a threshold-free analysis approach for sparse and highly heterogeneous datasets obtained from cultures of sensory neurons. This approach is based on population estimates and completely free of investigator-set parameters. With a quantitative automated microscope we measured the signaling state of single DRG neurons by immunofluorescently labeling phosphorylated, i.e., activated Erk1/2. The population density of sensory neurons with and without pain-sensitizing nerve growth factor (NGF) treatment was estimated using a kernel density estimator (KDE). By subtraction of both densities and integration of the positive part, a robust estimate for the size of the responsive subpopulations was obtained. To assure sufficiently large datasets, we determined the number of cells required for reliable estimates using a bootstrapping approach. The proposed methods were employed to analyze response kinetics and response amplitude of DRG neurons after NGF stimulation. We thereby determined the portion of NGF responsive cells on a true population basis. The analysis of the dose dependent NGF response unraveled a biphasic behavior, while the study of its time dependence showed a rapid response, which approached a steady state after less than five minutes. Analyzing two parameter correlations, we found that not only the number of responsive small-sized neurons exceeds the number of responsive large-sized neurons—which is commonly reported and could be explained by the excess of small-sized cells—but also the probability that small-sized cells respond to NGF is higher. In contrast, medium-sized and large-sized neurons showed a larger response amplitude in their mean Erk1/2 activity

Impaired mitochondrial biogenesis contributes to depletion of functional mitochondria in chronic MPP+ toxicity: dual roles for ERK1/2

The regulation of mitochondrial quality has emerged as a central issue in neurodegeneration, diabetes, and cancer. We utilized repeated low-dose applications of the complex I inhibitor 1-methyl-4-phenylpyridinium (MPP+) over 2 weeks to study cellular responses to chronic mitochondrial stress. Chronic MPP+ triggered depletion of functional mitochondria resulting in diminished capacities for aerobic respiration. Inhibiting autophagy/mitophagy only partially restored mitochondrial content. In contrast, inhibiting activation of extracellular signal-regulated protein kinases conferred complete cytoprotection with full restoration of mitochondrial functional and morphological parameters, enhancing spare respiratory capacity in MPP+ co-treated cells above that of control cells. Reversal of mitochondrial injury occurred when U0126 was added 1 week after MPP+, implicating enhanced repair mechanisms. Chronic MPP+ caused a >90% decrease in complex I subunits, along with decreases in complex III and IV subunits. Decreases in respiratory complex subunits were reversed by co-treatment with U0126, ERK1/2 RNAi or transfection of dominant-negative MEK1, but only partially restored by degradation inhibitors. Chronic MPP+ also suppressed the de novo synthesis of mitochondrial DNA-encoded proteins, accompanied by decreased expression of the mitochondrial transcription factor TFAM. U0126 completely reversed each of these deficits in mitochondrial translation and protein expression. These data indicate a key, limiting role for mitochondrial biogenesis in determining the outcome of injuries associated with elevated mitophagy

A comparison of approximation techniques for variance-based sensitivity analysis of biochemical reaction systems

<p>Abstract</p> <p>Background</p> <p>Sensitivity analysis is an indispensable tool for the analysis of complex systems. In a recent paper, we have introduced a thermodynamically consistent variance-based sensitivity analysis approach for studying the robustness and fragility properties of biochemical reaction systems under uncertainty in the standard chemical potentials of the activated complexes of the reactions and the standard chemical potentials of the molecular species. In that approach, key sensitivity indices were estimated by Monte Carlo sampling, which is computationally very demanding and impractical for large biochemical reaction systems. Computationally efficient algorithms are needed to make variance-based sensitivity analysis applicable to realistic cellular networks, modeled by biochemical reaction systems that consist of a large number of reactions and molecular species.</p> <p>Results</p> <p>We present four techniques, derivative approximation (DA), polynomial approximation (PA), Gauss-Hermite integration (GHI), and orthonormal Hermite approximation (OHA), for <it>analytically </it>approximating the variance-based sensitivity indices associated with a biochemical reaction system. By using a well-known model of the mitogen-activated protein kinase signaling cascade as a case study, we numerically compare the approximation quality of these techniques against traditional Monte Carlo sampling. Our results indicate that, although DA is computationally the most attractive technique, special care should be exercised when using it for sensitivity analysis, since it may only be accurate at low levels of uncertainty. On the other hand, PA, GHI, and OHA are computationally more demanding than DA but can work well at high levels of uncertainty. GHI results in a slightly better accuracy than PA, but it is more difficult to implement. OHA produces the most accurate approximation results and can be implemented in a straightforward manner. It turns out that the computational cost of the four approximation techniques considered in this paper is orders of magnitude smaller than traditional Monte Carlo estimation. Software, coded in MATLAB^®, which implements all sensitivity analysis techniques discussed in this paper, is available free of charge.</p> <p>Conclusions</p> <p>Estimating variance-based sensitivity indices of a large biochemical reaction system is a computationally challenging task that can only be addressed via approximations. Among the methods presented in this paper, a technique based on orthonormal Hermite polynomials seems to be an acceptable candidate for the job, producing very good approximation results for a wide range of uncertainty levels in a fraction of the time required by traditional Monte Carlo sampling.</p

Co-regulation map of the human proteome enables identification of protein functions

This is the author accepted manuscript. The final version is available from Nature Research via the DOI in this recordData availability: All mass spectrometry raw files generated in-house have been deposited in the ProteomeXchange Consortium (http://proteomecentral.proteomexchange.org) via the PRIDE partner repository36 with the dataset identifier PXD008888. The co-regulation map is hosted on our website at www.proteomeHD.net, and pair-wise co-regulation scores are available through STRING (https://string-db.org). A network of the top 0.5% co-regulated protein pairs can be explored interactively on NDEx (https://doi.org/10.18119/N9N30Q).Code availability: Data analysis was performed in R 3.5.1. R scripts and input files required to reproduce the results of this manuscript are available in the following GitHub repository: https://github.com/Rappsilber-Laboratory/ProteomeHD. R scripts related specifically to the benchmarking of the treeClust algorithm using synthetic data are available in the following GitHub repository: https://github.com/Rappsilber-Laboratory/treeClust-benchmarking. The R package data.table was used for fast data processing. Figures were prepared using ggplot2, gridExtra, cowplot and viridis.Note that the title of the AAM is different from the published versionThe annotation of protein function is a longstanding challenge of cell biology that suffers from the sheer magnitude of the task. Here we present ProteomeHD, which documents the response of 10,323 human proteins to 294 biological perturbations, measured by isotope-labelling mass spectrometry. We reveal functional associations between human proteins using the treeClust machine learning algorithm, which we show to improve protein co-regulation analysis due to robust selectivity for close linear relationships. Our co-regulation map identifies a functional context for many uncharacterized proteins, including microproteins that are difficult to study with traditional methods. Co-regulation also captures relationships between proteins which do not physically interact or co-localize. For example, co-regulation of the peroxisomal membrane protein PEX11β with mitochondrial respiration factors led us to discover a novel organelle interface between peroxisomes and mitochondria in mammalian cells. The co-regulation map can be explored at www.proteomeHD.net .Biotechnology & Biological Sciences Research Council (BBSRC)European Commissio

Genomic and Transcriptional Co-Localization of Protein-Coding and Long Non-Coding RNA Pairs in the Developing Brain

Besides protein-coding mRNAs, eukaryotic transcriptomes include many long non-protein-coding RNAs (ncRNAs) of unknown function that are transcribed away from protein-coding loci. Here, we have identified 659 intergenic long ncRNAs whose genomic sequences individually exhibit evolutionary constraint, a hallmark of functionality. Of this set, those expressed in the brain are more frequently conserved and are significantly enriched with predicted RNA secondary structures. Furthermore, brain-expressed long ncRNAs are preferentially located adjacent to protein-coding genes that are (1) also expressed in the brain and (2) involved in transcriptional regulation or in nervous system development. This led us to the hypothesis that spatiotemporal co-expression of ncRNAs and nearby protein-coding genes represents a general phenomenon, a prediction that was confirmed subsequently by in situ hybridisation in developing and adult mouse brain. We provide the full set of constrained long ncRNAs as an important experimental resource and present, for the first time, substantive and predictive criteria for prioritising long ncRNA and mRNA transcript pairs when investigating their biological functions and contributions to development and disease

Stochasticity in Protein Levels Drives Colinearity of Gene Order in Metabolic Operons of Escherichia coli

Gene order in some bacterial metabolic operons reflects ordering in the metabolic pathway. That this is true uniquely for operons expressed at low levels highlights the selective importance of fluctuations in protein levels