Extensive promoter-centered chromatin interactions provide a topological basis for transcription regulation.

Higher-order chromosomal organization for transcription regulation is poorly understood in eukaryotes. Using genome-wide Chromatin Interaction Analysis with Paired-End-Tag sequencing (ChIA-PET), we mapped long-range chromatin interactions associated with RNA polymerase II in human cells and uncovered widespread promoter-centered intragenic, extragenic, and intergenic interactions. These interactions further aggregated into higher-order clusters, wherein proximal and distal genes were engaged through promoter-promoter interactions. Most genes with promoter-promoter interactions were active and transcribed cooperatively, and some interacting promoters could influence each other implying combinatorial complexity of transcriptional controls. Comparative analyses of different cell lines showed that cell-specific chromatin interactions could provide structural frameworks for cell-specific transcription, and suggested significant enrichment of enhancer-promoter interactions for cell-specific functions. Furthermore, genetically-identified disease-associated noncoding elements were found to be spatially engaged with corresponding genes through long-range interactions. Overall, our study provides insights into transcription regulation by three-dimensional chromatin interactions for both housekeeping and cell-specific genes in human cells.

An AR-ERG transcriptional signature defined by long-range chromatin interactomes in prostate cancer cells.

The aberrant activities of transcription factors such as the androgen receptor (AR) underpin prostate cancer development. While the AR cis-regulation has been extensively studied in prostate cancer, information pertaining to the spatial architecture of the AR transcriptional circuitry remains limited. In this paper, we propose a novel framework to profile long-range chromatin interactions associated with AR and its collaborative transcription factor, erythroblast transformation-specific related gene (ERG), using chromatin interaction analysis by paired-end tag (ChIA-PET). We identified ERG-associated long-range chromatin interactions as a cooperative component in the AR-associated chromatin interactome, acting in concert to achieve coordinated regulation of a subset of AR target genes. Through multifaceted functional data analysis, we found that AR-ERG interaction hub regions are characterized by distinct functional signatures, including bidirectional transcription and cotranscription factor binding. In addition, cancer-associated long noncoding RNAs were found to be connected near protein-coding genes through AR-ERG looping. Finally, we found strong enrichment of prostate cancer genome-wide association study (GWAS) single nucleotide polymorphisms (SNPs) at AR-ERG co-binding sites participating in chromatin interactions and gene regulation, suggesting GWAS target genes identified from chromatin looping data provide more biologically relevant findings than using the nearest gene approach. Taken together, our results revealed an AR-ERG-centric higher-order chromatin structure that drives coordinated gene expression in prostate cancer progression and the identification of potential target genes for therapeutic intervention.

Concurrent Single-Cell RNA and Targeted DNA Sequencing on an Automated Platform for Comeasurement of Genomic and Transcriptomic Signatures.

BACKGROUND: The comeasurement of both genomic and transcriptomic signatures in single cells is of fundamental importance to accurately assess how the genetic information correlates with the transcriptomic phenotype. However, existing technologies have low throughput and laborious work flows. METHODS: We developed a new method for concurrent sequencing of the transcriptome and targeted genomic regions (CORTAD-seq) within the same single cell on an automated microfluidic platform. The method was compatible with the downstream library preparation, allowing easy integration into existing next-generation sequencing work flows. We incorporated a single-cell bioinformatics pipeline for transcriptome and mutation analysis. RESULTS: As proof of principle, we applied CORTAD-seq to lung cancer cell lines to dissect the cellular consequences of mutations that result in resistance to targeted therapy. We obtained a mean detection of 6000 expressed genes and an exonic rate of 50%. The targeted DNA-sequencing data achieved a 97.8% detection rate for mutations and allowed for the identification of copy number variations and haplotype construction. We detected expression signatures of tyrosine kinase inhibitor (TKI) resistance, epidermal growth factor receptor (EGFR) amplification, and expansion of the T790M mutation among resistant cells. We also identified characteristics for TKI resistance that were independent of EGFR T790M, indicating that other alterations are required for resistance in this context. CONCLUSIONS: CORTAD-seq allows assessment of the interconnection between genetic and transcriptomic changes in single cells. It is operated on an automated, commercially available single-cell isolation platform, making its implementation straightforward.

Chromatin connectivity maps reveal dynamic promoter-enhancer long-range associations.

In multicellular organisms, transcription regulation is one of the central mechanisms modelling lineage differentiation and cell-fate determination. Transcription requires dynamic chromatin configurations between promoters and their corresponding distal regulatory elements. It is believed that their communication occurs within large discrete foci of aggregated RNA polymerases termed transcription factories in three-dimensional nuclear space. However, the dynamic nature of chromatin connectivity has not been characterized at the genome-wide level. Here, through a chromatin interaction analysis with paired-end tagging approach using an antibody that primarily recognizes the pre-initiation complexes of RNA polymerase II, we explore the transcriptional interactomes of three mouse cells of progressive lineage commitment, including pluripotent embryonic stem cells, neural stem cells and neurosphere stem/progenitor cells. Our global chromatin connectivity maps reveal approximately 40,000 long-range interactions, suggest precise enhancer-promoter associations and delineate cell-type-specific chromatin structures. Analysis of the complex regulatory repertoire shows that there are extensive colocalizations among promoters and distal-acting enhancers. Most of the enhancers associate with promoters located beyond their nearest active genes, indicating that the linear juxtaposition is not the only guiding principle driving enhancer target selection. Although promoter-enhancer interactions exhibit high cell-type specificity, promoters involved in interactions are found to be generally common and mostly active among different cells. Chromatin connectivity networks reveal that the pivotal genes of reprogramming functions are transcribed within physical proximity to each other in embryonic stem cells, linking chromatin architecture to coordinated gene expression. Our study sets the stage for the full-scale dissection of spatial and temporal genome structures and their roles in orchestrating development.

TNFα signals through specialized factories where responsive coding and miRNA genes are transcribed.

Tumour necrosis factor alpha (TNFα) is a potent cytokine that signals through nuclear factor kappa B (NFκB) to activate a subset of human genes. It is usually assumed that this involves RNA polymerases transcribing responsive genes wherever they might be in the nucleus. Using primary human endothelial cells, variants of chromosome conformation capture (including 4C and chromatin interaction analysis with paired-end tag sequencing), and fluorescence in situ hybridization to detect single nascent transcripts, we show that TNFα induces responsive genes to congregate in discrete 'NFκB factories'. Some factories further specialize in transcribing responsive genes encoding micro-RNAs that target downregulated mRNAs. We expect all signalling pathways to contain this extra leg, where responding genes are transcribed in analogous specialized factories.

Integrative genomic analyses of European intrahepatic cholangiocarcinoma: Novel ROS1 fusion gene and PBX1 as prognostic marker.

BACKGROUND: Cholangiocarcinoma (CCA) is a fatal cancer of the bile duct with a poor prognosis owing to limited therapeutic options. The incidence of intrahepatic CCA (iCCA) is increasing worldwide, and its molecular basis is emerging. Environmental factors may contribute to regional differences in the mutation spectrum of European patients with iCCA, which are underrepresented in systematic genomic and transcriptomic studies of the disease. METHODS: We describe an integrated whole-exome sequencing and transcriptomic study of 37 iCCAs patients in Germany. RESULTS: We observed as most frequently mutated genes ARID1A (14%), IDH1, BAP1, TP53, KRAS, and ATM in 8% of patients. We identified FGFR2::BICC1 fusions in two tumours, and FGFR2::KCTD1 and TMEM106B::ROS1 as novel fusions with potential therapeutic implications in iCCA and confirmed oncogenic properties of TMEM106B::ROS1 in vitro. Using a data integration framework, we identified PBX1 as a novel central regulatory gene in iCCA. We performed extended screening by targeted sequencing of an additional 40 CCAs. In the joint analysis, IDH1 (13%), BAP1 (10%), TP53 (9%), KRAS (7%), ARID1A (7%), NF1 (5%), and ATM (5%) were the most frequently mutated genes, and we found PBX1 to show copy gain in 20% of the tumours. According to other studies, amplifications of PBX1 tend to occur in European iCCAs in contrast to liver fluke-associated Asian iCCAs. CONCLUSIONS: By analyzing an additional European cohort of iCCA patients, we found that PBX1 protein expression was a marker of poor prognosis. Overall, our findings provide insight into key molecular alterations in iCCA, reveal new targetable fusion genes, and suggest that PBX1 is a novel modulator of this disease.

ChIA-PET analysis of transcriptional chromatin interactions.

Long-range chromatin contacts between specific DNA regulatory elements play a pivotal role in gene expression regulation, and a global characterization of these interactions in the 3-dimensional (3D) chromatin structure is imperative in understanding signaling networks and cell states. Chromatin Interaction Analysis using Paired-End Tag sequencing (ChIA-PET) is a method which converts functional chromatin structure into millions of short tag sequences. Combining Chromatin Immunoprecipitation (ChIP), proximity ligation and high-throughput sequencing, ChIA-PET provides a global and unbiased interrogation of higher-order chromatin structures associated with specific protein factors. Here, we describe the detailed procedures of the ChIA-PET methodology, unraveling transcription-associated chromatin contacts in a model human cell line.

Experimental and bioinformatics considerations in cancer application of single cell genomics.

Single cell genomics offers an unprecedented resolution to interrogate genetic heterogeneity in a patient's tumour at the intercellular level. However, the DNA yield per cell is insufficient for today's sequencing library preparation protocols. This necessitates DNA amplification which is a key source of experimental noise. We provide an evaluation of two protocols using micro-fluidics based amplification for whole exome sequencing, which is an experimental scenario commonly used in single cell genomics. The results highlight their respective biases and relative strengths in identification of single nucleotide variations. Towards this end, we introduce a workflow SoVaTSiC, which allows for quality evaluation and somatic variant identification of single cell data. As proof of concept, the framework was applied to study a lung adenocarcinoma tumour. The analysis provides insights into tumour phylogeny by identifying key mutational events in lung adenocarcinoma evolution. The consequence of this inference is supported by the histology of the tumour and demonstrates usefulness of the approach.

Complementary Sequential Circulating Tumor Cell (CTC) and Cell-Free Tumor DNA (ctDNA) Profiling Reveals Metastatic Heterogeneity and Genomic Changes in Lung Cancer and Breast Cancer.

INTRODUCTION: Circulating tumor cells (CTCs) and cell-free tumor DNA (ctDNA) are tumor components present in circulation. Due to the limited access to both CTC enrichment platforms and ctDNA sequencing in most laboratories, they are rarely analyzed together. METHODS: Concurrent isolation of ctDNA and single CTCs were isolated from lung cancer and breast cancer patients using the combination of size-based and CD45-negative selection method via DropCell platform. We performed targeted amplicon sequencing to evaluate the genomic heterogeneity of CTCs and ctDNA in lung cancer and breast cancer patients. RESULTS: Higher degrees of genomic heterogeneity were observed in CTCs as compared to ctDNA. Several shared alterations present in CTCs and ctDNA were undetected in the primary tumor, highlighting the intra-tumoral heterogeneity of tumor components that were shed into systemic circulation. Accordingly, CTCs and ctDNA displayed higher degree of concordance with the metastatic tumor than the primary tumor. The alterations detected in circulation correlated with worse survival outcome for both lung and breast cancer patients emphasizing the impact of the metastatic phenotype. Notably, evolving genetic signatures were detected in the CTCs and ctDNA samples during the course of treatment and disease progression. CONCLUSIONS: A standardized sample processing and data analysis workflow for concurrent analysis of CTCs and ctDNA successfully dissected the heterogeneity of metastatic tumor in circulation as well as the progressive genomic changes that may potentially guide the selection of appropriate therapy against evolving tumor clonality.

Elucidating the genomic architecture of Asian EGFR-mutant lung adenocarcinoma through multi-region exome sequencing.

EGFR-mutant lung adenocarcinomas (LUAD) display diverse clinical trajectories and are characterized by rapid but short-lived responses to EGFR tyrosine kinase inhibitors (TKIs). Through sequencing of 79 spatially distinct regions from 16 early stage tumors, we show that despite low mutation burdens, EGFR-mutant Asian LUADs unexpectedly exhibit a complex genomic landscape with frequent and early whole-genome doubling, aneuploidy, and high clonal diversity. Multiple truncal alterations, including TP53 mutations and loss of CDKN2A and RB1, converge on cell cycle dysregulation, with late sector-specific high-amplitude amplifications and deletions that potentially beget drug resistant clones. We highlight the association between genomic architecture and clinical phenotypes, such as co-occurring truncal drivers and primary TKI resistance. Through comparative analysis with published smoking-related LUAD, we postulate that the high intra-tumor heterogeneity observed in Asian EGFR-mutant LUAD may be contributed by an early dominant driver, genomic instability, and low background mutation rates.

High resolution mapping of enhancer-promoter interactions.

RNA Polymerase II ChIA-PET data has revealed enhancers that are active in a profiled cell type and the genes that the enhancers regulate through chromatin interactions. The most commonly used computational method for analyzing ChIA-PET data, the ChIA-PET Tool, discovers interaction anchors at a spatial resolution that is insufficient to accurately identify individual enhancers. We introduce Germ, a computational method that estimates the likelihood that any two narrowly defined genomic locations are jointly occupied by RNA Polymerase II. Germ takes a blind deconvolution approach to simultaneously estimate the likelihood of RNA Polymerase II occupation as well as a model of the arrangement of read alignments relative to locations occupied by RNA Polymerase II. Both types of information are utilized to estimate the likelihood that RNA Polymerase II jointly occupies any two genomic locations. We apply Germ to RNA Polymerase II ChIA-PET data from embryonic stem cells to identify the genomic locations that are jointly occupied along with transcription start sites. We show that these genomic locations align more closely with features of active enhancers measured by ChIP-Seq than the locations identified using the ChIA-PET Tool. We also apply Germ to RNA Polymerase II ChIA-PET data from motor neuron progenitors. Based on the Germ results, we observe that a combination of cell type specific and cell type independent regulatory interactions are utilized by cells to regulate gene expression.

Recurrent Fusion Genes in Gastric Cancer: CLDN18-ARHGAP26 Induces Loss of Epithelial Integrity.

Genome rearrangements, a hallmark of cancer, can result in gene fusions with oncogenic properties. Using DNA paired-end-tag (DNA-PET) whole-genome sequencing, we analyzed 15 gastric cancers (GCs) from Southeast Asians. Rearrangements were enriched in open chromatin and shaped by chromatin structure. We identified seven rearrangement hot spots and 136 gene fusions. In three out of 100 GC cases, we found recurrent fusions between CLDN18, a tight junction gene, and ARHGAP26, a gene encoding a RHOA inhibitor. Epithelial cell lines expressing CLDN18-ARHGAP26 displayed a dramatic loss of epithelial phenotype and long protrusions indicative of epithelial-mesenchymal transition (EMT). Fusion-positive cell lines showed impaired barrier properties, reduced cell-cell and cell-extracellular matrix adhesion, retarded wound healing, and inhibition of RHOA. Gain of invasion was seen in cancer cell lines expressing the fusion. Thus, CLDN18-ARHGAP26 mediates epithelial disintegration, possibly leading to stomach H(+) leakage, and the fusion might contribute to invasiveness once a cell is transformed.

Cross-enhancement of ANGPTL4 transcription by HIF1 alpha and PPAR beta/delta is the result of the conformational proximity of two response elements.

BACKGROUND: Synergistic transcriptional activation by different stimuli has been reported along with a diverse array of mechanisms, but the full scope of these mechanisms has yet to be elucidated. RESULTS: We present a detailed investigation of hypoxia-inducible factor (HIF) 1 dependent gene expression in endothelial cells which suggests the importance of crosstalk between the peroxisome proliferator-activated receptor (PPAR) ß/δ and HIF signaling axes. A migration assay shows a synergistic interaction between these two stimuli, and we identify angiopoietin-like 4 (ANGPTL4) as a common target gene by using a combination of microarray and ChIP-seq analysis. We profile changes of histone marks at enhancers under hypoxia, PPARß/δ agonist and dual stimulations and these suggest that the spatial proximity of two response elements is the principal cause of the synergistic transcription induction. A newly developed quantitative chromosome conformation capture assay shows the quantitative change of the frequency of proximity of the two response elements. CONCLUSIONS: To the best of our knowledge, this is the first report that two different transcription factors cooperate in transcriptional regulation in a synergistic fashion through conformational change of their common target genes.

Direct evidence for pitavastatin induced chromatin structure change in the KLF4 gene in endothelial cells.

Statins exert atheroprotective effects through the induction of specific transcriptional factors in multiple organs. In endothelial cells, statin-dependent atheroprotective gene up-regulation is mediated by Kruppel-like factor (KLF) family transcription factors. To dissect the mechanism of gene regulation, we sought to determine molecular targets by performing microarray analyses of human umbilical vein endothelial cells (HUVECs) treated with pitavastatin, and KLF4 was determined to be the most highly induced gene. In addition, it was revealed that the atheroprotective genes induced with pitavastatin, such as nitric oxide synthase 3 (NOS3) and thrombomodulin (THBD), were suppressed by KLF4 knockdown. Myocyte enhancer factor-2 (MEF2) family activation is reported to be involved in pitavastatin-dependent KLF4 induction. We focused on MEF2C among the MEF2 family members and identified a novel functional MEF2C binding site 148 kb upstream of the KLF4 gene by chromatin immunoprecipitation along with deep sequencing (ChIP-seq) followed by luciferase assay. By applying whole genome and quantitative chromatin conformation analysis {chromatin interaction analysis with paired end tag sequencing (ChIA-PET), and real time chromosome conformation capture (3C) assay}, we observed that the MEF2C-bound enhancer and transcription start site (TSS) of KLF4 came into closer spatial proximity by pitavastatin treatment. 3D-Fluorescence in situ hybridization (FISH) imaging supported the conformational change in individual cells. Taken together, dynamic chromatin conformation change was shown to mediate pitavastatin-responsive gene induction in endothelial cells.

Chromatin Interaction Analysis with Paired-End Tag Sequencing (ChIA-PET) for mapping chromatin interactions and understanding transcription regulation.

Genomes are organized into three-dimensional structures, adopting higher-order conformations inside the micron-sized nuclear spaces (7, 2, 12). Such architectures are not random and involve interactions between gene promoters and regulatory elements (13). The binding of transcription factors to specific regulatory sequences brings about a network of transcription regulation and coordination (1, 14). Chromatin Interaction Analysis by Paired-End Tag Sequencing (ChIA-PET) was developed to identify these higher-order chromatin structures (5,6). Cells are fixed and interacting loci are captured by covalent DNA-protein cross-links. To minimize non-specific noise and reduce complexity, as well as to increase the specificity of the chromatin interaction analysis, chromatin immunoprecipitation (ChIP) is used against specific protein factors to enrich chromatin fragments of interest before proximity ligation. Ligation involving half-linkers subsequently forms covalent links between pairs of DNA fragments tethered together within individual chromatin complexes. The flanking MmeI restriction enzyme sites in the half-linkers allow extraction of paired end tag-linker-tag constructs (PETs) upon MmeI digestion. As the half-linkers are biotinylated, these PET constructs are purified using streptavidin-magnetic beads. The purified PETs are ligated with next-generation sequencing adaptors and a catalog of interacting fragments is generated via next-generation sequencers such as the Illumina Genome Analyzer. Mapping and bioinformatics analysis is then performed to identify ChIP-enriched binding sites and ChIP-enriched chromatin interactions (8). We have produced a video to demonstrate critical aspects of the ChIA-PET protocol, especially the preparation of ChIP as the quality of ChIP plays a major role in the outcome of a ChIA-PET library. As the protocols are very long, only the critical steps are shown in the video.

Large-scale functional organization of long-range chromatin interaction networks.

Chromatin interactions play important roles in transcription regulation. To better understand the underlying evolutionary and functional constraints of these interactions, we implemented a systems approach to examine RNA polymerase-II-associated chromatin interactions in human cells. We found that 40% of the total genomic elements involved in chromatin interactions converged to a giant, scale-free-like, hierarchical network organized into chromatin communities. The communities were enriched in specific functions and were syntenic through evolution. Disease-associated SNPs from genome-wide association studies were enriched among the nodes with fewer interactions, implying their selection against deleterious interactions by limiting the total number of interactions, a model that we further reconciled using somatic and germline cancer mutation data. The hubs lacked disease-associated SNPs, constituted a nonrandomly interconnected core of key cellular functions, and exhibited lethality in mouse mutants, supporting an evolutionary selection that favored the nonrandom spatial clustering of the least-evolving key genomic domains against random genetic or transcriptional errors in the genome. Altogether, our analyses reveal a systems-level evolutionary framework that shapes functionally compartmentalized and error-tolerant transcriptional regulation of human genome in three dimensions.

The Arabidopsis ARGOS-LIKE gene regulates cell expansion during organ growth.

Cell expansion, and its coordination with cell division, plays a critical role in the growth and development of plant organs. However, the genes controlling cell expansion during organogenesis are largely unknown. Here, we demonstrate that a novel Arabidopsis gene, ARGOS-LIKE (ARL), which has some sequence homology to the ARGOS gene, is involved in this process. Reduced expression or overexpression of ARL in Arabidopsis results in smaller or larger cotyledons and leaves as well as other lateral organs, respectively. Anatomical examination of cotyledons and leaves in ARL transgenic plants demonstrates that the alteration in size can be attributed to changes in cell size rather than cell number, indicating that ARL plays a role in cell expansion-dependent organ growth. ARL is upregulated by brassinosteroid (BR) and this induction is impaired in the BR-insensitive mutant bri1, but not in the BR-deficient mutant det2. Ectopic expression of ARL in bri1-119 partially restores cell growth in cotyledons and leaves. Our results suggest that ARL acts downstream of BRI1 and partially mediates BR-related cell expansion signals during organ growth.