Chromatin Profiling Research Articles

Single Cell Immune Atlas of Myalgic Encephalomyelitis/Chronic Fatigue Syndrome Identifies Substantial Transcriptional Dysregulation

Abstract Myalgic Encephalomyelitis/Chronic Fatigue Syndrome (ME/CFS) is a severe, debilitating disease with striking similarities to long-COVID/PASC, but its etiology remains unclear. Patients typically experience overwhleming fatigue, inflammatory symptoms and cognitive impairment, exacerbated by physical or mental exertion, also known as Post-Exertional Malaise (PEM). Mounting evidence points to immune system dysfunction being a key contributor to ME/CFS pathology, but our current understanding of gene regulation in ME/CFS leukocytes is hampered by a lack of quantitative and temporal resolution. To comprehensively characterize the gene regulatory state and precisely define the cell types underlying transcriptional dysregulation, we have generated single cell transcriptomes and chromatin profiles of peripheral blood mononuclear cells (PBMCs) using single cell RNA-seq and single cell ATAC-seq in matched patients and controls, at baseline and after exercise challenge. At baseline, ME/CFS patients displayed dysregulation of classical monocytes pointing to increased migration to tissues as well as inappropriate chemokine signaling, as well as dysregulation of certain subsets of effector CD8+ T cells, leading to a more inflammatory transcriptional program. Taken together, these data identify substantial immunological defects and serves as an important resource for understanding ME/CFS pathology.

Long-Read DNA Sequencing: Recent Advances and Remaining Challenges.

DNA sequencing has revolutionized medicine over recent decades. However, analysis of large structural variation and repetitive DNA, a hallmark of human genomes, has been limited by short-read technology, with read lengths of 100-300 bp. Long-read sequencing (LRS) permits routine sequencing of human DNA fragments tens to hundreds of kilobase pairs in size, using both real-time sequencing by synthesis and nanopore-based direct electronic sequencing. LRS permits analysis of large structural variation and haplotypic phasing in human genomes and has enabled the discovery and characterization of rare pathogenic structural variants and repeat expansions. It has also recently enabled the assembly of a complete, gapless human genome that includes previously intractable regions, such as highly repetitive centromeres and homologous acrocentric short arms. With the addition of protocols for targeted enrichment, direct epigenetic DNA modification detection, and long-range chromatin profiling, LRS promises to launch a new era of understanding of genetic diversity and pathogenic mutations in human populations.

OBERON3 and SUPPRESSOR OF MAX2 1-LIKE proteins form a regulatory module driving phloem development

Spatial specificity of cell fate decisions is central for organismal development. The phloem tissue mediates long-distance transport of energy metabolites along plant bodies and is characterized by an exceptional degree of cellular specialization. How a phloem-specific developmental program is implemented is, however, unknown. Here we reveal that the ubiquitously expressed PHD-finger protein OBE3 forms a central module with the phloem-specific SMXL5 protein for establishing the phloem developmental program in Arabidopsis thaliana. By protein interaction studies and phloem-specific ATAC-seq analyses, we show that OBE3 and SMXL5 proteins form a complex in nuclei of phloem stem cells where they promote a phloem-specific chromatin profile. This profile allows expression of OPS, BRX, BAM3, and CVP2 genes acting as mediators of phloem differentiation. Our findings demonstrate that OBE3/SMXL5 protein complexes establish nuclear features essential for determining phloem cell fate and highlight how a combination of ubiquitous and local regulators generate specificity of developmental decisions in plants.

Application of single cell multiomics points to changes in chromatin accessibility near calcitonin receptor like receptor and a possible role for adrenomedullin in the post-shock lung.

Acute lung injury (ALI)/acute respiratory distress syndrome (ARDS) is a commonly occurring sequelae of traumatic injury resulting from indirect insults like hypovolemic shock and/or extrapulmonary sepsis. The high lethality rate associated with these pathologies outlines the importance of clarifying the "priming" effects seen in the post-shock lung microenvironment, which are understood to bring about a dysregulated or overt immune response when triggered by a secondary systemic infectious/septic challenge culminating in ALI. In this pilot project, we test the hypothesis that application of a single cell multiomics approach can elucidate novel phenotype specific pathways potentially contributing to shock-induced ALI/ARDS. Hypovolemic shock was induced in C57BL/6 (wild-type), PD-1, PD-L1, or VISTA gene deficient male mice, 8-12 weeks old. Wild-type sham surgeries function as negative controls. A total of 24-h post-shock rodents were sacrificed, their lungs harvested and sectioned, with pools prepared from 2 mice per background, and flash frozen on liquid nitrogen. N = 2 biological replicates (representing 4 mice total) were achieved for all treatment groups across genetic backgrounds. Samples were received by the Boas Center for Genomics and Human Genetics, where single cell multiomics libraries were prepared for RNA/ATAC sequencing. The analysis pipeline Cell Ranger ARC was implemented to attain feature linkage assessments across genes of interest. Sham (pre-shock) results suggest high chromatin accessibility around calcitonin receptor like receptor (CALCRL) across cellular phenotypes with 17 and 18 feature links, exhibiting positive correlation with gene expression between biological replicates. Similarity between both sample chromatin profiles/linkage arcs is evident. Post-shock wild-type accessibility is starkly reduced across replicates where the number of feature links drops to 1 and 3, again presenting similar replicate profiles. Samples from shocked gene deficient backgrounds displayed high accessibility and similar profiles to the pre-shock lung microenvironment. High pre-shock availability of DNA segments and their positive correlation with CALCRL gene expression suggests an apparent regulatory capacity on transcription. Post-shock gene deficient chromatin profiles presented similar results to that of pre-shock wild-type samples, suggesting an influence on CALCRL accessibility. Key changes illustrated in the pre-ALI context of shock may allow for additional resolution of "priming" and "cellular pre-activation/pre-disposition" processes within the lung microenvironment.

Abstract 3523: A systematic comparison of molecular features shared by H3K27-altered diffuse midline gliomas and posterior fossa A ependymomas

Abstract In children, brain tumors account for the highest mortality rates of all malignancies, and in contrast to adult brain cancers, pediatric tumors often arise in infratentorial regions. Two tumors in particular, H3K27-altered diffuse midline gliomas (DMGs) and group-A posterior fossa ependymomas (PFAs) arise from structures around the fourth ventricle, and strikingly, both tumors exhibit a global decrease in the histone mark histone 3, lysine 27 tri-methylation (H3K27me3) that is associated with transcriptional repression. H3K27-altered DMGs primarily are found primarily in the pons and exhibit by lysine-to-methionine mutations at H3K27 (H3K27M). This missense mutation results in inhibition of the catalytic activity of the polycomb repressive complex 2 (PRC2) which mediates H3K27me3 deposition. PFAs are found in the lateral recess of the fourth ventricle and the cerebellum and demonstrate overexpression of EZH inhibitory protein (EZHIP). EZHIP contains a domain with striking resemblance to the mutant H3K27M histone tail and similarly blocks PRC2 activity. Although H3K27me3 levels are dramatically reduced in both areas, select genes retain H3K27me3 in both tumor types. Increasingly, data have emerged demonstrating that cases of PFAs actually harbor H3K27M mutations and cases of H3-wildtype DMGs in fact overexpress EZHIP. The neighboring origins and overlapping recurrent H3K27me3-modulating alterations observed in these two pediatric brain tumors suggest certain populations of cells in the developing hindbrain may be particularly prone to malignant transformation in the setting of PRC2 inhibition. A systematic comparison of genomic, transcriptomic, and epigenomic datasets for both tumors identified key molecular similarities between subsets of DMGs and PFAs and highlighted how developmental signatures may be preserved in tumor expression profiles. Specifically, these analyses identify shared copy-number profiles between PFAs and different subtypes of DMGs and demonstrate that the Activin A Receptor Type 1 (ACVR1), a gene with recurrent activating mutations in H3K27-altered DMGs, is overexpressed in a subset PFAs with worse survival outcomes. Additionally, interrogation of the H3K27me3 and enhancer landscapes identified critical genes with similar chromatin profiles and revealed heterogeneity in repression of H3K27me3-marked genes that correlated with patterns of inferred anatomic origins for each tumor. Moreover, one of these genes, cellular retinoic acid binding partner 1 (CRABP1) showed enrichment in progenitor cell types in single-cell transcriptomic datasets. Together, these analyses characterize the common and unique molecular features of H3K27-altered DMGs and PFAs and shed light on the cellular contexts of the developing hindbrain that lead to development of each tumor. Citation Format: Matthew Pun, Drew Pratt, Patricia R. Nano, Piyush K. Joshi, Li Jiang, Bernhard Englinger, Arvind Rao, Marcin Cieslik, Arul M. Chinnaiyan, Kenneth Aldape, Stefan Pfister, Mariella G. Filbin, Aparna Bhaduri, Sriram Venneti. A systematic comparison of molecular features shared by H3K27-altered diffuse midline gliomas and posterior fossa A ependymomas. [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 3523.

Abstract NG10: Chromatin profiles classify castration-resistant prostate cancers suggesting therapeutic targets

Abstract Untreated prostate cancers rely on androgen receptor (AR) signaling for growth and survival, forming the basis for the initial efficacy of androgen deprivation therapy (ADT). Yet the disease can relapse and progress to a lethal stage termed castration-resistant prostate cancer (CRPC). Reactivation of AR signaling represents the most common driver of CRPC growth, and next-generation AR signaling inhibitors (ARSIs) are now used in combination with ADT as first-line therapy. However, ARSIs can result in selective pressure, thereby generating AR-independent tumors. The transition from AR dependence frequently accompanies a change in a phenotype resembling developmental transdifferentiation or “lineage plasticity”. Neuroendocrine prostate cancer, which lacks a defined pathologic classification, is the most studied type of lineage plasticity. However, most AR-null tumors do not exhibit neuroendocrine features and are classified as “double-negative prostate cancer”, the drivers of which are poorly defined. Lineage plasticity studies in CRPC are limited by the lack of genetically defined patient-derived models that recapitulate the disease spectrum. To address this, we developed a biobank of organoids generated from patient biopsies to study the landscape of metastatic CRPC and allow for functional validation assays. Proteins called transcription factors (TFs) are drivers of tumor lineage plasticity. To identify the key TFs that drive the growth of AR-independent tumors, we integrated epigenetic and transcriptomic data generated from CRPC models. We generated ATAC-seq (assay for transposase-accessible chromatin sequencing) and RNA-seq data from 22 metastatic human prostate cancer organoids, six patient-derived xenografts (PDXs), and 12 derived or traditional cell lines. We classified the 40 models into four subtypes and predicted key TFs of each subtype. Besides the well-characterized AR-dependent (CRPC-AR) and neuroendocrine subtypes (CRPC-NE), we identified two novel AR-negative/low groups, including a Wnt-dependent subtype (CRPC-WNT), driven by TCF/LEF TFs, and a stem cell-like (SCL) subtype (CRPC-SCL), driven by the AP-1 family of TFs. To apply the subtype classification to patient samples, we derived RNA-seq signatures from the organoids and applied them to 366 patient samples from two independent CRPC cohorts. The generated signatures recapitulated the four-subtype classification and revealed that CRPC-SCL is the second most prevalent group. Patients from CRPC-SCL are also associated with shorter time under ARSI treatment compared to CRPC-AR, indicating that the ARSI treatments were less effective for CRPC-SCL patients. Additional chromatin immunoprecipitation sequencing (ChIP-seq) analysis indicated that AP-1 (FOSL1) collaboratively binds with TEAD and transcription coactivators, YAP and TAZ. Knocking down of AP-1 (FOSL1), YAP/TAZ decreased cell growth of CRPC-SCL and showed a decrease of chromatin accessibility at CRPC-SCL-specific open chromatin sites and down-regulation of YAP/TAZ target gene expression. In addition, the expression of AP-1 (FOSL1) decreased upon YAP/TAZ knockdown suggesting a positive feedback loop as well as YAP/TAZ as actional targets in CRPC-SCL. We used two small-molecule inhibitors, verteporfin and T-5224, that act on the YAP/TAZ/AP-1 pathway for their potential use as therapeutics for CRPC-SCL tumors, both inhibited the growth of samples from CRPC-SCL but not CRPC-AR. By overexpressing an AP-1 family gene (FOSL1) in AR-high cells, we observed an increase in chromatin accessibility at CRPC-SCL-specific open chromatin sites as well as significant up-regulation of CRPC-SCL signature genes, suggesting that AP-1 functions as a pioneering factor and master regulator for CRPC-SCL. All this work was recently published in Science (Tang, Xu et al. Science, 2022) where I am the co-first author. In summary, by using a diverse biobank of organoids, PDXs, and cell lines that recapitulate the heterogeneity of metastatic prostate cancer, we created a map of the chromatin accessibility and transcriptomic landscape of CRPC. We validated the CRPC-AR and CRPC-NE subtypes and report two novel subtypes of AR-negative/low samples, CRPC-SCL and CRPC-WNT, as well as their respective key TFs. Additional analysis revealed a model in which YAP, TAZ, TEAD, and AP-1 function together and drive oncogenic growth in CRPC-SCL samples. In addition, we proposed small inhibitors of YAP and TAZ that can potentially be used to treat CRPC-SCL patients. Overall, our results show how the stratification of CRPC patients into four subtypes using their transcriptomes can potentially inform appropriate clinical decisions. Citation Format: Fanying Tang, Duo Xu, Shangqian Wang, Chen Khuan Wong, Alexander Martinez-Fundichely, Cindy J. Lee, Sandra Cohen, Jane Park, Corinne E. Hill, Kenneth Eng, Rohan Bareja, Teng Han, Eric Minwei Liu, Ann Palladino, Wei Di, Dong Gao, Wassim Abida, Shaham Beg, Loredana Puca, Maximiliano Meneses, Elisa de Stanchina, Michael F. Berger, Anuradha Gopalan, Lukas E. Dow, Juan Miguel Mosquera, Himisha Beltran, Cora N. Sternberg, Ping Chi, Howard I. Scher, Andrea Sboner, Yu Chen, Ekta Khurana. Chromatin profiles classify castration-resistant prostate cancers suggesting therapeutic targets. [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr NG10.

Genome-Wide CRISPR-Cas9 Screen Identifies SMCHD1 as a Restriction Factor for Herpesviruses.

Intrinsic immunity is the frontline of host defense against invading pathogens. To combat viral infection, mammalian hosts deploy cell-intrinsic effectors to block viral replication prior to the onset of innate and adaptive immunity. In this study, SMCHD1 is identified as a pivotal cellular factor that restricts Kaposi's sarcoma-associated herpesvirus (KSHV) lytic reactivation through a genome-wide CRISPR-Cas9 knockout screen. Genome-wide chromatin profiling revealed that SMCHD1 associates with the KSHV genome, most prominently the origin of lytic DNA replication (ORI-Lyt). SMCHD1 mutants defective in DNA binding could not bind ORI-Lyt and failed to restrict KSHV lytic replication. Moreover, SMCHD1 functioned as a pan-herpesvirus restriction factor that potently suppressed a wide range of herpesviruses, including alpha, beta, and gamma subfamilies. SMCHD1 deficiency facilitated the replication of a murine herpesvirus in vivo. These findings uncovered SMCHD1 as a restriction factor against herpesviruses, and this could be harnessed for the development of antiviral therapies to limit viral infection. IMPORTANCE Intrinsic immunity represents the frontline of host defense against invading pathogens. However, our understanding of cell-intrinsic antiviral effectors remains limited. In this study, we identified SMCHD1 as a cell-intrinsic restriction factor that controlled KSHV lytic reactivation. Moreover, SMCHD1 restricted the replication of a wide range of herpesviruses by targeting the origins of viral DNA replication (ORIs), and SMCHD1 deficiency facilitated the replication of a murine herpesvirus in vivo. This study helps us to better understand intrinsic antiviral immunity, which may be harnessed to develop new therapeutics for the treatment of herpesvirus infection and the related diseases.

Cellcano: supervised cell type identification for single cell ATAC-seq data

Computational cell type identification is a fundamental step in single-cell omics data analysis. Supervised celltyping methods have gained increasing popularity in single-cell RNA-seq data because of the superior performance and the availability of high-quality reference datasets. Recent technological advances in profiling chromatin accessibility at single-cell resolution (scATAC-seq) have brought new insights to the understanding of epigenetic heterogeneity. With continuous accumulation of scATAC-seq datasets, supervised celltyping method specifically designed for scATAC-seq is in urgent need. Here we develop Cellcano, a computational method based on a two-round supervised learning algorithm to identify cell types from scATAC-seq data. The method alleviates the distributional shift between reference and target data and improves the prediction performance. After systematically benchmarking Cellcano on 50 well-designed celltyping tasks from various datasets, we show that Cellcano is accurate, robust, and computationally efficient. Cellcano is well-documented and freely available at https://marvinquiet.github.io/Cellcano/.

Conserved and variable heat stress responses of the Heat Shock Factor transcription factor family in maize and Setaria viridis.

The Heat Shock Factor (HSF) transcription factor family is a central and required component of plant heat stress responses and acquired thermotolerance. The HSF family has dramatically expanded in plant lineages, often including a repertoire of 20 or more genes. Here we assess and compare the composition, heat responsiveness, and chromatin profiles of the HSF families in maize and Setaria viridis (Setaria), two model C4 panicoid grasses. Both species encode a similar number of HSFs, and examples of both conserved and variable expression responses to a heat stress event were observed between the two species. Chromatin accessibility and genome-wide DNA-binding profiles were generated to assess the chromatin of HSF family members with distinct responses to heat stress. We observed significant variability for both chromatin accessibility and promoter occupancy within similarly regulated sets of HSFs between Setaria and maize, as well as between syntenic pairs of maize HSFs retained following its most recent genome duplication event. Additionally, we observed the widespread presence of TF binding at HSF promoters in control conditions, even at HSFs that are only expressed in response to heat stress. TF-binding peaks were typically near putative HSF-binding sites in HSFs upregulated in response to heat stress, but not in stable or not expressed HSFs. These observations collectively support a complex scenario of expansion and subfunctionalization within this transcription factor family and suggest that within-family HSF transcriptional regulation is a conserved, defining feature of the family.

Graded BMP signaling within intestinal crypt architecture directs self-organization of the Wnt-secreting stem cell niche.

Signals from the surrounding niche drive proliferation and suppress differentiation of intestinal stem cells (ISCs) at the bottom of intestinal crypts. Among sub-epithelial support cells, deep sub-cryptal CD81+ PDGFRAlo trophocytes capably sustain ISC functions exvivo. Here, we show that mRNA and chromatin profiles of abundant CD81- PDGFRAlo mouse stromal cells resemble those of trophocytes and that both populations provide crucial canonical Wnt ligands. Mesenchymal expression of key ISC-supportive factors extends along a spatial and molecular continuum from trophocytes into peri-cryptal CD81- CD55hi cells, which mimic trophocyte activity in organoid co-cultures. Graded expression of essential niche factors is not cell-autonomous but dictated by the distance from bone morphogenetic protein (BMP)-secreting PDGFRAhi myofibroblast aggregates. BMP signaling inhibits ISC-trophic genes in PDGFRAlo cells near high crypt tiers; that suppression is relieved in stromal cells near and below the crypt base, including trophocytes. Cell distances thus underlie a self-organized and polar ISC niche.

Chromatin activation profiling of stereotyped chronic lymphocytic leukemias reveals a subset 8–specific signature

AbstractThe chromatin activation landscape of chronic lymphocytic leukemia (CLL) with stereotyped B-cell receptor immunoglobulin is currently unknown. In this study, we report the results of a whole-genome chromatin profiling of histone 3 lysine 27 acetylation of 22 CLLs from major subsets, which were compared against nonstereotyped CLLs and normal B-cell subpopulations. Although subsets 1, 2, and 4 did not differ much from their nonstereotyped CLL counterparts, subset 8 displayed a remarkably distinct chromatin activation profile. In particular, we identified 209 de novo active regulatory elements in this subset, which showed similar patterns with U-CLLs undergoing Richter transformation. These regions were enriched for binding sites of 9 overexpressed transcription factors. In 78 of 209 regions, we identified 113 candidate overexpressed target genes, 11 regions being associated with more than 2 adjacent genes. These included blocks of up to 7 genes, suggesting local coupregulation within the same genome compartment. Our findings further underscore the uniqueness of subset 8 CLL, notable for the highest risk of Richter’s transformation among all CLLs and provide additional clues to decipher the molecular basis of its clinical behavior.

Transcriptomic and Chromatin Landscape Analysis Reveals That Involvement of Pituitary Level Transcription Factors Modulate Incubation Behaviors of Magang Geese.

The incubation behavior of geese seriously affects their egg production performance. Studies on incubation behavior have identified functional genes, but the regulatory architecture relationship between functional genes and chromatin accessibility remains poorly understood. Here, we present an integrated analysis of open chromatin profiles and transcriptome to identify the cis-regulatory element and their potential transcription factors involved in regulating incubation behavior in goose pituitary. Assay for transposase-accessible chromatin sequencing (ATAC-seq) revealed that open chromatin regions increased in the pituitary during the transition from incubation behavior to laying. We identified 920 significant differential accessible regions (DARs) in the pituitary. Compared to the laying stage, most DARs had higher chromatin accessibility in the brooding stage. Motif analysis of open DARs showed that the most significant transcription factor (TF) occupied sites predominantly enriched in motifs binding to the RFX family (RFX5, RFX2, and RFX1). While the majority of TF motifs enriched under sites of the nuclear receptor (NR) family (ARE, GRE, and PGR) in closed DARs at the incubation behavior stage. Footprint analysis indicated that the transcription factor RFX family exhibited higher binding on chromatin at the brooding stage. To further elucidate the effect of changes in chromatin accessibility on gene expression levels, a comparison of the transcriptome revealed 279 differentially expressed genes (DEGs). The transcriptome changes were associated with processes of steroid biosynthesis. By integrating ATAC-seq and RNA-seq, few DARs directly affect incubation behavior by regulating the transcription levels of genes. Five DAR-related DEGs were found to be closely related to maintaining the incubation behavior in geese. Footprinting analysis revealed a set of transcription factors (RFX1, RFX2, RFX3, RFX5, BHLHA15, SIX1, and DUX) which displayed the highest activity at the brooding stage. SREBF2 was predicted to be the unique differentially expressed transcription factor whose mRNA level was down-regulated and enriched in hyper-accessible regions of PRL in the broody stage. In the present study, we comprehensively profiled the transcriptome and chromatin accessibility in the pituitary related to incubation behavior. Our findings provided insight into the identification and analysis of regulatory elements in goose incubation behavior. The epigenetic alterations profiled here can help decipher the epigenetic mechanisms that contribute to the regulation of incubation behavior in birds.

Comprehensive molecular phenotyping of ARID1A-deficient gastric cancer reveals pervasive epigenomic reprogramming and therapeutic opportunities

ObjectiveGastric cancer (GC) is a leading cause of cancer mortality, with ARID1A being the second most frequently mutated driver gene in GC. We sought to decipher ARID1A-specific GC regulatory networks...

Distinct roles for canonical and variant histone H3 lysine-36 in Polycomb silencing.

Polycomb complexes regulate cell type-specific gene expression programs through heritable silencing of target genes. Trimethylation of histone H3 lysine 27 (H3K27me3) is essential for this process. Perturbation of H3K36 is thought to interfere with H3K27me3. We show that mutants of Drosophila replication-dependent (H3.2K36R) or replication-independent (H3.3K36R) histone H3 genes generally maintain Polycomb silencing and reach later stages of development. In contrast, combined (H3.3K36RH3.2K36R) mutants display widespread Hox gene misexpression and fail to develop past the first larval stage. Chromatin profiling revealed that the H3.2K36R mutation disrupts H3K27me3 levels broadly throughout silenced domains, whereas these regions are mostly unaffected in H3.3K36R animals. Analysis of H3.3 distributions showed that this histone is enriched at presumptive Polycomb response elements located outside of silenced domains but relatively depleted from those inside. We conclude that H3.2 and H3.3 K36 residues collaborate to repress Hox genes using different mechanisms.

Fluctuations in chromatin state at regulatory loci occur spontaneously under relaxed selection and are associated with epigenetically inherited variation in C. elegans gene expression.

Some epigenetic information can be transmitted between generations without changes in the underlying DNA sequence. Changes in epigenetic regulators, termed epimutations, can occur spontaneously and be propagated in populations in a manner reminiscent of DNA mutations. Small RNA-based epimutations occur in C. elegans and persist for around 3-5 generations on average. Here, we explored whether chromatin states also undergo spontaneous change and whether this could be a potential alternative mechanism for transgenerational inheritance of gene expression changes. We compared the chromatin and gene expression profiles at matched time points from three independent lineages of C. elegans propagated at minimal population size. Spontaneous changes in chromatin occurred in around 1% of regulatory regions each generation. Some were heritable epimutations and were significantly enriched for heritable changes in expression of nearby protein-coding genes. Most chromatin-based epimutations were short-lived but a subset had longer duration. Genes subject to long-lived epimutations were enriched for multiple components of xenobiotic response pathways. This points to a possible role for epimutations in adaptation to environmental stressors.

The molecular consequences of androgen activity in the human breast.

Estrogen and progesterone have been extensively studied in the mammary gland, but the molecular effects of androgen remain largely unexplored. Transgender men are recorded as female at birth but identify as male and may undergo gender-affirming androgen therapy to align their physical characteristics and gender identity. Here we perform single-cell-resolution transcriptome, chromatin, and spatial profiling of breast tissues from transgender men following androgen therapy. We find canonical androgen receptor gene targets are upregulated in cells expressing the androgen receptor and that paracrine signaling likely drives sex-relevant androgenic effects in other cell types. We also observe involution of the epithelium and a spatial reconfiguration of immune, fibroblast, and vascular cells, and identify a gene regulatory network associated with androgen-induced fat loss. This work elucidates the molecular consequences of androgen activity in the human breast at single-cell resolution.

Chromatin profiling identifies transcriptional readthrough as a conserved mechanism for piRNA biogenesis in mosquitoes

The piRNA pathway in mosquitoes differs substantially from other model organisms, with an expanded PIWI gene family and functions in antiviral defense. Here, we define core piRNA clusters as genomic loci that show ubiquitous piRNA expression in both somatic and germline tissues. These core piRNA clusters are enriched for non-retroviral endogenous viral elements (nrEVEs) in antisense orientation and depend on key biogenesis factors, Veneno, Tejas, Yb, and Shutdown. Combined transcriptome and chromatin state analyses identify transcriptional readthrough as a conserved mechanism for cluster-derived piRNA biogenesis in the vector mosquitoes Aedes aegypti, Aedes albopictus, Culex quinquefasciatus, and Anopheles gambiae. Comparative analyses between the two Aedes species suggest that piRNA clusters function as traps for nrEVEs, allowing adaptation to environmental challenges such as virus infection. Our systematic transcriptome and chromatin state analyses lay the foundation for studies of gene regulation, genome evolution, and piRNA function in these important vector species.

Mitochondrial single-cell ATAC-seq for high-throughput multi-omic detection of mitochondrial genotypes and chromatin accessibility.

Natural sequence variation within mitochondrial DNA (mtDNA) contributes to human phenotypes and may serve as natural genetic markers in human cells for clonal and lineage tracing. We recently developed a single-cell multi-omic approach, called 'mitochondrial single-cell assay for transposase-accessible chromatin with sequencing' (mtscATAC-seq), enabling concomitant high-throughput mtDNA genotyping and accessible chromatin profiling. Specifically, our technique allows the mitochondrial genome-wide inference of mtDNA variant heteroplasmy along with information on cell state and accessible chromatin variation in individual cells. Leveraging somatic mtDNA mutations, our method further enables inference of clonal relationships among native ex vivo-derived human cells not amenable to genetic engineering-based clonal tracing approaches. Here, we provide a step-by-step protocol for the use of mtscATAC-seq, including various cell-processing and flow cytometry workflows, by using primary hematopoietic cells, subsequent single-cell genomic library preparation and sequencing that collectively take ~3-4 days to complete. We discuss experimental and computational data quality control metrics and considerations for the extension to other mammalian tissues. Overall, mtscATAC-seq provides a broadly applicable platform to map clonal relationships between cells in human tissues, investigate fundamental aspects of mitochondrial genetics and enable additional modes of multi-omic discovery.

Aberrant JAK-STAT signaling-mediated chromatin remodeling impairs the sensitivity of NK/T-cell lymphoma to chidamide

BackgroundNatural killer/T-cell lymphoma (NKTL) is a rare type of aggressive and heterogeneous non-Hodgkin's lymphoma (NHL) with a poor prognosis and limited therapeutic options. Therefore, there is an urgent need to exploit potential novel therapeutic targets for the treatment of NKTL. Histone deacetylase (HDAC) inhibitor chidamide was recently approved for treating relapsed/refractory peripheral T-cell lymphoma (PTCL) patients. However, its therapeutic efficacy in NKTL remains unclear.MethodsWe performed a phase II clinical trial to evaluate the efficacy of chidamide in 28 relapsed/refractory NKTL patients. Integrative transcriptomic, chromatin profiling analysis and functional studies were performed to identify potential predictive biomarkers and unravel the mechanisms of resistance to chidamide. Immunohistochemistry (IHC) was used to validate the predictive biomarkers in tumors from the clinical trial.ResultsWe demonstrated that chidamide is effective in treating relapsed/refractory NKTL patients, achieving an overall response and complete response rate of 39 and 18%, respectively. In vitro studies showed that hyperactivity of JAK-STAT signaling in NKTL cell lines was associated with the resistance to chidamide. Mechanistically, our results revealed that aberrant JAK-STAT signaling remodels the chromatin and confers resistance to chidamide. Subsequently, inhibition of JAK-STAT activity could overcome resistance to chidamide by reprogramming the chromatin from a resistant to sensitive state, leading to synergistic anti-tumor effect in vitro and in vivo. More importantly, our clinical data demonstrated that combinatorial therapy with chidamide and JAK inhibitor ruxolitinib is effective against chidamide-resistant NKTL. In addition, we identified TNFRSF8 (CD30), a downstream target of the JAK-STAT pathway, as a potential biomarker that could predict NKTL sensitivity to chidamide.ConclusionsOur study suggests that chidamide, in combination with JAK-STAT inhibitors, can be a novel targeted therapy in the standard of care for NKTL.Trial registration: ClinicalTrials.gov, NCT02878278. Registered 25 August 2016, https://clinicaltrials.gov/ct2/show/NCT02878278

Cistrome Explorer: an interactive visual analysis tool for large-scale epigenomic data.

The regulation of genes by cis-regulatory elements (CREs) is complex and differs between cell types. Visual analysis of large collections of chromatin profiles across diverse cell types, integrated with computational methods, can reveal meaningful biological insights. We developed Cistrome Explorer, a web-based interactive visual analytics tool for exploring thousands of chromatin profiles in diverse cell types. Integrated with the Cistrome Data Browser database which contains thousands of ChIP-seq, DNase-seq and ATAC-seq samples, Cistrome Explorer enables the discovery of patterns of CREs across cell types and the identification of transcription factor binding underlying these patterns. Cistrome Explorer and its source code are available at http://cisvis.gehlenborglab.org/ and released under the MIT License. Documentation can be accessed via http://cisvis.gehlenborglab.org/docs/. Supplementary data are available at Bioinformatics online.