GOLGA8 increases bulk antisense oligonucleotide uptake and activity in mammalian cells.

Antisense oligonucleotides (ASOs) are short synthetic nucleic acids that recognize and bind to complementary RNA to modulate gene expression. It is well established that single-stranded, phosphorothioate-modified ASOs enter cells independent of carrier molecules, primarily via endocytic pathways, but that only a small portion of internalized ASO is released into the cytosol and/or nucleus, rendering the majority of ASO inaccessible to the targeted RNA. Identifying pathways that can increase the available ASO pool is valuable as a research tool and therapeutically. Here, we conducted a functional genomic screen for ASO activity by engineering GFP splice reporter cells and applying genome-wide CRISPR gene activation. The screen can identify factors that enhance ASO splice modulation activity. Characterization of hit genes uncovered GOLGA8, a largely uncharacterized protein, as a novel positive regulator enhancing ASO activity by ∼2-fold. Bulk ASO uptake is 2- to 5-fold higher in GOLGA8-overexpressing cells where GOLGA8 and ASOs are observed in the same intracellular compartments. We find GOLGA8 is highly localized to the trans-Golgi and readily detectable at the plasma membrane. Interestingly, overexpression of GOLGA8 increased activity for both splice modulation and RNase H1-dependent ASOs. Taken together, these results support a novel role for GOLGA8 in productive ASO uptake.

Chemical modification of PS-ASO therapeutics reduces cellular protein-binding and improves the therapeutic index.

The molecular mechanisms of toxicity of chemically modified phosphorothioate antisense oligonucleotides (PS-ASOs) are not fully understood. Here, we report that toxic gapmer PS-ASOs containing modifications such as constrained ethyl (cEt), locked nucleic acid (LNA) and 2'-O-methoxyethyl (2'-MOE) bind many cellular proteins with high avidity, altering their function, localization and stability. We show that RNase H1-dependent delocalization of paraspeckle proteins to nucleoli is an early event in PS-ASO toxicity, followed by nucleolar stress, p53 activation and apoptotic cell death. Introduction of a single 2'-O-methyl (2'-OMe) modification at gap position 2 reduced protein-binding, substantially decreasing hepatotoxicity and improving the therapeutic index with minimal impairment of antisense activity. We validated the ability of this modification to generally mitigate PS-ASO toxicity with more than 300 sequences. Our findings will guide the design of PS-ASOs with optimal therapeutic profiles.

A modular analysis of microglia gene expression, insights into the aged phenotype.

BACKGROUND: Microglia are multifunctional cells that are key players in brain development and homeostasis. Recent years have seen tremendous growth in our understanding of the role microglia play in neurodegeneration, CNS injury, and developmental disorders. Given that microglia show diverse functional phenotypes, there is a need for more precise tools to characterize microglial states. Here, we experimentally define gene modules as the foundation for describing microglial functional states. RESULTS: In an effort to develop a comprehensive classification scheme, we profiled transcriptomes of mouse microglia in a stimulus panel with 96 different conditions. Using the transcriptomic data, we generated fine-resolution gene modules that are robustly preserved across datasets. These modules served as the basis for a combinatorial code that we then used to characterize microglial activation under various inflammatory stimulus conditions. CONCLUSIONS: The microglial gene modules described here were robustly preserved, and could be applied to in vivo as well as in vitro conditions to dissociate the signaling pathways that distinguish acutely inflamed microglia from aged microglia. The microglial gene modules presented here are a novel resource for classifying and characterizing microglial states in health and disease.

Quantitative analysis of chromatin interaction changes upon a 4.3 Mb deletion at mouse 4E2.

BACKGROUND: Circular chromosome conformation capture (4C) has provided important insights into three dimensional (3D) genome organization and its critical impact on the regulation of gene expression. We developed a new quantitative framework based on polymer physics for the analysis of paired-end sequencing 4C (PE-4Cseq) data. We applied this strategy to the study of chromatin interaction changes upon a 4.3 Mb DNA deletion in mouse region 4E2. RESULTS: A significant number of differentially interacting regions (DIRs) and chromatin compaction changes were detected in the deletion chromosome compared to a wild-type (WT) control. Selected DIRs were validated by 3D DNA FISH experiments, demonstrating the robustness of our pipeline. Interestingly, significant overlaps of DIRs with CTCF/Smc1 binding sites and differentially expressed genes were observed. CONCLUSIONS: Altogether, our PE-4Cseq analysis pipeline provides a comprehensive characterization of DNA deletion effects on chromatin structure and function.

Breaking an epigenetic chromatin switch: curious features of hysteresis in Saccharomyces cerevisiae telomeric silencing.

In addition to gene network switches, local epigenetic modifications to DNA and histones play an important role in all-or-none cellular decision-making. Here, we study the dynamical design of a well-characterized epigenetic chromatin switch: the yeast SIR system, in order to understand the origin of the stability of epigenetic states. We study hysteresis in this system by perturbing it with a histone deacetylase inhibitor. We find that SIR silencing has many characteristics of a non-linear bistable system, as observed in conventional genetic switches, which are based on activities of a few promoters affecting each other through the abundance of their gene products. Quite remarkably, our experiments in yeast telomeric silencing show a very distinctive pattern when it comes to the transition from bistability to monostability. In particular, the loss of the stable silenced state, upon increasing the inhibitor concentration, does not seem to show the expected saddle node behavior, instead looking like a supercritical pitchfork bifurcation. In other words, the 'off' state merges with the 'on' state at a threshold concentration leading to a single state, as opposed to the two states remaining distinct up to the threshold and exhibiting a discontinuous jump from the 'off' to the 'on' state. We argue that this is an inevitable consequence of silenced and active regions coexisting with dynamic domain boundaries. The experimental observations in our study therefore have broad implications for the understanding of chromatin silencing in yeast and beyond.

The role of multiple marks in epigenetic silencing and the emergence of a stable bivalent chromatin state.

We introduce and analyze a minimal model of epigenetic silencing in budding yeast, built upon known biomolecular interactions in the system. Doing so, we identify the epigenetic marks essential for the bistability of epigenetic states. The model explicitly incorporates two key chromatin marks, namely H4K16 acetylation and H3K79 methylation, and explores whether the presence of multiple marks lead to a qualitatively different systems behavior. We find that having both modifications is important for the robustness of epigenetic silencing. Besides the silenced and transcriptionally active fate of chromatin, our model leads to a novel state with bivalent (i.e., both active and silencing) marks under certain perturbations (knock-out mutations, inhibition or enhancement of enzymatic activity). The bivalent state appears under several perturbations and is shown to result in patchy silencing. We also show that the titration effect, owing to a limited supply of silencing proteins, can result in counter-intuitive responses. The design principles of the silencing system is systematically investigated and disparate experimental observations are assessed within a single theoretical framework. Specifically, we discuss the behavior of Sir protein recruitment, spreading and stability of silenced regions in commonly-studied mutants (e.g., sas2[Formula: see text], dot1[Formula: see text]) illuminating the controversial role of Dot1 in the systems biology of yeast silencing.

Internucleosomal interactions mediated by histone tails allow distant communication in chromatin.

Action across long distances on chromatin is a hallmark of eukaryotic transcriptional regulation. Although chromatin structure per se can support long-range interactions, the mechanisms of efficient communication between widely spaced DNA modules in chromatin remain a mystery. The molecular simulations described herein suggest that transient binary internucleosomal interactions can mediate distant communication in chromatin. Electrostatic interactions between the N-terminal tails of the core histones and DNA enhance the computed probability of juxtaposition of sites that lie far apart along the DNA sequence. Experimental analysis of the rates of communication in chromatin constructs confirms that long-distance communication occurs efficiently and independently of distance on tail-containing, but not on tailless, chromatin. Taken together, our data suggest that internucleosomal interactions involving the histone tails are essential for highly efficient, long-range communication between regulatory elements and their targets in eukaryotic genomes.

Theoretical analysis of the role of chromatin interactions in long-range action of enhancers and insulators.

Long-distance regulatory interactions between enhancers and their target genes are commonplace in higher eukaryotes. Interposed boundaries or insulators are able to block these long-distance regulatory interactions. The mechanistic basis for insulator activity and how it relates to enhancer action-at-a-distance remains unclear. Here we explore the idea that topological loops could simultaneously account for regulatory interactions of distal enhancers and the insulating activity of boundary elements. We show that while loop formation is not in itself sufficient to explain action at a distance, incorporating transient nonspecific and moderate attractive interactions between the chromatin fibers strongly enhances long-distance regulatory interactions and is sufficient to generate a euchromatin-like state. Under these same conditions, the subdivision of the loop into two topologically independent loops by insulators inhibits interdomain interactions. The underlying cause of this effect is a suppression of crossings in the contact map at intermediate distances. Thus our model simultaneously accounts for regulatory interactions at a distance and the insulator activity of boundary elements. This unified model of the regulatory roles of chromatin loops makes several testable predictions that could be confronted with in vitro experiments, as well as genomic chromatin conformation capture and fluorescent microscopic approaches.

Locus dependence in epigenetic chromatin silencing.

Current biological models of epigenetic switches built on chromatin modifications lead to strong constraints on the repertoire of dynamic behaviors for the system. We use the structure of the bifurcation diagram of the underlying dynamical system to explain the existing single cell data in silencing by the SIR system in yeast.

Inheritance of epigenetic chromatin silencing.

Maintenance of alternative chromatin states through cell divisions pose some fundamental constraints on the dynamics of histone modifications. In this paper, we study the systems biology of epigenetic inheritance by defining and analyzing general classes of mathematical models. We discuss how the number of modification states involved plays an essential role in the stability of epigenetic states. In addition, DNA duplication and the consequent dilution of marked histones act as a large perturbation for a stable state of histone modifications. The requirement that this large perturbation falls into the basin of attraction of the original state sometimes leads to additional constraints on effective models. Two such models, inspired by two different biological systems, are compared in their fulfilling the requirements of multistability and of recovery after DNA duplication. We conclude that in the presence of multiple histone modifications that characterize alternative epigenetic stable states, these requirements are more easily fulfilled.

Nematic elastomers: from a microscopic model to macroscopic elasticity theory.

A Landau theory is constructed for the gelation transition in cross-linked polymer systems possessing spontaneous nematic ordering, based on symmetry principles and the concept of an order parameter for the amorphous solid state. This theory is substantiated with help of a simple microscopic model of cross-linked dimers. Minimization of the Landau free energy in the presence of nematic order yields the neoclassical theory of the elasticity of nematic elastomers and, in the isotropic limit, the classical theory of isotropic elasticity. These phenomenological theories of elasticity are thereby derived from a microscopic model, and it is furthermore demonstrated that they are universal mean-field descriptions of the elasticity for all chemical gels and vulcanized media.

Anomalous diffusion and stretched exponentials in heterogeneous glass-forming liquids: low-temperature behavior.

We propose a model of a heterogeneous glass-forming liquid and compute the low-temperature behavior of a tagged molecule moving within it. This model exhibits stretched-exponential decay of the wave-number-dependent, self-intermediate scattering function in the limit of long times. At temperatures close to the glass transition, where the heterogeneities are much larger in extent than the molecular spacing, the time dependence of the scattering function crosses over from stretched-exponential decay with an index b=1/2 at large wave numbers to normal, diffusive behavior with b=1 at small wave numbers. There is a clear separation between early-stage, cage-breaking beta relaxation and late-stage alpha relaxation. The spatial representation of the scattering function exhibits an anomalously broad exponential (non-Gaussian) tail for sufficiently large values of the molecular displacement at all finite times.

Scaling of entropic shear rigidity.

The scaling of shear modulus near the gelation-vulcanization transition is explored heuristically and analytically. It is found that in a dense melt the effective chains of the infinite cluster have sizes that scale sublinearly with their contour length. Consequently, each chain contributes k(B)T to the rigidity, which leads to a shear-modulus exponent dnu. In contrast, in phantom elastic networks the scaling is linear in the contour length, yielding an exponent identical to that of the random resistor network conductivity, as predicted by de Gennes. For nondense systems, the exponent should cross over to dnu when the percolation correlation length is much larger than the density-fluctuation length.