Core-binding factor fusion downregulation of ADAR2 RNA editing contributes to AML leukemogenesis

Adenosine-to-inosine RNA editing, which is catalyzed by adenosine deaminases acting on RNA (ADAR) family of enzymes, ADAR1 and ADAR2, has been shown to contribute to multiple cancers. However, other than the chronic myeloid leukemia blast crisis, relatively little is known about its role in other types of hematological malignancies. Here, we found that ADAR2, but not ADAR1 and ADAR3, was specifically downregulated in the core-binding factor (CBF) acute myeloid leukemia (AML) with t(8;21) or inv(16) translocations. In t(8;21) AML, RUNX1-driven transcription of ADAR2 was repressed by the RUNX1-ETO additional exon 9a fusion protein in a dominant-negative manner. Further functional studies confirmed that ADAR2 could suppress leukemogenesis specifically in t(8;21) and inv16 AML cells dependent on its RNA editing capability. Expression of 2 exemplary ADAR2-regulated RNA editing targets coatomer subunit α and component of oligomeric Golgi complex 3 inhibits the clonogenic growth of human t(8;21) AML cells. Our findings support a hitherto, unappreciated mechanism leading to ADAR2 dysregulation in CBF AML and highlight the functional relevance of loss of ADAR2-mediated RNA editing to CBF AML.

Using stable isotopes to measure the dietary responses of Costa Rican forest birds to agricultural countryside

How human modification of native habitats changes the feeding patterns and nutritional ecology of tropical birds is critical to conserving avian biodiversity, but tropical bird diets are laborious to investigate using the traditional methods of diet analysis. Stable isotope analysis provides a cost-effective and efficient proxy to identify general foraging patterns, especially when dietary shifts spanning multiple trophic levels have occurred due to ecosystem disturbance or transformation. To characterize the diets of forest bird species that persist in tropical agricultural countryside, we compared feather carbon (δ13C) and nitrogen (δ15N) isotope values of four species caught and radio-tracked in a 270 hectare forest reserve, smaller forest remnants (including mature forest, secondary forest, and riparian strips), and coffee plantations in mid-elevation (ca. 800–1,400 m) southern Costa Rica. Bird habitat choice had a significant effect on diet composition as revealed by δ13C and δ15N values. Three of the four species studied showed evidence of significantly reduced consumption of invertebrates in coffee plantations, with the isotope values of two species (Tangara icterocephala and Turdus assimilis) indicating, by comparison, nearly a doubling of invertebrate consumption in forest remnants. Our results suggest that coffee plantations are deficient in invertebrates preferred by forest generalist birds that forage in both native forest remnants and coffee plantations. In this region, typical of mountainous American tropics, small forest remnants and a larger forest reserve provide critical dietary resources for native forest birds that utilize the agricultural countryside.

MCI-frcnn: A deep learning method for topological micro-domain boundary detection

Chromatin structural domains, or topologically associated domains (TADs), are a general organizing principle in chromatin biology. RNA polymerase II (RNAPII) mediates multiple chromatin interactive loops, tethering together as RNAPII-associated chromatin interaction domains (RAIDs) to offer a framework for gene regulation. RAID and TAD alterations have been found to be associated with diseases. They can be further dissected as micro-domains (micro-TADs and micro-RAIDs) by clustering single-molecule chromatin-interactive complexes from next-generation three-dimensional (3D) genome techniques, such as ChIA-Drop. Currently, there are few tools available for micro-domain boundary identification. In this work, we developed the MCI-frcnn deep learning method to train a Faster Region-based Convolutional Neural Network (Faster R-CNN) for micro-domain boundary detection. At the training phase in MCI-frcnn, 50 images of RAIDs from Drosophila RNAPII ChIA-Drop data, containing 261 micro-RAIDs with ground truth boundaries, were trained for 7 days. Using this well-trained MCI-frcnn, we detected micro-RAID boundaries for the input new images, with a fast speed (5.26 fps), high recognition accuracy (AUROC = 0.85, mAP = 0.69), and high boundary region quantification (genomic IoU = 76%). We further applied MCI-frcnn to detect human micro-TADs boundaries using human GM12878 SPRITE data and obtained a high region quantification score (mean gIoU = 85%). In all, the MCI-frcnn deep learning method which we developed in this work is a general tool for micro-domain boundary detection.

MCIBox: a toolkit for single-molecule multi-way chromatin interaction visualization and micro-domains identification

The emerging ligation-free three-dimensional (3D) genome mapping technologies can identify multiplex chromatin interactions with single-molecule precision. These technologies not only offer new insight into high-dimensional chromatin organization and gene regulation, but also introduce new challenges in data visualization and analysis. To overcome these challenges, we developed MCIBox, a toolkit for multi-way chromatin interaction (MCI) analysis, including a visualization tool and a platform for identifying micro-domains with clustered single-molecule chromatin complexes. MCIBox is based on various clustering algorithms integrated with dimensionality reduction methods that can display multiplex chromatin interactions at single-molecule level, allowing users to explore chromatin extrusion patterns and super-enhancers regulation modes in transcription, and to identify single-molecule chromatin complexes that are clustered into micro-domains. Furthermore, MCIBox incorporates a two-dimensional kernel density estimation algorithm to identify micro-domains boundaries automatically. These micro-domains were stratified with distinctive signatures of transcription activity and contained different cell-cycle-associated genes. Taken together, MCIBox represents an invaluable tool for the study of multiple chromatin interactions and inaugurates a previously unappreciated view of 3D genome structure.

Accurate Identification of DNA Replication Origin by Fusing Epigenomics and Chromatin Interaction Information

DNA replication initiation is a complex process involving various genetic and epigenomic signatures. The correct identification of replication origins (ORIs) could provide important clues for the study of a variety of diseases caused by replication. Here, we design a computational approach named iORI-Epi to recognize ORIs by incorporating epigenome-based features, sequence-based features, and 3D genome-based features. The iORI-Epi displays excellent robustness and generalization ability on both training datasets and independent datasets of K562 cell line. Further experiments confirm that iORI-Epi is highly scalable in other cell lines (MCF7 and HCT116). We also analyze and clarify the regulatory role of epigenomic marks, DNA motifs, and chromatin interaction in DNA replication initiation of eukaryotic genomes. Finally, we discuss gene enrichment pathways from the perspective of ORIs in different replication timing states and heuristically dissect the effect of promoters on replication initiation. Our computational methodology is worth extending to ORI identification in other eukaryotic species.

Chromatin Rewiring by Mismatch Repair Protein MSH2 Alters Cell Adhesion Pathways and Sensitivity to BET Inhibition in Gastric Cancer

Mutations in the DNA mismatch repair gene MSH2 are causative of microsatellite instability (MSI) in multiple cancers. Here, we discovered that besides its well-established role in DNA repair, MSH2 exerts a novel epigenomic function in gastric cancer. Unbiased CRISPR-based mass spectrometry combined with genome-wide CRISPR functional screening revealed that in early-stage gastric cancer MSH2 genomic binding is not randomly distributed but rather is associated specifically with tumor-associated super-enhancers controlling the expression of cell adhesion genes. At these loci, MSH2 genomic binding was required for chromatin rewiring, de novo enhancer-promoter interactions, maintenance of histone acetylation levels, and regulation of cell adhesion pathway expression. The chromatin function of MSH2 was independent of its DNA repair catalytic activity but required MSH6, another DNA repair gene, and recruitment to gene loci by the SWI/SNF chromatin remodeler SMARCA4/BRG1. Loss of MSH2 in advanced gastric cancers was accompanied by deficient cell adhesion pathway expression, epithelial-mesenchymal transition, and enhanced tumorigenesis in vitro and in vivo. However, MSH2-deficient gastric cancers also displayed addiction to BAZ1B, a bromodomain-containing family member, and consequent synthetic lethality to bromodomain and extraterminal motif (BET) inhibition. Our results reveal a role for MSH2 in gastric cancer epigenomic regulation and identify BET inhibition as a potential therapy in MSH2-deficient gastric malignancies.

SIGNIFICANCE

DNA repair protein MSH2 binds and regulates cell adhesion genes by enabling enhancer-promoter interactions, and loss of MSH2 causes deficient cell adhesion and bromodomain and extraterminal motif inhibitor synthetic lethality in gastric cancer.

Pan-cancer pervasive upregulation of 3' UTR splicing drives tumourigenesis

Most mammalian genes generate messenger RNAs with variable untranslated regions (UTRs) that are important post-transcriptional regulators. In cancer, shortening at 3′ UTR ends via alternative polyadenylation can activate oncogenes. However, internal 3′ UTR splicing remains poorly understood as splicing studies have traditionally focused on protein-coding alterations. Here we systematically map the pan-cancer landscape of 3′ UTR splicing and present this in SpUR (http://www.cbrc.kaust.edu.sa/spur/home/). 3′ UTR splicing is widespread, upregulated in cancers, correlated with poor prognosis and more prevalent in oncogenes. We show that antisense oligonucleotide-mediated inhibition of 3′ UTR splicing efficiently reduces oncogene expression and impedes tumour progression. Notably, CTNNB1 3′ UTR splicing is the most consistently dysregulated event across cancers. We validate its upregulation in hepatocellular carcinoma and colon adenocarcinoma, and show that the spliced 3′ UTR variant is the predominant contributor to its oncogenic functions. Overall, our study highlights the importance of 3′ UTR splicing in cancer and may launch new avenues for RNA-based anti-cancer therapeutics.

Chan et al. report that 3′ UTR splicing is widespread and enhanced across different cancer types and is associated with more advanced tumour progression.

MYC overexpression leads to increased chromatin interactions at superenhancers and MYC binding sites

The MYC oncogene encodes for the MYC protein and is frequently dysregulated across multiple cancer cell types, making it an attractive target for cancer therapy. MYC overexpression leads to MYC binding at active enhancers, resulting in a global transcriptional amplification of active genes. Because super-enhancers are frequently dysregulated in cancer, we hypothesized that MYC preferentially invades into super-enhancers and alters the cancer genome organization. To that end, we performed ChIP-seq, RNA-seq, circular chromosome conformation capture (4C-seq), and Spike-in Quantitative Hi-C (SIQHiC) on the U2OS osteosarcoma cell line with tetracycline-inducible MYC. MYC overexpression in U2OS cells modulated histone acetylation and increased MYC binding at super-enhancers. SIQHiC analysis revealed increased global chromatin contact frequency, particularly at chromatin interactions connecting MYC binding sites at promoters and enhancers. Immunofluorescence staining showed that MYC molecules formed punctate foci at these transcriptionally active domains after MYC overexpression. These results demonstrate the accumulation of overexpressed MYC at promoter–enhancer hubs and suggest that MYC invades into enhancers through spatial proximity. At the same time, the increased protein–protein interactions may strengthen these chromatin interactions to increase chromatin contact frequency. CTCF siRNA knockdown in MYC-overexpressed U2OS cells demonstrated that removal of architectural proteins can disperse MYC and abrogate the increase in chromatin contacts. By elucidating the chromatin landscape of MYC-driven cancers, we can potentially target MYC-associated chromatin interactions for cancer therapy.

Super Enhancer-Mediated Upregulation of HJURP Promotes Growth and Survival of t(4;14)-Positive Multiple Myeloma

Multiple myeloma is an incurable malignancy with marked clinical and genetic heterogeneity. The cytogenetic abnormality t(4;14) (p16.3;q32.3) confers aggressive behavior in multiple myeloma. Recently, essential oncogenic drivers in a wide range of cancers have been shown to be controlled by super-enhancers (SE). We used chromatin immunoprecipitation sequencing of the active enhancer marker histone H3 lysine 27 acetylation (H3K27ac) to profile unique SEs in t(4;14)-translocated multiple myeloma. The histone chaperone HJURP was aberrantly overexpressed in t(4;14)-positive multiple myeloma due to transcriptional activation by a distal SE induced by the histone lysine methyltransferase NSD2. Silencing of HJURP with short hairpin RNA or CRISPR interference of SE function impaired cell viability and led to apoptosis. Conversely, HJURP overexpression promoted cell proliferation and abrogated apoptosis. Mechanistically, the NSD2/BRD4 complex positively coregulated HJURP transcription by binding the promoter and active elements of its SE. In summary, this study introduces SE profiling as an efficient approach to identify new targets and understand molecular pathogenesis in specific subtypes of cancer. Moreover, HJURP could be a valuable therapeutic target in patients with t(4;14)-positive myeloma. SIGNIFICANCE: A super-enhancer screen in t(4;14) multiple myeloma serves to identify genes that promote growth and survival of myeloma cells, which may be evaluated in future studies as therapeutic targets.

Chromatin interaction neural network (ChINN): a machine learning-based method for predicting chromatin interactions from DNA sequences

Chromatin interactions play important roles in regulating gene expression. However, the availability of genome-wide chromatin interaction data is limited. We develop a computational method, chromatin interaction neural network (ChINN), to predict chromatin interactions between open chromatin regions using only DNA sequences. ChINN predicts CTCF- and RNA polymerase II-associated and Hi-C chromatin interactions. ChINN shows good across-sample performances and captures various sequence features for chromatin interaction prediction. We apply ChINN to 6 chronic lymphocytic leukemia (CLL) patient samples and a published cohort of 84 CLL open chromatin samples. Our results demonstrate extensive heterogeneity in chromatin interactions among CLL patient samples.

Interplay and cooperation between SREBF1 and master transcription factors regulate lipid metabolism and tumor-promoting pathways in squamous cancer

Squamous cell carcinomas (SCCs) comprise one of the most common histologic types of human cancer. Transcriptional dysregulation of SCC cells is orchestrated by tumor protein p63 (TP63), a master transcription factor (TF) and a well-researched SCC-specific oncogene. In the present study, both Gene Set Enrichment Analysis (GSEA) of SCC patient samples and in vitro loss-of-function assays establish fatty-acid metabolism as a key pathway downstream of TP63. Further studies identify sterol regulatory element binding transcription factor 1 (SREBF1) as a central mediator linking TP63 with fatty-acid metabolism, which regulates the biosynthesis of fatty-acids, sphingolipids (SL), and glycerophospholipids (GPL), as revealed by liquid chromatography tandem mass spectrometry (LC-MS/MS)-based lipidomics. Moreover, a feedback co-regulatory loop consisting of SREBF1/TP63/Kruppel like factor 5 (KLF5) is identified, which promotes overexpression of all three TFs in SCCs. Downstream of SREBF1, a non-canonical, SCC-specific function is elucidated: SREBF1 cooperates with TP63/KLF5 to regulate hundreds of cis-regulatory elements across the SCC epigenome, which converge on activating cancer-promoting pathways. Indeed, SREBF1 is essential for SCC viability and migration, and its overexpression is associated with poor survival in SCC patients. Taken together, these data shed light on mechanisms of transcriptional dysregulation in cancer, identify specific epigenetic regulators of lipid metabolism, and uncover SREBF1 as a potential therapeutic target and prognostic marker in SCC.

Abstract LB189: Three-dimensional genome organization maps in normal haematopoietic stem cells and acute myeloid leukemia

Many oncogenes, such as MEIS1 and HOXA9, are overexpressed in some but not all cancers. We identified two key epigenetic mechanisms underlying this heterogeneity in oncogene expression in Acute Myeloid Leukemia. Acute Myeloid Leukemia (AML) is a highly lethal blood cancer arising due to aberrant differentiation of haematopoietic stem cells. MEIS1 and HOXA9 regulate stemness-related transcriptional programs in normal haematopoietic stem cells and AML. Here we obtained 3D genome organization maps in the CD34+ haematopoietic stem cells from healthy individuals and individuals with AML. The MEIS1 oncogenic transcription factor is regulated by a Frequently Interacting Region (FIRE). This FIRE is present in normal bone marrow samples, and an AML sample with high MEIS1 levels. However, it is absent in two AML samples that show low MEIS1 levels. CRISPR excision of the FIRE led to loss of MEIS1 and reduced cell growth. Moreover, MEIS1 can bind to the promoter of HOXA9. HOXA9 can also auto-regulate by binding to its own promoter as well as an Acute Myeloid Leukemia-specific super-enhancer that interacts with the HOXA9 promoter via chromatin interactions. The mechanisms elucidated here could be potentially exploited to utilize epigenetic inhibitors to specifically target oncogene expression in cancer.

Citation Format: Benny Wang, Lingshi Kong, Deepak Babu, Ruchi Choudhary, Winnie Fam, Jia Qi Tng, Yufen Goh, Xin Liu, Fang Fang Song, Priscella Chia, Ming Chun Chan, Omer An, Cheng Yong Tham, Touati Benoukraf, Henry Yang, Wilson Wang, Wee Joo Chng, Daniel Tenen, Melissa J. Fullwood. Three-dimensional genome organization maps in normal haematopoietic stem cells and acute myeloid leukemia [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr LB189.

Defining Essential Enhancers for Pluripotent Stem Cells Using a Features-Oriented CRISPR-Cas9 Screen

cis-regulatory elements (CREs) regulate the expression of genes in their genomic neighborhoods and influence cellular processes such as cell-fate maintenance and differentiation. To date, there remain major gaps in the functional characterization of CREs and the identification of their target genes in the cellular native environment. In this study, we perform a features-oriented CRISPR-utilized systematic (FOCUS) screen of OCT4-bound CREs using CRISPR-Cas9 to identify functional enhancers important for pluripotency maintenance in mESCs. From the initial 235 candidates tested, 16 CREs are identified to be essential stem cell enhancers. Using RNA-seq and genomic 4C-seq, we further uncover a complex network of candidate CREs and their downstream target genes, which supports the growth and self-renewal of mESCs. Notably, an essential enhancer, CRE111, and its target, Lrrc31, form the important switch to modulate the LIF-JAK1-STAT3 signaling pathway.

TP63, SOX2, and KLF5 Establish a Core Regulatory Circuitry That Controls Epigenetic and Transcription Patterns in Esophageal Squamous Cell Carcinoma Cell Lines

BACKGROUND & AIMS

We investigated the transcriptome of esophageal squamous cell carcinoma (ESCC) cells, activity of gene regulatory (enhancer and promoter regions), and the effects of blocking epigenetic regulatory proteins.

METHODS

We performed chromatin immunoprecipitation sequencing with antibodies against H3K4me1, H3K4me3, and H3K27ac and an assay for transposase-accessible chromatin to map the enhancer regions and accessible chromatin in 8 ESCC cell lines. We used the CRC_Mapper algorithm to identify core regulatory circuitry transcription factors in ESCC cell lines, and determined genome occupancy profiles for 3 of these factors. In ESCC cell lines, expression of transcription factors was knocked down with small hairpin RNAs, promoter and enhancer regions were disrupted by CRISPR/Cas9 genome editing, or bromodomains and extraterminal (BET) family proteins and histone deacetylases (HDACs) were inhibited with ARV-771 and romidepsin, respectively. ESCC cell lines were then analyzed by whole-transcriptome sequencing, immunoprecipitation, immunoblots, immunohistochemistry, and viability assays. Interactions between distal enhancers and promoters were identified and verified with circular chromosome conformation capture sequencing. NOD-SCID mice were given injections of modified ESCC cells, some mice where given injections of HDAC or BET inhibitors, and growth of xenograft tumors was measured.

RESULTS

We identified super-enhancer-regulated circuits and transcription factors TP63, SOX2, and KLF5 as core regulatory factors in ESCC cells. Super-enhancer regulation of ALDH3A1 mediated by core regulatory factors was required for ESCC viability. We observed direct interactions between the promoter region of TP63 and functional enhancers, mediated by the core regulatory circuitry transcription factors. Deletion of enhancer regions from ESCC cells decreased expression of the core regulatory circuitry transcription factors and reduced cell viability; these same results were observed with knockdown of each core regulatory circuitry transcription factor. Incubation of ESCC cells with BET and HDAC disrupted the core regulatory circuitry program and the epigenetic modifications observed in these cells; mice given injections of HDAC or BET inhibitors developed smaller xenograft tumors from the ESCC cell lines. Xenograft tumors grew more slowly in mice given the combination of ARV-771 and romidepsin than mice given either agent alone.

CONCLUSIONS

In epigenetic and transcriptional analyses of ESCC cell lines, we found the transcription factors TP63, SOX2, and KLF5 to be part of a core regulatory network that determines chromatin accessibility, epigenetic modifications, and gene expression patterns in these cells. A combination of epigenetic inhibitors slowed growth of xenograft tumors derived from ESCC cells in mice.

NanoVar: accurate characterization of patients' genomic structural variants using low-depth nanopore sequencing

The recent advent of third-generation sequencing technologies brings promise for better characterization of genomic structural variants by virtue of having longer reads. However, long-read applications are still constrained by their high sequencing error rates and low sequencing throughput. Here, we present NanoVar, an optimized structural variant caller utilizing low-depth (8X) whole-genome sequencing data generated by Oxford Nanopore Technologies. NanoVar exhibits higher structural variant calling accuracy when benchmarked against current tools using low-depth simulated datasets. In patient samples, we successfully validate structural variants characterized by NanoVar and uncover normal alternative sequences or alleles which are present in healthy individuals.

Super-Enhancers and Broad H3K4me3 Domains Form Complex Gene Regulatory Circuits Involving Chromatin Interactions

Stretched histone regions, such as super-enhancers and broad H3K4me3 domains, are associated with maintenance of cell identity and cancer. We connected super-enhancers and broad H3K4me3 domains in the K562 chronic myelogenous leukemia cell line as well as the MCF-7 breast cancer cell line with chromatin interactions. Super-enhancers and broad H3K4me3 domains showed higher association with chromatin interactions than their typical counterparts. Interestingly, we identified a subset of super-enhancers that overlap with broad H3K4me3 domains and show high association with cancer-associated genes including tumor suppressor genes. Besides cell lines, we could observe chromatin interactions by a Chromosome Conformation Capture (3C)-based method, in primary human samples. Several chromatin interactions involving super-enhancers and broad H3K4me3 domains are constitutive and can be found in both cancer and normal samples. Taken together, these results reveal a new layer of complexity in gene regulation by super-enhancers and broad H3K4me3 domains.

Deciphering Noncoding RNA and Chromatin Interactions: Multiplex Chromatin Interaction Analysis by Paired-End Tag Sequencing (mChIA-PET)

Genomic DNA is dynamically associated with protein factors and folded to form chromatin fibers. The 3-dimensional (3D) configuration of the chromatin will enable the distal genetic elements to come into close proximity, allowing transcriptional regulation. Noncoding RNA can mediate the 3D structure of chromatin. Chromatin Interaction Analysis by Paired-End Tag Sequencing (ChIA-PET) is a valuable and powerful technique in molecular biology which allows the study of unbiased, genome-wide de novo chromatin interactions with paired-end tags. Here, we describe the standard version of ChIA-PET and a Multiplex ChIA-PET version.

3D genome organization in health and disease: emerging opportunities in cancer translational medicine

Organizing the DNA to fit inside a spatially constrained nucleus is a challenging problem that has attracted the attention of scientists across all disciplines of science. Increasing evidence has demonstrated the importance of genome geometry in several cellular contexts that affect human health. Among several approaches, the application of sequencing technologies has substantially increased our understanding of this intricate organization, also known as chromatin interactions. These structures are involved in transcriptional control of gene expression by connecting distal regulatory elements with their target genes and regulating co-transcriptional splicing. In addition, chromatin interactions play pivotal roles in the organization of the genome, the formation of structural variants, recombination, DNA replication and cell division. Mutations in factors that regulate chromatin interactions lead to the development of pathological conditions, for example, cancer. In this review, we discuss key findings that have shed light on the importance of these structures in the context of cancers, and highlight the applicability of chromatin interactions as potential biomarkers in molecular medicine as well as therapeutic implications of chromatin interactions.

Roles, Functions, and Mechanisms of Long Non-coding RNAs in Cancer

Long non-coding RNAs (lncRNAs) play important roles in cancer. They are involved in chromatin remodeling, as well as transcriptional and post-transcriptional regulation, through a variety of chromatin-based mechanisms and via cross-talk with other RNA species. lncRNAs can function as decoys, scaffolds, and enhancer RNAs. This review summarizes the characteristics of lncRNAs, including their roles, functions, and working mechanisms, describes methods for identifying and annotating lncRNAs, and discusses future opportunities for lncRNA-based therapies using antisense oligonucleotides.

A resource of ribosomal RNA-depleted RNA-Seq data from different normal adult and fetal human tissues

Gene expression is the most fundamental level at which the genotype leads to the phenotype of the organism. Enabled by ultra-high-throughput next-generation DNA sequencing, RNA-Seq involves shotgun sequencing of fragmented RNA transcripts by next-generation sequencing followed by in silico assembly, and is rapidly becoming the most popular method for gene expression analysis. Poly[A]+ RNA-Seq analyses of normal human adult tissue samples such as Illumina's Human BodyMap 2.0 Project and the RNA-Seq atlas have provided a useful global resource and framework for comparisons with diseased tissues such as cancer. However, these analyses have failed to provide information on poly[A]-RNA, which is abundant in our cells. The most recent advances in RNA-Seq analyses use ribosomal RNA-depletion to provide information on both poly[A]+ and poly[A]-RNA. In this paper, we describe the use of Illumina's HiSeq 2000 to generate high quality rRNA-depleted RNA-Seq datasets from human fetal and adult tissues. The datasets reported here will be useful in understanding the different expression profiles in different tissues.

Landscape of transcription in human cells

Eukaryotic cells make many types of primary and processed RNAs that are found either in specific subcellular compartments or throughout the cells. A complete catalogue of these RNAs is not yet available and their characteristic subcellular localizations are also poorly understood. Because RNA represents the direct output of the genetic information encoded by genomes and a significant proportion of a cell's regulatory capabilities are focused on its synthesis, processing, transport, modification and translation, the generation of such a catalogue is crucial for understanding genome function. Here we report evidence that three-quarters of the human genome is capable of being transcribed, as well as observations about the range and levels of expression, localization, processing fates, regulatory regions and modifications of almost all currently annotated and thousands of previously unannotated RNAs. These observations, taken together, prompt a redefinition of the concept of a gene.

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.

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 ¹³. 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 ⁸.

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.

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.

Next-generation sequencing of apoptotic DNA breakpoints reveals association with actively transcribed genes and gene translocations

DNA fragmentation is a well-recognized hallmark of apoptosis. However, the precise DNA sequences cleaved during apoptosis triggered by distinct mechanisms remain unclear. We used next-generation sequencing of DNA fragments generated in Actinomycin D-treated human HL-60 leukemic cells to generate a high-throughput, global map of apoptotic DNA breakpoints. These data highlighted that DNA breaks are non-random and show a significant association with active genes and open chromatin regions. We noted that transcription factor binding sites were also enriched within a fraction of the apoptotic breakpoints. Interestingly, extensive apoptotic cleavage was noted within genes that are frequently translocated in human cancers. We speculate that the non-random fragmentation of DNA during apoptosis may contribute to gene translocations and the development of human cancers.

ChIA-PET tool for comprehensive chromatin interaction analysis with paired-end tag sequencing

ChIA-PET Tool can be used to process long-range chromatin interaction data. Results are visualized on a user-friendly genome browser.

Chromatin interaction analysis with paired-end tag sequencing (ChIA-PET) is a new technology to study genome-wide long-range chromatin interactions bound by protein factors. Here we present ChIA-PET Tool, a software package for automatic processing of ChIA-PET sequence data, including linker filtering, mapping tags to reference genomes, identifying protein binding sites and chromatin interactions, and displaying the results on a graphical genome browser. ChIA-PET Tool is fast, accurate, comprehensive, user-friendly, and open source (available at http://chiapet.gis.a-star.edu.sg).

Chromatin interaction analysis using paired-end tag sequencing

Chromatin Interaction Analysis using Paired-End Tag sequencing (ChIA-PET) is a technique developed for large-scale, de novo analysis of higher-order chromatin structures. Cells are treated with formaldehyde to cross-link chromatin interactions, DNA segments bound by protein factors are enriched by chromatin immunoprecipitation, and interacting DNA fragments are then captured by proximity ligation. The Paired-End Tag (PET) strategy is applied to the construction of ChIA-PET libraries, which are sequenced by high-throughput next-generation sequencing technologies. Finally, raw PET sequences are subjected to bioinformatics analysis, resulting in a genome-wide map of binding sites and chromatin interactions mediated by the protein factor under study. This unit describes ChIA-PET for genome-wide analysis of chromatin interactions in mammalian cells, with the application of Roche/454 and Illumina sequencing technologies.

An oestrogen-receptor-alpha-bound human chromatin interactome

Genomes are organized into high-level three-dimensional structures, and DNA elements separated by long genomic distances can in principle interact functionally. Many transcription factors bind to regulatory DNA elements distant from gene promoters. Although distal binding sites have been shown to regulate transcription by long-range chromatin interactions at a few loci, chromatin interactions and their impact on transcription regulation have not been investigated in a genome-wide manner. Here we describe the development of a new strategy, chromatin interaction analysis by paired-end tag sequencing (ChIA-PET) for the de novo detection of global chromatin interactions, with which we have comprehensively mapped the chromatin interaction network bound by oestrogen receptor alpha (ER-alpha) in the human genome. We found that most high-confidence remote ER-alpha-binding sites are anchored at gene promoters through long-range chromatin interactions, suggesting that ER-alpha functions by extensive chromatin looping to bring genes together for coordinated transcriptional regulation. We propose that chromatin interactions constitute a primary mechanism for regulating transcription in mammalian genomes.

ChIP-based methods for the identification of long-range chromatin interactions

Chromatin immunoprecipitation (ChIP) is an important technique for studying protein-DNA interactions. Whole genome ChIP methods have enjoyed much success, but are limited in that they cannot uncover important long-range chromatin interactions. Chromosome conformation capture (3C) and related methods are capable of detecting remote chromatin interactions, but are tedious, have low signal-to-noise ratios, and are not genome-wide. Although the addition of ChIP to 3C (ChIP-3C) would conceivably reduce noise and increase specificity for chromatin interaction detection, there are concerns that simple mixing of the ChIP and 3C protocols would lead to high levels of false positives. In this essay, we dissect current ChIP- and 3C-based methodologies, discuss the models of specific as opposed to non-specific chromatin interactions, and suggest approaches to separate specific chromatin complexes from non-specific chromatin fragments. We conclude that the combination of sonication-based chromatin fragmentation, ChIP-based enrichment, chromatin proximity ligation and Paired-End Tag ultra-high-throughput sequencing will be a winning implementation for genome-wide, unbiased and de novo discovery of long-range chromatin interactions, which will help to establish an emerging field for studying human chromatin interactomes and genome regulation networks in three-dimensional spaces.

Next-generation DNA sequencing of paired-end tags (PET) for transcriptome and genome analyses

Comprehensive understanding of functional elements in the human genome will require thorough interrogation and comparison of individual human genomes and genomic structures. Such an endeavor will require improvements in the throughputs and costs of DNA sequencing. Next-generation sequencing platforms have impressively low costs and high throughputs but are limited by short read lengths. An immediate and widely recognized solution to this critical limitation is the paired-end tag (PET) sequencing for various applications, collectively called the PET sequencing strategy, in which short and paired tags are extracted from the ends of long DNA fragments for ultra-high-throughput sequencing. The PET sequences can be accurately mapped to the reference genome, thus demarcating the genomic boundaries of PET-represented DNA fragments and revealing the identities of the target DNA elements. PET protocols have been developed for the analyses of transcriptomes, transcription factor binding sites, epigenetic sites such as histone modification sites, and genome structures. The exclusive advantage of the PET technology is its ability to uncover linkages between the two ends of DNA fragments. Using this unique feature, unconventional fusion transcripts, genome structural variations, and even molecular interactions between distant genomic elements can be unraveled by PET analysis. Extensive use of PET data could lead to efficient assembly of individual human genomes, transcriptomes, and interactomes, enabling new biological and clinical insights. With its versatile and powerful nature for DNA analysis, the PET sequencing strategy has a bright future ahead.

The use of multiple displacement amplification to amplify complex DNA libraries

Complex libraries for genomic DNA and cDNA sequencing analyses are typically amplified using bacterial propagation. To reduce biases, large numbers of colonies are plated and scraped from solid-surface agar. This process is time consuming, tedious and limits scaling up. At the same time, multiple displacement amplification (MDA) has been recently developed as a method for in vitro amplification of DNA. However, MDA has no selection function for the removal of ligation multimers. We developed a novel method of briefly introducing ligation reactions into bacteria to select single insert DNA clones followed by MDA to amplify. We applied these methods to a Gene Identification Signatures with Paired-End diTags (GIS-PET) library, which is a complex transcriptome library created by pairing short tags from the 5' and 3' ends of cDNA fragments together, and demonstrated that this selection and amplification strategy is unbiased and efficient.

Multiplex sequencing of paired-end ditags (MS-PET): a strategy for the ultra-high-throughput analysis of transcriptomes and genomes

The paired-end ditagging (PET) technique has been shown to be efficient and accurate for large-scale transcriptome and genome analysis. However, as with other DNA tag-based sequencing strategies, it is constrained by the current efficiency of Sanger technology. A recently developed multiplex sequencing method (454-sequencing™) using picolitre-scale reactions has achieved a remarkable advance in efficiency, but suffers from short-read lengths, and a lack of paired-end information. To further enhance the efficiency of PET analysis and at the same time overcome the drawbacks of the new sequencing method, we coupled multiplex sequencing with paired-end ditagging (MS-PET) using modified PET procedures to simultaneously sequence 200 000 to 300 000 dimerized PET (diPET) templates, with an output of nearly half-a-million PET sequences in a single 4 h machine run. We demonstrate the utility and robustness of MS-PET by analyzing the transcriptome of human breast carcinoma cells, and by mapping p53 binding sites in the genome of human colorectal carcinoma cells. This combined sequencing strategy achieved an approximate 100-fold efficiency increase over the current standard for PET analysis, and furthermore enables the short-read-length multiplex sequencing procedure to acquire paired-end information from large DNA fragments.