Defining super-enhancers by highly ranked histone H4 multi-acetylation levels identifies transcription factors associated with glioblastoma stem-like properties.

BACKGROUND: Super-enhancers (SEs), which activate genes involved in cell-type specificity, have mainly been defined as genomic regions with top-ranked enrichment(s) of histone H3 with acetylated K27 (H3K27ac) and/or transcription coactivator(s) including a bromodomain and extra-terminal domain (BET) family protein, BRD4. However, BRD4 preferentially binds to multi-acetylated histone H4, typically with acetylated K5 and K8 (H4K5acK8ac), leading us to hypothesize that SEs should be defined by high H4K5acK8ac enrichment at least as well as by that of H3K27ac. RESULTS: Here, we conducted genome-wide profiling of H4K5acK8ac and H3K27ac, BRD4 binding, and the transcriptome by using a BET inhibitor, JQ1, in three human glial cell lines. When SEs were defined as having the top ranks for H4K5acK8ac or H3K27ac signal, 43% of H4K5acK8ac-ranked SEs were distinct from H3K27ac-ranked SEs in a glioblastoma stem-like cell (GSC) line. CRISPR-Cas9-mediated deletion of the H4K5acK8ac-preferred SEs associated with MYCN and NFIC decreased the stem-like properties in GSCs. CONCLUSIONS: Collectively, our data highlights H4K5acK8ac's utility for identifying genes regulating cell-type specificity.

The ENCODE Imputation Challenge: a critical assessment of methods for cross-cell type imputation of epigenomic profiles.

A promising alternative to comprehensively performing genomics experiments is to, instead, perform a subset of experiments and use computational methods to impute the remainder. However, identifying the best imputation methods and what measures meaningfully evaluate performance are open questions. We address these questions by comprehensively analyzing 23 methods from the ENCODE Imputation Challenge. We find that imputation evaluations are challenging and confounded by distributional shifts from differences in data collection and processing over time, the amount of available data, and redundancy among performance measures. Our analyses suggest simple steps for overcoming these issues and promising directions for more robust research.

Antisense-oligonucleotide-mediated perturbation of long non-coding RNA reveals functional features in stem cells and across cell types.

Within the scope of the FANTOM6 consortium, we perform a large-scale knockdown of 200 long non-coding RNAs (lncRNAs) in human induced pluripotent stem cells (iPSCs) and systematically characterize their roles in self-renewal and pluripotency. We find 36 lncRNAs (18%) exhibiting cell growth inhibition. From the knockdown of 123 lncRNAs with transcriptome profiling, 36 lncRNAs (29.3%) show molecular phenotypes. Integrating the molecular phenotypes with chromatin-interaction assays further reveals cis- and trans-interacting partners as potential primary targets. Additionally, cell-type enrichment analysis identifies lncRNAs associated with pluripotency, while the knockdown of LINC02595, CATG00000090305.1, and RP11-148B6.2 modulates colony formation of iPSCs. We compare our results with previously published fibroblasts phenotyping data and find that 2.9% of the lncRNAs exhibit a consistent cell growth phenotype, whereas we observe 58.3% agreement in molecular phenotypes. This highlights that molecular phenotyping is more comprehensive in revealing affected pathways.

SCAFE: a software suite for analysis of transcribed cis-regulatory elements in single cells.

MOTIVATION: Cell type-specific activities of cis-regulatory elements (CRE) are central to understanding gene regulation and disease predisposition. Single-cell RNA 5'end sequencing (sc-end5-seq) captures the transcription start sites (TSS) which can be used as a proxy to measure the activity of transcribed CREs (tCREs). However, a substantial fraction of TSS identified from sc-end5-seq data may not be genuine due to various artifacts, hindering the use of sc-end5-seq for de novo discovery of tCREs. RESULTS: We developed SCAFE-Single-Cell Analysis of Five-prime Ends-a software suite that processes sc-end5-seq data to de novo identify TSS clusters based on multiple logistic regression. It annotates tCREs based on the identified TSS clusters and generates a tCRE-by-cell count matrix for downstream analyses. The software suite consists of a set of flexible tools that could either be run independently or as pre-configured workflows. AVAILABILITY AND IMPLEMENTATION: SCAFE is implemented in Perl and R. The source code and documentation are freely available for download under the MIT License from https://github.com/chung-lab/SCAFE. Docker images are available from https://hub.docker.com/r/cchon/scafe. The submitted software version and test data are archived at https://doi.org/10.5281/zenodo.7023163 and https://doi.org/10.5281/zenodo.7024060, respectively. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Functional annotation of human long noncoding RNAs via molecular phenotyping.

Long noncoding RNAs (lncRNAs) constitute the majority of transcripts in the mammalian genomes, and yet, their functions remain largely unknown. As part of the FANTOM6 project, we systematically knocked down the expression of 285 lncRNAs in human dermal fibroblasts and quantified cellular growth, morphological changes, and transcriptomic responses using Capped Analysis of Gene Expression (CAGE). Antisense oligonucleotides targeting the same lncRNAs exhibited global concordance, and the molecular phenotype, measured by CAGE, recapitulated the observed cellular phenotypes while providing additional insights on the affected genes and pathways. Here, we disseminate the largest-to-date lncRNA knockdown data set with molecular phenotyping (over 1000 CAGE deep-sequencing libraries) for further exploration and highlight functional roles for ZNF213-AS1 and lnc-KHDC3L-2.

OVOL2 induces mesenchymal-to-epithelial transition in fibroblasts and enhances cell-state reprogramming towards epithelial lineages.

Mesenchymal-to-epithelial transition (MET) is an important step in cell reprogramming from fibroblasts (a cell type frequently used for this purpose) to various epithelial cell types. However, the mechanism underlying MET induction in fibroblasts remains to be understood. The present study aimed to identify the transcription factors (TFs) that efficiently induce MET in dermal fibroblasts. OVOL2 was identified as a potent inducer of key epithelial genes, and OVOL2 cooperatively enhanced MET induced by HNF1A, TP63, and KLF4, which are known reprogramming TFs to epithelial lineages. In TP63/KLF4-induced keratinocyte-like cell-state reprogramming, OVOL2 greatly facilitated the activation of epithelial and keratinocyte-specific genes. This was accompanied by enhanced changes in chromatin accessibility across the genome. Mechanistically, motif enrichment analysis revealed that the target loci of KLF4 and TP63 become accessible upon induction of TFs, whereas the OVOL2 target loci become inaccessible. This indicates that KLF4 and TP63 positively regulate keratinocyte-associated genes whereas OVOL2 suppresses fibroblast-associated genes. The exogenous expression of OVOL2 therefore disrupts fibroblast lineage identity and facilitates fibroblast cell reprogramming into epithelial lineages cooperatively with tissue-specific reprogramming factors. Identification of OVOL2 as an MET inducer and an epithelial reprogramming enhancer in fibroblasts provides new insights into cellular reprogramming improvement for future applications.

Chromatin accessibility landscape of articular knee cartilage reveals aberrant enhancer regulation in osteoarthritis.

Osteoarthritis (OA) is a common joint disorder with increasing impact in an aging society. While genetic and transcriptomic analyses have revealed some genes and non-coding loci associated to OA, the pathogenesis remains incompletely understood. Chromatin profiling, which provides insight into gene regulation, has not been reported in OA mainly due to technical difficulties. Here, we employed Assay for Transposase-Accessible Chromatin with high throughput sequencing (ATAC-seq) to map the accessible chromatin landscape in articular knee cartilage of OA patients. We identified 109,215 accessible chromatin regions for cartilages, of which 71% were annotated as enhancers. By overlaying them with genetic and DNA methylation data, we have determined potential OA-relevant enhancers and their putative target genes. Furthermore, through integration with RNA-seq data, we characterized genes that are altered both at epigenomic and transcriptomic levels in OA. These genes are enriched in pathways regulating ossification and mesenchymal stem cell (MSC) differentiation. Consistently, the differentially accessible regions in OA are enriched for MSC-specific enhancers and motifs of transcription factor families involved in osteoblast differentiation. In conclusion, we demonstrate how direct chromatin profiling of clinical tissues can provide comprehensive epigenetic information for a disease and suggest candidate genes and enhancers of translational potential.

Fabrication of birefringent nanocylinders for single-molecule force and torque measurement.

Optically anisotropic subwavelength scale dielectric particles have been shown to enable studies of the mechanical properties of bio-molecules via optical trapping and manipulation. However, techniques emphasized to date for fabrication of such particles generally suffer from limited uniformity and control over particle dimensions, or low throughput and high cost. Here, an approach for rapid, low-cost, fabrication of large quantities of birefringent quartz nanocylinders with dimensions optimized for optical torque wrench experiments is described. For a typical process, 10(8) or more quartz cylinders with diameters of 500 nm and heights of 800 nm, with uniformity of ±5% in each dimension, can be fabricated over ∼10 cm(2) areas, for binding to a single bio-molecule, and harvested for use in optical trapping experiments. Use of these structures to measure extensional and torsional dynamics of single DNA molecules is demonstrated with measured forces and torques shown to be in very good agreement with previously reported results.

Single-molecule measurements of the CCR5 mRNA unfolding pathways.

Secondary or tertiary structure in an mRNA, such as a pseudoknot, can create a physical barrier that requires the ribosome to generate additional force to translocate. The presence of such a barrier can dramatically increase the probability that the ribosome will shift into an alternate reading frame, in which a different set of codons is recognized. The detailed biophysical mechanism by which frameshifting is induced remains unknown. Here we employ optical trapping techniques to investigate the structure of a -1 programmed ribosomal frameshift (-1 PRF) sequence element located in the CCR5 mRNA, which encodes a coreceptor for HIV-1 and is, to our knowledge, the first known human -1 PRF signal of nonviral origin. We begin by presenting a set of computationally predicted structures that include pseudoknots. We then employ what we believe to be new analytical techniques for measuring the effective free energy landscapes of biomolecules. We find that the -1 PRF element manifests several distinct unfolding pathways when subject to end-to-end force, one of which is consistent with a proposed pseudoknot conformation, and another of which we have identified as a folding intermediate. The dynamic ensemble of conformations that CCR5 mRNA exhibits in the single-molecule experiments may be a significant feature of the frameshifting mechanism.

Effect of handle length and microsphere size on transition kinetics in single-molecule experiments.

When subject to constant tension, a DNA or RNA hairpin will typically make abrupt transitions between the open and closed state. Although the transition kinetics are an intrinsic property of the molecule, the transition rates measured in single-molecule experiments can be influenced by the configuration of the measurement system. We investigate the transition kinetics for a DNA hairpin held under constant force by an optical trap as a function of microsphere size and double-stranded DNA handle length. We find the apparent transition lifetime cannot be expressed as a function of the drag coefficient of the microsphere alone or as a function of time scales relevant to the optical trap. The apparent transition lifetime is found to be a linear function of the factor ß(eff)·α(handle), where ß(eff) is the effective drag coefficient of the microsphere near the surface and α(handle) is the stiffness of the DNA tether. The results provide insight into the perturbation to the hairpin transition kinetics due to experimental configuration and guidance for designing single-molecule experiments which determine the intrinsic molecular kinetics.

Single-molecule study of G-quadruplex disruption using dynamic force spectroscopy.

Guanine-rich sequences in nucleic acids can fold into G quadruplexes, in which four guanines on a single strand combine to form G-tetrad planes stabilized by metallic ions. Sequence motifs which are predicted to form a G quadruplex are found throughout the genome and are believed to regulate a variety of biological processes. Detailed knowledge of the kinetics of G-quadruplex folding and unfolding would provide critical insight into these processes. To probe its structural stability, we used optical tweezers to disrupt single molecules of a single-stranded DNA G4 quadruplex. Dynamic force spectroscopy was employed, in which the distribution of rupture forces was measured for different loading rates and used to infer the nature of the transition state barrier for unfolding of the structure. The distance and height of the energy barriers were extracted for two observed conformations. The energy barrier was found to be close to the folded conformation, resulting in a high disruption force despite the relatively low energy barrier height.

FAST MONTE CARLO SIMULATION METHODS FOR BIOLOGICAL REACTION-DIFFUSION SYSTEMS IN SOLUTION AND ON SURFACES.

Many important physiological processes operate at time and space scales far beyond those accessible to atom-realistic simulations, and yet discrete stochastic rather than continuum methods may best represent finite numbers of molecules interacting in complex cellular spaces. We describe and validate new tools and algorithms developed for a new version of the MCell simulation program (MCell3), which supports generalized Monte Carlo modeling of diffusion and chemical reaction in solution, on surfaces representing membranes, and combinations thereof. A new syntax for describing the spatial directionality of surface reactions is introduced, along with optimizations and algorithms that can substantially reduce computational costs (e.g., event scheduling, variable time and space steps). Examples for simple reactions in simple spaces are validated by comparison to analytic solutions. Thus we show how spatially realistic Monte Carlo simulations of biological systems can be far more cost-effective than often is assumed, and provide a level of accuracy and insight beyond that of continuum methods.