Response to comment on "Human-specific gain of function in a developmental enhancer"

Duret and Galtier argue that human-specific sequence divergence and gain of function in the HACNS1 enhancer result from deleterious biased gene conversion (BGC) with no contribution from positive selection. We reinforce our previous conclusion by analyzing hypothesized BGC events genomewide and assessing the effect of recombination rates on human-accelerated conserved noncoding sequence ascertainment. We also provide evidence that AT → GC substitution bias can coexist with positive selection

Large-scale discovery of enhancers from human heart tissue.

Development and function of the human heart depend on the dynamic control of tissue-specific gene expression by distant-acting transcriptional enhancers. To generate an accurate genome-wide map of human heart enhancers, we used an epigenomic enhancer discovery approach and identified ∼6,200 candidate enhancer sequences directly from fetal and adult human heart tissue. Consistent with their predicted function, these elements were markedly enriched near genes implicated in heart development, function and disease. To further validate their in vivo enhancer activity, we tested 65 of these human sequences in a transgenic mouse enhancer assay and observed that 43 (66%) drove reproducible reporter gene expression in the heart. These results support the discovery of a genome-wide set of noncoding sequences highly enriched in human heart enhancers that is likely to facilitate downstream studies of the role of enhancers in development and pathological conditions of the heart

Functional autonomy of distant-acting human enhancers

Many human genes are associated with dispersed arrays of transcriptional enhancers that regulate their expression in time and space. Studies in invertebrate model systems have suggested that these elements function as discrete and independent regulatory units, but the in vivo combinatorial properties of vertebrate enhancers remain poorly understood. To explore the modularity and regulatory autonomy of human developmental enhancers, we experimentally concatenated up to four enhancers from different genes and used a transgenic mouse assay to compare the in vivo activity of these compound elements with that of the single modules. In all of the six different combinations of elements tested, the reporter gene activity patterns were additive without signs of interference between the individual modules, indicating that regulatory specificity was maintained despite the presence of closely-positioned heterologous enhancers. Even in cases where two elements drove expression in close anatomical proximity, such as within neighboring subregions of the developing limb bud, the compound patterns did not show signs of cross-inhibition between individual elements or novel expression sites. These data indicate that human developmental enhancers are highly modular and functionally autonomous and suggest that genomic enhancer shuffling may have contributed to the evolution of complex gene expression patterns in vertebrate

Deletion of Ultraconserved Elements Yields Viable Mice

Ultraconserved elements have been suggested to retain extended perfect sequence identity between the human, mouse, and rat genomes due to essential functional properties. To investigate the necessities of these elements in vivo, we removed four noncoding ultraconserved elements (ranging in length from 222 to 731 base pairs) from the mouse genome. To maximize the likelihood of observing a phenotype, we chose to delete elements that function as enhancers in a mouse transgenic assay and that are near genes that exhibit marked phenotypes both when completely inactivated in the mouse and when their expression is altered due to other genomic modifications. Remarkably, all four resulting lines of mice lacking these ultraconserved elements were viable and fertile, and failed to reveal any critical abnormalities when assayed for a variety of phenotypes including growth, longevity, pathology, and metabolism. In addition, more targeted screens, informed by the abnormalities observed in mice in which genes in proximity to the investigated elements had been altered, also failed to reveal notable abnormalities. These results, while not inclusive of all the possible phenotypic impact of the deleted sequences, indicate that extreme sequence constraint does not necessarily reflect crucial functions required for viability

ChIP-seq Accurately Predicts Tissue-Specific Activity of Enhancers

A major yet unresolved quest in decoding the human genome is the identification of the regulatory sequences that control the spatial and temporal expression of genes. Distant-acting transcriptional enhancers are particularly challenging to uncover since they are scattered amongst the vast non-coding portion of the genome. Evolutionary sequence constraint can facilitate the discovery of enhancers, but fails to predict when and where they are active in vivo. Here, we performed chromatin immunoprecipitation with the enhancer-associated protein p300, followed by massively-parallel sequencing, to map several thousand in vivo binding sites of p300 in mouse embryonic forebrain, midbrain, and limb tissue. We tested 86 of these sequences in a transgenic mouse assay, which in nearly all cases revealed reproducible enhancer activity in those tissues predicted by p300 binding. Our results indicate that in vivo mapping of p300 binding is a highly accurate means for identifying enhancers and their associated activities and suggest that such datasets will be useful to study the role of tissue-specific enhancers in human biology and disease on a genome-wide scale

Topologically associating domain boundaries are required for normal genome function

Topologically associating domain (TAD) boundaries partition the genome into distinct regulatory territories. Anecdotal evidence suggests that their disruption may interfere with normal gene expression and cause disease phenotypes1,2,3, but the overall extent to which this occurs remains unknown. Here we demonstrate that targeted deletions of TAD boundaries cause a range of disruptions to normal in vivo genome function and organismal development. We used CRISPR genome editing in mice to individually delete eight TAD boundaries (11–80 kb in size) from the genome. All deletions examined resulted in detectable molecular or organismal phenotypes, which included altered chromatin interactions or gene expression, reduced viability, and anatomical phenotypes. We observed changes in local 3D chromatin architecture in 7 of 8 (88%) cases, including the merging of TADs and altered contact frequencies within TADs adjacent to the deleted boundary. For 5 of 8 (63%) loci examined, boundary deletions were associated with increased embryonic lethality or other developmental phenotypes. For example, a TAD boundary deletion near Smad3/Smad6 caused complete embryonic lethality, while a deletion near Tbx5/Lhx5 resulted in a severe lung malformation. Our findings demonstrate the importance of TAD boundary sequences for in vivo genome function and reinforce the critical need to carefully consider the potential pathogenicity of noncoding deletions affecting TAD boundaries in clinical genetics screening.This work was supported by U.S. National Institutes of Health (NIH) grants to L.A.P. and A.V. (UM1HG009421). Research was conducted at the E.O. Lawrence Berkeley National Laboratory and performed under U.S. Department of Energy Contract DE-AC02-05CH11231, University of California (UC). Phenotyping performed by the UC Davis Mouse Biology Program (MBP) (www.mousebiology.org) was funded by an NIH administrative supplement to the KOMP2 grant, 3UM1OD023221-07S1, for phenotyping non-coding elements. Adyam Akeza was supported by the NIH Bridges to Baccalaureate Program Grant R25GM095401 via UC Berkeley. J.L.-R. is supported by the Spanish Ministerio de Ciencia e Innovacion (PID2020-113497GB-I00).Peer reviewe

Noncoding deletions reveal a gene that is critical for intestinal function.

Large-scale genome sequencing is poised to provide a substantial increase in the rate of discovery of disease-associated mutations, but the functional interpretation of such mutations remains challenging. Here we show that deletions of a sequence on human chromosome 16 that we term the intestine-critical region (ICR) cause intractable congenital diarrhoea in infants1,2. Reporter assays in transgenic mice show that the ICR contains a regulatory sequence that activates transcription during the development of the gastrointestinal system. Targeted deletion of the ICR in mice caused symptoms that recapitulated the human condition. Transcriptome analysis revealed that an unannotated open reading frame (Percc1) flanks the regulatory sequence, and the expression of this gene was lost in the developing gut of mice that lacked the ICR. Percc1-knockout mice displayed phenotypes similar to those observed upon ICR deletion in mice and patients, whereas an ICR-driven Percc1 transgene was sufficient to rescue the phenotypes found in mice that lacked the ICR. Together, our results identify a gene that is critical for intestinal function and underscore the need for targeted in vivo studies to interpret the growing number of clinical genetic findings that do not affect known protein-coding genes

Fine Tuning of Craniofacial Morphology by Distant-Acting Enhancers

The shape of the human face and skull is largely genetically determined, but the genetic drivers of craniofacial morphology remain poorly understood. Here we used a combination of epigenomic profiling, in vivo characterization of candidate enhancer sequences in transgenic mice, and targeted deletion experiments to examine the role of distant-acting enhancers in craniofacial development. We identified complex regulatory landscapes, consisting of enhancers that drive a remarkable spatial complexity of developmental expression patterns. Deletion of individual craniofacial enhancers from the mouse genome resulted in significant alterations of craniofacial shape, demonstrating their functional importance in defining face and skull morphology. These results demonstrate that enhancers play a pervasive role in mammalian craniofacial development and suggest that enhancer sequence variation contributes to human facial morphology

Ultraconservation identifies a small subset of extremely constrained developmental enhancers

While experimental studies have suggested that non-coding ultraconserved DNA elements are central nodes in the regulatory circuitry that specifies mammalian embryonic development, the possible functional relevance of their>200bp of perfect sequence conservation between human-mouse-rat remains obscure 1,2. Here we have compared the in vivo enhancer activity of a genome-wide set of 231 non-exonic sequences with ultraconserved cores to that of 206 sequences that are under equivalently severe human-rodent constraint (ultra-like), but lack perfect sequence conservation. In transgenic mouse assays, 50percent of the ultraconserved and 50percent of the ultra-like conserved elements reproducibly functioned as tissue-specific enhancers at embryonic day 11.5. In this in vivo assay, we observed that ultraconserved enhancers and constrained non-ultraconserved enhancers targeted expression to a similar spectrum of tissues with a particular enrichment in the developing central nervous system. A human genome-wide comparative screen uncovered ~;;2,600 non-coding elements that evolved under ultra-like human-rodent constraint and are similarly enriched near transcriptional regulators and developmental genes as the much smaller number of ultraconserved elements. These data indicate that ultraconserved elements possessing absolute human-rodent sequence conservation are not distinct from other non-coding elements that are under comparable purifying selection in mammals and suggest they are principal constituents of the cis-regulatory framework of mammalian development

Investigation of Genetic Variation Underlying Central Obesity amongst South Asians

The LOLIPOP study is supported by the National Institute for Health Research (NIHR) Comprehensive Biomedical Research Centre Imperial College Healthcare NHS Trust, the British Heart Foundation (SP/04/002), the Medical Research Council (G0601966,G0700931), the Wellcome Trust (084723/Z/08/Z), and the NIHR (RP-PG-0407-10371). The work was carried out in part at the NIHR/Wellcome Trust Imperial Clinical Research Facility. The Sikh Diabetes Study is supported by National Institute of Health grants KO1TW006087, funded by the Fogarty International Center, R01DK082766, funded by National Institute of Diabetes and Digestive and Kidney Diseases, and a seed grant from University of Oklahoma Health Sciences Center, Oklahoma City, USA. The Mauritius Family Study is supported by the Mauritius Ministry of Health and Quality of Life, Australian Government National Health and Medical Research Council NHMRC project grant numbers 1020285 and 1037916, the Victorian Government’s OIS Program, and partly funded by US National Institutes of Health Grant DK-25446. We thank the participants and research staff who made the study possible.South Asians are 1/4 of the world’s population and have increased susceptibility to central obesity and related cardiometabolic disease. Knowledge of genetic variants affecting risk of central obesity is largely based on genome-wide association studies of common SNPs in Europeans. To evaluate the contribution of DNA sequence variation to the higher levels of central obesity (defined as waist hip ratio adjusted for body mass index, WHR) among South Asians compared to Europeans we carried out: i) a genome-wide association analysis of >6M genetic variants in 10,318 South Asians with focused analysis of population-specific SNPs; ii) an exome-wide association analysis of ~250K SNPs in protein-coding regions in 2,637 South Asians; iii) a comparison of risk allele frequencies and effect sizes of 48 known WHR SNPs in 12,240 South Asians compared to Europeans. In genome-wide analyses, we found no novel associations between common genetic variants and WHR in South Asians at P<5x10-8; variants showing equivocal association with WHR (P<1x10-5) did not replicate at P<0.05 in an independent cohort of South Asians (N = 1,922) or in published, predominantly European meta-analysis data. In the targeted analyses of 122,391 population-specific SNPs we also found no associations with WHR in South Asians at P<0.05 after multiple testing correction. Exome-wide analyses showed no new associations between genetic variants and WHR in South Asians, either individually at P<1.5x10-6 or grouped by gene locus at P<2.5x10−6. At known WHR loci, risk allele frequencies were not higher in South Asians compared to Europeans (P = 0.77), while effect sizes were unexpectedly smaller in South Asians than Europeans (P<5.0x10-8). Our findings argue against an important contribution for population-specific or cosmopolitan genetic variants underlying the increased risk of central obesity in South Asians compared to Europeans.Yeshttp://www.plosone.org/static/editorial#pee