Mapping networks of physical interactions between genomic elements using 5C technology

Long-read ChIA-PET for base-pair-resolution mapping of haplotype-specific chromatin interactions

Chromatin interaction analysis by paired-end tag sequencing (ChIA-PET) is a robust method for capturing genome-wide chromatin interactions. Unlike other 3C-based methods, it includes a chromatin immunoprecipitation (ChIP) step that enriches for interactions mediated by specific target proteins. This unique feature allows ChIA-PET to provide the functional specificity and higher resolution needed to detect chromatin interactions, which chromosome conformation capture (3C)/Hi-C approaches have not achieved. The original ChIA-PET protocol generates short paired-end tags (2 × 20 base pairs (bp)) to detect two genomic loci that are far apart on linear chromosomes but are in spatial proximity in the folded genome. We have improved the original approach by developing long-read ChIA-PET, in which the length of the paired-end tags is increased (up to 2 × 250 bp). The longer PET reads not only improve the tag-mapping efficiency but also increase the probability of covering phased single-nucleotide polymorphisms (SNPs), which allows haplotype-specific chromatin interactions to be identified. Here, we provide the detailed protocol for long-read ChIA-PET that includes cell fixation and lysis, chromatin fragmentation by sonication, ChIP, proximity ligation with a bridge linker, Tn5 tagmentation, PCR amplification and high-throughput sequencing. For a well-trained molecular biologist, it typically takes 6 d from cell harvesting to the completion of library construction, up to a further 36 h for DNA sequencing and <20 h for processing of raw sequencing reads.

Decoding transcriptional states in cancer

Gene regulatory networks determine cellular identity. In cancer, aberrations of gene networks are caused by driver mutations that often affect transcription factors and chromatin modifiers. Nevertheless, gene transcription in cancer follows the same cis-regulatory rules as normal cells, and cancer cells have served as convenient model systems to study transcriptional regulation. Tumours often show regulatory heterogeneity, with subpopulations of cells in different transcriptional states, which has important therapeutic implications. Here, we review recent experimental and computational techniques to reverse engineer cancer gene networks using transcriptome and epigenome data. New algorithms, data integration strategies, and increasing amounts of single cell genomics data provide exciting opportunities to model dynamic regulatory states at unprecedented resolution.

The identification of hematopoietic-specific regulatory elements for WASp gene expression

Chromosome Conformation Capture (3C) technology was used to identify physical interactions between the proximal Wiskott-Aldrich Syndrome protein (WASp) promoter and its distant DNA segments in Jurkat-T cells. We found that two hematopoietic specific DNase I hypersensitive (DHS) sites (proximal DHS-A, and distal DHS-B) which had high interaction frequencies with the proximal WASp promoter indicating potential regulatory activity for these DHS sites. Subsequently, we cloned several DNA fragments around the proximal DHS-A site into a luciferase reporter vector. Interestingly, no fragments showed enhancer activity, but two fragments exhibited strong silencing activity in Jurkat-T cells. After aligning the chromatin state profiling for hematopoietic and nonhematopoietic cells using the human genome browser (UCSC), we found a 5 kb putative hematopoietic specific enhancer region located 250 kb downstream of the WAS gene. This putative enhancer region contains two hematopoietic cell specific DHS sites. Subsequently, the hematopoietic specific DHS sites enhanced luciferase expression from the proximal WASp promoter in all hematopoietic cells we tested. Finally, using a lentiviral vector stable expression system, the hematopoietic specific-enhancer(s) increased GFP reporter gene expression in hematopoietic cells, and increased WASp gene expression in WASp deficient cells. This enhancer may have the potential to be used in gene therapy for hematological diseases.

The detailed 3D multi-loop aggregate/rosette chromatin architecture and functional dynamic organization of the human and mouse genomes

BACKGROUND

The dynamic three-dimensional chromatin architecture of genomes and its co-evolutionary connection to its function-the storage, expression, and replication of genetic information-is still one of the central issues in biology. Here, we describe the much debated 3D architecture of the human and mouse genomes from the nucleosomal to the megabase pair level by a novel approach combining selective high-throughput high-resolution chromosomal interaction capture (T2C), polymer simulations, and scaling analysis of the 3D architecture and the DNA sequence.

RESULTS

The genome is compacted into a chromatin quasi-fibre with ~5 ± 1 nucleosomes/11 nm, folded into stable ~30-100 kbp loops forming stable loop aggregates/rosettes connected by similar sized linkers. Minor but significant variations in the architecture are seen between cell types and functional states. The architecture and the DNA sequence show very similar fine-structured multi-scaling behaviour confirming their co-evolution and the above.

CONCLUSIONS

This architecture, its dynamics, and accessibility, balance stability and flexibility ensuring genome integrity and variation enabling gene expression/regulation by self-organization of (in)active units already in proximity. Our results agree with the heuristics of the field and allow "architectural sequencing" at a genome mechanics level to understand the inseparable systems genomic properties.

Dynamic properties of independent chromatin domains measured by correlation spectroscopy in living cells

Background

Genome organization into subchromosomal topologically associating domains (TADs) is linked to cell-type-specific gene expression programs. However, dynamic properties of such domains remain elusive, and it is unclear how domain plasticity modulates genomic accessibility for soluble factors.

Results

Here, we combine and compare a high-resolution topology analysis of interacting chromatin loci with fluorescence correlation spectroscopy measurements of domain dynamics in single living cells. We identify topologically and dynamically independent chromatin domains of ~1 Mb in size that are best described by a loop-cluster polymer model. Hydrodynamic relaxation times and gyration radii of domains are larger for open (161 ± 15 ms, 297 ± 9 nm) than for dense chromatin (88 ± 7 ms, 243 ± 6 nm) and increase globally upon chromatin hyperacetylation or ATP depletion.

Conclusions

Based on the domain structure and dynamics measurements, we propose a loop-cluster model for chromatin domains. It suggests that the regulation of chromatin accessibility for soluble factors displays a significantly stronger dependence on factor concentration than search processes within a static network.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-016-0093-1) contains supplementary material, which is available to authorized users.

Shh and ZRS enhancer colocalisation is specific to the zone of polarising activity

Limb-specific Shh expression is regulated by the (∼1 Mb distant) ZRS enhancer. In the mouse, limb bud-restricted spatiotemporal Shh expression occurs from ∼E10 to E11.5 at the distal posterior margin and is essential for correct autopod formation. Here, we have analysed the higher-order chromatin conformation of Shh in expressing and non-expressing tissues, both by fluorescence in situ hybridisation (FISH) and by chromosome conformation capture (5C). Conventional and super-resolution light microscopy identified significantly elevated frequencies of Shh/ZRS colocalisation only in the Shh-expressing regions of the limb bud, in a conformation consistent with enhancer-promoter loop formation. However, in all tissues and at all developmental stages analysed, Shh-ZRS spatial distances were still consistently shorter than those to a neural enhancer located between Shh and ZRS in the genome. 5C identified a topologically associating domain (TAD) over the Shh/ZRS genomic region and enriched interactions between Shh and ZRS throughout E11.5 embryos. Shh/ZRS colocalisation, therefore, correlates with the spatiotemporal domain of limb bud-specific Shh expression, but close Shh and ZRS proximity in the nucleus occurs regardless of whether the gene or enhancer is active. We suggest that this constrained chromatin configuration optimises the opportunity for the active enhancer to locate and instigate the expression of Shh.

Summary: Super-resolution microscopy reveals that, during mouse limb development, enhancer-driven gene expression results in the juxtaposition of Shh and its limb bud-specific enhancer only within cells of the distal posterior limb bud.

Metagenomics and Bioinformatics in Microbial Ecology: Current Status and Beyond

Metagenomic approaches are now commonly used in microbial ecology to study microbial communities in more detail, including many strains that cannot be cultivated in the laboratory. Bioinformatic analyses make it possible to mine huge metagenomic datasets and discover general patterns that govern microbial ecosystems. However, the findings of typical metagenomic and bioinformatic analyses still do not completely describe the ecology and evolution of microbes in their environments. Most analyses still depend on straightforward sequence similarity searches against reference databases. We herein review the current state of metagenomics and bioinformatics in microbial ecology and discuss future directions for the field. New techniques will allow us to go beyond routine analyses and broaden our knowledge of microbial ecosystems. We need to enrich reference databases, promote platforms that enable meta- or comprehensive analyses of diverse metagenomic datasets, devise methods that utilize long-read sequence information, and develop more powerful bioinformatic methods to analyze data from diverse perspectives.

Decoding transcriptional enhancers: Evolving from annotation to functional interpretation

Deciphering the intricate molecular processes that orchestrate the spatial and temporal regulation of genes has become an increasingly major focus of biological research. The differential expression of genes by diverse cell types with a common genome is a hallmark of complex cellular functions, as well as the basis for multicellular life. Importantly, a more coherent understanding of gene regulation is critical for defining developmental processes, evolutionary principles and disease etiologies. Here we present our current understanding of gene regulation by focusing on the role of enhancer elements in these complex processes. Although functional genomic methods have provided considerable advances to our understanding of gene regulation, these assays, which are usually performed on a genome-wide scale, typically provide correlative observations that lack functional interpretation. Recent innovations in genome editing technologies have placed gene regulatory studies at an exciting crossroads, as systematic, functional evaluation of enhancers and other transcriptional regulatory elements can now be performed in a coordinated, high-throughput manner across the entire genome. This review provides insights on transcriptional enhancer function, their role in development and disease, and catalogues experimental tools commonly used to study these elements. Additionally, we discuss the crucial role of novel techniques in deciphering the complex gene regulatory landscape and how these studies will shape future research.

EVI1 Interferes with Myeloid Maturation via Transcriptional Repression of Cebpa, via Binding to Two Far Downstream Regulatory Elements

One mechanism by which oncoproteins work is through perturbation of cellular maturation; understanding the mechanisms by which this occurs can lead to the development of targeted therapies. EVI1 is a zinc finger oncoprotein involved in the development of acute myeloid leukemia; previous work has shown it to interfere with the maturation of granulocytes from immature precursors. Here we investigate the mechanism by which that occurs, using an immortalized hematopoietic progenitor cell line, EML-C1, as a model system. We document that overexpression of EVI1 abrogates retinoic acid-induced maturation of EML cells into committed myeloid cells, a process that can be documented by the down-regulation of stem cell antigen-1 and acquisition of responsiveness to granulocyte-macrophage colony-stimulating factor. We show that this requires DNA binding capacity of EVI1, suggesting that downstream target genes are involved. We identify the myeloid regulator Cebpa as a target gene and identify two EVI1 binding regions within evolutionarily conserved enhancer elements at +35 and +37 kb relative to the gene. EVI1 can strongly suppress Cebpa transcription, and add-back of Cebpa into EVI1-expressing EML cells partially corrects the block in maturation. We identify the DNA sequences to which EVI1 binds at +35 and +37 kb and show that mutation of one of these releases Cebpa from EVI1-induced suppression. We observe a more complex picture in primary bone marrow cells, where EVI1 suppresses Cebpa in stem cells but not in more committed progenitors. Our data thus identify a regulatory node by which EVI1 contributes to leukemia, and this represents a possible therapeutic target for treatment of EVI1-expressing leukemia.

Genomic Locations of Conserved Noncoding Sequences and Their Proximal Protein-Coding Genes in Mammalian Expression Dynamics

Experimental studies have found the involvement of certain conserved noncoding sequences (CNSs) in the regulation of the proximal protein-coding genes in mammals. However, reported cases of long range enhancer activities and inter-chromosomal regulation suggest that proximity of CNSs to protein-coding genes might not be important for regulation. To test the importance of the CNS genomic location, we extracted the CNSs conserved between chicken and four mammalian species (human, mouse, dog, and cattle). These CNSs were confirmed to be under purifying selection. The intergenic CNSs are often found in clusters in gene deserts, where protein-coding genes are in paucity. The distribution pattern, ChIP-Seq, and RNA-Seq data suggested that the CNSs are more likely to be regulatory elements and not corresponding to long intergenic noncoding RNAs. Physical distances between CNS and their nearest protein coding genes were well conserved between human and mouse genomes, and CNS-flanking genes were often found in evolutionarily conserved genomic neighborhoods. ChIP-Seq signal and gene expression patterns also suggested that CNSs regulate nearby genes. Interestingly, genes with more CNSs have more evolutionarily conserved expression than those with fewer CNSs. These computationally obtained results suggest that the genomic locations of CNSs are important for their regulatory functions. In fact, various kinds of evolutionary constraints may be acting to maintain the genomic locations of CNSs and protein-coding genes in mammals to ensure proper regulation.

Biochemical Methods To Investigate lncRNA and the Influence of lncRNA:Protein Complexes on Chromatin

Long noncoding RNAs (lncRNAs), defined as nontranslated transcripts greater than 200 nucleotides in length, are often differentially expressed throughout developmental stages, tissue types, and disease states. The identification, visualization, and suppression/overexpression of these sequences have revealed impacts on a wide range of biological processes, including epigenetic regulation. Biochemical investigations on select systems have revealed striking insight into the biological roles of lncRNAs and lncRNA:protein complexes, which in turn prompt even more unanswered questions. To begin, multiple protein- and RNA-centric technologies have been employed to isolate lncRNA:protein and lncRNA:chromatin complexes. LncRNA interactions with the multi-subunit protein complex PRC2, which acts as a transcriptional silencer, represent some of the few cases where the binding affinity, selectivity, and activity of a lncRNA:protein complex have been investigated. At the same time, recent reports of full-length lncRNA secondary structures suggest the formation of complex structures with multiple independent folding domains and pave the way for more detailed structural investigations and predictions of lncRNA three-dimensional structure. This review will provide an overview of the methods and progress made to date as well as highlight new methods that promise to further inform the molecular recognition, specificity, and function of lncRNAs.

Ensembl regulation resources

New experimental techniques in epigenomics allow researchers to assay a diversity of highly dynamic features such as histone marks, DNA modifications or chromatin structure. The study of their fluctuations should provide insights into gene expression regulation, cell differentiation and disease. The Ensembl project collects and maintains the Ensembl regulation data resources on epigenetic marks, transcription factor binding and DNA methylation for human and mouse, as well as microarray probe mappings and annotations for a variety of chordate genomes. From this data, we produce a functional annotation of the regulatory elements along the human and mouse genomes with plans to expand to other species as data becomes available. Starting from well-studied cell lines, we will progressively expand our library of measurements to a greater variety of samples. Ensembl’s regulation resources provide a central and easy-to-query repository for reference epigenomes. As with all Ensembl data, it is freely available at http://www.ensembl.org, from the Perl and REST APIs and from the public Ensembl MySQL database server at ensembldb.ensembl.org.

Database URL: http://www.ensembl.org

Extremely Long-Range Chromatin Loops Link Topological Domains to Facilitate a Diverse Antibody Repertoire

Early B cell development is characterized by large-scale Igh locus contraction prior to V(D)J recombination to facilitate a highly diverse Ig repertoire. However, an understanding of the molecular architecture that mediates locus contraction remains unclear. We have combined high-resolution chromosome conformation capture (3C) techniques with 3D DNA FISH to identify three conserved topological subdomains. Each of these topological folds encompasses a major VH gene family that become juxtaposed in pro-B cells via megabase-scale chromatin looping. The transcription factor Pax5 organizes the subdomain that spans the VHJ558 gene family. In its absence, the J558 VH genes fail to associate with the proximal VH genes, thereby providing a plausible explanation for reduced VHJ558 gene rearrangements in Pax5-deficient pro-B cells. We propose that Igh locus contraction is the cumulative effect of several independently controlled chromatin subdomains that provide the structural infrastructure to coordinate optimal antigen receptor assembly.

In Silico Functional Annotation of Genomic Variation

This unit describes the concepts and practical techniques for annotating genomic variants in the human genome to estimate their functional significance. With the rapid increase of available whole exome and whole genome sequencing information for human studies, annotation techniques have become progressively more important for highlighting and prioritizing nucleotide variants and their potential impact on genes and other genetic constructs. Here, we present an overview of different types of variant annotation approaches and elaborate on their foundations, assumptions, and the downstream consequences of their use. Computational approaches and tools to assign annotations and to identify variants are reviewed. Further, the general philosophy of assigning potential function to a genetic change within the biological context of a disease is discussed.

Analysis of long-range interactions in primary human cells identifies cooperative CFTR regulatory elements

A mechanism by which control DNA elements regulate transcription over large linear genomic distances is by achieving close physical proximity with genes, and looping of the intervening chromatin paths. Alterations of such regulatory 'chromatin looping' systems are likely to play a critical role in human genetic disease at large. Here, we studied the spatial organization of a ≈790 kb locus encompassing the cystic fibrosis transmembrane conductance regulator (CFTR) gene. Dysregulation of CFTR is responsible for cystic fibrosis, which is the most common lethal genetic disorder in Caucasian populations. CFTR is a relatively large gene of 189 kb with a rather complex tissue-specific and temporal expression profile. We used chromatin conformation at the CFTR locus to identify new DNA sequences that regulate its transcription. By comparing 5C chromatin interaction maps of the CFTR locus in expressing and non-expressing human primary cells, we identified several new contact points between the CFTR promoter and its surroundings, in addition to regions featuring previously described regulatory elements. We demonstrate that two of these novel interacting regions cooperatively increase CFTR expression, and suggest that the new enhancer elements located on either side of the gene are brought together through chromatin looping via CTCF.

Demystifying the secret mission of enhancers: linking distal regulatory elements to target genes

Enhancers are short regulatory sequences bound by sequence-specific transcription factors and play a major role in the spatiotemporal specificity of gene expression patterns in development and disease. While it is now possible to identify enhancer regions genomewide in both cultured cells and primary tissues using epigenomic approaches, it has been more challenging to develop methods to understand the function of individual enhancers because enhancers are located far from the gene(s) that they regulate. However, it is essential to identify target genes of enhancers not only so that we can understand the role of enhancers in disease but also because this information will assist in the development of future therapeutic options. After reviewing models of enhancer function, we discuss recent methods for identifying target genes of enhancers. First, we describe chromatin structure-based approaches for directly mapping interactions between enhancers and promoters. Second, we describe the use of correlation-based approaches to link enhancer state with the activity of nearby promoters and/or gene expression. Third, we describe how to test the function of specific enhancers experimentally by perturbing enhancer–target relationships using high-throughput reporter assays and genome editing. Finally, we conclude by discussing as yet unanswered questions concerning how enhancers function, how target genes can be identified, and how to distinguish direct from indirect changes in gene expression mediated by individual enhancers.

An Overview of Genome Organization and How We Got There: from FISH to Hi-C

In humans, nearly two meters of genomic material must be folded to fit inside each micrometer-scale cell nucleus while remaining accessible for gene transcription, DNA replication, and DNA repair. This fact highlights the need for mechanisms governing genome organization during any activity and to maintain the physical organization of chromosomes at all times. Insight into the functions and three-dimensional structures of genomes comes mostly from the application of visual techniques such as fluorescence in situ hybridization (FISH) and molecular approaches including chromosome conformation capture (3C) technologies. Recent developments in both types of approaches now offer the possibility of exploring the folded state of an entire genome and maybe even the identification of how complex molecular machines govern its shape. In this review, we present key methodologies used to study genome organization and discuss what they reveal about chromosome conformation as it relates to transcription regulation across genomic scales in mammals.

Boundary Associated Long Noncoding RNA Mediates Long-Range Chromosomal Interactions

CCCTC binding factor (CTCF) is involved in organizing chromosomes into mega base-sized, topologically associated domains (TADs) along with other factors that define sub-TAD organization. CTCF-Cohesin interactions have been shown to be critical for transcription insulation activity as it stabilizes long-range interactions to promote proper gene expression. Previous studies suggest that heterochromatin boundary activity of CTCF may be independent of Cohesin, and there may be additional mechanisms for defining topological domains. Here, we show that a boundary site we previously identified known as CTCF binding site 5 (CBS5) from the homeotic gene cluster A (HOXA) locus exhibits robust promoter activity. This promoter activity from the CBS5 boundary element generates a long noncoding RNA that we designate as boundary associated long noncoding RNA-1 (blncRNA1). Functional characterization of this RNA suggests that the transcript stabilizes long-range interactions at the HOXA locus and promotes proper expression of HOXA genes. Additionally, our functional analysis also shows that this RNA is not needed in the stabilization of CTCF-Cohesin interactions however CTCF-Cohesin interactions are critical in the transcription of blncRNA1. Thus, the CTCF-associated boundary element, CBS5, employs both Cohesin and noncoding RNA to establish and maintain topologically associated domains at the HOXA locus.

Prediction of Causal Candidate Genes in Coronary Artery Disease Loci

OBJECTIVE

Genome-wide association studies have to date identified 159 significant and suggestive loci for coronary artery disease (CAD). We now report comprehensive bioinformatics analyses of sequence variation in these loci to predict candidate causal genes.

APPROACH AND RESULTS

All annotated genes in the loci were evaluated with respect to protein-coding single-nucleotide polymorphism and gene expression parameters. The latter included expression quantitative trait loci, tissue specificity, and miRNA binding. High priority candidate genes were further identified based on literature searches and our experimental data. We conclude that the great majority of causal variations affecting CAD risk occur in noncoding regions, with 41% affecting gene expression robustly versus 6% leading to amino acid changes. Many of these genes differed from the traditionally annotated genes, which was usually based on proximity to the lead single-nucleotide polymorphism. Indeed, we obtained evidence that genetic variants at CAD loci affect 98 genes which had not been linked to CAD previously.

CONCLUSIONS

Our results substantially revise the list of likely candidates for CAD and suggest that genome-wide association studies efforts in other diseases may benefit from similar bioinformatics analyses.

Inferring 3D chromatin structure using a multiscale approach based on quaternions

Background

The knowledge of the spatial organisation of the chromatin fibre in cell nuclei helps researchers to understand the nuclear machinery that regulates dna activity. Recent experimental techniques of the type Chromosome Conformation Capture (3c, or similar) provide high-resolution, high-throughput data consisting in the number of times any possible pair of dna fragments is found to be in contact, in a certain population of cells. As these data carry information on the structure of the chromatin fibre, several attempts have been made to use them to obtain high-resolution 3d reconstructions of entire chromosomes, or even an entire genome. The techniques proposed treat the data in different ways, possibly exploiting physical-geometric chromatin models. One popular strategy is to transform contact data into Euclidean distances between pairs of fragments, and then solve a classical distance-to-geometry problem.

Results

We developed and tested a reconstruction technique that does not require translating contacts into distances, thus avoiding a number of related drawbacks. Also, we introduce a geometrical chromatin chain model that allows us to include sound biochemical and biological constraints in the problem. This model can be scaled at different genomic resolutions, where the structures of the coarser models are influenced by the reconstructions at finer resolutions. The search in the solution space is then performed by a classical simulated annealing, where the model is evolved efficiently through quaternion operators. The presence of appropriate constraints permits the less reliable data to be overlooked, so the result is a set of plausible chromatin configurations compatible with both the data and the prior knowledge.

Conclusions

To test our method, we obtained a number of 3d chromatin configurations from hi-c data available in the literature for the long arm of human chromosome 1, and validated their features against known properties of gene density and transcriptional activity. Our results are compatible with biological features not introduced a priori in the problem: structurally different regions in our reconstructions highly correlate with functionally different regions as known from literature and genomic repositories.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0667-0) contains supplementary material, which is available to authorized users.

The ensembl regulatory build

Most genomic variants associated with phenotypic traits or disease do not fall within gene coding regions, but in regulatory regions, rendering their interpretation difficult. We collected public data on epigenetic marks and transcription factor binding in human cell types and used it to construct an intuitive summary of regulatory regions in the human genome. We verified it against independent assays for sensitivity. The Ensembl Regulatory Build will be progressively enriched when more data is made available. It is freely available on the Ensembl browser, from the Ensembl Regulation MySQL database server and in a dedicated track hub.

Spatial genome organization: contrasting views from chromosome conformation capture and fluorescence in situ hybridization

Although important for gene regulation, most studies of genome organization use either fluorescence in situ hybridization (FISH) or chromosome conformation capture (3C) methods. FISH directly visualizes the spatial relationship of sequences but is usually applied to a few loci at a time. The frequency at which sequences are ligated together by formaldehyde cross-linking can be measured genome-wide by 3C methods, with higher frequencies thought to reflect shorter distances. FISH and 3C should therefore give the same views of genome organization, but this has not been tested extensively. We investigated the murine HoxD locus with 3C carbon copy (5C) and FISH in different developmental and activity states and in the presence or absence of epigenetic regulators. We identified situations in which the two data sets are concordant but found other conditions under which chromatin topographies extrapolated from 5C or FISH data are not compatible. We suggest that products captured by 3C do not always reflect spatial proximity, with ligation occurring between sequences located hundreds of nanometers apart, influenced by nuclear environment and chromatin composition. We conclude that results obtained at high resolution with either 3C methods or FISH alone must be interpreted with caution and that views about genome organization should be validated by independent methods.

Chromatin associations in Arabidopsis interphase nuclei

The arrangement of chromatin within interphase nuclei seems to be caused by topological constraints and related to gene expression depending on tissue and developmental stage. In yeast and animals it was found that homologous and heterologous chromatin association are required to realize faithful expression and DNA repair. To test whether such associations are present in plants we analyzed Arabidopsis thaliana interphase nuclei by FISH using probes from different chromosomes. We found that chromatin fiber movement and variable associations, although in general relatively seldom, may occur between euchromatin segments along chromosomes, sometimes even over large distances. The combination of euchromatin segments bearing high or low co-expressing genes did not reveal different association frequencies probably due to adjacent genes of deviating expression patterns. Based on previous data and on FISH analyses presented here, we conclude that the global interphase chromatin organization in A. thaliana is relatively stable, due to the location of its 10 centromeres at the nuclear periphery and of the telomeres mainly at the centrally localized nucleolus. Nevertheless, chromatin movement enables a flexible spatial genome arrangement in plant nuclei.

Genome-wide analysis of local chromatin packing in Arabidopsis thaliana

The spatial arrangement of interphase chromosomes in the nucleus is important for gene expression and genome function in animals and in plants. The recently developed Hi-C technology is an efficacious method to investigate genome packing. Here we present a detailed Hi-C map of the three-dimensional genome organization of the plant Arabidopsis thaliana. We find that local chromatin packing differs from the patterns seen in animals, with kilobasepair-sized segments that have much higher intrachromosome interaction rates than neighboring regions, representing a dominant local structural feature of genome conformation in A. thaliana. These regions, which appear as positive strips on two-dimensional representations of chromatin interaction, are enriched in epigenetic marks H3K27me3, H3.1, and H3.3. We also identify more than 400 insulator-like regions. Furthermore, although topologically associating domains (TADs), which are prominent in animals, are not an obvious feature of A. thaliana genome packing, we found more than 1000 regions that have properties of TAD boundaries, and a similar number of regions analogous to the interior of TADs. The insulator-like, TAD-boundary-like, and TAD-interior-like regions are each enriched for distinct epigenetic marks and are each correlated with different gene expression levels. We conclude that epigenetic modifications, gene density, and transcriptional activity combine to shape the local packing of the A. thaliana nuclear genome.

The Runx-PU.1 pathway preserves normal and AML/ETO9a leukemic stem cells

Runx transcription factors contribute to hematopoiesis and are frequently implicated in hematologic malignancies. All three Runx isoforms are expressed at the earliest stages of hematopoiesis; however, their function in hematopoietic stem cells (HSCs) is not fully elucidated. Here, we show that Runx factors are essential in HSCs by driving the expression of the hematopoietic transcription factor PU.1. Mechanistically, by using a knockin mouse model in which all three Runx binding sites in the -14kb enhancer of PU.1 are disrupted, we observed failure to form chromosomal interactions between the PU.1 enhancer and its proximal promoter. Consequently, decreased PU.1 levels resulted in diminished long-term HSC function through HSC exhaustion, which could be rescued by reintroducing a PU.1 transgene. Similarly, in a mouse model of AML/ETO9a leukemia, disrupting the Runx binding sites resulted in decreased PU.1 levels. Leukemia onset was delayed, and limiting dilution transplantation experiments demonstrated functional loss of leukemia-initiating cells. This is surprising, because low PU.1 levels have been considered a hallmark of AML/ETO leukemia, as indicated in mouse models and as shown here in samples from leukemic patients. Our data demonstrate that Runx-dependent PU.1 chromatin interaction and transcription of PU.1 are essential for both normal and leukemia stem cells.

Targeted Chromatin Capture (T2C): a novel high resolution high throughput method to detect genomic interactions and regulatory elements

Background

Significant efforts have recently been put into the investigation of the spatial organization and the chromatin-interaction networks of genomes. Chromosome conformation capture (3C) technology and its derivatives are important tools used in this effort. However, many of these have limitations, such as being limited to one viewpoint, expensive with moderate to low resolution, and/or requiring a large sequencing effort. Techniques like Hi-C provide a genome-wide analysis. However, it requires massive sequencing effort with considerable costs. Here we describe a new technique termed Targeted Chromatin Capture (T2C), to interrogate large selected regions of the genome. T2C provides an unbiased view of the spatial organization of selected loci at superior resolution (single restriction fragment resolution, from 2 to 6 kbp) at much lower costs than Hi-C due to the lower sequencing effort.

Results

We applied T2C on well-known model regions, the mouse β-globin locus and the human H19/IGF2 locus. In both cases we identified all known chromatin interactions. Furthermore, we compared the human H19/IGF2 locus data obtained from different chromatin conformation capturing methods with T2C data. We observed the same compartmentalization of the locus, but at a much higher resolution (single restriction fragments vs. the common 40 kbp bins) and higher coverage. Moreover, we compared the β-globin locus in two different biological samples (mouse primary erythroid cells and mouse fetal brain), where it is either actively transcribed or not, to identify possible transcriptional dependent interactions. We identified the known interactions in the β-globin locus and the same topological domains in both mouse primary erythroid cells and in mouse fetal brain with the latter having fewer interactions probably due to the inactivity of the locus. Furthermore, we show that interactions due to the important chromatin proteins, Ldb1 and Ctcf, in both tissues can be analyzed easily to reveal their role on transcriptional interactions and genome folding.

Conclusions

T2C is an efficient, easy, and affordable with high (restriction fragment) resolution tool to address both genome compartmentalization and chromatin-interaction networks for specific genomic regions at high resolution for both clinical and non-clinical research.

The genome in three dimensions: a new frontier in human brain research

Less than 1.5% of the human genome encodes protein. However, vast portions of the human genome are subject to transcriptional and epigenetic regulation, and many noncoding regulatory DNA elements are thought to regulate the spatial organization of interphase chromosomes. For example, chromosomal "loopings" are pivotal for the orderly process of gene expression, by enabling distal regulatory enhancer or silencer elements to directly interact with proximal promoter and transcription start sites, potentially bypassing hundreds of kilobases of interspersed sequence on the linear genome. To date, however, epigenetic studies in the human brain are mostly limited to the exploration of DNA methylation and posttranslational modifications of the nucleosome core histones. In contrast, very little is known about the regulation of supranucleosomal structures. Here, we show that chromosome conformation capture, a widely used approach to study higher-order chromatin, is applicable to tissue collected postmortem, thereby informing about genome organization in the human brain. We introduce chromosome conformation capture protocols for brain and compare higher-order chromatin structures at the chromosome 6p22.2-22.1 schizophrenia and bipolar disorder susceptibility locus, and additional neurodevelopmental risk genes, (DPP10, MCPH1) in adult prefrontal cortex and various cell culture systems, including neurons derived from reprogrammed skin cells. We predict that the exploration of three-dimensional genome architectures and function will open up new frontiers in human brain research and psychiatric genetics and provide novel insights into the epigenetic risk architectures of regulatory noncoding DNA.

3D chromosome modeling with semi-definite programming and Hi-C data

For a long period of time, scientists studied genomes while assuming they are linear. Recently, chromosome conformation capture (3C)-based technologies, such as Hi-C, have been developed that provide the loci contact frequencies among loci pairs in a genome-wide scale. The technology unveiled that two far-apart loci can interact in the tested genome. It indicated that the tested genome forms a three-dimensional (3D) chromosomal structure within the nucleus. With the available Hi-C data, our next challenge is to model the 3D chromosomal structure from the 3C-derived data computationally. This article presents a deterministic method called ChromSDE, which applies semi-definite programming techniques to find the best structure fitting the observed data and uses golden section search to find the correct parameter for converting the contact frequency to spatial distance. Further, we develop a measure called consensus index to indicate if the Hi-C data corresponds to a single structure or a mixture of structures. To the best of our knowledge, ChromSDE is the only method that can guarantee recovering the correct structure in the noise-free case. In addition, we prove that the parameter of conversion from contact frequency to spatial distance will change under different resolutions theoretically and empirically. Using simulation data and real Hi-C data, we showed that ChromSDE is much more accurate and robust than existing methods. Finally, we demonstrated that interesting biological findings can be uncovered from our predicted 3D structure.

Human colorectal cancer-specific CCAT1-L lncRNA regulates long-range chromatin interactions at the MYC locus

The human 8q24 gene desert contains multiple enhancers that form tissue-specific long-range chromatin loops with the MYC oncogene, but how chromatin looping at the MYC locus is regulated remains poorly understood. Here we demonstrate that a long noncoding RNA (lncRNA), CCAT1-L, is transcribed specifically in human colorectal cancers from a locus 515 kb upstream of MYC. This lncRNA plays a role in MYC transcriptional regulation and promotes long-range chromatin looping. Importantly, the CCAT1-L locus is located within a strong super-enhancer and is spatially close to MYC. Knockdown of CCAT1-L reduced long-range interactions between the MYC promoter and its enhancers. In addition, CCAT1-L interacts with CTCF and modulates chromatin conformation at these loop regions. These results reveal an important role of a previously unannotated lncRNA in gene regulation at the MYC locus.

Hox in motion: tracking HoxA cluster conformation during differentiation

Three-dimensional genome organization is an important higher order transcription regulation mechanism that can be studied with the chromosome conformation capture techniques. Here, we combined chromatin organization analysis by chromosome conformation capture-carbon copy, computational modeling and epigenomics to achieve the first integrated view, through time, of a connection between chromatin state and its architecture. We used this approach to examine the chromatin dynamics of the HoxA cluster in a human myeloid leukemia cell line at various stages of differentiation. We found that cellular differentiation involves a transient activation of the 5'-end HoxA genes coinciding with a loss of contacts throughout the cluster, and by specific silencing at the 3'-end with H3K27 methylation. The 3D modeling of the data revealed an extensive reorganization of the cluster between the two previously reported topologically associated domains in differentiated cells. Our results support a model whereby silencing by polycomb group proteins and reconfiguration of CTCF interactions at a topologically associated domain boundary participate in changing the HoxA cluster topology, which compartmentalizes the genes following differentiation.

The grammar of transcriptional regulation

Eukaryotes employ combinatorial strategies to generate a variety of expression patterns from a relatively small set of regulatory DNA elements. As in any other language, deciphering the mapping between DNA and expression requires an understanding of the set of rules that govern basic principles in transcriptional regulation, the functional elements involved, and the ways in which they combine to orchestrate a transcriptional output. Here, we review the current understanding of various grammatical rules, including the effect on expression of the number of transcription factor binding sites, their location, orientation, affinity and activity; co-association with different factors; and intrinsic nucleosome organization. We review different methods that are used to study the grammar of transcription regulation, highlight gaps in current understanding, and discuss how recent technological advances may be utilized to bridge them.

Detection of short-range chromatin interactions by chromosome conformation capture (3C) in yeast

We describe a modified 3C ("chromosome conformation capture") protocol for detection of transient, short-range chromatin interactions in the yeast Saccharomyces cerevisiae. 3C was initially described by Job Dekker and involves formaldehyde cross-linking to stabilize transient chromatin interactions, followed by restriction digestion, ligation, and locus-specific PCR. As such, 3C reveals complex three-dimensional interactions between distal genetic elements within intact cells at high resolution. Using a modified version of Dekker's protocol, we are able to detect gene loops that juxtapose promoter and terminator regions of yeast genes with ORFs as short as 1 kb. We are using this technique to define the cis- and trans-acting requirements for the formation and maintenance of gene loops, and to elucidate their physiological consequences. We anticipate that this method will be generally applicable to detect dynamic, short-range chromatin interactions, not limited to gene loops.

Clustering of tissue-specific sub-TADs accompanies the regulation of HoxA genes in developing limbs

HoxA genes exhibit central roles during development and causal mutations have been found in several human syndromes including limb malformation. Despite their importance, information on how these genes are regulated is lacking. Here, we report on the first identification of bona fide transcriptional enhancers controlling HoxA genes in developing limbs and show that these enhancers are grouped into distinct topological domains at the sub-megabase scale (sub-TADs). We provide evidence that target genes and regulatory elements physically interact with each other through contacts between sub-TADs rather than by the formation of discreet “DNA loops”. Interestingly, there is no obvious relationship between the functional domains of the enhancers within the limb and how they are partitioned among the topological domains, suggesting that sub-TAD formation does not rely on enhancer activity. Moreover, we show that suppressing the transcriptional activity of enhancers does not abrogate their contacts with HoxA genes. Based on these data, we propose a model whereby chromatin architecture defines the functional landscapes of enhancers. From an evolutionary standpoint, our data points to the convergent evolution of HoxA and HoxD regulation in the fin-to-limb transition, one of the major morphological innovations in vertebrates.

Hox genes encode transcription factors with crucial roles during development. These genes are grouped in four different clusters names HoxA, B, C, and D. Mutations in genes of the HoxA and D clusters have been found in several human syndromes, affecting in some cases limb development. Despite their essential role and contrary to the genes of the HoxD cluster, little is known about how the HoxA genes are regulated. Here, we identified a large set of regulatory elements controlling HoxA genes during limb development. By studying spatial chromatin organization at the HoxA region, we found that the regulatory elements are spatially clustered regardless of their activity. Clustering of enhancers define tissue-specific chromatin domains that interact specifically with each other and with active genes in the limb. Our findings give support to the emerging concept that chromatin architecture defines the functional properties of genomes. Additionally, our study suggests a common constraint of the chromatin topology in the evolution of HoxA and HoxD regulation in the emergence of the hand/foot, which is one of the major morphological innovations in vertebrates.

Chromosome conformation capture-on-chip analysis of long-range cis-interactions of the SOX9 promoter

Evolutionarily conserved transcription factor SOX9 is essential for the differentiation of chondrocytes and the development of testes. Heterozygous point mutations and genomic deletions involving SOX9 lead to campomelic dysplasia (CD), a skeletal malformation syndrome often associated with sex reversal. Chromosomal rearrangements with breakpoints mapping up to 1.6 Mb up- and downstream to SOX9, and likely disrupting its distant cis-regulatory elements, have been described in patients with milder forms of CD. Based on the location of these aberration breakpoints, four clusters upstream of SOX9 have been defined. Interestingly, we found that each of these intervals overlaps a gene encoding long noncoding RNA (lncRNA), suggesting that lncRNAs may contribute to long-range regulation of SOX9 expression. One of the four upstream regions, RevSex (517-595 kb 5' to SOX9), is associated with sex reversal, and was suggested to harbor a testis-specific and sex-determining enhancer. Another sex-determining interval was mapped to a gene desert >1.3 Mb downstream of SOX9. We have performed chromosome conformation capture-on-chip (4C) analysis in Sertoli cells and lymphoblasts to verify the proposed long-range interactions of the SOX9 promoter and to identify potential novel regulatory elements that might be responsible for sex reversal in patients with CD. We identified several novel potentially cis-interacting regions both up- and downstream to SOX9, with some of them overlapping lncRNA genes. Our data point to lncRNAs as likely mediators of some of these regulatory interactions.

The interactomes of POU5F1 and SOX2 enhancers in human embryonic stem cells

The genes POU5F1 and SOX2 are critical for pluripotency and reprogramming, yet the chromosomal organization around these genes remains poorly understood. We assayed long-range chromosomal interactions on putative enhancers of POU5F1 and SOX2 genes in human embryonic stem cells (hESCs) using 4C-Seq technique. We discovered that their frequent interacting regions mainly overlap with early DNA replication domains. The interactomes are associated with active histone marks and enriched with 5-hydroxymethylcytosine sites. In hESCs, genes within the interactomes have elevated expression. Additionally, some genes associated with the POU5F1 enhancer contribute to pluripotency. Binding sites for multiple DNA binding proteins, including ATF3, CTCF, GABPA, JUND, NANOG, RAD21 and YY1, are enriched in both interactomes. The RARG locus, frequently interacting with the POU5F1 locus, has abundant RAD21 binding sites co-localized with other protein binding sites. Thus the interactomes of these two pluripotency genes could be an important part of the regulatory network in hESCs.

The statistical-mechanics of chromosome conformation capture

Since Jacob and Monod’s characterization of the role of DNA elements in gene control, it has been recognized that the linear organization of genome structure is important for the regulation of gene transcription and hence the manifestation of phenotypes. Similarly, it has long been hypothesized that the spatial organization (in three dimensions evolving through time), as part of the epigenome, makes a significant contribution to the genotype-phenotype transition. Proximity ligation assays commonly known as chromosome conformation capture (3C) and 3C based methodologies (e.g., GCC, HiC, and ChIA-Pet) are increasingly being incorporated into empirical studies to investigate the role that three-dimensional genome structure plays in the regulation of phenotype. The apparent simplicity of these methodologies—crosslink chromatin, digest, dilute, ligate, detect interactions—belies the complexity of the data and the considerations that should be taken into account to ensure the generation and accurate interpretation of reliable data. Here we discuss the probabilistic nature of these methodologies and how this contributes to their endogenous limitations.

Role of CTCF protein in regulating FMR1 locus transcription

Fragile X syndrome (FXS), the leading cause of inherited intellectual disability, is caused by epigenetic silencing of the FMR1 gene, through expansion and methylation of a CGG triplet repeat (methylated full mutation). An antisense transcript (FMR1-AS1), starting from both promoter and intron 2 of the FMR1 gene, was demonstrated in transcriptionally active alleles, but not in silent FXS alleles. Moreover, a DNA methylation boundary, which is lost in FXS, was recently identified upstream of the FMR1 gene. Several nuclear proteins bind to this region, like the insulator protein CTCF. Here we demonstrate for the first time that rare unmethylated full mutation (UFM) alleles present the same boundary described in wild type (WT) alleles and that CTCF binds to this region, as well as to the FMR1 gene promoter, exon 1 and intron 2 binding sites. Contrariwise, DNA methylation prevents CTCF binding to FXS alleles. Drug-induced CpGs demethylation does not restore this binding. CTCF knock-down experiments clearly established that CTCF does not act as insulator at the active FMR1 locus, despite the presence of a CGG expansion. CTCF depletion induces heterochromatinic histone configuration of the FMR1 locus and results in reduction of FMR1 transcription, which however is not accompanied by spreading of DNA methylation towards the FMR1 promoter. CTCF depletion is also associated with FMR1-AS1 mRNA reduction. Antisense RNA, like sense transcript, is upregulated in UFM and absent in FXS cells and its splicing is correlated to that of the FMR1-mRNA. We conclude that CTCF has a complex role in regulating FMR1 expression, probably through the organization of chromatin loops between sense/antisense transcriptional regulatory regions, as suggested by bioinformatics analysis.

Understanding spatial organizations of chromosomes via statistical analysis of Hi-C data

Understanding how chromosomes fold provides insights into the transcription regulation, hence, the functional state of the cell. Using the next generation sequencing technology, the recently developed Hi-C approach enables a global view of spatial chromatin organization in the nucleus, which substantially expands our knowledge about genome organization and function. However, due to multiple layers of biases, noises and uncertainties buried in the protocol of Hi-C experiments, analyzing and interpreting Hi-C data poses great challenges, and requires novel statistical methods to be developed. This article provides an overview of recent Hi-C studies and their impacts on biomedical research, describes major challenges in statistical analysis of Hi-C data, and discusses some perspectives for future research.

Architectural protein subclasses shape 3D organization of genomes during lineage commitment

Understanding the topological configurations of chromatin may reveal valuable insights into how the genome and epigenome act in concert to control cell fate during development. Here, we generate high-resolution architecture maps across seven genomic loci in embryonic stem cells and neural progenitor cells. We observe a hierarchy of 3D interactions that undergo marked reorganization at the submegabase scale during differentiation. Distinct combinations of CCCTC-binding factor (CTCF), Mediator, and cohesin show widespread enrichment in chromatin interactions at different length scales. CTCF/cohesin anchor long-range constitutive interactions that might form the topological basis for invariant subdomains. Conversely, Mediator/cohesin bridge short-range enhancer-promoter interactions within and between larger subdomains. Knockdown of Smc1 or Med12 in embryonic stem cells results in disruption of spatial architecture and downregulation of genes found in cohesin-mediated interactions. We conclude that cell-type-specific chromatin organization occurs at the submegabase scale and that architectural proteins shape the genome in hierarchical length scales.

Multiplexed chromosome conformation capture sequencing for rapid genome-scale high-resolution detection of long-range chromatin interactions

Chromosome conformation capture (3C) technology is a powerful and increasingly popular tool for analyzing the spatial organization of genomes. Several 3C variants have been developed (e.g., 4C, 5C, ChIA-PET, Hi-C), allowing large-scale mapping of long-range genomic interactions. Here we describe multiplexed 3C sequencing (3C-seq), a 4C variant coupled to next-generation sequencing, allowing genome-scale detection of long-range interactions with candidate regions. Compared with several other available techniques, 3C-seq offers a superior resolution (typically single restriction fragment resolution; approximately 1-8 kb on average) and can be applied in a semi-high-throughput fashion. It allows the assessment of long-range interactions of up to 192 genes or regions of interest in parallel by multiplexing library sequencing. This renders multiplexed 3C-seq an inexpensive, quick (total hands-on time of 2 weeks) and efficient method that is ideal for the in-depth analysis of complex genetic loci. The preparation of multiplexed 3C-seq libraries can be performed by any investigator with basic skills in molecular biology techniques. Data analysis requires basic expertise in bioinformatics and in Linux and Python environments. The protocol describes all materials, critical steps and bioinformatics tools required for successful application of 3C-seq technology.

Sustained PU.1 levels balance cell-cycle regulators to prevent exhaustion of adult hematopoietic stem cells

To provide a lifelong supply of blood cells, hematopoietic stem cells (HSCs) need to carefully balance both self-renewing cell divisions and quiescence. Although several regulators that control this mechanism have been identified, we demonstrate that the transcription factor PU.1 acts upstream of these regulators. So far, attempts to uncover PU.1's role in HSC biology have failed because of the technical limitations of complete loss-of-function models. With the use of hypomorphic mice with decreased PU.1 levels specifically in phenotypic HSCs, we found reduced HSC long-term repopulation potential that could be rescued completely by restoring PU.1 levels. PU.1 prevented excessive HSC division and exhaustion by controlling the transcription of multiple cell-cycle regulators. Levels of PU.1 were sustained through autoregulatory PU.1 binding to an upstream enhancer that formed an active looped chromosome architecture in HSCs. These results establish that PU.1 mediates chromosome looping and functions as a master regulator of HSC proliferation.

High-throughput sequencing for biology and medicine

Advances in genome sequencing have progressed at a rapid pace, with increased throughput accompanied by plunging costs. But these advances go far beyond faster and cheaper. High-throughput sequencing technologies are now routinely being applied to a wide range of important topics in biology and medicine, often allowing researchers to address important biological questions that were not possible before. In this review, we discuss these innovative new approaches-including ever finer analyses of transcriptome dynamics, genome structure and genomic variation-and provide an overview of the new insights into complex biological systems catalyzed by these technologies. We also assess the impact of genotyping, genome sequencing and personal omics profiling on medical applications, including diagnosis and disease monitoring. Finally, we review recent developments in single-cell sequencing, and conclude with a discussion of possible future advances and obstacles for sequencing in biology and health.

ENCODE data in the UCSC Genome Browser: year 5 update

The Encyclopedia of DNA Elements (ENCODE), http://encodeproject.org, has completed its fifth year of scientific collaboration to create a comprehensive catalog of functional elements in the human genome, and its third year of investigations in the mouse genome. Since the last report in this journal, the ENCODE human data repertoire has grown by 898 new experiments (totaling 2886), accompanied by a major integrative analysis. In the mouse genome, results from 404 new experiments became available this year, increasing the total to 583, collected during the course of the project. The University of California, Santa Cruz, makes this data available on the public Genome Browser http://genome.ucsc.edu for visual browsing and data mining. Download of raw and processed data files are all supported. The ENCODE portal provides specialized tools and information about the ENCODE data sets.

From cells to chromatin: capturing snapshots of genome organization with 5C technology

In eukaryotes, genome organization can be observed on many levels and at different scales. This organization is important not only to reduce chromosome length but also for the proper execution of various biological processes. High-resolution mapping of spatial chromatin structure was made possible by the development of the chromosome conformation capture (3C) technique. 3C uses chemical cross-linking followed by proximity-based ligation of fragmented DNA to capture frequently interacting chromatin segments in cell populations. Several 3C-related methods capable of higher chromosome conformation mapping throughput were reported afterwards. These techniques include the 3C-carbon copy (5C) approach, which offers the advantage of being highly quantitative and reproducible. We provide here an updated reference protocol for the production of 5C libraries analyzed by next-generation sequencing or onto microarrays. A procedure used to verify that 3C library templates bear the high quality required to produce superior 5C libraries is also described. We believe that this detailed protocol will help guide researchers in probing spatial genome organization and its role in various biological processes.

Small noncoding differentially methylated copy-number variants, including lncRNA genes, cause a lethal lung developmental disorder

An unanticipated and tremendous amount of the noncoding sequence of the human genome is transcribed. Long noncoding RNAs (lncRNAs) constitute a significant fraction of non-protein-coding transcripts; however, their functions remain enigmatic. We demonstrate that deletions of a small noncoding differentially methylated region at 16q24.1, including lncRNA genes, cause a lethal lung developmental disorder, alveolar capillary dysplasia with misalignment of pulmonary veins (ACD/MPV), with parent-of-origin effects. We identify overlapping deletions 250 kb upstream of FOXF1 in nine patients with ACD/MPV that arose de novo specifically on the maternally inherited chromosome and delete lung-specific lncRNA genes. These deletions define a distant cis-regulatory region that harbors, besides lncRNA genes, also a differentially methylated CpG island, binds GLI2 depending on the methylation status of this CpG island, and physically interacts with and up-regulates the FOXF1 promoter. We suggest that lung-transcribed 16q24.1 lncRNAs may contribute to long-range regulation of FOXF1 by GLI2 and other transcription factors. Perturbation of lncRNA-mediated chromatin interactions may, in general, be responsible for position effect phenomena and potentially cause many disorders of human development.

Cis-acting transcriptional repression establishes a sharp boundary in chordate embryos

The function of bone morphogenetic protein (BMP) signaling in dorsoventral (DV) patterning of animal embryos is conserved among Bilateria. In vertebrates, the BMP ligand antidorsalizing morphogenetic protein (Admp) is expressed dorsally and moves to the opposite side to specify the ventral fate. Here, we show that Pinhead is an antagonist specific for Admp with a role in establishing the DV axis of the trunk epidermis in embryos of the ascidian Ciona intestinalis. Pinhead and Admp exist in tandem in the genomes of various animals from arthropods to vertebrates. This genomic configuration is important for mutually exclusive expression of these genes, because Pinhead transcription directly disturbs the action of the Admp enhancer. Our data suggest that this dual negative regulatory mechanism is widely conserved in animals.