chromVAR: inferring transcription-factor-associated accessibility from single-cell epigenomic data

Transposable Elements Are Co-opted as Oncogenic Regulatory Elements by Lineage-Specific Transcription Factors in Prostate Cancer

Transposable elements hold regulatory functions that impact cell fate determination by controlling gene expression. However, little is known about the transcriptional machinery engaged at transposable elements in pluripotent and mature versus oncogenic cell states. Through positional analysis over repetitive DNA sequences of H3K27ac chromatin immunoprecipitation sequencing data from 32 normal cell states, we report pluripotent/stem and mature cell state-specific "regulatory transposable elements." Pluripotent/stem elements are binding sites for pluripotency factors (e.g., NANOG, SOX2, OCT4). Mature cell elements are docking sites for lineage-specific transcription factors, including AR and FOXA1 in prostate epithelium. Expanding the analysis to prostate tumors, we identify a subset of regulatory transposable elements shared with pluripotent/stem cells, including Tigger3a. Using chromatin editing technology, we show how such elements promote prostate cancer growth by regulating AR transcriptional activity. Collectively, our results suggest that oncogenesis arises from lineage-specific transcription factors hijacking pluripotent/stem cell regulatory transposable elements.

SIGNIFICANCE

We show that oncogenesis relies on co-opting transposable elements from pluripotent stem cells as regulatory elements altering the recruitment of lineage-specific transcription factors. We further discover how co-option is dependent on active chromatin states with important implications for developing treatment options against drivers of oncogenesis across the repetitive DNA. This article is featured in Selected Articles from This Issue, p. 2293.

Retrospective analysis of enhancer activity and transcriptome history

Cell state changes in development and disease are controlled by gene regulatory networks, the dynamics of which are difficult to track in real time. In this study, we used an inducible DCM-RNA polymerase subunit b fusion protein which labels active genes and enhancers with a bacterial methylation mark that does not affect gene transcription and is propagated in S-phase. This DCM-RNA polymerase fusion protein enables transcribed genes and active enhancers to be tagged and then examined at later stages of development or differentiation. We apply this DCM-time machine (DCM-TM) technology to study intestinal homeostasis, revealing rapid and coordinated activation of enhancers and nearby genes during enterocyte differentiation. We provide new insights in absorptive-secretory lineage decision-making in intestinal stem cell (ISC) differentiation and show that ISCs retain a unique chromatin landscape required to maintain ISC identity and delineate future expression of differentiation-associated genes. DCM-TM has wide applicability in tracking cell states, providing new insights in the regulatory networks underlying cell state changes.

Gene regulatory network inference in the era of single-cell multi-omics

The interplay between chromatin, transcription factors and genes generates complex regulatory circuits that can be represented as gene regulatory networks (GRNs). The study of GRNs is useful to understand how cellular identity is established, maintained and disrupted in disease. GRNs can be inferred from experimental data - historically, bulk omics data - and/or from the literature. The advent of single-cell multi-omics technologies has led to the development of novel computational methods that leverage genomic, transcriptomic and chromatin accessibility information to infer GRNs at an unprecedented resolution. Here, we review the key principles of inferring GRNs that encompass transcription factor-gene interactions from transcriptomics and chromatin accessibility data. We focus on the comparison and classification of methods that use single-cell multimodal data. We highlight challenges in GRN inference, in particular with respect to benchmarking, and potential further developments using additional data modalities.

Multi-modal skin atlas identifies a multicellular immune-stromal community associated with altered cornification and specific T cell expansion in atopic dermatitis

In healthy skin, a cutaneous immune system maintains the balance between tolerance towards innocuous environmental antigens and immune responses against pathological agents. In atopic dermatitis (AD), barrier and immune dysfunction result in chronic tissue inflammation. Our understanding of the skin tissue ecosystem in AD remains incomplete with regard to the hallmarks of pathological barrier formation, and cellular state and clonal composition of disease-promoting cells. Here, we generated a multi-modal cell census of 310,691 cells spanning 86 cell subsets from whole skin tissue of 19 adult individuals, including non-lesional and lesional skin from 11 AD patients, and integrated it with 396,321 cells from four studies into a comprehensive human skin cell atlas in health and disease. Reconstruction of human keratinocyte differentiation from basal to cornified layers revealed a disrupted cornification trajectory in AD. This disrupted epithelial differentiation was associated with signals from a unique immune and stromal multicellular community comprised of MMP12 + dendritic cells (DCs), mature migratory DCs, cycling ILCs, NK cells, inflammatory CCL19 + IL4I1 + fibroblasts, and clonally expanded IL13 + IL22 + IL26 + T cells with overlapping type 2 and type 17 characteristics. Cell subsets within this immune and stromal multicellular community were connected by multiple inter-cellular positive feedback loops predicted to impact community assembly and maintenance. AD GWAS gene expression was enriched both in disrupted cornified keratinocytes and in cell subsets from the lesional immune and stromal multicellular community including IL13 + IL22 + IL26 + T cells and ILCs, suggesting that epithelial or immune dysfunction in the context of the observed cellular communication network can initiate and then converge towards AD. Our work highlights specific, disease-associated cell subsets and interactions as potential targets in progression and resolution of chronic inflammation.

RORγt+ c-Maf+ Vγ4+ γδ T cells are generated in the adult thymus but do not reach the periphery

T cell receptor (TCR) Vγ4-expressing γδ T cells comprise interferon γ (IFNγ)- and interleukin-17 (IL-17)-producing effector subsets, with a preference for IL-17 effector fate decisions during early ontogeny. The existence of adult-thymus-derived IL-17+ T cells (γδ17) remains controversial. Here, we use a mouse model in which T cells are generated exclusively in the adult thymus and employ single-cell chromatin state analysis to study their development. We identify adult-thymus-derived Vγ4 T cells that have all the molecular programs to become IL-17 producers. However, they have reduced IL-17 production capabilities and rarely reach the periphery. Moreover, this study provides high-resolution profiles of Vγ4 T cells in the adult thymus and lymph nodes and identifies Zeb1 as a potential γδ17 cell regulator. Together, this study provides valuable insights into the developmental traits of Vγ4 T cells during adulthood and supports the idea of age-specific signals required for thymic export and/or peripheral maturation of γδ17 cells.

Ikaros is a principal regulator of Aire+ mTEC homeostasis, thymic mimetic cell diversity, and central tolerance

Mutations in the gene encoding the zinc-finger transcription factor Ikaros (IKZF1) are found in patients with immunodeficiency, leukemia, and autoimmunity. Although Ikaros has a well-established function in modulating gene expression programs important for hematopoietic development, its role in other cell types is less well defined. Here, we uncover functions for Ikaros in thymic epithelial lineage development in mice and show that Ikzf1 expression in medullary thymic epithelial cells (mTECs) is required for both autoimmune regulator-positive (Aire+) mTEC development and tissue-specific antigen (TSA) gene expression. Accordingly, TEC-specific deletion of Ikzf1 in mice results in a profound decrease in Aire+ mTECs, a global loss of TSA gene expression, and the development of autoimmunity. Moreover, Ikaros shapes thymic mimetic cell diversity, and its deletion results in a marked expansion of thymic tuft cells and muscle-like mTECs and a loss of other Aire-dependent mimetic populations. Single-cell analysis reveals that Ikaros modulates core transcriptional programs in TECs that correlate with the observed cellular changes. Our findings highlight a previously undescribed role for Ikaros in regulating epithelial lineage development and function and suggest that failed thymic central tolerance could contribute to the autoimmunity seen in humans with IKZF1 mutations.

A Germline Point Mutation in the MYC-FBW7 Phosphodegron Initiates Hematopoietic Malignancies

Oncogenic activation of MYC in cancers predominantly involves increased transcription rather than coding region mutations. However, MYC-dependent lymphomas frequently contain point mutations in the MYC phospho-degron, including at threonine-58 (T58), where phosphorylation permits binding by the FBW7 ubiquitin ligase triggering MYC degradation. To understand how T58 phosphorylation functions in normal cell physiology, we introduced an alanine mutation at T58 (T58A) into the endogenous c-Myc locus in the mouse germline. While MYC-T58A mice develop normally, lymphomas and myeloid leukemias emerge in ~60% of adult homozygous T58A mice. We find that primitive hematopoietic progenitor cells from MYC-T58A mice exhibit aberrant self-renewal normally associated with hematopoietic stem cells (HSCs) and upregulate a subset of Myc target genes important in maintaining stem/progenitor cell balance. Genomic occupancy by MYC-T58A was increased at all promoters, compared to WT MYC, while genes differentially expressed in a T58A-dependent manner were significantly more proximal to MYC-bound enhancers. MYC-T58A lymphocyte progenitors exhibited metabolic alterations and decreased activation of inflammatory and apoptotic pathways. Our data demonstrate that a single point mutation in Myc is sufficient to produce a profound gain of function in multipotential hematopoietic progenitors associated with self-renewal and initiation of lymphomas and leukemias.

Chromatin accessibility profiling of targeted cell populations with laser capture microdissection coupled to ATAC-seq

Spatially resolved omics technologies reveal context-dependent cellular regulatory networks in tissues of interest. Beyond transcriptome analysis, information on epigenetic traits and chromatin accessibility can provide further insights on gene regulation in health and disease. Nevertheless, compared to the enormous advancements in spatial transcriptomics technologies, the field of spatial epigenomics is much younger and still underexplored. In this study, we report laser capture microdissection coupled to ATAC-seq (LCM-ATAC-seq) applied to fresh frozen samples for the spatial characterization of chromatin accessibility. We first demonstrate the efficient use of LCM coupled to in situ tagmentation and evaluate its performance as a function of cell number, microdissected areas, and tissue type. Further, we demonstrate its use for the targeted chromatin accessibility analysis of discrete contiguous or scattered cell populations in tissues via single-nuclei capture based on immunostaining for specific cellular markers.

Comprehensive single cell aging atlas of mammary tissues reveals shared epigenomic and transcriptomic signatures of aging and cancer

Aging is the greatest risk factor for breast cancer; however, how age-related cellular and molecular events impact cancer initiation is unknown. We investigate how aging rewires transcriptomic and epigenomic programs of mouse mammary glands at single cell resolution, yielding a comprehensive resource for aging and cancer biology. Aged epithelial cells exhibit epigenetic and transcriptional changes in metabolic, pro-inflammatory, or cancer-associated genes. Aged stromal cells downregulate fibroblast marker genes and upregulate markers of senescence and cancer-associated fibroblasts. Among immune cells, distinct T cell subsets (Gzmk+, memory CD4+, γδ) and M2-like macrophages expand with age. Spatial transcriptomics reveal co-localization of aged immune and epithelial cells in situ. Lastly, transcriptional signatures of aging mammary cells are found in human breast tumors, suggesting mechanistic links between aging and cancer. Together, these data uncover that epithelial, immune, and stromal cells shift in proportions and cell identity, potentially impacting cell plasticity, aged microenvironment, and neoplasia risk.

Transcription factor stoichiometry, motif affinity and syntax regulate single-cell chromatin dynamics during fibroblast reprogramming to pluripotency

Ectopic expression of OCT4, SOX2, KLF4 and MYC (OSKM) transforms differentiated cells into induced pluripotent stem cells. To refine our mechanistic understanding of reprogramming, especially during the early stages, we profiled chromatin accessibility and gene expression at single-cell resolution across a densely sampled time course of human fibroblast reprogramming. Using neural networks that map DNA sequence to ATAC-seq profiles at base-resolution, we annotated cell-state-specific predictive transcription factor (TF) motif syntax in regulatory elements, inferred affinity- and concentration-dependent dynamics of Tn5-bias corrected TF footprints, linked peaks to putative target genes, and elucidated rewiring of TF-to-gene cis-regulatory networks. Our models reveal that early in reprogramming, OSK, at supraphysiological concentrations, rapidly open transient regulatory elements by occupying non-canonical low-affinity binding sites. As OSK concentration falls, the accessibility of these transient elements decays as a function of motif affinity. We find that these OSK-dependent transient elements sequester the somatic TF AP-1. This redistribution is strongly associated with the silencing of fibroblast-specific genes within individual nuclei. Together, our integrated single-cell resource and models reveal insights into the cis-regulatory code of reprogramming at unprecedented resolution, connect TF stoichiometry and motif syntax to diversification of cell fate trajectories, and provide new perspectives on the dynamics and role of transient regulatory elements in somatic silencing.

Reprogramming of cis-regulatory networks during skeletal muscle atrophy in male mice

A comprehensive atlas of cis-regulatory elements and their dynamic activity is necessary to understand the transcriptional basis of cellular structure maintenance, metabolism, and responses to the environment. Here we show, using matched single-nucleus chromatin accessibility and RNA-sequencing from juvenile male C57BL6 mice, an atlas of accessible chromatin regions in both normal and denervated skeletal muscles. We identified cell-type-specific cis-regulatory networks, highlighting the dynamic regulatory circuits mediating transitions between myonuclear types. Through comparison of normal and perturbed muscle, we delineated the reprogramming of cis-regulatory networks in response to denervation, described the interplay of promoters/enhancers and target genes. We further unveil a hierarchical structure of transcription factors that delineate a regulatory network in atrophic muscle, identifying ELK4 as a key atrophy-related transcription factor that instigates muscle atrophy through TGF-β1 regulation. This study furnishes a rich genomic resource, essential for decoding the regulatory dynamics of skeletal muscle in both physiological and pathological states.

Single-cell chromatin accessibility and transcriptomic characterization of Behcet's disease

Behect's disease is a chronic vasculitis characterized by complex multi-organ immune aberrations. However, a comprehensive understanding of the gene-regulatory profile of peripheral autoimmunity and the diverse immune responses across distinct cell types in Behcet's disease (BD) is still lacking. Here, we present a multi-omic single-cell study of 424,817 cells in BD patients and non-BD individuals. This study maps chromatin accessibility and gene expression in the same biological samples, unraveling vast cellular heterogeneity. We identify widespread cell-type-specific, disease-associated active and pro-inflammatory immunity in both transcript and epigenomic aspects. Notably, integrative multi-omic analysis reveals putative TF regulators that might contribute to chromatin accessibility and gene expression in BD. Moreover, we predicted gene-regulatory networks within nominated TF activators, including AP-1, NF-kB, and ETS transcript factor families, which may regulate cellular interaction and govern inflammation. Our study illustrates the epigenetic and transcriptional landscape in BD peripheral blood and expands understanding of potential epigenomic immunopathology in this disease.

Using single-cell chromatin accessibility sequencing to characterize CD4+ T cells from murine tissues

The Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq) is a cutting-edge technology that enables researchers to assess genome-wide chromatin accessibility and to characterize cell type specific gene-regulatory programs. Recent technological progress allows for using this technology also on the single-cell level. In this article, we describe the whole value chain from the isolation of T cells from murine tissues to a complete bioinformatic analysis workflow. We start with methods for isolating scATAC-seq-ready CD4+ T cells from murine tissues such as visceral adipose tissue, skin, colon, and secondary lymphoid tissues such as the spleen. We describe the preparation of nuclei and quality control parameters during library preparation. Based on publicly available sequencing data that was generated using these protocols, we describe a step-by-step bioinformatic analysis pipeline for data pre-processing and downstream analysis. Our analysis workflow will follow the R-based bioinformatics framework ArchR, which is currently well established for scATAC-seq datasets. All in all, this work serves as a one-stop shop for generating and analyzing chromatin accessibility landscapes in T cells.

Inflammatory stress-mediated chromatin changes underlie dysfunction in endothelial cells

Inflammatory stresses underlie endothelial dysfunction and contribute to the development of chronic cardiovascular disorders such as atherosclerosis and vascular fibrosis. The initial transcriptional response of endothelial cells to pro-inflammatory cytokines such as TNF-alpha is well established. However, very few studies uncover the effects of inflammatory stresses on chromatin architecture. We used integrative analysis of ATAC-seq and RNA-seq data to investigate chromatin alterations in human endothelial cells in response to TNF-alpha and febrile-range heat stress exposure. Multi-omics data analysis suggests a correlation between the transcription of stress-related genes and endothelial dysfunction drivers with chromatin regions exhibiting differential accessibility. Moreover, microscopy identified the dynamics in the nuclear organization, specifically, the changes in a subset of heterochromatic nucleoli-associated chromatin domains, the centromeres. Upon inflammatory stress exposure, the centromeres decreased association with nucleoli in a p38-dependent manner and increased the number of transcripts from pericentromeric regions. Overall, we provide two lines of evidence that suggest chromatin alterations in vascular endothelial cells during inflammatory stresses.

Synaptic Activity Causes Minute-scale Changes in BAF Complex Composition and Function

Genes encoding subunits of the SWI/SNF or BAF ATP-dependent chromatin remodeling complex are among the most enriched for deleterious de novo mutations in intellectual disabilities and autism spectrum disorder, but the causative molecular pathways are not fully known 1,2 . Synaptic activity in neurons is critical for learning and memory and proper neural development 3 . Neural activity prompts calcium influx and transcription within minutes, facilitated in the nucleus by various transcription factors (TFs) and chromatin modifiers 4 . While BAF is required for activity-dependent developmental processes such as dendritic outgrowth 5-7 , the immediate molecular consequences of neural activity on BAF complexes and their functions are unknown. Here we mapped minute-scale biochemical consequences of neural activity, modeled by membrane depolarization of embryonic mouse primary cortical neurons, on BAF complexes. We used acute chemical perturbations of BAF ATPase activity and kinase signaling to define the activity-dependent effects on BAF complexes and activity-dependent BAF functions. Our studies found that BAF complexes change in subunit composition and are selectively phosphorylated within 10 minutes of depolarization. Increased levels of the core PBAF subunit Baf200/ Arid2 , uniquely containing an RFX-like DNA-binding domain, are concurrent with ATPase-dependent opening of chromatin at RFX/X-box motifs. Changes in BAF composition and phosphorylation lead to the regulation of chromatin accessibility for critical neurogenesis TFs. These biochemical effects are a convergent phenomenon downstream of multiple growth factor signaling pathways in mouse neurons and fibroblasts suggesting that BAF integrates signaling information from the membrane. In support of such a membrane-to-nucleus signaling cascade, we also identified a BAF-interacting kinase, Dclk2, whose inhibition attenuates BAF phosphorylation selectively. Our findings support a direct role of BAF complexes in responding to synaptic activity to regulate TF binding and transcription.

A comparative atlas of single-cell chromatin accessibility in the human brain

Recent advances in single-cell transcriptomics have illuminated the diverse neuronal and glial cell types within the human brain. However, the regulatory programs governing cell identity and function remain unclear. Using a single-nucleus assay for transposase-accessible chromatin using sequencing (snATAC-seq), we explored open chromatin landscapes across 1.1 million cells in 42 brain regions from three adults. Integrating this data unveiled 107 distinct cell types and their specific utilization of 544,735 candidate cis-regulatory DNA elements (cCREs) in the human genome. Nearly a third of the cCREs demonstrated conservation and chromatin accessibility in the mouse brain cells. We reveal strong links between specific brain cell types and neuropsychiatric disorders including schizophrenia, bipolar disorder, Alzheimer's disease (AD), and major depression, and have developed deep learning models to predict the regulatory roles of noncoding risk variants in these disorders.

Multi-omics analysis of mucosal and systemic immunity to SARS-CoV-2 after birth

The dynamics of immunity to infection in infants remain obscure. Here, we used a multi-omics approach to perform a longitudinal analysis of immunity to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in infants and young children by analyzing blood samples and weekly nasal swabs collected before, during, and after infection with Omicron and non-Omicron variants. Infection stimulated robust antibody titers that, unlike in adults, showed no sign of decay for up to 300 days. Infants mounted a robust mucosal immune response characterized by inflammatory cytokines, interferon (IFN) α, and T helper (Th) 17 and neutrophil markers (interleukin [IL]-17, IL-8, and CXCL1). The immune response in blood was characterized by upregulation of activation markers on innate cells, no inflammatory cytokines, but several chemokines and IFNα. The latter correlated with viral load and expression of interferon-stimulated genes (ISGs) in myeloid cells measured by single-cell multi-omics. Together, these data provide a snapshot of immunity to infection during the initial weeks and months of life.

Analysis of potential biomarkers for diabetic kidney disease based on single-cell RNA-sequencing integrated with a single-cell sequencing assay for transposase-accessible chromatin

Diabetic kidney disease (DKD) is a renal microvascular disease caused by hyperglycemia that involves metabolic remodeling, oxidative stress, inflammation, and other factors. The mechanism is complex and not fully unraveled. We performed an integrated single-cell sequencing assay for transposase-accessible chromatin (scATAC-seq) and single-cell RNA-sequencing (scRNA-seq) analyses of kidneys from db/db and db/m mice to identify differential open chromatin regions and gene expression, particularly in genes related to proximal tubular reabsorption and secretion. We identified 9,776 differentially expressed genes (DEGs) and 884 cell type-specific transcription factors (TFs) across 15 cell types. Glucose and lipid transporters, and TFs related to the circadian rhythm in the proximal tubules had significantly higher expression in db/db mice than in db/m mice (P<0.01). Crosstalk between podocytes and tubular cells in the proximal tubules was enhanced, and renal inflammation, oxidative stress, and fibrosis pathways were activated in db/db mice. Western blotting and immunohistochemical staining results showed that Wfdc2 expression in the urine and kidneys of DKD patients was higher than that in non-diabetic kidney disease (NDKD) controls. The revealed landscape of chromatin accessibility and transcriptional profiles in db/db mice provide insights into the pathological mechanism of DKD.

Integrated single-cell multiome analysis reveals muscle fiber-type gene regulatory circuitry modulated by endurance exercise

Endurance exercise is an important health modifier. We studied cell-type specific adaptations of human skeletal muscle to acute endurance exercise using single-nucleus (sn) multiome sequencing in human vastus lateralis samples collected before and 3.5 hours after 40 min exercise at 70% VO2max in four subjects, as well as in matched time of day samples from two supine resting circadian controls. High quality same-cell RNA-seq and ATAC-seq data were obtained from 37,154 nuclei comprising 14 cell types. Among muscle fiber types, both shared and fiber-type specific regulatory programs were identified. Single-cell circuit analysis identified distinct adaptations in fast, slow and intermediate fibers as well as LUM-expressing FAP cells, involving a total of 328 transcription factors (TFs) acting at altered accessibility sites regulating 2,025 genes. These data and circuit mapping provide single-cell insight into the processes underlying tissue and metabolic remodeling responses to exercise.

Systematic identification of gene combinations to target in innate immune cells to enhance T cell activation

Genetic engineering of immune cells has opened new avenues for improving their functionality but it remains a challenge to pinpoint which genes or combination of genes are the most beneficial to target. Here, we conduct High Multiplicity of Perturbations and Cellular Indexing of Transcriptomes and Epitopes (HMPCITE-seq) to find combinations of genes whose joint targeting improves antigen-presenting cell activity and enhances their ability to activate T cells. Specifically, we perform two genome-wide CRISPR screens in bone marrow dendritic cells and identify negative regulators of CD86, that participate in the co-stimulation programs, including Chd4, Stat5b, Egr2, Med12, and positive regulators of PD-L1, that participate in the co-inhibitory programs, including Sptlc2, Nckap1l, and Pi4kb. To identify the genetic interactions between top-ranked genes and find superior combinations to target, we perform high-order Perturb-Seq experiments and we show that targeting both Cebpb and Med12 results in a better phenotype compared to the single perturbations or other combinations of perturbations.

A transcriptional and regulatory map of mouse somite maturation

The mammalian body plan is shaped by rhythmic segmentation of mesoderm into somites, which are transient embryonic structures that form down each side of the neural tube. We have analyzed the genome-wide transcriptional and chromatin dynamics occurring within nascent somites, from early inception of somitogenesis to the latest stages of body plan establishment. We created matched gene expression and open chromatin maps for the three leading pairs of somites at six time points during mouse embryonic development. We show that the rate of somite differentiation accelerates as development progresses. We identified a conserved maturation program followed by all somites, but somites from more developed embryos concomitantly switch on differentiation programs from derivative cell lineages soon after segmentation. Integrated analysis of the somitic transcriptional and chromatin activities identified opposing regulatory modules controlling the onset of differentiation. Our results provide a powerful, high-resolution view of the molecular genetics underlying somitic development in mammals.

Integrome signatures of lentiviral gene therapy for SCID-X1 patients

Lentiviral vector (LV)-based gene therapy holds promise for a broad range of diseases. Analyzing more than 280,000 vector integration sites (VISs) in 273 samples from 10 patients with X-linked severe combined immunodeficiency (SCID-X1), we discovered shared LV integrome signatures in 9 of 10 patients in relation to the genomics, epigenomics, and 3D structure of the human genome. VISs were enriched in the nuclear subcompartment A1 and integrated into super-enhancers close to nuclear pore complexes. These signatures were validated in T cells transduced with an LV encoding a CD19-specific chimeric antigen receptor. Intriguingly, the one patient whose VISs deviated from the identified integrome signatures had a distinct clinical course. Comparison of LV and gamma retrovirus integromes regarding their 3D genome signatures identified differences that might explain the lower risk of insertional mutagenesis in LV-based gene therapy. Our findings suggest that LV integrome signatures, shaped by common features such as genome organization, may affect the efficacy of LV-based cellular therapies.

Chromatin Accessibility Landscape of Human Triple-negative Breast Cancer Cell Lines Reveals Variation by Patient Donor Ancestry

African American (AA) women have an excessive risk of developing triple-negative breast cancer (TNBC). We employed Assay for Transposase-Accessible Chromatin using sequencing to characterize differences in chromatin accessibility between nine commonly used TNBC cell lines derived from patients of European and African ancestry. Principal component and chromosome mapping analyses of accessibility peaks with the most variance revealed separation of chromatin profiles by patient group. Motif enrichment and footprinting analyses of disparate open chromatin regions revealed differences in transcription factor activity, identifying 79 with ancestry-associated binding patterns (FDR < 0.01). AA TNBC cell lines exhibited increased accessibility for 62 transcription factors associated with epithelial-to-mesenchymal transition, cancer stemness/chemotherapeutic resistance, proliferation, and aberrant p53 regulation, as well as KAISO, which has been previously linked to aggressive tumor characteristics in AA patients with cancer. Differential Assay for Transposase-Accessible Chromatin signal analysis identified 1,596 genes located within promoters of differentially open chromatin regions in AA-derived TNBC, identifying DNA methyltransferase 1 as the top upregulated gene associated with African ancestry. Pathway analyses with these genes revealed enrichment in several pathways, including hypoxia. Culturing cells under hypoxia showed ancestry-specific stress responses that led to the identification of a core set of AA-associated transcription factors, which included members of the Kruppel-like factor and Sp subfamilies, as well as KAISO, and identified ZDHHC1, a gene previously implicated in immunity and STING activation, as the top upregulated AA-specific gene under hypoxia. Together, these data reveal a differential chromatin landscape in TNBC associated with donor ancestry. The open chromatin structure of AA TNBC may contribute to a more lethal disease.

SIGNIFICANCE

We identify an ancestry-associated open chromatin landscape and related transcription factors that may contribute to aggressive TNBC in AA women. Furthermore, this study advocates for the inclusion of diversely sourced cell lines in experimental in vitro studies to advance health equity at all levels of scientific research.

Hnf4 activates mimetic-cell enhancers to recapitulate gut and liver development within the thymus

Mimetic cells are medullary thymic epithelial cells (mTECs) that mimic extra-thymic cell types to tolerize T cells to self-antigens. Here, we dissected the biology of entero-hepato mTECs, mimetic cells expressing gut- and liver-associated transcripts. Entero-hepato mTECs conserved their thymic identity yet accessed wide swaths of enterocyte chromatin and transcriptional programs via the transcription factors Hnf4α and Hnf4γ. Deletion of Hnf4α and Hnf4γ in TECs ablated entero-hepato mTECs and downregulated numerous gut- and liver-associated transcripts, with a primary contribution from Hnf4γ. Loss of Hnf4 impaired enhancer activation and CTCF redistribution in mTECs but did not impact Polycomb-mediated repression or promoter-proximal histone marks. By single-cell RNA sequencing, Hnf4 loss produced three distinct effects on mimetic cell state, fate, and accumulation. Serendipitously, a requirement for Hnf4 in microfold mTECs was discovered, which exposed a requirement for Hnf4γ in gut microfold cells and the IgA response. Study of Hnf4 in entero-hepato mTECs thus revealed mechanisms of gene control in the thymus and periphery alike.

RUNX1 loss renders hematopoietic and leukemic cells dependent on IL-3 and sensitive to JAK inhibition

Disease-initiating mutations in the transcription factor RUNX1 occur as germline and somatic events that cause leukemias with particularly poor prognosis. However, the role of RUNX1 in leukemogenesis is not fully understood, and effective therapies for RUNX1-mutant leukemias remain elusive. Here, we used primary patient samples and a RUNX1-KO model in primary human hematopoietic cells to investigate how RUNX1 loss contributes to leukemic progression and to identify targetable vulnerabilities. Surprisingly, we found that RUNX1 loss decreased proliferative capacity and stem cell function. However, RUNX1-deficient cells selectively upregulated the IL-3 receptor. Exposure to IL-3, but not other JAK/STAT cytokines, rescued RUNX1-KO proliferative and competitive defects. Further, we demonstrated that RUNX1 loss repressed JAK/STAT signaling and rendered RUNX1-deficient cells sensitive to JAK inhibitors. Our study identifies a dependency of RUNX1-mutant leukemias on IL-3/JAK/STAT signaling, which may enable targeting of these aggressive blood cancers with existing agents.

PD-1 instructs a tumor-suppressive metabolic program that restricts glycolysis and restrains AP-1 activity in T cell lymphoma

The PDCD1-encoded immune checkpoint receptor PD-1 is a key tumor suppressor in T cells that is recurrently inactivated in T cell non-Hodgkin lymphomas (T-NHLs). The highest frequencies of PDCD1 deletions are detected in advanced disease, predicting inferior prognosis. However, the tumor-suppressive mechanisms of PD-1 signaling remain unknown. Here, using tractable mouse models for T-NHL and primary patient samples, we demonstrate that PD-1 signaling suppresses T cell malignancy by restricting glycolytic energy and acetyl coenzyme A (CoA) production. In addition, PD-1 inactivation enforces ATP citrate lyase (ACLY) activity, which generates extramitochondrial acetyl-CoA for histone acetylation to enable hyperactivity of activating protein 1 (AP-1) transcription factors. Conversely, pharmacological ACLY inhibition impedes aberrant AP-1 signaling in PD-1-deficient T-NHLs and is toxic to these cancers. Our data uncover genotype-specific vulnerabilities in PDCD1-mutated T-NHL and identify PD-1 as regulator of AP-1 activity.

Transcriptional and epigenetic decoding of the microglial aging process

As important immune cells, microglia undergo a series of alterations during aging that increase the susceptibility to brain dysfunctions. However, the longitudinal characteristics of microglia remain poorly understood. In this study, we mapped the transcriptional and epigenetic profiles of microglia from 3- to 24-month-old mice. We first discovered unexpected sex differences and identified age-dependent microglia (ADEM) genes during the aging process. We then compared the features of aging and reactivity in female microglia at single-cell resolution and epigenetic level. To dissect functions of aged microglia excluding the influence from other aged brain cells, we established an accelerated microglial turnover model without directly affecting other brain cells. By this model, we achieved aged-like microglia in non-aged brains and confirmed that aged-like microglia per se contribute to cognitive decline. Collectively, our work provides a comprehensive resource for decoding the aging process of microglia, shedding light on how microglia maintain brain functions.

Thymic mimetic cells function beyond self-tolerance

Development of immunocompetent T cells in the thymus is required for effective defence against all types of pathogens, including viruses, bacteria and fungi. To this end, T cells undergo a very strict educational program in the thymus, during which both non-functional and self-reactive T cell clones are eliminated by means of positive and negative selection1.Thymic epithelial cells (TECs) have an indispensable role in these processes, and previous studies have shown the notable heterogeneity of these cells2-7. Here, using multiomic analysis, we provide further insights into the functional and developmental diversity of TECs in mice, and reveal a detailed atlas of the TEC compartment according to cell transcriptional states and chromatin landscapes. Our analysis highlights unconventional TEC subsets that are similar to functionally well-defined parenchymal populations, including endocrine cells, microfold cells and myocytes. By focusing on the endocrine and microfold TEC populations, we show that endocrine TECs require Insm1 for their development and are crucial to maintaining thymus cellularity in a ghrelin-dependent manner; by contrast, microfold TECs require Spib for their development and are essential for the generation of thymic IgA+ plasma cells. Collectively, our study reveals that medullary TECs have the potential to differentiate into various types of molecularly distinct and functionally defined cells, which not only contribute to the induction of central tolerance, but also regulate the homeostasis of other thymus-resident populations.

Pluripotent stem cell-derived model of the post-implantation human embryo

The human embryo undergoes morphogenetic transformations following implantation into the uterus, but our knowledge of this crucial stage is limited by the inability to observe the embryo in vivo. Models of the embryo derived from stem cells are important tools for interrogating developmental events and tissue-tissue crosstalk during these stages1. Here we establish a model of the human post-implantation embryo, a human embryoid, comprising embryonic and extraembryonic tissues. We combine two types of extraembryonic-like cell generated by overexpression of transcription factors with wild-type embryonic stem cells and promote their self-organization into structures that mimic several aspects of the post-implantation human embryo. These self-organized aggregates contain a pluripotent epiblast-like domain surrounded by extraembryonic-like tissues. Our functional studies demonstrate that the epiblast-like domain robustly differentiates into amnion, extraembryonic mesenchyme and primordial germ cell-like cells in response to bone morphogenetic protein cues. In addition, we identify an inhibitory role for SOX17 in the specification of anterior hypoblast-like cells2. Modulation of the subpopulations in the hypoblast-like compartment demonstrates that extraembryonic-like cells influence epiblast-like domain differentiation, highlighting functional tissue-tissue crosstalk. In conclusion, we present a modular, tractable, integrated3 model of the human embryo that will enable us to probe key questions of human post-implantation development, a critical window during which substantial numbers of pregnancies fail.

Gene-regulation modules in nonalcoholic fatty liver disease revealed by single-nucleus ATAC-seq

We investigated the progression of nonalcoholic fatty liver disease from fatty liver to steatohepatitis using single-nucleus and bulk ATAC-seq on the livers of rats fed a high-fat diet (HFD). Rats fed HFD for 4 wk developed fatty liver, and those fed HFD for 8 wk further progressed to steatohepatitis. We observed an increase in the proportion of inflammatory macrophages, consistent with the pathological progression. Utilizing machine learning, we divided global gene regulation into modules, wherein transcription factors within a module could regulate genes within the same module, reaffirming known regulatory relationships between transcription factors and biological processes. We identified core genes-central to co-expression and protein-protein interaction-for the biological processes discovered. Notably, a large part of the core genes overlapped with genes previously implicated in nonalcoholic fatty liver disease. Single-nucleus ATAC-seq, combined with data-driven statistical analysis, offers insight into in vivo global gene regulation as a combination of modules and assists in identifying core genes of relevant biological processes.

Tracking cell-type-specific temporal dynamics in human and mouse brains

Progenitor cells are critical in preserving organismal homeostasis, yet their diversity and dynamics in the aged brain remain underexplored. We introduced TrackerSci, a single-cell genomic method that combines newborn cell labeling and combinatorial indexing to characterize the transcriptome and chromatin landscape of proliferating progenitor cells in vivo. Using TrackerSci, we investigated the dynamics of newborn cells in mouse brains across various ages and in a mouse model of Alzheimer's disease. Our dataset revealed diverse progenitor cell types in the brain and their epigenetic signatures. We further quantified aging-associated shifts in cell-type-specific proliferation and differentiation and deciphered the associated molecular programs. Extending our study to the progenitor cells in the aged human brain, we identified conserved genetic signatures across species and pinpointed region-specific cellular dynamics, such as the reduced oligodendrogenesis in the cerebellum. We anticipate that TrackerSci will be broadly applicable to unveil cell-type-specific temporal dynamics in diverse systems.

Osteogenic Differentiation Potential of Mesenchymal Stem Cells Using Single Cell Multiomic Analysis

Mesenchymal stem cells (MSC) are multipotent stem cells that can differentiate into multiple cell types, including osteoblasts, chondrocytes, and adipocytes. Osteoblast differentiation is reduced during osteoporosis development, resulting in reduced bone formation. Further, MSC isolated from different donors possess distinct osteogenic capacity. In this study, we used single-cell multiomic analysis to profile the transcriptome and epigenome of MSC from four healthy donors. Data were obtained from ~1300 to 1600 cells for each donor. These cells were clustered into four groups, indicating that MSC from different donors have distinct chromatin accessible regulatory elements for regulating gene expression. To investigate the mechanism by which MSC undergo osteogenic differentiation, we used the chromatin accessibility data from the single-cell multiome data to identify individual-specific enhancer-promoter pairs and evaluated the expression levels and activities of the transcriptional regulators. The MSC from four donors showed distinct differentiation potential into osteoblasts. MSC of donor 1 showed the largest average motif activities, indicating that MSC from donor 1 was most likely to differentiate into osteoblasts. The results of our validation experiments were consistent with the bioinformatics prediction. We also tested the enrichment of genome-wide association study (GWAS) signals of several musculoskeletal disease traits in the patient-specific chromatin accessible regions identified in the single-cell multiome data, including osteoporosis, osteopenia, and osteoarthritis. We found that osteoarthritis-associated variants were only enriched in the regions identified from donor 4. In contrast, osteoporosis and osteopenia variants were enriched in regions from donor 1 and least enriched in donor 4. Since osteoporosis and osteopenia are related to the density of bone cells, the enrichment of variants from these traits should be correlated with the osteogenic potential of MSC. In summary, this study provides large-scale data to link regulatory elements with their target genes to study the regulatory relationships during the differentiation of mesenchymal stem cells and provide a deeper insight into the gene regulatory mechanism.

Single-nuclei multiomic analyses identify human cardiac lymphatic endothelial cells associated with coronary arteries in the epicardium

Cardiac lymphatic vessels play important roles in fluid homeostasis, inflammation, disease, and regeneration of the heart. The developing cardiac lymphatics in human fetal hearts are closely associated with coronary arteries, similar to those in zebrafish hearts. We identify a population of cardiac lymphatic endothelial cells (LECs) that reside in the epicardium. Single-nuclei multiomic analysis of the human fetal heart reveals the plasticity and heterogeneity of the cardiac endothelium. Furthermore, we find that VEGFC is highly expressed in arterial endothelial cells and epicardium-derived cells, providing a molecular basis for the arterial association of cardiac lymphatic development. Using a cell-type-specific integrative analysis, we identify a population of cardiac lymphatic endothelial cells marked by the PROX1 and the lymphangiocrine RELN and enriched in binding motifs of erythroblast transformation specific (ETS) variant (ETV) transcription factors. We report the in vivo molecular characterization of human cardiac lymphatics and provide a valuable resource to understand fetal heart development.

A branching model of lineage differentiation underpinning the neurogenic potential of enteric glia

Glial cells have been proposed as a source of neural progenitors, but the mechanisms underpinning the neurogenic potential of adult glia are not known. Using single cell transcriptomic profiling, we show that enteric glial cells represent a cell state attained by autonomic neural crest cells as they transition along a linear differentiation trajectory that allows them to retain neurogenic potential while acquiring mature glial functions. Key neurogenic loci in early enteric nervous system progenitors remain in open chromatin configuration in mature enteric glia, thus facilitating neuronal differentiation under appropriate conditions. Molecular profiling and gene targeting of enteric glial cells in a cell culture model of enteric neurogenesis and a gut injury model demonstrate that neuronal differentiation of glia is driven by transcriptional programs employed in vivo by early progenitors. Our work provides mechanistic insight into the regulatory landscape underpinning the development of intestinal neural circuits and generates a platform for advancing glial cells as therapeutic agents for the treatment of neural deficits.

Evaluation of the determinants for improved pluripotency induction and maintenance by engineered SOX17

An engineered SOX17 variant with point mutations within its DNA binding domain termed SOX17FNV is a more potent pluripotency inducer than SOX2, yet the underlying mechanism remains unclear. Although wild-type SOX17 was incapable of inducing pluripotency, SOX17FNV outperformed SOX2 in mouse and human pluripotency reprogramming. In embryonic stem cells, SOX17FNV could replace SOX2 to maintain pluripotency despite considerable sequence differences and upregulated genes expressed in cleavage-stage embryos. Mechanistically, SOX17FNV co-bound OCT4 more cooperatively than SOX2 in the context of the canonical SoxOct DNA element. SOX2, SOX17, and SOX17FNV were all able to bind nucleosome core particles in vitro, which is a prerequisite for pioneer transcription factors. Experiments using purified proteins and in cellular contexts showed that SOX17 variants phase-separated more efficiently than SOX2, suggesting an enhanced ability to self-organise. Systematic deletion analyses showed that the N-terminus of SOX17FNV was dispensable for its reprogramming activity. However, the C-terminus encodes essential domains indicating multivalent interactions that drive transactivation and reprogramming. We defined a minimal SOX17FNV (miniSOX) that can support reprogramming with high activity, reducing the payload of reprogramming cassettes. This study uncovers the mechanisms behind SOX17FNV-induced pluripotency and establishes engineered SOX factors as powerful cell engineering tools.

Toward the functional interpretation of somatic structural variations: bulk- and single-cell approaches

Structural variants (SVs) are genomic rearrangements that can take many different forms such as copy number alterations, inversions and translocations. During cell development and aging, somatic SVs accumulate in the genome with potentially neutral, deleterious or pathological effects. Generation of somatic SVs is a key mutational process in cancer development and progression. Despite their importance, the detection of somatic SVs is challenging, making them less studied than somatic single-nucleotide variants. In this review, we summarize recent advances in whole-genome sequencing (WGS)-based approaches for detecting somatic SVs at the tissue and single-cell levels and discuss their advantages and limitations. First, we describe the state-of-the-art computational algorithms for somatic SV calling using bulk WGS data and compare the performance of somatic SV detectors in the presence or absence of a matched-normal control. We then discuss the unique features of cutting-edge single-cell-based techniques for analyzing somatic SVs. The advantages and disadvantages of bulk and single-cell approaches are highlighted, along with a discussion of their sensitivity to copy-neutral SVs, usefulness for functional inferences and experimental and computational costs. Finally, computational approaches for linking somatic SVs to their functional readouts, such as those obtained from single-cell transcriptome and epigenome analyses, are illustrated, with a discussion of the promise of these approaches in health and diseases.

Epigenetic Induction of Smooth Muscle Cell Phenotypic Alterations in Aortic Aneurysms and Dissections

BACKGROUND

Smooth muscle cell (SMC) phenotypic switching has been increasingly detected in aortic aneurysm and dissection (AAD) tissues. However, the diverse SMC phenotypes in AAD tissues and the mechanisms driving SMC phenotypic alterations remain to be identified.

METHODS

We examined the transcriptomic and epigenomic dynamics of aortic SMC phenotypic changes in mice with angiotensin II-induced AAD by using single-cell RNA sequencing and single-cell sequencing assay for transposase-accessible chromatin. SMC phenotypic alteration in aortas from patients with ascending thoracic AAD was examined by using single-cell RNA sequencing analysis.

RESULTS

Single-cell RNA sequencing analysis revealed that aortic stress induced the transition of SMCs from a primary contractile phenotype to proliferative, extracellular matrix-producing, and inflammatory phenotypes. Lineage tracing showed the complete transformation of SMCs to fibroblasts and macrophages. Single-cell sequencing assay for transposase-accessible chromatin analysis indicated that these phenotypic alterations were controlled by chromatin remodeling marked by the reduced chromatin accessibility of contractile genes and the induced chromatin accessibility of genes involved in proliferation, extracellular matrix, and inflammation. IRF3 (interferon regulatory factor 3), a proinflammatory transcription factor activated by cytosolic DNA, was identified as a key driver of the transition of aortic SMCs from a contractile phenotype to an inflammatory phenotype. In cultured SMCs, cytosolic DNA signaled through its sensor STING (stimulator of interferon genes)-TBK1 (tank-binding kinase 1) to activate IRF3, which bound and recruited EZH2 (enhancer of zeste homolog 2) to contractile genes to induce repressive H3K27me3 modification and gene suppression. In contrast, double-stranded DNA-STING-IRF3 signaling induced inflammatory gene expression in SMCs. In Sting-/- mice, the aortic stress-induced transition of SMCs into an inflammatory phenotype was prevented, and SMC populations were preserved. Finally, profound SMC phenotypic alterations toward diverse directions were detected in human ascending thoracic AAD tissues.

CONCLUSIONS

Our study reveals the dynamic epigenetic induction of SMC phenotypic alterations in AAD. DNA damage and cytosolic leakage drive SMCs from a contractile phenotype to an inflammatory phenotype.

Common and divergent gene regulatory networks control injury-induced and developmental neurogenesis in zebrafish retina

Following acute retinal damage, zebrafish possess the ability to regenerate all neuronal subtypes. This regeneration requires Müller glia (MG) to reprogram and divide asymmetrically to produce a multipotent Müller glia-derived neuronal progenitor cell (MGPC). This raises three key questions. First, does loss of different retinal cell subtypes induce unique MG regeneration responses? Second, do MG reprogram to a developmental retinal progenitor cell state? And finally, to what extent does regeneration recapitulate retinal development? We examined these questions by performing single-nuclear and single-cell RNA-Seq and ATAC-Seq in both developing and regenerating retinas. While MG reprogram to a state similar to late-stage retinal progenitors in developing retinas, there are transcriptional differences between reprogrammed MG/MGPCs and late progenitors, as well as reprogrammed MG in outer and inner retinal damage models. Validation of candidate genes confirmed that loss of different subtypes induces differences in transcription factor gene expression and regeneration outcomes. This work identifies major differences between gene regulatory networks activated following the selective loss of different subtypes of retina neurons, as well as between retinal regeneration and development.

Divergent single cell transcriptome and epigenome alterations in ALS and FTD patients with C9orf72 mutation

A repeat expansion in the C9orf72 (C9) gene is the most common genetic cause of amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). Here we investigate single nucleus transcriptomics (snRNA-seq) and epigenomics (snATAC-seq) in postmortem motor and frontal cortices from C9-ALS, C9-FTD, and control donors. C9-ALS donors present pervasive alterations of gene expression with concordant changes in chromatin accessibility and histone modifications. The greatest alterations occur in upper and deep layer excitatory neurons, as well as in astrocytes. In neurons, the changes imply an increase in proteostasis, metabolism, and protein expression pathways, alongside a decrease in neuronal function. In astrocytes, the alterations suggest activation and structural remodeling. Conversely, C9-FTD donors have fewer high-quality neuronal nuclei in the frontal cortex and numerous gene expression changes in glial cells. These findings highlight a context-dependent molecular disruption in C9-ALS and C9-FTD, indicating unique effects across cell types, brain regions, and diseases.

Dynamic single-cell regulomes characterize human peripheral blood innate lymphoid cell subpopulations

Innate lymphoid cells (ILCs) are plastic immune cells divided into 3 main subsets, characterized by distinct phenotypic and functional profiles. Using single cell approaches, heightened heterogeneity of mouse ILCs has been appreciated, imprinted by tissue signals that shape their transcriptome and epigenome. Intra-subset diversity has also been observed in human ILCs. However, combined transcriptomic and epigenetic analyses of single ILCs in humans are lacking. Here, we show high transcriptional and epigenetic heterogeneity among human circulating ILCs in healthy individuals. We describe phenotypically distinct subclusters and diverse chromatin accessibility within main ILC populations, compatible with differentially poised states. We validate the use of this healthy donor-based analysis as resource dataset to help inferring ILC changes occurring in disease conditions. Overall, our work provides insights in the complex human ILC biology. We anticipate it to facilitate hypothesis-driven studies in patients, without the need to perform single cell OMICs using precious patients' material.

Human microglia maturation is underpinned by specific gene regulatory networks

Microglia phenotypes are highly regulated by the brain environment, but the transcriptional networks that specify the maturation of human microglia are poorly understood. Here, we characterized stage-specific transcriptomes and epigenetic landscapes of fetal and postnatal human microglia and acquired corresponding data in induced pluripotent stem cell (iPSC)-derived microglia, in cerebral organoids, and following engraftment into humanized mice. Parallel development of computational approaches that considered transcription factor (TF) co-occurrence and enhancer activity allowed prediction of shared and state-specific gene regulatory networks associated with fetal and postnatal microglia. Additionally, many features of the human fetal-to-postnatal transition were recapitulated in a time-dependent manner following the engraftment of iPSC cells into humanized mice. These data and accompanying computational approaches will facilitate further efforts to elucidate mechanisms by which human microglia acquire stage- and disease-specific phenotypes.

Integration of spatial and single-cell data across modalities with weakly linked features

Although single-cell and spatial sequencing methods enable simultaneous measurement of more than one biological modality, no technology can capture all modalities within the same cell. For current data integration methods, the feasibility of cross-modal integration relies on the existence of highly correlated, a priori 'linked' features. We describe matching X-modality via fuzzy smoothed embedding (MaxFuse), a cross-modal data integration method that, through iterative coembedding, data smoothing and cell matching, uses all information in each modality to obtain high-quality integration even when features are weakly linked. MaxFuse is modality-agnostic and demonstrates high robustness and accuracy in the weak linkage scenario, achieving 20~70% relative improvement over existing methods under key evaluation metrics on benchmarking datasets. A prototypical example of weak linkage is the integration of spatial proteomic data with single-cell sequencing data. On two example analyses of this type, MaxFuse enabled the spatial consolidation of proteomic, transcriptomic and epigenomic information at single-cell resolution on the same tissue section.

Single cell multiomic analysis reveals diabetes-associated β-cell heterogeneity driven by HNF1A

Broad heterogeneity in pancreatic β-cell function and morphology has been widely reported. However, determining which components of this cellular heterogeneity serve a diabetes-relevant function remains challenging. Here, we integrate single-cell transcriptome, single-nuclei chromatin accessibility, and cell-type specific 3D genome profiles from human islets and identify Type II Diabetes (T2D)-associated β-cell heterogeneity at both transcriptomic and epigenomic levels. We develop a computational method to explicitly dissect the intra-donor and inter-donor heterogeneity between single β-cells, which reflect distinct mechanisms of T2D pathogenesis. Integrative transcriptomic and epigenomic analysis identifies HNF1A as a principal driver of intra-donor heterogeneity between β-cells from the same donors; HNF1A expression is also reduced in β-cells from T2D donors. Interestingly, HNF1A activity in single β-cells is significantly associated with lower Na+ currents and we nominate a HNF1A target, FXYD2, as the primary mitigator. Our study demonstrates the value of investigating disease-associated single-cell heterogeneity and provides new insights into the pathogenesis of T2D.

Deciphering early human pancreas development at the single-cell level

Understanding pancreas development can provide clues for better treatments of pancreatic diseases. However, the molecular heterogeneity and developmental trajectory of the early human pancreas are poorly explored. Here, we performed large-scale single-cell RNA sequencing and single-cell assay for transposase accessible chromatin sequencing of human embryonic pancreas tissue obtained from first-trimester embryos. We unraveled the molecular heterogeneity, developmental trajectories and regulatory networks of the major cell types. The results reveal that dorsal pancreatic multipotent cells in humans exhibit different gene expression patterns than ventral multipotent cells. Pancreato-biliary progenitors that generate ventral multipotent cells in humans were identified. Notch and MAPK signals from mesenchymal cells regulate the differentiation of multipotent cells into trunk and duct cells. Notably, we identified endocrine progenitor subclusters with different differentiation potentials. Although the developmental trajectories are largely conserved between humans and mice, some distinct gene expression patterns have also been identified. Overall, we provide a comprehensive landscape of early human pancreas development to understand its lineage transitions and molecular complexity.

Single-cell multi-omics sequencing reveals the immunological disturbance underlying STAT3-V637M Hyper-IgE syndrome

Hyper-IgE syndrome (HIES) is a primary immunodeficiency characterized by, among others, the excessive production of IgE and repetitive bacterial/fungal infections. Mutations in STAT3, a transcription factor that orchestrates immune responses, may cause HIES, but the underlying mechanisms are not fully understood. Here, we used multi-omic approaches to comprehensively decipher the immune disturbance in a male HIES patient harboring STAT3-V637M. In his peripheral blood mononuclear cell (PBMC) we found significant clonal expansion of CD8 T cells (with increased CD8 subunits expression, potentially enhancing responsiveness to MHC I molecules), but not in his CD4 T cells and B cells. Although his B cells exhibited a higher potential in producing immunoglobulin, elevated SPIC binding might bias the products toward IgE isotype. Immune checkpoint inhibitors, including CTLA4, LAG3, were overexpressed in his PBMC-CD4 T cells, accompanied by reduced CD28 and IL6ST (gp130) expression. In his CD4 T cells, integrative analyses predicted upstream transcription factors (including ETV6, KLF13, and RORA) for LAG3, IL6ST, and CD28, respectively. The down-regulation of phagocytosis and nitric oxide synthesis-related genes in his PBMC-monocytes seem to be the culprit of his disseminated bacterial/fungal infection. Counterintuitively, in his PBMC we predicted increased STAT3 binding in both naïve and mature CD4 compartments, although this was not observed in most of his PBMC. In his bronchoalveolar lavage fluid (BALF), we found two macrophage subtypes with anti-bacterial properties, which were identified by CXCL8/S100A8/S100A9, or SOD2, respectively. Together, we described how the immune cell landscape was disturbed in STAT3-V637M HIES, providing a resource for further studies.

Chromatin accessibility dynamics in colorectal cancer liver metastasis: Uncovering the liver tropism at single cell resolution

Tumor metastasis causes over 90% of cancer related death and no currently available therapies target it. However, there is limited understanding regarding the epigenetic regulation of genes during this complex process. Here by integrating single-cell ATAC-seq (scATAC-seq), single-cell RNA-seq (scRNA-seq), microarray, bulk RNA-seq, immunohistochemistry (IHC) staining, as well as proteomics datasets from paired primary and liver metastatic colorectal cancer (CRC) patient-derived xenograft (PDX) model and patients, we discovered that liver metastatic CRC cells lose their colon-specific chromatin accessible sites yet gain liver-specific ones. Importantly, we observed elevated accessibility of HNF4A, a liver-specific transcription factor, in liver metastatic CRC cells. Subsequently, we performed clustering analysis of liver metastatic CRC cells together with cells involved in liver development, revealing significant heterogeneity among the liver metastatic CRC cells. Over 50% of the liver metastatic CRC cells exhibited characteristics similar to those of erythroid progenitors and hepatocytes, showing increased expression of genes involved in oxidative phosphorylation and glycolysis. Moreover, our discovery further revealed that the MHC and IFN response genes in these cells exhibit moderate epigenetic activity, which is significantly associated with the low objective response rates in checkpoint blockade immunotherapy. Our findings uncovered the critical roles of HNF4A and the cell populations within liver metastatic CRC cells might serve as crucial therapeutic targets for addressing liver metastasis and improving the immunotherapy response in patients with CRC.

Differences in syncytia formation by SARS-CoV-2 variants modify host chromatin accessibility and cellular senescence via TP53

COVID-19 remains a significant public health threat due to the ability of SARS-CoV-2 variants to evade the immune system and cause breakthrough infections. Although pathogenic coronaviruses such as SARS-CoV-2 and MERS-CoV lead to severe respiratory infections, how these viruses affect the chromatin proteomic composition upon infection remains largely uncharacterized. Here we used our recently developed integrative DNA And Protein Tagging (iDAPT) methodology to identify changes in host chromatin accessibility states and chromatin proteomic composition upon infection with pathogenic coronaviruses. SARS-CoV-2 infection induces TP53 stabilization on chromatin, which contributes to its host cytopathic effect. We mapped this TP53 stabilization to the SARS-CoV-2 spike and its propensity to form syncytia, a consequence of cell-cell fusion. Differences in SARS-CoV-2 spike variant-induced syncytia formation modify chromatin accessibility, cellular senescence, and inflammatory cytokine release via TP53. Our findings suggest that differences in syncytia formation alter senescence-associated inflammation, which varies among SARS-CoV-2 variants.

Epigenetic regulation of HBV-specific tumor-infiltrating T cells in HBV-related HCC

BACKGROUND AND AIMS

HBV shapes the T-cell immune responses in HBV-related HCC. T cells can be recruited to the nidus, but limited T cells participate specifically in response to the HBV-related tumor microenvironment and HBV antigens. How epigenomic programs regulate T-cell compartments in virus-specific immune processes is unclear.

APPROACH AND RESULTS

We developed Ti-ATAC-seq. 2 to map the T-cell receptor repertoire, epigenomic, and transcriptomic landscape of αβ T cells at both the bulk-cell and single-cell levels in 54 patients with HCC. We deeply investigated HBV-specific T cells and HBV-related T-cell subsets that specifically responded to HBV antigens and the HBV + tumor microenvironment, respectively, characterizing their T-cell receptor clonality and specificity and performing epigenomic profiling. A shared program comprising NFKB1/2-, Proto-Oncogene, NF-KB Sub unit, NFATC2-, and NR4A1-associated unique T-cell receptor-downstream core epigenomic and transcriptomic regulome commonly regulated the differentiation of HBV-specific regulatory T-cell (Treg) cells and CD8 + exhausted T cells; this program was also selectively enriched in the HBV-related Treg-CTLA4 and CD8-exhausted T cell-thymocyte selection associated high mobility subsets and drove greater clonal expansion in HBV-related Treg-CTLA4 subset. Overall, 54% of the effector and memory HBV-specific T cells are governed by transcription factor motifs of activator protein 1, NFE2, and BACH1/2, which have been reported to be associated with prolonged patient relapse-free survival. Moreover, HBV-related tumor-infiltrating Tregs correlated with both increased viral titer and poor prognosis in patients.

CONCLUSIONS

This study provides insight into the cellular and molecular basis of the epigenomic programs that regulate the differentiation and generation of HBV-related T cells from viral infection and HBV + HCC unique immune exhaustion.

Hallmarks of CD8+ T cell dysfunction are established within hours of tumor antigen encounter before cell division

Tumor-specific CD8+ T cells (TST) in patients with cancer are dysfunctional and unable to halt cancer progression. TST dysfunction, also known as exhaustion, is thought to be driven by chronic T cell antigen receptor (TCR) stimulation over days to weeks. However, we know little about the interplay between CD8+ T cell function, cell division and epigenetic remodeling within hours of activation. Here, we assessed early CD8+ T cell differentiation, cell division, chromatin accessibility and transcription in tumor-bearing mice and acutely infected mice. Surprisingly, despite robust activation and proliferation, TST had near complete effector function impairment even before undergoing cell division and had acquired hallmark chromatin accessibility features previously associated with later dysfunction/exhaustion. Moreover, continued tumor/antigen exposure drove progressive epigenetic remodeling, 'imprinting' the dysfunctional state. Our study reveals the rapid divergence of T cell fate choice before cell division in the context of tumors versus infection.

Characterizing control of memory CD8 T cell differentiation by BTB-ZF transcription factor Zbtb20

Members of the BTB-ZF transcription factor family regulate the immune system. Our laboratory identified that family member Zbtb20 contributes to the differentiation, recall responses, and metabolism of CD8 T cells. Here, we report a characterization of the transcriptional and epigenetic signatures controlled by Zbtb20 at single-cell resolution during the effector and memory phases of the CD8 T cell response. Without Zbtb20, transcriptional programs associated with memory CD8 T cell formation were up-regulated throughout the CD8 T response. A signature of open chromatin was associated with genes controlling T cell activation, consistent with the known impact on differentiation. In addition, memory CD8 T cells lacking Zbtb20 were characterized by open chromatin regions with overrepresentation of AP-1 transcription factor motifs and elevated RNA- and protein-level expressions of the corresponding AP-1 components. Finally, we describe motifs and genomic annotations from the DNA targets of Zbtb20 in CD8 T cells identified by cleavage under targets and release under nuclease (CUT&RUN). Together, these data establish the transcriptional and epigenetic networks contributing to the control of CD8 T cell responses by Zbtb20.