Three-dimensional genome architecture persists in a 52,000-year-old woolly mammoth skin sample.

Analyses of ancient DNA typically involve sequencing the surviving short oligonucleotides and aligning to genome assemblies from related, modern species. Here, we report that skin from a female woolly mammoth (†Mammuthus primigenius) that died 52,000 years ago retained its ancient genome architecture. We use PaleoHi-C to map chromatin contacts and assemble its genome, yielding 28 chromosome-length scaffolds. Chromosome territories, compartments, loops, Barr bodies, and inactive X chromosome (Xi) superdomains persist. The active and inactive genome compartments in mammoth skin more closely resemble Asian elephant skin than other elephant tissues. Our analyses uncover new biology. Differences in compartmentalization reveal genes whose transcription was potentially altered in mammoths vs. elephants. Mammoth Xi has a tetradic architecture, not bipartite like human and mouse. We hypothesize that, shortly after this mammoth's death, the sample spontaneously freeze-dried in the Siberian cold, leading to a glass transition that preserved subfossils of ancient chromosomes at nanometer scale.

Perivascular neurons instruct 3D vascular lattice formation via neurovascular contact.

The vasculature of the central nervous system is a 3D lattice composed of laminar vascular beds interconnected by penetrating vessels. The mechanisms controlling 3D lattice network formation remain largely unknown. Combining viral labeling, genetic marking, and single-cell profiling in the mouse retina, we discovered a perivascular neuronal subset, annotated as Fam19a4/Nts-positive retinal ganglion cells (Fam19a4/Nts-RGCs), directly contacting the vasculature with perisomatic endfeet. Developmental ablation of Fam19a4/Nts-RGCs led to disoriented growth of penetrating vessels near the ganglion cell layer (GCL), leading to a disorganized 3D vascular lattice. We identified enriched PIEZO2 expression in Fam19a4/Nts-RGCs. Piezo2 loss from all retinal neurons or Fam19a4/Nts-RGCs abolished the direct neurovascular contacts and phenocopied the Fam19a4/Nts-RGC ablation deficits. The defective vascular structure led to reduced capillary perfusion and sensitized the retina to ischemic insults. Furthermore, we uncovered a Piezo2-dependent perivascular granule cell subset for cerebellar vascular patterning, indicating neuronal Piezo2-dependent 3D vascular patterning in the brain.

Structure-guided discovery of protein and glycan components in native mastigonemes.

Mastigonemes, the hair-like lateral appendages lining cilia or flagella, participate in mechanosensation and cellular motion, but their constituents and structure have remained unclear. Here, we report the cryo-EM structure of native mastigonemes isolated from Chlamydomonas at 3.0 Å resolution. The long stem assembles as a super spiral, with each helical turn comprising four pairs of anti-parallel mastigoneme-like protein 1 (Mst1). A large array of arabinoglycans, which represents a common class of glycosylation in plants and algae, is resolved surrounding the type II poly-hydroxyproline (Hyp) helix in Mst1. The EM map unveils a mastigoneme axial protein (Mstax) that is rich in heavily glycosylated Hyp and contains a PKD2-like transmembrane domain (TMD). Mstax, with nearly 8,000 residues spanning from the intracellular region to the distal end of the mastigoneme, provides the framework for Mst1 assembly. Our study provides insights into the complexity of protein and glycan interactions in native bio-architectures.

Chromosome-level organization of the regulatory genome in the Drosophila nervous system.

Previous studies have identified topologically associating domains (TADs) as basic units of genome organization. We present evidence of a previously unreported level of genome folding, where distant TAD pairs, megabases apart, interact to form meta-domains. Within meta-domains, gene promoters and structural intergenic elements present in distant TADs are specifically paired. The associated genes encode neuronal determinants, including those engaged in axonal guidance and adhesion. These long-range associations occur in a large fraction of neurons but support transcription in only a subset of neurons. Meta-domains are formed by diverse transcription factors that are able to pair over long and flexible distances. We present evidence that two such factors, GAF and CTCF, play direct roles in this process. The relative simplicity of higher-order meta-domain interactions in Drosophila, compared with those previously described in mammals, allowed the demonstration that genomes can fold into highly specialized cell-type-specific scaffolds that enable megabase-scale regulatory associations.

Modeling epigenetic lesions that cause gliomas.

Epigenetic lesions that disrupt regulatory elements represent potential cancer drivers. However, we lack experimental models for validating their tumorigenic impact. Here, we model aberrations arising in isocitrate dehydrogenase-mutant gliomas, which exhibit DNA hypermethylation. We focus on a CTCF insulator near the PDGFRA oncogene that is recurrently disrupted by methylation in these tumors. We demonstrate that disruption of the syntenic insulator in mouse oligodendrocyte progenitor cells (OPCs) allows an OPC-specific enhancer to contact and induce Pdgfra, thereby increasing proliferation. We show that a second lesion, methylation-dependent silencing of the Cdkn2a tumor suppressor, cooperates with insulator loss in OPCs. Coordinate inactivation of the Pdgfra insulator and Cdkn2a drives gliomagenesis in vivo. Despite locus synteny, the insulator is CpG-rich only in humans, a feature that may confer human glioma risk but complicates mouse modeling. Our study demonstrates the capacity of recurrent epigenetic lesions to drive OPC proliferation in vitro and gliomagenesis in vivo.

Non-cell-autonomous disruption of nuclear architecture as a potential cause of COVID-19-induced anosmia.

SARS-CoV-2 infects less than 1% of cells in the human body, yet it can cause severe damage in a variety of organs. Thus, deciphering the non-cell-autonomous effects of SARS-CoV-2 infection is imperative for understanding the cellular and molecular disruption it elicits. Neurological and cognitive defects are among the least understood symptoms of COVID-19 patients, with olfactory dysfunction being their most common sensory deficit. Here, we show that both in humans and hamsters, SARS-CoV-2 infection causes widespread downregulation of olfactory receptors (ORs) and of their signaling components. This non-cell-autonomous effect is preceded by a dramatic reorganization of the neuronal nuclear architecture, which results in dissipation of genomic compartments harboring OR genes. Our data provide a potential mechanism by which SARS-CoV-2 infection alters the cellular morphology and the transcriptome of cells it cannot infect, offering insight to its systemic effects in olfaction and beyond.

Cohesin mediates DNA loop extrusion by a "swing and clamp" mechanism.

Structural maintenance of chromosomes (SMC) complexes organize genome topology in all kingdoms of life and have been proposed to perform this function by DNA loop extrusion. How this process works is unknown. Here, we have analyzed how loop extrusion is mediated by human cohesin-NIPBL complexes, which enable chromatin folding in interphase cells. We have identified DNA binding sites and large-scale conformational changes that are required for loop extrusion and have determined how these are coordinated. Our results suggest that DNA is translocated by a spontaneous 50 nm-swing of cohesin's hinge, which hands DNA over to the ATPase head of SMC3, where upon binding of ATP, DNA is clamped by NIPBL. During this process, NIPBL "jumps ship" from the hinge toward the SMC3 head and might thereby couple the spontaneous hinge swing to ATP-dependent DNA clamping. These results reveal mechanistic principles of how cohesin-NIPBL and possibly other SMC complexes mediate loop extrusion.

Deciphering osteoarthritis genetics across 826,690 individuals from 9 populations.

Osteoarthritis affects over 300 million people worldwide. Here, we conduct a genome-wide association study meta-analysis across 826,690 individuals (177,517 with osteoarthritis) and identify 100 independently associated risk variants across 11 osteoarthritis phenotypes, 52 of which have not been associated with the disease before. We report thumb and spine osteoarthritis risk variants and identify differences in genetic effects between weight-bearing and non-weight-bearing joints. We identify sex-specific and early age-at-onset osteoarthritis risk loci. We integrate functional genomics data from primary patient tissues (including articular cartilage, subchondral bone, and osteophytic cartilage) and identify high-confidence effector genes. We provide evidence for genetic correlation with phenotypes related to pain, the main disease symptom, and identify likely causal genes linked to neuronal processes. Our results provide insights into key molecular players in disease processes and highlight attractive drug targets to accelerate translation.

Senescence: An Identity Crisis Originating from Deep Within the Nucleus.

Cellular senescence is implicated in a wide range of physiological and pathological conditions throughout an organism's entire lifetime. In particular, it has become evident that senescence plays a causative role in aging and age-associated disorders. This is not due simply to the loss of function of senescent cells. Instead, the substantial alterations of the cellular activities of senescent cells, especially the array of secretory factors, impact the surrounding tissues or even entire organisms. Such non-cell-autonomous functionality is largely coordinated by tissue-specific genes, constituting a cell fate-determining state. Senescence can be viewed as a gain-of-function phenotype or a process of cell identity shift. Cellular functionality or lineage-specific gene expression is tightly linked to the cell type-specific epigenetic landscape, reinforcing the heterogeneity of senescence across cell types. Here, we aim to define the senescence cellular functionality and epigenetic features that may contribute to the gain-of-function phenotype.

Large-Scale Topological Changes Restrain Malignant Progression in Colorectal Cancer.

Widespread changes to DNA methylation and chromatin are well documented in cancer, but the fate of higher-order chromosomal structure remains obscure. Here we integrated topological maps for colon tumors and normal colons with epigenetic, transcriptional, and imaging data to characterize alterations to chromatin loops, topologically associated domains, and large-scale compartments. We found that spatial partitioning of the open and closed genome compartments is profoundly compromised in tumors. This reorganization is accompanied by compartment-specific hypomethylation and chromatin changes. Additionally, we identify a compartment at the interface between the canonical A and B compartments that is reorganized in tumors. Remarkably, similar shifts were evident in non-malignant cells that have accumulated excess divisions. Our analyses suggest that these topological changes repress stemness and invasion programs while inducing anti-tumor immunity genes and may therefore restrain malignant progression. Our findings call into question the conventional view that tumor-associated epigenomic alterations are primarily oncogenic.

The Self-Organizing Genome: Principles of Genome Architecture and Function.

Genomes have complex three-dimensional architectures. The recent convergence of genetic, biochemical, biophysical, and cell biological methods has uncovered several fundamental principles of genome organization. They highlight that genome function is a major driver of genome architecture and that structural features of chromatin act as modulators, rather than binary determinants, of genome activity. The interplay of these principles in the context of self-organization can account for the emergence of structural chromatin features, the diversity and single-cell heterogeneity of nuclear architecture in cell types and tissues, and explains evolutionarily conserved functional features of genomes, including plasticity and robustness.

Trans-ethnic and Ancestry-Specific Blood-Cell Genetics in 746,667 Individuals from 5 Global Populations.

Most loci identified by GWASs have been found in populations of European ancestry (EUR). In trans-ethnic meta-analyses for 15 hematological traits in 746,667 participants, including 184,535 non-EUR individuals, we identified 5,552 trait-variant associations at p < 5 × 10-9, including 71 novel associations not found in EUR populations. We also identified 28 additional novel variants in ancestry-specific, non-EUR meta-analyses, including an IL7 missense variant in South Asians associated with lymphocyte count in vivo and IL-7 secretion levels in vitro. Fine-mapping prioritized variants annotated as functional and generated 95% credible sets that were 30% smaller when using the trans-ethnic as opposed to the EUR-only results. We explored the clinical significance and predictive value of trans-ethnic variants in multiple populations and compared genetic architecture and the effect of natural selection on these blood phenotypes between populations. Altogether, our results for hematological traits highlight the value of a more global representation of populations in genetic studies.

Heterochromatin-Driven Nuclear Softening Protects the Genome against Mechanical Stress-Induced Damage.

Tissue homeostasis requires maintenance of functional integrity under stress. A central source of stress is mechanical force that acts on cells, their nuclei, and chromatin, but how the genome is protected against mechanical stress is unclear. We show that mechanical stretch deforms the nucleus, which cells initially counteract via a calcium-dependent nuclear softening driven by loss of H3K9me3-marked heterochromatin. The resulting changes in chromatin rheology and architecture are required to insulate genetic material from mechanical force. Failure to mount this nuclear mechanoresponse results in DNA damage. Persistent, high-amplitude stretch induces supracellular alignment of tissue to redistribute mechanical energy before it reaches the nucleus. This tissue-scale mechanoadaptation functions through a separate pathway mediated by cell-cell contacts and allows cells/tissues to switch off nuclear mechanotransduction to restore initial chromatin state. Our work identifies an unconventional role of chromatin in altering its own mechanical state to maintain genome integrity in response to deformation.

Coordinated Cellular Neighborhoods Orchestrate Antitumoral Immunity at the Colorectal Cancer Invasive Front.

Antitumoral immunity requires organized, spatially nuanced interactions between components of the immune tumor microenvironment (iTME). Understanding this coordinated behavior in effective versus ineffective tumor control will advance immunotherapies. We re-engineered co-detection by indexing (CODEX) for paraffin-embedded tissue microarrays, enabling simultaneous profiling of 140 tissue regions from 35 advanced-stage colorectal cancer (CRC) patients with 56 protein markers. We identified nine conserved, distinct cellular neighborhoods (CNs)-a collection of components characteristic of the CRC iTME. Enrichment of PD-1+CD4+ T cells only within a granulocyte CN positively correlated with survival in a high-risk patient subset. Coupling of tumor and immune CNs, fragmentation of T cell and macrophage CNs, and disruption of inter-CN communication was associated with inferior outcomes. This study provides a framework for interrogating how complex biological processes, such as antitumoral immunity, occur through concerted actions of cells and spatial domains.

Identifying cis Elements for Spatiotemporal Control of Mammalian DNA Replication.

The temporal order of DNA replication (replication timing [RT]) is highly coupled with genome architecture, but cis-elements regulating either remain elusive. We created a series of CRISPR-mediated deletions and inversions of a pluripotency-associated topologically associating domain (TAD) in mouse ESCs. CTCF-associated domain boundaries were dispensable for RT. CTCF protein depletion weakened most TAD boundaries but had no effect on RT or A/B compartmentalization genome-wide. By contrast, deletion of three intra-TAD CTCF-independent 3D contact sites caused a domain-wide early-to-late RT shift, an A-to-B compartment switch, weakening of TAD architecture, and loss of transcription. The dispensability of TAD boundaries and the necessity of these "early replication control elements" (ERCEs) was validated by deletions and inversions at additional domains. Our results demonstrate that discrete cis-regulatory elements orchestrate domain-wide RT, A/B compartmentalization, TAD architecture, and transcription, revealing fundamental principles linking genome structure and function.

Stromal Microenvironment Shapes the Intratumoral Architecture of Pancreatic Cancer.

Single-cell technologies have described heterogeneity across tissues, but the spatial distribution and forces that drive single-cell phenotypes have not been well defined. Combining single-cell RNA and protein analytics in studying the role of stromal cancer-associated fibroblasts (CAFs) in modulating heterogeneity in pancreatic cancer (pancreatic ductal adenocarcinoma [PDAC]) model systems, we have identified significant single-cell population shifts toward invasive epithelial-to-mesenchymal transition (EMT) and proliferative (PRO) phenotypes linked with mitogen-activated protein kinase (MAPK) and signal transducer and activator of transcription 3 (STAT3) signaling. Using high-content digital imaging of RNA in situ hybridization in 195 PDAC tumors, we quantified these EMT and PRO subpopulations in 319,626 individual cancer cells that can be classified within the context of distinct tumor gland "units." Tumor gland typing provided an additional layer of intratumoral heterogeneity that was associated with differences in stromal abundance and clinical outcomes. This demonstrates the impact of the stroma in shaping tumor architecture by altering inherent patterns of tumor glands in human PDAC.

Trans Effects on Gene Expression Can Drive Omnigenic Inheritance.

Early genome-wide association studies (GWASs) led to the surprising discovery that, for typical complex traits, most of the heritability is due to huge numbers of common variants with tiny effect sizes. Previously, we argued that new models are needed to understand these patterns. Here, we provide a formal model in which genetic contributions to complex traits are partitioned into direct effects from core genes and indirect effects from peripheral genes acting in trans. We propose that most heritability is driven by weak trans-eQTL SNPs, whose effects are mediated through peripheral genes to impact the expression of core genes. In particular, if the core genes for a trait tend to be co-regulated, then the effects of peripheral variation can be amplified such that nearly all of the genetic variance is driven by weak trans effects. Thus, our model proposes a framework for understanding key features of the architecture of complex traits.

Genomic Relationships, Novel Loci, and Pleiotropic Mechanisms across Eight Psychiatric Disorders.

Genetic influences on psychiatric disorders transcend diagnostic boundaries, suggesting substantial pleiotropy of contributing loci. However, the nature and mechanisms of these pleiotropic effects remain unclear. We performed analyses of 232,964 cases and 494,162 controls from genome-wide studies of anorexia nervosa, attention-deficit/hyperactivity disorder, autism spectrum disorder, bipolar disorder, major depression, obsessive-compulsive disorder, schizophrenia, and Tourette syndrome. Genetic correlation analyses revealed a meaningful structure within the eight disorders, identifying three groups of inter-related disorders. Meta-analysis across these eight disorders detected 109 loci associated with at least two psychiatric disorders, including 23 loci with pleiotropic effects on four or more disorders and 11 loci with antagonistic effects on multiple disorders. The pleiotropic loci are located within genes that show heightened expression in the brain throughout the lifespan, beginning prenatally in the second trimester, and play prominent roles in neurodevelopmental processes. These findings have important implications for psychiatric nosology, drug development, and risk prediction.

Root System Depth in Arabidopsis Is Shaped by EXOCYST70A3 via the Dynamic Modulation of Auxin Transport.

Root system architecture (RSA), the distribution of roots in soil, plays a major role in plant survival. RSA is shaped by multiple developmental processes that are largely governed by the phytohormone auxin, suggesting that auxin regulates responses of roots that are important for local adaptation. However, auxin has a central role in numerous processes, and it is unclear which molecular mechanisms contribute to the variation in RSA for environmental adaptation. Using natural variation in Arabidopsis, we identify EXOCYST70A3 as a modulator of the auxin system that causes variation in RSA by acting on PIN4 protein distribution. Allelic variation and genetic perturbation of EXOCYST70A3 lead to alteration of root gravitropic responses, resulting in a different RSA depth profile and drought resistance. Overall our findings suggest that the local modulation of the pleiotropic auxin pathway can gives rise to distinct RSAs that can be adaptive in specific environments.

The ChAHP Complex Counteracts Chromatin Looping at CTCF Sites that Emerged from SINE Expansions in Mouse.

CCCTC-binding factor (CTCF) and cohesin are key players in three-dimensional chromatin organization. The topologically associating domains (TADs) demarcated by CTCF are remarkably well conserved between species, although genome-wide CTCF binding has diverged substantially following transposon-mediated motif expansions. Therefore, the CTCF consensus motif poorly predicts TADs, and additional factors must modulate CTCF binding and subsequent TAD formation. Here, we demonstrate that the ChAHP complex (CHD4, ADNP, HP1) competes with CTCF for a common set of binding motifs. In Adnp knockout cells, novel insulated regions are formed at sites normally bound by ChAHP, whereas proximal canonical boundaries are weakened. These data reveal that CTCF-mediated loop formation is modulated by a distinct zinc-finger protein complex. Strikingly, ChAHP-bound loci are mainly situated within less diverged SINE B2 transposable elements. This implicates ChAHP in maintenance of evolutionarily conserved spatial chromatin organization by buffering novel CTCF binding sites that emerged through SINE expansions.