Disruptions of topological chromatin domains cause pathogenic rewiring of gene-enhancer interactions.

Mammalian genomes are organized into megabase-scale topologically associated domains (TADs). We demonstrate that disruption of TADs can rewire long-range regulatory architecture and result in pathogenic phenotypes. We show that distinct human limb malformations are caused by deletions, inversions, or duplications altering the structure of the TAD-spanning WNT6/IHH/EPHA4/PAX3 locus. Using CRISPR/Cas genome editing, we generated mice with corresponding rearrangements. Both in mouse limb tissue and patient-derived fibroblasts, disease-relevant structural changes cause ectopic interactions between promoters and non-coding DNA, and a cluster of limb enhancers normally associated with Epha4 is misplaced relative to TAD boundaries and drives ectopic limb expression of another gene in the locus. This rewiring occurred only if the variant disrupted a CTCF-associated boundary domain. Our results demonstrate the functional importance of TADs for orchestrating gene expression via genome architecture and indicate criteria for predicting the pathogenicity of human structural variants, particularly in non-coding regions of the human genome.

The little skate genome and the evolutionary emergence of wing-like fins.

Skates are cartilaginous fish whose body plan features enlarged wing-like pectoral fins, enabling them to thrive in benthic environments1,2. However, the molecular underpinnings of this unique trait remain unclear. Here we investigate the origin of this phenotypic innovation by developing the little skate Leucoraja erinacea as a genomically enabled model. Analysis of a high-quality chromosome-scale genome sequence for the little skate shows that it preserves many ancestral jawed vertebrate features compared with other sequenced genomes, including numerous ancient microchromosomes. Combining genome comparisons with extensive regulatory datasets in developing fins-including gene expression, chromatin occupancy and three-dimensional conformation-we find skate-specific genomic rearrangements that alter the three-dimensional regulatory landscape of genes that are involved in the planar cell polarity pathway. Functional inhibition of planar cell polarity signalling resulted in a reduction in anterior fin size, confirming that this pathway is a major contributor to batoid fin morphology. We also identified a fin-specific enhancer that interacts with several hoxa genes, consistent with the redeployment of hox gene expression in anterior pectoral fins, and confirmed its potential to activate transcription in the anterior fin using zebrafish reporter assays. Our findings underscore the central role of genome reorganization and regulatory variation in the evolution of phenotypes, shedding light on the molecular origin of an enigmatic trait.

Structural variation in the 3D genome.

Structural and quantitative chromosomal rearrangements, collectively referred to as structural variation (SV), contribute to a large extent to the genetic diversity of the human genome and thus are of high relevance for cancer genetics, rare diseases and evolutionary genetics. Recent studies have shown that SVs can not only affect gene dosage but also modulate basic mechanisms of gene regulation. SVs can alter the copy number of regulatory elements or modify the 3D genome by disrupting higher-order chromatin organization such as topologically associating domains. As a result of these position effects, SVs can influence the expression of genes distant from the SV breakpoints, thereby causing disease. The impact of SVs on the 3D genome and on gene expression regulation has to be considered when interpreting the pathogenic potential of these variant types.

Cell adhesion and immune response, two main functions altered in the transcriptome of seasonally regressed testes of two mammalian species.

In species with seasonal breeding, male specimens undergo substantial testicular regression during the nonbreeding period of the year. However, the molecular mechanisms that control this biological process are largely unknown. Here, we report a transcriptomic analysis on the Iberian mole, Talpa occidentalis, in which the desquamation of live, nonapoptotic germ cells is the major cellular event responsible for testis regression. By comparing testes at different reproductive states (active, regressing, and inactive), we demonstrate that the molecular pathways controlling the cell adhesion function in the seminiferous epithelium, such as the MAPK, ERK, and TGF-ß signaling, are altered during the regression process. In addition, inactive testes display a global upregulation of genes associated with immune response, indicating a selective loss of the "immune privilege" that normally operates in sexually active testes. Interspecies comparative analyses using analogous data from the Mediterranean pine vole, a rodent species where testis regression is controlled by halting meiosis entry, revealed a common gene expression signature in the regressed testes of these two evolutionary distant species. Our study advances in the knowledge of the molecular mechanisms associated to gonadal seasonal breeding, highlighting the existence of a conserved transcriptional program of testis involution across mammalian clades.

Unraveling the transcriptional regulation of TWIST1 in limb development.

The transcription factor TWIST1 plays a vital role in mesoderm development, particularly in limb and craniofacial formation. Accordingly, haploinsufficiency of TWIST1 can cause limb and craniofacial malformations as part of Saethre-Chotzen syndrome. However, the molecular basis of TWIST1 transcriptional regulation during development has yet to be elucidated. Here, we characterized active enhancers in the TWIST1-HDAC9 locus that drive transcription in the developing limb and branchial arches. Using available p300 and H3K27ac ChIP-seq data, we identified 12 enhancer candidates, located both within and outside the coding sequences of the neighboring gene, Histone deacetyase 9 (HDAC9). Using zebrafish and mouse enhancer assays, we showed that eight of these candidates have limb/fin and branchial arch enhancer activity that resemble Twist1 expression. Using 4C-seq, we showed that the Twist1 promoter region interacts with three enhancers (eTw-5, 6, 7) in the limb bud and branchial arch of mouse embryos at day 11.5. Furthermore, we found that two transcription factors, LMX1B and TFAP2, bind these enhancers and modulate their enhancer activity. Finally, using CRISPR/Cas9 genome editing, we showed that homozygous deletion of eTw5-7 enhancers reduced Twist1 expression in the limb bud and caused pre-axial polydactyly, a phenotype observed in Twist1+/- mice. Taken together, our findings reveal that each enhancer has a discrete activity pattern, and together comprise a spatiotemporal regulatory network of Twist1 transcription in the developing limbs/fins and branchial arches. Our study suggests that mutations in TWIST1 enhancers could lead to reduced TWIST1 expression, resulting in phenotypic outcome as seen with TWIST1 coding mutations.

GOPHER: Generator Of Probes for capture Hi-C Experiments at high Resolution.

BACKGROUND: Target enrichment combined with chromosome conformation capturing methodologies such as capture Hi-C (CHC) can be used to investigate spatial layouts of genomic regions with high resolution and at scalable costs. A common application of CHC is the investigation of regulatory elements that are in contact with promoters, but CHC can be used for a range of other applications. Therefore, probe design for CHC needs to be adapted to experimental needs, but no flexible tool is currently available for this purpose. RESULTS: We present a Java desktop application called GOPHER (Generator Of Probes for capture Hi-C Experiments at high Resolution) that implements three strategies for CHC probe design. GOPHER's simple approach is similar to the probe design of previous approaches that employ CHC to investigate all promoters, with one probe being placed at each margin of a single digest that overlaps the transcription start site (TSS) of each promoter. GOPHER's simple-patched approach extends this methodology with a heuristic that improves coverage of viewpoints in which the TSS is located near to one of the boundaries of the digest. GOPHER's extended approach is intended mainly for focused investigations of smaller gene sets. GOPHER can also be used to design probes for regions other than TSS such as GWAS hits or large blocks of genomic sequence. GOPHER additionally provides a number of features that allow users to visualize and edit viewpoints, and outputs a range of files useful for documentation, ordering probes, and downstream analysis. CONCLUSION: GOPHER is an easy-to-use and robust desktop application for CHC probe design. Source code and a precompiled executable can be downloaded from the GOPHER GitHub page at https://github.com/TheJacksonLaboratory/Gopher .

Breaking TADs: How Alterations of Chromatin Domains Result in Disease.

Spatial organization is an inherent property of the vertebrate genome to accommodate the roughly 2m of DNA in the nucleus of a cell. In this nonrandom organization, topologically associating domains (TADs) emerge as a fundamental structural unit that is thought to guide regulatory elements to their cognate promoters. In this review we summarize the most recent findings about TADs and the boundary regions separating them. We discuss how the disruption of these structures by genomic rearrangements can result in gene misexpression and disease.

Exome sequencing and CRISPR/Cas genome editing identify mutations of ZAK as a cause of limb defects in humans and mice.

The CRISPR/Cas technology enables targeted genome editing and the rapid generation of transgenic animal models for the study of human genetic disorders. Here we describe an autosomal recessive human disease in two unrelated families characterized by a split-foot defect, nail abnormalities of the hands, and hearing loss, due to mutations disrupting the SAM domain of the protein kinase ZAK. ZAK is a member of the MAPKKK family with no known role in limb development. We show that Zak is expressed in the developing limbs and that a CRISPR/Cas-mediated knockout of the two Zak isoforms is embryonically lethal in mice. In contrast, a deletion of the SAM domain induces a complex hindlimb defect associated with down-regulation of Trp63, a known split-hand/split-foot malformation disease gene. Our results identify ZAK as a key player in mammalian limb patterning and demonstrate the rapid utility of CRISPR/Cas genome editing to assign causality to human mutations in the mouse in <10 wk.

Complete male-to-female sex reversal in XY mice lacking the miR-17~92 cluster.

Mammalian sex determination is controlled by antagonistic gene cascades operating in embryonic undifferentiated gonads. The expression of the Y-linked gene SRY is sufficient to trigger the testicular pathway, whereas its absence in XX embryos leads to ovarian differentiation. Yet, the potential involvement of non-coding regulation in this process remains unclear. Here we show that the deletion of a single microRNA cluster, miR-17~92, induces complete primary male-to-female sex reversal in XY mice. Sry expression is delayed in XY knockout gonads, which develop as ovaries. Sertoli cell differentiation is reduced, delayed and unable to sustain testicular development. Pre-supporting cells in mutant gonads undergo a transient state of sex ambiguity which is subsequently resolved towards the ovarian fate. The miR-17~92 predicted target genes are upregulated, affecting the fine regulation of gene networks controlling gonad development. Thus, microRNAs emerge as key components for mammalian sex determination, controlling Sry expression timing and Sertoli cell differentiation.

Cancer Evolution: A Multifaceted Affair.

Cancer cells adapt and survive through the acquisition and selection of molecular modifications. This process defines cancer evolution. Building on a theoretical framework based on heritable genetic changes has provided insights into the mechanisms supporting cancer evolution. However, cancer hallmarks also emerge via heritable nongenetic mechanisms, including epigenetic and chromatin topological changes, and interactions between tumor cells and the tumor microenvironment. Recent findings on tumor evolutionary mechanisms draw a multifaceted picture where heterogeneous forces interact and influence each other while shaping tumor progression. A comprehensive characterization of the cancer evolutionary toolkit is required to improve personalized medicine and biomarker discovery. SIGNIFICANCE: Tumor evolution is fueled by multiple enabling mechanisms. Importantly, genetic instability, epigenetic reprogramming, and interactions with the tumor microenvironment are neither alternative nor independent evolutionary mechanisms. As demonstrated by findings highlighted in this perspective, experimental and theoretical approaches must account for multiple evolutionary mechanisms and their interactions to ultimately understand, predict, and steer tumor evolution.

A microRNA (mmu-miR-124) prevents Sox9 expression in developing mouse ovarian cells.

In mammals, sex differentiation depends on gonad development, which is controlled by two groups of sex-determining genes that promote one gonadal sex and antagonize the opposite one. SOX9 plays a key role during testis development in all studied vertebrates, whereas it is kept inactive in the XX gonad at the critical time of sex determination, otherwise, ovary-to-testis gonadal sex reversal occurs. However, molecular mechanisms underlying repression of Sox9 at the beginning of ovarian development, as well as other important aspects of gonad organogenesis, remain largely unknown. Because there is indirect evidence that micro-RNAs (miRNA) are necessary for testicular function, the possible involvement of miRNAs in mammalian sex determination deserved further research. Using microarray technology, we have identified 22 miRNAs showing sex-specific expression in the developing gonads during the critical period of sex determination. Bioinformatics analyses led to the identification of miR-124 as the candidate gene for ovarian development. We knocked down or overexpressed miR-124 in primary gonadal cell cultures and observed that miR-124 is sufficient to induce the repression of both SOX9 translation and transcription in ovarian cells. Our results provide the first evidence of the involvement of a miRNA in the regulation of the gene controlling gonad development and sex determination. The miRNA microarray data reported here will help promote further research in this field, to unravel the role of other miRNAs in the genetic control of mammalian sex determination.

Identification of live germ-cell desquamation as a major mechanism of seasonal testis regression in mammals: a study in the Iberian mole (Talpa occidentalis).

In males of seasonally breeding species, testes undergo a severe involution at the end of the breeding season, with a major volume decrease due to massive germ-cell depletion associated with photoperiod-dependent reduced levels of testosterone and gonadotropins. Although it has been repeatedly suggested that apoptosis is the principal effector of testicular regression in vertebrates, recent studies do not support this hypothesis in some mammals. The purpose of our work is to discover alternative mechanisms of testis regression in these species. In this paper, we have performed a morphological, hormonal, ultrastructural, molecular, and functional study of the mechanism of testicular regression and the role that cell junctions play in the cell-content dynamics of the testis of the Iberian mole, Talpa occidentalis, throughout the seasonal breeding cycle. Desquamation of live, nonapoptotic germ cells has been identified here as a new mechanism for seasonal testis involution in mammals, indicating that testis regression is regulated by modulating the expression and distribution of the cell-adhesion molecules in the seminiferous epithelium. During this process, which is mediated by low intratesticular testosterone levels, Sertoli cells lose their nursing and supporting function, as well as the impermeability of the blood-testis barrier. Our results contradict the current paradigm that apoptosis is the major testis regression effector in vertebrates, as it is clearly not true in all mammals. The new testis regression mechanism described here for the mole could then be generalized to other mammalian species. Available data from some previously studied mammals should be reevaluated.

Evolution of 3D chromatin organization at different scales.

Most animal genomes fold in 3D chromatin domains called topologically associated domains (TADs) that facilitate interactions between cis-regulatory elements (CREs) and promoters. Owing to their critical role in the control of developmental gene expression, we explore how TADs have shaped animal evolution. In the light of recent studies that profile TADs in disparate animal lineages, we discuss their phylogenetic distribution and the mechanisms that underlie their formation. We present evidence indicating that TADs are plastic entities composed of genomic strata of different ages: ancient cores are combined with newer regions and brought into extant TADs through genomic rearrangements. We highlight that newly incorporated TAD strata enable the establishment of new CRE-promoter interactions and in turn new expression patterns that can drive phenotypical innovation. We further highlight how subtle changes in chromatin folding may fine-tune the expression levels of developmental genes and hold a potential for evolutionary significance.

Dynamic antagonism between key repressive pathways maintains the placental epigenome.

DNA and Histone 3 Lysine 27 methylation typically function as repressive modifications and operate within distinct genomic compartments. In mammals, the majority of the genome is kept in a DNA methylated state, whereas the Polycomb repressive complexes regulate the unmethylated CpG-rich promoters of developmental genes. In contrast to this general framework, the extra-embryonic lineages display non-canonical, globally intermediate DNA methylation levels, including disruption of local Polycomb domains. Here, to better understand this unusual landscape's molecular properties, we genetically and chemically perturbed major epigenetic pathways in mouse trophoblast stem cells. We find that the extra-embryonic epigenome reflects ongoing and dynamic de novo methyltransferase recruitment, which is continuously antagonized by Polycomb to maintain intermediate, locally disordered methylation. Despite its disorganized molecular appearance, our data point to a highly controlled equilibrium between counteracting repressors within extra-embryonic cells, one that can seemingly persist indefinitely without bistable features typically seen for embryonic forms of epigenetic regulation.

Pattern and density of vascularization in mammalian testes, ovaries, and ovotestes.

According to the classical paradigm, the vasculature of the embryonic testis is more dense and complex than that of the ovary, but recent studies based on whole-mount detection of Caveolin-1 (CAV1) as an endothelial cell marker, have suggested that the level of ovarian vascularization is higher than previously assumed. However, this new hypothesis has been neither tested using alternative methodology nor investigated in other mammalian species. In this paper, we have studied the vascularization process in the gonads of males and females of two mammalian species, the mouse (Mus musculus) and the Iberian mole (Talpa occidentalis). Our results show that the pattern of testis vascularization is very well conserved among mammals, including both pre- and postnatal stages of development and, at least in the mole, it is conserved irrespectively of whether the testicular tissue is XY or XX. We have shown that CAV1 is present not only in endothelial cells but also in prefollicular oocytes and in an ovarian population of somatic cortical cells. These data clearly establish that: (1) according to the classical hypothesis, the degree of vascularization of the developing ovary is lower than that of the testis, (2) ovarian vascularization is also evolutionarily conserved as it occurs similarly both in moles and in mice, and (3) that the degree of vascular development of the mammalian ovary is age-dependent increasing significatively at puberty. The expression of CAV1 in the ovary of most animal taxa, from nematodes to mammals, strongly suggests a role for this gene in the female meiosis.

Morphology and ultrastructure of the chorioallantoic placenta of the Iberian mole (Talpa occidentalis) with special reference to heterophagous areolas and the nature of interhaemal barrier.

This study provides a contribution to the reconstruction of the eulipotyphlan placental morphotype and also may help resolving a long-standing conflict about the interhaemal barrier in moles. As detailed descriptions of talpid placentation, only available for Talpa europaea and Scalopus aquaticus, led to a controversial debate on the nature of interhaemal barrier, the collection of more placental data of further mole species was strongly desired. Hence, the placentas of six gestational stages of Talpa occidentalis have been studied concerning their morphogenesis and ultrastructure with special focus on the structure of the interhaemal barrier and heterophagous regions. Generally, the mode of placentation in T. occidentalis resembles that of T. europaea, including a broad, discoid, antimesometrial, definitive chorioallantoic placenta of labyrinthine type being still villous in earlier stages. Within the labyrinth, the zona intima shows an endotheliochorial interhaemal barrier with a two-layered trophoblast. This clearly contradicts former statements on the S. aquaticus placenta made by Prasad et al. (1979), although their findings cannot exclude a totally different interpretation. Regardless, the placenta of moles represents the least invasive mode of placentation among Eulipotyphla, which otherwise have highly invasive placentas. Although the phagocytic areolas situated above uterine gland openings are heterophagous, they mainly seem to serve fetal histiotrophic nutrition, at least early in pregnancy. In later stages the number of glands and areolas decreases. This special type of additional phagocytic region is usually most common in species with noninvasive, epitheliochorial placentation, which suggests a correlation between placental invasiveness and the occurrence and type of phagocytic placental structures. The compact and invasive mode of placentation of Talpidae and all other Eulipotyphla seems to be plesiomorphic within Laurasiatheria and is always correlated with an altricial neonate.

In vivo dissection of a clustered-CTCF domain boundary reveals developmental principles of regulatory insulation.

Vertebrate genomes organize into topologically associating domains, delimited by boundaries that insulate regulatory elements from nontarget genes. However, how boundary function is established is not well understood. Here, we combine genome-wide analyses and transgenic mouse assays to dissect the regulatory logic of clustered-CCCTC-binding factor (CTCF) boundaries in vivo, interrogating their function at multiple levels: chromatin interactions, transcription and phenotypes. Individual CTCF binding site (CBS) deletions revealed that the characteristics of specific sites can outweigh other factors such as CBS number and orientation. Combined deletions demonstrated that CBSs cooperate redundantly and provide boundary robustness. We show that divergent CBS signatures are not strictly required for effective insulation and that chromatin loops formed by nonconvergently oriented sites could be mediated by a loop interference mechanism. Further, we observe that insulation strength constitutes a quantitative modulator of gene expression and phenotypes. Our results highlight the modular nature of boundaries and their control over developmental processes.

C/EBPß regulates lipid metabolism and Pparg isoform 2 expression in alveolar macrophages.

Pulmonary alveolar proteinosis (PAP) is a syndrome characterized by accumulation of surfactant lipoproteins within the lung alveoli. Alveolar macrophages (AMs) are crucial for surfactant clearance, and their differentiation depends on colony-stimulating factor 2 (CSF2), which regulates the establishment of an AM-characteristic gene regulatory network. Here, we report that the transcription factor CCAAT/enhancer binding protein ß (C/EBPß) is essential for the development of the AM identity, as demonstrated by transcriptome and chromatin accessibility analysis. Furthermore, C/EBPß-deficient AMs showed severe defects in proliferation, phagocytosis, and lipid metabolism, collectively resulting in a PAP-like syndrome. Mechanistically, the long C/EBPß protein variants LAP* and LAP together with CSF2 signaling induced the expression of Pparg isoform 2 but not Pparg isoform 1, a molecular regulatory mechanism that was also observed in other CSF2-primed macrophages. These results uncover C/EBPß as a key regulator of AM cell fate and shed light on the molecular networks controlling lipid metabolism in macrophages.

Order and disorder: abnormal 3D chromatin organization in human disease.

A precise three-dimensional (3D) organization of chromatin is central to achieve the intricate transcriptional patterns that are required to form complex organisms. Growing evidence supports an important role of 3D chromatin architecture in development and delineates its alterations as prominent causes of disease. In this review, we discuss emerging concepts on the fundamental forces shaping genomes in space and on how their disruption can lead to pathogenic phenotypes. We describe the molecular mechanisms underlying a wide range of diseases, from the systemic effects of coding mutations on 3D architectural factors, to the more tissue-specific phenotypes resulting from genetic and epigenetic modifications at specific loci. Understanding the connection between the 3D organization of the genome and its underlying biological function will allow a better interpretation of human pathogenesis.