We developed the Genomic Regions Enrichment of Annotations Tool (GREAT) to analyze the functional significance of cis-regulatory regions identified by localized measurements of DNA binding events across an entire genome. Whereas previous methods took into account only binding proximal to genes, GREAT is able to properly incorporate distal binding sites and control for false positives using a binomial test over the input genomic regions. GREAT incorporates annotations from 20 ontologies and is available as a web application. Applying GREAT to data sets from chromatin immunoprecipitation coupled with massively parallel sequencing (ChIP-seq) of multiple transcription-associated factors, including SRF, NRSF, GABP, Stat3 and p300 in different developmental contexts, we recover many functions of these factors that are missed by existing gene-based tools, and we generate testable hypotheses. The utility of GREAT is not limited to ChIP-seq, as it could also be applied to open chromatin, localized epigenomic markers and similar functional data sets, as well as comparative genomics sets.

/pdf/great-improves-functional-interpretation-of-cis-regulatory-47gi1yqk9r.pdf

GREAT improves functional interpretation of cis-regulatory regions

Genome-wide association studies have identified many noncoding variants associated with common diseases and traits. We show that these variants are concentrated in regulatory DNA marked by deoxyribonuclease I (DNase I) hypersensitive sites (DHSs). Eighty-eight percent of such DHSs are active during fetal development and are enriched in variants associated with gestational exposure–related phenotypes. We identified distant gene targets for hundreds of variant-containing DHSs that may explain phenotype associations. Disease-associated variants systematically perturb transcription factor recognition sequences, frequently alter allelic chromatin states, and form regulatory networks. We also demonstrated tissue-selective enrichment of more weakly disease-associated variants within DHSs and the de novo identification of pathogenic cell types for Crohn’s disease, multiple sclerosis, and an electrocardiogram trait, without prior knowledge of physiological mechanisms. Our results suggest pervasive involvement of regulatory DNA variation in common human disease and provide pathogenic insights into diverse disorders.

/pdf/systematic-localization-of-common-disease-associated-2h35ncsl4r.pdf

Systematic localization of common disease-associated variation in regulatory DNA.

Traditional approaches to measuring the complexity of biological signals fail to account for the multiple time scales inherent in such time series. These algorithms have yielded contradictory findings when applied to real-world datasets obtained in health and disease states. We describe in detail the basis and implementation of the multiscale entropy (MSE) method. We extend and elaborate previous findings showing its applicability to the fluctuations of the human heartbeat under physiologic and pathologic conditions. The method consistently indicates a loss of complexity with aging, with an erratic cardiac arrhythmia (atrial fibrillation), and with a life-threatening syndrome (congestive heart failure). Further, these different conditions have distinct MSE curve profiles, suggesting diagnostic uses. The results support a general "complexity-loss" theory of aging and disease. We also apply the method to the analysis of coding and noncoding DNA sequences and find that the latter have higher multiscale entropy, consistent with the emerging view that so-called "junk DNA" sequences contain important biological information.

/pdf/multiscale-entropy-analysis-of-biological-signals-ho0qlfaxad.pdf

Multiscale entropy analysis of biological signals

https://www.cell.com/cell/pdf/S0092-8674(15)00377-3.pdf

Disruptions of Topological Chromatin Domains Cause Pathogenic Rewiring of Gene-Enhancer Interactions

Gene targeting in embryonic stem cells has become the principal technology for manipulation of the mouse genome, offering unrivalled accuracy in allele design and access to conditional mutagenesis. To bring these advantages to the wider research community, large-scale mouse knockout programmes are producing a permanent resource of targeted mutations in all protein-coding genes. Here we report the establishment of a high-throughput gene-targeting pipeline for the generation of reporter-tagged, conditional alleles. Computational allele design, 96-well modular vector construction and high-efficiency gene-targeting strategies have been combined to mutate genes on an unprecedented scale. So far, more than 12,000 vectors and 9,000 conditional targeted alleles have been produced in highly germline-competent C57BL/6N embryonic stem cells. High-throughput genome engineering highlighted by this study is broadly applicable to rat and human stem cells and provides a foundation for future genome-wide efforts aimed at deciphering the function of all genes encoded by the mammalian genome.

/pdf/a-conditional-knockout-resource-for-the-genome-wide-study-of-1inxfz9xq3.pdf

A conditional knockout resource for the genome-wide study of mouse gene function.

Unequivocal identification of the full composition of a gene is made difficult by the cryptic nature of regulatory elements. Regulatory elements are notoriously difficult to locate and may reside at considerable distances from the transcription units on which they operate and, moreover, may be incorporated into the structure of neighbouring genes. The importance of regulatory mutations as the basis of human abnormalities remains obscure. Here, we show that the chromosome 7q36 associated preaxial polydactyly, a frequently observed congenital limb malformation, results from point mutations in a Shh regulatory element. Shh, normally expressed in the ZPA posteriorly in the limb bud, is expressed in an additional ectopic site at the anterior margin in mouse models of PPD. Our investigations into the basis of the ectopic Shh expression identified the enhancer element that drives normal Shh expression in the ZPA. The regulator, designated ZRS, lies within intron 5 of the Lmbr1 gene 1 Mb from the target gene Shh. The ZRS drives the early spatio-temporal expression pattern in the limb of tetrapods. Despite the morphological differences between limbs and fins, an equivalent regulatory element is found in fish. The ZRS contains point mutations that segregate with polydactyly in four unrelated families with PPD and in the Hx mouse mutant. Thus point mutations residing in long-range regulatory elements are capable of causing congenital abnormalities, and possess the capacity to modify gene activity such that a novel gamut of abnormalities is detected.

/pdf/a-long-range-shh-enhancer-regulates-expression-in-the-1lf7f9g0ii.pdf

A long-range Shh enhancer regulates expression in the developing limb and fin and is associated with preaxial polydactyly

Preaxial polydactyly (PPD) is a common limb malformation in human. A number of polydactylous mouse mutants indicate that misexpression of Shh is a common requirement for generating extra digits. Here we identify a translocation breakpoint in a PPD patient and a transgenic insertion site in the polydactylous mouse mutant sasquatch (Ssq). The genetic lesions in both lie within the same respective intron of the LMBR1/Lmbr1 gene, which resides ≈1 Mb away from Shh. Genetic analysis of Ssq reveals that the Lmbr1 gene is incidental to the phenotype and that the mutation directly interrupts a cis-acting regulator of Shh. This regulator is most likely the target for generating PPD mutations in human.

https://www.pnas.org/content/pnas/99/11/7548.full.pdf

Disruption of a long-range cis-acting regulator for Shh causes preaxial polydactyly.

Sonic hedgehog, SHH, is required for patterning the limb. The array of skeletal elements that compose the hands and feet, and the ordered arrangement of these bones to form the pattern of fingers and toes are dependent on SHH. The mechanism of action of SHH in the limb is not fully understood; however, an aspect that appears to be important is the localized, asymmetric expression of Shh. Shh is expressed in the posterior margin of the limb bud in a region defined as the zone of polarizing activity (ZPA). Analysis of mouse mutants which have polydactyly (extra toes) shows that asymmetric expression of Shh is lost due to the appearance of an ectopic domain of expression in the anterior limb margin. One such polydactylous mouse mutant, sasquatch (Ssq), maps to the corresponding chromosomal region of the human condition pre-axial polydactyly (PPD) and thus represents a model for this condition. The mutation responsible for Ssq is located 1 Mb away from the Shh gene; however, the mutation disrupts a long-range cis-acting regulator of Shh expression. By inference, human pre-axial polydactyly results from a similar disruption of Shh expression. Other human congenital abnormalities also map near the pre-axial polydactyly locus, suggesting a major chromosomal region for limb dysmorphologies. The distinct phenotypes range from loss of all bones of the hands and feet to syndactyly of the soft tissue and fusion of the digits. We discuss the role played by Shh expression in mouse mutant phenotypes and the human limb dysmorphologies.

Sonic hedgehog: restricted expression and limb dysmorphologies.

In mouse, a number of different polydactylous mutants indicate that the ultimate misexpression of Shh is a common requirement for generating extra digits. We show that one of these mutants, called sasquatch (Ssq), is the counterpart to human preaxial polydactyly (PPD). We identify a translocation breakpoint in a PPD patient and the transgenic insertion site in the Ssq mutant. The genetic lesion in both mouse and human are located within the same intron of a gene termed LMBR1/Lmbr1 (Clark et al. 2000) situated approximately 1 Mb from Shh. We previously showed (Sharpe et al. 1999) that the reporter gene contained within the Ssq transgene insertion has acquired a limb-specific Shh-like expression pattern, not only in the appropriate posterior ZPA but also at an ectopic site at the anterior margin of the limb bud. It was unclear whether the mechanism for Shh misexpression involves disruption of the Lmbr1 gene and is therefore indirect or direct due to disruption of a long range Shh regulatory element. In order to distinguish the 2 possibilities we established a genetic cis trans test which revealed that the Ssq mutation interrupts a long-range cis-acting regulator of Shh ˜1 Mb away. Thus PPD in human most likely results from similar acting mutations and supports the prospect that Shh regulatory elements underlie this complex locus. Subsequently we have shown that Shh regulatory information resides within this locus and we present a model to explain other limb-related defects that map to this locus.

Simon J. H. Heaney

Papers

A long-range Shh enhancer regulates expression in the developing limb and fin and is associated with preaxial polydactyly

Disruption of a long-range cis-acting regulator for Shh causes preaxial polydactyly.

Sonic hedgehog: restricted expression and limb dysmorphologies.

2 Preaxial polydactyly In human and mouse: regulatory anomalies in digit patterning.