A parallelized, automated platform enabling individual or sequential ChIP of histone marks and transcription factors

Multiplexed spatial mapping of chromatin features, transcriptome, and proteins in tissues

The phenotypic and functional states of a cell are modulated by a complex interactive molecular hierarchy of multiple omics layers, involving the genome, epigenome, transcriptome, proteome, and metabolome. Spatial omics approaches have enabled the capture of information from different molecular layers directly in the tissue context. However, current technologies are limited to map one to two modalities at the same time, providing an incomplete representation of cellular identity. Such data is inadequate to fully understand complex biological systems and their underlying regulatory mechanisms. Here we present spatial-Mux-seq, a multi-modal spatial technology that allows simultaneous profiling of five different modalities, including genome-wide profiles of two histone modifications and open chromatin, whole transcriptome, and a panel of proteins at tissue scale and cellular level in a spatially resolved manner. We applied this technology to generate multi-modal tissue maps in mouse embryos and mouse brains, which discriminated more cell types and states than unimodal data. We investigated the spatiotemporal relationship between histone modifications, chromatin accessibility, gene and protein expression in neuron differentiation revealing the relationship between tissue organization, function, and gene regulatory networks. We were able to identify a radial glia spatial niche and revealed spatially changing gradient of epigenetic signals in this region. Moreover, we revealed previously unappreciated involvement of repressive histone marks in the mouse hippocampus. Collectively, the spatial multi-omics approach heralds a new era for characterizing tissue and cellular heterogeneity that single modality studies alone could not reveal.

Metabolic stress-induced long ncRNA transcription governs the formation of meiotic DNA breaks in the fission yeast fbp1 gene

Meiotic recombination is a pivotal process that ensures faithful chromosome segregation and contributes to the generation of genetic diversity in offspring, which is initiated by the formation of double-strand breaks (DSBs). The distribution of meiotic DSBs is not uniform and is clustered at hotspots, which can be affected by environmental conditions. Here, we show that non-coding RNA (ncRNA) transcription creates meiotic DSBs through local chromatin remodeling in the fission yeast fbp1 gene. The fbp1 gene is activated upon glucose starvation stress, in which a cascade of ncRNA-transcription in the fbp1 upstream region converts the chromatin configuration into an open structure, leading to the subsequent binding of transcription factors. We examined the distribution of meiotic DSBs around the fbp1 upstream region in the presence and absence of glucose and observed several new DSBs after chromatin conversion under glucose starvation conditions. Moreover, these DSBs disappeared when cis-elements required for ncRNA transcription were mutated. These results indicate that ncRNA transcription creates meiotic DSBs in response to stress conditions in the fbp1 upstream region. This study addressed part of a long-standing unresolved mechanism underlying meiotic recombination plasticity in response to environmental fluctuation.

Patterns of transcription factor binding and epigenome at promoters allow interpretable predictability of multiple functions of non-coding and coding genes

Understanding the biological roles of all genes only through experimental methods is challenging. A computational approach with reliable interpretability is needed to infer the function of genes, particularly for non-coding RNAs. We have analyzed genomic features that are present across both coding and non-coding genes like transcription factor (TF) and cofactor ChIP-seq (823), histone modifications ChIP-seq (n = 621), cap analysis gene expression (CAGE) tags (n = 255), and DNase hypersensitivity profiles (n = 255) to predict ontology-based functions of genes. Our approach for gene function prediction was reliable (>90% balanced accuracy) for 486 gene-sets. PubMed abstract mining and CRISPR screens supported the inferred association of genes with biological functions, for which our method had high accuracy. Further analysis revealed that TF-binding patterns at promoters have high predictive strength for multiple functions. TF-binding patterns at the promoter add an unexplored dimension of explainable regulatory aspects of genes and their functions. Therefore, we performed a comprehensive analysis for the functional-specificity of TF-binding patterns at promoters and used them for clustering functions to reveal many latent groups of gene-sets involved in common major cellular processes. We also showed how our approach could be used to infer the functions of non-coding genes using the CRISPR screens of coding genes, which were validated using a long non-coding RNA CRISPR screen. Thus our results demonstrated the generality of our approach by using gene-sets from CRISPR screens. Overall, our approach opens an avenue for predicting the involvement of non-coding genes in various functions.

A microfluidics-enabled automated workflow of sample preparation for MS-based immunopeptidomics

Mass spectrometry (MS)-based immunopeptidomics is an attractive antigen discovery method with growing clinical implications. However, the current experimental approach to extract HLA-restricted peptides requires a bulky sample source, which remains a challenge for obtaining clinical specimens. We present an innovative workflow that requires a low sample volume, which streamlines the immunoaffinity purification (IP) and C18 peptide cleanup on a single microfluidics platform with automated liquid handling and minimal sample transfers, resulting in higher assay sensitivity. We also demonstrate how the state-of-the-art data-independent acquisition (DIA) method further enhances the depth of tandem MS spectra-based peptide sequencing. Consequently, over 4,000 and 5,000 HLA-I-restricted peptides were identified from as few as 0.2 million RA957 cells and a melanoma tissue of merely 5 mg, respectively. We also identified multiple immunogenic tumor-associated antigens and hundreds of peptides derived from non-canonical protein sources. This workflow represents a powerful tool for identifying the immunopeptidome of sparse samples.

Natural and Experimental Rewiring of Gene Regulatory Regions

The successful development and ongoing functioning of complex organisms depend on the faithful execution of the genetic code. A critical step in this process is the correct spatial and temporal expression of genes. The highly orchestrated transcription of genes is controlled primarily by cis-regulatory elements: promoters, enhancers, and insulators. The medical importance of this key biological process can be seen by the frequency with which mutations and inherited variants that alter cis-regulatory elements lead to monogenic and complex diseases and cancer. Here, we provide an overview of the methods available to characterize and perturb gene regulatory circuits. We then highlight mechanisms through which regulatory rewiring contributes to disease, and conclude with a perspective on how our understanding of gene regulation can be used to improve human health.

Single-cell joint detection of chromatin occupancy and transcriptome enables higher-dimensional epigenomic reconstructions

Deciphering mechanisms in cell-fate decisions requires single-cell holistic reconstructions of multidimensional epigenomic states in transcriptional regulation. Here we develop CoTECH, a combinatorial barcoding method allowing high-throughput single-cell joint detection of chromatin occupancy and transcriptome. We used CoTECH to examine bivalent histone marks (H3K4me3 and H3K27me3) with transcription from naive to primed mouse embryonic stem cells. We also derived concurrent bivalent marks in pseudosingle cells using transcriptome as an anchor for resolving pseudotemporal bivalency trajectories and disentangling a context-specific interplay between H3K4me3/H3K27me3 and transcription level. Next, we revealed the regulatory basis of endothelial-to-hematopoietic transition in two waves of hematopoietic cells and distinctive enhancer-gene-linking schemes guiding hemogenic endothelial cell emergence, indicating a unique epigenetic control of transcriptional regulation for hematopoietic stem cell priming. CoTECH provides an efficient framework for single-cell coassay of chromatin occupancy and transcription, thus enabling higher-dimensional epigenomic reconstructions.