In situ Chromatin Interaction Analysis Using Paired-End Tag Sequencing.

Chromatin Interaction Analysis Using Paired-End Tag Sequencing (ChIA-PET) is an established method to map protein-mediated chromatin interactions. A limitation, however, is that it requires a hundred million cells per experiment, which hampers its broad application in biomedical research, particularly in studies in which it is impractical to obtain a large number of cells from rare samples. To reduce the required input cell number while retaining high data quality, we developed an in situ ChIA-PET protocol, which requires as few as 1 million cells. Here, we describe detailed step-by-step procedures for performing in situ ChIA-PET from cultured cells, including both an experimental protocol for sample preparation and data generation and a computational protocol for data processing and visualization using the ChIA-PIPE pipeline. As the protocol significantly simplifies the experimental procedure, reduces ligation noise, and decreases the required input of cells compared to previous versions of ChIA-PET protocols, it can be applied to generate high-resolution chromatin contact maps mediated by various protein factors for a wide range of human and mouse primary cells. © 2021 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Sample preparation and data generation Support Protocol: Bridge linker preparation Basic Protocol 2: Data processing and visualization.

Sequencing of E. coli strain UTI89 on multiple sequencing platforms.

OBJECTIVES: The availability of matched sequencing data for the same sample across different sequencing platforms is a necessity for validation and effective comparison of sequencing platforms. A commonly sequenced sample is the lab-adapted MG1655 strain of Escherichia coli; however, this strain is not fully representative of more complex and dynamic genomes of pathogenic E. coli strains. DATA DESCRIPTION: We present six new sequencing data sets for another E. coli strain, UTI89, which is an extraintestinal pathogenic strain isolated from a patient suffering from a urinary tract infection. We now provide matched whole genome sequencing data generated using the PacBio RSII, Oxford Nanopore MinION R9.4, Ion Torrent, ABI SOLiD, and Illumina NextSeq sequencers. Together with other publically available datasets, UTI89 has a nearly complete suite of data generated on most second- and third-generation sequencers. These data can be used as an additional validation set for new sequencing technologies and analytical methods. More than being another E. coli strain, however, UTI89 is pathogenic, with a 10% larger genome, additional pathogenicity islands, and a large plasmid, features that are common among other naturally occurring and disease-causing E. coli isolates. These data therefore provide a more medically relevant test set for development of algorithms.

Chromatin topology reorganization and transcription repression by PML-RARα in acute promyeloid leukemia.

BACKGROUND: Acute promyeloid leukemia (APL) is characterized by the oncogenic fusion protein PML-RARα, a major etiological agent in APL. However, the molecular mechanisms underlying the role of PML-RARα in leukemogenesis remain largely unknown. RESULTS: Using an inducible system, we comprehensively analyze the 3D genome organization in myeloid cells and its reorganization after PML-RARα induction and perform additional analyses in patient-derived APL cells with native PML-RARα. We discover that PML-RARα mediates extensive chromatin interactions genome-wide. Globally, it redefines the chromatin topology of the myeloid genome toward a more condensed configuration in APL cells; locally, it intrudes RNAPII-associated interaction domains, interrupts myeloid-specific transcription factors binding at enhancers and super-enhancers, and leads to transcriptional repression of genes critical for myeloid differentiation and maturation. CONCLUSIONS: Our results not only provide novel topological insights for the roles of PML-RARα in transforming myeloid cells into leukemia cells, but further uncover a topological framework of a molecular mechanism for oncogenic fusion proteins in cancers.

Coordinated demethylation of H3K9 and H3K27 is required for rapid inflammatory responses of endothelial cells.

Histone H3 lysine-9 di-methylation (H3K9me2) and lysine-27 tri-methylation (H3K27me3) are linked to repression of gene expression, but the functions of repressive histone methylation dynamics during inflammatory responses remain enigmatic. Here, we report that lysine demethylases 7A (KDM7A) and 6A (UTX) play crucial roles in tumor necrosis factor (TNF)-α signaling in endothelial cells (ECs), where they are regulated by a novel TNF-α-responsive microRNA, miR-3679-5p. TNF-α rapidly induces co-occupancy of KDM7A and UTX at nuclear factor kappa-B (NF-κB)-associated elements in human ECs. KDM7A and UTX demethylate H3K9me2 and H3K27me3, respectively, and are both required for activation of NF-κB-dependent inflammatory genes. Chromosome conformation capture-based methods furthermore uncover increased interactions between TNF-α-induced super enhancers at NF-κB-relevant loci, coinciding with KDM7A and UTX recruitments. Simultaneous pharmacological inhibition of KDM7A and UTX significantly reduces leukocyte adhesion in mice, establishing the biological and potential translational relevance of this mechanism. Collectively, these findings suggest that rapid erasure of repressive histone marks by KDM7A and UTX is essential for NF-κB-dependent regulation of genes that control inflammatory responses of ECs.

The Genome of C57BL/6J "Eve", the Mother of the Laboratory Mouse Genome Reference Strain.

Isogenic laboratory mouse strains enhance reproducibility because individual animals are genetically identical. For the most widely used isogenic strain, C57BL/6, there exists a wealth of genetic, phenotypic, and genomic data, including a high-quality reference genome (GRCm38.p6). Now 20 years after the first release of the mouse reference genome, C57BL/6J mice are at least 26 inbreeding generations removed from GRCm38 and the strain is now maintained with periodic reintroduction of cryorecovered mice derived from a single breeder pair, aptly named Adam and Eve. To provide an update to the mouse reference genome that more accurately represents the genome of today's C57BL/6J mice, we took advantage of long read, short read, and optical mapping technologies to generate a de novo assembly of the C57BL/6J Eve genome (B6Eve). Using these data, we have addressed recurring variants observed in previous mouse genomic studies. We have also identified structural variations, closed gaps in the mouse reference assembly, and revealed previously unannotated coding sequences. This B6Eve assembly explains discrepant observations that have been associated with GRCm38-based analyses, and will inform a reference genome that is more representative of the C57BL/6J mice that are in use today.

Multiplex chromatin interactions with single-molecule precision.

The genomes of multicellular organisms are extensively folded into 3D chromosome territories within the nucleus1. Advanced 3D genome-mapping methods that combine proximity ligation and high-throughput sequencing (such as chromosome conformation capture, Hi-C)2, and chromatin immunoprecipitation techniques (such as chromatin interaction analysis by paired-end tag sequencing, ChIA-PET)3, have revealed topologically associating domains4 with frequent chromatin contacts, and have identified chromatin loops mediated by specific protein factors for insulation and regulation of transcription5-7. However, these methods rely on pairwise proximity ligation and reflect population-level views, and thus cannot reveal the detailed nature of chromatin interactions. Although single-cell Hi-C8 potentially overcomes this issue, this method may be limited by the sparsity of data that is inherent to current single-cell assays. Recent advances in microfluidics have opened opportunities for droplet-based genomic analysis9 but this approach has not yet been adapted for chromatin interaction analysis. Here we describe a strategy for multiplex chromatin-interaction analysis via droplet-based and barcode-linked sequencing, which we name ChIA-Drop. We demonstrate the robustness of ChIA-Drop in capturing complex chromatin interactions with single-molecule precision, which has not been possible using methods based on population-level pairwise contacts. By applying ChIA-Drop to Drosophila cells, we show that chromatin topological structures predominantly consist of multiplex chromatin interactions with high heterogeneity; ChIA-Drop also reveals promoter-centred multivalent interactions, which provide topological insights into transcription.

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.

Chromatin Immunoprecipitation for Chromatin Interaction Analysis Using Paired-End-Tag (ChIA-PET) Sequencing in Tadpole Tissues.

Proper gene expression involves communication between the regulatory elements and promoters of genes. Because regulatory elements can be located over a large range of genomic distances (from as close as a few hundred bp to as much as several Mb away), contact and communication between regulators and the core transcriptional machinery at promoters are mediated through DNA looping. Today, chromosome conformation capture (3C)-based methods efficiently probe chromosome folding in the nucleus and thus provide a molecular description of physical proximity between enhancer(s) and their target promoter(s). One such method, chromatin interaction analysis using paired-end-tag (ChIA-PET) sequencing, is a leading high-throughput method for detection of genome wide chromatin interactions. Briefly, the method involves cross-linkage of chromatin (-DNA) fibers in cells in situ, fragmentation of the fixed chromatin-DNA complexes by sonication, followed by enrichment of the chromatin complexes with a dedicated antibody through the process of immunoprecipitation (IP). Next, application of the ChIA-PET protocol followed by deep sequencing and mapping of reads to the reference genome reveals both binding sites and remote chromatin interactions mediated by the protein factors of interest. The method detailed here focuses on ChIP sample preparation and can be completed in ∼5 d. The ChIA-PET method is detailed in an associated protocol. Because not all chromatin immunoprecipitation protocols are suitable for ChIA-PET, it is important to strictly follow this procedure before performing the ChIA-PET protocol.

Chromatin Interaction Analysis Using Paired-End-Tag (ChIA-PET) Sequencing in Tadpole Tissues.

Proper gene expression involves communication between the regulatory elements and promoters of genes. Today, chromosome conformation capture (3C)-based methods efficiently probe chromosome folding in the nucleus and thus provide a molecular description of physical proximity through DNA looping between enhancer(s) and their target promoter(s). One such method, chromatin interaction analysis using paired-end-tag (ChIA-PET) sequencing is a powerful high-throughput method for detection of genome-wide chromatin interactions. Following enrichment of the chromatin complexes with a dedicated antibody, through a process of immunoprecipitation (IP), DNA fragments are end-joined with specifically designed DNA-linkers through proximity ligation. The DNA-linkers contain the binding site for the type II restriction enzyme MmeI, which cleaves 20 bp from each end of the ligated fragments, thus releasing a "paired end tag" (PET): [20 bp tag]-[linker]-[20 bp tag]. The PETs are then deep-sequenced and reads are mapped to the reference genome, revealing both binding sites, as well as remote chromatin interactions mediated by the protein factors of interest. The method detailed here focuses on ChIA-PET library construction and can be completed in 2 wk.

Mammalian Systems Biotechnology Reveals Global Cellular Adaptations in a Recombinant CHO Cell Line.

Effective development of host cells for therapeutic protein production is hampered by the poor characterization of cellular transfection. Here, we employed a multi-omics-based systems biotechnology approach to elucidate the genotypic and phenotypic differences between a wild-type and recombinant antibody-producing Chinese hamster ovary (CHO) cell line. At the genomic level, we observed extensive rearrangements in specific targeted loci linked to transgene integration sites. Transcriptional re-wiring of DNA damage repair and cellular metabolism in the antibody producer, via changes in gene copy numbers, was also detected. Subsequent integration of transcriptomic data with a genome-scale metabolic model showed a substantial increase in energy metabolism in the antibody producer. Metabolomics, lipidomics, and glycomics analyses revealed an elevation in long-chain lipid species, potentially associated with protein transport and secretion requirements, and a surprising stability of N-glycosylation profiles between both cell lines. Overall, the proposed knowledge-based systems biotechnology framework can further accelerate mammalian cell-line engineering in a targeted manner.

CTCF-Mediated Human 3D Genome Architecture Reveals Chromatin Topology for Transcription.

Spatial genome organization and its effect on transcription remains a fundamental question. We applied an advanced chromatin interaction analysis by paired-end tag sequencing (ChIA-PET) strategy to comprehensively map higher-order chromosome folding and specific chromatin interactions mediated by CCCTC-binding factor (CTCF) and RNA polymerase II (RNAPII) with haplotype specificity and nucleotide resolution in different human cell lineages. We find that CTCF/cohesin-mediated interaction anchors serve as structural foci for spatial organization of constitutive genes concordant with CTCF-motif orientation, whereas RNAPII interacts within these structures by selectively drawing cell-type-specific genes toward CTCF foci for coordinated transcription. Furthermore, we show that haplotype variants and allelic interactions have differential effects on chromosome configuration, influencing gene expression, and may provide mechanistic insights into functions associated with disease susceptibility. 3D genome simulation suggests a model of chromatin folding around chromosomal axes, where CTCF is involved in defining the interface between condensed and open compartments for structural regulation. Our 3D genome strategy thus provides unique insights in the topological mechanism of human variations and diseases.

The genome sequence of Sea-Island cotton (Gossypium barbadense) provides insights into the allopolyploidization and development of superior spinnable fibres.

Gossypium hirsutum contributes the most production of cotton fibre, but G. barbadense is valued for its better comprehensive resistance and superior fibre properties. However, the allotetraploid genome of G. barbadense has not been comprehensively analysed. Here we present a high-quality assembly of the 2.57 gigabase genome of G. barbadense, including 80,876 protein-coding genes. The double-sized genome of the A (or At) (1.50 Gb) against D (or Dt) (853 Mb) primarily resulted from the expansion of Gypsy elements, including Peabody and Retrosat2 subclades in the Del clade, and the Athila subclade in the Athila/Tat clade. Substantial gene expansion and contraction were observed and rich homoeologous gene pairs with biased expression patterns were identified, suggesting abundant gene sub-functionalization occurred by allopolyploidization. More specifically, the CesA gene family has adapted differentially temporal expression patterns, suggesting an integrated regulatory mechanism of CesA genes from At and Dt subgenomes for the primary and secondary cellulose biosynthesis of cotton fibre in a "relay race"-like fashion. We anticipate that the G. barbadense genome sequence will advance our understanding the mechanism of genome polyploidization and underpin genome-wide comparison research in this genus.

Xenopus tropicalis Genome Re-Scaffolding and Re-Annotation Reach the Resolution Required for In Vivo ChIA-PET Analysis.

Genome-wide functional analyses require high-resolution genome assembly and annotation. We applied ChIA-PET to analyze gene regulatory networks, including 3D chromosome interactions, underlying thyroid hormone (TH) signaling in the frog Xenopus tropicalis. As the available versions of Xenopus tropicalis assembly and annotation lacked the resolution required for ChIA-PET we improve the genome assembly version 4.1 and annotations using data derived from the paired end tag (PET) sequencing technologies and approaches (e.g., DNA-PET [gPET], RNA-PET etc.). The large insert (~10 Kb, ~17 Kb) paired end DNA-PET with high throughput NGS sequencing not only significantly improved genome assembly quality, but also strongly reduced genome "fragmentation", reducing total scaffold numbers by ~60%. Next, RNA-PET technology, designed and developed for the detection of full-length transcripts and fusion mRNA in whole transcriptome studies (ENCODE consortia), was applied to capture the 5' and 3' ends of transcripts. These amendments in assembly and annotation were essential prerequisites for the ChIA-PET analysis of TH transcription regulation. Their application revealed complex regulatory configurations of target genes and the structures of the regulatory networks underlying physiological responses. Our work allowed us to improve the quality of Xenopus tropicalis genomic resources, reaching the standard required for ChIA-PET analysis of transcriptional networks. We consider that the workflow proposed offers useful conceptual and methodological guidance and can readily be applied to other non-conventional models that have low-resolution genome data.

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.

Coherence analysis discriminates between retroviral integration patterns in CD34(+) cells transduced under differing clinical trial conditions.

Unequivocal demonstration of the therapeutic utility of γ-retroviral vectors for gene therapy applications targeting the hematopoietic system was accompanied by instances of insertional mutagenesis. These events stimulated the ongoing development of putatively safer integrating vector systems and analysis methods to characterize and compare integration site (IS) biosafety profiles. Continuing advances in next-generation sequencing technologies are driving the generation of ever-more complex IS datasets. Available bioinformatic tools to compare such datasets focus on the association of integration sites (ISs) with selected genomic and epigenetic features, and the choice of these features determines the ability to discriminate between datasets. We describe the scalable application of point-process coherence analysis (CA) to compare patterns produced by vector ISs across genomic intervals, uncoupled from association with genomic features. To explore the utility of CA in the context of an unresolved question, we asked whether the differing transduction conditions used in the initial Paris and London SCID-X1 gene therapy trials result in divergent genome-wide integration profiles. We tested a transduction carried out under each condition, and showed that CA could indeed resolve differences in IS distributions. Existence of these differences was confirmed by the application of established methods to compare integration datasets.

Deciphering the regulatory logic of an ancient, ultraconserved nuclear receptor enhancer module.

Cooperative, synergistic gene regulation by nuclear hormone receptors can increase sensitivity and amplify cellular responses to hormones. We investigated thyroid hormone (TH) and glucocorticoid (GC) synergy on the Krüppel-like factor 9 (Klf9) gene, which codes for a zinc finger transcription factor involved in development and homeostasis of diverse tissues. We identified regions of the Xenopus and mouse Klf9 genes 5-6 kb upstream of the transcription start sites that supported synergistic transactivation by TH plus GC. Within these regions, we found an orthologous sequence of approximately 180 bp that is highly conserved among tetrapods, but absent in other chordates, and possesses chromatin marks characteristic of an enhancer element. The Xenopus and mouse approximately 180-bp DNA element conferred synergistic transactivation by hormones in transient transfection assays, so we designate this the Klf9 synergy module (KSM). We identified binding sites within the mouse KSM for TH receptor, GC receptor, and nuclear factor κB. TH strongly increased recruitment of liganded GC receptor and serine 5 phosphorylated (initiating) RNA polymerase II to chromatin at the KSM, suggesting a mechanism for transcriptional synergy. The KSM is transcribed to generate long noncoding RNAs, which are also synergistically induced by combined hormone treatment, and the KSM interacts with the Klf9 promoter and a far upstream region through chromosomal looping. Our findings support that the KSM plays a central role in hormone regulation of vertebrate Klf9 genes, it evolved in the tetrapod lineage, and has been maintained by strong stabilizing selection.

Systems consequences of amplicon formation in human breast cancer.

Chromosomal structural variations play an important role in determining the transcriptional landscape of human breast cancers. To assess the nature of these structural variations, we analyzed eight breast tumor samples with a focus on regions of gene amplification using mate-pair sequencing of long-insert genomic DNA with matched transcriptome profiling. We found that tandem duplications appear to be early events in tumor evolution, especially in the genesis of amplicons. In a detailed reconstruction of events on chromosome 17, we found large unpaired inversions and deletions connect a tandemly duplicated ERBB2 with neighboring 17q21.3 amplicons while simultaneously deleting the intervening BRCA1 tumor suppressor locus. This series of events appeared to be unusually common when examined in larger genomic data sets of breast cancers albeit using approaches with lesser resolution. Using siRNAs in breast cancer cell lines, we showed that the 17q21.3 amplicon harbored a significant number of weak oncogenes that appeared consistently coamplified in primary tumors. Down-regulation of BRCA1 expression augmented the cell proliferation in ERBB2-transfected human normal mammary epithelial cells. Coamplification of other functionally tested oncogenic elements in other breast tumors examined, such as RIPK2 and MYC on chromosome 8, also parallel these findings. Our analyses suggest that structural variations efficiently orchestrate the gain and loss of cancer gene cassettes that engage many oncogenic pathways simultaneously and that such oncogenic cassettes are favored during the evolution of a cancer.

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.

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.