LIGHT signaling through LTßR and HVEM in keratinocytes promotes psoriasis and atopic dermatitis-like skin inflammation.

Psoriasis (PS) and atopic dermatitis (AD) are common skin inflammatory diseases characterized by hyper-responsive keratinocytes. Although, some cytokines have been suggested to be specific for each disease, other cytokines might be central to both diseases. Here, we show that Tumor necrosis factor superfamily member 14 (TNFSF14), known as LIGHT, is required for experimental PS, similar to its requirement in experimental AD. Mice devoid of LIGHT, or deletion of either of its receptors, lymphotoxin ß receptor (LTßR) and herpesvirus entry mediator (HVEM), in keratinocytes, were protected from developing imiquimod-induced psoriatic features, including epidermal thickening and hyperplasia, and expression of PS-related genes. Correspondingly, in single cell RNA-seq analysis of PS patient biopsies, LTßR transcripts were found strongly expressed with HVEM in keratinocytes, and LIGHT was upregulated in T cells. Similar transcript expression profiles were also seen in AD biopsies, and LTßR deletion in keratinocytes also protected mice from allergen-induced AD features. Moreover, in vitro, LIGHT upregulated a broad spectrum of genes in human keratinocytes that are clinical features of both PS and AD skin lesions. Our data suggest that agents blocking LIGHT activity might be useful for therapeutic intervention in PS as well as in AD.

Rewiring of the 3D genome during acquisition of carboplatin resistance in a triple-negative breast cancer patient-derived xenograft.

Changes in the three-dimensional (3D) structure of the genome are an emerging hallmark of cancer. Cancer-associated copy number variants and single nucleotide polymorphisms promote rewiring of chromatin loops, disruption of topologically associating domains (TADs), active/inactive chromatin state switching, leading to oncogene expression and silencing of tumor suppressors. However, little is known about 3D changes during cancer progression to a chemotherapy-resistant state. We integrated chromatin conformation capture (Hi-C), RNA-seq, and whole-genome sequencing obtained from triple-negative breast cancer patient-derived xenograft primary tumors (UCD52) and carboplatin-resistant samples and found increased short-range (< 2 Mb) interactions, chromatin looping, formation of TAD, chromatin state switching into a more active state, and amplification of ATP-binding cassette transporters. Transcriptome changes suggested the role of long-noncoding RNAs in carboplatin resistance. Rewiring of the 3D genome was associated with TP53, TP63, BATF, FOS-JUN family of transcription factors and led to activation of aggressiveness-, metastasis- and other cancer-related pathways. Integrative analysis highlighted increased ribosome biogenesis and oxidative phosphorylation, suggesting the role of mitochondrial energy metabolism. Our results suggest that 3D genome remodeling may be a key mechanism underlying carboplatin resistance.

3D genome mapping identifies subgroup-specific chromosome conformations and tumor-dependency genes in ependymoma.

Ependymoma is a tumor of the brain or spinal cord. The two most common and aggressive molecular groups of ependymoma are the supratentorial ZFTA-fusion associated and the posterior fossa ependymoma group A. In both groups, tumors occur mainly in young children and frequently recur after treatment. Although molecular mechanisms underlying these diseases have recently been uncovered, they remain difficult to target and innovative therapeutic approaches are urgently needed. Here, we use genome-wide chromosome conformation capture (Hi-C), complemented with CTCF and H3K27ac ChIP-seq, as well as gene expression and DNA methylation analysis in primary and relapsed ependymoma tumors, to identify chromosomal conformations and regulatory mechanisms associated with aberrant gene expression. In particular, we observe the formation of new topologically associating domains ('neo-TADs') caused by structural variants, group-specific 3D chromatin loops, and the replacement of CTCF insulators by DNA hyper-methylation. Through inhibition experiments, we validate that genes implicated by these 3D genome conformations are essential for the survival of patient-derived ependymoma models in a group-specific manner. Thus, this study extends our ability to reveal tumor-dependency genes by 3D genome conformations even in tumors that lack targetable genetic alterations.

Lymphotoxin beta receptor signaling directly controls airway smooth muscle deregulation and asthmatic lung dysfunction.

BACKGROUND: Dysregulation of airway smooth muscle cells (ASM) is central to the severity of asthma. Which molecules dominantly control ASM in asthma is unclear. High levels of the cytokine LIGHT (aka TNFSF14) have been linked to asthma severity and lower baseline predicted FEV1 percentage, implying that signals through its receptors might directly control ASM dysfunction. OBJECTIVE: Our study sought to determine whether signaling via lymphotoxin beta receptor (LTßR) or herpesvirus entry mediator from LIGHT dominantly drives ASM hyperreactivity induced by allergen. METHODS: Conditional knockout mice deficient for LTßR or herpesvirus entry mediator in smooth muscle cells were used to determine their role in ASM deregulation and airway hyperresponsiveness (AHR) in vivo. Human ASM were used to study signals induced by LTßR. RESULTS: LTßR was strongly expressed in ASM from normal and asthmatic subjects compared to several other receptors implicated in smooth muscle deregulation. Correspondingly, conditional deletion of LTßR only in smooth muscle cells in smMHCCreLTßRfl/fl mice minimized changes in their numbers and mass as well as AHR induced by house dust mite allergen in a model of severe asthma. Intratracheal LIGHT administration independently induced ASM hypertrophy and AHR in vivo dependent on direct LTßR signals to ASM. LIGHT promoted contractility, hypertrophy, and hyperplasia of human ASM in vitro. Distinguishing LTßR from the receptors for IL-13, TNF, and IL-17, which have also been implicated in smooth muscle dysregulation, LIGHT promoted NF-κB-inducing kinase-dependent noncanonical nuclear factor kappa-light-chain enhancer of activated B cells in ASM in vitro, leading to sustained accumulation of F-actin, phosphorylation of myosin light chain kinase, and contractile activity. CONCLUSIONS: LTßR signals directly and dominantly drive airway smooth muscle hyperresponsiveness relevant for pathogenesis of airway remodeling in severe asthma.

Intermittent PI3Kδ inhibition sustains anti-tumour immunity and curbs irAEs.

Phosphoinositide 3-kinase δ (PI3Kδ) has a key role in lymphocytes, and inhibitors that target this PI3K have been approved for treatment of B cell malignancies1-3. Although studies in mouse models of solid tumours have demonstrated that PI3Kδ inhibitors (PI3Kδi) can induce anti-tumour immunity4,5, its effect on solid tumours in humans remains unclear. Here we assessed the effects of the PI3Kδi AMG319 in human patients with head and neck cancer in a neoadjuvant, double-blind, placebo-controlled randomized phase II trial (EudraCT no. 2014-004388-20). PI3Kδ inhibition decreased the number of tumour-infiltrating regulatory T (Treg) cells and enhanced the cytotoxic potential of tumour-infiltrating T cells. At the tested doses of AMG319, immune-related adverse events (irAEs) required treatment to be discontinued in 12 out of 21 of patients treated with AMG319, suggestive of systemic effects on Treg cells. Accordingly, in mouse models, PI3Kδi decreased the number of Treg cells systemically and caused colitis. Single-cell RNA-sequencing analysis revealed a PI3Kδi-driven loss of tissue-resident colonic ST2 Treg cells, accompanied by expansion of pathogenic T helper 17 (TH17) and type 17 CD8+ T (TC17) cells, which probably contributed to toxicity; this points towards a specific mode of action for the emergence of irAEs. A modified treatment regimen with intermittent dosing of PI3Kδi in mouse models led to a significant decrease in tumour growth without inducing pathogenic T cells in colonic tissue, indicating that alternative dosing regimens might limit toxicity.

TWEAK functions with TNF and IL-17 on keratinocytes and is a potential target for psoriasis therapy.

TNF and IL-17 are two cytokines that drive dysregulated keratinocyte activity, and their targeting is highly efficacious in patients with psoriasis, but whether these molecules act with other inflammatory factors is not clear. Here, we show that mice having a keratinocyte-specific deletion of Fn14 (Tnfrsf12a), the receptor for the TNF superfamily cytokine TWEAK (Tnfsf12), displayed reduced imiquimod-induced skin inflammation, including diminished epidermal hyperplasia and less expression of psoriasis signature genes. This corresponded with Fn14 being expressed in keratinocytes in human psoriasis lesions and TWEAK being found in several subsets of skin cells. Transcriptomic studies in human keratinocytes revealed that TWEAK strongly overlaps with IL-17A and TNF in up-regulating the expression of CXC chemokines, along with cytokines such as IL-23 and inflammation-associated proteins like S100A8/9 and SERPINB1/B9, all previously found to be highly expressed in the lesional skin of patients with psoriasis. TWEAK displayed strong synergism with TNF or IL-17A in up-regulating messenger RNA for many psoriasis-associated genes in human keratinocytes, including IL23A, IL36G, and multiple chemokines, implying that TWEAK acts with TNF and IL-17 to enhance feedback inflammatory activity. Correspondingly, therapeutic treatment of mice with anti-TWEAK was equally as effective as antibodies to IL-17A or TNF in reducing clinical and immunological features of psoriasis-like skin inflammation and combination targeting of TWEAK with either cytokine had no greater inhibitory effect, reinforcing the conclusion that all three cytokines function together. Thus, blocking TWEAK could be comparable to targeting TNF or IL-17 and might be considered as an alternate therapeutic treatment for psoriasis.

Intratumoral follicular regulatory T cells curtail anti-PD-1 treatment efficacy.

Immune-checkpoint blockade (ICB) has shown remarkable clinical success in boosting antitumor immunity. However, the breadth of its cellular targets and specific mode of action remain elusive. We find that tumor-infiltrating follicular regulatory T (TFR) cells are prevalent in tumor tissues of several cancer types. They are primarily located within tertiary lymphoid structures and exhibit superior suppressive capacity and in vivo persistence as compared with regulatory T cells, with which they share a clonal and developmental relationship. In syngeneic tumor models, anti-PD-1 treatment increases the number of tumor-infiltrating TFR cells. Both TFR cell deficiency and the depletion of TFR cells with anti-CTLA-4 before anti-PD-1 treatment improve tumor control in mice. Notably, in a cohort of 271 patients with melanoma, treatment with anti-CTLA-4 followed by anti-PD-1 at progression was associated with better a survival outcome than monotherapy with anti-PD-1 or anti-CTLA-4, anti-PD-1 followed by anti-CTLA-4 at progression or concomitant combination therapy.

The tumor suppressor kinase DAPK3 drives tumor-intrinsic immunity through the STING-IFN-ß pathway.

Evasion of host immunity is a hallmark of cancer; however, mechanisms linking oncogenic mutations and immune escape are incompletely understood. Through loss-of-function screening of 1,001 tumor suppressor genes, we identified death-associated protein kinase 3 (DAPK3) as a previously unrecognized driver of anti-tumor immunity through the stimulator of interferon genes (STING) pathway of cytosolic DNA sensing. Loss of DAPK3 expression or kinase activity impaired STING activation and interferon (IFN)-ß-stimulated gene induction. DAPK3 deficiency in IFN-ß-producing tumors drove rapid growth and reduced infiltration of CD103+CD8α+ dendritic cells and cytotoxic lymphocytes, attenuating the response to cancer chemo-immunotherapy. Mechanistically, DAPK3 coordinated post-translational modification of STING. In unstimulated cells, DAPK3 inhibited STING K48-linked poly-ubiquitination and proteasome-mediated degradation. After cGAMP stimulation, DAPK3 was required for STING K63-linked poly-ubiquitination and STING-TANK-binding kinase 1 interaction. Comprehensive phospho-proteomics uncovered a DAPK3-specific phospho-site on the E3 ligase LMO7, critical for LMO7-STING interaction and STING K63-linked poly-ubiquitination. Thus, DAPK3 is an essential kinase for STING activation that drives tumor-intrinsic innate immunity and tumor immune surveillance.

NSD2 overexpression drives clustered chromatin and transcriptional changes in a subset of insulated domains.

CTCF and cohesin play a key role in organizing chromatin into topologically associating domain (TAD) structures. Disruption of a single CTCF binding site is sufficient to change chromosomal interactions leading to alterations in chromatin modifications and gene regulation. However, the extent to which alterations in chromatin modifications can disrupt 3D chromosome organization leading to transcriptional changes is unknown. In multiple myeloma, a 4;14 translocation induces overexpression of the histone methyltransferase, NSD2, resulting in expansion of H3K36me2 and shrinkage of antagonistic H3K27me3 domains. Using isogenic cell lines producing high and low levels of NSD2, here we find oncogene activation is linked to alterations in H3K27ac and CTCF within H3K36me2 enriched chromatin. A logistic regression model reveals that differentially expressed genes are significantly enriched within the same insulated domain as altered H3K27ac and CTCF peaks. These results identify a bidirectional relationship between 2D chromatin and 3D genome organization in gene regulation.

Single-cell transcriptomic analysis of tissue-resident memory T cells in human lung cancer.

High numbers of tissue-resident memory T (TRM) cells are associated with better clinical outcomes in cancer patients. However, the molecular characteristics that drive their efficient immune response to tumors are poorly understood. Here, single-cell and bulk transcriptomic analysis of TRM and non-TRM cells present in tumor and normal lung tissue from patients with lung cancer revealed that PD-1-expressing TRM cells in tumors were clonally expanded and enriched for transcripts linked to cell proliferation and cytotoxicity when compared with PD-1-expressing non-TRM cells. This feature was more prominent in the TRM cell subset coexpressing PD-1 and TIM-3, and it was validated by functional assays ex vivo and also reflected in their chromatin accessibility profile. This PD-1+TIM-3+ TRM cell subset was enriched in responders to PD-1 inhibitors and in tumors with a greater magnitude of CTL responses. These data highlight that not all CTLs expressing PD-1 are dysfunctional; on the contrary, TRM cells with PD-1 expression were enriched for features suggestive of superior functionality.

Human Eosinophils Express a Distinct Gene Expression Program in Response to IL-3 Compared with Common ß-Chain Cytokines IL-5 and GM-CSF.

Despite recent advances in asthma management with anti-IL-5 therapies, many patients have eosinophilic asthma that remains poorly controlled. IL-3 shares a common ß subunit receptor with both IL-5 and GM-CSF but, through α-subunit-specific properties, uniquely influences eosinophil biology and may serve as a potential therapeutic target. We aimed to globally characterize the transcriptomic profiles of GM-CSF, IL-3, and IL-5 stimulation on human circulating eosinophils and identify differences in gene expression using advanced statistical modeling. Human eosinophils were isolated from the peripheral blood of healthy volunteers and stimulated with either GM-CSF, IL-3, or IL-5 for 48 h. RNA was then extracted and bulk sequencing performed. DESeq analysis identified differentially expressed genes and weighted gene coexpression network analysis independently defined modules of genes that are highly coexpressed. GM-CSF, IL-3, and IL-5 commonly upregulated 252 genes and downregulated 553 genes, producing a proinflammatory and survival phenotype that was predominantly mediated through TWEAK signaling. IL-3 stimulation yielded the most numbers of differentially expressed genes that were also highly coexpressed (n = 119). These genes were enriched in pathways involving JAK/STAT signaling. GM-CSF and IL-5 stimulation demonstrated redundancy in eosinophil gene expression. In conclusion, IL-3 produces a distinct eosinophil gene expression program among the ß-chain receptor cytokines. IL-3-upregulated genes may provide a foundation for research into therapeutics for patients with eosinophilic asthma who do not respond to anti-IL-5 therapies.

TET enzymes augment activation-induced deaminase (AID) expression via 5-hydroxymethylcytosine modifications at the Aicda superenhancer.

TET enzymes are dioxygenases that promote DNA demethylation by oxidizing the methyl group of 5-methylcytosine to 5-hydroxymethylcytosine (5hmC). Here, we report a close correspondence between 5hmC-marked regions, chromatin accessibility and enhancer activity in B cells, and a strong enrichment for consensus binding motifs for basic region-leucine zipper (bZIP) transcription factors at TET-responsive genomic regions. Functionally, Tet2 and Tet3 regulate class switch recombination (CSR) in murine B cells by enhancing expression of Aicda, which encodes the activation-induced cytidine deaminase (AID) enzyme essential for CSR. TET enzymes deposit 5hmC, facilitate DNA demethylation, and maintain chromatin accessibility at two TET-responsive enhancer elements, TetE1 and TetE2, located within a superenhancer in the Aicda locus. Our data identify the bZIP transcription factor, ATF-like (BATF) as a key transcription factor involved in TET-dependent Aicda expression. 5hmC is not deposited at TetE1 in activated Batf-deficient B cells, indicating that BATF facilitates TET recruitment to this Aicda enhancer. Our study emphasizes the importance of TET enzymes for bolstering AID expression and highlights 5hmC as an epigenetic mark that captures enhancer dynamics during cell activation.

Measuring the reproducibility and quality of Hi-C data.

BACKGROUND: Hi-C is currently the most widely used assay to investigate the 3D organization of the genome and to study its role in gene regulation, DNA replication, and disease. However, Hi-C experiments are costly to perform and involve multiple complex experimental steps; thus, accurate methods for measuring the quality and reproducibility of Hi-C data are essential to determine whether the output should be used further in a study. RESULTS: Using real and simulated data, we profile the performance of several recently proposed methods for assessing reproducibility of population Hi-C data, including HiCRep, GenomeDISCO, HiC-Spector, and QuASAR-Rep. By explicitly controlling noise and sparsity through simulations, we demonstrate the deficiencies of performing simple correlation analysis on pairs of matrices, and we show that methods developed specifically for Hi-C data produce better measures of reproducibility. We also show how to use established measures, such as the ratio of intra- to interchromosomal interactions, and novel ones, such as QuASAR-QC, to identify low-quality experiments. CONCLUSIONS: In this work, we assess reproducibility and quality measures by varying sequencing depth, resolution and noise levels in Hi-C data from 13 cell lines, with two biological replicates each, as well as 176 simulated matrices. Through this extensive validation and benchmarking of Hi-C data, we describe best practices for reproducibility and quality assessment of Hi-C experiments. We make all software publicly available at http://github.com/kundajelab/3DChromatin_ReplicateQC to facilitate adoption in the community.

Reduced expression of phosphatase PTPN2 promotes pathogenic conversion of Tregs in autoimmunity.

Genetic variants at the PTPN2 locus, which encodes the tyrosine phosphatase PTPN2, cause reduced gene expression and are linked to rheumatoid arthritis (RA) and other autoimmune diseases. PTPN2 inhibits signaling through the T cell and cytokine receptors, and loss of PTPN2 promotes T cell expansion and CD4- and CD8-driven autoimmunity. However, it remains unknown whether loss of PTPN2 in FoxP3+ regulatory T cells (Tregs) plays a role in autoimmunity. Here we aimed to model human autoimmune-predisposing PTPN2 variants, the presence of which results in a partial loss of PTPN2 expression, in mouse models of RA. We identified that reduced expression of Ptpn2 enhanced the severity of autoimmune arthritis in the T cell-dependent SKG mouse model and demonstrated that this phenotype was mediated through a Treg-intrinsic mechanism. Mechanistically, we found that through dephosphorylation of STAT3, PTPN2 inhibits IL-6-driven pathogenic loss of FoxP3 after Tregs have acquired RORγt expression, at a stage when chromatin accessibility for STAT3-targeted IL-17-associated transcription factors is maximized. We conclude that PTPN2 promotes FoxP3 stability in mouse RORγt+ Tregs and that loss of function of PTPN2 in Tregs contributes to the association between PTPN2 and autoimmunity.

Integrative detection and analysis of structural variation in cancer genomes.

Structural variants (SVs) can contribute to oncogenesis through a variety of mechanisms. Despite their importance, the identification of SVs in cancer genomes remains challenging. Here, we present a framework that integrates optical mapping, high-throughput chromosome conformation capture (Hi-C), and whole-genome sequencing to systematically detect SVs in a variety of normal or cancer samples and cell lines. We identify the unique strengths of each method and demonstrate that only integrative approaches can comprehensively identify SVs in the genome. By combining Hi-C and optical mapping, we resolve complex SVs and phase multiple SV events to a single haplotype. Furthermore, we observe widespread structural variation events affecting the functions of noncoding sequences, including the deletion of distal regulatory sequences, alteration of DNA replication timing, and the creation of novel three-dimensional chromatin structural domains. Our results indicate that noncoding SVs may be underappreciated mutational drivers in cancer genomes.

Identification of copy number variations and translocations in cancer cells from Hi-C data.

MOTIVATION: Eukaryotic chromosomes adapt a complex and highly dynamic three-dimensional (3D) structure, which profoundly affects different cellular functions and outcomes including changes in epigenetic landscape and in gene expression. Making the scenario even more complex, cancer cells harbor chromosomal abnormalities [e.g. copy number variations (CNVs) and translocations] altering their genomes both at the sequence level and at the level of 3D organization. High-throughput chromosome conformation capture techniques (e.g. Hi-C), which are originally developed for decoding the 3D structure of the chromatin, provide a great opportunity to simultaneously identify the locations of genomic rearrangements and to investigate the 3D genome organization in cancer cells. Even though Hi-C data has been used for validating known rearrangements, computational methods that can distinguish rearrangement signals from the inherent biases of Hi-C data and from the actual 3D conformation of chromatin, and can precisely detect rearrangement locations de novo have been missing. RESULTS: In this work, we characterize how intra and inter-chromosomal Hi-C contacts are distributed for normal and rearranged chromosomes to devise a new set of algorithms (i) to identify genomic segments that correspond to CNV regions such as amplifications and deletions (HiCnv), (ii) to call inter-chromosomal translocations and their boundaries (HiCtrans) from Hi-C experiments and (iii) to simulate Hi-C data from genomes with desired rearrangements and abnormalities (AveSim) in order to select optimal parameters for and to benchmark the accuracy of our methods. Our results on 10 different cancer cell lines with Hi-C data show that we identify a total number of 105 amplifications and 45 deletions together with 90 translocations, whereas we identify virtually no such events for two karyotypically normal cell lines. Our CNV predictions correlate very well with whole genome sequencing data among chromosomes with CNV events for a breast cancer cell line (r = 0.89) and capture most of the CNVs we simulate using Avesim. For HiCtrans predictions, we report evidence from the literature for 30 out of 90 translocations for eight of our cancer cell lines. Furthermore, we show that our tools identify and correctly classify relatively understudied rearrangements such as double minutes and homogeneously staining regions. Considering the inherent limitations of existing techniques for karyotyping (i.e. missing balanced rearrangements and those near repetitive regions), the accurate identification of CNVs and translocations in a cost-effective and high-throughput setting is still a challenge. Our results show that the set of tools we develop effectively utilize moderately sequenced Hi-C libraries (100-300 million reads) to identify known and de novo chromosomal rearrangements/abnormalities in well-established cancer cell lines. With the decrease in required number of cells and the increase in attainable resolution, we believe that our framework will pave the way towards comprehensive mapping of genomic rearrangements in primary cells from cancer patients using Hi-C. AVAILABILITY AND IMPLEMENTATION: CNV calling: https://github.com/ay-lab/HiCnv, Translocation calling: https://github.com/ay-lab/HiCtrans and Hi-C simulation: https://github.com/ay-lab/AveSim. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Topologically associating domains and their long-range contacts are established during early G1 coincident with the establishment of the replication-timing program.

Mammalian genomes are partitioned into domains that replicate in a defined temporal order. These domains can replicate at similar times in all cell types (constitutive) or at cell type-specific times (developmental). Genome-wide chromatin conformation capture (Hi-C) has revealed sub-megabase topologically associating domains (TADs), which are the structural counterparts of replication domains. Hi-C also segregates inter-TAD contacts into defined 3D spatial compartments that align precisely to genome-wide replication timing profiles. Determinants of the replication-timing program are re-established during early G1 phase of each cell cycle and lost in G2 phase, but it is not known when TAD structure and inter-TAD contacts are re-established after their elimination during mitosis. Here, we use multiplexed 4C-seq to study dynamic changes in chromatin organization during early G1. We find that both establishment of TADs and their compartmentalization occur during early G1, within the same time frame as establishment of the replication-timing program. Once established, this 3D organization is preserved either after withdrawal into quiescence or for the remainder of interphase including G2 phase, implying 3D structure is not sufficient to maintain replication timing. Finally, we find that developmental domains are less well compartmentalized than constitutive domains and display chromatin properties that distinguish them from early and late constitutive domains. Overall, this study uncovers a strong connection between chromatin re-organization during G1, establishment of replication timing, and its developmental control.

Identifying multi-locus chromatin contacts in human cells using tethered multiple 3C.

BACKGROUND: Several recently developed experimental methods, each an extension of the chromatin conformation capture (3C) assay, have enabled the genome-wide profiling of chromatin contacts between pairs of genomic loci in 3D. Especially in complex eukaryotes, data generated by these methods, coupled with other genome-wide datasets, demonstrated that non-random chromatin folding correlates strongly with cellular processes such as gene expression and DNA replication. RESULTS: We describe a genome architecture assay, tethered multiple 3C (TM3C), that maps genome-wide chromatin contacts via a simple protocol of restriction enzyme digestion and religation of fragments upon agarose gel beads followed by paired-end sequencing. In addition to identifying contacts between pairs of loci, TM3C enables identification of contacts among more than two loci simultaneously. We use TM3C to assay the genome architectures of two human cell lines: KBM7, a near-haploid chronic leukemia cell line, and NHEK, a normal diploid human epidermal keratinocyte cell line. We confirm that the contact frequency maps produced by TM3C exhibit features characteristic of existing genome architecture datasets, including the expected scaling of contact probabilities with genomic distance, megabase scale chromosomal compartments and sub-megabase scale topological domains. We also confirm that TM3C captures several known cell type-specific contacts, ploidy shifts and translocations, such as Philadelphia chromosome formation (Ph+) in KBM7. We confirm a subset of the triple contacts involving the IGF2-H19 imprinting control region (ICR) using PCR analysis for KBM7 cells. Our genome-wide analysis of pairwise and triple contacts demonstrates their preference for linking open chromatin regions to each other and for linking regions with higher numbers of DNase hypersensitive sites (DHSs) to each other. For near-haploid KBM7 cells, we infer whole genome 3D models that exhibit clustering of small chromosomes with each other and large chromosomes with each other, consistent with previous studies of the genome architectures of other human cell lines. CONCLUSION: TM3C is a simple protocol for ascertaining genome architecture and can be used to identify simultaneous contacts among three or four loci. Application of TM3C to a near-haploid human cell line revealed large-scale features of chromosomal organization and multi-way chromatin contacts that preferentially link regions of open chromatin.

Joint annotation of chromatin state and chromatin conformation reveals relationships among domain types and identifies domains of cell-type-specific expression.

The genomic neighborhood of a gene influences its activity, a behavior that is attributable in part to domain-scale regulation. Previous genomic studies have identified many types of regulatory domains. However, due to the difficulty of integrating genomics data sets, the relationships among these domain types are poorly understood. Semi-automated genome annotation (SAGA) algorithms facilitate human interpretation of heterogeneous collections of genomics data by simultaneously partitioning the human genome and assigning labels to the resulting genomic segments. However, existing SAGA methods cannot integrate inherently pairwise chromatin conformation data. We developed a new computational method, called graph-based regularization (GBR), for expressing a pairwise prior that encourages certain pairs of genomic loci to receive the same label in a genome annotation. We used GBR to exploit chromatin conformation information during genome annotation by encouraging positions that are close in 3D to occupy the same type of domain. Using this approach, we produced a model of chromatin domains in eight human cell types, thereby revealing the relationships among known domain types. Through this model, we identified clusters of tightly regulated genes expressed in only a small number of cell types, which we term "specific expression domains." We found that domain boundaries marked by promoters and CTCF motifs are consistent between cell types even when domain activity changes. Finally, we showed that GBR can be used to transfer information from well-studied cell types to less well-characterized cell types during genome annotation, making it possible to produce high-quality annotations of the hundreds of cell types with limited available data.