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

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

Synthetic STARR-seq reveals how DNA shape and sequence modulate transcriptional output and noise.

The binding of transcription factors to short recognition sequences plays a pivotal role in controlling the expression of genes. The sequence and shape characteristics of binding sites influence DNA binding specificity and have also been implicated in modulating the activity of transcription factors downstream of binding. To quantitatively assess the transcriptional activity of tens of thousands of designed synthetic sites in parallel, we developed a synthetic version of STARR-seq (synSTARR-seq). We used the approach to systematically analyze how variations in the recognition sequence of the glucocorticoid receptor (GR) affect transcriptional regulation. Our approach resulted in the identification of a novel highly active functional GR binding sequence and revealed that sequence variation both within and flanking GR's core binding site can modulate GR activity without apparent changes in DNA binding affinity. Notably, we found that the sequence composition of variants with similar activity profiles was highly diverse. In contrast, groups of variants with similar activity profiles showed specific DNA shape characteristics indicating that DNA shape may be a better predictor of activity than DNA sequence. Finally, using single cell experiments with individual enhancer variants, we obtained clues indicating that the architecture of the response element can independently tune expression mean and cell-to cell variability in gene expression (noise). Together, our studies establish synSTARR as a powerful method to systematically study how DNA sequence and shape modulate transcriptional output and noise.

Expanding the repertoire of glucocorticoid receptor target genes by engineering genomic response elements.

The glucocorticoid receptor (GR), a hormone-activated transcription factor, binds to a myriad of genomic binding sites yet seems to regulate a much smaller number of genes. Genome-wide analysis of GR binding and gene regulation has shown that the likelihood of GR-dependent regulation increases with decreased distance of its binding to the transcriptional start site of a gene. To test if we can adopt this knowledge to expand the repertoire of GR target genes, we used CRISPR/Cas-mediated homology-directed repair to add a single GR-binding site directly upstream of the transcriptional start site of each of four genes. To our surprise, we found that the addition of a single GR-binding site can be enough to convert a gene into a GR target. The gain of GR-dependent regulation was observed for two of four genes analyzed and coincided with acquired GR binding at the introduced binding site. However, the gene-specific gain of GR-dependent regulation could not be explained by obvious differences in chromatin accessibility between converted genes and their non-converted counterparts. Furthermore, by introducing GR-binding sequences with different nucleotide compositions, we show that activation can be facilitated by distinct sequences without obvious differences in activity between the GR-binding sequence variants we tested. The approach to use genome engineering to build genomic response elements facilitates the generation of cell lines with tailored repertoires of GR-responsive genes and a framework to test and refine our understanding of the cis-regulatory logic of gene regulation by testing if engineered response elements behave as predicted.

Glucocorticoid resistance of allogeneic T cells alters the gene expression profile in the inflamed small intestine of mice suffering from acute graft-versus-host disease.

Glucocorticoids (GCs) play an important role in controlling acute graft-versus-host disease (aGvHD), a frequent complication of allogeneic hematopoietic stem cell transplantation. The anti-inflammatory activity of GCs is mainly ascribed to the modulation of T cells and macrophages, for which reason a genetically induced GC resistance of either of these cell types causes aggravated aGvHD. Since only a few genes are currently known that are differentially regulated under these conditions, we analyzed the expression of 54 candidate genes in the inflamed small intestine of mice suffering from aGvHD when either allogeneic T cells or host myeloid cells were GC resistant using a microfluidic dynamic array platform for high-throughput quantitative PCR. The majority of genes categorized as cytokines (e.g. Il2, Il6), chemokines (e.g. Ccl2, Cxcl1), cell surface receptors (e.g. Fasl, Ctla4) and intracellular molecules (e.g. Dusp1, Arg1) were upregulated in mice transplanted with GC resistant allogeneic T cells. Moreover, the expression of several genes linked to energy metabolism (e.g. Glut1) was altered. Surprisingly, mice harboring GC resistant myeloid cells showed almost no changes in gene expression despite their fatal disease course after aGvHD induction. To identify additional genes in the inflamed small intestine that were affected by a GC resistance of allogeneic T cells, we performed an RNAseq analysis, which uncovered more than 500 differentially expressed transcripts (e.g. Cxcr6, Glut3, Otc, Aoc1, Il1r1, Sphk1) that were enriched for biological processes associated with inflammation and tissue disassembly. The changes in gene expression could be confirmed during full-blown disease but hardly any of them in the preclinical phase using high-throughput quantitative PCR. Further analysis of some of these genes revealed a highly selective expression pattern in T cells, intestinal epithelial cells and macrophages, which correlated with their regulation during disease progression. Collectively, we identified an altered gene expression profile caused by GC resistance of transplanted allogeneic T cells, which could help to define new targets for aGvHD therapy.

Deletions, Inversions, Duplications: Engineering of Structural Variants using CRISPR/Cas in Mice.

Structural variations (SVs) contribute to the variability of our genome and are often associated with disease. Their study in model systems was hampered until now by labor-intensive genetic targeting procedures and multiple mouse crossing steps. Here we present the use of CRISPR/Cas for the fast (10 weeks) and efficient generation of SVs in mice. We specifically produced deletions, inversions, and also duplications at six different genomic loci ranging from 1.1 kb to 1.6 Mb with efficiencies up to 42%. After PCR-based selection, clones were successfully used to create mice via aggregation. To test the practicability of the method, we reproduced a human 500 kb disease-associated deletion and were able to recapitulate the human phenotype in mice. Furthermore, we evaluated the regulatory potential of a large genomic interval by deleting a 1.5 Mb fragment. The method presented permits rapid in vivo modeling of genomic rearrangements.