Pairtools: from sequencing data to chromosome contacts

The brittle star genome illuminates the genetic basis of animal appendage regeneration

Species within nearly all extant animal lineages are capable of regenerating body parts. However, it remains unclear whether the gene expression programme controlling regeneration is evolutionarily conserved. Brittle stars are a species-rich class of echinoderms with outstanding regenerative abilities, but investigations into the genetic bases of regeneration in this group have been hindered by the limited genomic resources. Here we report a chromosome-scale genome assembly for the brittle star Amphiura filiformis. We show that the brittle star genome is the most rearranged among echinoderms sequenced so far, featuring a reorganized Hox cluster reminiscent of the rearrangements observed in sea urchins. In addition, we performed an extensive profiling of gene expression during brittle star adult arm regeneration and identified sequential waves of gene expression governing wound healing, proliferation and differentiation. We conducted comparative transcriptomic analyses with other invertebrate and vertebrate models for appendage regeneration and uncovered hundreds of genes with conserved expression dynamics, particularly during the proliferative phase of regeneration. Our findings emphasize the crucial importance of echinoderms to detect long-range expression conservation between vertebrates and classical invertebrate regeneration model systems.

Extensive long-range polycomb interactions and weak compartmentalization are hallmarks of human neuronal 3D genome

Chromatin architecture regulates gene expression and shapes cellular identity, particularly in neuronal cells. Specifically, polycomb group (PcG) proteins enable establishment and maintenance of neuronal cell type by reorganizing chromatin into repressive domains that limit the expression of fate-determining genes and sustain distinct gene expression patterns in neurons. Here, we map the 3D genome architecture in neuronal and non-neuronal cells isolated from the Wernicke's area of four human brains and comprehensively analyze neuron-specific aspects of chromatin organization. We find that genome segregation into active and inactive compartments is greatly reduced in neurons compared to other brain cells. Furthermore, neuronal Hi-C maps reveal strong long-range interactions, forming a specific network of PcG-mediated contacts in neurons that is nearly absent in other brain cells. These interacting loci contain developmental transcription factors with repressed expression in neurons and other mature brain cells. But only in neurons, they are rich in bivalent promoters occupied by H3K4me3 histone modification together with H3K27me3, which points to a possible functional role of PcG contacts in neurons. Importantly, other layers of chromatin organization also exhibit a distinct structure in neurons, characterized by an increase in short-range interactions and a decrease in long-range ones.

Chromosomal-scale genome assembly of the Mediterranean mussel Mytilus galloprovincialis

The Mediterranean mussel, Mytilus galloprovincialis, is a significant marine bivalve species that has ecological and economic importance. This species is robustly resilient and highly invasive. Despite the scientific and commercial interest in studying its biology and aquaculture, there remains a need for a high-quality, chromosome-scale reference genome. In this study, we have assembled a high-quality chromosome-scale reference genome for M. galloprovincialis. The total length of our reference genome is 1.41 Gb, with a scaffold N50 sequence length of 96.9 Mb. BUSCO analysis revealed a 97.5% completeness based on complete BUSCOs. Compared to the four other available M. galloprovincialis assemblies, the assembly described here is dramatically improved in both contiguity and completeness. This new reference genome will greatly contribute to a deeper understanding of the resilience and invasiveness of M. galloprovincialis.

Systematic assessment of ISWI subunits shows that NURF creates local accessibility for CTCF

Catalytic activity of the imitation switch (ISWI) family of remodelers is critical for nucleosomal organization and DNA binding of certain transcription factors, including the insulator protein CTCF. Here we define the contribution of individual subcomplexes by deriving a panel of isogenic mouse stem cell lines, each lacking one of six ISWI accessory subunits. Individual deletions of subunits of either CERF, RSF, ACF, WICH or NoRC subcomplexes only moderately affect the chromatin landscape, while removal of the NURF-specific subunit BPTF leads to a strong reduction in chromatin accessibility and SNF2H ATPase localization around CTCF sites. This affects adjacent nucleosome occupancy and CTCF binding. At a group of sites with reduced chromatin accessibility, CTCF binding persists but cohesin occupancy is reduced, resulting in decreased insulation. These results suggest that CTCF binding can be separated from its function as an insulator in nuclear organization and identify a specific role for NURF in mediating SNF2H localization and chromatin opening at bound CTCF sites.

Accelerating 3D genomics data analysis with Microcket

The three-dimensional (3D) organization of genome is fundamental to cell biology. To explore 3D genome, emerging high-throughput approaches have produced billions of sequencing reads, which is challenging and time-consuming to analyze. Here we present Microcket, a package for mapping and extracting interacting pairs from 3D genomics data, including Hi-C, Micro-C, and derivant protocols. Microcket utilizes a unique read-stitch strategy that takes advantage of the long read cycles in modern DNA sequencers; benchmark evaluations reveal that Microcket runs much faster than the current tools along with improved mapping efficiency, and thus shows high potential in accelerating and enhancing the biological investigations into 3D genome. Microcket is freely available at https://github.com/hellosunking/Microcket .

CRISPRi screens identify the lncRNA, LOUP, as a multifunctional locus regulating macrophage differentiation and inflammatory signaling

Long noncoding RNAs (lncRNAs) account for the largest portion of RNA from the transcriptome, yet most of their functions remain unknown. Here, we performed two independent high-throughput CRISPRi screens to understand the role of lncRNAs in monocyte function and differentiation. The first was a reporter-based screen to identify lncRNAs that regulate TLR4-NFkB signaling in human monocytes and the second screen identified lncRNAs involved in monocyte to macrophage differentiation. We successfully identified numerous noncoding and protein-coding genes that can positively or negatively regulate inflammation and differentiation. To understand the functional roles of lncRNAs in both processes, we chose to further study the lncRNA LOUP [lncRNA originating from upstream regulatory element of SPI1 (also known as PU.1)], as it emerged as a top hit in both screens. Not only does LOUP regulate its neighboring gene, the myeloid fate-determining factor SPI1, thereby affecting monocyte to macrophage differentiation, but knockdown of LOUP leads to a broad upregulation of NFkB-targeted genes at baseline and upon TLR4-NFkB activation. LOUP also harbors three small open reading frames capable of being translated and are responsible for LOUP's ability to negatively regulate TLR4/NFkB signaling. This work emphasizes the value of high-throughput screening to rapidly identify functional lncRNAs in the innate immune system.

Genome of the endangered eastern quoll (Dasyurus viverrinus) reveals signatures of historical decline and pelage color evolution

The eastern quoll (Dasyurus viverrinus) is an endangered marsupial native to Australia. Since the extirpation of its mainland populations in the 20th century, wild eastern quolls have been restricted to two islands at the southern end of their historical range. Eastern quolls are the subject of captive breeding programs and attempts have been made to re-establish a population in mainland Australia. However, few resources currently exist to guide the genetic management of this species. Here, we generated a reference genome for the eastern quoll with gene annotations supported by multi-tissue transcriptomes. Our assembly is among the most complete marsupial genomes currently available. Using this assembly, we infer the species' demographic history, identifying potential evidence of a long-term decline beginning in the late Pleistocene. Finally, we identify a deletion at the ASIP locus that likely underpins pelage color differences between the eastern quoll and the closely related Tasmanian devil (Sarcophilus harrisii).

Rules of engagement for condensins and cohesins guide mitotic chromosome formation

During mitosis, interphase chromatin is rapidly converted into rod-shaped mitotic chromosomes. Using Hi-C, imaging, proteomics and polymer modeling, we determine how the activity and interplay between loop-extruding SMC motors accomplishes this dramatic transition. Our work reveals rules of engagement for SMC complexes that are critical for allowing cells to refold interphase chromatin into mitotic chromosomes. We find that condensin disassembles interphase chromatin loop organization by evicting or displacing extrusive cohesin. In contrast, condensin bypasses cohesive cohesins, thereby maintaining sister chromatid cohesion while separating the sisters. Studies of mitotic chromosomes formed by cohesin, condensin II and condensin I alone or in combination allow us to develop new models of mitotic chromosome conformation. In these models, loops are consecutive and not overlapping, implying that condensins do not freely pass one another but stall upon encountering each other. The dynamics of Hi-C interactions and chromosome morphology reveal that during prophase loops are extruded in vivo at ~1-3 kb/sec by condensins as they form a disordered discontinuous helical scaffold within individual chromatids.

Limited genomic signatures of population collapse in the critically endangered black abalone (Haliotis cracherodii)

The black abalone, Haliotis cracherodii, is a large, long-lived marine mollusc that inhabits rocky intertidal habitats along the coast of California and Mexico. In 1985, populations were impacted by a bacterial disease known as withering syndrome (WS) that wiped out >90% of individuals, leading to the closure of all U.S. black abalone fisheries since 1993. Current conservation strategies include restoring diminished populations by translocating healthy individuals. However, population collapse on this scale may have dramatically lowered genetic diversity and strengthened geographic differentiation, making translocation-based recovery contentious. Additionally, the current prevalence of WS remains unknown. To address these uncertainties, we sequenced and analysed the genomes of 133 black abalone individuals from across their present range. We observed no spatial genetic structure among black abalone, with the exception of a single chromosomal inversion that increases in frequency with latitude. Outside the inversion, genetic differentiation between sites is minimal and does not scale with either geographic distance or environmental dissimilarity. Genetic diversity appears uniformly high across the range. Demographic inference does indicate a severe population bottleneck beginning just 15 generations in the past, but this decline is short lived, with present-day size far exceeding the pre-bottleneck status quo. Finally, we find the bacterial agent of WS is equally present across the sampled range, but only in 10% of individuals. The lack of population genetic structure, uniform diversity and prevalence of WS bacteria indicates that translocation could be a valid and low-risk means of population restoration for black abalone species' recovery.

High resolution maps of chromatin reorganization through mouse meiosis reveal novel features of the 3D meiotic structure

When germ cells transition from the mitotic cycle into meiotic prophase I (MPI), chromosomes condense into an array of chromatin loops that are required to promote homolog pairing and genetic recombination. To identify the changes in chromosomal conformation, we isolated nuclei on a trajectory from spermatogonia to the end of MPI. At each stage along this trajectory, we built genomic interaction maps with the highest temporal and spatial resolution to date. The changes in chromatin folding coincided with a concurrent decline in mitotic cohesion and a rise in meiotic cohesin complexes. We found that the stereotypical large-scale A and B compartmentalization was lost during meiotic prophase I alongside the loss of topological associating domains (TADs). Still, local subcompartments were detected and maintained throughout meiosis. The enhanced Micro-C resolution revealed that, despite the loss of TADs, higher frequency contact sites between two loci were detectable during meiotic prophase I coinciding with CTCF bound sites. The pattern of interactions around these CTCF sites with their neighboring loci showed that CTCF sites were often anchoring the meiotic loops. Additionally, the localization of CTCF to the meiotic axes indicated that these anchors were at the base of loops. Strikingly, even in the face of the dramatic reconfiguration of interphase chromatin into a condensed loop-array, the interactions between regulatory elements remained well preserved. This establishes a potential mechanism for how the meiotic chromatin maintains active transcription within a highly structured genome. In summary, the high temporal and spatial resolution of these data revealed previously unappreciated aspects of mammalian meiotic chromatin organization.

KDM3B inhibitors disrupt the oncogenic activity of PAX3-FOXO1 in fusion-positive rhabdomyosarcoma

Fusion-positive rhabdomyosarcoma (FP-RMS) is an aggressive pediatric sarcoma driven primarily by the PAX3-FOXO1 fusion oncogene, for which therapies targeting PAX3-FOXO1 are lacking. Here, we screen 62,643 compounds using an engineered cell line that monitors PAX3-FOXO1 transcriptional activity identifying a hitherto uncharacterized compound, P3FI-63. RNA-seq, ATAC-seq, and docking analyses implicate histone lysine demethylases (KDMs) as its targets. Enzymatic assays confirm the inhibition of multiple KDMs with the highest selectivity for KDM3B. Structural similarity search of P3FI-63 identifies P3FI-90 with improved solubility and potency. Biophysical binding of P3FI-90 to KDM3B is demonstrated using NMR and SPR. P3FI-90 suppresses the growth of FP-RMS in vitro and in vivo through downregulating PAX3-FOXO1 activity, and combined knockdown of KDM3B and KDM1A phenocopies P3FI-90 effects. Thus, we report KDM inhibitors P3FI-63 and P3FI-90 with the highest specificity for KDM3B. Their potent suppression of PAX3-FOXO1 activity indicates a possible therapeutic approach for FP-RMS and other transcriptionally addicted cancers.

YY1-dependent transcriptional regulation manifests at the morula stage

YY1 plays multifaceted roles in various cell types. We recently reported that YY1 regulates nucleosome organization in early mouse embryos. However, despite the impaired nucleosome organization in the absence of YY1, the transcriptome was minimally affected in eight-cell embryos. We then hypothesized that YY1 might prepare a chromatin environment to regulate gene expression at later stages. To test this possibility, we performed a transcriptome analysis at the morula stage. We found that a substantial number of genes are aberrantly expressed in the absence of YY1. Furthermore, our analysis revealed that YY1 is required for the transcription of LINE-1 retrotransposons.

Histone binding of ASF1 is required for fruiting body development but not for genome stability in the filamentous fungus Sordaria macrospora

IMPORTANCE

Histone chaperones are proteins that are involved in nucleosome assembly and disassembly and can therefore influence all DNA-dependent processes including transcription, DNA replication, and repair. ASF1 is a histone chaperone that is conserved throughout eukaryotes. In contrast to most other multicellular organisms, a deletion mutant of asf1 in the fungus Sordaria macrospora is viable; however, the mutant is sterile. In this study, we could show that the histone-binding ability of ASF1 is required for fertility in S. macrospora, whereas the function of ASF1 in maintenance of genome stability does not require histone binding. We also showed that the histone modifications H3K27me3 and H3K56ac are misregulated in the Δasf1 mutant. Furthermore, we identified a large duplication on chromosome 2 of the mutant strain that is genetically linked to the Δasf1 allele present on chromosome 6, suggesting that viability of the mutant might depend on the presence of the duplicated region.

A chromosome-level reference genome for the common octopus, Octopus vulgaris (Cuvier, 1797)

Cephalopods are emerging animal models and include iconic species for studying the link between genomic innovations and physiological and behavioral complexities. Coleoid cephalopods possess the largest nervous system among invertebrates, both for cell counts and brain-to-body ratio. Octopus vulgaris has been at the center of a long-standing tradition of research into diverse aspects of cephalopod biology, including behavioral and neural plasticity, learning and memory recall, regeneration, and sophisticated cognition. However, no chromosome-scale genome assembly was available for O. vulgaris to aid in functional studies. To fill this gap, we sequenced and assembled a chromosome-scale genome of the common octopus, O. vulgaris. The final assembly spans 2.8 billion basepairs, 99.34% of which are in 30 chromosome-scale scaffolds. Hi-C heatmaps support a karyotype of 1n = 30 chromosomes. Comparisons with other octopus species' genomes show a conserved octopus karyotype and a pattern of local genome rearrangements between species. This new chromosome-scale genome of O. vulgaris will further facilitate research in all aspects of cephalopod biology, including various forms of plasticity and the neural machinery underlying sophisticated cognition, as well as an understanding of cephalopod evolution.

High-resolution CTCF footprinting reveals impact of chromatin state on cohesin extrusion dynamics

DNA looping is vital for establishing many enhancer-promoter interactions. While CTCF is known to anchor many cohesin-mediated loops, the looped chromatin fiber appears to predominantly exist in a poorly characterized actively extruding state. To better characterize extruding chromatin loop structures, we used CTCF MNase HiChIP data to determine both CTCF binding at high resolution and 3D contact information. Here we present FactorFinder, a tool that identifies CTCF binding sites at near base-pair resolution. We leverage this substantial advance in resolution to determine that the fully extruded (CTCF-CTCF) state is rare genome-wide with locus-specific variation from ~1-10%. We further investigate the impact of chromatin state on loop extrusion dynamics, and find that active enhancers and RNA Pol II impede cohesin extrusion, facilitating an enrichment of enhancer-promoter contacts in the partially extruded loop state. We propose a model of topological regulation whereby the transient, partially extruded states play active roles in transcription.

Droplet-based single-cell joint profiling of histone modifications and transcriptomes

We previously reported Paired-Tag, a combinatorial indexing-based method that can simultaneously map histone modifications and gene expression at single-cell resolution at scale. However, the lengthy procedure of Paired-Tag has hindered its general adoption in the community. To address this bottleneck, we developed a droplet-based Paired-Tag protocol that is faster and more accessible than the previous method. Using cultured mammalian cells and primary brain tissues, we demonstrate its superior performance at identifying candidate cis-regulatory elements and associating their dynamic chromatin state to target gene expression in each constituent cell type in a complex tissue.

Integrative PTEN Enhancer Discovery Reveals a New Model of Enhancer Organization

Enhancers possess both structural elements mediating promoter looping and functional elements mediating gene expression. Traditional models of enhancer-mediated gene regulation imply genomic overlap or immediate adjacency of these elements. We test this model by combining densely-tiled CRISPRa screening with nucleosome-resolution Region Capture Micro-C topology analysis. Using this integrated approach, we comprehensively define the cis-regulatory landscape for the tumor suppressor PTEN, identifying and validating 10 distinct enhancers and defining their 3D spatial organization. Unexpectedly, we identify several long-range functional enhancers whose promoter proximity is facilitated by chromatin loop anchors several kilobases away, and demonstrate that accounting for this spatial separation improves the computational prediction of validated enhancers. Thus, we propose a new model of enhancer organization incorporating spatial separation of essential functional and structural components.

The 3D chromatin landscape of rhabdomyosarcoma

Rhabdomyosarcoma (RMS) is a pediatric soft tissue cancer with a lack of precision therapy options for patients. We hypothesized that with a general paucity of known mutations in RMS, chromatin structural driving mechanisms are essential for tumor proliferation. Thus, we carried out high-depth in situ Hi-C in representative cell lines and patient-derived xenografts (PDXs) to define chromatin architecture in each major RMS subtype. We report a comprehensive 3D chromatin structural analysis and characterization of fusion-positive (FP-RMS) and fusion-negative RMS (FN-RMS). We have generated spike-in in situ Hi-C chromatin interaction maps for the most common FP-RMS and FN-RMS cell lines and compared our data with PDX models. In our studies, we uncover common and distinct structural elements in large Mb-scale chromatin compartments, tumor-essential genes within variable topologically associating domains and unique patterns of structural variation. Our high-depth chromatin interactivity maps and comprehensive analyses provide context for gene regulatory events and reveal functional chromatin domains in RMS.

MacroH2A restricts inflammatory gene expression in melanoma cancer-associated fibroblasts by coordinating chromatin looping

MacroH2A has established tumour suppressive functions in melanoma and other cancers, but an unappreciated role in the tumour microenvironment. Using an autochthonous, immunocompetent mouse model of melanoma, we demonstrate that mice devoid of macroH2A variants exhibit increased tumour burden compared with wild-type counterparts. MacroH2A-deficient tumours accumulate immunosuppressive monocytes and are depleted of functional cytotoxic T cells, characteristics consistent with a compromised anti-tumour response. Single cell and spatial transcriptomics identify increased dedifferentiation along the neural crest lineage of the tumour compartment and increased frequency and activation of cancer-associated fibroblasts following macroH2A loss. Mechanistically, macroH2A-deficient cancer-associated fibroblasts display increased myeloid chemoattractant activity as a consequence of hyperinducible expression of inflammatory genes, which is enforced by increased chromatin looping of their promoters to enhancers that gain H3K27ac. In summary, we reveal a tumour suppressive role for macroH2A variants through the regulation of chromatin architecture in the tumour stroma with potential implications for human melanoma.

Variant-to-function analysis of the childhood obesity chr12q13 locus implicates rs7132908 as a causal variant within the 3' UTR of FAIM2

The ch12q13 obesity locus is among the most significant childhood obesity loci identified in genome-wide association studies. This locus resides in a non-coding region within FAIM2; thus, the underlying causal variant(s) presumably influence disease susceptibility via an influence on cis-regulation within the genomic region. We implicated rs7132908 as a putative causal variant at this locus leveraging a combination of our inhouse 3D genomic data, public domain datasets, and several computational approaches. Using a luciferase reporter assay in human primary astrocytes, we observed allele-specific cis-regulatory activity of the immediate region harboring rs7132908. Motivated by this finding, we went on to generate isogenic human embryonic stem cell lines homozygous for either rs7132908 allele with CRISPR-Cas9 homology-directed repair to assess changes in gene expression due to genotype and chromatin accessibility throughout a differentiation to hypothalamic neurons, a key cell type known to regulate feeding behavior. We observed that the rs7132908 obesity risk allele influenced the expression of FAIM2 along with other genes, decreased the proportion of neurons produced during differentiation, up-regulated cell death gene sets, and conversely down-regulated neuron differentiation gene sets. We have therefore functionally validated rs7132908 as a causal obesity variant which temporally regulates nearby effector genes at the ch12q13 locus and influences neurodevelopment and survival.

Extrusion fountains are hallmarks of chromosome organization emerging upon zygotic genome activation

The first activation of gene expression during development (zygotic genome activation, ZGA) is accompanied by massive changes in chromosome organization. The connection between these two processes remains unknown. Using Hi-C for zebrafish embryos, we found that chromosome folding starts by establishing "fountains", novel elements of chromosome organization, emerging selectively at enhancers upon ZGA. Using polymer simulations, we demonstrate that fountains can emerge as sites of targeted cohesin loading and require two-sided, yet desynchronized, loop extrusion. Specific loss of fountains upon loss of pioneer transcription factors that drive ZGA reveals a causal connection between enhancer activity and fountain formation. Finally, we show that fountains emerge in early Medaka and Xenopus embryos; moreover, we found cohesin-dependent fountain pattern on enhancers of mouse embryonic stem cells. Taken together, fountains are the first enhancer-specific elements of chromosome organization; they constitute starting points of chromosome folding during early development, likely serving as sites of targeted cohesin loading.

Ancient gene linkages support ctenophores as sister to other animals

A central question in evolutionary biology is whether sponges or ctenophores (comb jellies) are the sister group to all other animals. These alternative phylogenetic hypotheses imply different scenarios for the evolution of complex neural systems and other animal-specific traits1-6. Conventional phylogenetic approaches based on morphological characters and increasingly extensive gene sequence collections have not been able to definitively answer this question7-11. Here we develop chromosome-scale gene linkage, also known as synteny, as a phylogenetic character for resolving this question12. We report new chromosome-scale genomes for a ctenophore and two marine sponges, and for three unicellular relatives of animals (a choanoflagellate, a filasterean amoeba and an ichthyosporean) that serve as outgroups for phylogenetic analysis. We find ancient syntenies that are conserved between animals and their close unicellular relatives. Ctenophores and unicellular eukaryotes share ancestral metazoan patterns, whereas sponges, bilaterians, and cnidarians share derived chromosomal rearrangements. Conserved syntenic characters unite sponges with bilaterians, cnidarians, and placozoans in a monophyletic clade to the exclusion of ctenophores, placing ctenophores as the sister group to all other animals. The patterns of synteny shared by sponges, bilaterians, and cnidarians are the result of rare and irreversible chromosome fusion-and-mixing events that provide robust and unambiguous phylogenetic support for the ctenophore-sister hypothesis. These findings provide a new framework for resolving deep, recalcitrant phylogenetic problems and have implications for our understanding of animal evolution.