Genomic context sensitizes regulatory elements to genetic disruption.

Genomic context critically modulates regulatory function but is difficult to manipulate systematically. The murine insulin-like growth factor 2 (Igf2)/H19 locus is a paradigmatic model of enhancer selectivity, whereby CTCF occupancy at an imprinting control region directs downstream enhancers to activate either H19 or Igf2. We used synthetic regulatory genomics to repeatedly replace the native locus with 157-kb payloads, and we systematically dissected its architecture. Enhancer deletion and ectopic delivery revealed previously uncharacterized long-range regulatory dependencies at the native locus. Exchanging the H19 enhancer cluster with the Sox2 locus control region (LCR) showed that the H19 enhancers relied on their native surroundings while the Sox2 LCR functioned autonomously. Analysis of regulatory DNA actuation across cell types revealed that these enhancer clusters typify broader classes of context sensitivity genome wide. These results show that unexpected dependencies influence even well-studied loci, and our approach permits large-scale manipulation of complete loci to investigate the relationship between regulatory architecture and function.

Synthetic regulatory genomics uncovers enhancer context dependence at the Sox2 locus.

Sox2 expression in mouse embryonic stem cells (mESCs) depends on a distal cluster of DNase I hypersensitive sites (DHSs), but their individual contributions and degree of interdependence remain a mystery. We analyzed the endogenous Sox2 locus using Big-IN to scarlessly integrate large DNA payloads incorporating deletions, rearrangements, and inversions affecting single or multiple DHSs, as well as surgical alterations to transcription factor (TF) recognition sequences. Multiple mESC clones were derived for each payload, sequence-verified, and analyzed for Sox2 expression. We found that two DHSs comprising a handful of key TF recognition sequences were each sufficient for long-range activation of Sox2 expression. By contrast, three nearby DHSs were entirely context dependent, showing no activity alone but dramatically augmenting the activity of the autonomous DHSs. Our results highlight the role of context in modulating genomic regulatory element function, and our synthetic regulatory genomics approach provides a roadmap for the dissection of other genomic loci.

Synthetic reversed sequences reveal default genomic states.

Pervasive transcriptional activity is observed across diverse species. The genomes of extant organisms have undergone billions of years of evolution, making it unclear whether these genomic activities represent effects of selection or 'noise'1-4. Characterizing default genome states could help understand whether pervasive transcriptional activity has biological meaning. Here we addressed this question by introducing a synthetic 101-kb locus into the genomes of Saccharomyces cerevisiae and Mus musculus and characterizing genomic activity. The locus was designed by reversing but not complementing human HPRT1, including its flanking regions, thus retaining basic features of the natural sequence but ablating evolved coding or regulatory information. We observed widespread activity of both reversed and native HPRT1 loci in yeast, despite the lack of evolved yeast promoters. By contrast, the reversed locus displayed no activity at all in mouse embryonic stem cells, and instead exhibited repressive chromatin signatures. The repressive signature was alleviated in a locus variant lacking CpG dinucleotides; nevertheless, this variant was also transcriptionally inactive. These results show that synthetic genomic sequences that lack coding information are active in yeast, but inactive in mouse embryonic stem cells, consistent with a major difference in 'default genomic states' between these two divergent eukaryotic cell types, with implications for understanding pervasive transcription, horizontal transfer of genetic information and the birth of new genes.

On the genetic basis of tail-loss evolution in humans and apes.

The loss of the tail is among the most notable anatomical changes to have occurred along the evolutionary lineage leading to humans and to the 'anthropomorphous apes'1-3, with a proposed role in contributing to human bipedalism4-6. Yet, the genetic mechanism that facilitated tail-loss evolution in hominoids remains unknown. Here we present evidence that an individual insertion of an Alu element in the genome of the hominoid ancestor may have contributed to tail-loss evolution. We demonstrate that this Alu element-inserted into an intron of the TBXT gene7-9-pairs with a neighbouring ancestral Alu element encoded in the reverse genomic orientation and leads to a hominoid-specific alternative splicing event. To study the effect of this splicing event, we generated multiple mouse models that express both full-length and exon-skipped isoforms of Tbxt, mimicking the expression pattern of its hominoid orthologue TBXT. Mice expressing both Tbxt isoforms exhibit a complete absence of the tail or a shortened tail depending on the relative abundance of Tbxt isoforms expressed at the embryonic tail bud. These results support the notion that the exon-skipped transcript is sufficient to induce a tail-loss phenotype. Moreover, mice expressing the exon-skipped Tbxt isoform develop neural tube defects, a condition that affects approximately 1 in 1,000 neonates in humans10. Thus, tail-loss evolution may have been associated with an adaptive cost of the potential for neural tube defects, which continue to affect human health today.

Mouse genome rewriting and tailoring of three important disease loci.

Genetically engineered mouse models (GEMMs) help us to understand human pathologies and develop new therapies, yet faithfully recapitulating human diseases in mice is challenging. Advances in genomics have highlighted the importance of non-coding regulatory genome sequences, which control spatiotemporal gene expression patterns and splicing in many human diseases1,2. Including regulatory extensive genomic regions, which requires large-scale genome engineering, should enhance the quality of disease modelling. Existing methods set limits on the size and efficiency of DNA delivery, hampering the routine creation of highly informative models that we call genomically rewritten and tailored GEMMs (GREAT-GEMMs). Here we describe 'mammalian switching antibiotic resistance markers progressively for integration' (mSwAP-In), a method for efficient genome rewriting in mouse embryonic stem cells. We demonstrate the use of mSwAP-In for iterative genome rewriting of up to 115 kb of a tailored Trp53 locus, as well as for humanization of mice using 116 kb and 180 kb human ACE2 loci. The ACE2 model recapitulated human ACE2 expression patterns and splicing, and notably, presented milder symptoms when challenged with SARS-CoV-2 compared with the existing K18-hACE2 model, thus representing a more human-like model of infection. Finally, we demonstrated serial genome writing by humanizing mouse Tmprss2 biallelically in the ACE2 GREAT-GEMM, highlighting the versatility of mSwAP-In in genome writing.

Genomic context sensitivity of insulator function.

The specificity of interactions between genomic regulatory elements and potential target genes is influenced by the binding of insulator proteins such as CTCF, which can act as potent enhancer blockers when interposed between an enhancer and a promoter in a reporter assay. But not all CTCF sites genome-wide function as insulator elements, depending on cellular and genomic context. To dissect the influence of genomic context on enhancer blocker activity, we integrated reporter constructs with promoter-only, promoter and enhancer, and enhancer blocker configurations at hundreds of thousands of genomic sites using the Sleeping Beauty transposase. Deconvolution of reporter activity by genomic position reveals distinct expression patterns subject to genomic context, including a compartment of enhancer blocker reporter integrations with robust expression. The high density of integration sites permits quantitative delineation of characteristic genomic context sensitivity profiles and their decomposition into sensitivity to both local and distant DNase I hypersensitive sites. Furthermore, using a single-cell expression approach to test the effect of integrated reporters for differential expression of nearby endogenous genes reveals that CTCF insulator elements do not completely abrogate reporter effects on endogenous gene expression. Collectively, our results lend new insight into genomic regulatory compartmentalization and its influence on the determinants of promoter-enhancer specificity.

CREEPY: CRISPR-mediated editing of synthetic episomes in yeast.

Use of synthetic genomics to design and build 'big' DNA has revolutionized our ability to answer fundamental biological questions by employing a bottom-up approach. Saccharomyces cerevisiae, or budding yeast, has become the major platform to assemble large synthetic constructs thanks to its powerful homologous recombination machinery and the availability of well-established molecular biology techniques. However, introducing designer variations to episomal assemblies with high efficiency and fidelity remains challenging. Here we describe CRISPR Engineering of EPisomes in Yeast, or CREEPY, a method for rapid engineering of large synthetic episomal DNA constructs. We demonstrate that CRISPR editing of circular episomes presents unique challenges compared to modifying native yeast chromosomes. We optimize CREEPY for efficient and precise multiplex editing of >100 kb yeast episomes, providing an expanded toolkit for synthetic genomics.

A versatile platform for locus-scale genome rewriting and verification.

Routine rewriting of loci associated with human traits and diseases would facilitate their functional analysis. However, existing DNA integration approaches are limited in terms of scalability and portability across genomic loci and cellular contexts. We describe Big-IN, a versatile platform for targeted integration of large DNAs into mammalian cells. CRISPR/Cas9-mediated targeting of a landing pad enables subsequent recombinase-mediated delivery of variant payloads and efficient positive/negative selection for correct clones in mammalian stem cells. We demonstrate integration of constructs up to 143 kb, and an approach for one-step scarless delivery. We developed a staged pipeline combining PCR genotyping and targeted capture sequencing for economical and comprehensive verification of engineered stem cells. Our approach should enable combinatorial interrogation of genomic functional elements and systematic locus-scale analysis of genome function.

Inactivation of myocardin and p16 during malignant transformation contributes to a differentiation defect.

Myocardin is known as an important transcriptional regulator in smooth and cardiac muscle development. Here we found that myocardin is frequently repressed during human malignant transformation, contributing to a differentiation defect. We demonstrate that myocardin is a transcriptional target of TGFbeta required for TGFbeta-mediated differentiation of human fibroblasts. Serum deprivation, intact contact inhibition response, and the p16ink4a/Rb pathway contribute to myocardin induction and differentiation. Restoration of myocardin expression in sarcoma cells results in differentiation and inhibition of malignant growth, whereas inactivation of myocardin in normal fibroblasts increases their proliferative potential. Myocardin expression is reduced in multiple types of human tumors. Collectively, our results demonstrate that myocardin is an important suppressive modifier of the malignant transformation process.

Various p53 mutant proteins differently regulate the Ras circuit to induce a cancer-related gene signature.

Concomitant expression of mutant p53 and oncogenic Ras, leading to cellular transformation, is well documented. However, the mechanisms by which the various mutant p53 categories cooperate with Ras remain largely obscure. From this study we suggest that different mutant p53 categories cooperate with H-Ras in different ways to induce a unique expression pattern of a cancer-related gene signature (CGS). The DNA-contact p53 mutants (p53(R248Q) and p53(R273H)) exhibited the highest level of CGS expression by cooperating with NFκB. Furthermore, the Zn(+2) region conformational p53 mutants (p53(R175H) and p53(H179R)) induced the CGS by elevating H-Ras activity. This elevation in H-Ras activity stemmed from a perturbed function of the p53 transcription target gene, BTG2. By contrast, the L3 loop region conformational mutant (p53(G245S)) did not affect CGS expression. Our findings were further corroborated in human tumor-derived cell lines expressing Ras and the aforementioned mutated p53 proteins. These data might assist in future tailor-made therapy targeting the mutant p53-Ras axis in cancer.

Genomic context sensitizes regulatory elements to genetic disruption.

Enhancer function is frequently investigated piecemeal using truncated reporter assays or single deletion analysis. Thus it remains unclear to what extent enhancer function at native loci relies on surrounding genomic context. Using the Big-IN technology for targeted integration of large DNAs, we analyzed the regulatory architecture of the murine Igf2/H19 locus, a paradigmatic model of enhancer selectivity. We assembled payloads containing a 157-kb functional Igf2/H19 locus and engineered mutations to genetically direct CTCF occupancy at the imprinting control region (ICR) that switches the target gene of the H19 enhancer cluster. Contrasting activity of payloads delivered at the endogenous Igf2/H19 locus or ectopically at Hprt revealed that the Igf2/H19 locus includes additional, previously unknown long-range regulatory elements. Exchanging components of the Igf2/H19 locus with the well-studied Sox2 locus showed that the H19 enhancer cluster functioned poorly out of context, and required its native surroundings to activate Sox2 expression. Conversely, the Sox2 locus control region (LCR) could activate both Igf2 and H19 outside its native context, but its activity was only partially modulated by CTCF occupancy at the ICR. Analysis of regulatory DNA actuation across different cell types revealed that, while the H19 enhancers are tightly coordinated within their native locus, the Sox2 LCR acts more independently. We show that these enhancer clusters typify broader classes of loci genome-wide. Our results show that unexpected dependencies may influence even the most studied functional elements, and our synthetic regulatory genomics approach permits large-scale manipulation of complete loci to investigate the relationship between locus architecture and function.

Longitudinal scRNA-seq analysis in mouse and human informs optimization of rapid mouse astrocyte differentiation protocols.

Macroglia (astrocytes and oligodendrocytes) are required for normal development and function of the central nervous system, yet many questions remain about their emergence during the development of the brain and spinal cord. Here we used single-cell/single-nucleus RNA sequencing (scRNA-seq/snRNA-seq) to analyze over 298,000 cells and nuclei during macroglia differentiation from mouse embryonic and human-induced pluripotent stem cells. We computationally identify candidate genes involved in the fate specification of glia in both species and report heterogeneous expression of astrocyte surface markers across differentiating cells. We then used our transcriptomic data to optimize a previous mouse astrocyte differentiation protocol, decreasing the overall protocol length and complexity. Finally, we used multi-omic, dual single-nuclei (sn)RNA-seq/snATAC-seq analysis to uncover potential genomic regulatory sites mediating glial differentiation. These datasets will enable future optimization of glial differentiation protocols and provide insight into human glial differentiation.

Regulatory architecture of housekeeping genes is driven by promoter assemblies.

Genes that are key to cell identity are generally regulated by cell-type-specific enhancer elements bound by transcription factors, some of which facilitate looping to distant gene promoters. In contrast, genes that encode housekeeping functions, whose regulation is essential for normal cell metabolism and growth, generally lack interactions with distal enhancers. We find that Ronin (Thap11) assembles multiple promoters of housekeeping and metabolic genes to regulate gene expression. This behavior is analogous to how enhancers are brought together with promoters to regulate cell identity genes. Thus, Ronin-dependent promoter assemblies provide a mechanism to explain why housekeeping genes can forgo distal enhancer elements and why Ronin is important for cellular metabolism and growth control. We propose that clustering of regulatory elements is a mechanism common to cell identity and housekeeping genes but is accomplished by different factors binding distinct control elements to establish enhancer-promoter or promoter-promoter interactions, respectively.

A conditional counterselectable Piga knockout in mouse embryonic stem cells for advanced genome writing applications.

Overwriting counterselectable markers is an efficient strategy for removing wild-type DNA or replacing it with payload DNA of interest. Currently, one bottleneck of efficient genome engineering in mammals is the shortage of counterselectable (negative selection) markers that work robustly without affecting organismal developmental potential. Here, we report a conditional Piga knockout strategy that enables efficient proaerolysin-based counterselection in mouse embryonic stem cells. The conditional Piga knockout cells show similar proaerolysin resistance as full (non-conditional) Piga deletion cells, which enables the use of a PIGA transgene as a counterselectable marker for genome engineering purposes. Native Piga function is readily restored in conditional Piga knockout cells to facilitate subsequent mouse development. We also demonstrate the generality of our strategy by engineering a conditional knockout of endogenous Hprt. Taken together, our work provides a new tool for advanced mouse genome writing and mouse model establishment.

Synthetic regulatory reconstitution reveals principles of mammalian Hox cluster regulation.

Precise Hox gene expression is crucial for embryonic patterning. Intra-Hox transcription factor binding and distal enhancer elements have emerged as the major regulatory modules controlling Hox gene expression. However, quantifying their relative contributions has remained elusive. Here, we introduce "synthetic regulatory reconstitution," a conceptual framework for studying gene regulation, and apply it to the HoxA cluster. We synthesized and delivered variant rat HoxA clusters (130 to 170 kilobases) to an ectopic location in the mouse genome. We found that a minimal HoxA cluster recapitulated correct patterns of chromatin remodeling and transcription in response to patterning signals, whereas the addition of distal enhancers was needed for full transcriptional output. Synthetic regulatory reconstitution could provide a generalizable strategy for deciphering the regulatory logic of gene expression in complex genomes.

Transcriptional activity of ATF3 in the stromal compartment of tumors promotes cancer progression.

Compelling evidences have rendered the tumor microenvironment a crucial determinant in cancer outcome. Activating transcription factor 3 (ATF3), a stress response transcription factor, is known to have a dichotomous role in tumor cells, acting either as a tumor suppressor or an oncogene in a context-dependent manner. However, its expression and possible role in the tumor microenvironment are hitherto unknown. Here we show that ATF3 is upregulated in the stromal compartment of several types of cancer. Accordingly, Cancer-associated fibroblasts (CAFs) ectopically expressing ATF3 proliferated faster as indicated by increased colony-forming capacity and promoted the growth of adjacent tumor cells when co-injected into nude mice. Utilizing a genome-wide profiling approach, we unraveled a robust gene expression program induced by ATF3 in CAFs. Focusing on a specific subset of genes, we found that the ability of stromal ATF3 to promote cancer progression is mediated by transcriptional repression of CLDN1 and induction of CXCL12 and RGS4. In addition, regulation of LIF, CLDN1, SERPINE2, HSD17B2, ITGA7 and PODXL by ATF3 mediated the increased proliferation capacity of CAFs. In sum, our findings implicate ATF3 as a novel stromal tumor promoter and suggest that targeting ATF3 pathway might be beneficial for anticancer therapy.

De novo assembly and delivery to mouse cells of a 101 kb functional human gene.

Design and large-scale synthesis of DNA has been applied to the functional study of viral and microbial genomes. New and expanded technology development is required to unlock the transformative potential of such bottom-up approaches to the study of larger mammalian genomes. Two major challenges include assembling and delivering long DNA sequences. Here, we describe a workflow for de novo DNA assembly and delivery that enables functional evaluation of mammalian genes on the length scale of 100 kilobase pairs (kb). The DNA assembly step is supported by an integrated robotic workcell. We demonstrate assembly of the 101 kb human HPRT1 gene in yeast from 3 kb building blocks, precision delivery of the resulting construct to mouse embryonic stem cells, and subsequent expression of the human protein from its full-length human gene in mouse cells. This workflow provides a framework for mammalian genome writing. We envision utility in producing designer variants of human genes linked to disease and their delivery and functional analysis in cell culture or animal models.

p53 is balancing development, differentiation and de-differentiation to assure cancer prevention.

Many of the roles played by the tumor suppressor p53 in restraining cancer initiation and progression are well established. These include the ability of p53 to induce cell-cycle arrest, DNA repair, senescence and apoptosis. In addition, during the 30 years of p53 research, numerous studies have implicated p53 in the regulation of differentiation and developmental pathways. Here, we summarize the data on these relatively less-characterized functions of p53, including its involvement in embryogenesis and various differentiation programs, as well as its function in restraining de-differentiation of mature somatic cells. Besides the well-known functions of p53 as a cell-cycle regulator and a mediator of apoptosis, both coincide with differentiation processes, p53 was shown to exert its effects on various differentiation programs via direct regulation of specific key factors controlling these programs. The complex regulation by p53, which acts to suppress or to induce differentiation, is mainly the result of the specific cell type and fate. We argue that regulation of differentiation is pivotal for the tumor-suppressive activity of p53, which act to maintain the proper cellular state, preventing improper maturation or reprogramming. This conclusion is further supporting the notion that aberrant differentiation is associated with malignant transformation.

A critical role of PRDM14 in human primordial germ cell fate revealed by inducible degrons.

PRDM14 is a crucial regulator of mouse primordial germ cells (mPGCs), epigenetic reprogramming and pluripotency, but its role in the evolutionarily divergent regulatory network of human PGCs (hPGCs) remains unclear. Besides, a previous knockdown study indicated that PRDM14 might be dispensable for human germ cell fate. Here, we decided to use inducible degrons for a more rapid and comprehensive PRDM14 depletion. We show that PRDM14 loss results in significantly reduced specification efficiency and an aberrant transcriptome of hPGC-like cells (hPGCLCs) obtained in vitro from human embryonic stem cells (hESCs). Chromatin immunoprecipitation and transcriptomic analyses suggest that PRDM14 cooperates with TFAP2C and BLIMP1 to upregulate germ cell and pluripotency genes, while repressing WNT signalling and somatic markers. Notably, PRDM14 targets are not conserved between mouse and human, emphasising the divergent molecular mechanisms of PGC specification. The effectiveness of degrons for acute protein depletion is widely applicable in various developmental contexts.

Memory of Divisional History Directs the Continuous Process of Primitive Hematopoietic Lineage Commitment.

Hematopoietic stem cells (HSCs) exist in a dormant state and progressively lose regenerative potency as they undergo successive divisions. Why this functional decline occurs and how this information is encoded is unclear. To better understand how this information is stored, we performed RNA sequencing on HSC populations differing only in their divisional history. Comparative analysis revealed that genes upregulated with divisions are enriched for lineage genes and regulated by cell-cycle-associated transcription factors, suggesting that proliferation itself drives lineage priming. Downregulated genes are, however, associated with an HSC signature and targeted by the Polycomb Repressive Complex 2 (PRC2). The PRC2 catalytic subunits Ezh1 and Ezh2 promote and suppress the HSC state, respectively, and successive divisions cause a switch from Ezh1 to Ezh2 dominance. We propose that cell divisions drive lineage priming and Ezh2 accumulation, which represses HSC signature genes to consolidate information on divisional history into memory.