Verena Heinrich

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.

Efficient CRISPR/Cas9 genome editing with low off-target effects in zebrafish

Gene modifications in animal models have been greatly facilitated through the application of targeted genome editing tools. The prokaryotic CRISPR/Cas9 type II genome editing system has recently been applied in cell lines and vertebrates. However, we still have very limited information about the efficiency of mutagenesis, germline transmission rates and off-target effects in genomes of model organisms. We now demonstrate that CRISPR/Cas9 mutagenesis in zebrafish is highly efficient, reaching up to 86.0%, and is heritable. The efficiency of the CRISPR/Cas9 system further facilitated the targeted knock-in of a protein tag provided by a donor oligonucleotide with knock-in efficiencies of 3.5-15.6%. Mutation rates at potential off-target sites are only 1.1-2.5%, demonstrating the specificity of the CRISPR/Cas9 system. The ease and efficiency of the CRISPR/Cas9 system with limited off-target effects make it a powerful genome engineering tool for in vivo studies.

Formation of new chromatin domains determines pathogenicity of genomic duplications

Chromosome conformation capture methods have identified subchromosomal structures of higher-order chromatin interactions called topologically associated domains (TADs) that are separated from each other by boundary regions. By subdividing the genome into discrete regulatory units, TADs restrict the contacts that enhancers establish with their target genes. However, the mechanisms that underlie partitioning of the genome into TADs remain poorly understood. Here we show by chromosome conformation capture (capture Hi-C and 4C-seq methods) that genomic duplications in patient cells and genetically modified mice can result in the formation of new chromatin domains (neo-TADs) and that this process determines their molecular pathology. Duplications of non-coding DNA within the mouse Sox9 TAD (intra-TAD) that cause female to male sex reversal in humans, showed increased contact of the duplicated regions within the TAD, but no change in the overall TAD structure. In contrast, overlapping duplications that extended over the next boundary into the neighbouring TAD (inter-TAD), resulted in the formation of a new chromatin domain (neo-TAD) that was isolated from the rest of the genome. As a consequence of this insulation, inter-TAD duplications had no phenotypic effect. However, incorporation of the next flanking gene, Kcnj2, in the neo-TAD resulted in ectopic contacts of Kcnj2 with the duplicated part of the Sox9 regulatory region, consecutive misexpression of Kcnj2, and a limb malformation phenotype. Our findings provide evidence that TADs are genomic regulatory units with a high degree of internal stability that can be sculptured by structural genomic variations. This process is important for the interpretation of copy number variations, as these variations are routinely detected in diagnostic tests for genetic disease and cancer. This finding also has relevance in an evolutionary setting because copy-number differences are thought to have a crucial role in the evolution of genome complexity.

Filtering for Compound Heterozygous Sequence Variants in Non-Consanguineous Pedigrees

The identification of disease-causing mutations in next-generation sequencing (NGS) data requires efficient filtering techniques. In patients with rare recessive diseases, compound heterozygosity of pathogenic mutations is the most likely inheritance model if the parents are non-consanguineous. We developed a web-based compound heterozygous filter that is suited for data from NGS projects and that is easy to use for non-bioinformaticians. We analyzed the power of compound heterozygous mutation filtering by deriving background distributions for healthy individuals from different ethnicities and studied the effectiveness in trios as well as more complex pedigree structures. While usually more then 30 genes harbor potential compound heterozygotes in single exomes, this number can be markedly reduced with every additional member of the pedigree that is included in the analysis. In a real data set with exomes of four family members, two sisters affected by Mabry syndrome and their healthy parents, the disease-causing gene PIGO, which harbors the pathogenic compound heterozygous variants, could be readily identified. Compound heterozygous filtering is an efficient means to reduce the number of candidate mutations in studies aiming at identifying recessive disease genes in non-consanguineous families. A web-server is provided to make this filtering strategy available at www.gene-talk.de.

The allele distribution in next-generation sequencing data sets is accurately described as the result of a stochastic branching process

With the availability of next-generation sequencing (NGS) technology, it is expected that sequence variants may be called on a genomic scale. Here, we demonstrate that a deeper understanding of the distribution of the variant call frequencies at heterozygous loci in NGS data sets is a prerequisite for sensitive variant detection. We model the crucial steps in an NGS protocol as a stochastic branching process and derive a mathematical framework for the expected distribution of alleles at heterozygous loci before measurement that is sequencing. We confirm our theoretical results by analyzing technical replicates of human exome data and demonstrate that the variance of allele frequencies at heterozygous loci is higher than expected by a simple binomial distribution. Due to this high variance, mutation callers relying on binomial distributed priors are less sensitive for heterozygous variants that deviate strongly from the expected mean frequency. Our results also indicate that error rates can be reduced to a greater degree by technical replicates than by increasing sequencing depth.

Characterization of hundreds of regulatory landscapes in developing limbs reveals two regimes of chromatin folding

Complex regulatory landscapes control the pleiotropic transcriptional activities of developmental genes. For most genes, the number, location, and dynamics of their associated regulatory elements are unknown. In this work, we characterized the three-dimensional chromatin microarchitecture and regulatory landscape of 446 limb-associated gene loci in mouse using Capture-C, ChIP-seq, and RNA-seq in forelimb, hindlimb at three developmental stages, and midbrain. The fine mapping of chromatin interactions revealed a strong preference for functional genomic regions such as repressed or active domains. By combining chromatin marks and interaction peaks, we annotated more than 1000 putative limb enhancers and their associated genes. Moreover, the analysis of chromatin interactions revealed two regimes of chromatin folding, one producing interactions stable across tissues and stages and another one associated with tissue and/or stage-specific interactions. Whereas stable interactions associate strongly with CTCF/RAD21 binding, the intensity of variable interactions correlates with changes in underlying chromatin modifications, specifically at the viewpoint and at the interaction site. In conclusion, this comprehensive data set provides a resource for the characterization of hundreds of limb-associated regulatory landscapes and a framework to interpret the chromatin folding dynamics observed during embryogenesis.

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

Copy number variations (CNVs) often include noncoding sequences and putative enhancers, but how these rearrangements induce disease is poorly understood. Here we investigate CNVs involving the regulatory landscape of IHH (encoding Indian hedgehog), which cause multiple, highly localized phenotypes including craniosynostosis and synpolydactyly. We show through transgenic reporter and genome-editing studies in mice that Ihh is regulated by a constellation of at least nine enhancers with individual tissue specificities in the digit anlagen, growth plates, skull sutures and fingertips. Consecutive deletions, resulting in growth defects of the skull and long bones, showed that these enhancers function in an additive manner. Duplications, in contrast, caused not only dose-dependent upregulation but also misexpression of Ihh, leading to abnormal phalanges, fusion of sutures and syndactyly. Thus, precise spatiotemporal 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.

Estimating exome genotyping accuracy by comparing to data from large scale sequencing projects

With exome sequencing becoming a tool for mutation detection in routine diagnostics there is an increasing need for platform-independent methods of quality control. We present a genotype-weighted metric that allows comparison of all the variant calls of an exome to a high-quality reference dataset of an ethnically matched population. The exome-wide genotyping accuracy is estimated from the distance to this reference set, and does not require any further knowledge about data generation or the bioinformatics involved. The distances of our metric are visualized by non-metric multidimensional scaling and serve as an intuitive, standardizable score for the quality assessment of exome data.

Screening for single nucleotide variants, small indels and exon deletions with a next-generation sequencing based gene panel approach for Usher syndrome

Usher syndrome is an autosomal recessive disorder characterized both by deafness and blindness. For the three clinical subtypes of Usher syndrome causal mutations in altogether 12 genes and a modifier gene have been identified. Due to the genetic heterogeneity of Usher syndrome, the molecular analysis is predestined for a comprehensive and parallelized analysis of all known genes by next‐generation sequencing (NGS) approaches. We describe here the targeted enrichment and deep sequencing for exons of Usher genes and compare the costs and workload of this approach compared to Sanger sequencing. We also present a bioinformatics analysis pipeline that allows us to detect single‐nucleotide variants, short insertions and deletions, as well as copy number variations of one or more exons on the same sequence data. Additionally, we present a flexible in silico gene panel for the analysis of sequence variants, in which newly identified genes can easily be included. We applied this approach to a cohort of 44 Usher patients and detected biallelic pathogenic mutations in 35 individuals and monoallelic mutations in eight individuals of our cohort. Thirty‐nine of the sequence variants, including two heterozygous deletions comprising several exons of USH2A, have not been reported so far. Our NGS‐based approach allowed us to assess single‐nucleotide variants, small indels, and whole exon deletions in a single test. The described diagnostic approach is fast and cost‐effective with a high molecular diagnostic yield.

Homozygous and Compound-Heterozygous Mutations in TGDS Cause Catel-Manzke Syndrome

Catel-Manzke syndrome is characterized by Pierre Robin sequence and a unique form of bilateral hyperphalangy causing a clinodactyly of the index finger. We describe the identification of homozygous and compound heterozygous mutations in TGDS in seven unrelated individuals with typical Catel-Manzke syndrome by exome sequencing. Six different TGDS mutations were detected: c.892A>G (p.Asn298Asp), c.270_271del (p.Lys91Asnfs(∗)22), c.298G>T (p.Ala100Ser), c.294T>G (p.Phe98Leu), c.269A>G (p.Glu90Gly), and c.700T>C (p.Tyr234His), all predicted to be disease causing. By using haplotype reconstruction we showed that the mutation c.298G>T is probably a founder mutation. Due to the spectrum of the amino acid changes, we suggest that loss of function in TGDS is the underlying mechanism of Catel-Manzke syndrome. TGDS (dTDP-D-glucose 4,6-dehydrogenase) is a conserved protein belonging to the SDR family and probably plays a role in nucleotide sugar metabolism.