Large-scale analysis of the regulatory architecture of the mouse genome with a transposon-associated sensor

Analysis of the transcriptional activity of model piggyBac transgenes stably integrated into different loci of the genome of CHO cells in the absence of selection pressure

CHO cells are most commonly used for the synthesis of recombinant proteins in biopharmaceutical production. When stable producer cell lines are obtained, the locus of transgene integration into the genome has a great influence on the level of its expression. Therefore, the identification of genomic loci ensuring a high level of protein production is very important. Here, we used the TRIP assay to study the influence of the local chromatin environment on the activity of transgenes in CHO cells. For this purpose, reporter constructs encoding eGFP under the control of four promoters were stably integrated into the genome of CHO cells using the piggyBac transposon. Each individual transgene contained a unique tag, a DNA barcode, and the resulting polyclonal cell population was cultured for almost a month without any selection. Next, using the high-throughput sequencing, genomic localizations of barcodes, as well as their abundances in the population and transcriptional activities were identified. In total, ~640 transgenes more or less evenly distributed across all chromosomes of CHO cells were characterized. More than half of the transgenes were completely silent. The most active transgenes were identified to be inserted in gene promoters and 5' UTRs. Transgenes carrying Chinese hamster full-length promoter of the EF-1α gene showed the highest activity. Transgenes with a truncated version of the same promoter and with the mouse PGK gene promoter were on average 10 and 19 times less active, respectively. In total, combinations of genomic loci of CHO cells and transgene promoters that together provide different levels of transcriptional activity of the model reporter construct were described.

Advances in transposable elements: from mechanisms to applications in mammalian genomics

It has been 70 years since Barbara McClintock discovered transposable elements (TE), and the mechanistic studies and functional applications of transposable elements have been at the forefront of life science research. As an essential part of the genome, TEs have been discovered in most species of prokaryotes and eukaryotes, and the relative proportion of the total genetic sequence they comprise gradually increases with the expansion of the genome. In humans, TEs account for about 40% of the genome and are deeply involved in gene regulation, chromosome structure maintenance, inflammatory response, and the etiology of genetic and non-genetic diseases. In-depth functional studies of TEs in mammalian cells and the human body have led to a greater understanding of these fundamental biological processes. At the same time, as a potent mutagen and efficient genome editing tool, TEs have been transformed into biological tools critical for developing new techniques. By controlling the random insertion of TEs into the genome to change the phenotype in cells and model organisms, critical proteins of many diseases have been systematically identified. Exploiting the TE's highly efficient in vitro insertion activity has driven the development of cutting-edge sequencing technologies. Recently, a new technology combining CRISPR with TEs was reported, which provides a novel targeted insertion system to both academia and industry. We suggest that interrogating biological processes that generally depend on the actions of TEs with TEs-derived genetic tools is a very efficient strategy. For example, excessive activation of TEs is an essential factor in the occurrence of cancer in humans. As potent mutagens, TEs have also been used to unravel the key regulatory elements and mechanisms of carcinogenesis. Through this review, we aim to effectively combine the traditional views of TEs with recent research progress, systematically link the mechanistic discoveries of TEs with the technological developments of TE-based tools, and provide a comprehensive approach and understanding for researchers in different fields.

Disruption of regulatory domains and novel transcripts as disease-causing mechanisms

Deletions, duplications, insertions, inversions, and translocations, collectively called structural variations (SVs), affect more base pairs of the genome than any other sequence variant. The recent technological advancements in genome sequencing have enabled the discovery of tens of thousands of SVs per human genome. These SVs primarily affect non-coding DNA sequences, but the difficulties in interpreting their impact limit our understanding of human disease etiology. The functional annotation of non-coding DNA sequences and methodologies to characterize their three-dimensional (3D) organization in the nucleus have greatly expanded our understanding of the basic mechanisms underlying gene regulation, thereby improving the interpretation of SVs for their pathogenic impact. Here, we discuss the various mechanisms by which SVs can result in altered gene regulation and how these mechanisms can result in rare genetic disorders. Beyond changing gene expression, SVs can produce novel gene-intergenic fusion transcripts at the SV breakpoints.

Integrative approaches to study enhancer-promoter communication

The spatiotemporal control of gene expression in complex multicellular organisms relies on noncoding regulatory sequences such as enhancers, which activate transcription of target genes often over large genomic distances. Despite the advances in the identification and characterization of enhancers, the principles and mechanisms by which enhancers select and control their target genes remain largely unknown. Here, we review recent interdisciplinary and quantitative approaches based on emerging techniques that aim to address open questions in the field, notably how regulatory information is encoded in the DNA sequence, how this information is transferred from enhancers to promoters, and how these processes are regulated in time.

CRISMERE Chromosome Engineering in Mouse and Rat

CRISPR/Cas9 technology is a versatile tool for engineering biology that has dramatically transformed our ability to manipulate genomes. In this protocol, we use its capacity to generate two double-strand breaks simultaneously, at precise positions in the genome, to generate mouse or rat lines with deletion, inversion, and duplication of a specific genomic segment. The technic is called CRISMERE for CRISpr-MEdiated REarrangement. This protocol describes the different steps to generate and validate the different chromosomal rearrangements that can be obtained with the technology. These new genetic configurations can be useful to model rare diseases with copy number variation, understand the genomic organization, or provide genetic tools (like balancer chromosome) to keep lethal mutations.

Genetic and Molecular Quality Control of Genetically Engineered Mice

Genetically engineered mice are used as avatars to understand mammalian gene function and develop therapies for human disease. During genetic modification, unintended changes can occur, and these changes may result in misassigned gene-phenotype relationships leading to incorrect or incomplete experimental interpretations. The types of unintended changes that may occur depend on the allele type being made and the genetic engineering approach used. Here we broadly categorize allele types as deletions, insertions, base changes, and transgenes derived from engineered embryonic stem (ES) cells or edited mouse embryos. However, the methods we describe can be adapted to other allele types and engineering strategies. We describe the sources and consequ ences of common unintended changes and best practices for detecting both intended and unintended changes by screening and genetic and molecular quality control (QC) of chimeras, founders, and their progeny. Employing these practices, along with careful allele design and good colony management, will increase the chance that investigations using genetically engineered mice will produce high-quality reproducible results, to enable a robust understanding of gene function, human disease etiology, and therapeutic development.

Enhancer-gene specificity in development and disease

Enhancers control the establishment of spatiotemporal gene expression patterns throughout development. Over the past decade, the development of new technologies has improved our capacity to link enhancers with their target genes based on their colocalization within the same topological domains. However, the mechanisms that regulate how enhancers specifically activate some genes but not others within a given domain remain unclear. In this Review, we discuss recent insights into the factors controlling enhancer specificity, including the genetic composition of enhancers and promoters, the linear and 3D distance between enhancers and their target genes, and cell-type specific chromatin landscapes. We also discuss how elucidating the molecular principles of enhancer specificity might help us to better understand and predict the pathological consequences of human genetic, epigenetic and structural variants.

Strand asymmetry influences mismatch resolution during a single-strand annealing

Background

Biases of DNA repair can shape the nucleotide landscape of genomes at evolutionary timescales. The molecular mechanisms of those biases are still poorly understood because it is difficult to isolate the contributions of DNA repair from those of DNA damage.

Results

Here, we develop a genome-wide assay whereby the same DNA lesion is repaired in different genomic contexts. We insert thousands of barcoded transposons carrying a reporter of DNA mismatch repair in the genome of mouse embryonic stem cells. Upon inducing a double-strand break between tandem repeats, a mismatch is generated if the break is repaired through single-strand annealing. The resolution of the mismatch showed a 60–80% bias in favor of the strand with the longest 3′ flap. The location of the lesion in the genome and the type of mismatch had little influence on the bias. Instead, we observe a complete reversal of the bias when the longest 3′ flap is moved to the opposite strand by changing the position of the double-strand break in the reporter.

Conclusions

These results suggest that the processing of the double-strand break has a major influence on the repair of mismatches during a single-strand annealing.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02665-3.

Cell-specific alterations in Pitx1 regulatory landscape activation caused by the loss of a single enhancer

Developmental genes are frequently controlled by multiple enhancers sharing similar specificities. As a result, deletions of such regulatory elements have often failed to reveal their full function. Here, we use the Pitx1 testbed locus to characterize in detail the regulatory and cellular identity alterations following the deletion of one of its enhancers (Pen). By combining single cell transcriptomics and an in-embryo cell tracing approach, we observe an increased fraction of Pitx1 non/low-expressing cells and a decreased fraction of Pitx1 high-expressing cells. We find that the over-representation of Pitx1 non/low-expressing cells originates from a failure of the Pitx1 locus to coordinate enhancer activities and 3D chromatin changes. This locus mis-activation induces a localized heterochrony and a concurrent loss of irregular connective tissue, eventually leading to a clubfoot phenotype. This data suggests that, in some cases, redundant enhancers may be used to locally enforce a robust activation of their host regulatory landscapes.

The Stochastic Genome and Its Role in Gene Expression

Mammalian genomes have distinct levels of spatial organization and structure that have been hypothesized to play important roles in transcription regulation. Although much has been learned about these architectural features with ensemble techniques, single-cell studies are showing a new universal trend: Genomes are stochastic and dynamic at every level of organization. Stochastic gene expression, on the other hand, has been studied for years. In this review, we probe whether there is a causative link between the two phenomena. We specifically discuss the functionality of chromatin state, topologically associating domains (TADs), and enhancer biology in light of their stochastic nature and their specific roles in stochastic gene expression. We highlight persistent fundamental questions in this area of research.

Mesomelic dysplasias associated with the HOXD locus are caused by regulatory reallocations

Human families with chromosomal rearrangements at 2q31, where the human HOXD locus maps, display mesomelic dysplasia, a severe shortening and bending of the limb. In mice, the dominant Ulnaless inversion of the HoxD cluster produces a similar phenotype suggesting the same origin for these malformations in humans and mice. Here we engineer 1 Mb inversion including the HoxD gene cluster, which positioned Hoxd13 close to proximal limb enhancers. Using this model, we show that these enhancers contact and activate Hoxd13 in proximal cells, inducing the formation of mesomelic dysplasia. We show that a secondary Hoxd13 null mutation in-cis with the inversion completely rescues the alterations, demonstrating that ectopic HOXD13 is directly responsible for this bone anomaly. Single-cell expression analysis and evaluation of HOXD13 binding sites suggests that the phenotype arises primarily by acting through genes normally controlled by HOXD13 in distal limb cells. Altogether, these results provide a conceptual and mechanistic framework to understand and unify the molecular origins of human mesomelic dysplasia associated with 2q31.

Murine allele and transgene symbols: ensuring unique, concise, and informative nomenclature

In addition to naturally occurring sequence variation and spontaneous mutations, a wide array of technologies exist for modifying the mouse genome. Standardized nomenclature, including allele, transgene, and other mutation nomenclature, as well as persistent unique identifiers (PUID) are critical for effective scientific communication, comparison of results, and integration of data into knowledgebases such as Mouse Genome Informatics (MGI), Alliance for Genome Resources, and International Mouse Strain Resource (IMSR). As well as being the authoritative source for mouse gene, allele, and strain nomenclature, MGI integrates published and unpublished genomic, phenotypic, and expression data while linking to other online resources for a complete view of the mouse as a valuable model organism. The International Committee on Standardized Genetic Nomenclature for Mice has developed allele nomenclature rules and guidelines that take into account the number of genes impacted, the method of allele generation, and the nature of the sequence alteration. To capture details that cannot be included in allele symbols, MGI has further developed allele to gene relationships using sequence ontology (SO) definitions for mutations that provide links between alleles and the genes affected. MGI is also using (HGVS) variant nomenclature for variants associated with alleles that will enhance searching for mutations and will improve cross-species comparison. With the ability to assign unique and informative symbols as well as to link alleles with more than one gene, allele and transgene nomenclature rules and guidelines provide an unambiguous way to represent alterations in the mouse genome and facilitate data integration among multiple resources such the Alliance of Genome Resources and International Mouse Strain Resource.

Genetic dissection identifies Necdin as a driver gene in a mouse model of paternal 15q duplications

Maternally inherited duplication of chromosome 15q11-q13 (Dup15q) is a pathogenic copy number variation (CNV) associated with autism spectrum disorder (ASD). Recently, paternally derived duplication has also been shown to contribute to the development of ASD. The molecular mechanism underlying paternal Dup15q remains unclear. Here, we conduct genetic and overexpression-based screening and identify Necdin (Ndn) as a driver gene for paternal Dup15q resulting in the development of ASD-like phenotypes in mice. An excess amount of Ndn results in enhanced spine formation and density as well as hyperexcitability of cortical pyramidal neurons. We generate 15q dupΔNdn mice with a normalized copy number of Ndn by excising its one copy from Dup15q mice using a CRISPR-Cas9 system. 15q dupΔNdn mice do not show ASD-like phenotypes and show dendritic spine dynamics and cortical excitatory-inhibitory balance similar to wild type animals. Our study provides an insight into the role of Ndn in paternal 15q duplication and a mouse model of paternal Dup15q syndrome.

Engineering three-dimensional genome folding

Animal genomes are partitioned and folded at various scales that contribute distinctly to nuclear processes. While structural features have been disrupted either globally or at select loci in loss-of-function studies, gain-of-function studies that probe the role of genome architecture have lagged behind. Here we examine recent advances in experimentally creating chromatin loops, contact domains, boundaries and compartments. Furthermore, we explore parallels between this emerging theme and natural evolution of mammalian genomes with increasing architectural complexity. Finally, we provide a perspective on how insights arising from recent gain-of-function studies may inform future endeavors toward engineering the three-dimensional genome.

A native, highly active Tc1/mariner transposon from zebrafish (ZB) offers an efficient genetic manipulation tool for vertebrates

New genetic tools and strategies are currently under development to facilitate functional genomics analyses. Here, we describe an active member of the Tc1/mariner transposon superfamily, named ZB, which invaded the zebrafish genome very recently. ZB exhibits high activity in vertebrate cells, in the range of those of the widely used transposons piggyBac (PB), Sleeping Beauty (SB) and Tol2. ZB has a similar structural organization and target site sequence preference to SB, but a different integration profile with respect to genome-wide preference among mammalian functional annotation features. Namely, ZB displays a preference for integration into transcriptional regulatory regions of genes. Accordingly, we demonstrate the utility of ZB for enhancer trapping in zebrafish embryos and in the mouse germline. These results indicate that ZB may be a powerful tool for genetic manipulation in vertebrate model species.

Jumping Ahead with Sleeping Beauty: Mechanistic Insights into Cut-and-Paste Transposition

Sleeping Beauty (SB) is a transposon system that has been widely used as a genetic engineering tool. Central to the development of any transposon as a research tool is the ability to integrate a foreign piece of DNA into the cellular genome. Driven by the need for efficient transposon-based gene vector systems, extensive studies have largely elucidated the molecular actors and actions taking place during SB transposition. Close transposon relatives and other recombination enzymes, including retroviral integrases, have served as useful models to infer functional information relevant to SB. Recently obtained structural data on the SB transposase enable a direct insight into the workings of this enzyme. These efforts cumulatively allowed the development of novel variants of SB that offer advanced possibilities for genetic engineering due to their hyperactivity, integration deficiency, or targeting capacity. However, many aspects of the process of transposition remain poorly understood and require further investigation. We anticipate that continued investigations into the structure-function relationships of SB transposition will enable the development of new generations of transposition-based vector systems, thereby facilitating the use of SB in preclinical studies and clinical trials.

Gene regulation gravitates toward either addition or multiplication when combining the effects of two signals

Two different cell signals often affect transcription of the same gene. In such cases, it is natural to ask how the combined transcriptional response compares to the individual responses. The most commonly used mechanistic models predict additive or multiplicative combined responses, but a systematic genome-wide evaluation of these predictions is not available. Here, we analyzed the transcriptional response of human MCF-7 cells to retinoic acid and TGF-β, applied individually and in combination. The combined transcriptional responses of induced genes exhibited a range of behaviors, but clearly favored both additive and multiplicative outcomes. We performed paired chromatin accessibility measurements and found that increases in accessibility were largely additive. There was some association between super-additivity of accessibility and multiplicative or super-multiplicative combined transcriptional responses, while sub-additivity of accessibility associated with additive transcriptional responses. Our findings suggest that mechanistic models of combined transcriptional regulation must be able to reproduce a range of behaviors.

Identifying chromatin features that regulate gene expression distribution

Gene expression variability, differences in the number of mRNA per cell across a population of cells, is ubiquitous across diverse organisms with broad impacts on cellular phenotypes. The role of chromatin in regulating average gene expression has been extensively studied. However, what aspects of the chromatin contribute to gene expression variability is still underexplored. Here we addressed this problem by leveraging chromatin diversity and using a systematic investigation of randomly integrated expression reporters to identify what aspects of chromatin microenvironment contribute to gene expression variability. Using DNA barcoding and split-pool decoding, we created a large library of isogenic reporter clones and identified reporter integration sites in a massive and parallel manner. By mapping our measurements of reporter expression at different genomic loci with multiple epigenetic profiles including the enrichment of transcription factors and the distance to different chromatin states, we identified new factors that impact the regulation of gene expression distributions.

Evaluating Enhancer Function and Transcription

Cell-type- and condition-specific profiles of gene expression require coordination between protein-coding gene promoters and cis-regulatory sequences called enhancers. Enhancers can stimulate gene activity at great genomic distances from their targets, raising questions about how enhancers communicate with specific gene promoters and what molecular mechanisms underlie enhancer function. Characterization of enhancer loci has identified the molecular features of active enhancers that accompany the binding of transcription factors and local opening of chromatin. These characteristics include coactivator recruitment, histone modifications, and noncoding RNA transcription. However, it remains unclear which of these features functionally contribute to enhancer activity. Here, we discuss what is known about how enhancers regulate their target genes and how enhancers and promoters communicate. Further, we describe recent data demonstrating many similarities between enhancers and the gene promoters they control, and we highlight unanswered questions in the field, such as the potential roles of transcription at enhancers.

Comprehensive In Vivo Interrogation Reveals Phenotypic Impact of Human Enhancer Variants

Establishing causal links between non-coding variants and human phenotypes is an increasing challenge. Here, we introduce a high-throughput mouse reporter assay for assessing the pathogenic potential of human enhancer variants in vivo and examine nearly a thousand variants in an enhancer repeatedly linked to polydactyly. We show that 71% of all rare non-coding variants previously proposed as causal lead to reporter gene expression in a pattern consistent with their pathogenic role. Variants observed to alter enhancer activity were further confirmed to cause polydactyly in knockin mice. We also used combinatorial and single-nucleotide mutagenesis to evaluate the in vivo impact of mutations affecting all positions of the enhancer and identified additional functional substitutions, including potentially pathogenic variants hitherto not observed in humans. Our results uncover the functional consequences of hundreds of mutations in a phenotype-associated enhancer and establish a widely applicable strategy for systematic in vivo evaluation of human enhancer variants.

CRISPR/Cas9: targeted genome editing for the treatment of hereditary hearing loss

Hereditary hearing loss (HHL) is a neurosensory disorder that affects every 1/500 newborns worldwide and nearly 1/3 people over the age of 65. Congenital deafness is inherited as monogenetic or polygenic disorder. The delicacy, tissue heterogeneity, deep location of the inner ear down the brainstem, and minute quantity of cells present in cochlea are the major challenges for current therapeutic approaches to cure deafness. Targeted genome editing is considered a suitable approach to treat HHL since it can target defective molecular components of auditory transduction to restore normal cochlear function. With the advent of CRISPR/Cas9 technique, targeted genome editing and biomedical research have been revolutionized. The robustness and simplicity of this technology lie in its design and delivery methods. It can directly deliver a complex of Cas9 endonuclease and single guide RNA (sgRNA) into zygote using either vector-mediated stable transfection or transient delivery of ribonucleoproteins complexes. This strategy induces DNA double strand breaks (DSBs) at target site followed by endogenous DNA repairing mechanisms of the cell. CRISPR/Cas9 has been successfully used in model animals to edit hearing genes like calcium and integrin-binding protein 2, myosin VIIA, Xin-actin binding repeat containing 2, leucine-zipper and sterile-alpha motif kinase Zak, epiphycan, transmembrane channel-like protein 1, and cadherin 23. This review discusses the utility of lipid-mediated transient delivery of Cas9/sgRNA complexes, an efficient way to restore hearing in humans, suffering from HHL. Notwithstanding, challenges like PAM requirement, HDR efficiency, off-target activity, and optimized delivery systems need to be addressed.

A highly soluble Sleeping Beauty transposase improves control of gene insertion

The Sleeping Beauty (SB) transposon system is an efficient non-viral gene transfer tool in mammalian cells but its broad use has been hampered by uncontrolled transposase gene activity from DNA vectors, posing a risk for genome instability, and by the inability to use transposase protein directly. Here, we used rational protein design based on the crystal structure of the hyperactive SB100X variant to create an SB transposase (hsSB) with enhanced solubility and stability. We demonstrate that hsSB can be delivered with transposon DNA to genetically modify cell lines and embryonic, hematopoietic and induced pluripotent stem cells (iPSCs), overcoming uncontrolled transposase activity. We used hsSB to generate chimeric antigen receptor (CAR) T-cells, which exhibit potent anti-tumor activity in vitro and in xenograft mice. We found that hsSB spontaneously penetrates cells, enabling modification of iPSCs and generation of CAR-T cells without the use of transfection reagents. Titration of hsSB to modulate genomic integration frequency achieved as few as two integrations per genome.

Zebrafish: A Powerful Model for Understanding the Functional Relevance of Noncoding Region Mutations in Human Genetic Diseases

Determining aetiology of genetic disorders caused by damaging mutations in protein-coding genes is well established. However, understanding how mutations in the vast stretches of the noncoding genome contribute to genetic abnormalities remains a huge challenge. Cis-regulatory elements (CREs) or enhancers are an important class of noncoding elements. CREs function as the primary determinants of precise spatial and temporal regulation of their target genes during development by serving as docking sites for tissue-specific transcription factors. Although a large number of potential disease-associated CRE mutations are being identified in patients, lack of robust methods for mechanistically linking these mutations to disease phenotype is currently hampering the understanding of their roles in disease aetiology. Here, we have described the various systems available for testing the CRE potential of stretches of noncoding regions harbouring mutations implicated in human disease. We highlight advances in the field leading to the establishment of zebrafish as a powerful system for robust and cost-effective functional assays of CRE activity, enabling rapid identification of causal variants in regulatory regions and the validation of their role in disruption of appropriate gene expression.

Functional dissection of the Sox9-Kcnj2 locus identifies nonessential and instructive roles of TAD architecture

The genome is organized in three-dimensional units called topologically associating domains (TADs), through a process dependent on the cooperative action of cohesin and the DNA-binding factor CTCF. Genomic rearrangements of TADs have been shown to cause gene misexpression and disease, but genome-wide depletion of CTCF has no drastic effects on transcription. Here, we investigate TAD function in vivo in mouse limb buds at the Sox9-Kcnj2 locus. We show that the removal of all major CTCF sites at the boundary and within the TAD resulted in a fusion of neighboring TADs, without major effects on gene expression. Gene misexpression and disease phenotypes, however, were achieved by redirecting regulatory activity through inversions and/or the repositioning of boundaries. Thus, TAD structures provide robustness and precision but are not essential for developmental gene regulation. Aberrant disease-related gene activation is not induced by a mere loss of insulation but requires CTCF-dependent redirection of enhancer-promoter contacts.

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.

Developmental dynamics of lncRNAs across mammalian organs and species

While many long noncoding RNAs (lncRNAs) have been identified in human and other mammalian genomes, there has been limited systematic functional characterization. In particular, the contribution of lncRNAs to organ development remains largely unexplored. Here we analyze the expression patterns of lncRNAs across developmental timepoints in seven major organs, from early organogenesis to adulthood, across seven species (human, macaque, mouse, rat, rabbit, opossum, and chicken). Our analyses identified ~15,000-35,000 candidate lncRNAs in each species, most of which show species specificity. We characterized expression patterns of lncRNAs across developmental stages, and found many with dynamic expression patterns across time that show signatures of enrichment for functionality. During development, there is a transition from broadly expressed and conserved lncRNAs towards an increasing number of lineage- and organ-specific lncRNAs. Our study provides a resource of candidate lncRNAs and their patterns of expression and evolutionary conservation across mammalian organ development.

A tour of 3D genome with a focus on CTCF

The complex three-dimensional (3D) structure of the genome plays critical roles in the maintenance of genome stability, organization, and dynamics and in regulation of gene expression for understanding molecular mechanisms and diseases. Chromatin maintains biological functions and transcriptional activities through long distance interaction and interactions between loops and enhancers-promoters. We firstly overview the architecture and biology of chromatin and loops, topologically associated domains (TADs) and interactions, and compartments and functions. We specifically focus on CCCTC-binding factor (CTCF) in 3D genome organization and function to furthermore understand the significance of CTCF biology, transcriptional regulations, interactions with cohesin, roles in DNA binding, influences of CTCF degradation, and communication with wings-apart like (Wapl) protein. We also summarize the advanced single cell approaches to further monitor dynamics of CTCF functions and structures in the maintenance of 3D genome organization and function at single cell level.

Generation of barcoded plasmid libraries for massively parallel analysis of chromatin position effects

The discovery of the position effect variegation phenomenon and the subsequent comprehensive analysis of its molecular mechanisms led to understanding that the local chromatin composition has a dramatic effect on gene activity. To study this effect in a high-throughput mode and at the genome-wide level, the Thousands of Reporters Integrated in Parallel (TRIP) approach based on the usage of barcoded reporter gene constructs was recently developed. Here we describe the construction and quality checks of high-diversity barcoded plasmid libraries supposed to be used for high-throughput analysis of chromatin position effects in Drosophila cells. First, we highlight the critical parameters that should be considered in the generation of barcoded plasmid libraries and introduce a simple method to assess the diversity of random sequences (barcodes) of synthetic oligonucleotides using PCR amplification followed by Sanger sequencing. Second, we compare the conventional restriction-ligation method with the Gibson assembly approach for cloning barcodes into the same plasmid vector. Third, we provide optimized parameters for the construction of barcoded plasmid libraries, such as the vector : insert ratio in the Gibson assembly reaction and the voltage used for electroporation of bacterial cells with ligation products. We also compare different approaches to check the quality of barcoded plasmid libraries. Finally, we briefly describe alternative approaches that can be used for the generation of such libraries. Importantly, all improvements and modifications of the techniques described here can be applied to a wide range of experiments involving barcoded plasmid libraries.

Functional characteristics of novel pancreatic Pax6 regulatory elements

Complex diseases, such as diabetes, are influenced by comprehensive transcriptional networks. Genome-wide association studies have revealed that variants located in regulatory elements for pancreatic transcription factors are linked to diabetes, including those functionally linked to the paired box transcription factor Pax6. Pax6 deletions in adult mice cause rapid onset of classic diabetes, but the full spectrum of pancreatic Pax6 regulators is unknown. Using a regulatory element discovery approach, we identified two novel Pax6 pancreatic cis-regulatory elements in a poorly characterized regulatory desert. Both new elements, Pax6 pancreas cis-regulatory element 3 (PE3) and PE4, are located 50 and 100 kb upstream and interact with different parts of the Pax6 promoter and nearby non-coding RNAs. They drive expression in the developing pancreas and brain and code for multiple pancreas-related transcription factor-binding sites. PE3 binds CCCTC-binding factor (CTCF) and is marked by stem cell identity markers in embryonic stem cells, whilst a common variant located in the PE4 element affects binding of Pax4, a known pancreatic regulator, altering Pax6 gene expression. To determine the ability of these elements to regulate gene expression, synthetic transcriptional activators and repressors were targeted to PE3 and PE4, modulating Pax6 gene expression, as well as influencing neighbouring genes and long non-coding RNAs, implicating the Pax6 locus in pancreas function and diabetes.

Enhancer Trapping and Annotation in Zebrafish Mediated with Sleeping Beauty, piggyBac and Tol2 Transposons

Although transposon-mediated enhancer trapping (ET) is successfully applied in diverse models, the efficiency of various transposon systems varies significantly, and little information is available regarding efficiency of enhancer trapping by various transposons in zebrafish. Most potential enhancers (Ens) still lack evidence of actual En activity. Here, we compared the differences in ET efficiency between sleeping beauty (SB), piggyBac (PB) and Tol2 transposons. Tol2 represented the highest germline transfer efficiencies at 55.56% (N_F0 = 165), followed by SB (38.36%, N_F0 = 151) and PB (32.65%, N_F0 = 149). ET lines generated by the Tol2 transposon tended to produce offspring with a single expression pattern per line, while PB and SB tended to generate embryos with multiple expression patterns. In our tests, 10 putative Ens (En1–10) were identified by splinkerette PCR and comparative genomic analysis. Combining the GFP expression profiles and mRNA expression patterns revealed that En1 and En2 may be involved in regulation of the expression of dlx1a and dlx2a, while En6 may be involved in regulation of the expression of line TK4 transgene and rps26, and En7 may be involved in the regulation of the expression of wnt1 and wnt10b. Most identified Ens were found to be transcribed in zebrafish embryos, and their regulatory function may involve eRNAs.

Three-dimensional organization and dynamics of the genome

Genome is a complex hierarchical structure, and its spatial organization plays an important role in its function. Chromatin loops and topological domains form the basic structural units of this multiscale organization and are essential to orchestrate complex regulatory networks and transcription mechanisms. They also form higher-order structures such as chromosomal compartments and chromosome territories. Each level of this intrinsic architecture is governed by principles and mechanisms that we only start to understand. In this review, we summarize the current view of the genome architecture on the scales ranging from chromatin loops to the whole genome. We describe cell-to-cell variability, links between genome reorganization and various genomic processes, such as chromosome X inactivation and cell differentiation, and the interplay between different experimental techniques.

Dynamic 3D chromatin architecture contributes to enhancer specificity and limb morphogenesis

The regulatory specificity of enhancers and their interaction with gene promoters is thought to be controlled by their sequence and the binding of transcription factors. By studying Pitx1, a regulator of hindlimb development, we show that dynamic changes in chromatin conformation can restrict the activity of enhancers. Inconsistent with its hindlimb-restricted expression, Pitx1 is controlled by an enhancer (Pen) that shows activity in forelimbs and hindlimbs. By Capture Hi-C and three-dimensional modeling of the locus, we demonstrate that forelimbs and hindlimbs have fundamentally different chromatin configurations, whereby Pen and Pitx1 interact in hindlimbs and are physically separated in forelimbs. Structural variants can convert the inactive into the active conformation, thereby inducing Pitx1 misexpression in forelimbs, causing partial arm-to-leg transformation in mice and humans. Thus, tissue-specific three-dimensional chromatin conformation can contribute to enhancer activity and specificity in vivo and its disturbance can result in gene misexpression and disease.

Rodent models in Down syndrome research: impact and future opportunities

Down syndrome is caused by trisomy of chromosome 21. To date, a multiplicity of mouse models with Down-syndrome-related features has been developed to understand this complex human chromosomal disorder. These mouse models have been important for determining genotype-phenotype relationships and identification of dosage-sensitive genes involved in the pathophysiology of the condition, and in exploring the impact of the additional chromosome on the whole genome. Mouse models of Down syndrome have also been used to test therapeutic strategies. Here, we provide an overview of research in the last 15 years dedicated to the development and application of rodent models for Down syndrome. We also speculate on possible and probable future directions of research in this fast-moving field. As our understanding of the syndrome improves and genome engineering technologies evolve, it is necessary to coordinate efforts to make all Down syndrome models available to the community, to test therapeutics in models that replicate the whole trisomy and design new animal models to promote further discovery of potential therapeutic targets.

Summary: Mouse models have boosted therapeutic options for Down syndrome, and improved models are being developed to better understand the pathophysiology of this genetic condition.

An evolutionary perspective of regulatory landscape dynamics in development and disease

The organization of animal genomes into topologically associating domains (TADs) provides a structural scaffold in which cis-regulatory elements (CREs) operate on their target genes. Determining the position of CREs and genes relative to TADs has become instrumental to trace gene expression changes during evolution and in diseases. Here we will review recent studies and discuss TADs as structural units with respect to their conservation and stability during genome reorganization. Furthermore, we describe how TAD restructuring contributed to morphological novelties during evolution but also their deleterious effects associated with disease. Despite considering TADs as structural units, the nested and dynamic scaffold within TADs contributes to tissue-specific gene expression, implying that such changes can also account for gene expression differences during evolution.

Developmentally regulated higher-order chromatin interactions orchestrate B cell fate commitment

Genome organization in 3D nuclear-space is important for regulation of gene expression. However, the alterations of chromatin architecture that impinge on the B cell-fate choice of multi-potent progenitors are still unclear. By integrating in situ Hi-C analyses with epigenetic landscapes and genome-wide expression profiles, we tracked the changes in genome architecture as the cells transit from a progenitor to a committed state. We identified the genomic loci that undergo developmental switch between A and B compartments during B-cell fate determination. Furthermore, although, topologically associating domains (TADs) are stable, a significant number of TADs display structural alterations that are associated with changes in cis-regulatory interaction landscape. Finally, we demonstrate the potential roles for Ebf1 and its downstream factor, Pax5, in chromatin reorganization and transcription regulation. Collectively, our studies provide a general paradigm of the dynamic relationship between chromatin reorganization and lineage-specific gene expression pattern that dictates cell-fate determination.

The three-dimensional genome: regulating gene expression during pluripotency and development

The precise expression of genes in time and space during embryogenesis is largely influenced by communication between enhancers and promoters, which is propagated and governed by the physical proximity of these elements in the nucleus. Here, we review how chromatin domains organize the genome by guiding enhancers to their target genes thereby preventing non-specific interactions with other neighboring regions. We also discuss the dynamics of chromatin interactions between enhancers and promoters, as well as the consequent changes in gene expression, that occur in pluripotent cells and during development. Finally, we evaluate how genomic changes such as deletions, inversions and duplications affect 3D chromatin configuration overall and lead to ectopic enhancer-promoter contacts, and thus gene misexpression, which can contribute to abnormal development and disease.

Between form and function: the complexity of genome folding

It has been known for over a century that chromatin is not randomly distributed within the nucleus. However, the question of how DNA is folded and the influence of such folding on nuclear processes remain topics of intensive current research. A longstanding, unanswered question is whether nuclear organization is simply a reflection of nuclear processes such as transcription and replication, or whether chromatin is folded by independent mechanisms and this per se encodes function? Evidence is emerging that both may be true. Here, using the α-globin gene cluster as an illustrative model, we provide an overview of the most recent insights into the layers of genome organization across different scales and how this relates to gene activity.

The functions of long noncoding RNAs in development and stem cells

Eukaryotic genomes are pervasively transcribed, with tens of thousands of RNAs emanating from uni- and bi-directional promoters and from active enhancers. In vertebrates, thousands of loci in each species produce a class of transcripts called long noncoding RNAs (lncRNAs) that are typically expressed at low levels and do not appear to give rise to functional proteins. Substantial numbers of lncRNAs are expressed at specific stages of embryonic development, in many cases from regions flanking key developmental regulators. Here, we review the known biological functions of such lncRNAs and the emerging paradigms of their modes of action. We also provide an overview of the growing arsenal of methods for lncRNA identification, perturbation and functional characterization.

Topological Domains, Metagenes, and the Emergence of Pleiotropic Regulations at Hox Loci

Hox gene clusters of jaw vertebrates display a tight genomic organization, which has no equivalent in any other bilateria genomes sequenced thus far. It was previously argued that such a topological consolidation toward a condensed, metagenic structure was due to the accumulation of long-range regulations flanking Hox loci, a phenomenon made possible by the successive genome duplications that occurred at the root of the vertebrate lineage, similar to gene neofunctionalization but applied to a coordinated multigenic system. Here, we propose that the emergence of such large vertebrate regulatory landscapes containing a range of global enhancers was greatly facilitated by the presence of topologically associating domains (TADs). These chromatin domains, mostly constitutive, may have been used as genomic niches where novel regulations could evolve due to both the preexistence of a structural backbone poised to integrate novel regulatory inputs, and a highly adaptive transcriptional readout. We propose a scenario for the coevolution of such TADs and the emergence of pleiotropy at ancestral vertebrate Hox loci.

Efficient chromosomal-scale DNA looping in Escherichia coli using multiple DNA-looping elements

Genes are frequently regulated by interactions between proteins that bind to the DNA near the gene and proteins that bind to DNA sites located far away, with the intervening DNA looped out. But it is not understood how efficient looping can occur when the sites are very far apart. We develop a simple theoretical framework that relates looping efficiency to the energetic cost and benefit of looping, allowing prediction of the efficiency of single or multiple nested loops at different distances. Measurements of absolute loop efficiencies for Lac repressor and λ CI using gene expression reporters in Escherichia coli cells show that, as predicted by the model, long-range DNA looping between a pair of sites can be strongly enhanced by the use of nested DNA loops or by the use of additional protein-binding sequences. A combination of these approaches was able to generate efficient DNA looping at a 200 kb distance.

Composition and dosage of a multipartite enhancer cluster control developmental expression of Ihh (Indian hedgehog)

Copy number variations (CNVs) often include non-coding sequence and putative enhancers but how these rearrangements induce disease is poorly understood. Here we investigate CNVs involving the regulatory landscape of Indian hedgehog (IHH), causing multiple, highly localised phenotypes including craniosynostosis and synpolydactyly^1,2. We show through transgenic reporter and genome editing studies in mice that Ihh is regulated by a constellation of at least 9 enhancers with individual tissue specificities in the digit anlagen, growth plates, skull sutures and fingertips. Consecutive deletions show that they function in an additive manner resulting in growth defects of the skull and long bones. Duplications, in contrast, cause not only dose-dependent upregulation but also misexpression of Ihh, leading to abnormal phalanges, fusion of sutures and syndactyly. Thus, precise spatio-temporal control of developmental gene expression is achieved by complex multipartite enhancer ensembles. Alterations in the composition of such clusters can result in gene misexpression and disease.

Modeling human disease in rodents by CRISPR/Cas9 genome editing

Modeling human disease has proven to be a challenge for the scientific community. For years, generating an animal model was complicated and restricted to very few species. With the rise of CRISPR/Cas9, it is now possible to generate more or less any animal model. In this review, we will show how this technology is and will change our way to obtain relevant disease animal models and how it should impact human health.

Krox20 hindbrain regulation incorporates multiple modes of cooperation between cis-acting elements

Developmental genes can harbour multiple transcriptional enhancers that act simultaneously or in succession to achieve robust and precise spatiotemporal expression. However, the mechanisms underlying cooperation between cis-acting elements are poorly documented, notably in vertebrates. The mouse gene Krox20 encodes a transcription factor required for the specification of two segments (rhombomeres) of the developing hindbrain. In rhombomere 3, Krox20 is subject to direct positive feedback governed by an autoregulatory enhancer, element A. In contrast, a second enhancer, element C, distant by 70 kb, is active from the initiation of transcription independent of the presence of the KROX20 protein. Here, using both enhancer knock-outs and investigations of chromatin organisation, we show that element C possesses a dual activity: besides its classical enhancer function, it is also permanently required in cis to potentiate the autoregulatory activity of element A, by increasing its chromatin accessibility. This work uncovers a novel, asymmetrical, long-range mode of cooperation between cis-acting elements that might be essential to avoid promiscuous activation of positive autoregulatory elements.

The formation of multicellular organisms from the egg to the adult stage is largely under genetic control. The activation of specific genes is governed by regulatory DNA sequences present nearby on the chromosome. Most of these sequences promote activation and are called enhancers. In this paper, we study two enhancers governing the expression of a gene involved in the formation of the posterior brain in vertebrates. One of these enhancers is involved in a positive feedback loop: it is itself activated by the protein product of the gene that it regulates. The other enhancer was thought to be only involved in the initial accumulation of the protein, necessary for the subsequent activation of the feedback loop. Here we show that the second enhancer directly cooperates with the autoregulatory enhancer to increase its accessibility and its activity. Our work uncovers a novel, long-range mode of cooperation between enhancers that restricts the domain of action of autoregulatory enhancers within embryos and might be essential to avoid their inappropriate activation.

The Role of Noncoding Genetic Variation in Isolated Orofacial Clefts

In the past decade, medical genetic research has generated multiple discoveries, many of which were obtained via genome-wide association studies (GWASs). A major GWAS finding is that the majority of risk variants for complex traits map to noncoding regions. This has resulted in a paradigm shift in terms of the interpretation of human genomic sequence variation, with more attention now being paid to what was previously termed “junk DNA.” Translation of genetic findings into biologically meaningful results requires 1) large-scale and cell-specific efforts to annotate non-protein–coding regions and 2) the integration of comprehensive genomic data sets. However, this represents an enormous challenge, particularly in the case of human traits that arise during embryonic development, such as orofacial clefts (OFCs). OFC is a multifactorial trait and ranks among the most common of all human congenital malformations. These 2 attributes apply in particular to its isolated forms (nonsyndromic OFC [nsOFC]). Although genetic studies (including GWASs) have yielded novel insights into the genetic architecture of nsOFC, few data are available concerning causality and affected biological pathways. Reasons for this deficiency include the complex genetic architecture at risk loci and the limited availability of functional data sets from human tissues that represent relevant embryonic sites and time points. The present review summarizes current knowledge of the role of noncoding regions in nsOFC etiology. We describe the identification of genetic risk factors for nsOFC and several of the approaches used to identify causal variants at these loci. These strategies include the use of biological and genetic information from public databases, the assessment of the full spectrum of genetic variability within 1 locus, and comprehensive in vitro and in vivo experiments. This review also highlights the role of the emerging research field “functional genomics” and its increasing contribution to our biological understanding of nsOFC.

Recent advances in high-throughput approaches to dissect enhancer function

The regulation of gene transcription in higher eukaryotes is accomplished through the involvement of transcription start site (TSS)-proximal (promoters) and -distal (enhancers) regulatory elements. It is now well acknowledged that enhancer elements play an essential role during development and cell differentiation, while genetic alterations in these elements are a major cause of human disease. Many strategies have been developed to identify and characterize enhancers. Here, we discuss recent advances in high-throughput approaches to assess enhancer activity, from the well-established massively parallel reporter assays to the recent clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9-based technologies. We highlight how these approaches contribute toward a better understanding of enhancer function, eventually leading to the discovery of new types of regulatory sequences, and how the alteration of enhancers can affect transcriptional regulation.

Global reorganisation of cis-regulatory units upon lineage commitment of human embryonic stem cells

Long-range cis-regulatory elements such as enhancers coordinate cell-specific transcriptional programmes by engaging in DNA looping interactions with target promoters. Deciphering the interplay between the promoter connectivity and activity of cis-regulatory elements during lineage commitment is crucial for understanding developmental transcriptional control. Here, we use Promoter Capture Hi-C to generate a high-resolution atlas of chromosomal interactions involving ~22,000 gene promoters in human pluripotent and lineage-committed cells, identifying putative target genes for known and predicted enhancer elements. We reveal extensive dynamics of cis-regulatory contacts upon lineage commitment, including the acquisition and loss of promoter interactions. This spatial rewiring occurs preferentially with predicted changes in the activity of cis-regulatory elements and is associated with changes in target gene expression. Our results provide a global and integrated view of promoter interactome dynamics during lineage commitment of human pluripotent cells.

DOI:http://dx.doi.org/10.7554/eLife.21926.001

Virtually every cell in the body contains the same set of DNA, which encodes thousands of genes. The activities of these genes vary between different types of cells and at different points in time. As a result, our DNA contains a complex array of molecular switches that instruct genes to switch on and off at the right time and in the right cells. These molecular switches, termed regulatory elements, are often a long way away from the genes that they control, and this can make it difficult to find out which switch controls which genes.

DNA is made up of several different building blocks known as bases and the order of these bases encodes specific information about the gene. Every human cell contains approximately two meters of DNA, which is highly folded in the cell nucleus. This three-dimensional folding allows regions that are far apart on the DNA thread to physically contact each other. To reach the genes they control, regulatory elements form loops on the DNA that are near-impossible to predict from looking at the sequence of bases alone. Mapping the locations of these loops can reveal the hidden circuitry within our DNA and help us to understand how unwanted changes (mutations) within regulatory elements may cause disease.

Freire-Pritchett, Schoenfelder et al. set out to reveal how loops between genes and their regulatory elements change as the stem cells specialise into immature brain cells. The experiments show that the pattern of DNA loops is extensively altered after the stem cells specialise into brain cells, that is, some loops are lost and new ones form. Furthermore, the regulatory elements themselves were often toggled between “on” and “off” states. These two processes tend to occur together, so that new loops are formed at the same time as the switch is activated.

Freire-Pritchett, Schoenfelder et al. also show that individual genes are often connected to many different regulatory elements. The next step is to understand how these multiple connections work together to coordinate gene activity, and whether this information could be used to control how stem cells specialise. This knowledge may lead to the development of stem cell-based therapies for stroke, Parkinson’s disease and other conditions.

DOI:http://dx.doi.org/10.7554/eLife.21926.002

Efficient and rapid generation of large genomic variants in rats and mice using CRISMERE

Modelling Down syndrome (DS) in mouse has been crucial for the understanding of the disease and the evaluation of therapeutic targets. Nevertheless, the modelling so far has been limited to the mouse and, even in this model, generating duplication of genomic regions has been labour intensive and time consuming. We developed the CRISpr MEdiated REarrangement (CRISMERE) strategy, which takes advantage of the CRISPR/Cas9 system, to generate most of the desired rearrangements from a single experiment at much lower expenses and in less than 9 months. Deletions, duplications, and inversions of genomic regions as large as 24.4 Mb in rat and mouse founders were observed and germ line transmission was confirmed for fragment as large as 3.6 Mb. Interestingly we have been able to recover duplicated regions from founders in which we only detected deletions. CRISMERE is even more powerful than anticipated it allows the scientific community to manipulate the rodent and probably other genomes in a fast and efficient manner which was not possible before.

[High-throughput approaches to study cis-regulating elements]

Gene expression in higher eukaryotes is regulated through the involvement of transcription start site (TSS)-proximal (promoters) and -distal (enhancers) regulatory elements. Enhancer elements play an essential role during development and cell differentiation, while genetic alterations in these elements are a major cause of human disease. Here, we discuss recent advances in high-throughput approaches to identify and characterize enhancer elements, from the well-established massively parallel reporter assays to the recent clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9-based technologies. We discuss how these approaches contribute toward a better understanding of enhancer function in normal and pathological conditions.