A versatile platform for locus-scale genome rewriting and verification

Identification of the Molecular Components of Enhancer-Mediated Gene Expression Variation in Multiple Tissues Regulating Blood Pressure

BACKGROUND

Inter-individual variation in blood pressure (BP) arises in part from sequence variants within enhancers modulating the expression of causal genes. We propose that these genes, active in tissues relevant to BP physiology, can be identified from tissue-level epigenomic data and genotypes of BP-phenotyped individuals.

METHODS

We used chromatin accessibility data from the heart, adrenal, kidney, and artery to identify cis-regulatory elements (CREs) in these tissues and estimate the impact of common human single-nucleotide variants within these CREs on gene expression, using machine learning methods. To identify causal genes, we performed a gene-wise association test. We conducted analyses in 2 separate large-scale cohorts: 77 822 individuals from the Genetic Epidemiology Research on Adult Health and Aging and 315 270 individuals from the UK Biobank.

RESULTS

We identified 309, 259, 331, and 367 genes (false discovery rate <0.05) for diastolic BP and 191, 184, 204, and 204 genes for systolic BP in the artery, kidney, heart, and adrenal, respectively, in Genetic Epidemiology Research on Adult Health and Aging; 50% to 70% of these genes were replicated in the UK Biobank, significantly higher than the 12% to 15% expected by chance (P<0.0001). These results enabled tissue expression prediction of these 988 to 2875 putative BP genes in individuals of both cohorts to construct an expression polygenic score. This score explained ≈27% of the reported single-nucleotide variant heritability, substantially higher than expected from prior studies.

CONCLUSIONS

Our work demonstrates the power of tissue-restricted comprehensive CRE analysis, followed by CRE-based expression prediction, for understanding BP regulation in relevant tissues and provides dual-modality supporting evidence, CRE and expression, for the causality genes.

Hepatic organoids move from adolescence to maturity

Since organoids were developed 15 years ago, they are now in their adolescence as a research tool. The ability to generate 'tissue in a dish' has created enormous opportunities for biomedical research. We examine the contributions that hepatic organoids have made to three areas of liver research: as a source of cells and tissue for basic research, for drug discovery and drug safety testing, and for understanding disease pathobiology. We discuss the features that enable hepatic organoids to provide useful models for human liver diseases and identify four types of advances that will enable them to become a mature (i.e., adult) research tool over the next 5 years. During this period, advances in single-cell RNA sequencing and CRISPR technologies coupled with improved hepatic organoid methodology, which enables them to have a wider range of cell types that are present in liver and to be grown in microwells, will generate discoveries that will dramatically advance our understanding of liver development and the pathogenesis of liver diseases. It will generate also new approaches for treating liver fibrosis, which remains a major public health problem with few treatment options.

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 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.

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.

Multiplex generation and single cell analysis of structural variants in a mammalian genome

The functional consequences of structural variants (SVs) in mammalian genomes are challenging to study. This is due to several factors, including: 1) their numerical paucity relative to other forms of standing genetic variation such as single nucleotide variants (SNVs) and short insertions or deletions (indels); 2) the fact that a single SV can involve and potentially impact the function of more than one gene and/or cis regulatory element; and 3) the relative immaturity of methods to generate and map SVs, either randomly or in targeted fashion, in in vitro or in vivo model systems. Towards addressing these challenges, we developed Genome-Shuffle-seq, a straightforward method that enables the multiplex generation and mapping of several major forms of SVs (deletions, inversions, translocations) throughout a mammalian genome. Genome-Shuffle-seq is based on the integration of "shuffle cassettes" to the genome, wherein each shuffle cassette contains components that facilitate its site-specific recombination (SSR) with other integrated shuffle cassettes (via Cre-loxP), its mapping to a specific genomic location (via T7-mediated in vitro transcription or IVT), and its identification in single-cell RNA-seq (scRNA-seq) data (via T7-mediated in situ transcription or IST). In this proof-of-concept, we apply Genome-Shuffle-seq to induce and map thousands of genomic SVs in mouse embryonic stem cells (mESCs) in a single experiment. Induced SVs are rapidly depleted from the cellular population over time, possibly due to Cre-mediated toxicity and/or negative selection on the rearrangements themselves. Leveraging T7 IST of barcodes whose positions are already mapped, we further demonstrate that we can efficiently genotype which SVs are present in association with each of many single cell transcriptomes in scRNA-seq data. Finally, preliminary evidence suggests our method may be a powerful means of generating extrachromosomal circular DNAs (ecDNAs). Looking forward, we anticipate that Genome-Shuffle-seq may be broadly useful for the systematic exploration of the functional consequences of SVs on gene expression, the chromatin landscape, and 3D nuclear architecture. We further anticipate potential uses for in vitro modeling of ecDNAs, as well as in paving the path to a minimal mammalian genome.

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.

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.

Cell-type-specific prediction of 3D chromatin organization enables high-throughput in silico genetic screening

Investigating how chromatin organization determines cell-type-specific gene expression remains challenging. Experimental methods for measuring three-dimensional chromatin organization, such as Hi-C, are costly and have technical limitations, restricting their broad application particularly in high-throughput genetic perturbations. We present C.Origami, a multimodal deep neural network that performs de novo prediction of cell-type-specific chromatin organization using DNA sequence and two cell-type-specific genomic features-CTCF binding and chromatin accessibility. C.Origami enables in silico experiments to examine the impact of genetic changes on chromatin interactions. We further developed an in silico genetic screening approach to assess how individual DNA elements may contribute to chromatin organization and to identify putative cell-type-specific trans-acting regulators that collectively determine chromatin architecture. Applying this approach to leukemia cells and normal T cells, we demonstrate that cell-type-specific in silico genetic screening, enabled by C.Origami, can be used to systematically discover novel chromatin regulation circuits in both normal and disease-related biological systems.

Functional roles of reactive astrocytes in neuroinflammation and neurodegeneration

Despite advances in uncovering the mechanisms that underlie neuroinflammation and neurodegenerative disease, therapies that prevent neuronal loss remain elusive. Targeting of disease-defining markers in conditions such as Alzheimer disease (amyloid-β and tau) or Parkinson disease (α-synuclein) has been met with limited success, suggesting that these proteins do not act in isolation but form part of a pathological network. This network could involve phenotypic alteration of multiple cell types in the CNS, including astrocytes, which have a major neurosupportive, homeostatic role in the healthy CNS but adopt reactive states under acute or chronic adverse conditions. Transcriptomic studies in human patients and disease models have revealed the co-existence of many putative reactive sub-states of astrocytes. Inter-disease and even intra-disease heterogeneity of reactive astrocytic sub-states are well established, but the extent to which specific sub-states are shared across different diseases is unclear. In this Review, we highlight how single-cell and single-nuclei RNA sequencing and other 'omics' technologies can enable the functional characterization of defined reactive astrocyte states in various pathological scenarios. We provide an integrated perspective, advocating cross-modal validation of key findings to define functionally important sub-states of astrocytes and their triggers as tractable therapeutic targets with cross-disease relevance.

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.

CRaTER enrichment for on-target gene editing enables generation of variant libraries in hiPSCs

Standard transgenic cell line generation requires screening 100-1000s of colonies to isolate correctly edited cells. We describe CRISPRa On-Target Editing Retrieval (CRaTER) which enriches for cells with on-target knock-in of a cDNA-fluorescent reporter transgene by transient activation of the targeted locus followed by flow sorting to recover edited cells. We show CRaTER recovers rare cells with heterozygous, biallelic-editing of the transcriptionally-inactive MYH7 locus in human induced pluripotent stem cells (hiPSCs), enriching on average 25-fold compared to standard antibiotic selection. We leveraged CRaTER to enrich for heterozygous knock-in of a library of variants in MYH7, a gene in which missense mutations cause cardiomyopathies, and recovered hiPSCs with 113 different variants. We differentiated these hiPSCs to cardiomyocytes and show MHC-β fusion proteins can localize as expected. Additionally, single-cell contractility analyses revealed cardiomyocytes with a pathogenic, hypertrophic cardiomyopathy-associated MYH7 variant exhibit salient HCM physiology relative to isogenic controls. Thus, CRaTER substantially reduces screening required for isolation of gene-edited cells, enabling generation of functional transgenic cell lines at unprecedented scale.

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.

CRaTER enrichment for on-target gene-editing enables generation of variant libraries in hiPSCs

Standard transgenic cell line generation requires screening 100-1000s of colonies to isolate correctly edited cells. We describe CR ISPR a On- T arget E diting R etrieval (CRaTER) which enriches for cells with on-target knock-in of a cDNA-fluorescent reporter transgene by transient activation of the targeted locus followed by flow sorting to recover edited cells. We show CRaTER recovers rare cells with heterozygous, biallelic-editing of the transcriptionally-inactive MYH7 locus in human induced pluripotent stem cells (hiPSCs), enriching on average 25-fold compared to standard antibiotic selection. We leveraged CRaTER to enrich for heterozygous knock-in of a library of single nucleotide variants (SNVs) in MYH7 , a gene in which missense mutations cause cardiomyopathies, and recovered hiPSCs with 113 different MYH7 SNVs. We differentiated these hiPSCs to cardiomyocytes and show MYH7 fusion proteins can localize as expected. Thus, CRaTER substantially reduces screening required for isolation of gene-edited cells, enabling generation of transgenic cell lines at unprecedented scale.

Systems genomics in age-related macular degeneration

Genomic studies in age-related macular degeneration (AMD) have identified genetic variants that account for the majority of AMD risk. An important next step is to understand the functional consequences and downstream effects of the identified AMD-associated genetic variants. Instrumental for this next step are 'omics' technologies, which enable high-throughput characterization and quantification of biological molecules, and subsequent integration of genomics with these omics datasets, a field referred to as systems genomics. Single cell sequencing studies of the retina and choroid demonstrated that the majority of candidate AMD genes identified through genomic studies are expressed in non-neuronal cells, such as the retinal pigment epithelium (RPE), glia, myeloid and choroidal cells, highlighting that many different retinal and choroidal cell types contribute to the pathogenesis of AMD. Expression quantitative trait locus (eQTL) studies in retinal tissue have identified putative causal genes by demonstrating a genetic overlap between gene regulation and AMD risk. Linking genetic data to complement measurements in the systemic circulation has aided in understanding the effect of AMD-associated genetic variants in the complement system, and supports that protein QTL (pQTL) studies in plasma or serum samples may aid in understanding the effect of genetic variants and pinpointing causal genes in AMD. A recent epigenomic study fine-mapped AMD causal variants by determing regulatory regions in RPE cells differentiated from induced pluripotent stem cells (iPSC-RPE). Another approach that is being employed to pinpoint causal AMD genes is to produce synthetic DNA assemblons representing risk and protective haplotypes, which are then delivered to cellular or animal model systems. Pinpointing causal genes and understanding disease mechanisms is crucial for the next step towards clinical translation. Clinical trials targeting proteins encoded by the AMD-associated genomic loci C3, CFB, CFI, CFH, and ARMS2/HTRA1 are currently ongoing, and a phase III clinical trial for C3 inhibition recently showed a modest reduction of lesion growth in geographic atrophy. The EYERISK consortium recently developed a genetic test for AMD that allows genotyping of common and rare variants in AMD-associated genes. Polygenic risk scores (PRS) were applied to quantify AMD genetic risk, and may aid in predicting AMD progression. In conclusion, genomic studies represent a turning point in our exploration of AMD. The results of those studies now serve as a driving force for several clinical trials. Expanding to omics and systems genomics will further decipher function and causality from the associations that have been reported, and will enable the development of therapies that will lessen the burden of AMD.

STRAIGHT-IN enables high-throughput targeting of large DNA payloads in human pluripotent stem cells

Inserting large DNA payloads (>10 kb) into specific genomic sites of mammalian cells remains challenging. Applications ranging from synthetic biology to evaluating the pathogenicity of disease-associated variants for precision medicine initiatives would greatly benefit from tools that facilitate this process. Here, we merge the strengths of different classes of site-specific recombinases and combine these with CRISPR-Cas9-mediated homologous recombination to develop a strategy for stringent site-specific replacement of genomic fragments at least 50 kb in size in human induced pluripotent stem cells (hiPSCs). We demonstrate the versatility of STRAIGHT-IN (serine and tyrosine recombinase-assisted integration of genes for high-throughput investigation) by (1) inserting various combinations of fluorescent reporters into hiPSCs to assess the excitation-contraction coupling cascade in derivative cardiomyocytes and (2) simultaneously targeting multiple variants associated with inherited cardiac arrhythmic disorders into a pool of hiPSCs. STRAIGHT-IN offers a precise approach to generate genetically matched panels of hiPSC lines efficiently and cost effectively.

Natural and Experimental Rewiring of Gene Regulatory Regions

The successful development and ongoing functioning of complex organisms depend on the faithful execution of the genetic code. A critical step in this process is the correct spatial and temporal expression of genes. The highly orchestrated transcription of genes is controlled primarily by cis-regulatory elements: promoters, enhancers, and insulators. The medical importance of this key biological process can be seen by the frequency with which mutations and inherited variants that alter cis-regulatory elements lead to monogenic and complex diseases and cancer. Here, we provide an overview of the methods available to characterize and perturb gene regulatory circuits. We then highlight mechanisms through which regulatory rewiring contributes to disease, and conclude with a perspective on how our understanding of gene regulation can be used to improve human health.

In vivo processing of digital information molecularly with targeted specificity and robust reliability

DNA has attracted increasing interest as an appealing medium for information storage. However, target-specific rewriting of the digital data stored in intracellular DNA remains a grand challenge because the highly repetitive nature and uneven guanine-cytosine content render the encoded DNA sequences poorly compatible with endogenous ones. In this study, a dual-plasmid system based on gene editing tools was introduced into Escherichia coli to process information accurately. Digital data containing large repeat units in binary codes, such as text, codebook, or image, were involved in the realization of target-specific rewriting in vivo, yielding up to 94% rewriting reliability. An optical reporter was introduced as an advanced tool for presenting data processing at the molecular level. Rewritten information was stored stably and amplified over hundreds of generations. Our work demonstrates a digital-to-biological information processing approach for highly efficient data storage, amplification, and rewriting, thus robustly promoting the application of DNA-based information technology.

Programmable RNA-Guided Large DNA Transgenesis by CRISPR/Cas9 and Site-Specific Integrase Bxb1

The inability to insert large DNA constructs into the genome efficiently and precisely is a key challenge in genomic engineering. Random transgenesis, which is widely used, lacks precision, and comes with a slew of drawbacks. Lentiviral and adeno-associated viral methods are plagued by, respectively, DNA toxicity and a payload capacity of less than 5 kb. Homology-directed repair (HDR) techniques based on CRISPR-Cas9 can be effective, but only in the 1–5 kb range. In addition, long homology arms—DNA sequences that permit construct insertion—of lengths ranging from 0.5 to 5 kb are required by currently known HDR-based techniques. A potential new method that uses Cas9-guided transposases to insert DNA structures up to 10 kb in length works well in bacteria, but only in bacteria. Surmounting these roadblocks, a new toolkit has recently been developed that combines RNA-guided Cas9 and the site-specific integrase Bxb1 to integrate DNA constructs ranging in length from 5 to 43 kb into mouse zygotes with germline transmission and into human cells. This ground-breaking toolkit will give researchers a valuable resource for developing novel, urgently needed mouse and human induced pluripotent stem cell (hiPSC) models of cancer and other genetic diseases, as well as therapeutic gene integration and biopharmaceutical applications, such as the development of stable cell lines to produce therapeutic protein products.

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.

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.

Efficient targeted transgenesis of large donor DNA into multiple mouse genetic backgrounds using bacteriophage Bxb1 integrase

The development of mouse models of human disease and synthetic biology research by targeted transgenesis of large DNA constructs represent a significant genetic engineering hurdle. We developed an efficient, precise, single-copy integration of large transgenes directly into zygotes using multiple mouse genetic backgrounds. We used in vivo Bxb1 mediated recombinase-mediated cassette exchange (RMCE) with a transgene “landing pad” composed of dual heterologous Bxb1 attachment (att) sites in cis, within the Gt(ROSA)26Sor safe harbor locus. RMCE of donor was achieved by microinjection of vector DNA carrying cognate attachment sites flanking the donor transgene with Bxb1-integrase mRNA. This approach achieves perfect vector-free integration of donor constructs at efficiencies > 40% with up to ~ 43 kb transgenes. Coupled with a nanopore-based Cas9-targeted sequencing (nCATS), complete verification of precise insertion sequence was achieved. As a proof-of-concept we describe the development of C57BL/6J and NSG Krt18-ACE2 models for SARS-CoV2 research with verified heterozygous N1 animals within ~ 4 months. Additionally, we created a series of mice with diverse backgrounds carrying a single att site including FVB/NJ, PWK/PhJ, NOD/ShiLtJ, CAST/EiJ and DBA/2J allowing for rapid transgene insertion. Combined, this system enables predictable, rapid development with simplified characterization of precisely targeted transgenic animals across multiple genetic backgrounds.

Darwinian genomics and diversity in the tree of life

Genomics encompasses the entire tree of life, both extinct and extant, and the evolutionary processes that shape this diversity. To date, genomic research has focused on humans, a small number of agricultural species, and established laboratory models. Fewer than 18,000 of ∼2,000,000 eukaryotic species (<1%) have a representative genome sequence in GenBank, and only a fraction of these have ancillary information on genome structure, genetic variation, gene expression, epigenetic modifications, and population diversity. This imbalance reflects a perception that human studies are paramount in disease research. Yet understanding how genomes work, and how genetic variation shapes phenotypes, requires a broad view that embraces the vast diversity of life. We have the technology to collect massive and exquisitely detailed datasets about the world, but expertise is siloed into distinct fields. A new approach, integrating comparative genomics with cell and evolutionary biology, ecology, archaeology, anthropology, and conservation biology, is essential for understanding and protecting ourselves and our world. Here, we describe potential for scientific discovery when comparative genomics works in close collaboration with a broad range of fields as well as the technical, scientific, and social constraints that must be addressed.

Tissue context determines the penetrance of regulatory DNA variation

Functional assessment of disease-associated sequence variation at non-coding regulatory elements is complicated by their high degree of context sensitivity to both the local chromatin and nuclear environments. Allelic profiling of DNA accessibility across individuals has shown that only a select minority of sequence variation affects transcription factor (TF) occupancy, yet low sequence diversity in human populations means that no experimental assessment is available for the majority of disease-associated variants. Here we describe high-resolution in vivo maps of allelic DNA accessibility in liver, kidney, lung and B cells from 5 increasingly diverged strains of F1 hybrid mice. The high density of heterozygous sites in these hybrids enables precise quantification of effect size and cell-type specificity for hundreds of thousands of variants throughout the mouse genome. We show that chromatin-altering variants delineate characteristic sensitivity profiles for hundreds of TF motifs. We develop a compendium of TF-specific sensitivity profiles accounting for genomic context effects. Finally, we link maps of allelic accessibility to allelic transcript levels in the same samples. This work provides a foundation for quantitative prediction of cell-type specific effects of non-coding variation on TF activity, which will facilitate both fine-mapping and systems-level analyses of common disease-associated variation in human genomes.

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.