Genome-wide modeling of transcription kinetics reveals patterns of RNA production delays

Kinetic pathway of HIV-1 TAR cotranscriptional folding

The Trans-Activator Receptor (TAR) RNA, located at the 5'-end untranslated region (5' UTR) of the human immunodeficiency virus type 1 (HIV-1), is pivotal in the virus's life cycle. As the initial functional domain, it folds during the transcription of viral mRNA. Although TAR's role in recruiting the Tat protein for trans-activation is established, the detailed kinetic mechanisms at play during early transcription, especially at points of temporary transcriptional pausing, remain elusive. Moreover, the precise physical processes of transcriptional pause and subsequent escape are not fully elucidated. This study focuses on the folding kinetics of TAR and the biological implications by integrating computer simulations of RNA folding during transcription with nuclear magnetic resonance (NMR) spectroscopy data. The findings reveal insights into the folding mechanism of a non-native intermediate that triggers transcriptional pause, along with different folding pathways leading to transcriptional pause and readthrough. The profiling of the cotranscriptional folding pathway and identification of kinetic structural intermediates reveal a novel mechanism for viral transcriptional regulation, which could pave the way for new antiviral drug designs targeting kinetic cotranscriptional folding pathways in viral RNAs.

Deconstructing cancer with precision genome editing

Recent advances in genome editing technologies are allowing investigators to engineer and study cancer-associated mutations in their endogenous genetic contexts with high precision and efficiency. Of these, base editing and prime editing are quickly becoming gold-standards in the field due to their versatility and scalability. Here, we review the merits and limitations of these precision genome editing technologies, their application to modern cancer research, and speculate how these could be integrated to address future directions in the field.

Post-transcriptional splicing can occur in a slow-moving zone around the gene

Splicing is the stepwise molecular process by which introns are removed from pre-mRNA and exons are joined together to form mature mRNA sequences. The ordering and spatial distribution of these steps remain controversial, with opposing models suggesting splicing occurs either during or after transcription. We used single-molecule RNA FISH, expansion microscopy, and live-cell imaging to reveal the spatiotemporal distribution of nascent transcripts in mammalian cells. At super-resolution levels, we found that pre-mRNA formed clouds around the transcription site. These clouds indicate the existence of a transcription-site-proximal zone through which RNA move more slowly than in the nucleoplasm. Full-length pre-mRNA undergo continuous splicing as they move through this zone following transcription, suggesting a model in which splicing can occur post-transcriptionally but still within the proximity of the transcription site, thus seeming co-transcriptional by most assays. These results may unify conflicting reports of co-transcriptional versus post-transcriptional splicing.

Enhancer-driven 3D chromatin domain folding modulates transcription in human mammary tumor cells

The genome is organized in functional compartments and structural domains at the sub-megabase scale. How within these domains interactions between numerous cis-acting enhancers and promoters regulate transcription remains an open question. Here, we determined chromatin folding and composition over several hundred kb around estrogen-responsive genes in human breast cancer cell lines after hormone stimulation. Modeling of 5C data at 1.8 kb resolution was combined with quantitative 3D analysis of multicolor FISH measurements at 100 nm resolution and integrated with ChIP-seq data on transcription factor binding and histone modifications. We found that rapid estradiol induction of the progesterone gene expression occurs in the context of preexisting, cell type-specific chromosomal architectures encompassing the 90 kb progesterone gene coding region and an enhancer-spiked 5' 300 kb upstream genomic region. In response to estradiol, interactions between estrogen receptor α (ERα) bound regulatory elements are reinforced. Whereas initial enhancer-gene contacts coincide with RNA Pol 2 binding and transcription initiation, sustained hormone stimulation promotes ERα accumulation creating a regulatory hub stimulating transcript synthesis. In addition to implications for estrogen receptor signaling, we uncover that preestablished chromatin architectures efficiently regulate gene expression upon stimulation without the need for de novo extensive rewiring of long-range chromatin interactions.

Circadian ribosome profiling reveals a role for the Period2 upstream open reading frame in sleep

Many mammalian proteins have circadian cycles of production and degradation, and many of these rhythms are altered posttranscriptionally. We used ribosome profiling to examine posttranscriptional control of circadian rhythms by quantifying RNA translation in the liver over a 24-h period from circadian-entrained mice transferred to constant darkness conditions and by comparing ribosome binding levels to protein levels for 16 circadian proteins. We observed large differences in ribosome binding levels compared to protein levels, and we observed delays between peak ribosome binding and peak protein abundance. We found extensive binding of ribosomes to upstream open reading frames (uORFs) in circadian mRNAs, including the core clock gene Period2 (Per2). An increase in the number of uORFs in the 5'UTR was associated with a decrease in ribosome binding in the main coding sequence and a reduction in expression of synthetic reporter constructs. Mutation of the Per2 uORF increased luciferase and fluorescence reporter expression in 3T3 cells and increased luciferase expression in PER2:LUC MEF cells. Mutation of the Per2 uORF in mice increased Per2 mRNA expression, enhanced ribosome binding on Per2, and reduced total sleep time compared to that in wild-type mice. These results suggest that uORFs affect mRNA posttranscriptionally, which can impact physiological rhythms and sleep.

Meta-analysis of integrated ChIP-seq and transcriptome data revealed genomic regions affected by estrogen receptor alpha in breast cancer

BACKGROUND

The largest group of patients with breast cancer are estrogen receptor-positive (ER+) type. The estrogen receptor acts as a transcription factor and triggers cell proliferation and differentiation. Hence, investigating ER-DNA interaction genomic regions can help identify genes directly regulated by ER and understand the mechanism of ER action in cancer progression.

METHODS

In the present study, we employed a workflow to do a meta-analysis of ChIP-seq data of ER+ cell lines stimulated with 10 nM and 100 nM of E2. All publicly available data sets were re-analyzed with the same platform. Then, the known and unknown batch effects were removed. Finally, the meta-analysis was performed to obtain meta-differentially bound sites in estrogen-treated MCF7 cell lines compared to vehicles (as control). Also, the meta-analysis results were compared with the results of T47D cell lines for more precision. Enrichment analyses were also employed to find the functional importance of common meta-differentially bound sites and associated genes among both cell lines.

RESULTS

Remarkably, POU5F1B, ZNF662, ZNF442, KIN, ZNF410, and SGSM2 transcription factors were recognized in the meta-analysis but not in individual studies. Enrichment of the meta-differentially bound sites resulted in the candidacy of pathways not previously reported in breast cancer. PCGF2, HNF1B, and ZBED6 transcription factors were also predicted through the enrichment analysis of associated genes. In addition, comparing the meta-analysis results of both ChIP-seq and RNA-seq data showed that many transcription factors affected by ER were up-regulated.

CONCLUSION

The meta-analysis of ChIP-seq data of estrogen-treated MCF7 cell line leads to the identification of new binding sites of ER that have not been previously reported. Also, enrichment of the meta-differentially bound sites and their associated genes revealed new terms and pathways involved in the development of breast cancer which should be examined in future in vitro and in vivo studies.

Bayesian inference of polymerase dynamics over the exclusion process

Transcription is a complex phenomenon that permits the conversion of genetic information into phenotype by means of an enzyme called RNA polymerase, which erratically moves along and scans the DNA template. We perform Bayesian inference over a paradigmatic mechanistic model of non-equilibrium statistical physics, i.e. the asymmetric exclusion processes in the hydrodynamic limit, assuming a Gaussian process prior for the polymerase progression rate as a latent variable. Our framework allows us to infer the speed of polymerases during transcription given their spatial distribution, while avoiding the explicit inversion of the system's dynamics. The results, which show processing rates strongly varying with genomic position and minor role of traffic-like congestion, may have strong implications for the understanding of gene expression.

Combined modelling of mRNA decay dynamics and single-molecule imaging in the Drosophila embryo uncovers a role for P-bodies in 5' to 3' degradation

Regulation of mRNA degradation is critical for a diverse array of cellular processes and developmental cell fate decisions. Many methods for determining mRNA half-lives rely on transcriptional inhibition or metabolic labelling. Here, we use a non-invasive method for estimating half-lives for hundreds of mRNAs in the early Drosophila embryo. This approach uses the intronic and exonic reads from a total RNA-seq time series and Gaussian process regression to model the dynamics of premature and mature mRNAs. We show how regulation of mRNA stability is used to establish a range of mature mRNA dynamics during embryogenesis, despite shared transcription profiles. Using single-molecule imaging, we provide evidence that, for the mRNAs tested, there is a correlation between short half-life and mRNA association with P-bodies. Moreover, we detect an enrichment of mRNA 3' ends in P-bodies in the early embryo, consistent with 5' to 3' degradation occurring in P-bodies for at least a subset of mRNAs. We discuss our findings in relation to recently published data suggesting that the primary function of P-bodies in other biological contexts is mRNA storage.

Constrained Fourier estimation of short-term time-series gene expression data reduces noise and improves clustering and gene regulatory network predictions

BACKGROUND

Biological data suffers from noise that is inherent in the measurements. This is particularly true for time-series gene expression measurements. Nevertheless, in order to to explore cellular dynamics, scientists employ such noisy measurements in predictive and clustering tools. However, noisy data can not only obscure the genes temporal patterns, but applying predictive and clustering tools on noisy data may yield inconsistent, and potentially incorrect, results.

RESULTS

To reduce the noise of short-term (< 48 h) time-series expression data, we relied on the three basic temporal patterns of gene expression: waves, impulses and sustained responses. We constrained the estimation of the true signals to these patterns by estimating the parameters of first and second-order Fourier functions and using the nonlinear least-squares trust-region optimization technique. Our approach lowered the noise in at least 85% of synthetic time-series expression data, significantly more than the spline method ([Formula: see text]). When the data contained a higher signal-to-noise ratio, our method allowed downstream network component analyses to calculate consistent and accurate predictions, particularly when the noise variance was high. Conversely, these tools led to erroneous results from untreated noisy data. Our results suggest that at least 5-7 time points are required to efficiently de-noise logarithmic scaled time-series expression data. Investing in sampling additional time points provides little benefit to clustering and prediction accuracy.

CONCLUSIONS

Our constrained Fourier de-noising method helps to cluster noisy gene expression and interpret dynamic gene networks more accurately. The benefit of noise reduction is large and can constitute the difference between a successful application and a failing one.

Landscape Zooming toward the Prediction of RNA Cotranscriptional Folding

RNA molecules fold as they are transcribed. Cotranscriptional folding of RNA plays a critical role in RNA functions in vivo. Present computational strategies focus on simulations where large structural changes may not be completely sampled. Here, we describe an alternative approach to predicting cotranscriptional RNA folding by zooming in and out of the RNA folding energy landscape. By classifying the RNA structural ensemble into "partitions" based on long, stable helices, we zoom out of the landscape and predict the overall slow folding kinetics from the interpartition kinetic network, and for each interpartition transition, we zoom in on the landscape to simulate the kinetics. Applications of the model to the 117-nucleotide E. coli SRP RNA and the 59-nucleotide HIV-1 TAR RNA show agreements with the experimental data and new structural and kinetic insights into biologically significant conformational switches and pathways for these important systems. This approach, by zooming in/out of an RNA folding landscape at different resolutions, might allow us to treat large RNAs in vivo with transcriptional pause, transcription speed, and other in vivo effects.

Translation complex stabilization on messenger RNA and footprint profiling to study the RNA responses and dynamics of protein biosynthesis in the cells

During protein biosynthesis, ribosomes bind to messenger (m)RNA, locate its protein-coding information, and translate the nucleotide triplets sequentially as codons into the corresponding sequence of amino acids, forming proteins. Non-coding mRNA features, such as 5' and 3' untranslated regions (UTRs), start sites or stop codons of different efficiency, stretches of slower or faster code and nascent polypeptide interactions can alter the translation rates transcript-wise. Most of the homeostatic and signal response pathways of the cells converge on individual mRNA control, as well as alter the global translation output. Among the multitude of approaches to study translational control, one of the most powerful is to infer the locations of translational complexes on mRNA based on the mRNA fragments protected by these complexes from endonucleolytic hydrolysis, or footprints. Translation complex profiling by high-throughput sequencing of the footprints allows to quantify the transcript-wise, as well as global, alterations of translation, and uncover the underlying control mechanisms by attributing footprint locations and sizes to different configurations of the translational complexes. The accuracy of all footprint profiling approaches critically depends on the fidelity of footprint generation and many methods have emerged to preserve certain or multiple configurations of the translational complexes, often in challenging biological material. In this review, a systematic summary of approaches to stabilize translational complexes on mRNA for footprinting is presented and major findings are discussed. Future directions of translation footprint profiling are outlined, focusing on the fidelity and accuracy of inference of the native in vivo translation complex distribution on mRNA.

Inferring kinetic parameters of oscillatory gene regulation from single cell time-series data

Gene expression dynamics, such as stochastic oscillations and aperiodic fluctuations, have been associated with cell fate changes in multiple contexts, including development and cancer. Single cell live imaging of protein expression with endogenous reporters is widely used to observe such gene expression dynamics. However, the experimental investigation of regulatory mechanisms underlying the observed dynamics is challenging, since these mechanisms include complex interactions of multiple processes, including transcription, translation and protein degradation. Here, we present a Bayesian method to infer kinetic parameters of oscillatory gene expression regulation using an auto-negative feedback motif with delay. Specifically, we use a delay-adapted nonlinear Kalman filter within a Metropolis-adjusted Langevin algorithm to identify posterior probability distributions. Our method can be applied to time-series data on gene expression from single cells and is able to infer multiple parameters simultaneously. We apply it to published data on murine neural progenitor cells and show that it outperforms alternative methods. We further analyse how parameter uncertainty depends on the duration and time resolution of an imaging experiment, to make experimental design recommendations. This work demonstrates the utility of parameter inference on time course data from single cells and enables new studies on cell fate changes and population heterogeneity.

Enhancer RNAs Mediate Estrogen-Induced Decommissioning of Selective Enhancers by Recruiting ERα and Its Cofactor

The function of enhancer RNAs (eRNAs) in transcriptional regulation remains obscure. By analyzing the genome-wide nascent transcript profiles in breast cancer cells, we identify a special group of eRNAs that are essential for estrogen-induced transcriptional repression. Using eRNAs of TM4SF1 and EFEMP1 as the paradigms, we find that these RNA molecules not only stabilize promoter-enhancer interactions but also recruit liganded estrogen receptor α (ERα) to particular enhancer regions, facilitate the formation of a functional transcriptional complex, and cause gene silencing. Interestingly, ERα is shown to directly bind with eRNAs by its DNA-binding domain. These eRNAs help with the formation of a specific ERα-centered transcriptional complex and promote the association of the histone demethylase KDM2A, which dismisses RNA polymerase II from designated enhancers and suppresses the transcription of target genes. Our work demonstrates a complete mechanism underlying the action of eRNAs in modulating and refining the locus-specific transcriptional program.

Yang et al. identified a group of eRNAs that are essential for estrogen-induced transcriptional repression by assisting with the chromatin recruitment of ERα through binding to its DNA-binding domain and facilitating the interaction of ERα with its cofactors, which leads to the dismissal of RNA polymerase II.

Physically-Inspired Gaussian Process Models for Post-Transcriptional Regulation in Drosophila

The regulatory process of Drosophila is thoroughly studied for understanding a great variety of biological principles. While pattern-forming gene networks are analyzed in the transcription step, post-transcriptional events (e.g., translation, protein processing) play an important role in establishing protein expression patterns and levels. Since the post-transcriptional regulation of Drosophila depends on spatiotemporal interactions between mRNAs and gap proteins, proper physically-inspired stochastic models are required to study the link between both quantities. Previous research attempts have shown that using Gaussian processes (GPs) and differential equations lead to promising predictions when analyzing regulatory networks. Here, we aim at further investigating two types of physically-inspired GP models based on a reaction-diffusion equation where the main difference lies in where the prior is placed. While one of them has been studied previously using protein data only, the other is novel and yields a simple approach requiring only the differentiation of kernel functions. In contrast to other stochastic frameworks, discretizing the spatial space is not required here. Both GP models are tested under different conditions depending on the availability of gap gene mRNA expression data. Finally, their performances are assessed on a high-resolution dataset describing the blastoderm stage of the early embryo of Drosophila melanogaster.

Kinetic modeling of stem cell transcriptome dynamics to identify regulatory modules of normal and disturbed neuroectodermal differentiation

Thousands of transcriptome data sets are available, but approaches for their use in dynamic cell response modelling are few, especially for processes affected simultaneously by two orthogonal influencing variables. We approached this problem for neuroepithelial development of human pluripotent stem cells (differentiation variable), in the presence or absence of valproic acid (signaling variable). Using few basic assumptions (sequential differentiation states of cells; discrete on/off states for individual genes in these states), and time-resolved transcriptome data, a comprehensive model of spontaneous and perturbed gene expression dynamics was developed. The model made reliable predictions (average correlation of 0.85 between predicted and subsequently tested expression values). Even regulations predicted to be non-monotonic were successfully validated by PCR in new sets of experiments. Transient patterns of gene regulation were identified from model predictions. They pointed towards activation of Wnt signaling as a candidate pathway leading to a redirection of differentiation away from neuroepithelial cells towards neural crest. Intervention experiments, using a Wnt/beta-catenin antagonist, led to a phenotypic rescue of this disturbed differentiation. Thus, our broadly applicable model allows the analysis of transcriptome changes in complex time/perturbation matrices.

The common misuse of noise decomposition as applied to genetic systems

The noise-decomposition technique is applied in several fields, including genetic systems, optical images, recording, and navigation. In genetic systems, noise decomposition is usually achieved by using two reporters [Elowitz M.B., Levine A.J., Siggia E.D., Swain P·S., 2002. Stochastic gene expression in a single cell. Science 297, 1183-6.]. A reporter is a protein with fluorescence, an RNA hybridized with a fluorescent probe, or any other detectable intracellular component. If a reporter is constructed in addition to the original reporter, the system's stochasticity may change. Such phenomena became severe for genes in plasmids with a high copy number. By SSA (stochastic simulation algorithm), we observed an approximately 50% increment in the coefficient of variation while introducing additional reporters. Besides, if two reporters respond to the upstream element at a different time, the trunk noise (or extrinsic noise) cannot be accurately determined. This is because the "calculative trunk noise" changes along with the delay, though the real trunk noise does not. For RNA reporters, a 5-min transcriptional delay caused a calculative trunk noise that was 90% less than the real trunk noise. Fortunately, this problem is negligible when the degradation rate constant is low, and it is usually true in the case of the protein reporters. One can check the lifespan of the reporter before applying the noise-decomposition technique.

Long-Chain Non-Coding SOX21-AS1 Promotes Proliferation and Migration of Breast Cancer Cells Through the PI3K/AKT Signaling Pathway

Aim

This study aimed to investigate the effect of long-chain non-coding SOX21-AS1 on the proliferation and migration of breast cancer (BC) cells through the PI3K/AKT signaling pathway.

Methods

Eighty-eight BC and adjacent tissues were collected, and BC cells and normal breast epithelial cells were purchased. LncRNA SOX21-AS1 expression in tissues and cells was detected by RT-PCR. miR-NC, si-SOX21-AS1, and Sh-SOX21-AS1 were transfected into BC cells. The PI3K/AKT signaling pathway was interfered with L740Y-P (activator of the PI3K/AKT pathway) and LY294002 (inhibitor of the PI3K/AKT pathway) in BC cells. The SOX21-AS1 expression in BC tissues and cells was tested by qRT-PCR, and the expression levels of p-PI3K, p-AKT, N-cadherin, E-cadherin, and vimentin were detected by WB.

Results

SOX21-AS1 was highly expressed in BC, and the p-PI3K and p-AKT levels were also high. Cell experiments showed that inhibiting SOX21-AS1 expression could inhibit the proliferation, invasion, migration, and EMT of BC cells, and up-regulating its expression could promote the proliferation, invasion, migration, and EMT of ovarian cancer cells. The tumor-forming experiment in nude mice was consistent with the results in vitro. 740Y-P intervention could reverse the inhibition effect of Si-SOX21-AS1 on BC cell proliferation, invasion, migration, and EMT, while LY294002 intervention could reverse the promotion effect of Sh-SOX21-AS1 on BC cell proliferation, invasion, migration, and EMT.

Conclusion

SOX21-AS1 is highly expressed in BC tissues. Silencing BC expression can inhibit the proliferation, invasion, migration, and EMT of cells by inhibiting the PI3K/AKT signaling pathway, which may be a new target for diagnosis and treatment.

Combinatorial dynamics of protein synthesis time delay and negative feedback loop in NF-κB signalling pathway

The transcription factor NF-κB links immune response and inflammatory reaction and its different oscillation patterns determine different cell fates. In this study, a mathematical model with IκBα protein synthesis time delay is developed based on the experimental evidences. The results show that time delay has the ability to drive oscillation of NF-κB via Hopf bifurcation. Meanwhile, the amplitude and period are sensitive to the time delay. Moreover, the time delay threshold is a function of four parameters characterising the negative feedback loop. Likewise, the parameters also have effects on the amplitude and period of NF-κB oscillation induced by time delay. Therefore, the oscillation patterns of NF-κB are collaborative results of time delay coupled with the negative feedback loop. These results not only enhance the understanding of NF-κB biological oscillation but also provide clues for the development of anti-inflammatory or anti-cancer drugs.

Cell type specific gene expression profiling reveals a role for complement component C3 in neutrophil responses to tissue damage

Tissue damage induces rapid recruitment of leukocytes and changes in the transcriptional landscape that influence wound healing. However, the cell-type specific transcriptional changes that influence leukocyte function and tissue repair have not been well characterized. Here, we employed translating ribosome affinity purification (TRAP) and RNA sequencing, TRAP-seq, in larval zebrafish to identify genes differentially expressed in neutrophils, macrophages, and epithelial cells in response to wounding. We identified the complement pathway and c3a.1, homologous to the C3 component of human complement, as significantly increased in neutrophils in response to wounds. c3a.1^−/− zebrafish larvae have impaired neutrophil directed migration to tail wounds with an initial lag in recruitment early after wounding. Moreover, c3a.1^−/− zebrafish larvae have impaired recruitment to localized bacterial infections and reduced survival that is, at least in part, neutrophil mediated. Together, our findings support the power of TRAP-seq to identify cell type specific changes in gene expression that influence neutrophil behavior in response to tissue damage.

Oscillatory Dynamics of p53-Mdm2 Circuit in Response to DNA Damage Caused by Ionizing Radiation

Although the dynamical behavior of the p53-Mdm2 loop has been extensively studied, the understanding of the mechanism underlying the regulation of this pathway still remains limited. Herein, we developed an integrated model with five basic components and three ubiquitous time delays for the p53-Mdm2 interaction in response to DNA damage following ionizing radiation (IR). We showed that a sufficient amount of activated ATM level can initiate the p53 oscillations with nearly the same amplitude over a wide range of the ATM level; a proper range of p53 level is also required for generating the oscillations, for too high or too low levels it would fail to generate the oscillations; and increased Mdm2 level leads to decreased amplitude of the p53 oscillation and reduced expression of the p53 activity. Moreover, we found that the negative feedback loop formed between p53 and nuclear Mdm2 plays a dominant role in determining the p53 dynamics, whereas when interaction strength of the negative feedback loop becomes weaker, the positive feedback loop formed between p53 and cytoplasmatic Mdm2 can induce different types of dynamics. Furthermore, we demonstrated that the total time delay required for protein production and nuclear translocation of Mdm2 can induce p53 oscillations even when the p53 level is at a certain stable high steady state or at a certain stable low steady state. In addition, the two important features of the oscillatory dynamics-amplitude and period-can be controlled by such time delay. These results are in agreement with multiple experimental observations and may enrich our understanding of the dynamics of the p53 network.

Oscillatory behavior of p53-Mdm2 system driven by transcriptional and translational time delays

Biological experiments clarify that p53-Mdm2 module is the core of tumor network and p53 oscillation plays an important role in determining the tumor cell fate. In this paper, we investigate the effect of time delay on the oscillatory behavior induced by Hopf bifurcation in p53-Mdm2 system. First, the stability of the unique positive equilibrium point and the existence of Hopf bifurcation are investigated by using the time delay as the bifurcation parameter and by applying the bifurcation theory. Second, the explicit criteria determining the direction of Hopf bifurcation and the stability of bifurcating periodic solutions are developed based on the normal form theory and the center manifold theorem. In addition, the combination of numerical simulation results and theoretical calculation results indicates that time delays in p53-Mdm2 system are critical for p53 oscillations. The results may help us to better understand the biological functions of p53 pathway and provide clues for treatment of cancer.

A Gaussian process model and Bayesian variable selection for mapping function-valued quantitative traits with incomplete phenotypic data

Motivation

Recent advances in high dimensional phenotyping bring time as an extra dimension into the phenotypes. This promotes the quantitative trait locus (QTL) studies of function-valued traits such as those related to growth and development. Existing approaches for analyzing functional traits utilize either parametric methods or semi-parametric approaches based on splines and wavelets. However, very limited choices of software tools are currently available for practical implementation of functional QTL mapping and variable selection.

Results

We propose a Bayesian Gaussian process (GP) approach for functional QTL mapping. We use GPs to model the continuously varying coefficients which describe how the effects of molecular markers on the quantitative trait are changing over time. We use an efficient gradient based algorithm to estimate the tuning parameters of GPs. Notably, the GP approach is directly applicable to the incomplete datasets having even larger than 50% missing data rate (among phenotypes). We further develop a stepwise algorithm to search through the model space in terms of genetic variants, and use a minimal increase of Bayesian posterior probability as a stopping rule to focus on only a small set of putative QTL. We also discuss the connection between GP and penalized B-splines and wavelets. On two simulated and three real datasets, our GP approach demonstrates great flexibility for modeling different types of phenotypic trajectories with low computational cost. The proposed model selection approach finds the most likely QTL reliably in tested datasets.

Availability and implementation

Software and simulated data are available as a MATLAB package ‘GPQTLmapping’, and they can be downloaded from GitHub (https://github.com/jpvanhat/GPQTLmapping). Real datasets used in case studies are publicly available at QTL Archive.

Supplementary information

Supplementary data are available at Bioinformatics online.

The role of time delays in P53 gene regulatory network stimulated by growth factor

In this paper, a delayed mathematical model for the P53-Mdm2 network is developed. The P53-Mdm2 network we study is triggered by growth factor instead of DNA damage and the amount of DNA damage is regarded as zero. We study the influences of time delays, growth factor and other important chemical reaction rates on the dynamic behaviors in the system. It is shown that the time delay is a critical factor and its length determines the period, amplitude and stability of the P53 oscillation. Furthermore, as for some important chemical reaction rates, we also obtain some interesting results through numerical simulation. Especially, S (growth factor), k₃ (rate constant for Mdm2_p dephosphorylation), k₁₀ (basal expression of PTEN) and k₁₄ (Rate constant for PTEN-induced Akt dephosphorylation) could undermine the dynamic behavior of the system in different degree. These findings are expected to understand the mechanisms of action of several carcinogenic and tumor suppressor factors in humans under normal conditions.

Nutrient dose-responsive transcriptome changes driven by Michaelis-Menten kinetics underlie plant growth rates

An increase in nutrient dose leads to proportional increases in crop biomass and agricultural yield. However, the molecular underpinnings of this nutrient dose-response are largely unknown. To investigate, we assayed changes in the Arabidopsis root transcriptome to different doses of nitrogen (N)-a key plant nutrient-as a function of time. By these means, we found that rate changes of genome-wide transcript levels in response to N-dose could be explained by a simple kinetic principle: the Michaelis-Menten (MM) model. Fitting the MM model allowed us to estimate the maximum rate of transcript change (V _max), as well as the N-dose at which one-half of V _max was achieved (K _m) for 1,153 N-dose-responsive genes. Since transcription factors (TFs) can act in part as the catalytic agents that determine the rates of transcript change, we investigated their role in regulating N-dose-responsive MM-modeled genes. We found that altering the abundance of TGA1, an early N-responsive TF, perturbed the maximum rates of N-dose transcriptomic responses (V _max), K _m, as well as the rate of N-dose-responsive plant growth. We experimentally validated that MM-modeled N-dose-responsive genes included both direct and indirect TGA1 targets, using a root cell TF assay to detect TF binding and/or TF regulation genome-wide. Taken together, our results support a molecular mechanism of transcriptional control that allows an increase in N-dose to lead to a proportional change in the rate of genome-wide expression and plant growth.

BICORN: An R package for integrative inference of de novo cis-regulatory modules

Genome-wide transcription factor (TF) binding signal analyses reveal co-localization of TF binding sites based on inferred cis-regulatory modules (CRMs). CRMs play a key role in understanding the cooperation of multiple TFs under specific conditions. However, the functions of CRMs and their effects on nearby gene transcription are highly dynamic and context-specific and therefore are challenging to characterize. BICORN (Bayesian Inference of COoperative Regulatory Network) builds a hierarchical Bayesian model and infers context-specific CRMs based on TF-gene binding events and gene expression data for a particular cell type. BICORN automatically searches for a list of candidate CRMs based on the input TF bindings at regulatory regions associated with genes of interest. Applying Gibbs sampling, BICORN iteratively estimates model parameters of CRMs, TF activities, and corresponding regulation on gene transcription, which it models as a sparse network of functional CRMs regulating target genes. The BICORN package is implemented in R (version 3.4 or later) and is publicly available on the CRAN server at https://cran.r-project.org/web/packages/BICORN/index.html.

Stochastic simulation and statistical inference platform for visualization and estimation of transcriptional kinetics

Recent advances in single-molecule fluorescent imaging have enabled quantitative measurements of transcription at a single gene copy, yet an accurate understanding of transcriptional kinetics is still lacking due to the difficulty of solving detailed biophysical models. Here we introduce a stochastic simulation and statistical inference platform for modeling detailed transcriptional kinetics in prokaryotic systems, which has not been solved analytically. The model includes stochastic two-state gene activation, mRNA synthesis initiation and stepwise elongation, release to the cytoplasm, and stepwise co-transcriptional degradation. Using the Gillespie algorithm, the platform simulates nascent and mature mRNA kinetics of a single gene copy and predicts fluorescent signals measurable by time-lapse single-cell mRNA imaging, for different experimental conditions. To approach the inverse problem of estimating the kinetic parameters of the model from experimental data, we develop a heuristic optimization method based on the genetic algorithm and the empirical distribution of mRNA generated by simulation. As a demonstration, we show that the optimization algorithm can successfully recover the transcriptional kinetics of simulated and experimental gene expression data. The platform is available as a MATLAB software package at https://data.caltech.edu/records/1287.

A Practicable Method of Tuning the Noise Intensity at Protein Level

The expression of genes is inevitably subject to intracellular noise. Noise, for some regulatory networks, is constructive but detrimental to many others. The intensity of the noise is a determinant factor and the method of tuning it is of great value. In this study, we illustrated that the transcriptional delay in an incoherent feedforward loop (FFL) grants the target protein modulation the intensity of noise. Remarkably, for a wide range, the coefficient of variation (COV) of the target protein appeared to be about linear to the time span of the transcriptional delay. Without a noise-buffering method, the COV of the target protein is 0.455. While applying incoherent FFL, the COV reduced to 0.236. Then, it changed from 0.236 to 0.630 as the transcriptional delay raised from 0 to 1000 seconds. If we further increased the delay out of the linear range, the COV finally reached 0.779. In addition, we incorporated the distribution of the transcriptional delay in the delay stochastic simulation algorithm. This distribution is based on the experimental observation in the literature. The outcome suggested that the distributed delay slightly improved the ability of tuning noise. In conclusion, we demonstrated a noise-tuning method that altered only the intensity of noise without changing the deterministic steady-state behaviors. It is ready to be applied to various systems in the field of synthetic biology.

Osteoglycin (OGN) Inhibits Cell Proliferation and Invasiveness in Breast Cancer via PI3K/Akt/mTOR Signaling Pathway

Introduction

Previous studies have indicated that the small leucine-rich proteoglycan (SLR) osteoglycin (OGN) is downregulated in various cancers, including squamous cervical carcinoma, gastric cancer, and colorectal adenoma, indicating that OGN is a putative tumor suppressor. However, its exact role in the pathology of human cancers, especially breast cancer (BC), is not clear.

Methods

The expression of OGN in BC tissues was examined using qRT-PCR. Online databases were employed to analyze the correlation between OGN expression and clinicopathological characteristics. CCK-8 assay, colony formation assay, transwell migration and invasion assays were applied to detect cell proliferation, colony formation, migration and invasion of BC cells, respectively. Xenograft tumor models were constructed to explore the role of OGN on tumor growth in vivo.

Results

OGN expression was reduced in 24 paired BC samples compared with normal tissue. Decreased expression of OGN was correlated with greater pathological grade, a more aggressive tumor subtype, and poor overall survival. In vitro experiments showed that OGN overexpressed by plasmid transfection significantly inhibited cell proliferation, colony formation, migration, and invasion of BC cell lines. In xenograft tumor models, overexpression of OGN repressed the growth of MCF-7 cells in vivo and alleviated the compression of the tumor on surrounding structures. We also observed that OGN expression reversed EMT via repressing the PI3K/Akt/mTOR pathway.

Conclusion

This study revealed that OGN could function as a tumor suppressor during breast carcinogenesis, and we contribute new evidence to the body of research on the SLRP family.

Oscillation induced by Hopf bifurcation in the p53-Mdm2 feedback module

This study develops an integrated model of the p53-Mdm2 interaction composed of five basic components and time delay in the DNA damage response based on the existing research work. Some critical factors, including time delay, system parameters, and their interactions in the p53-Mdm2 system are investigated to examine their effects on the oscillatory behaviour induced by Hopf bifurcation. It is shown that the positive feedback formed between p53 and the activity of Mdm2 in the cytoplasm can cause a slight decrease in the amplitude of the p53 oscillation. The length of the time delay plays an important role in determining the amplitude and period of the oscillation and can significantly extend the parameter range for the system to demonstrate oscillatory behaviour. The numerical simulation results are found to be in good agreement with the published experimental observation. It is expected that the results of this research would be helpful to better understand the biological functions of p53 pathway and provide some clues in the treatment of cancer.

RNA components of the spliceosome regulate tissue- and cancer-specific alternative splicing

Alternative splicing of pre-mRNAs plays a pivotal role during the establishment and maintenance of human cell types. Characterizing the trans-acting regulatory proteins that control alternative splicing has therefore been the focus of much research. Recent work has established that even core protein components of the spliceosome, which are required for splicing to proceed, can nonetheless contribute to splicing regulation by modulating splice site choice. We here show that the RNA components of the spliceosome likewise influence alternative splicing decisions. Although these small nuclear RNAs (snRNAs), termed U1, U2, U4, U5, and U6 snRNA, are present in equal stoichiometry within the spliceosome, we found that their relative levels vary by an order of magnitude during development, across tissues, and across cancer samples. Physiologically relevant perturbation of individual snRNAs drove widespread gene-specific differences in alternative splicing but not transcriptome-wide splicing failure. Genes that were particularly sensitive to variations in snRNA abundance in a breast cancer cell line model were likewise preferentially misspliced within a clinically diverse cohort of invasive breast ductal carcinomas. As aberrant mRNA splicing is prevalent in many cancers, we propose that a full understanding of such dysregulated pre-mRNA processing requires study of snRNAs, as well as protein splicing factors. Together, our data show that the RNA components of the spliceosome are not merely basal factors, as has long been assumed. Instead, these noncoding RNAs constitute a previously uncharacterized layer of regulation of alternative splicing, and contribute to the establishment of global splicing programs in both healthy and malignant cells.

Temporal dynamic reorganization of 3D chromatin architecture in hormone-induced breast cancer and endocrine resistance

Recent studies have demonstrated that chromatin architecture is linked to the progression of cancers. However, the roles of 3D structure and its dynamics in hormone-dependent breast cancer and endocrine resistance are largely unknown. Here we report the dynamics of 3D chromatin structure across a time course of estradiol (E2) stimulation in human estrogen receptor α (ERα)-positive breast cancer cells. We identified subsets of temporally highly dynamic compartments predominantly associated with active open chromatin and found that these highly dynamic compartments showed higher alteration in tamoxifen-resistant breast cancer cells. Remarkably, these compartments are characterized by active chromatin states, and enhanced ERα binding but decreased transcription factor CCCTC-binding factor (CTCF) binding. We finally identified a set of ERα-bound promoter-enhancer looping genes enclosed within altered domains that are enriched with cancer invasion, aggressiveness or metabolism signaling pathways. This large-scale analysis expands our understanding of high-order temporal chromatin reorganization underlying hormone-dependent breast cancer.

Theoretical study on the oscillation mechanism of p53-Mdm2 network

In this paper, a delayed mathematical model was developed based on experimental data to understand how the time delays required for transcription and translation in Mdm2 gene expression affect the kinetic behavior of the p53-Mdm2 network. Taking the time delays as the main research parameters, the stability of the system at the positive equilibrium was studied by using the theoretical method of delay differential equation. We found that such delays can induce oscillations by undergoing a supercritical Hopf bifurcation. Then, we used the normal form theory and the center manifold reduction to study the direction and stability of the bifurcation in detail. Furthermore, we also studied the effects of the length of time delays and the model parameters by numerical simulations. We found that time delays in Mdm2 synthesis are required for p53 oscillations and the length of such delays can determine the amplitude and period of the oscillations. In addition, the model parameters can also change the stability of the system. These results illustrate that the repair process after DNA damage can be regulated by varying time delays and the model parameters.

Genome-wide Estrogen Receptor-α activation is sustained, not cyclical

Estrogen Receptor-alpha (ER) drives 75% of breast cancers. Stimulation of the ER by estra-2-diol forms a transcriptionally-active chromatin-bound complex. Previous studies reported that ER binding follows a cyclical pattern. However, most studies have been limited to individual ER target genes and without replicates. Thus, the robustness and generality of ER cycling are not well understood. We present a comprehensive genome-wide analysis of the ER after activation, based on 6 replicates at 10 time-points, using our method for precise quantification of binding, Parallel-Factor ChIP-seq. In contrast to previous studies, we identified a sustained increase in affinity, alongside a class of estra-2-diol independent binding sites. Our results are corroborated by quantitative re-analysis of multiple independent studies. Our new model reconciles the conflicting studies into the ER at the TFF1 promoter and provides a detailed understanding in the context of the ER's role as both the driver and therapeutic target of breast cancer.

GPrank: an R package for detecting dynamic elements from genome-wide time series

Background

Genome-wide high-throughput sequencing (HTS) time series experiments are a powerful tool for monitoring various genomic elements over time. They can be used to monitor, for example, gene or transcript expression with RNA sequencing (RNA-seq), DNA methylation levels with bisulfite sequencing (BS-seq), or abundances of genetic variants in populations with pooled sequencing (Pool-seq). However, because of high experimental costs, the time series data sets often consist of a very limited number of time points with very few or no biological replicates, posing challenges in the data analysis.

Results

Here we present the GPrank R package for modelling genome-wide time series by incorporating variance information obtained during pre-processing of the HTS data using probabilistic quantification methods or from a beta-binomial model using sequencing depth. GPrank is well-suited for analysing both short and irregularly sampled time series. It is based on modelling each time series by two Gaussian process (GP) models, namely, time-dependent and time-independent GP models, and comparing the evidence provided by data under two models by computing their Bayes factor (BF). Genomic elements are then ranked by their BFs, and temporally most dynamic elements can be identified.

Conclusions

Incorporating the variance information helps GPrank avoid false positives without compromising computational efficiency. Fitted models can be easily further explored in a browser. Detection and visualisation of temporally most active dynamic elements in the genome can provide a good starting point for further downstream analyses for increasing our understanding of the studied processes.

SparseIso: a novel Bayesian approach to identify alternatively spliced isoforms from RNA-seq data

Motivation

Recent advances in high-throughput RNA sequencing (RNA-seq) technologies have made it possible to reconstruct the full transcriptome of various types of cells. It is important to accurately assemble transcripts or identify isoforms for an improved understanding of molecular mechanisms in biological systems.

Results

We have developed a novel Bayesian method, SparseIso, to reliably identify spliced isoforms from RNA-seq data. A spike-and-slab prior is incorporated into the Bayesian model to enforce the sparsity for isoform identification, effectively alleviating the problem of overfitting. A Gibbs sampling procedure is further developed to simultaneously identify and quantify transcripts from RNA-seq data. With the sampling approach, SparseIso estimates the joint distribution of all candidate transcripts, resulting in a significantly improved performance in detecting lowly expressed transcripts and multiple expressed isoforms of genes. Both simulation study and real data analysis have demonstrated that the proposed SparseIso method significantly outperforms existing methods for improved transcript assembly and isoform identification.

Availability and implementation

The SparseIso package is available at http://github.com/henryxushi/SparseIso.

Contact

xuan@vt.edu.

Supplementary information

Supplementary data are available at Bioinformatics online.

Hamiltonian Monte Carlo Inversion of Seismic Sources in Complex Media

We present a probabilistic seismic point source inversion, taking into account 3-D heterogeneous Earth structure. Our method rests on (1) reciprocity and numerical wavefield simulations in complex media and (2) Hamiltonian Monte Carlo sampling that requires only a small amount of test models to provide reliable uncertainty information on the timing, location, and mechanism of the source. Using spectral element simulations of 3-D, viscoelastic, anisotropic wave propagation, we precompute receiver side strain tensors in time and space. This enables the fast computation of synthetic seismograms for any hypothetical source within the volume of interest, and thus a Bayesian solution of the inverse problem. To improve efficiency, we developed a variant of Hamiltonian Monte Carlo sampling. Taking advantage of easily computable derivatives, numerical examples indicate that Hamiltonian Monte Carlo can converge to the posterior probability density with orders of magnitude less samples than the derivative-free Metropolis-Hastings algorithm, which we use for benchmarking. Exact numbers depend on observational errors and the quality of the prior. We apply our method to the Japanese Islands region where we previously constrained 3-D structure of the crust and upper mantle using full-waveform inversion with a minimum period of 15 s.

Hypomethylation of the miRNA-34a gene promoter is associated with Severe Preeclampsia

PURPOSE

PE is a pregnancy-specific complication, which genetic and epigenetic factors play key roles in its pathogenesis. DNA methylation is a main epigenetic alteration with important roles in gene regulation. Micro RNAs (miRNAs) as another member of epigenetic machinery regulate the gene expression and involve in different biological pathways including apoptosis and placental development. Therefore, the present study performed to assess the association between miRNA-34a promoter methylation and PE susceptibility.

METHODS

The placenta of 104 PE pregnant women and 119 normotensive pregnant women were collected after delivery. The methylation status of the miRNA-34a promoter was assessed using Methylation Specific PCR (MSP).

RESULTS

The frequency of the hemi-methylated (MU) miR-34a promoter was significantly lower in PE women compared to the controls (17.3 vs. 29.4%) (OR, 0.45 [95% CI, 0.2-0.9], P = 0.016). The overall methylation rate was 23.1% in PE women and 41.2% in the control group and was significantly lower in PE women (OR, 0.4 [95% CI, 0.2-0.8], P = 0.004). The frequency of hemi-methylated (MU) and overall methylated (MU+MM) promoter of miR-34a gene was significantly lower in severe PE but not in mild PE women compared to the controls [(OR, 0.3 [95% CI, 0.1-0.8], P = 0.02) and (OR, 0.3 [95% CI, 0.1-0.7], P = 0.009), respectively]. There was an association between hemi-methylated (MU) and overall methylated (MU+MM) promoter and late onset PE [(OR, 0.4 [95% CI, 0.2-0.9], P = 0.03) and (OR, 0.4 [95% CI, 0.2-0.8], P = 0.01), respectively].

CONCLUSIONS

An association was found between hypo-methylation of the miR-34a promoter and PE and PE severity.

Luminal lncRNAs Regulation by ERα-Controlled Enhancers in a Ligand-Independent Manner in Breast Cancer Cells

Estrogen receptor-α (ERα) is a ligand-inducible protein which mediates estrogenic hormones signaling and defines the luminal BC phenotype. Recently, we demonstrated that even in absence of ligands ERα (apoERα) binds chromatin sites where it regulates transcription of several protein-coding and lncRNA genes. Noteworthy, apoERα-regulated lncRNAs marginally overlap estrogen-induced transcripts, thus representing a new signature of luminal BC genes. By the analysis of H3K27ac enrichment in hormone-deprived MCF-7 cells, we defined a set of Super Enhancers (SEs) occupied by apoERα, including one mapped in proximity of the DSCAM-AS1 lncRNA gene. This represents a paradigm of apoERα activity since its expression is largely unaffected by estrogenic treatment, despite the fact that E2 increases ERα binding on DSCAM-AS1 promoter. We validated the enrichment of apoERα, p300, GATA3, FoxM1 and CTCF at both DSCAM-AS1 TSS and at its associated SE by ChIP-qPCR. Furthermore, by analyzing MCF-7 ChIA-PET data and by 3C assays, we confirmed long range chromatin interaction between the SE and the DSCAM-AS1 TSS. Interestingly, CTCF and p300 binding showed an enrichment in hormone-depleted medium and in the presence of ERα, elucidating the dynamics of the estrogen-independent regulation of DSCAM-AS1 expression. The analysis of this lncRNA provides a paradigm of transcriptional regulation of a luminal specific apoERα regulated lncRNA.

Predicting stimulation-dependent enhancer-promoter interactions from ChIP-Seq time course data

We have developed a machine learning approach to predict stimulation-dependent enhancer-promoter interactions using evidence from changes in genomic protein occupancy over time. The occupancy of estrogen receptor alpha (ERα), RNA polymerase (Pol II) and histone marks H2AZ and H3K4me3 were measured over time using ChIP-Seq experiments in MCF7 cells stimulated with estrogen. A Bayesian classifier was developed which uses the correlation of temporal binding patterns at enhancers and promoters and genomic proximity as features to predict interactions. This method was trained using experimentally determined interactions from the same system and was shown to achieve much higher precision than predictions based on the genomic proximity of nearest ERα binding. We use the method to identify a genome-wide confident set of ERα target genes and their regulatory enhancers genome-wide. Validation with publicly available GRO-Seq data demonstrates that our predicted targets are much more likely to show early nascent transcription than predictions based on genomic ERα binding proximity alone.

Oscillatory dynamics of p38 activity with transcriptional and translational time delays

Recent experimental evidence reports that oscillations of p38 MAPK (p38) activity would efficiently induce pro-inflammatory gene expression, which might be deleterious to immune systems and may even cause cellular damage and apoptosis. It is widely accepted now that transcriptional and translational delays are ubiquitous in gene expression, which can typically result in oscillatory responses of gene regulations. Consequently, delay-driven sustained oscillations in p38 activity (p38*) could in principle be commonplace. Nevertheless, so far the studies of the impact of such delays on p38* have been lacking both experimentally and theoretically. Here, we use experimental data to develop a delayed mathematical model, with the aim of understanding how such delays affect oscillatory behaviour on p38*. We analyze the stability and oscillation of the model with and without explicit time delays. We show that a sufficiently input stimulation strength is prerequisite for generating p38* oscillations, and that an optimal rate of model parameters is also essential to these oscillations. Moreover, we find that the time delays required for transcription and translation in mitogen-activated protein kinase phosphatase-1 (MKP-1) gene expression can drive p38* to be oscillatory even when the concentration of p38* level is at a stable state. Furthermore, the length of these delays can determine the amplitude and period of the oscillations and can enormously extend the oscillatory ranges of model parameters. These results indicate that time delays in MKP-1 synthesis are required, albeit not sufficient, for p38* oscillations, which may lead to new insights related to p38 oscillations.

Analysis of differential splicing suggests different modes of short-term splicing regulation

Motivation: Alternative splicing is an important mechanism in which the regions of pre-mRNAs are differentially joined in order to form different transcript isoforms. Alternative splicing is involved in the regulation of normal physiological functions but also linked to the development of diseases such as cancer. We analyse differential expression and splicing using RNA-sequencing time series in three different settings: overall gene expression levels, absolute transcript expression levels and relative transcript expression levels.

Results: Using estrogen receptor α signaling response as a model system, our Gaussian process-based test identifies genes with differential splicing and/or differentially expressed transcripts. We discover genes with consistent changes in alternative splicing independent of changes in absolute expression and genes where some transcripts change whereas others stay constant in absolute level. The results suggest classes of genes with different modes of alternative splicing regulation during the experiment.

Availability and Implementation: R and Matlab codes implementing the method are available at https://github.com/PROBIC/diffsplicing. An interactive browser for viewing all model fits is available at http://users.ics.aalto.fi/hande/splicingGP/

Contact:hande.topa@helsinki.fi or antti.honkela@helsinki.fi

Supplementary information:Supplementary data are available at Bioinformatics online.

Associating transcription factors and conserved RNA structures with gene regulation in the human brain

Anatomical subdivisions of the human brain can be associated with different neuronal functions. This functional diversification is reflected by differences in gene expression. By analyzing post-mortem gene expression data from the Allen Brain Atlas, we investigated the impact of transcription factors (TF) and RNA secondary structures on the regulation of gene expression in the human brain. First, we modeled the expression of a gene as a linear combination of the expression of TFs. We devised an approach to select robust TF-gene interactions and to determine localized contributions to gene expression of TFs. Among the TFs with the most localized contributions, we identified EZH2 in the cerebellum, NR3C1 in the cerebral cortex and SRF in the basal forebrain. Our results suggest that EZH2 is involved in regulating ZIC2 and SHANK1 which have been linked to neurological diseases such as autism spectrum disorder. Second, we associated enriched regulatory elements inside differentially expressed mRNAs with RNA secondary structure motifs. We found a group of purine-uracil repeat RNA secondary structure motifs plus other motifs in neuron related genes such as ACSL4 and ERLIN2.

Kinetic CRAC uncovers a role for Nab3 in determining gene expression profiles during stress

RNA-binding proteins play a key role in shaping gene expression profiles during stress, however, little is known about the dynamic nature of these interactions and how this influences the kinetics of gene expression. To address this, we developed kinetic cross-linking and analysis of cDNAs (χCRAC), an ultraviolet cross-linking method that enabled us to quantitatively measure the dynamics of protein-RNA interactions in vivo on a minute time-scale. Here, using χCRAC we measure the global RNA-binding dynamics of the yeast transcription termination factor Nab3 in response to glucose starvation. These measurements reveal rapid changes in protein-RNA interactions within 1 min following stress imposition. Changes in Nab3 binding are largely independent of alterations in transcription rate during the early stages of stress response, indicating orthogonal transcriptional control mechanisms. We also uncover a function for Nab3 in dampening expression of stress-responsive genes. χCRAC has the potential to greatly enhance our understanding of in vivo dynamics of protein-RNA interactions.Protein RNA interactions are dynamic and regulated in response to environmental changes. Here the authors describe 'kinetic CRAC', an approach that allows time resolved analyses of protein RNA interactions with minute time point resolution and apply it to gain insight into the function of the RNA-binding protein Nab3.

Transient dynamics and their control in time-delay autonomous Boolean ring networks

Biochemical systems with switch-like interactions, such as gene regulatory networks, are well modeled by autonomous Boolean networks. Specifically, the topology and logic of gene interactions can be described by systems of continuous piecewise-linear differential equations, enabling analytical predictions of the dynamics of specific networks. However, most models do not account for time delays along links associated with spatial transport, mRNA transcription, and translation. To address this issue, we have developed an experimental test bed to realize a time-delay autonomous Boolean network with three inhibitory nodes, known as a repressilator, and use it to study the dynamics that arise as time delays along the links vary. We observe various nearly periodic oscillatory transient patterns with extremely long lifetime, which emerge in small network motifs due to the delay, and which are distinct from the eventual asymptotically stable periodic attractors. For repeated experiments with a given network, we find that stochastic processes give rise to a broad distribution of transient times with an exponential tail. In some cases, the transients are so long that it is doubtful the attractors will ever be approached in a biological system that has a finite lifetime. To counteract the long transients, we show experimentally that small, occasional perturbations applied to the time delays can force the trajectories to rapidly approach the attractors.

Respectful Modeling: Addressing Uncertainty in Dynamic System Models for Molecular Biology

Although there is still some skepticism in the biological community regarding the value and significance of quantitative computational modeling, important steps are continually being taken to enhance its accessibility and predictive power. We view these developments as essential components of an emerging 'respectful modeling' framework which has two key aims: (i) respecting the models themselves and facilitating the reproduction and update of modeling results by other scientists, and (ii) respecting the predictions of the models and rigorously quantifying the confidence associated with the modeling results. This respectful attitude will guide the design of higher-quality models and facilitate the use of models in modern applications such as engineering and manipulating microbial metabolism by synthetic biology.

Genes, epigenetics and miRNA regulation in the placenta

This text reviews briefly the context in which epigenetics regulate gene expression in trophoblast development and function. It is an attempt to focus on a limited number of recent papers that, according to the author, shed new light on placental development, and constitute possible trails for improving knowledge and women follow-up in pathological pregnancies.

Perfect timing: splicing and transcription rates in living cells

An important step toward understanding gene regulation is the elucidation of the time necessary for the completion of individual steps. Measurement of reaction rates can reveal potential nodes for regulation. For example, measurements of in vivo transcription elongation rates reveal regulation by DNA sequence, gene architecture, and chromatin. Pre-mRNA splicing is regulated by transcription elongation rates and vice versa, yet the rates of RNA processing reactions remain largely elusive. Since the 1980s, numerous model systems and approaches have been used to determine the precise timing of splicing in vivo. Because splicing can be co-transcriptional, the position of Pol II when splicing is detected has been used as a proxy for time by some investigators. In addition to these 'distance-based' measurements, 'time-based' measurements have been possible through live cell imaging, metabolic labeling of RNA, and gene induction. Yet splicing rates can be convolved by the time it takes for transcription, spliceosome assembly and spliceosome disassembly. The variety of assays and systems used has, perhaps not surprisingly, led to reports of widely differing splicing rates in vivo. Recently, single molecule RNA-seq has indicated that splicing occurs more quickly than previously deduced. Here we comprehensively review these findings and discuss evidence that splicing and transcription rates are closely coordinated, facilitating the efficiency of gene expression. On the other hand, introduction of splicing delays through as yet unknown mechanisms provide opportunity for regulation. More work is needed to understand how cells optimize the rates of gene expression for a range of biological conditions. WIREs RNA 2017, 8:e1401. doi: 10.1002/wrna.1401 For further resources related to this article, please visit the WIREs website.