BET inhibition reforms the immune microenvironment and alleviates T cell dysfunction in chronic lymphocytic leukemia.

Redundant tumor microenvironment (TME) immunosuppressive mechanisms and epigenetic maintenance of terminal T cell exhaustion greatly hinder functional antitumor immune responses in chronic lymphocytic leukemia (CLL). Bromodomain and extraterminal (BET) proteins regulate key pathways contributing to CLL pathogenesis and TME interactions, including T cell function and differentiation. Herein, we report that blocking BET protein function alleviates immunosuppressive networks in the CLL TME and repairs inherent CLL T cell defects. The pan-BET inhibitor OPN-51107 reduced exhaustion-associated cell signatures resulting in improved T cell proliferation and effector function in the Eµ-TCL1 splenic TME. Following BET inhibition (BET-i), TME T cells coexpressed significantly fewer inhibitory receptors (IRs) (e.g., PD-1, CD160, CD244, LAG3, VISTA). Complementary results were witnessed in primary CLL cultures, wherein OPN-51107 exerted proinflammatory effects on T cells, regardless of leukemic cell burden. BET-i additionally promotes a progenitor T cell phenotype through reduced expression of transcription factors that maintain terminal differentiation and increased expression of TCF-1, at least in part through altered chromatin accessibility. Moreover, direct T cell effects of BET-i were unmatched by common targeted therapies in CLL. This study demonstrates the immunomodulatory action of BET-i on CLL T cells and supports the inclusion of BET inhibitors in the management of CLL to alleviate terminal T cell dysfunction and potentially enhance tumoricidal T cell activity.

Novel Spirocyclic Dimer, SpiD3, Targets Chronic Lymphocytic Leukemia Survival Pathways with Potent Preclinical Effects.

Chronic lymphocytic leukemia (CLL) cell survival and growth is fueled by the induction of B-cell receptor (BCR) signaling within the tumor microenvironment (TME) driving activation of NFκB signaling and the unfolded protein response (UPR). Malignant cells have higher basal levels of UPR posing a unique therapeutic window to combat CLL cell growth using pharmacologic agents that induce accumulation of misfolded proteins. Frontline CLL therapeutics that directly target BCR signaling such as Bruton tyrosine kinase (BTK) inhibitors (e.g., ibrutinib) have enhanced patient survival. However, resistance mechanisms wherein tumor cells bypass BTK inhibition through acquired BTK mutations, and/or activation of alternative survival mechanisms have rendered ibrutinib ineffective, imposing the need for novel therapeutics. We evaluated SpiD3, a novel spirocyclic dimer, in CLL cell lines, patient-derived CLL samples, ibrutinib-resistant CLL cells, and in the Eµ-TCL1 mouse model. Our integrated multi-omics and functional analyses revealed BCR signaling, NFκB signaling, and endoplasmic reticulum stress among the top pathways modulated by SpiD3. This was accompanied by marked upregulation of the UPR and inhibition of global protein synthesis in CLL cell lines and patient-derived CLL cells. In ibrutinib-resistant CLL cells, SpiD3 retained its antileukemic effects, mirrored in reduced activation of key proliferative pathways (e.g., PRAS, ERK, MYC). Translationally, we observed reduced tumor burden in SpiD3-treated Eµ-TCL1 mice. Our findings reveal that SpiD3 exploits critical vulnerabilities in CLL cells including NFκB signaling and the UPR, culminating in profound antitumor properties independent of TME stimuli. SIGNIFICANCE: SpiD3 demonstrates cytotoxicity in CLL partially through inhibition of NFκB signaling independent of tumor-supportive stimuli. By inducing the accumulation of unfolded proteins, SpiD3 activates the UPR and hinders protein synthesis in CLL cells. Overall, SpiD3 exploits critical CLL vulnerabilities (i.e., the NFκB pathway and UPR) highlighting its use in drug-resistant CLL.

Predicting A/B compartments from histone modifications using deep learning.

The three-dimensional organization of genomes plays a crucial role in essential biological processes. The segregation of chromatin into A and B compartments highlights regions of activity and inactivity, providing a window into the genomic activities specific to each cell type. Yet, the steep costs associated with acquiring Hi-C data, necessary for studying this compartmentalization across various cell types, pose a significant barrier in studying cell type specific genome organization. To address this, we present a prediction tool called compartment prediction using recurrent neural networks (CoRNN), which predicts compartmentalization of 3D genome using histone modification enrichment. CoRNN demonstrates robust cross-cell-type prediction of A/B compartments with an average AuROC of 90.9%. Cell-type-specific predictions align well with known functional elements, with H3K27ac and H3K36me3 identified as highly predictive histone marks. We further investigate our mispredictions and found that they are located in regions with ambiguous compartmental status. Furthermore, our model's generalizability is validated by predicting compartments in independent tissue samples, which underscores its broad applicability.

Mechanistic drivers of chromatin organization into compartments.

The human genome is not just a simple string of DNA, it is a complex and dynamic entity intricately folded within the cell's nucleus. This three-dimensional organization of chromatin, the combination of DNA and proteins in the nucleus, is crucial for many biological processes and has been prominently studied for its intricate relationship to gene expression. Indeed, the transcriptional machinery does not operate in isolation but interacts intimately with the folded chromatin structure. Techniques for chromatin conformation capture, including genome-wide sequencing approaches, have revealed key organizational features of chromatin, such as the formation of loops by CCCTC-binding factor (CTCF) and the division of loci into chromatin compartments. While much of the recent research and reviews have focused on CTCF loops, we discuss several new revelations that have emerged concerning chromatin compartments, with a particular focus on what is known about mechanistic drivers of compartmentalization. These insights challenge the traditional views of chromatin organization and reveal the complexity behind the formation and maintenance of chromatin compartments.

An environmentally relevant mixture of Perfluoroalkyl Substances (PFAS) impacts proliferation, steroid hormone synthesis, and gene transcription in primary human granulosa cells.

Perfluoroalkyl substances (PFAS) are a group of synthetic chemicals that are resistant to biodegradation and are environmentally persistent. PFAS are found in many consumer products and are a major source of water and soil contamination. This study investigated the effects of an environmentally relevant PFAS mixture [perfluorooctanoic acid (PFOA), perfluorooctanesulfonic acid (PFOS), perfluorohexanesulfonic acid (PFHxS)] on the transcriptome and function of human granulosa cells (hGCs). Primary hGCs were harvested from follicular aspirates of healthy, reproductive-age women who were undergoing oocyte retrieval for in vitro fertilization. LC/MS-MS was performed to identify PFAS compounds in pure follicular fluid. Cells were cultured with vehicle control or a PFAS mixture (2 nM PFHxS, 7 nM PFOA, 10 nM PFOS) for 96h. Analyses of cell proliferation/apoptosis, steroidogenesis, and gene expression were measured via MTT assays/immunofluorescence, ELISA/western blotting, and RNA sequencing/bioinformatics, respectively. PFOA, PFOS, and PFHxS were detected in 100% of follicle fluid samples. Increased cell proliferation was observed in hGCs treated with the PFAS mixture with no impacts on cellular apoptosis. The PFAS mixture also altered steroid hormone synthesis, increasing both FSH-stimulated and basal progesterone secretion and concomitant upregulation of STAR protein. RNA sequencing revealed inherent differences in transcriptomic profiles in hGCs after PFAS exposure. This study demonstrates functional and transcriptomic changes in hGCs after exposure to a PFAS mixture, improving our knowledge about the impacts of PFAS exposures and female reproductive health. These findings suggest that PFAS compounds can disrupt normal granulosa cell function with possible long-term consequences on overall reproductive health.

The zinc-finger protein Z4 cooperates with condensin II to regulate somatic chromosome pairing and 3D chromatin organization.

Chromosome pairing constitutes an important level of genome organization, yet the mechanisms that regulate pairing in somatic cells and the impact on 3D chromatin organization are still poorly understood. Here, we address these questions in Drosophila, an organism with robust somatic pairing. In Drosophila, pairing preferentially occurs at loci consisting of numerous architectural protein binding sites (APBSs), suggesting a role of architectural proteins (APs) in pairing regulation. Amongst these, the anti-pairing function of the condensin II subunit CAP-H2 is well established. However, the factors that regulate CAP-H2 localization and action at APBSs remain largely unknown. Here, we identify two factors that control CAP-H2 occupancy at APBSs and, therefore, regulate pairing. We show that Z4, interacts with CAP-H2 and is required for its localization at APBSs. We also show that hyperosmotic cellular stress induces fast and reversible unpairing in a Z4/CAP-H2 dependent manner. Moreover, by combining the opposite effects of Z4 depletion and osmostress, we show that pairing correlates with the strength of intrachromosomal 3D interactions, such as active (A) compartment interactions, intragenic gene-loops, and polycomb (Pc)-mediated chromatin loops. Altogether, our results reveal new players in CAP-H2-mediated pairing regulation and the intimate interplay between inter-chromosomal and intra-chromosomal 3D interactions.

HiCrayon reveals distinct layers of multi-state 3D chromatin organization.

The co-visualization of chromatin conformation with 1D 'omics data is key to the multi-omics driven data analysis of 3D genome organization. Chromatin contact maps are often shown as 2D heatmaps and visually compared to 1D genomic data by simple juxtaposition. While common, this strategy is imprecise, placing the onus on the reader to align features with each other. To remedy this, we developed HiCrayon, an interactive tool that facilitates the integration of 3D chromatin organization maps and 1D datasets. This visualization method integrates data from genomic assays directly into the chromatin contact map by coloring interactions according to 1D signal. HiCrayon is implemented using R shiny and python to create a graphical user interface (GUI) application, available in both web or containerized format to promote accessibility. HiCrayon is implemented in R, and includes a graphical user interface (GUI), as well as a slimmed-down web-based version that lets users quickly produce publication-ready images. We demonstrate the utility of HiCrayon in visualizing the effectiveness of compartment calling and the relationship between ChIP-seq and various features of chromatin organization. We also demonstrate the improved visualization of other 3D genomic phenomena, such as differences between loops associated with CTCF/cohesin vs. those associated with H3K27ac. We then demonstrate HiCrayon's visualization of organizational changes that occur during differentiation and use HiCrayon to detect compartment patterns that cannot be assigned to either A or B compartments, revealing a distinct 3rd chromatin compartment. Overall, we demonstrate the utility of co-visualizing 2D chromatin conformation with 1D genomic signals within the same matrix to reveal fundamental aspects of genome organization. Local version: https://github.com/JRowleyLab/HiCrayon Web version: https://jrowleylab.com/HiCrayon.

Regulation of CTCF loop formation during pancreatic cell differentiation.

Transcription reprogramming during cell differentiation involves targeting enhancers to genes responsible for establishment of cell fates. To understand the contribution of CTCF-mediated chromatin organization to cell lineage commitment, we analyzed 3D chromatin architecture during the differentiation of human embryonic stem cells into pancreatic islet organoids. We find that CTCF loops are formed and disassembled at different stages of the differentiation process by either recruitment of CTCF to new anchor sites or use of pre-existing sites not previously involved in loop formation. Recruitment of CTCF to new sites in the genome involves demethylation of H3K9me3 to H3K9me2, demethylation of DNA, recruitment of pioneer factors, and positioning of nucleosomes flanking the new CTCF sites. Existing CTCF sites not involved in loop formation become functional loop anchors via the establishment of new cohesin loading sites containing NIPBL and YY1 at sites between the new anchors. In both cases, formation of new CTCF loops leads to strengthening of enhancer promoter interactions and increased transcription of genes adjacent to loop anchors. These results suggest an important role for CTCF and cohesin in controlling gene expression during cell differentiation.

Characterization of m6A modifiers and RNA modifications in uterine fibroids.

Uterine leiomyoma or fibroids are the most common prevalent noncancerous tumors of the uterine muscle layer. Common symptoms associated with fibroids include pelvic pain, heavy menstrual bleeding, anemia, and pelvic pressure. These tumors are a leading cause of gynecological care but lack long-term therapy as the origin and development of fibroids are not well understood. Several next-generation sequencing technologies have been performed to identify the underlying genetic and epigenetic basis of fibroids. However, there remains a systemic gap in our understanding of molecular and biological process that define uterine fibroids. Recent epitranscriptomics studies have unraveled RNA modifications that are associated with all forms of RNA and are thought to influence both normal physiological functions and the progression of diseases. We quantified RNA expression profiles by analyzing publicly available RNA-seq data for 15 known epigenetic mediators to identify their expression profile in uterine fibroids compared to myometrium. To validate our findings, we performed RT-qPCR on a separate cohort of uterine fibroids targeting these modifiers confirming our RNA-seq data. We then examined protein profiles of key m6A modifiers in fibroids and their matched myometrium. In concordance with our RNA expression profiles, no significant differences were observed in these proteins in uterine fibroids compared to myometrium. To determine abundance of RNA modifications, mRNA and small RNA from fibroids and matched myometrium were analyzed by UHPLC MS/MS. In addition to the prevalent N6-methyladenosine (m6A), we identified 11 other known modifiers but did not identify any aberrant expression in fibroids. We then mined a previously published dataset and identified differential expression of m6A modifiers that were specific to fibroid genetic sub-type. Our analysis also identified m6A consensus motifs on genes previously identified to be dysregulated in uterine fibroids. Overall, using state-of-the-art mass spectrometry, RNA expression and protein profiles, we characterized and identified differentially expressed m6A modifiers in relation to driver mutations. Despite the use of several different approaches, we identified limited differential expression of RNA modifiers and associated modifications in uterine fibroids. However, considering the highly heterogenous genomic and cellular nature of fibroids, and the possible contribution of single molecule m6A modifications to fibroid pathology, there is a need for greater in-depth characterization of m6A marks and modifiers in a larger and varied patient cohort.

Chromatin alternates between A and B compartments at kilobase scale for subgenic organization.

Nuclear compartments are prominent features of 3D chromatin organization, but sequencing depth limitations have impeded investigation at ultra fine-scale. CTCF loops are generally studied at a finer scale, but the impact of looping on proximal interactions remains enigmatic. Here, we critically examine nuclear compartments and CTCF loop-proximal interactions using a combination of in situ Hi-C at unparalleled depth, algorithm development, and biophysical modeling. Producing a large Hi-C map with 33 billion contacts in conjunction with an algorithm for performing principal component analysis on sparse, super massive matrices (POSSUMM), we resolve compartments to 500 bp. Our results demonstrate that essentially all active promoters and distal enhancers localize in the A compartment, even when flanking sequences do not. Furthermore, we find that the TSS and TTS of paused genes are often segregated into separate compartments. We then identify diffuse interactions that radiate from CTCF loop anchors, which correlate with strong enhancer-promoter interactions and proximal transcription. We also find that these diffuse interactions depend on CTCF's RNA binding domains. In this work, we demonstrate features of fine-scale chromatin organization consistent with a revised model in which compartments are more precise than commonly thought while CTCF loops are more protracted.

Considerations and caveats for analyzing chromatin compartments.

Genomes are organized into nuclear compartments, separating active from inactive chromatin. Chromatin compartments are readily visible in a large number of species by experiments that map chromatin conformation genome-wide. When analyzing these maps, a common step is the identification of genomic intervals that interact within A (active) and B (inactive) compartments. It has also become increasingly common to identify and analyze subcompartments. We review different strategies to identify A/B and subcompartment intervals, including a discussion of various machine-learning approaches to predict these features. We then discuss the strengths and limitations of current strategies and examine how these aspects of analysis may have impacted our understanding of chromatin compartments.

Somatic chromosome pairing has a determinant impact on 3D chromatin organization.

In the nucleus, chromatin is intricately structured into multiple layers of 3D organization important for genome activity. How distinct layers influence each other is not well understood. In particular, the contribution of chromosome pairing to 3D chromatin organization has been largely neglected. Here, we address this question in Drosophila, an organism that shows robust chromosome pairing in interphasic somatic cells. The extent of chromosome pairing depends on the balance between pairing and anti-pairing factors, with the anti-pairing activity of the CAP-H2 condensin II subunit being the best documented. Here, we identify the zinc-finger protein Z4 as a strong anti-pairer that interacts with and mediates the chromatin binding of CAP-H2. We also report that hyperosmotic cellular stress induces fast and reversible chromosome unpairing that depends on Z4/CAP-H2. And, most important, by combining Z4 depletion and osmostress, we show that chromosome pairing reinforces intrachromosomal 3D interactions. On the one hand, pairing facilitates RNAPII occupancy that correlates with enhanced intragenic gene-loop interactions. In addition, acting at a distance, pairing reinforces chromatin-loop interactions mediated by Polycomb (Pc). In contrast, chromosome pairing does not affect which genomic intervals segregate to active (A) and inactive (B) compartments, with only minimal effects on the strength of A-A compartmental interactions. Altogether, our results unveil the intimate interplay between inter-chromosomal and intra-chromosomal 3D interactions, unraveling the interwoven relationship between different layers of chromatin organization and the essential contribution of chromosome pairing.

Defining a cohort of anemia-activated cis elements reveals a mechanism promoting erythroid precursor function.

Acute anemia elicits broad transcriptional changes in erythroid progenitors and precursors. We previously discovered a cis-regulatory transcriptional enhancer at the sterile alpha motif domain-14 enhancer locus (S14E), defined by a CANNTG-spacer-AGATAA composite motif and occupied by GATA1 and TAL1 transcription factors, is required for survival in severe anemia. However, S14E is only 1 of dozens of anemia-activated genes containing similar motifs. In a mouse model of acute anemia, we identified populations of expanding erythroid precursors, which increased expression of genes that contain S14E-like cis elements. We reveal that several S14E-like cis elements provide important transcriptional control of newly identified anemia-inducing genes, including the Ssx-2 interacting protein (Ssx2ip). Ssx2ip expression was determined to play an important role in erythroid progenitor/precursor cell activities, cell cycle regulation, and cell proliferation. Over a weeklong course of acute anemia recovery, we observed that erythroid gene activation mediated by S14E-like cis elements occurs during a phase coincident with low hematocrit and high progenitor activities, with distinct transcriptional programs activated at earlier and later time points. Our results define a genome-wide mechanism in which S14E-like enhancers control transcriptional responses during erythroid regeneration. These findings provide a framework to understand anemia-specific transcriptional mechanisms, ineffective erythropoiesis, anemia recovery, and phenotypic variability within human populations.

Sex-specific multi-level 3D genome dynamics in the mouse brain.

The female mammalian brain exhibits sex hormone-driven plasticity during the reproductive period. Recent evidence implicates chromatin dynamics in gene regulation underlying this plasticity. However, whether ovarian hormones impact higher-order chromatin organization in post-mitotic neurons in vivo is unknown. Here, we mapped the 3D genome of ventral hippocampal neurons across the oestrous cycle and by sex in mice. In females, we find cycle-driven dynamism in 3D chromatin organization, including in oestrogen response elements-enriched X chromosome compartments, autosomal CTCF loops, and enhancer-promoter interactions. With rising oestrogen levels, the female 3D genome becomes more similar to the male 3D genome. Cyclical enhancer-promoter interactions are partially associated with gene expression and enriched for brain disorder-relevant genes and pathways. Our study reveals unique 3D genome dynamics in the female brain relevant to female-specific gene regulation, neuroplasticity, and disease risk.

Elevating SOX2 Downregulates MYC through a SOX2:MYC Signaling Axis and Induces a Slowly Cycling Proliferative State in Human Tumor Cells.

Slowly cycling/infrequently proliferating tumor cells present a clinical challenge due to their ability to evade treatment. Previous studies established that high levels of SOX2 in both fetal and tumor cells restrict cell proliferation and induce a slowly cycling state. However, the mechanisms through which elevated SOX2 levels inhibit tumor cell proliferation have not been identified. To identify common mechanisms through which SOX2 elevation restricts tumor cell proliferation, we initially performed RNA-seq using two diverse tumor cell types. SOX2 elevation in both cell types downregulated MYC target genes. Consistent with these findings, elevating SOX2 in five cell lines representing three different human cancer types decreased MYC expression. Importantly, the expression of a dominant-negative MYC variant, omomyc, recapitulated many of the effects of SOX2 on proliferation, cell cycle, gene expression, and biosynthetic activity. We also demonstrated that rescuing MYC activity in the context of elevated SOX2 induces cell death, indicating that the downregulation of MYC is a critical mechanistic step necessary to maintain survival in the slowly cycling state induced by elevated SOX2. Altogether, our findings uncover a novel SOX2:MYC signaling axis and provide important insights into the molecular mechanisms through which SOX2 elevation induces a slowly cycling proliferative state.

Implications of Dosage Deficiencies in CTCF and Cohesin on Genome Organization, Gene Expression, and Human Neurodevelopment.

Properly organizing DNA within the nucleus is critical to ensure normal downstream nuclear functions. CTCF and cohesin act as major architectural proteins, working in concert to generate thousands of high-intensity chromatin loops. Due to their central role in loop formation, a massive research effort has been dedicated to investigating the mechanism by which CTCF and cohesin create these loops. Recent results lead to questioning the direct impact of CTCF loops on gene expression. Additionally, results of controlled depletion experiments in cell lines has indicated that genome architecture may be somewhat resistant to incomplete deficiencies in CTCF or cohesin. However, heterozygous human genetic deficiencies in CTCF and cohesin have illustrated the importance of their dosage in genome architecture, cellular processes, animal behavior, and disease phenotypes. Thus, the importance of considering CTCF or cohesin levels is especially made clear by these heterozygous germline variants that characterize genetic syndromes, which are increasingly recognized in clinical practice. Defined primarily by developmental delay and intellectual disability, the phenotypes of CTCF and cohesin deficiency illustrate the importance of architectural proteins particularly in neurodevelopment. We discuss the distinct roles of CTCF and cohesin in forming chromatin loops, highlight the major role that dosage of each protein plays in the amplitude of observed effects on gene expression, and contrast these results to heterozygous mutation phenotypes in murine models and clinical patients. Insights highlighted by this comparison have implications for future research into these newly emerging genetic syndromes.

Ecdysoneless Protein Regulates Viral and Cellular mRNA Splicing to Promote Cervical Oncogenesis.

High-risk human papillomaviruses (HPV), exemplified by HPV16/18, are causally linked to human cancers of the anogenital tract, skin, and upper aerodigestive tract. Previously, we identified Ecdysoneless (ECD) protein, the human homolog of the Drosophila ecdysoneless gene, as a novel HPV16 E6-interacting protein. Here, we show that ECD, through its C-terminal region, selectively binds to high-risk but not to low-risk HPV E6 proteins. We demonstrate that ECD is overexpressed in cervical and head and neck squamous cell carcinoma (HNSCC) cell lines as well as in tumor tissues. Using The Cancer Genome Atlas dataset, we show that ECD mRNA overexpression predicts shorter survival in patients with cervical and HNSCC. We demonstrate that ECD knockdown in cervical cancer cell lines led to impaired oncogenic behavior, and ECD co-overexpression with E7 immortalized primary human keratinocytes. RNA-sequencing analyses of SiHa cells upon ECD knockdown showed to aberrations in E6/E7 RNA splicing, as well as RNA splicing of several HPV oncogenesis-linked cellular genes, including splicing of components of mRNA splicing machinery itself. Taken together, our results support a novel role of ECD in viral and cellular mRNA splicing to support HPV-driven oncogenesis. IMPLICATIONS: This study links ECD overexpression to poor prognosis and shorter survival in HNSCC and cervical cancers and identifies a critical role of ECD in cervical oncogenesis through regulation of viral and cellular mRNA splicing.

The inhibition of LSD1 via sequestration contributes to tau-mediated neurodegeneration.

Tauopathies are a class of neurodegenerative diseases associated with pathological tau. Despite many advances in our understanding of these diseases, the direct mechanism through which tau contributes to neurodegeneration remains poorly understood. Previously, our laboratory implicated the histone demethylase LSD1 in tau-induced neurodegeneration by showing that LSD1 localizes to pathological tau aggregates in Alzheimer's disease cases, and that it is continuously required for the survival of hippocampal and cortical neurons in mice. Here, we utilize the P301S tauopathy mouse model to demonstrate that pathological tau can exclude LSD1 from the nucleus in neurons. In addition, we show that reducing LSD1 in these mice is sufficient to highly exacerbate tau-mediated neurodegeneration and tau-induced gene expression changes. Finally, we find that overexpressing LSD1 in the hippocampus of tauopathy mice, even after pathology has formed, is sufficient to significantly delay neurodegeneration and counteract tau-induced expression changes. These results suggest that inhibiting LSD1 via sequestration contributes to tau-mediated neurodegeneration. Thus, LSD1 is a promising therapeutic target for tauopathies such as Alzheimer's disease.

Analysis of Hi-C data using SIP effectively identifies loops in organisms from C. elegans to mammals.

Chromatin loops are a major component of 3D nuclear organization, visually apparent as intense point-to-point interactions in Hi-C maps. Identification of these loops is a critical part of most Hi-C analyses. However, current methods often miss visually evident CTCF loops in Hi-C data sets from mammals, and they completely fail to identify high intensity loops in other organisms. We present SIP, Significant Interaction Peak caller, and SIPMeta, which are platform independent programs to identify and characterize these loops in a time- and memory-efficient manner. We show that SIP is resistant to noise and sequencing depth, and can be used to detect loops that were previously missed in human cells as well as loops in other organisms. SIPMeta corrects for a common visualization artifact by accounting for Manhattan distance to create average plots of Hi-C and HiChIP data. We then demonstrate that the use of SIP and SIPMeta can lead to biological insights by characterizing the contribution of several transcription factors to CTCF loop stability in human cells. We also annotate loops associated with the SMC component of the dosage compensation complex (DCC) in Caenorhabditis elegans and demonstrate that loop anchors represent bidirectional blocks for symmetrical loop extrusion. This is in contrast to the asymmetrical extrusion until unidirectional blockage by CTCF that is presumed to occur in mammals. Using HiChIP and multiway ligation events, we then show that DCC loops form a network of strong interactions that may contribute to X Chromosome-wide condensation in C. elegans hermaphrodites.

Evolutionary History and Activity of RNase H1-Like Proteins in Arabidopsis thaliana.

RNase H1 is an endonuclease specific toward the RNA strand of RNA:DNA hybrids. Members of this protein family are present in most living organisms and are essential for removing RNA that base pairs with DNA. It prevents detrimental effects of RNA:DNA hybrids and is involved in several biological processes. Arabidopsis thaliana has been previously shown to contain three genes encoding RNase H1 proteins that localize to three distinct cellular compartments. We show that these genes originate from two gene duplication events. One occurred in the common ancestor of dicots and produced nuclear and organellar RNase H1 paralogs. Second duplication occurred in the common ancestor of Brassicaceae and produced mitochondrial- and plastid-localized proteins. These proteins have the canonical RNase H1 activity, which requires at least four ribonucleotides for endonucleolytic digestion. Analysis of mutants in the RNase H1 genes revealed that the nuclear RNH1A and mitochondrial RNH1B are dispensable for development under normal growth conditions. However, the presence of at least one organellar RNase H1 (RNH1B or RNH1C) is required for embryonic development. The plastid-localized RNH1C affects plastid DNA copy number and sensitivity to replicative stress. Our results present the evolutionary history of RNH1 proteins in A. thaliana, demonstrate their canonical RNase H1 activity and indicate their role in early embryonic development.