Protocol variations in run-on transcription dataset preparation produce detectable signatures in sequencing libraries.

BACKGROUND: A variety of protocols exist for producing whole genome run-on transcription datasets. However, little is known about how differences between these protocols affect the signal within the resulting libraries. RESULTS: Using run-on transcription datasets generated from the same biological system, we show that a variety of GRO- and PRO-seq preparation methods leave identifiable signatures within each library. Specifically we show that the library preparation method results in differences in quality control metrics, as well as differences in the signal distribution at the 5 ' end of transcribed regions. These shifts lead to disparities in eRNA identification, but do not impact analyses aimed at inferring the key regulators involved in changes to transcription. CONCLUSIONS: Run-on sequencing protocol variations result in technical signatures that can be used to identify both the enrichment and library preparation method of a particular data set. These technical signatures are batch effects that limit detailed comparisons of pausing ratios and eRNAs identified across protocols. However, these batch effects have only limited impact on our ability to infer which regulators underlie the observed transcriptional changes.

Determination of steady-state transcriptome modifications associated with repeated homotypic stress in the rat rostral posterior hypothalamic region.

Chronic stress is epidemiologically correlated with physical and psychiatric disorders. Whereas many animal models of chronic stress induce symptoms of psychopathology, repeated homotypic stressors to moderate intensity stimuli typically reduce stress-related responses with fewer, if any, pathological symptoms. Recent results indicate that the rostral posterior hypothalamic (rPH) region is a significant component of the brain circuitry underlying response reductions (habituation) associated with repeated homotypic stress. To test whether posterior hypothalamic transcriptional regulation associates with the neuroendocrine modifications induced by repeated homotypic stress, RNA-seq was performed in the rPH dissected from adult male rats that experienced either no stress, 1, 3, or 7 stressful loud noise exposures. Plasma samples displayed reliable increases of corticosterone in all stressed groups, with the smallest increase in the group exposed to 7 loud noises, indicating significant habituation compared to the other stressed groups. While few or no differentially expressed genes were detected 24-h after one or three loud noise exposures, relatively large numbers of transcripts were differentially expressed between the group exposed to 7 loud noises when compared to the control or 3-stress groups, respectively, which correlated with the corticosterone response habituation observed. Gene ontology analyses indicated multiple significant functional terms related to neuron differentiation, neural membrane potential, pre- and post-synaptic elements, chemical synaptic transmission, vesicles, axon guidance and projection, glutamatergic and GABAergic neurotransmission. Some of the differentially expressed genes (Myt1l, Zmat4, Dlx6, Csrnp3) encode transcription factors that were independently predicted by transcription factor enrichment analysis to target other differentially regulated genes in this study. A similar experiment employing in situ hybridization histochemical analysis in additional animals validated the direction of change of the 5 transcripts investigated (Camk4, Gabrb2, Gad1, Grin2a and Slc32a) with a high level of temporal and regional specificity for the rPH. In aggregate, the results suggest that distinct patterns of gene regulation are obtained in response to a repeated homotypic stress regimen; they also point to a significant reorganization of the rPH region that may critically contribute to the phenotypic modifications associated with repeated homotypic stress habituation.

Atlas of nascent RNA transcripts reveals enhancer to gene linkages.

Gene transcription is controlled and modulated by regulatory regions, including enhancers and promoters. These regions are abundant in unstable, non-coding bidirectional transcription. Using nascent RNA transcription data across hundreds of human samples, we identified over 800,000 regions containing bidirectional transcription. We then identify highly correlated transcription between bidirectional and gene regions. The identified correlated pairs, a bidirectional region and a gene, are enriched for disease associated SNPs and often supported by independent 3D data. We present these resources as an SQL database which serves as a resource for future studies into gene regulation, enhancer associated RNAs, and transcription factors.

Aging disrupts gene expression timing during muscle regeneration.

Skeletal muscle function and regenerative capacity decline during aging, yet factors driving these changes are incompletely understood. Muscle regeneration requires temporally coordinated transcriptional programs to drive myogenic stem cells to activate, proliferate, fuse to form myofibers, and to mature as myonuclei, restoring muscle function after injury. We assessed global changes in myogenic transcription programs distinguishing muscle regeneration in aged mice from young mice by comparing pseudotime trajectories from single-nucleus RNA sequencing of myogenic nuclei. Aging-specific differences in coordinating myogenic transcription programs necessary for restoring muscle function occur following muscle injury, likely contributing to compromised regeneration in aged mice. Differences in pseudotime alignment of myogenic nuclei when comparing aged with young mice via dynamic time warping revealed pseudotemporal differences becoming progressively more severe as regeneration proceeds. Disruptions in timing of myogenic gene expression programs may contribute to incomplete skeletal muscle regeneration and declines in muscle function as organisms age.

Transcription factor enrichment analysis (TFEA) quantifies the activity of multiple transcription factors from a single experiment.

Detecting changes in the activity of a transcription factor (TF) in response to a perturbation provides insights into the underlying cellular process. Transcription Factor Enrichment Analysis (TFEA) is a robust and reliable computational method that detects positional motif enrichment associated with changes in transcription observed in response to a perturbation. TFEA detects positional motif enrichment within a list of ranked regions of interest (ROIs), typically sites of RNA polymerase initiation inferred from regulatory data such as nascent transcription. Therefore, we also introduce muMerge, a statistically principled method of generating a consensus list of ROIs from multiple replicates and conditions. TFEA is broadly applicable to data that informs on transcriptional regulation including nascent transcription (eg. PRO-Seq), CAGE, histone ChIP-Seq, and accessibility data (e.g., ATAC-Seq). TFEA not only identifies the key regulators responding to a perturbation, but also temporally unravels regulatory networks with time series data. Consequently, TFEA serves as a hypothesis-generating tool that provides an easy, rigorous, and cost-effective means to broadly assess TF activity yielding new biological insights.

Two-stage ML Classifier for Identifying Host Protein Targets of the Dengue Protease.

Flaviviruses such as dengue encode a protease that is essential for viral replication. The protease functions by cleaving well-conserved positions in the viral polyprotein. In addition to the viral polyprotein, the dengue protease cleaves at least one host protein involved in immune response. This raises the question, what other host proteins are targeted and cleaved? Here we present a new computational method for identifying putative host protein targets of the dengue virus protease. Our method relies on biochemical and secondary structure features at the known cleavage sites in the viral polyprotein in a two-stage classification process to identify putative cleavage targets. The accuracy of our predictions scaled inversely with evolutionary distance when we applied it to the known cleavage sites of several other flaviviruses-a good indication of the validity of our predictions. Ultimately, our classifier identified 257 human protein sites possessing both a similar target motif and accessible local structure. These proteins are promising candidates for further investigation. As the number of viral sequences expands, our method could be adopted to predict host targets of other flaviviruses.

Stress relaxation in quasi-two-dimensional self-assembled nanoparticle monolayers.

We experimentally probed the stress relaxation of a monolayer of iron oxide nanoparticles at the water-air interface. Upon drop-casting onto a water surface, the nanoparticles self-assembled into islands of two-dimensional hexagonally close packed crystalline domains surrounded by large voids. When compressed laterally, the voids gradually disappeared as the surface pressure increased. After the compression was stopped, the surface pressure (as measured by a Wilhelmy plate) evolved as a function of the film aging time with three distinct timescales. These aging dynamics were intrinsic to the stressed state built up during the non-equilibrium compression of the film. Utilizing x-ray photon correlation spectroscopy, we measured the characteristic relaxation time (τ) of in-plane nanoparticle motion as a function of the aging time through both second-order and two-time autocorrelation analysis. Compressed and stretched exponential fitting of the intermediate scattering function yielded exponents (ß) indicating different relaxation mechanisms of the films under different compression stresses. For a monolayer compressed to a lower surface pressure (between 20 mN/m and 30 mN/m), the relaxation time (τ) decreased continuously as a function of the aging time, as did the fitted exponent, which transitioned from being compressed (>1) to stretched (<1), indicating that the monolayer underwent a stress release through crystalline domain reorganization. However, for a monolayer compressed to a higher surface pressure (around 40 mN/m), the relaxation time increased continuously and the compressed exponent varied very little from a value of 1.6, suggesting that the system may have been highly stressed and jammed. Despite the interesting stress relaxation signatures seen in these samples, the structural ordering of the monolayer remained the same over the sample lifetime, as revealed by grazing incidence x-ray diffraction.