WiggleTools: parallel processing of large collections of genome-wide datasets for visualization and statistical analysis

CRWN nuclear lamina components maintain the H3K27me3 landscape and promote successful reproduction in Arabidopsis

Arabidopsis lamin analogs CROWDED NUCLEIs (CRWNs) are necessary to maintain nuclear structure, genome function, and proper plant growth. However, whether and how CRWNs impact reproduction and genome-wide epigenetic modifications is unknown. Here, we investigate the role of CRWNs during the development of gametophytes, seeds, and endosperm, using genomic and epigenomic profiling methods. We observed defects in crwn mutant seeds including seed abortion and reduced germination rate. Quadruple crwn null genotypes were rarely transmitted through gametophytes. Because defects in seeds often stem from abnormal endosperm development, we focused on crwn1 crwn2 (crwn1/2) endosperm. These mutant seeds exhibited enlarged chalazal endosperm cysts and increased expression of stress-related genes and the MADS-box transcription factor PHERES1 and its targets. Previously, it was shown that PHERES1 expression is regulated by H3K27me3 and that CRWN1 interacts with the PRC2 interactor PWO1. Thus, we tested whether crwn1/2 alters H3K27me3 patterns. We observed a mild loss of H3K27me3 at several hundred loci, which differed between endosperm and leaves. These data indicate that CRWNs are necessary to maintain the H3K27me3 landscape, with tissue-specific chromatin and transcriptional consequences.

The Rtf1/Prf1-dependent histone modification axis counteracts multi-drug resistance in fission yeast

RNA polymerase II transcription elongation directs an intricate pattern of histone modifications. This pattern includes a regulatory cascade initiated by the elongation factor Rtf1, leading to monoubiquitylation of histone H2B, and subsequent methylation of histone H3 on lysine 4. Previous studies have defined the molecular basis for these regulatory relationships, but it remains unclear how they regulate gene expression. To address this question, we investigated a drug resistance phenotype that characterizes defects in this axis in the model eukaryote Schizosaccharomyces pombe (fission yeast). The mutations caused resistance to the ribonucleotide reductase inhibitor hydroxyurea (HU) that correlated with a reduced effect of HU on dNTP pools, reduced requirement for the S-phase checkpoint, and blunting of the transcriptional response to HU treatment. Mutations in the C-terminal repeat domain of the RNA polymerase II large subunit Rpb1 led to similar phenotypes. Moreover, all the HU-resistant mutants also exhibited resistance to several azole-class antifungal agents. Our results suggest a novel, shared gene regulatory function of the Rtf1-H2Bub1-H3K4me axis and the Rpb1 C-terminal repeat domain in controlling fungal drug tolerance.

SIV infection and ARV treatment reshape the transcriptional and epigenetic profile of naïve and memory T cells in vivo

Human and simian immunodeficiency viruses (HIV and SIV) are lentiviruses that reverse transcribe their RNA genome with subsequent integration into the genome of the target cell. How progressive infection and administration of antiretrovirals (ARVs) longitudinally influence the transcriptomic and epigenetic landscape of particular T cell subsets, and how these may influence the genetic location of integration are unclear. Here, we use RNAseq and ATACseq to study the transcriptomics and epigenetic landscape of longitudinally sampled naïve and memory CD4+ and CD8+ T cells in two species of non-human primates prior to SIV infection, during chronic SIV infection, and after administration of ARVs. We find that SIV infection leads to significant alteration to the transcriptomic profile of all T cell subsets that are only partially reversed by administration of ARVs. Epigenetic changes were more apparent in animals with longer periods of untreated SIV infection and correlated well with changes in corresponding gene expression. Known SIV integration sites did not vary due to SIV status but did contain more open chromatin in rhesus macaque memory T cells, and the expression of proteasome-related genes at the pre-SIV timepoint correlated with subsequent viremia.IMPORTANCEChronic inflammation during progressive human and simian immunodeficiency virus (HIV and SIV) infections leads to significant co-morbidities in infected individuals with significant consequences. Antiretroviral (ARV)-treated individuals also manifest increased levels of inflammation which are associated with increased mortalities. These data will help guide rational development of modalities to reduce inflammation observed in people living with HIV and suggest mechanisms underlying lentiviral integration site preferences.

The plant POLYMERASE-ASSOCIATED FACTOR1 complex links transcription and H2B monoubiquitination genome wide

The evolutionarily conserved POLYMERASE-ASSOCIATED FACTOR1 complex (Paf1C) participates in transcription, and research in animals and fungi suggests that it facilitates RNA POLYMERASE II (RNAPII) progression through chromatin. We examined the genomic distribution of the EARLY FLOWERING7 (ELF7) and VERNALIZATION INDEPENDENCE3 subunits of Paf1C in Arabidopsis (Arabidopsis thaliana). The occupancy of both subunits was confined to thousands of gene bodies and positively associated with RNAPII occupancy and the level of gene expression, supporting a role as a transcription elongation factor. We found that monoubiquitinated histone H2B, which marks most transcribed genes, was strongly reduced genome wide in elf7 seedlings. Genome-wide profiling of RNAPII revealed that in elf7 mutants, RNAPII occupancy was reduced throughout the gene body and at the transcription end site of Paf1C-targeted genes, suggesting a direct role for the complex in transcription elongation. Overall, our observations suggest a direct functional link between Paf1C activity, monoubiquitination of histone H2B, and the transition of RNAPII to productive elongation. However, for several genes, Paf1C may also act independently of H2Bub deposition or occupy these genes more stable than H2Bub marking, possibly reflecting the dynamic nature of Paf1C association and H2Bub turnover during transcription.

Dynamic co-evolution of transposable elements and the piRNA pathway in African cichlid fishes

East African cichlid fishes have diversified in an explosive fashion, but the (epi)genetic basis of the phenotypic diversity of these fishes remains largely unknown. Although transposable elements (TEs) have been associated with phenotypic variation in cichlids, little is known about their transcriptional activity and epigenetic silencing. Here, we describe dynamic patterns of TE expression in African cichlid gonads and during early development. Orthology inference revealed an expansion of piwil1 genes in Lake Malawi cichlids, likely driven by PiggyBac TEs. The expanded piwil1 copies have signatures of positive selection and retain amino acid residues essential for catalytic activity. Furthermore, the gonads of African cichlids express a Piwi-interacting RNA (piRNA) pathway that target TEs. We define the genomic sites of piRNA production in African cichlids and find divergence in closely related species, in line with fast evolution of piRNA-producing loci. Our findings suggest dynamic co-evolution of TEs and host silencing pathways in the African cichlid radiations. We propose that this co-evolution has contributed to cichlid genomic diversity.

Thiolutin has complex effects in vivo but is a direct inhibitor of RNA polymerase II in vitro

Thiolutin is a natural product transcription inhibitor with an unresolved mode of action. Thiolutin and the related dithiolopyrrolone holomycin chelate Zn2+ and previous studies have concluded that RNA Polymerase II (Pol II) inhibition in vivo is indirect. Here, we present chemicogenetic and biochemical approaches to investigate thiolutin's mode of action in Saccharomyces cerevisiae. We identify mutants that alter sensitivity to thiolutin. We provide genetic evidence that thiolutin causes oxidation of thioredoxins in vivo and that thiolutin both induces oxidative stress and interacts functionally with multiple metals including Mn2+ and Cu2+, and not just Zn2+. Finally, we show direct inhibition of RNA polymerase II (Pol II) transcription initiation by thiolutin in vitro in support of classical studies that thiolutin can directly inhibit transcription in vitro. Inhibition requires both Mn2+ and appropriate reduction of thiolutin as excess DTT abrogates its effects. Pause prone, defective elongation can be observed in vitro if inhibition is bypassed. Thiolutin effects on Pol II occupancy in vivo are widespread but major effects are consistent with prior observations for Tor pathway inhibition and stress induction, suggesting that thiolutin use in vivo should be restricted to studies on its modes of action and not as an experimental tool.

Oncogenic enhancers prime quiescent metastatic cells to escape NK immune surveillance by eliciting transcriptional memory

Metastasis arises from disseminated tumour cells (DTCs) that are characterized by intrinsic phenotypic plasticity and the capability of seeding to secondary organs. DTCs can remain latent for years before giving rise to symptomatic overt metastasis. In this context, DTCs fluctuate between a quiescent and proliferative state in response to systemic and microenvironmental signals including immune-mediated surveillance. Despite its relevance, how intrinsic mechanisms sustain DTCs plasticity has not been addressed. By interrogating the epigenetic state of metastatic cells, we find that tumour progression is coupled with the activation of oncogenic enhancers that are organized in variable interconnected chromatin domains. This spatial chromatin context leads to the activation of a robust transcriptional response upon repeated exposure to retinoic acid (RA). We show that this adaptive mechanism sustains the quiescence of DTCs through the activation of the master regulator SOX9. Finally, we determine that RA-stimulated transcriptional memory increases the fitness of metastatic cells by supporting the escape of quiescent DTCs from NK-mediated immune surveillance. Overall, these findings highlight the contribution of oncogenic enhancers in establishing transcriptional memories as an adaptive mechanism to reinforce cancer dormancy and immune escape, thus amenable for therapeutic intervention.

Transcriptome analysis of Edwardsiella piscicida during intracellular infection reveals excludons are involved with the activation of a mitochondrion-like energy generation program

Phylogenetic evidence suggests a shared ancestry between mitochondria and modern Proteobacteria, a phylum including several genera of intracellular pathogens. Studying these diverse pathogens, particularly during intracellular infection of their hosts, can reveal characteristics potentially representative of the mitochondrial-Proteobacterial ancestor by identifying traits shared with mitochondria. While transcriptomic approaches can provide global insights into intracellular acclimatization by pathogens, they are often limited by excess host RNAs in extracts. Here, we developed a method employing magnetic nanoparticles to enrich RNA from an intracellular Gammaproteobacterium, Edwardsiella piscicida, within zebrafish, Danio rerio, fin fibroblasts, enabling comprehensive exploration of the bacterial transcriptome. Our findings revealed that the intracellular E. piscicida transcriptome reflects a mitochondrion-like energy generation program characterized by the suppression of glycolysis and sugar transport, coupled with upregulation of the tricarboxylic acid (TCA) cycle and alternative import of simple organic acids that directly flux into TCA cycle intermediates or electron transport chain donors. Additionally, genes predicted to be members of excludons, loci of gene pairs antagonistically co-regulated by overlapping antisense transcription, are significantly enriched in the set of all genes with perturbed sense and antisense transcription, suggesting a general but important involvement of excludons with intracellular acclimatization. Notably, genes involved with the activation of the mitochondrion-like energy generation program, specifically with metabolite import and glycolysis, are also members of predicted excludons. Other intracellular Proteobacterial pathogens appear to employ a similar mitochondrion-like energy generation program, suggesting a potentially conserved mechanism for optimized energy acquisition from hosts centered around the TCA cycle.IMPORTANCEPhylogenetic evidence suggests that mitochondria and Proteobacteria, a phylum encompassing various intracellular pathogens, share a common ancestral lineage. In this study, we developed a novel method employing magnetic nanoparticles to explore the transcriptome of an aquatic Gammaproteobacterium, Edwardsiella piscicida, during intracellular infection of host cells. We show that the strategy E. piscicida uses to generate energy strikingly mirrors the function of mitochondria-energy generators devoid of glycolytic processes. Notably, several implicated genes are members of excludons-gene pairs antagonistically co-regulated by overlapping antisense transcription. Other intracellular Proteobacterial pathogens appear to adopt a similar mitochondrion-like energy generation program, indicating a possibly conserved strategy for optimized energy acquisition from hosts centered around the tricarboxylic acid cycle.

Rarγ -Foxa1 signaling promotes luminal identity in prostate progenitors and is disrupted in prostate cancer

Retinoic acid (RA) signaling is a master regulator of vertebrate development with crucial roles in directing body axis orientation and tissue differentiation, including in the reproductive system. However, a mechanistic understanding of how RA signaling promotes cell lineage identity in different tissues is often missing. Here, leveraging prostate organoid technology, we demonstrated that RA signaling orchestrates the commitment of adult mouse prostate progenitors to glandular identity, epithelial barrier integrity, and ultimately, proper specification of the prostatic lumen. Mechanistically, RA-dependent RARγ activation promotes the expression of the pioneer factor Foxa1, which synergizes with the androgen pathway for proper luminal expansion, cytoarchitecture and function. FOXA1 nucleotide variants are common in human prostate and breast cancers and considered driver mutations, though their pathogenic mechanism is incompletely understood. Combining functional genetics experiments with structural modeling of FOXA1 folding and chromatin binding analyses, we discovered that FOXA1 F254E255 is a loss-of-function mutation leading to compromised transcriptional function and lack of luminal fate commitment of prostate progenitors. Overall, we define RA as a crucial instructive signal for glandular identity in adult prostate progenitors. We propose deregulation of vitamin A metabolism as a risk factor for benign and malignant prostate disease, and identified cancer associated FOXA1 indels affecting residue F254 as loss-of-function mutations promoting dedifferentiation of adult prostate progenitors. Summary: Retinoic acid signaling orchestrates luminal differentiation of adult prostate progenitors.

A GREB1-steroid receptor feedforward mechanism governs differential GREB1 action in endometrial function and endometriosis

Cellular responses to the steroid hormones, estrogen (E2), and progesterone (P4) are governed by their cognate receptor's transcriptional output. However, the feed-forward mechanisms that shape cell-type-specific transcriptional fulcrums for steroid receptors are unidentified. Herein, we found that a common feed-forward mechanism between GREB1 and steroid receptors regulates the differential effect of GREB1 on steroid hormones in a physiological or pathological context. In physiological (receptive) endometrium, GREB1 controls P4-responses in uterine stroma, affecting endometrial receptivity and decidualization, while not affecting E2-mediated epithelial proliferation. Of mechanism, progesterone-induced GREB1 physically interacts with the progesterone receptor, acting as a cofactor in a positive feedback mechanism to regulate P4-responsive genes. Conversely, in endometrial pathology (endometriosis), E2-induced GREB1 modulates E2-dependent gene expression to promote the growth of endometriotic lesions in mice. This differential action of GREB1 exerted by a common feed-forward mechanism with steroid receptors advances our understanding of mechanisms that underlie cell- and tissue-specific steroid hormone actions.

ARID1A governs the silencing of sex-linked transcription during male meiosis in the mouse

We present evidence implicating the BAF (BRG1/BRM Associated Factor) chromatin remodeler in meiotic sex chromosome inactivation (MSCI). By immunofluorescence (IF), the putative BAF DNA binding subunit, ARID1A (AT-rich Interaction Domain 1a), appeared enriched on the male sex chromosomes during diplonema of meiosis I. Those germ cells showing a Cre-induced loss of ARID1A were arrested in pachynema and failed to repress sex-linked genes, indicating a defective MSCI. Consistent with this defect, mutant sex chromosomes displayed an abnormal presence of elongating RNA polymerase II coupled with an overall increase in chromatin accessibility detectable by ATAC-seq. By investigating potential mechanisms underlying these anomalies, we identified a role for ARID1A in promoting the preferential enrichment of the histone variant, H3.3, on the sex chromosomes, a known hallmark of MSCI. Without ARID1A, the sex chromosomes appeared depleted of H3.3 at levels resembling autosomes. Higher resolution analyses by CUT&RUN revealed shifts in sex-linked H3.3 associations from discrete intergenic sites and broader gene-body domains to promoters in response to the loss of ARID1A. Several sex-linked sites displayed ectopic H3.3 occupancy that did not co-localize with DMC1 (DNA Meiotic Recombinase 1). This observation suggests a requirement for ARID1A in DMC1 localization to the asynapsed sex chromatids. We conclude that ARID1A-directed H3.3 localization influences meiotic sex chromosome gene regulation and DNA repair.

Stat5 opposes the transcription factor Tox and rewires exhausted CD8+ T cells toward durable effector-like states during chronic antigen exposure

Rewiring exhausted CD8+ T (Tex) cells toward functional states remains a therapeutic challenge. Tex cells are epigenetically programmed by the transcription factor Tox. However, epigenetic remodeling occurs as Tex cells transition from progenitor (Texprog) to intermediate (Texint) and terminal (Texterm) subsets, suggesting development flexibility. We examined epigenetic transitions between Tex cell subsets and revealed a reciprocally antagonistic circuit between Stat5a and Tox. Stat5 directed Texint cell formation and re-instigated partial effector biology during this Texprog-to-Texint cell transition. Constitutive Stat5a activity antagonized Tox and rewired CD8+ T cells from exhaustion to a durable effector and/or natural killer (NK)-like state with superior anti-tumor potential. Temporal induction of Stat5 activity in Tex cells using an orthogonal IL-2:IL2Rβ-pair fostered Texint cell accumulation, particularly upon PD-L1 blockade. Re-engaging Stat5 also partially reprogrammed the epigenetic landscape of exhaustion and restored polyfunctionality. These data highlight therapeutic opportunities of manipulating the IL-2-Stat5 axis to rewire Tex cells toward more durably protective states.

Temporal evolution reveals bifurcated lineages in aggressive neuroendocrine small cell prostate cancer trans-differentiation

Trans-differentiation from an adenocarcinoma to a small cell neuroendocrine state is associated with therapy resistance in multiple cancer types. To gain insight into the underlying molecular events of the trans-differentiation, we perform a multi-omics time course analysis of a pan-small cell neuroendocrine cancer model (termed PARCB), a forward genetic transformation using human prostate basal cells and identify a shared developmental, arc-like, and entropy-high trajectory among all transformation model replicates. Further mapping with single cell resolution reveals two distinct lineages defined by mutually exclusive expression of ASCL1 or ASCL2. Temporal regulation by groups of transcription factors across developmental stages reveals that cellular reprogramming precedes the induction of neuronal programs. TFAP4 and ASCL1/2 feedback are identified as potential regulators of ASCL1 and ASCL2 expression. Our study provides temporal transcriptional patterns and uncovers pan-tissue parallels between prostate and lung cancers, as well as connections to normal neuroendocrine cell states.

The BRD4-NUT Fusion Alone Drives Malignant Transformation of NUT Carcinoma

NUT carcinoma (NC) is an aggressive squamous carcinoma defined by the BRD4-NUT fusion oncoprotein. Routinely effective systemic treatments are unavailable for most NC patients. The lack of an adequate animal model precludes identifying and leveraging cell-extrinsic factors therapeutically in NC. Here, we created a genetically engineered mouse model (GEMM) of NC that forms a Brd4::NUTM1 fusion gene upon tamoxifen induction of Sox2-driven Cre. The model displayed complete disease penetrance, with tumors arising from the squamous epithelium weeks after induction and all mice succumbing to the disease shortly thereafter. Closely resembling human NC (hNC), GEMM tumors (mNC) were poorly differentiated squamous carcinomas with high expression of MYC that metastasized to solid organs and regional lymph nodes. Two GEMM-derived cell lines were developed whose transcriptomic and epigenetic landscapes harbored key features of primary GEMM tumors. Importantly, GEMM tumor and cell line transcriptomes co-classified with those of human NC. BRD4-NUT also blocked differentiation and maintained the growth of mNC as in hNC. Mechanistically, GEMM primary tumors and cell lines formed large histone H3K27ac-enriched domains, termed megadomains, that were invariably associated with the expression of key NC-defining proto-oncogenes, Myc and Trp63. Small-molecule BET bromodomain inhibition (BETi) of mNC induced differentiation and growth arrest and prolonged survival of NC GEMMs, as it does in hNC models. Overall, tumor formation in the NC GEMM is definitive evidence that BRD4-NUT alone can potently drive the malignant transformation of squamous progenitor cells into NC.

SIGNIFICANCE

The development of an immunocompetent model of NUT carcinoma that closely mimics the human disease provides a valuable global resource for mechanistic and preclinical studies to improve treatment of this incurable disease.

MSL2 ensures biallelic gene expression in mammals

In diploid organisms, biallelic gene expression enables the production of adequate levels of mRNA1,2. This is essential for haploinsufficient genes, which require biallelic expression for optimal function to prevent the onset of developmental disorders1,3. Whether and how a biallelic or monoallelic state is determined in a cell-type-specific manner at individual loci remains unclear. MSL2 is known for dosage compensation of the male X chromosome in flies. Here we identify a role of MSL2 in regulating allelic expression in mammals. Allele-specific bulk and single-cell analyses in mouse neural progenitor cells revealed that, in addition to the targets showing biallelic downregulation, a class of genes transitions from biallelic to monoallelic expression after MSL2 loss. Many of these genes are haploinsufficient. In the absence of MSL2, one allele remains active, retaining active histone modifications and transcription factor binding, whereas the other allele is silenced, exhibiting loss of promoter-enhancer contacts and the acquisition of DNA methylation. Msl2-knockout mice show perinatal lethality and heterogeneous phenotypes during embryonic development, supporting a role for MSL2 in regulating gene dosage. The role of MSL2 in preserving biallelic expression of specific dosage-sensitive genes sets the stage for further investigation of other factors that are involved in allelic dosage compensation in mammalian cells, with considerable implications for human disease.

Age- and caste-independent piRNAs in the germline and miRNA profiles linked to caste and fecundity in the ant Temnothorax rugatulus

Social insects are models for studies of phenotypic plasticity. Ant queens and workers vary in fecundity and lifespan, which are enhanced and extended in queens. Yet, the regulatory mechanisms underlying this variation are not well understood. Ant queens live and reproduce for years, so that they need to protect their germline from transposable element (TE) activity, which may be redundant in short-lived, often sterile workers. We analysed the expression of two protective classes of small RNAs, microRNAs (miRNAs) and Piwi-interacting RNAs (piRNAs), in various tissues, castes and age classes of the ant Temnothorax rugatulus. In queens, piRNAs were highly abundant in ovaries with TEs being their clear targets, with reduced but still detectable piRNA-specific ping-pong signatures in thorax and brains. piRNA pathway activity varied little with age in queens. Moreover, the reduced ovaries of workers also exhibited similar piRNA activity and this not only in young, fertile workers, but also in older foragers with regressed ovaries. Therefore, these ants protect their germline through piRNA activity, regardless of ovarian development, age or caste, even in sterile workers often considered the soma of the superorganism. Our tissue-specific miRNA analysis detected the expression of 304 miRNAs, of which 105 were expressed in all tissues, 10 enriched in the brain, three in the thorax, whereas 83 were ovarian-specific. We identified ovarian miRNAs whose expression was related to caste, fecundity and age, and which likely regulate group-specific gene expression. sRNA shifts in young- to middle-aged queens were minor, suggesting delayed senescence in this reproductive caste.

Structural Evolution of Gene Promoters Driven by Primate-Specific KRAB Zinc Finger Proteins

Krüppel-associated box (KRAB) zinc finger proteins (KZNFs) recognize and repress transposable elements (TEs); TEs are DNA elements that are capable of replicating themselves throughout our genomes with potentially harmful consequences. However, genes from this family of transcription factors have a much wider potential for genomic regulation. KZNFs have become integrated into gene-regulatory networks through the control of TEs that function as enhancers and gene promoters; some KZNFs also bind directly to gene promoters, suggesting an additional, more direct layer of KZNF co-option into gene-regulatory networks. Binding site analysis of ZNF519, ZNF441, and ZNF468 suggests the structural evolution of KZNFs to recognize TEs can result in coincidental binding to gene promoters independent of TE sequences. We show a higher rate of sequence turnover in gene promoter KZNF binding sites than neighboring regions, implying a selective pressure is being applied by the binding of a KZNF. Through CRISPR/Cas9 mediated genetic deletion of ZNF519, ZNF441, and ZNF468, we provide further evidence for genome-wide co-option of the KZNF-mediated gene-regulatory functions; KZNF knockout leads to changes in expression of KZNF-bound genes in neuronal lineages. Finally, we show that the opposite can be established upon KZNF overexpression, further strengthening the support for the role of KZNFs as bona-fide gene regulators. With no eminent role for ZNF519 in controlling its TE target, our study may provide a snapshot into the early stages of the completed co-option of a KZNF, showing the lasting, multilayered impact that retrovirus invasions and host response mechanisms can have upon the evolution of our genomes.

Isoform-resolved transcriptome of the human preimplantation embryo

Human preimplantation development involves extensive remodeling of RNA expression and splicing. However, its transcriptome has been compiled using short-read sequencing data, which fails to capture most full-length mRNAs. Here, we generate an isoform-resolved transcriptome of early human development by performing long- and short-read RNA sequencing on 73 embryos spanning the zygote to blastocyst stages. We identify 110,212 unannotated isoforms transcribed from known genes, including highly conserved protein-coding loci and key developmental regulators. We further identify 17,964 isoforms from 5,239 unannotated genes, which are largely non-coding, primate-specific, and highly associated with transposable elements. These isoforms are widely supported by the integration of published multi-omics datasets, including single-cell 8CLC and blastoid studies. Alternative splicing and gene co-expression network analyses further reveal that embryonic genome activation is associated with splicing disruption and transient upregulation of gene modules. Together, these findings show that the human embryo transcriptome is far more complex than currently known, and will act as a valuable resource to empower future studies exploring development.

The BAF chromatin remodeling complex licenses planarian stem cells access to ectodermal and mesodermal cell fates

BACKGROUND

The flatworm planarian, Schmidtea mediterranea, has a large population of adult stem cells (ASCs) that replace any cell type during tissue turnover or regeneration. How planarian ASCs (called neoblasts) manage self-renewal with the ability to produce daughter cells of different cell lineages (multipotency) is not well understood. Chromatin remodeling complexes ultimately control access to DNA regions of chromosomes and together with specific transcription factors determine whether a gene is transcribed in a given cell type. Previous work in planarians determined that RNAi of core components of the BAF chromatin remodeling complex, brg1 and smarcc2, caused increased ASCs and failed regeneration, but how these cellular defects arise at the level of gene regulation in neoblasts is unknown.

RESULTS

Here, we perform ATAC and RNA sequencing on purified neoblasts, deficient for the BAF complex subunits brg-1 and smarcc2. The data demonstrate that the BAF complex promotes chromatin accessibility and facilitates transcription at target loci, as in other systems. Interestingly, we find that the BAF complex enables access to genes known to be required for the generation of mesoderm- and ectoderm-derived lineages, including muscle, parenchymal cathepsin, neural, and epithelial lineages. BAF complex knockdowns result in disrupted differentiation into these cell lineages and functional consequences on planarian regeneration and tissue turnover. Notably, we did not detect a role for the BAF complex in neoblasts making endodermal lineages.

CONCLUSIONS

Our study provides functional insights into how the BAF complex contributes to cell fate decisions in planarian ASCs in vivo.

Palbociclib releases the latent differentiation capacity of neuroblastoma cells

Neuroblastoma is the most common extracranial solid tumor in infants, arising from developmentally stalled neural crest-derived cells. Driving tumor differentiation is a promising therapeutic approach for this devastating disease. Here, we show that the CDK4/6 inhibitor palbociclib not only inhibits proliferation but induces extensive neuronal differentiation of adrenergic neuroblastoma cells. Palbociclib-mediated differentiation is manifested by extensive phenotypic and transcriptional changes accompanied by the establishment of an epigenetic program driving expression of mature neuronal features. In vivo palbociclib significantly inhibits tumor growth in mouse neuroblastoma models. Furthermore, dual treatment with retinoic acid resets the oncogenic adrenergic core regulatory circuit of neuroblastoma cells, further suppresses proliferation, and can enhance differentiation, altering gene expression in ways that significantly correlate with improved patient survival. We therefore identify palbociclib as a therapeutic approach to dramatically enhance neuroblastoma differentiation efficacy that could be used in combination with retinoic acid to improve patient outcomes.

CRWN nuclear lamina components maintain the H3K27me3 landscape and promote successful reproduction in Arabidopsis

The nuclear lamina, a sub-nuclear protein matrix, maintains nuclear structure and genome function. Here, we investigate the role of Arabidopsis lamin analogs CROWDED NUCLEIs during gametophyte and seed development. We observed defects in crwn mutant seeds, including seed abortion and reduced germination rate. Quadruple crwn null genotypes were rarely transmitted through gametophytes. We focused on the crwn1 crwn2 (crwn1/2) endosperm, which exhibited enlarged chalazal cysts and increased expression of stress-related genes and the MADS-box transcription factor PHERES1 and its targets. Previously, it was shown that PHERES1 is regulated by H3K27me3 and that CRWN1 interacts with the PRC2 interactor PWO1. Thus, we tested whether crwn1/2 alters H3K27me3 patterns. We observed a mild loss of H3K27me3 at several hundred loci, which differed between endosperm and leaves. These data indicate that CRWNs are necessary to maintain the H3K27me3 landscape, with tissue-specific chromatin and transcriptional consequences.

Protocol for quantitative analysis of RNA 3'-end processing induced by disassociated subunits using chromatin-associated RNA-seq data

Sequencing chromatin-associated RNA using libraries from the chromatin fraction makes it possible to characterize RNA processing driven by disassociated subunits. Here, we present an experimental strategy and computational pipeline for processing chromatin-associated RNA-seq data to detect and quantify readthrough transcripts. We describe steps for constructing degron mouse embryonic stem cells, detecting readthrough genes, data processing, and data analysis. This protocol can be adapted to various biological scenarios and other types of nascent RNA-seq, such as TT-seq. For complete details on the use and execution of this protocol, please refer to Li et al. (2023).1.

Genome-wide mapping of i-motifs reveals their association with transcription regulation in live human cells

i-Motifs (iMs) are four-stranded DNA structures that form at cytosine (C)-rich sequences in acidic conditions in vitro. Their formation in cells is still under debate. We performed CUT&Tag sequencing using the anti-iM antibody iMab and showed that iMs form within the human genome in live cells. We mapped iMs in two human cell lines and recovered C-rich sequences that were confirmed to fold into iMs in vitro. We found that iMs in cells are mainly present at actively transcribing gene promoters, in open chromatin regions, they overlap with R-loops, and their abundance and distribution are specific to each cell type. iMs with both long and short C-tracts were recovered, further extending the relevance of iMs. By simultaneously mapping G-quadruplexes (G4s), which form at guanine-rich regions, and comparing the results with iMs, we proved that the two structures can form in independent regions; however, when both iMs and G4s are present in the same genomic tract, their formation is enhanced. iMs and G4s were mainly found at genes with low and high transcription rates, respectively. Our findings support the in vivo formation of iM structures and provide new insights into their interplay with G4s as new regulatory elements in the human genome.

KDM6A epigenetically regulates subtype plasticity in small cell lung cancer

Small cell lung cancer (SCLC) exists broadly in four molecular subtypes: ASCL1, NEUROD1, POU2F3 and Inflammatory. Initially, SCLC subtypes were thought to be mutually exclusive, but recent evidence shows intra-tumoural subtype heterogeneity and plasticity between subtypes. Here, using a CRISPR-based autochthonous SCLC genetically engineered mouse model to study the consequences of KDM6A/UTX inactivation, we show that KDM6A inactivation induced plasticity from ASCL1 to NEUROD1 resulting in SCLC tumours that express both ASCL1 and NEUROD1. Mechanistically, KDM6A normally maintains an active chromatin state that favours the ASCL1 subtype with its loss decreasing H3K4me1 and increasing H3K27me3 at enhancers of neuroendocrine genes leading to a cell state that is primed for ASCL1-to-NEUROD1 subtype switching. This work identifies KDM6A as an epigenetic regulator that controls ASCL1 to NEUROD1 subtype plasticity and provides an autochthonous SCLC genetically engineered mouse model to model ASCL1 and NEUROD1 subtype heterogeneity and plasticity, which is found in 35-40% of human SCLCs.

Cellular and Genomic Features of Muscle Differentiation from Isogenic Fibroblasts and Myoblasts

The ability to recapitulate muscle differentiation in vitro enables the exploration of mechanisms underlying myogenesis and muscle diseases. However, obtaining myoblasts from patients with neuromuscular diseases or from healthy subjects poses ethical and procedural challenges that limit such investigations. An alternative consists in converting skin fibroblasts into myogenic cells by forcing the expression of the myogenic regulator MYOD. Here, we directly compared cellular phenotype, transcriptome, and nuclear lamina-associated domains (LADs) in myo-converted human fibroblasts and myotubes differentiated from myoblasts. We used isogenic cells from a 16-year-old donor, ruling out, for the first time to our knowledge, genetic factors as a source of variations between the two myogenic models. We show that myo-conversion of fibroblasts upregulates genes controlling myogenic pathways leading to multinucleated cells expressing muscle cell markers. However, myotubes are more advanced in myogenesis than myo-converted fibroblasts at the phenotypic and transcriptomic levels. While most LADs are shared between the two cell types, each also displays unique domains of lamin A/C interactions. Furthermore, myotube-specific LADs are more gene-rich and less heterochromatic than shared LADs or LADs unique to myo-converted fibroblasts, and they uniquely sequester developmental genes. Thus, myo-converted fibroblasts and myotubes retain cell type-specific features of radial and functional genome organization. Our results favor a view of myo-converted fibroblasts as a practical model to investigate the phenotypic and genomic properties of muscle cell differentiation in normal and pathological contexts, but also highlight current limitations in using fibroblasts as a source of myogenic cells.

Mutant FUS induces chromatin reorganization in the hippocampus and alters memory processes

Cytoplasmic mislocalization of the nuclear Fused in Sarcoma (FUS) protein is associated to amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). Cytoplasmic FUS accumulation is recapitulated in the frontal cortex and spinal cord of heterozygous Fus∆NLS/+ mice. Yet, the mechanisms linking FUS mislocalization to hippocampal function and memory formation are still not characterized. Herein, we show that in these mice, the hippocampus paradoxically displays nuclear FUS accumulation. Multi-omic analyses showed that FUS binds to a set of genes characterized by the presence of an ETS/ELK-binding motifs, and involved in RNA metabolism, transcription, ribosome/mitochondria and chromatin organization. Importantly, hippocampal nuclei showed a decompaction of the neuronal chromatin at highly expressed genes and an inappropriate transcriptomic response was observed after spatial training of Fus∆NLS/+ mice. Furthermore, these mice lacked precision in a hippocampal-dependent spatial memory task and displayed decreased dendritic spine density. These studies shows that mutated FUS affects epigenetic regulation of the chromatin landscape in hippocampal neurons, which could participate in FTD/ALS pathogenic events. These data call for further investigation in the neurological phenotype of FUS-related diseases and open therapeutic strategies towards epigenetic drugs.

Predicting prime editing efficiency and product purity by deep learning

Prime editing is a versatile genome editing tool but requires experimental optimization of the prime editing guide RNA (pegRNA) to achieve high editing efficiency. Here we conducted a high-throughput screen to analyze prime editing outcomes of 92,423 pegRNAs on a highly diverse set of 13,349 human pathogenic mutations that include base substitutions, insertions and deletions. Based on this dataset, we identified sequence context features that influence prime editing and trained PRIDICT (prime editing guide prediction), an attention-based bidirectional recurrent neural network. PRIDICT reliably predicts editing rates for all small-sized genetic changes with a Spearman's R of 0.85 and 0.78 for intended and unintended edits, respectively. We validated PRIDICT on endogenous editing sites as well as an external dataset and showed that pegRNAs with high (>70) versus low (<70) PRIDICT scores showed substantially increased prime editing efficiencies in different cell types in vitro (12-fold) and in hepatocytes in vivo (tenfold), highlighting the value of PRIDICT for basic and for translational research applications.

The SWI/SNF chromatin remodeling complexes BAF and PBAF differentially regulate epigenetic transitions in exhausted CD8+ T cells

CD8⁺ T cell exhaustion (Tex) limits disease control during chronic viral infections and cancer. Here, we investigated the epigenetic factors mediating major chromatin-remodeling events in Tex-cell development. A protein-domain-focused in vivo CRISPR screen identified distinct functions for two versions of the SWI/SNF chromatin-remodeling complex in Tex-cell differentiation. Depletion of the canonical SWI/SNF form, BAF, impaired initial CD8⁺ T cell responses in acute and chronic infection. In contrast, disruption of PBAF enhanced Tex-cell proliferation and survival. Mechanistically, PBAF regulated the epigenetic and transcriptional transition from TCF-1⁺ progenitor Tex cells to more differentiated TCF-1^- Tex subsets. Whereas PBAF acted to preserve Tex progenitor biology, BAF was required to generate effector-like Tex cells, suggesting that the balance of these factors coordinates Tex-cell subset differentiation. Targeting PBAF improved tumor control both alone and in combination with anti-PD-L1 immunotherapy. Thus, PBAF may present a therapeutic target in cancer immunotherapy.

EZH2 inhibition remodels the inflammatory senescence-associated secretory phenotype to potentiate pancreatic cancer immune surveillance

Immunotherapies that produce durable responses in some malignancies have failed in pancreatic ductal adenocarcinoma (PDAC) due to rampant immune suppression and poor tumor immunogenicity. We and others have demonstrated that induction of the senescence-associated secretory phenotype (SASP) can be an effective approach to activate anti-tumor natural killer (NK) cell and T cell immunity. In the present study, we found that the pancreas tumor microenvironment suppresses NK cell and T cell surveillance after therapy-induced senescence through enhancer of zeste homolog 2 (EZH2)-mediated epigenetic repression of proinflammatory SASP genes. EZH2 blockade stimulated production of SASP chemokines CCL2 and CXCL9/10, leading to enhanced NK cell and T cell infiltration and PDAC eradication in mouse models. EZH2 activity was also associated with suppression of chemokine signaling and cytotoxic lymphocytes and reduced survival in patients with PDAC. These results demonstrate that EZH2 represses the proinflammatory SASP and that EZH2 inhibition combined with senescence-inducing therapy could be a powerful means to achieve immune-mediated tumor control in PDAC.

Transcriptome sequencing suggests that pre-mRNA splicing counteracts widespread intronic cleavage and polyadenylation

Alternative splicing (AS) and alternative polyadenylation (APA) are two crucial steps in the post-transcriptional regulation of eukaryotic gene expression. Protocols capturing and sequencing RNA 3'-ends have uncovered widespread intronic polyadenylation (IPA) in normal and disease conditions, where it is currently attributed to stochastic variations in the pre-mRNA processing. Here, we took advantage of the massive amount of RNA-seq data generated by the Genotype Tissue Expression project (GTEx) to simultaneously identify and match tissue-specific expression of intronic polyadenylation sites with tissue-specific splicing. A combination of computational methods including the analysis of short reads with non-templated adenines revealed that APA events are more abundant in introns than in exons. While the rate of IPA in composite terminal exons and skipped terminal exons expectedly correlates with splicing, we observed a considerable fraction of IPA events that lack AS support and attributed them to spliced polyadenylated introns (SPI). We hypothesize that SPIs represent transient byproducts of a dynamic coupling between APA and AS, in which the spliceosome removes the intron while it is being cleaved and polyadenylated. These findings indicate that cotranscriptional pre-mRNA splicing could serve as a rescue mechanism to suppress premature transcription termination at intronic polyadenylation sites.

Optimal HSF1 activation in response to acute cold stress in BAT requires nuclear TXNIP

While TXNIP (thioredoxin interacting protein) in the plasma membrane and vesicular location is known to negatively regulate cellular glucose uptake by facilitating glucose transporter endocytosis, the function of TXNIP in the nucleus is far less understood. Herein, we sought to determine the function of nuclear TXNIP in vivo, using a new HA-tagged TXNIP knock-in mouse model. We observed that TXNIP can be found in the nucleus of a variety of cells from different tissues including hepatocytes (liver), enterocytes (small intestine), exocrine cells (pancreas), and brown adipocytes (BAT). Further investigations into the role of nuclear TXNIP in BAT revealed that cold stress rapidly and transiently activated HSF1 (heat shock factor 1). HSF1 interaction with TXNIP during its activation is required for optimal HSF1 directed cold shock response in BAT.

Transition from Transient DNA Rereplication to Inherited Gene Amplification Following Prolonged Environmental Stress

Cells require the ability to adapt to changing environmental conditions, however, it is unclear how these changes elicit stable permanent changes in genomes. We demonstrate that, in response to environmental metal exposure, the metallothionein (MT) locus undergoes DNA rereplication generating transient site-specific gene amplifications (TSSGs). Chronic metal exposure allows transition from MT TSSG to inherited MT gene amplification through homologous recombination within and outside of the MT locus. DNA rereplication of the MT locus is suppressed by H3K27me3 and EZH2. Long-term ablation of EZH2 activity eventually leads to integration and inheritance of MT gene amplifications without the selective pressure of metal exposure. The rereplication and inheritance of MT gene amplification is an evolutionarily conserved response to environmental metal from yeast to human. Our results describe a new paradigm for adaptation to environmental stress where targeted, transient DNA rereplication precedes stable inherited gene amplification.

RNA Pol II preferentially regulates ribosomal protein expression by trapping disassociated subunits

RNA polymerase II (RNA Pol II) has been recognized as a passively regulated multi-subunit holoenzyme. However, the extent to which RNA Pol II subunits might be important beyond the RNA Pol II complex remains unclear. Here, fractions containing disassociated RPB3 (dRPB3) were identified by size exclusion chromatography in various cells. Through a unique strategy, i.e., "specific degradation of disassociated subunits (SDDS)," we demonstrated that dRPB3 functions as a regulatory component of RNA Pol II to enable the preferential control of 3' end processing of ribosomal protein genes directly through its N-terminal domain. Machine learning analysis of large-scale genomic features revealed that the little elongation complex (LEC) helps to specialize the functions of dRPB3. Mechanistically, dRPB3 facilitates CBC-PCF11 axis activity to increase the efficiency of 3' end processing. Furthermore, RPB3 is dynamically regulated during development and diseases. These findings suggest that RNA Pol II gains specific regulatory functions by trapping disassociated subunits in mammalian cells.

High level of novelty under the hood of convergent evolution

Little is known about the extent to which species use homologous regulatory architectures to achieve phenotypic convergence. By characterizing chromatin accessibility and gene expression in developing wing tissues, we compared the regulatory architecture of convergence between a pair of mimetic butterfly species. Although a handful of color pattern genes are known to be involved in their convergence, our data suggest that different mutational paths underlie the integration of these genes into wing pattern development. This is supported by a large fraction of accessible chromatin being exclusive to each species, including the de novo lineage-specific evolution of a modular optix enhancer. These findings may be explained by a high level of developmental drift and evolutionary contingency that occurs during the independent evolution of mimicry.

Single-strand mismatch and damage patterns revealed by single-molecule DNA sequencing

Mutations accumulate in the genome of every cell of the body throughout life, causing cancer and other genetic diseases1-4. Almost all of these mosaic mutations begin as nucleotide mismatches or damage in only one of the two strands of the DNA prior to becoming double-strand mutations if unrepaired or misrepaired5. However, current DNA sequencing technologies cannot resolve these initial single-strand events. Here, we developed a single-molecule, long-read sequencing method that achieves single-molecule fidelity for single-base substitutions when present in either one or both strands of the DNA. It also detects single-strand cytosine deamination events, a common type of DNA damage. We profiled 110 samples from diverse tissues, including from individuals with cancer-predisposition syndromes, and define the first single-strand mismatch and damage signatures. We find correspondences between these single-strand signatures and known double-strand mutational signatures, which resolves the identity of the initiating lesions. Tumors deficient in both mismatch repair and replicative polymerase proofreading show distinct single-strand mismatch patterns compared to samples deficient in only polymerase proofreading. In the mitochondrial genome, our findings support a mutagenic mechanism occurring primarily during replication. Since the double-strand DNA mutations interrogated by prior studies are only the endpoint of the mutation process, our approach to detect the initiating single-strand events at single-molecule resolution will enable new studies of how mutations arise in a variety of contexts, especially in cancer and aging.

Gene Regulatory Interactions at Lamina-Associated Domains

The nuclear lamina provides a repressive chromatin environment at the nuclear periphery. However, whereas most genes in lamina-associated domains (LADs) are inactive, over ten percent reside in local euchromatic contexts and are expressed. How these genes are regulated and whether they are able to interact with regulatory elements remain unclear. Here, we integrate publicly available enhancer-capture Hi-C data with our own chromatin state and transcriptomic datasets to show that inferred enhancers of active genes in LADs are able to form connections with other enhancers within LADs and outside LADs. Fluorescence in situ hybridization analyses show proximity changes between differentially expressed genes in LADs and distant enhancers upon the induction of adipogenic differentiation. We also provide evidence of involvement of lamin A/C, but not lamin B1, in repressing genes at the border of an in-LAD active region within a topological domain. Our data favor a model where the spatial topology of chromatin at the nuclear lamina is compatible with gene expression in this dynamic nuclear compartment.

Cell type-specific DNA methylome signatures reveal epigenetic mechanisms for neuronal diversity and neurodevelopmental disorder

DNA methylation plays a critical function in establishing and maintaining cell identity in brain. Disruption of DNA methylation-related processes leads to diverse neurological disorders. However, the role of DNA methylation characteristics in neuronal diversity remains underexplored. Here, we report detailed context-specific DNA methylation maps for GABAergic, glutamatergic (Glu) and Purkinje neurons, together with matched transcriptome profiles. Genome-wide mCH levels are distinguishable, while the mCG levels are similar among the three cell types. Substantial CG-differentially methylated regions (DMRs) are also seen, with Glu neurons experiencing substantial hypomethylation events. The relationship between mCG levels and gene expression displays cell type-specific patterns, while genic CH methylation exhibits a negative effect on transcriptional abundance. We found that cell type-specific CG-DMRs are informative in terms of represented neuronal function. Furthermore, we observed that the identified Glu-specific hypo-DMRs have a high level of consistency with the chromatin accessibility of excitatory neurons and the regions enriched for histone modifications (H3K27ac and H3K4me1) of active enhancers, suggesting their regulatory potential. Hypomethylation regions specific to each cell type are predicted to bind neuron type-specific transcription factors. Finally, we show that the DNA methylation changes in a mouse model of Rett syndrome, a neurodevelopmental disorder caused by the de novo mutations in MECP2, are cell type- and brain region-specific. Our results suggest that cell type-specific DNA methylation signatures are associated with the functional characteristics of the neuronal subtypes. The presented results emphasize the importance of DNA methylation-mediated epigenetic regulation in neuronal diversity and disease.

Spurious transcription causing innate immune responses is prevented by 5-hydroxymethylcytosine

Generation of functional transcripts requires transcriptional initiation at regular start sites, avoiding production of aberrant and potentially hazardous aberrant RNAs. The mechanisms maintaining transcriptional fidelity and the impact of spurious transcripts on cellular physiology and organ function have not been fully elucidated. Here we show that TET3, which successively oxidizes 5-methylcytosine to 5-hydroxymethylcytosine (5hmC) and other derivatives, prevents aberrant intragenic entry of RNA polymerase II pSer5 into highly expressed genes of airway smooth muscle cells, assuring faithful transcriptional initiation at canonical start sites. Loss of TET3-dependent 5hmC production in SMCs results in accumulation of spurious transcripts, which stimulate the endosomal nucleic-acid-sensing TLR7/8 signaling pathway, thereby provoking massive inflammation and airway remodeling resembling human bronchial asthma. Furthermore, we found that 5hmC levels are substantially lower in human asthma airways compared with control samples. Suppression of spurious transcription might be important to prevent chronic inflammation in asthma.

Circulating cell-free DNA methylation mirrors alterations in cerebral patterns in epilepsy

BACKGROUND

DNA methylation profiling of circulating cell-free DNA (cfDNA) has rapidly become a promising strategy for biomarker identification and development. The cell-type-specific nature of DNA methylation patterns and the direct relationship between cfDNA and apoptosis can potentially be used non-invasively to predict local alterations. In addition, direct detection of altered DNA methylation patterns performs well as a biomarker. In a previous study, we demonstrated marked DNA methylation alterations in brain tissue from patients with mesial temporal lobe epilepsy with hippocampal sclerosis (MTLE-HS).

RESULTS

We performed DNA methylation profiling in cfDNA isolated from the serum of MTLE patients and healthy controls using BeadChip arrays followed by systematic bioinformatic analysis including deconvolution analysis and integration with DNase accessibility data sets. Differential cfDNA methylation analysis showed an overrepresentation of gene ontology terms and transcription factors related to central nervous system function and regulation. Deconvolution analysis of the DNA methylation data sets ruled out the possibility that the observed differences were due to changes in the proportional contribution of cortical neurons in cfDNA. Moreover, we found no overrepresentation of neuron- or glia-specific patterns in the described cfDNA methylation patterns. However, the MTLE-HS cfDNA methylation patterns featured a significant overrepresentation of the epileptic DNA methylation alterations previously observed in the hippocampus.

CONCLUSIONS

Our results support the use of cfDNA methylation profiling as a rational approach to seeking non-invasive and reproducible epilepsy biomarkers.

PPARγ and C/EBPα response to acute cold stress in brown adipose tissue

Brown adipose tissue (BAT) has the ability to burn calories as heat. Utilizing BAT thermogenesis is thus an attractive way to combat obesity. However, the transcriptional network resulting in the lipid synthesis to oxidation shift during thermogenesis is not completely understood. Here, we report the regulation of two master regulators of adipogenesis, peroxisome proliferator-activated receptor gamma (PPARγ) and CCAAT/enhancer-binding protein alpha (C/EBPα), during acute cold stress in BAT. We found PPARγ dissociates from DNA in a fifth of its binding sites and these include Cebpa enhancers, leading to decreased C/EBPα expression. This dissociation requires PPARγ binding to activating ligands and is thus modulated by diet. Meanwhile, PPARα also detaches from DNA, and co-activator PGC1α associates with ERRα as part of a transcriptional network regulating lipid metabolism. Subsequent global replacement of C/EBPα by C/EBPβ and its associated transcriptional machinery is required for upregulation of structural lipid synthesis despite general upregulation of fatty acid oxidation.

Integration of 3D genome topology and local chromatin features uncovers enhancers underlying craniofacial-specific cartilage defects

Aberrations in tissue-specific enhancers underlie many developmental defects. Disrupting a noncoding region distal from the human SOX9 gene causes the Pierre Robin sequence (PRS) characterized by the undersized lower jaw. Such a craniofacial-specific defect has been previously linked to enhancers transiently active in cranial neural crest cells (CNCCs). We demonstrate that the PRS region also strongly regulates Sox9 in CNCC-derived Meckel's cartilage (MC), but not in limb cartilages, even after decommissioning of CNCC enhancers. Such an MC-specific regulatory effect correlates with the MC-specific chromatin contacts between the PRS region and Sox9, highlighting the importance of lineage-dependent chromatin topology in instructing enhancer usage. By integrating the enhancer signatures and chromatin topology, we uncovered >10,000 enhancers that function differentially between MC and limb cartilages and demonstrated their association with human diseases. Our findings provide critical insights for understanding the choreography of gene regulation during development and interpreting the genetic basis of craniofacial pathologies.

EBF1 primes B-lymphoid enhancers and limits the myeloid bias in murine multipotent progenitors

Hematopoietic stem cells (HSCs) and multipotent progenitors (MPPs) generate all cells of the blood system. Despite their multipotency, MPPs display poorly understood lineage bias. Here, we examine whether lineage-specifying transcription factors, such as the B-lineage determinant EBF1, regulate lineage preference in early progenitors. We detect low-level EBF1 expression in myeloid-biased MPP3 and lymphoid-biased MPP4 cells, coinciding with expression of the myeloid determinant C/EBPα. Hematopoietic deletion of Ebf1 results in enhanced myelopoiesis and reduced HSC repopulation capacity. Ebf1-deficient MPP3 and MPP4 cells exhibit an augmented myeloid differentiation potential and a transcriptome with an enriched C/EBPα signature. Correspondingly, EBF1 binds the Cebpa enhancer, and the deficiency and overexpression of Ebf1 in MPP3 and MPP4 cells lead to an up- and downregulation of Cebpa expression, respectively. In addition, EBF1 primes the chromatin of B-lymphoid enhancers specifically in MPP3 cells. Thus, our study implicates EBF1 in regulating myeloid/lymphoid fate bias in MPPs by constraining C/EBPα-driven myelopoiesis and priming the B-lymphoid fate.

Vitamin C enhances NF-κB-driven epigenomic reprogramming and boosts the immunogenic properties of dendritic cells

Dendritic cells (DCs), the most potent antigen-presenting cells, are necessary for effective activation of naïve T cells. DCs’ immunological properties are modulated in response to various stimuli. Active DNA demethylation is crucial for DC differentiation and function. Vitamin C, a known cofactor of ten-eleven translocation (TET) enzymes, drives active demethylation. Vitamin C has recently emerged as a promising adjuvant for several types of cancer; however, its effects on human immune cells are poorly understood. In this study, we investigate the epigenomic and transcriptomic reprogramming orchestrated by vitamin C in monocyte-derived DC differentiation and maturation. Vitamin C triggers extensive demethylation at NF-κB/p65 binding sites, together with concordant upregulation of antigen-presentation and immune response-related genes during DC maturation. p65 interacts with TET2 and mediates the aforementioned vitamin C-mediated changes, as demonstrated by pharmacological inhibition. Moreover, vitamin C increases TNFβ production in DCs through NF-κB, in concordance with the upregulation of its coding gene and the demethylation of adjacent CpGs. Finally, vitamin C enhances DC’s ability to stimulate the proliferation of autologous antigen-specific T cells. We propose that vitamin C could potentially improve monocyte-derived DC-based cell therapies.

The TRIPLE PHD FINGERS proteins are required for SWI/SNF complex-mediated +1 nucleosome positioning and transcription start site determination in Arabidopsis

Eukaryotes have evolved multiple ATP-dependent chromatin remodelers to shape the nucleosome landscape. We recently uncovered an evolutionarily conserved SWItch/Sucrose Non-Fermentable (SWI/SNF) chromatin remodeler complex in plants reminiscent of the mammalian BAF subclass, which specifically incorporates the MINUSCULE (MINU) catalytic subunits and the TRIPLE PHD FINGERS (TPF) signature subunits. Here we report experimental evidence that establishes the functional relevance of TPF proteins for the complex activity. Our results show that depletion of TPF triggers similar pleiotropic phenotypes and molecular defects to those found in minu mutants. Moreover, we report the genomic location of MINU2 and TPF proteins as representative members of this SWI/SNF complex and their impact on nucleosome positioning and transcription. These analyses unravel the binding of the complex to thousands of genes where it modulates the position of the +1 nucleosome. These targets tend to produce 5′-shifted transcripts in the tpf and minu mutants pointing to the participation of the complex in alternative transcription start site usage. Interestingly, there is a remarkable correlation between +1 nucleosome shift and 5′ transcript length change suggesting their functional connection. In summary, this study unravels the function of a plant SWI/SNF complex involved in +1 nucleosome positioning and transcription start site determination.

A butterfly pan-genome reveals that a large amount of structural variation underlies the evolution of chromatin accessibility

Despite insertions and deletions being the most common structural variants (SVs) found across genomes, not much is known about how much these SVs vary within populations and between closely related species, nor their significance in evolution. To address these questions, we characterized the evolution of indel SVs using genome assemblies of three closely related Heliconius butterfly species. Over the relatively short evolutionary timescales investigated, up to 18.0% of the genome was composed of indels between two haplotypes of an individual Heliconius charithonia butterfly and up to 62.7% included lineage-specific SVs between the genomes of the most distant species (11 Mya). Lineage-specific sequences were mostly characterized as transposable elements (TEs) inserted at random throughout the genome and their overall distribution was similarly affected by linked selection as single nucleotide substitutions. Using chromatin accessibility profiles (i.e., ATAC-seq) of head tissue in caterpillars to identify sequences with potential cis-regulatory function, we found that out of the 31,066 identified differences in chromatin accessibility between species, 30.4% were within lineage-specific SVs and 9.4% were characterized as TE insertions. These TE insertions were localized closer to gene transcription start sites than expected at random and were enriched for sites with significant resemblance to several transcription factor binding sites with known function in neuron development in Drosophila We also identified 24 TE insertions with head-specific chromatin accessibility. Our results show high rates of structural genome evolution that were previously overlooked in comparative genomic studies and suggest a high potential for structural variation to serve as raw material for adaptive evolution.

Progression of prostate cancer reprograms MYC-mediated lipid metabolism via lysine methyltransferase 2A

BACKGROUND

The activities of MYC, the androgen receptor, and its associated pioneer factors demonstrate substantial reprogramming between early and advanced prostate cancer. Although previous studies have shown a shift in cellular metabolic requirements associated with prostate cancer progression, the epigenetic regulation of these processes is incompletely described. Here, we have integrated chromatin immunoprecipitation sequencing (ChIP-seq) and whole-transcriptome sequencing to identify novel regulators of metabolism in advanced prostate tumors characterized by elevated MYC activity.

RESULTS

Using ChIP-seq against MYC, HOXB13, and AR in LNCaP cells, we observed redistribution of co-bound sites suggestive of differential KMT2A activity as a function of MYC expression. In a cohort of 177 laser-capture microdissected foci of prostate tumors, KMT2A expression was positively correlated with MYC activity, AR activity, and HOXB13 expression, but decreased with tumor grade severity. However, KMT2A expression was negatively correlated with these factors in 25 LuCaP patient-derived xenograft models of advanced prostate cancer and 99 laser-capture microdissected foci of metastatic castration-resistant prostate cancer. Stratified by KMT2A expression, ChIP-seq against AR and HOXB13 in 15 LuCaP patient-derived xenografts showed an inverse association with sites involving genes implicated in lipid metabolism, including the arachidonic acid metabolic enzyme PLA2G4F. LuCaP patient-derived xenograft models grown as organoids recapitulated the inverse association between KMT2A expression and fluorine-18 labeled arachidonic acid uptake in vitro.

CONCLUSIONS

Our study demonstrates that the epigenetic activity of transcription factor oncogenes exhibits a shift during prostate cancer progression with distinctive phenotypic effects on metabolism. These epigenetically driven changes in lipid metabolism may serve as novel targets for the development of novel imaging agents and therapeutics.

Targeted protein degradation reveals RNA Pol II heterogeneity and functional diversity

RNA polymerase II (RNA Pol II) subunits are thought to be involved in various transcription-associated processes, but it is unclear whether they play different regulatory roles in modulating gene expression. Here, we performed nascent and mature transcript sequencing after the acute degradation of 12 mammalian RNA Pol II subunits and profiled their genomic binding sites and protein interactomes to dissect their molecular functions. We found that RNA Pol II subunits contribute differently to RNA Pol II cellular localization and transcription processes and preferentially regulate RNA processing (such as RNA splicing and 3' end maturation). Genes sensitive to the depletion of different RNA Pol II subunits tend to be involved in diverse biological functions and show different RNA half-lives. Sequences, associated protein factors, and RNA structures are correlated with RNA Pol II subunit-mediated differential gene expression. These findings collectively suggest that the heterogeneity of RNA Pol II and different genes appear to depend on some of the subunits.

Identification of putative enhancer-like elements predicts regulatory networks active in planarian adult stem cells

Planarians have become an established model system to study regeneration and stem cells, but the regulatory elements in the genome remain almost entirely undescribed. Here, by integrating epigenetic and expression data we use multiple sources of evidence to predict enhancer elements active in the adult stem cell populations that drive regeneration. We have used ChIP-seq data to identify genomic regions with histone modifications consistent with enhancer activity, and ATAC-seq data to identify accessible chromatin. Overlapping these signals allowed for the identification of a set of high-confidence candidate enhancers predicted to be active in planarian adult stem cells. These enhancers are enriched for predicted transcription factor (TF) binding sites for TFs and TF families expressed in planarian adult stem cells. Footprinting analyses provided further evidence that these potential TF binding sites are likely to be occupied in adult stem cells. We integrated these analyses to build testable hypotheses for the regulatory function of TFs in stem cells, both with respect to how pluripotency might be regulated, and to how lineage differentiation programs are controlled. We found that our predicted GRNs were independently supported by existing TF RNAi/RNA-seq datasets, providing further evidence that our work predicts active enhancers that regulate adult stem cells and regenerative mechanisms.

Maternal H3K36 and H3K27 HMTs protect germline development via regulation of the transcription factor LIN-15B

Maternally synthesized products play critical roles in the development of offspring. A premier example is the Caenorhabditis elegans H3K36 methyltransferase MES-4, which is essential for germline survival and development in offspring. How maternal MES-4 protects the germline is not well understood, but its role in H3K36 methylation hinted that it may regulate gene expression in primordial germ cells (PGCs). We tested this hypothesis by profiling transcripts from nascent germlines (PGCs and their descendants) dissected from wild-type and mes-4 mutant (lacking maternal and zygotic MES-4) larvae. mes-4 nascent germlines displayed downregulation of some germline genes, upregulation of some somatic genes, and dramatic upregulation of hundreds of genes on the X chromosome. We demonstrated that upregulation of one or more genes on the X is the cause of germline death by generating and analyzing mes-4 mutants that inherited different endowments of X chromosome(s). Intriguingly, removal of the THAP transcription factor LIN-15B from mes-4 mutants reduced X misexpression and prevented germline death. lin-15B is X-linked and misexpressed in mes-4 PGCs, identifying it as a critical target for MES-4 repression. The above findings extend to the H3K27 methyltransferase MES-2/3/6, the C. elegans version of polycomb repressive complex 2. We propose that maternal MES-4 and PRC2 cooperate to protect germline survival by preventing synthesis of germline-toxic products encoded by genes on the X chromosome, including the key transcription factor LIN-15B.

High enhancer activity is an epigenetic feature of HPV negative atypical head and neck squamous cell carcinoma

Head and neck squamous cell carcinoma (HNSCC) is a heterogeneous disease with significant mortality and frequent recurrence. Prior efforts to transcriptionally classify HNSCC into groups of varying prognoses have identified four accepted molecular subtypes of the disease: Atypical (AT), Basal (BA), Classical (CL), and Mesenchymal (MS). Here, we investigate the active enhancer landscapes of these subtypes using representative HNSCC cell lines and identify samples belonging to the AT subtype as having increased enhancer activity compared to the other 3 HNSCC subtypes. Cell lines belonging to the AT subtype are more resistant to enhancer-blocking bromodomain inhibitors (BETi). Examination of nascent transcripts reveals that both AT TCGA tumors and cell lines express higher levels of enhancer RNA (eRNA) transcripts for enhancers controlling BETi resistance pathways, such as lipid metabolism and MAPK signaling. Additionally, investigation of higher-order chromatin structure suggests more enhancer-promoter (E-P) contacts in the AT subtype, including on genes identified in the eRNA analysis. Consistently, known BETi resistance pathways are upregulated upon exposure to these inhibitors. Together, our results identify that the AT subtype of HNSCC is associated with higher enhancer activity, resistance to enhancer blockade, and increased signaling through pathways that could serve as future targets for sensitizing HNSCC to BET inhibition.