Hi-C detects novel structural variants in HL-60 and HL-60/S4 cell lines

Evaluation of Hi-C sequencing for the detection of gene fusions in hematologic and solid pediatric cancer samples

HiC sequencing is a DNA-based next-generation sequencing method that preserves the 3D conformation of the genome and has shown promise in detecting genomic rearrangements in translational research studies. To evaluate HiC as a potential clinical diagnostic platform, analytical concordance with routine laboratory testing was assessed using primary pediatric leukemia and sarcoma specimens previously positive for clinically significant genomic rearrangements. Archived specimen types tested included viable and nonviable frozen leukemic cells, as well as formalin-fixed paraffin-embedded (FFPE) tumor tissues. Initially, pediatric acute myeloid leukemia (AML) and alveolar rhabdomyosarcoma (A-RMS) specimens with known genomic rearrangements were subjected to HiC analysis to assess analytical concordance. Subsequently, a discovery cohort consisting of AML and acute lymphoblastic leukemia (ALL) cases with no known genomic rearrangements based on prior clinical diagnostic testing were evaluated to determine whether HiC could detect rearrangements. Using a standard sequencing depth of 50 million raw read-pairs per sample, or approximately 5X raw genomic coverage, 100% concordance was observed between HiC and previous clinical cytogenetic and molecular testing. In the discovery cohort, a clinically relevant gene fusion was detected in 45% of leukemia cases (5/11). This study demonstrates the value of HiC sequencing to medical diagnostic testing as it identified several clinically significant rearrangements, including those that might have been missed by current clinical testing workflows.

Key points

HiC sequencing is a DNA-based next-generation sequencing method that preserves the 3D conformation of the genome, facilitating detection of genomic rearrangements.HiC was 100% concordant with clinical diagnostic testing workflows for detecting clinically significant genomic rearrangements in pediatric leukemia and rhabdomyosarcoma specimens.HiC detected clinically significant genomic rearrangements not previously detected by prior clinical cytogenetic and molecular testing.HiC performed well with archived non-viable and viable frozen leukemic cell samples, as well as archived formalin-fixed paraffin-embedded tumor tissue specimens.

Hi-C analysis of genomic contacts revealed karyotype abnormalities in chicken HD3 cell line

BACKGROUND

Karyotype abnormalities are frequent in immortalized continuous cell lines either transformed or derived from primary tumors. Chromosomal rearrangements can cause dramatic changes in gene expression and affect cellular phenotype and behavior during in vitro culture. Structural variations of chromosomes in many continuous mammalian cell lines are well documented, but chromosome aberrations in cell lines from other vertebrate models often remain understudied. The chicken LSCC-HD3 cell line (HD3), generated from erythroid precursors, was used as an avian model for erythroid differentiation and lineage-specific gene expression. However, karyotype abnormalities in the HD3 cell line were not assessed. In the present study, we applied high-throughput chromosome conformation capture to analyze 3D genome organization and to detect chromosome rearrangements in the HD3 cell line.

RESULTS

We obtained Hi-C maps of genomic interactions for the HD3 cell line and compared A/B compartments and topologically associating domains between HD3 and several other cell types. By analysis of contact patterns in the Hi-C maps of HD3 cells, we identified more than 25 interchromosomal translocations of regions ≥ 200 kb on both micro- and macrochromosomes. We classified most of the observed translocations as unbalanced, leading to the formation of heteromorphic chromosomes. In many cases of microchromosome rearrangements, an entire microchromosome together with other macro- and microchromosomes participated in the emergence of a derivative chromosome, resembling "chromosomal fusions'' between acrocentric microchromosomes. Intrachromosomal inversions, deletions and duplications were also detected in HD3 cells. Several of the identified simple and complex chromosomal rearrangements, such as between GGA2 and GGA1qter; GGA5, GGA4p and GGA7p; GGA4q, GGA6 and GGA19; and duplication of the sex chromosome GGAW, were confirmed by FISH.

CONCLUSIONS

In the erythroid progenitor HD3 cell line, in contrast to mature and immature erythrocytes, the genome is organized into distinct topologically associating domains. The HD3 cell line has a severely rearranged karyotype with most of the chromosomes engaged in translocations and can be used in studies of genome structure-function relationships. Hi-C proved to be a reliable tool for simultaneous assessment of the spatial genome organization and chromosomal aberrations in karyotypes of birds with a large number of microchromosomes.

A global high-density chromatin interaction network reveals functional long-range and trans-chromosomal relationships

BACKGROUND

Chromatin contacts are essential for gene-expression regulation; however, obtaining a high-resolution genome-wide chromatin contact map is still prohibitively expensive owing to large genome sizes and the quadratic scale of pairwise data. Chromosome conformation capture (3C)-based methods such as Hi-C have been extensively used to obtain chromatin contacts. However, since the sparsity of these maps increases with an increase in genomic distance between contacts, long-range or trans-chromatin contacts are especially challenging to sample.

RESULTS

Here, we create a high-density reference genome-wide chromatin contact map using a meta-analytic approach. We integrate 3600 human, 6700 mouse, and 500 fly Hi-C experiments to create species-specific meta-Hi-C chromatin contact maps with 304 billion, 193 billion, and 19 billion contacts in respective species. We validate that meta-Hi-C contact maps are uniquely powered to capture functional chromatin contacts in both cis and trans. We find that while individual dataset Hi-C networks are largely unable to predict any long-range coexpression (median 0.54 AUC), meta-Hi-C networks perform comparably in both cis and trans (0.65 AUC vs 0.64 AUC). Similarly, for long-range expression quantitative trait loci (eQTL), meta-Hi-C contacts outperform all individual Hi-C experiments, providing an improvement over the conventionally used linear genomic distance-based association. Assessing between species, we find patterns of chromatin contact conservation in both cis and trans and strong associations with coexpression even in species for which Hi-C data is lacking.

CONCLUSIONS

We have generated an integrated chromatin interaction network which complements a large number of methodological and analytic approaches focused on improved specificity or interpretation. This high-depth "super-experiment" is surprisingly powerful in capturing long-range functional relationships of chromatin interactions, which are now able to predict coexpression, eQTLs, and cross-species relationships. The meta-Hi-C networks are available at https://labshare.cshl.edu/shares/gillislab/resource/HiC/ .

Chromosomal Translocations Detection in Cancer Cells Using Chromosomal Conformation Capture Data

Complex chromosomal rearrangements such as translocations play a critical role in oncogenesis. Translocation detection is vital to decipher their biological role in activating cancer-associated mechanisms. High-throughput chromosomal conformations capture (Hi-C) data have shown promising progress in unveiling the genome variations in a disease condition. Until now, multiple structural data (Hi-C)-based methods are available that can detect translocations in cancer genomes. However, the consistency and specificity of Hi-C-based translocation results still need to be validated with conventional methods. This study used Hi-C data of cancerous cell lines, namely lung cancer (A549), Chronic Myelogenous Leukemia (K562), and Acute Monocytic Leukemia (THP-1), to detect the translocations. The results were cross-validated through whole-genome sequencing (WGS) and paired-read analysis. Moreover, PCR amplification validated the presence of translocated reads in different chromosomes. By integrating different data types, we showed that the results of Hi-C data are as reliable as WGS and can be utilized as an assistive method for detecting translocations in the diseased genome. Our findings support the utility of Hi-C technology to detect the translocations and study their effects on the three-dimensional architecture of the genome in cancer condition.

Epigenetic Silencing of PTEN and Epi-Transcriptional Silencing of MDM2 Underlied Progression to Secondary Acute Myeloid Leukemia in Myelodysplastic Syndrome Treated with Hypomethylating Agents

In myelodysplastic syndrome (MDS), resistance to hypomethylating agents (HMA) portends a poor prognosis, underscoring the importance of understanding the molecular mechanisms leading to HMA-resistance. In this study, P39 and Kasumi-1 cells and their azacitidine-resistant and decitabine-resistant sublines were evaluated comparatively with transcriptomic and methylomic analyses. Expression profiling and genome-wide methylation microarray showed downregulation of PTEN associated with DNA hypermethylation in P39 cell lines resistant to azacitidine and decitabine. This pattern of PTEN dysregulation was also confirmed in a cohort of patients failing treatment with HMA. DNA hypomethylation of MDM2 was detected with downregulation of MDM2 in HMA resistant cell lines. Long-read sequencing revealed significant RNA hypomethylation of MDM2 resulting in alternative splicing and production of a truncated MDM2 transcript in azacitidine-resistant P39 cells. The expression of this MDM2 truncated transcript was also significantly increased in HMA-resistant patients compared with HMA-responsive patients. In conclusion, epigenetic and epi-transcriptomic dysregulation of PTEN and MDM2 were associated with resistance to hypomethylating agents.

Neutrophils Promote Larynx Squamous Cell Carcinoma Progression via Activating the IL-17/JAK/STAT3 Pathway

Laryngeal squamous cell carcinoma (LSCC) is the main type of laryngeal cancer with poor prognosis. Incidence of LSCC increases every year, posing a great threat to human health. The underlying mechanism needs further study. Neutrophils are the most prevalent type of immune cells, which play vital roles in crosstalk between the microenvironment and cancer cells. In our study, we aim to figure out the complex regulation between neutrophils and LSCC. Our experiments showed that LSCC cells could promote the activation and mobility of neutrophils. And, in return, neutrophils enhanced the proliferation, migration, and invasion of LSCC. The subsequent results showed that IL-17 was highly expressed in neutrophil conditioned medium. Block of IL-17 could effectively inhibit the progression of LSCC induced by neutrophils. What is more, the results showed that IL-17 activated the JAK/STAT3 pathway in LSCC. Inhibition of the JAK/STAT3 pathway could significantly block neutrophil-induced LSCC progression. Our research reveals the complex interaction between neutrophils and LSCC cells, providing new ideas for the treatment of LSCC.

Unravelling the Shared Genetic Mechanisms Underlying 18 Autoimmune Diseases Using a Systems Approach

Autoimmune diseases (AiDs) are complex heterogeneous diseases characterized by hyperactive immune responses against self. Genome-wide association studies have identified thousands of single nucleotide polymorphisms (SNPs) associated with several AiDs. While these studies have identified a handful of pleiotropic loci that confer risk to multiple AiDs, they lack the power to detect shared genetic factors residing outside of these loci. Here, we integrated chromatin contact, expression quantitative trait loci and protein-protein interaction (PPI) data to identify genes that are regulated by both pleiotropic and non-pleiotropic SNPs. The PPI analysis revealed complex interactions between the shared and disease-specific genes. Furthermore, pathway enrichment analysis demonstrated that the shared genes co-occur with disease-specific genes within the same biological pathways. In conclusion, our results are consistent with the hypothesis that genetic risk loci associated with multiple AiDs converge on a core set of biological processes that potentially contribute to the emergence of polyautoimmunity.

Transcriptional drug repositioning and cheminformatics approach for differentiation therapy of leukaemia cells

Differentiation therapy is attracting increasing interest in cancer as it can be more specific than conventional chemotherapy approaches, and it has offered new treatment options for some cancer types, such as treating acute promyelocytic leukaemia (APL) by retinoic acid. However, there is a pressing need to identify additional molecules which act in this way, both in leukaemia and other cancer types. In this work, we hence developed a novel transcriptional drug repositioning approach, based on both bioinformatics and cheminformatics components, that enables selecting such compounds in a more informed manner. We have validated the approach for leukaemia cells, and retrospectively retinoic acid was successfully identified using our method. Prospectively, the anti-parasitic compound fenbendazole was tested in leukaemia cells, and we were able to show that it can induce the differentiation of leukaemia cells to granulocytes in low concentrations of 0.1 μM and within as short a time period as 3 days. This work hence provides a systematic and validated approach for identifying small molecules for differentiation therapy in cancer.

Structural variation in the sequencing era

Identifying structural variation (SV) is essential for genome interpretation but has been historically difficult due to limitations inherent to available genome technologies. Detection methods that use ensemble algorithms and emerging sequencing technologies have enabled the discovery of thousands of SVs, uncovering information about their ubiquity, relationship to disease and possible effects on biological mechanisms. Given the variability in SV type and size, along with unique detection biases of emerging genomic platforms, multiplatform discovery is necessary to resolve the full spectrum of variation. Here, we review modern approaches for investigating SVs and proffer that, moving forwards, studies integrating biological information with detection will be necessary to comprehensively understand the impact of SV in the human genome.

Neutrophil-like HL-60 cells expressing only GFP-tagged β-actin exhibit nearly normal motility

Observations of actin dynamics in living cells using fluorescence microscopy have been foundational in the exploration of the mechanisms underlying cell migration. We used CRISPR/Cas9 gene editing to generate neutrophil‐like HL‐60 cell lines expressing GFP‐β‐actin from the endogenous locus (ACTB). In light of many previous reports outlining functional deficiencies of labeled actin, we anticipated that HL‐60 cells would only tolerate a monoallelic edit, as biallelic edited cells would produce no normal β‐actin. Surprisingly, we recovered viable monoallelic GFP‐β‐actin cells as well as biallelic edited GFP‐β‐actin cells, in which one copy of the ACTB gene is silenced and the other contains the GFP tag. Furthermore, the edited cells migrate with similar speeds and persistence as unmodified cells in a variety of motility assays, and have nearly normal cell shapes. These results might partially be explained by our observation that GFP‐β‐actin incorporates into the F‐actin network in biallelic edited cells at similar efficiencies as normal β‐actin in unedited cells. Additionally, the edited cells significantly upregulate γ‐actin, perhaps helping to compensate for the loss of normal β‐actin. Interestingly, biallelic edited cells have only modest changes in global gene expression relative to the monoallelic line, as measured by RNA sequencing. While monoallelic edited cells downregulate expression of the tagged allele and are thus only weakly fluorescent, biallelic edited cells are quite bright and well‐suited for live cell microscopy. The nondisruptive phenotype and direct interpretability of this fluorescent tagging approach make it a promising tool for studying actin dynamics in these rapidly migrating and highly phagocytic cells.

Challenges in identifying large germline structural variants for clinical use by long read sequencing

Genomic structural variations, previously considered rare events, are widely recognized as a major source of inter-individual variability and hence, a major hurdle in optimum patient stratification and disease management. Herein, we focus on large complex germline structural variations and present challenges towards target treatment via the synergy of state-of-the-art approaches and information technology tools. A complex structural variation detection remains challenging, as there is no gold standard for identifying such genomic variations with long reads, especially when the chromosomal rearrangement in question is a few Mb in length. A clinical case with a large complex chromosomal rearrangement serves as a paradigm. We feel that functional validation and data interpretation are of outmost importance for information growth to be translated into knowledge growth and hence, new working practices are highlighted.

TNF-α Differentially Regulates Cell Cycle Genes in Promyelocytic and Granulocytic HL-60/S4 Cells

Tumor necrosis factor alpha (TNF-α) is a potent cytokine involved in systemic inflammation and immune modulation. Signaling responses that involve TNF-α are context dependent and capable of stimulating pathways promoting both cell death and survival. TNF-α treatment has been investigated as part of a combined therapy for acute myeloid leukemia due to its modifying effects on all-trans retinoic acid (ATRA) mediated differentiation into granulocytes. To investigate the interaction between cellular differentiation and TNF-α, we performed RNA-sequencing on two forms of the human HL-60/S4 promyelocytic leukemia cell line treated with TNF-α. The ATRA-differentiated granulocytic form of HL-60/S4 cells had an enhanced transcriptional response to TNF-α treatment compared to the undifferentiated promyelocytes. The observed TNF-α responses included differential expression of cell cycle gene sets, which were generally upregulated in TNF-α treated promyelocytes, and downregulated in TNF-α treated granulocytes. This is consistent with TNF-α induced cell cycle repression in granulocytes and cell cycle progression in promyelocytes. Moreover, we found evidence that TNF-α treatment of granulocytes shifts the transcriptome toward that of a macrophage. We conclude that TNF-α treatment promotes a divergent transcriptional program in promyelocytes and granulocytes. TNF-α promotes cell cycle associated gene expression in promyelocytes. In contrast, TNF-α stimulated granulocytes have reduced cell cycle gene expression, and a macrophage-like transcriptional program.

Characterization of Structural Variations in the Context of 3D Chromatin Structure

Chromosomes located in the nucleus form discrete units of genetic material composed of DNA and protein complexes. The genetic information is encoded in linear DNA sequences, but its interpretation requires an understanding of threedimensional (3D) structure of the chromosome, in which distant DNA sequences can be juxtaposed by highly condensed chromatin packing in the space of nucleus to precisely control gene expression. Recent technological innovations in exploring higher-order chromatin structure have uncovered organizational principles of the 3D genome and its various biological implications. Very recently, it has been reported that large-scale genomic variations may disrupt higher-order chromatin organization and as a consequence, greatly contribute to disease-specific gene regulation for a range of human diseases. Here, we review recent developments in studying the effect of structural variation in gene regulation, and the detection and the interpretation of structural variations in the context of 3D chromatin structure.