Comprehensive, integrated, and phased whole-genome analysis of the primary ENCODE cell line K562

Interface-guided phenotyping of coding variants in the transcription factor RUNX1

Single-gene missense mutations remain challenging to interpret. Here, we deploy scalable functional screening by sequencing (SEUSS), a Perturb-seq method, to generate mutations at protein interfaces of RUNX1 and quantify their effect on activities of downstream cellular programs. We evaluate single-cell RNA profiles of 115 mutations in myelogenous leukemia cells and categorize them into three functionally distinct groups, wild-type (WT)-like, loss-of-function (LoF)-like, and hypomorphic, that we validate in orthogonal assays. LoF-like variants dominate the DNA-binding site and are recurrent in cancer; however, recurrence alone does not predict functional impact. Hypomorphic variants share characteristics with LoF-like but favor protein interactions, promoting gene expression indicative of nerve growth factor (NGF) response and cytokine recruitment of neutrophils. Accessible DNA near differentially expressed genes frequently contains RUNX1-binding motifs. Finally, we reclassify 16 variants of uncertain significance and train a classifier to predict 103 more. Our work demonstrates the potential of targeting protein interactions to better define the landscape of phenotypes reachable by missense mutations.

Simultaneous assessment of human genome and methylome data in a single experiment using limited deamination of methylated cytosine

Multiomics require concerted recording of independent information, ideally from a single experiment. In this study, we introduce RIMS-seq2, a high-throughput technique to simultaneously sequence genomes and overlay methylation information while requiring only a small modification of the experimental protocol for high-throughput DNA sequencing to include a controlled deamination step. Importantly, the rate of deamination of 5-methylcytosine is negligible and thus does not interfere with standard DNA sequencing and data processing. Thus, RIMS-seq2 libraries from whole- or targeted-genome sequencing show the same germline variation calling accuracy and sensitivity compared with standard DNA-seq. Additionally, regional methylation levels provide an accurate map of the human methylome.

mLiftOver: harmonizing data across Infinium DNA methylation platforms

MOTIVATION

Infinium DNA methylation BeadChips are widely used for genome-wide DNA methylation profiling at the population scale. Recent updates to probe content and naming conventions in the EPIC version 2 (EPICv2) arrays have complicated integrating new data with previous Infinium array platforms, such as the MethylationEPIC (EPIC) and the HumanMethylation450 (HM450) BeadChip.

RESULTS

We present mLiftOver, a user-friendly tool that harmonizes probe ID, methylation level, and signal intensity data across different Infinium platforms. It manages probe replicates, missing data imputation, and platform-specific bias for accurate data conversion. We validated the tool by applying HM450-based cancer classifiers to EPICv2 cancer data, achieving high accuracy. Additionally, we successfully integrated EPICv2 healthy tissue data with legacy HM450 data for tissue identity analysis and produced consistent copy number profiles in cancer cells.

AVAILABILITY AND IMPLEMENTATION

mLiftOver is implemented R and available in the Bioconductor package SeSAMe (version 1.21.13+): https://bioconductor.org/packages/release/bioc/html/sesame.html. Analysis of EPIC and EPICv2 platform-specific bias and high-confidence mapping is available at https://github.com/zhou-lab/InfiniumAnnotationV1/raw/main/Anno/EPICv2/EPICv2ToEPIC_conversion.tsv.gz. The source code is available at https://github.com/zwdzwd/sesame/blob/devel/R/mLiftOver.R under the MIT license.

Circumvention of Topoisomerase IIα Intron 19 Intronic Polyadenylation in Acquired Etoposide-Resistant Human Leukemia K562 Cells

DNA topoisomerase IIα (TOP2α; 170 kDa, TOP2α/170) is an essential enzyme for proper chromosome dysjunction by producing transient DNA double-stranded breaks and is an important target for DNA damage-stabilizing anticancer agents, such as etoposide. Therapeutic effects of TOP2α poisons can be limited due to acquired drug resistance. We previously demonstrated decreased TOP2α/170 levels in an etoposide-resistant human leukemia K562 subline, designated K/VP.5, accompanied by increased expression of a C-terminal truncated TOP2α isoform (90 kDa; TOP2α/90), which heterodimerized with TOP2α/170 and was a determinant of resistance by exhibiting dominant-negative effects against etoposide activity. Based on 3'-rapid amplification of cDNA ends, we confirmed TOP2α/90 as the translation product of a TOP2α mRNA in which a cryptic polyadenylation site (PAS) harbored in intron 19 (I19) was used. In this report, we investigated whether the resultant intronic polyadenylation (IPA) would be attenuated by blocking or mutating the I19 PAS, thereby circumventing acquired drug resistance. An antisense morpholino oligonucleotide was used to hybridize/block the PAS in TOP2α pre-mRNA in K/VP.5 cells, resulting in decreased TOP2α/90 mRNA/protein levels in K/VP.5 cells and partially circumventing drug resistance. Subsequently, CRISPR/CRISPR-associated protein 9 with homology-directed repair was used to mutate the cryptic I19 PAS (AATAAA→ACCCAA) to prevent IPA. Gene-edited clones exhibited increased TOP2α/170 and decreased TOP2α/90 mRNA/protein and demonstrated restored sensitivity to etoposide and other TOP2α-targeted drugs. Together, results indicated that blocking/mutating a cryptic I19 PAS in K/VP.5 cells reduced IPA and restored sensitivity to TOP2α-targeting drugs. SIGNIFICANCE STATEMENT: The results presented in this study indicate that CRISPR/CRISPR-associated protein 9 gene editing of a cryptic polyadenylation site (PAS) within I19 of the TOP2α gene results in the reversal of acquired resistance to etoposide and other TOP2-targeted drugs. An antisense morpholino oligonucleotide targeting the PAS also partially circumvented resistance.

Iterative deep learning-design of human enhancers exploits condensed sequence grammar to achieve cell type-specificity

An important and largely unsolved problem in synthetic biology is how to target gene expression to specific cell types. Here, we apply iterative deep learning to design synthetic enhancers with strong differential activity between two human cell lines. We initially train models on published datasets of enhancer activity and chromatin accessibility and use them to guide the design of synthetic enhancers that maximize predicted specificity. We experimentally validate these sequences, use the measurements to re-optimize the predictor, and design a second generation of enhancers with improved specificity. Our design methods embed relevant transcription factor binding site (TFBS) motifs with higher frequencies than comparable endogenous enhancers while using a more selective motif vocabulary, and we show that enhancer activity is correlated with transcription factor expression at the single cell level. Finally, we characterize causal features of top enhancers via perturbation experiments and show enhancers as short as 50bp can maintain specificity.

Linking regulatory variants to target genes by integrating single-cell multiome methods and genomic distance

Methods that analyze single-cell paired RNA-seq and ATAC-seq multiome data have shown great promise in linking regulatory elements to genes. However, existing methods differ in their modeling assumptions and approaches to account for biological and technical noise-leading to low concordance in their linking scores-and do not capture the effects of genomic distance. We propose pgBoost, an integrative modeling framework that trains a non-linear combination of existing linking strategies (including genomic distance) on fine-mapped eQTL data to assign a probabilistic score to each candidate SNP-gene link. We applied pgBoost to single-cell multiome data from 85k cells representing 6 major immune/blood cell types. pgBoost attained higher enrichment for fine-mapped eSNP-eGene pairs (e.g. 21x at distance >10kb) than existing methods (1.2-10x; p-value for difference = 5e-13 vs. distance-based method and < 4e-35 for each other method), with larger improvements at larger distances (e.g. 35x vs. 0.89-6.6x at distance >100kb; p-value for difference < 0.002 vs. each other method). pgBoost also outperformed existing methods in enrichment for CRISPR-validated links (e.g. 4.8x vs. 1.6-4.1x at distance >10kb; p-value for difference = 0.25 vs. distance-based method and < 2e-5 for each other method), with larger improvements at larger distances (e.g. 15x vs. 1.6-2.5x at distance >100kb; p-value for difference < 0.009 for each other method). Similar improvements in enrichment were observed for links derived from Activity-By-Contact (ABC) scores and GWAS data. We further determined that restricting pgBoost to features from a focal cell type improved the identification of SNP-gene links relevant to that cell type. We highlight several examples where pgBoost linked fine-mapped GWAS variants to experimentally validated or biologically plausible target genes that were not implicated by other methods. In conclusion, a non-linear combination of linking strategies, including genomic distance, improves power to identify target genes underlying GWAS associations.

Engineered CRISPR-Cas12a for higher-order combinatorial chromatin perturbations

Multiplexed genetic perturbations are critical for testing functional interactions among coding or non-coding genetic elements. Compared to double-stranded DNA cutting, repressive chromatin formation using CRISPR interference (CRISPRi) avoids genotoxicity and is more effective for perturbing non-coding regulatory elements in pooled assays. However, current CRISPRi pooled screening approaches are limited to targeting one to three genomic sites per cell. We engineer an Acidaminococcus Cas12a (AsCas12a) variant, multiplexed transcriptional interference AsCas12a (multiAsCas12a), that incorporates R1226A, a mutation that stabilizes the ribonucleoprotein-DNA complex via DNA nicking. The multiAsCas12a-KRAB fusion improves CRISPRi activity over DNase-dead AsCas12a-KRAB fusions, often rescuing the activities of lentivirally delivered CRISPR RNAs (crRNA) that are inactive when used with the latter. multiAsCas12a-KRAB supports CRISPRi using 6-plex crRNA arrays in high-throughput pooled screens. Using multiAsCas12a-KRAB, we discover enhancer elements and dissect the combinatorial function of cis-regulatory elements in human cells. These results instantiate a group testing framework for efficiently surveying numerous combinations of chromatin perturbations for biological discovery and engineering.

Multiplex, single-cell CRISPRa screening for cell type specific regulatory elements

CRISPR-based gene activation (CRISPRa) is a promising therapeutic approach for gene therapy, upregulating gene expression by targeting promoters or enhancers in a tissue/cell-type specific manner. Here, we describe an experimental framework that combines highly multiplexed perturbations with single-cell RNA sequencing (sc-RNA-seq) to identify cell-type-specific, CRISPRa-responsive cis- regulatory elements and the gene(s) they regulate. Random combinations of many gRNAs are introduced to each of many cells, which are then profiled and partitioned into test and control groups to test for effect(s) of CRISPRa perturbations of both enhancers and promoters on the expression of neighboring genes. Applying this method to a library of 493 gRNAs targeting candidate cis- regulatory elements in both K562 cells and iPSC-derived excitatory neurons, we identify gRNAs capable of specifically upregulating intended target genes and no other neighboring genes within 1 Mb, including gRNAs yielding upregulation of six autism spectrum disorder (ASD) and neurodevelopmental disorder (NDD) risk genes in neurons. A consistent pattern is that the responsiveness of individual enhancers to CRISPRa is restricted by cell type, implying a dependency on either chromatin landscape and/or additional trans- acting factors for successful gene activation. The approach outlined here may facilitate large-scale screens for gRNAs that activate therapeutically relevant genes in a cell type-specific manner.

Variable characteristics overlooked in human K-562 leukemia cell lines with a common signature

K-562 is a well-known in vitro cellular model that represents human leukemia cell lines. Although the K-562 cells have been extensively characterized, there are inconsistencies in the data across publications, showing the presence of multiple K-562 cell lines. This suggests that analyzing a single K-562 cell line is insufficient to provide reliable reference data. In this study, we compared three K-562 cell lines with different IDs (RCB0027, RCB1635, and RCB1897) to investigate the fundamental characteristics of K-562 cells. Amplifications of the BCR-ABL1 fusion gene and at 13q31 were detected in all three cell lines, whereas each genome exhibited distinctive features of sequence variants and loss of heterozygosity. This implies that each K-562 cell line can be characterized by common and unique features through a comparison of multiple K-562 cell lines. Variations in transcriptome profiles and hemoglobin synthesis were also observed among the three cell lines, indicating that they should be considered sublines that have diverged from the common ancestral K-562 despite no changes from the original cell name. This leads to unintentional differences in genotypes and/or phenotypes among cell lines that share the same name. These data show that characterizing a single K-562 cell line does not necessarily provide data that are applicable to other K-562 cells. In this context, it is essential to modify cell names in accordance with changes in characteristics during cell culture. Furthermore, our data could serve as a reference for evaluating other K-562 sublines, facilitating the discovery of new K-562 sublines with distinct characteristics. This approach results in the accumulation of K-562 sublines with diverged characteristics and expands the options available, which may help in selecting the most suitable K-562 subline for each experiment.

Use of CRISPR/Cas9 with Homology-Directed Repair to Gene-Edit Topoisomerase IIβ in Human Leukemia K562 Cells: Generation of a Resistance Phenotype

DNA topoisomerase IIβ (TOP2β/180; 180 kDa) is a nuclear enzyme that regulates DNA topology by generation of short-lived DNA double-strand breaks, primarily during transcription. TOP2β/180 can be a target for DNA damage-stabilizing anticancer drugs, whose efficacy is often limited by chemoresistance. Our laboratory previously demonstrated reduced levels of TOP2β/180 (and the paralog TOP2α/170) in an acquired etoposide-resistant human leukemia (K562) clonal cell line, K/VP.5, in part due to overexpression of microRNA-9-3p/5p impacting post-transcriptional events. To evaluate the effect on drug sensitivity upon reduction/elimination of TOP2β/180, a premature stop codon was generated at the TOP2β/180 gene exon 19/intron 19 boundary (AGAA//GTAA→ATAG//GTAA) in parental K562 cells (which contain four TOP2β/180 alleles) by CRISPR/Cas9 editing with homology-directed repair to disrupt production of full-length TOP2β/180. Gene-edited clones were identified and verified by quantitative polymerase chain reaction and Sanger sequencing, respectively. Characterization of TOP2β/180 gene-edited clones, with one or all four TOP2β/180 alleles mutated, revealed partial or complete loss of TOP2β mRNA/protein, respectively. The loss of TOP2β/180 protein correlated with decreased (2-{4-[(7-chloro-2-quinoxalinyl)oxy]phenoxy}propionic acid)-induced DNA damage and partial resistance in growth inhibition assays. Partial resistance to mitoxantrone was also noted in the gene-edited clone with all four TOP2β/180 alleles modified. No cross-resistance to etoposide or mAMSA was noted in the gene-edited clones. Results demonstrated the role of TOP2β/180 in drug sensitivity/resistance in K562 cells and revealed differential paralog activity of TOP2-targeted agents. SIGNIFICANCE STATEMENT: Data indicated that CRISPR/Cas9 editing of the exon 19/intron 19 boundary in the TOP2β/180 gene to introduce a premature stop codon resulted in partial to complete disruption of TOP2β/180 expression in human leukemia (K562) cells depending on the number of edited alleles. Edited clones were partially resistant to mitoxantrone and XK469, while lacking resistance to etoposide and mAMSA. Results demonstrated the import of TOP2β/180 in drug sensitivity/resistance in K562 cells and revealed differential paralog activity of TOP2-targeted agents.

Higher-order combinatorial chromatin perturbations by engineered CRISPR-Cas12a for functional genomics

Multiplexed genetic perturbations are critical for testing functional interactions among coding or non-coding genetic elements. Compared to double-stranded DNA cutting, repressive chromatin formation using CRISPR interference (CRISPRi) avoids genotoxicity and is more effective for perturbing non-coding regulatory elements in pooled assays. However, current CRISPRi pooled screening approaches are limited to targeting 1-3 genomic sites per cell. To develop a tool for higher-order ( > 3) combinatorial targeting of genomic sites with CRISPRi in functional genomics screens, we engineered an Acidaminococcus Cas12a variant -- referred to as mul tiplexed transcriptional interference AsCas12a (multiAsCas12a). multiAsCas12a incorporates a key mutation, R1226A, motivated by the hypothesis of nicking-induced stabilization of the ribonucleoprotein:DNA complex for improving CRISPRi activity. multiAsCas12a significantly outperforms prior state-of-the-art Cas12a variants in combinatorial CRISPRi targeting using high-order multiplexed arrays of lentivirally transduced CRISPR RNAs (crRNA), including in high-throughput pooled screens using 6-plex crRNA array libraries. Using multiAsCas12a CRISPRi, we discover new enhancer elements and dissect the combinatorial function of cis-regulatory elements. These results instantiate a group testing framework for efficiently surveying potentially numerous combinations of chromatin perturbations for biological discovery and engineering.

Small-molecule correctors divert CFTR-F508del from ERAD by stabilizing sequential folding states

Over 80% of people with cystic fibrosis (CF) carry the F508del mutation in the cystic fibrosis transmembrane conductance regulator (CFTR), a chloride ion channel at the apical plasma membrane (PM) of epithelial cells. F508del impairs CFTR folding causing it to be destroyed by endoplasmic reticulum associated degradation (ERAD). Small-molecule correctors, which act as pharmacological chaperones to divert CFTR-F508del from ERAD, are the primary strategy for treating CF, yet corrector development continues with only a rudimentary understanding of how ERAD targets CFTR-F508del. We conducted genome-wide CRISPR/Cas9 knockout screens to systematically identify the molecular machinery that underlies CFTR-F508del ERAD. Although the ER-resident ubiquitin ligase, RNF5 was the top E3 hit, knocking out RNF5 only modestly reduced CFTR-F508del degradation. Sublibrary screens in an RNF5 knockout background identified RNF185 as a redundant ligase and demonstrated that CFTR-F508del ERAD is robust. Gene-drug interaction experiments illustrated that correctors tezacaftor (VX-661) and elexacaftor (VX-445) stabilize sequential, RNF5-resistant folding states. We propose that binding of correctors to nascent CFTR-F508del alters its folding landscape by stabilizing folding states that are not substrates for RNF5-mediated ubiquitylation.

Evidence for widespread existence of functional novel and non-canonical human transcripts

BACKGROUND

Fraction of functional sequence in the human genome remains a key unresolved question in Biology and the subject of vigorous debate. While a plethora of studies have connected a significant fraction of human DNA to various biochemical processes, the classical definition of function requires evidence of effects on cellular or organismal fitness that such studies do not provide. Although multiple high-throughput reverse genetics screens have been developed to address this issue, they are limited to annotated genomic elements and suffer from non-specific effects, arguing for a strong need to develop additional functional genomics approaches.

RESULTS

In this work, we established a high-throughput lentivirus-based insertional mutagenesis strategy as a forward genetics screen tool in aneuploid cells. Application of this approach to human cell lines in multiple phenotypic screens suggested the presence of many yet uncharacterized functional elements in the human genome, represented at least in part by novel exons of known and novel genes. The novel transcripts containing these exons can be massively, up to thousands-fold, induced by specific stresses, and at least some can represent bi-cistronic protein-coding mRNAs.

CONCLUSIONS

Altogether, these results argue that many unannotated and non-canonical human transcripts, including those that appear as aberrant splice products, have biological relevance under specific biological conditions.

Improving Hi-C contact matrices using genome graphs

Three-dimensional chromosome structure plays an important role in fundamental genomic functions. Hi-C, a high-throughput, sequencing-based technique, has drastically expanded our comprehension of 3D chromosome structures. The first step of Hi-C analysis pipeline involves mapping sequencing reads from Hi-C to linear reference genomes. However, the linear reference genome does not incorporate genetic variation information, which can lead to incorrect read alignments, especially when analyzing samples with substantial genomic differences from the reference such as cancer samples. Using genome graphs as the reference facilitates more accurate mapping of reads, however, new algorithms are required for inferring linear genomes from Hi-C reads mapped on genome graphs and constructing corresponding Hi-C contact matrices, which is a prerequisite for the subsequent steps of the Hi-C analysis such as identifying topologically associated domains and calling chromatin loops. We introduce the problem of genome sequence inference from Hi-C data mediated by genome graphs. We formalize this problem, show the hardness of solving this problem, and introduce a novel heuristic algorithm specifically tailored to this problem. We provide a theoretical analysis to evaluate the efficacy of our algorithm. Finally, our empirical experiments indicate that the linear genomes inferred from our method lead to the creation of improved Hi-C contact matrices. These enhanced matrices show a reduction in erroneous patterns caused by structural variations and are more effective in accurately capturing the structures of topologically associated domains.

Tracing cancer evolution and heterogeneity using Hi-C

Chromosomal rearrangements can initiate and drive cancer progression, yet it has been challenging to evaluate their impact, especially in genetically heterogeneous solid cancers. To address this problem we developed HiDENSEC, a new computational framework for analyzing chromatin conformation capture in heterogeneous samples that can infer somatic copy number alterations, characterize large-scale chromosomal rearrangements, and estimate cancer cell fractions. After validating HiDENSEC with in silico and in vitro controls, we used it to characterize chromosome-scale evolution during melanoma progression in formalin-fixed tumor samples from three patients. The resulting comprehensive annotation of the genomic events includes copy number neutral translocations that disrupt tumor suppressor genes such as NF1, whole chromosome arm exchanges that result in loss of CDKN2A, and whole-arm copy-number neutral loss of homozygosity involving PTEN. These findings show that large-scale chromosomal rearrangements occur throughout cancer evolution and that characterizing these events yields insights into drivers of melanoma progression.

The chromatin network helps prevent cancer-associated mutagenesis at transcription-replication conflicts

Genome instability is a feature of cancer cells, transcription being an important source of DNA damage. This is in large part associated with R-loops, which hamper replication, especially at head-on transcription-replication conflicts (TRCs). Here we show that TRCs trigger a DNA Damage Response (DDR) involving the chromatin network to prevent genome instability. Depletion of the key chromatin factors INO80, SMARCA5 and MTA2 results in TRCs, fork stalling and R-loop-mediated DNA damage which mostly accumulates at S/G2, while histone H3 Ser10 phosphorylation, a mark of chromatin compaction, is enriched at TRCs. Strikingly, TRC regions show increased mutagenesis in cancer cells with signatures of homologous recombination deficiency, transcription-coupled nucleotide excision repair (TC-NER) and of the AID/APOBEC cytidine deaminases, being predominant at head-on collisions. Thus, our results support that the chromatin network prevents R-loops and TRCs from genomic instability and mutagenic signatures frequently associated with cancer.

scNanoHi-C: a single-cell long-read concatemer sequencing method to reveal high-order chromatin structures within individual cells

The high-order three-dimensional (3D) organization of regulatory genomic elements provides a topological basis for gene regulation, but it remains unclear how multiple regulatory elements across the mammalian genome interact within an individual cell. To address this, herein, we developed scNanoHi-C, which applies Nanopore long-read sequencing to explore genome-wide proximal high-order chromatin contacts within individual cells. We show that scNanoHi-C can reliably and effectively profile 3D chromatin structures and distinguish structure subtypes among individual cells. This method could also be used to detect genomic variations, including copy-number variations and structural variations, as well as to scaffold the de novo assembly of single-cell genomes. Notably, our results suggest that extensive high-order chromatin structures exist in active chromatin regions across the genome, and multiway interactions between enhancers and their target promoters were systematically identified within individual cells. Altogether, scNanoHi-C offers new opportunities to investigate high-order 3D genome structures at the single-cell level.

Small molecule correctors divert CFTR-F508del from ERAD by stabilizing sequential folding states

Over 80% of people with cystic fibrosis (CF) carry the F508del mutation in the cystic fibrosis transmembrane conductance regulator (CFTR), a chloride ion channel at the apical plasma membrane (PM) of epithelial cells. F508del impairs CFTR folding causing it to be destroyed by endoplasmic reticulum associated degradation (ERAD). Small molecule correctors, which act as pharmacological chaperones to divert CFTR-F508del from ERAD, are the primary strategy for treating CF, yet corrector development continues with only a rudimentary understanding of how ERAD targets CFTR-F508del. We conducted genome-wide CRISPR/Cas9 knockout screens to systematically identify the molecular machinery that underlies CFTR-F508del ERAD. Although the ER-resident ubiquitin ligase, RNF5 was the top E3 hit, knocking out RNF5 only modestly reduced CFTR-F508del degradation. Sublibrary screens in an RNF5 knockout background identified RNF185 as a redundant ligase, demonstrating that CFTR-F508del ERAD is highly buffered. Gene-drug interaction experiments demonstrated that correctors tezacaftor (VX-661) and elexacaftor (VX-445) stabilize sequential, RNF5-resistant folding states. We propose that binding of correctors to nascent CFTR-F508del alters its folding landscape by stabilizing folding states that are not substrates for RNF5-mediated ubiquitylation.

Interface-guided phenotyping of coding variants in the transcription factor RUNX1 with SEUSS

Understanding the consequences of single amino acid substitutions in cancer driver genes remains an unmet need. Perturb-seq provides a tool to investigate the effects of individual mutations on cellular programs. Here we deploy SEUSS, a Perturb-seq like approach, to generate and assay mutations at physical interfaces of the RUNX1 Runt domain. We measured the impact of 115 mutations on RNA profiles in single myelogenous leukemia cells and used the profiles to categorize mutations into three functionally distinct groups: wild-type (WT)-like, loss-of-function (LOF)-like and hypomorphic. Notably, the largest concentration of functional mutations (non-WT-like) clustered at the DNA binding site and contained many of the more frequently observed mutations in human cancers. Hypomorphic variants shared characteristics with loss of function variants but had gene expression profiles indicative of response to neural growth factor and cytokine recruitment of neutrophils. Additionally, DNA accessibility changes upon perturbations were enriched for RUNX1 binding motifs, particularly near differentially expressed genes. Overall, our work demonstrates the potential of targeting protein interaction interfaces to better define the landscape of prospective phenotypes reachable by amino acid substitutions.

Simulated Microgravity-Induced Changes to Drug Response in Cancer Cells Quantified Using Fluorescence Morphometry

Unlike plants that have special gravity-sensing cells, such special cells in animals are yet to be discovered. However, microgravity, the condition of apparent weightlessness, causes bone, muscular and immune system dysfunctions in astronauts following spaceflights. Decades of investigations show correlations between these organ and system-level dysfunctions with changes induced at the cellular level both by simulated microgravity as well as microgravity conditions in outer space. Changes in single bone, muscle and immune cells include morphological abnormalities, altered gene expression, protein expression, metabolic pathways and signaling pathways. These suggest that human cells mount some response to microgravity. However, the implications of such adjustments on many cellular functions and responses are not clear. Here, we addressed the question whether microgravity induces alterations to drug response in cancer cells. We used both adherent cancer cells (T98G) and cancer cells in suspension (K562) to confirm the known effects of simulated microgravity and then treated the K562 cells with common cancer drugs (hydroxyurea and paclitaxel) following 48 h of exposure to simulated microgravity via a NASA-developed rotary cell culture system. Through fluorescence-guided morphometry, we found that microgravity abolished a significant reduction (p < 0.01) in the nuclear-to-cytoplasm ratio of cancer cells treated with hydroxyurea. Our results call for more studies on the impact of microgravity on cellular drug response, in light of the growing need for space medicine, as space exploration grows.

RNA polymerase II dynamics shape enhancer-promoter interactions

How enhancers control target gene expression over long genomic distances remains an important unsolved problem. Here we investigated enhancer-promoter communication by integrating data from nucleosome-resolution genomic contact maps, nascent transcription and perturbations affecting either RNA polymerase II (Pol II) dynamics or the activity of thousands of candidate enhancers. Integration of new Micro-C experiments with published CRISPRi data demonstrated that enhancers spend more time in close proximity to their target promoters in functional enhancer-promoter pairs compared to nonfunctional pairs, which can be attributed in part to factors unrelated to genomic position. Manipulation of the transcription cycle demonstrated a key role for Pol II in enhancer-promoter interactions. Notably, promoter-proximal paused Pol II itself partially stabilized interactions. We propose an updated model in which elements of transcriptional dynamics shape the duration or frequency of interactions to facilitate enhancer-promoter communication.

Linked-read based analysis of the medulloblastoma genome

Introduction

Medulloblastoma is the most common type of malignant pediatric brain tumor with group 4 medulloblastomas (G4 MBs) accounting for 40% of cases. However, the molecular mechanisms that underlie this subgroup are still poorly understood. Point mutations are detected in a large number of genes at low incidence per gene while the detection of complex structural variants in recurrently affected genes typically requires the application of long-read technologies.

Methods

Here, we applied linked-read sequencing, which combines the long-range genome information of long-read sequencing with the high base pair accuracy of short read sequencing and very low sample input requirements.

Results

We demonstrate the detection of complex structural variants and point mutations in these tumors, and, for the first time, the detection of extrachromosomal DNA (ecDNA) with linked-reads. We provide further evidence for the high heterogeneity of somatic mutations in G4 MBs and add new complex events associated with it.

Discussion

We detected several enhancer-hijacking events, an ecDNA containing the MYCN gene, and rare structural rearrangements, such a chromothripsis in a G4 medulloblastoma, chromoplexy involving 8 different chromosomes, a TERT gene rearrangement, and a PRDM6 duplication.

Comparison and benchmark of structural variants detected from long read and long-read assembly

Structural variant (SV) detection is essential for genomic studies, and long-read sequencing technologies have advanced our capacity to detect SVs directly from read or de novo assembly, also known as read-based and assembly-based strategy. However, to date, no independent studies have compared and benchmarked the two strategies. Here, on the basis of SVs detected by 20 read-based and eight assembly-based detection pipelines from six datasets of HG002 genome, we investigated the factors that influence the two strategies and assessed their performance with well-curated SVs. We found that up to 80% of the SVs could be detected by both strategies among different long-read datasets, whereas variant type, size, and breakpoint detected by read-based strategy were greatly affected by aligners. For the high-confident insertions and deletions at non-tandem repeat regions, a remarkable subset of them (82% in assembly-based calls and 93% in read-based calls), accounting for around 4000 SVs, could be captured by both reads and assemblies. However, discordance between two strategies was largely caused by complex SVs and inversions, which resulted from inconsistent alignment of reads and assemblies at these loci. Finally, benchmarking with SVs at medically relevant genes, the recall of read-based strategy reached 77% on 5X coverage data, whereas assembly-based strategy required 20X coverage data to achieve similar performance. Therefore, integrating SVs from read and assembly is suggested for general-purpose detection because of inconsistently detected complex SVs and inversions, whereas assembly-based strategy is optional for applications with limited resources.

HiSV: A control-free method for structural variation detection from Hi-C data

Structural variations (SVs) play an essential role in the evolution of human genomes and are associated with cancer genetics and rare disease. High-throughput chromosome capture (Hi-C) technology probed all genome-wide crosslinked chromatin to study the spatial architecture of chromosomes. Hi-C read pairs can span megabases, making the technology useful for detecting large-scale SVs. So far, the identification of SVs from Hi-C data is still in the early stages with only a few methods available. Especially, no algorithm has been developed that can detect SVs without control samples. Therefore, we developed HiSV (Hi-C for Structural Variation), a control-free method for identifying large-scale SVs from a Hi-C sample. Inspired by the single image saliency detection model, HiSV constructed a saliency map of interaction frequencies and extracted saliency segments as large-scale SVs. By evaluating both simulated and real data, HiSV not only detected all variant types, but also achieved a higher level of accuracy and sensitivity than existing methods. Moreover, our results on cancer cell lines showed that HiSV effectively detected eight complex SV events and identified two novel SVs of key factors associated with cancer development. Finally, we found that integrating the result of HiSV helped the WGS method to identify a total number of 94 novel SVs in two cancer cell lines.

Local assembly of long reads enables phylogenomics of transposable elements in a polyploid cell line

Animal cell lines often undergo extreme genome restructuring events, including polyploidy and segmental aneuploidy that can impede de novo whole-genome assembly (WGA). In some species like Drosophila, cell lines also exhibit massive proliferation of transposable elements (TEs). To better understand the role of transposition during animal cell culture, we sequenced the genome of the tetraploid Drosophila S2R+ cell line using long-read and linked-read technologies. WGAs for S2R+ were highly fragmented and generated variable estimates of TE content across sequencing and assembly technologies. We therefore developed a novel WGA-independent bioinformatics method called TELR that identifies, locally assembles, and estimates allele frequency of TEs from long-read sequence data (https://github.com/bergmanlab/telr). Application of TELR to a ∼130x PacBio dataset for S2R+ revealed many haplotype-specific TE insertions that arose by transposition after initial cell line establishment and subsequent tetraploidization. Local assemblies from TELR also allowed phylogenetic analysis of paralogous TEs, which revealed that proliferation of TE families in vitro can be driven by single or multiple source lineages. Our work provides a model for the analysis of TEs in complex heterozygous or polyploid genomes that are recalcitrant to WGA and yields new insights into the mechanisms of genome evolution in animal cell culture.

Integration of multi-omics data reveals cis-regulatory variants that are associated with phenotypic differentiation of eastern from western pigs

Background

The genetic mechanisms that underlie phenotypic differentiation in breeding animals have important implications in evolutionary biology and agriculture. However, the contribution of cis-regulatory variants to pig phenotypes is poorly understood. Therefore, our aim was to elucidate the molecular mechanisms by which non-coding variants cause phenotypic differences in pigs by combining evolutionary biology analyses and functional genomics.

Results

We obtained a high-resolution phased chromosome-scale reference genome with a contig N50 of 18.03 Mb for the Luchuan pig breed (a representative eastern breed) and profiled potential selective sweeps in eastern and western pigs by resequencing the genomes of 234 pigs. Multi-tissue transcriptome and chromatin accessibility analyses of these regions suggest that tissue-specific selection pressure is mediated by promoters and distal cis-regulatory elements. Promoter variants that are associated with increased expression of the lysozyme (LYZ) gene in the small intestine might enhance the immunity of the gastrointestinal tract and roughage tolerance in pigs. In skeletal muscle, an enhancer-modulating single-nucleotide polymorphism that is associated with up-regulation of the expression of the troponin C1, slow skeletal and cardiac type (TNNC1) gene might increase the proportion of slow muscle fibers and affect meat quality.

Conclusions

Our work sheds light on the molecular mechanisms by which non-coding variants shape phenotypic differences in pigs and provides valuable resources and novel perspectives to dissect the role of gene regulatory evolution in animal domestication and breeding.

The genomic landscape of canine osteosarcoma cell lines reveals conserved structural complexity and pathway alterations

The characterization of immortalized canine osteosarcoma (OS) cell lines used for research has historically been based on phenotypic features such as cellular morphology and expression of bone specific markers. With the increasing use of these cell lines to investigate novel therapeutic approaches prior to in vivo translation, a much more detailed understanding regarding the genomic landscape of these lines is required to ensure accurate interpretation of findings. Here we report the first whole genome characterization of eight canine OS cell lines, including single nucleotide variants, copy number variants and other structural variants. Many alterations previously characterized in primary canine OS tissue were observed in these cell lines, including TP53 mutations, MYC copy number gains, loss of CDKN2A, PTEN, DLG2, MAGI2, and RB1 and structural variants involving SETD2, DLG2 and DMD. These data provide a new framework for understanding how best to incorporate in vitro findings generated using these cell lines into the design of future clinical studies involving dogs with spontaneous OS.

Intronic Polyadenylation in Acquired Cancer Drug Resistance Circumvented by Utilizing CRISPR/Cas9 with Homology-Directed Repair: The Tale of Human DNA Topoisomerase IIα

Simple Summary

DNA topoisomerase IIα (170 kDa, TOP2α/170) resolves nucleic acid topological entanglements by generating transient double-strand DNA breaks. TOP2α inhibitors/poisons stabilize TOP2α-DNA covalent complexes resulting in persistent DNA damage and are frequently utilized to treat a variety of cancers. Acquired resistance to these chemotherapeutic agents is often associated with decreased TOP2α/170 expression levels. Studies have demonstrated that a reduction in TOP2α/170 results from a type of alternative polyadenylation designated intronic polyadenylation (IPA). As a consequence of IPA, variant TOP2α mRNA transcripts have been characterized that have resulted in the translation of C-terminal truncated TOP2α isoforms with altered biological activities. In this paper, an example is discussed where circumvention of acquired TOP2α-mediated drug resistance was achieved by utilizing CRISPR/Cas9 specific gene editing of an exon/intron boundary through homology directed repair (HDR) to reduce TOP2α IPA. These results illustrate the therapeutic potential of CRISPR/Cas9/HDR to impact drug resistance associated with aberrant IPA.

Abstract

Intronic polyadenylation (IPA) plays a critical role in malignant transformation, development, progression, and cancer chemoresistance by contributing to transcriptome/proteome alterations. DNA topoisomerase IIα (170 kDa, TOP2α/170) is an established clinical target for anticancer agents whose efficacy is compromised by drug resistance often associated with a reduction of nuclear TOP2α/170 levels. In leukemia cell lines with acquired resistance to TOP2α-targeted drugs and reduced TOP2α/170 expression, variant TOP2α mRNA transcripts have been reported due to IPA that resulted in the translation of C-terminal truncated isoforms with altered nuclear-cytoplasmic distribution or heterodimerization with wild-type TOP2α/170. This review provides an overview of the various mechanisms regulating pre-mRNA processing and alternative polyadenylation, as well as the utilization of CRISPR/Cas9 specific gene editing through homology directed repair (HDR) to decrease IPA when splice sites are intrinsically weak or potentially mutated. The specific case of TOP2α exon 19/intron 19 splice site editing is discussed in etoposide-resistant human leukemia K562 cells as a tractable strategy to circumvent acquired TOP2α-mediated drug resistance. This example supports the importance of aberrant IPA in acquired drug resistance to TOP2α-targeted drugs. In addition, these results demonstrate the therapeutic potential of CRISPR/Cas9/HDR to impact drug resistance associated with aberrant splicing/polyadenylation.

Building regulatory landscapes reveals that an enhancer can recruit cohesin to create contact domains, engage CTCF sites and activate distant genes

Developmental gene expression is often controlled by distal regulatory DNA elements called enhancers. Distant enhancer action is restricted to structural chromosomal domains that are flanked by CTCF-associated boundaries and formed through cohesin chromatin loop extrusion. To better understand how enhancers, genes and CTCF boundaries together form structural domains and control expression, we used a bottom-up approach, building series of active regulatory landscapes in inactive chromatin. We demonstrate here that gene transcription levels and activity over time reduce with increased enhancer distance. The enhancer recruits cohesin to stimulate domain formation and engage flanking CTCF sites in loop formation. It requires cohesin exclusively for the activation of distant genes, not of proximal genes, with nearby CTCF boundaries supporting efficient long-range enhancer action. Our work supports a dual activity model for enhancers: its classic role of stimulating transcription initiation and elongation from target gene promoters and a role of recruiting cohesin for the creation of chromosomal domains, the engagement of CTCF sites in chromatin looping and the activation of distal target genes.

Use of CRISPR/Cas9 with homology-directed repair to silence the human topoisomerase IIα intron-19 5’ splice site: Generation of etoposide resistance in human leukemia K562 cells

DNA Topoisomerase IIα (TOP2α/170) is an enzyme essential for proliferating cells. For rapidly multiplying malignancies, this has made TOP2α/170 an important target for etoposide and other clinically active anticancer drugs. Efficacy of these agents is often limited by chemoresistance related to alterations in TOP2α/170 expression levels. Our laboratory recently demonstrated reduced levels of TOP2α/170 and overexpression of a C-terminal truncated 90-kDa isoform, TOP2α/90, due to intronic polyadenylation (IPA; within intron 19) in an acquired etoposide-resistant K562 clonal cell line, K/VP.5. We previously reported that this isoform heterodimerized with TOP2α/170 and was a determinant of acquired resistance to etoposide. Optimization of the weak TOP2α exon 19/intron 19 5′ splice site in drug-resistant K/VP.5 cells by gene-editing restored TOP2α/170 levels, diminished TOP2α/90 expression, and circumvented drug resistance. Conversely, in the present study, silencing of the exon 19/intron 19 5′ splice site in parental K562 cells by CRISPR/Cas9 with homology-directed repair (HDR), and thereby forcing intron 19 retention, was used to induce resistance by disrupting normal RNA processing (i.e., gene knockout), and to further evaluate the role of TOP2α/170 and TOP2α/90 isoforms as resistance determinants. Gene-edited clones were identified by quantitative polymerase chain reaction (qPCR) and verified by Sanger sequencing. TOP2α/170 mRNA/protein expression levels were attenuated in the TOP2α gene-edited clones which resulted in resistance to etoposide as assessed by reduced etoposide-induced DNA damage (γH2AX, Comet assays) and growth inhibition. RNA-seq and qPCR studies suggested that intron 19 retention leads to decreased TOP2α/170 expression by degradation of the TOP2α edited mRNA transcripts. Forced expression of TOP2α/90 in the gene-edited K562 cells further decreased etoposide-induced DNA damage in support of a dominant negative role for this truncated isoform. Together results support the important role of both TOP2α/170 and TOP2α/90 as determinants of sensitivity/resistance to TOP2α-targeting agents.

Predicting genotype-specific gene regulatory networks

Understanding how each person's unique genotype influences their individual patterns of gene regulation has the potential to improve our understanding of human health and development, and to refine genotype-specific disease risk assessments and treatments. However, the effects of genetic variants are not typically considered when constructing gene regulatory networks, despite the fact that many disease-associated genetic variants are thought to have regulatory effects, including the disruption of transcription factor (TF) binding. We developed EGRET (Estimating the Genetic Regulatory Effect on TFs), which infers a genotype-specific gene regulatory network for each individual in a study population. EGRET begins by constructing a genotype-informed TF-gene prior network derived using TF motif predictions, expression quantitative trait locus (eQTL) data, individual genotypes, and the predicted effects of genetic variants on TF binding. It then uses a technique known as message passing to integrate this prior network with gene expression and TF protein–protein interaction data to produce a refined, genotype-specific regulatory network. We used EGRET to infer gene regulatory networks for two blood-derived cell lines and identified genotype-associated, cell line–specific regulatory differences that we subsequently validated using allele-specific expression, chromatin accessibility QTLs, and differential ChIP-seq TF binding. We also inferred EGRET networks for three cell types from each of 119 individuals and identified cell type–specific regulatory differences associated with diseases related to those cell types. EGRET is, to our knowledge, the first method that infers networks reflective of individual genetic variation in a way that provides insight into the genetic regulatory associations driving complex phenotypes.

Protocol: A Multiplexed Reporter Assay to Study Effects of Chromatin Context on DNA Double-Strand Break Repair

DNA double-strand breaks (DSBs) can be repaired through various pathways. Understanding how these pathways are regulated is of great interest for cancer research and optimization of gene editing. The local chromatin environment can affect the balance between repair pathways, but this is still poorly understood. Here we provide a detailed protocol for DSB-TRIP, a technique that utilizes the specific DNA scars left by DSB repair pathways to study pathway usage throughout the genome. DSB-TRIP randomly integrates a repair reporter into many genomic locations, followed by the induction of DSBs in the reporter. Multiplexed sequencing of the resulting scars at all integration sites then reveals the balance between several repair pathways, which can be linked to the local chromatin state of the integration sites. Here we present a step-by-step protocol to perform DSB-TRIP in K562 cells and to analyse the data by a dedicated computational pipeline. We discuss strengths and limitations of the technique, as well as potential additional applications to study DNA repair.

Genomic context sensitivity of insulator function

The specificity of interactions between genomic regulatory elements and potential target genes is influenced by the binding of insulator proteins such as CTCF, which can act as potent enhancer blockers when interposed between an enhancer and a promoter in a reporter assay. But not all CTCF sites genome-wide function as insulator elements, depending on cellular and genomic context. To dissect the influence of genomic context on enhancer blocker activity, we integrated reporter constructs with promoter-only, promoter and enhancer, and enhancer blocker configurations at hundreds of thousands of genomic sites using the Sleeping Beauty transposase. Deconvolution of reporter activity by genomic position reveals distinct expression patterns subject to genomic context, including a compartment of enhancer blocker reporter integrations with robust expression. The high density of integration sites permits quantitative delineation of characteristic genomic context sensitivity profiles and their decomposition into sensitivity to both local and distant DNase I hypersensitive sites. Furthermore, using a single-cell expression approach to test the effect of integrated reporters for differential expression of nearby endogenous genes reveals that CTCF insulator elements do not completely abrogate reporter effects on endogenous gene expression. Collectively, our results lend new insight into genomic regulatory compartmentalization and its influence on the determinants of promoter–enhancer specificity.

New approaches to moderate CRISPR-Cas9 activity: Addressing issues of cellular uptake and endosomal escape

CRISPR-Cas9 is rapidly entering molecular biology and biomedicine as a promising gene-editing tool. A unique feature of CRISPR-Cas9 is a single-guide RNA directing a Cas9 nuclease toward its genomic target. Herein, we highlight new approaches for improving cellular uptake and endosomal escape of CRISPR-Cas9. As opposed to other recently published works, this review is focused on non-viral carriers as a means to facilitate the cellular uptake of CRISPR-Cas9 through endocytosis. The majority of non-viral carriers, such as gold nanoparticles, polymer nanoparticles, lipid nanoparticles, and nanoscale zeolitic imidazole frameworks, is developed with a focus toward optimizing the endosomal escape of CRISPR-Cas9 by taking advantage of the acidic environment in the late endosomes. Among the most broadly used methods for in vitro and ex vivo ribonucleotide protein transfection are electroporation and microinjection. Thus, other delivery formats are warranted for in vivo delivery of CRISPR-Cas9. Herein, we specifically revise the use of peptide and nanoparticle-based systems as platforms for CRISPR-Cas9 delivery in vivo. Finally, we highlight future perspectives of the CRISPR-Cas9 gene-editing tool and the prospects of using non-viral vectors to improve its bioavailability and therapeutic potential.

The FORCE Panel: An All-in-One SNP Marker Set for Confirming Investigative Genetic Genealogy Leads and for General Forensic Applications

The FORensic Capture Enrichment (FORCE) panel is an all-in-one SNP panel for forensic applications. This panel of 5422 markers encompasses common, forensically relevant SNPs (identity, ancestry, phenotype, X- and Y-chromosomal SNPs), a novel set of 3931 autosomal SNPs for extended kinship analysis, and no clinically relevant/disease markers. The FORCE panel was developed as a custom hybridization capture assay utilizing ~20,000 baits to target the selected SNPs. Five non-probative, previously identified World War II (WWII) cases were used to assess the kinship panel. Each case included one bone sample and associated family reference DNA samples. Additionally, seven reference quality samples, two 200-year-old bone samples, and four control DNAs were processed for kit performance and concordance assessments. SNP recovery after capture resulted in a mean of ~99% SNPs exceeding 10X coverage for reference and control samples, and 44.4% SNPs for bone samples. The WWII case results showed that the FORCE panel could predict first to fifth degree relationships with strong statistical support (likelihood ratios over 10,000 and posterior probabilities over 99.99%). To conclude, SNPs will be important for further advances in forensic DNA analysis. The FORCE panel shows promising results and demonstrates the utility of a 5000 SNP panel for forensic applications.

Transposable elements that have recently been mobile in the human genome

Background

Transposable elements (TE) comprise nearly half of the human genome and their insertions have profound effects to human genetic diversification and as well as disease. Despite their abovementioned significance, there is no consensus on the TE subfamilies that remain active in the human genome. In this study, we therefore developed a novel statistical test for recently mobile subfamilies (RMSs), based on patterns of overlap with > 100,000 polymorphic indels.

Results

Our analysis produced a catalogue of 20 high-confidence RMSs, which excludes many false positives in public databases. Intriguingly though, it includes HERV-K, an LTR subfamily previously thought to be extinct. The RMS catalogue is strongly enriched for contributions to germline genetic disorders (P = 1.1e-10), and thus constitutes a valuable resource for diagnosing disorders of unknown aetiology using targeted TE-insertion screens. Remarkably, RMSs are also highly enriched for somatic insertions in diverse cancers (P = 2.8e-17), thus indicating strong correlations between germline and somatic TE mobility. Using CRISPR/Cas9 deletion, we show that an RMS-derived polymorphic TE insertion increased the expression of RPL17, a gene associated with lower survival in liver cancer. More broadly, polymorphic TE insertions from RMSs were enriched near genes with allele-specific expression, suggesting widespread effects on gene regulation.

Conclusions

By using a novel statistical test we have defined a catalogue of 20 recently mobile transposable element subfamilies. We illustrate the gene regulatory potential of RMS-derived polymorphic TE insertions, using CRISPR/Cas9 deletion in vitro on a specific candidate, as well as by genome wide analysis of allele-specific expression. Our study presents novel insights into TE mobility and regulatory potential and provides a key resource for human disease genetics and population history studies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08085-0.

Simple oligonucleotide-based multiplexing of single-cell chromatin accessibility

Microdroplet single-cell ATAC-seq is widely used to measure chromatin accessibility, however, highly scalable and simple sample multiplexing procedures are not available. Here, we present a transposome-assisted single nucleus barcoding approach for ATAC-seq (SNuBar-ATAC) that utilizes a single oligonucleotide adaptor for multiplexing samples during the existing tagmentation step and does not require a pre-labeling procedure. The accuracy and scalability of SNuBar-ATAC was evaluated using cell line mixture experiments. We applied SNuBar-ATAC to investigate treatment-induced chromatin accessibility dynamics by multiplexing 28 mice with lung tumors that received different combinations of chemo, radiation, and targeted immunotherapy. We also applied SNuBar-ATAC to study spatial epigenetic heterogeneity by multiplexing 32 regions from a human breast tissue. Additionally, we show that SNuBar can multiplex single cell ATAC and RNA multiomic assays in cell lines and human breast tissue samples. Our data show that SNuBar is a highly accurate, easy-to-use, and scalable system for multiplexing scATAC-seq and scATAC and RNA co-assay experiments.

The SWI/SNF chromatin remodeling complex helps resolve R-loop-mediated transcription-replication conflicts

ATP-dependent chromatin remodelers are commonly mutated in human cancer. Mammalian SWI/SNF complexes comprise three conserved multisubunit chromatin remodelers (cBAF, ncBAF and PBAF) that share the BRG1 (also known as SMARCA4) subunit responsible for the main ATPase activity. BRG1 is the most frequently mutated Snf2-like ATPase in cancer. In the present study, we have investigated the role of SWI/SNF in genome instability, a hallmark of cancer cells, given its role in transcription, DNA replication and DNA-damage repair. We show that depletion of BRG1 increases R-loops and R-loop-dependent DNA breaks, as well as transcription-replication (T-R) conflicts. BRG1 colocalizes with R-loops and replication fork blocks, as determined by FANCD2 foci, with BRG1 depletion being epistatic to FANCD2 silencing. Our study, extended to other components of SWI/SNF, uncovers a key role of the SWI/SNF complex, in particular cBAF, in helping resolve R-loop-mediated T-R conflicts, thus, unveiling a new mechanism by which chromatin remodeling protects genome integrity.

SMOOTH-seq: single-cell genome sequencing of human cells on a third-generation sequencing platform

There is no effective way to detect structure variations (SVs) and extra-chromosomal circular DNAs (ecDNAs) at single-cell whole-genome level. Here, we develop a novel third-generation sequencing platform-based single-cell whole-genome sequencing (scWGS) method named SMOOTH-seq (single-molecule real-time sequencing of long fragments amplified through transposon insertion). We evaluate the method for detecting CNVs, SVs, and SNVs in human cancer cell lines and a colorectal cancer sample and show that SMOOTH-seq reliably and effectively detects SVs and ecDNAs in individual cells, but shows relatively limited accuracy in detection of CNVs and SNVs. SMOOTH-seq opens a new chapter in scWGS as it generates high fidelity reads of kilobases long.

FaNDOM: Fast nested distance-based seeding of optical maps

Optical mapping (OM) provides single-molecule readouts of fluorescently labeled sequence motifs on long fragments of DNA, resolved to nucleotide-level coordinates. With the advent of microfluidic technologies for analysis of DNA molecules, it is possible to inexpensively generate long OM data ( > 150 kbp) at high coverage. In addition to scaffolding for de novo assembly, OM data can be aligned to a reference genome for identification of genomic structural variants. We introduce FaNDOM (Fast Nested Distance Seeding of Optical Maps)-an optical map alignment tool that greatly reduces the search space of the alignment process. On four benchmark human datasets, FaNDOM was significantly (4-14×) faster than competing tools while maintaining comparable sensitivity and specificity. We used FaNDOM to map variants in three cancer cell lines and identified many biologically interesting structural variants, including deletions, duplications, gene fusions and gene-disrupting rearrangements. FaNDOM is publicly available at https://github.com/jluebeck/FaNDOM.

Hi-C as a molecular rangefinder to examine genomic rearrangements

The mammalian genome is highly packed into the nucleus. Over the past decade, the development of Hi-C has contributed significantly to our understanding of the three-dimensional (3D) chromatin structure, uncovering the principles and functions of higher-order chromatin organizations. Recent studies have repositioned its property in spatial proximity measurement to address challenging problems in genome analyses including genome assembly, haplotype phasing, and the detection of genomic rearrangements. In particular, the power of Hi-C in detecting large-scale structural variations (SVs) in the cancer genome has been demonstrated, which is challenging to be addressed solely with short-read-based whole-genome sequencing analyses. In this review, we first provide a comprehensive view of Hi-C as an intuitive and effective SV detection tool. Then, we introduce recently developed bioinformatics tools utilizing Hi-C to investigate genomic rearrangements. Finally, we discuss the potential application of single-cell Hi-C to address the heterogeneity of genomic rearrangements and sub-population identification in the cancer genome.

Determination of complete chromosomal haplotypes by bulk DNA sequencing

Haplotype phase represents the collective genetic variation between homologous chromosomes and is an essential feature of non-haploid genomes. Here we describe a computational strategy to reliably determine complete whole-chromosome haplotypes using a combination of bulk long-range sequencing and Hi-C sequencing. We demonstrate that this strategy can resolve the haplotypes of parental chromosomes in diploid human genomes with high precision (>99%) and completeness (>98%) and assemble the syntenic structure of rearranged chromosomes in aneuploid cancer genomes at base pair level resolution. Our work enables direct interrogation of chromosome-specific alterations and chromatin reorganization using bulk DNA sequencing.

HICANCER: accurate and complete cancer genome phasing with Hi-C reads

Due to the high complexity of cancer genome, it is too difficult to generate complete cancer genome map which contains the sequence of every DNA molecule until now. Nevertheless, phasing each chromosome in cancer genome into two haplotypes according to germline mutations provides a suboptimal solution to understand cancer genome. However, phasing cancer genome is also a challenging problem, due to the limit in experimental and computational technologies. Hi-C data is widely used in phasing in recent years due to its long-range linkage information and provides an opportunity for solving the problem of phasing cancer genome. The existing Hi-C based phasing methods can not be applied to cancer genome directly, because the somatic mutations in cancer genome such as somatic SNPs, copy number variations and structural variations greatly reduce the correctness and completeness. Here, we propose a new Hi-C based pipeline for phasing cancer genome called HiCancer. HiCancer solves different kinds of somatic mutations and variations, and take advantage of allelic copy number imbalance and linkage disequilibrium to improve the correctness and completeness of phasing. According to our experiments in K562 and KBM-7 cell lines, HiCancer is able to generate very high-quality chromosome-level haplotypes for cancer genome with only Hi-C data.

CRISPR/Cas9 Genome Editing of the Human Topoisomerase IIα Intron 19 5' Splice Site Circumvents Etoposide Resistance in Human Leukemia K562 Cells

An essential function of DNA topoisomerase IIα (TOP2α; 170 kDa, TOP2α/170) is to resolve DNA topologic entanglements during chromosome disjunction by introducing transient DNA double-stranded breaks. TOP2α/170 is an important target for DNA damage-stabilizing anticancer drugs, whose clinical efficacy is compromised by drug resistance often associated with decreased TOP2α/170 expression. We recently demonstrated that an etoposide-resistant K562 clonal subline, K/VP.5, with reduced levels of TOP2α/170, expresses high levels of a novel C-terminal truncated TOP2α isoform (90 kDa, TOP2α/90). TOP2α/90, the translation product of a TOP2α mRNA that retains a processed intron 19 (I19), heterodimerizes with TOP2α/170 and is a resistance determinant through a dominant-negative effect on drug activity. We hypothesized that genome editing to enhance I19 removal would provide a tractable strategy to circumvent acquired TOP2α-mediated drug resistance. To enhance I19 removal in K/VP.5 cells, CRISPR/Cas9 was used to make changes (GAG//GTAA AC →GAG//GTAA GT ) in the TOP2α gene's suboptimal exon 19/intron 19 5' splice site (E19/I19 5' SS). Gene-edited clones were identified by quantitative polymerase chain reaction and verified by sequencing. Characterization of a clone with all TOP2α alleles edited revealed improved I19 removal, decreased TOP2α/90 mRNA/protein, and increased TOP2α/170 mRNA/protein. Sensitivity to etoposide-induced DNA damage (γH2AX, Comet assays) and growth inhibition was restored to levels comparable to those in parental K562 cells. Together, the results indicate that our gene-editing strategy for optimizing the TOP2α E19/I19 5' SS in K/VP.5 cells circumvents resistance to etoposide and other TOP2α-targeted drugs. SIGNIFICANCE STATEMENT: Results presented here indicate that CRISPR/Cas9 gene editing of a suboptimal exon 19/intron 19 5' splice site in the DNA topoisomerase IIα (TOP2α) gene results in circumvention of acquired drug resistance to etoposide and other TOP2α-targeted drugs in a clonal K562 cell line by enhancing removal of intron 19 and thereby decreasing formation of a truncated TOP2α 90 kDa isoform and increasing expression of full-length TOP2α 170 kDa in these resistant cells. Results demonstrate the importance of RNA processing in acquired drug resistance to TOP2α-targeted drugs.

A homology independent sequence replacement strategy in human cells using a CRISPR nuclease

Precision genomic alterations largely rely on homology directed repair (HDR), but targeting without homology using the non-homologous end-joining (NHEJ) pathway has gained attention as a promising alternative. Previous studies demonstrated precise insertions formed by the ligation of donor DNA into a targeted genomic double-strand break in both dividing and non-dividing cells. Here, we demonstrate the use of NHEJ repair to replace genomic segments with donor sequences; we name this method 'Replace' editing (Rational end-joining protocol delivering a targeted sequence exchange). Using CRISPR/Cas9, we create two genomic breaks and ligate a donor sequence in-between. This exchange of a genomic for a donor sequence uses neither microhomology nor homology arms. We target four loci in cell lines and show successful exchange of exons in 16-54% of human cells. Using linear amplification methods and deep sequencing, we quantify the diversity of outcomes following Replace editing and profile the ligated interfaces. The ability to replace exons or other genomic sequences in cells not efficiently modified by HDR holds promise for both basic research and medicine.

Reduced thermal tolerance in a coral carrying CRISPR-induced mutations in the gene for a heat-shock transcription factor

Reef-building corals are keystone species that are threatened by anthropogenic stresses including climate change. To investigate corals' responses to stress and other aspects of their biology, numerous genomic and transcriptomic studies have been performed, generating many hypotheses about the roles of particular genes and molecular pathways. However, it has not generally been possible to test these hypotheses rigorously because of the lack of genetic tools for corals or closely related cnidarians. CRISPR technology seems likely to alleviate this problem. Indeed, we show here that microinjection of single-guide RNA/Cas9 ribonucleoprotein complexes into fertilized eggs of the coral Acropora millepora can produce a sufficiently high frequency of mutations to detect a clear phenotype in the injected generation. Based in part on experiments in a sea-anemone model system, we targeted the gene encoding Heat Shock Transcription Factor 1 (HSF1) and obtained larvae in which >90% of the gene copies were mutant. The mutant larvae survived well at 27 °C but died rapidly at 34 °C, a temperature that did not produce detectable mortality over the duration of the experiment in wild-type (WT) larvae or larvae injected with Cas9 alone. We conclude that HSF1 function (presumably its induction of genes in response to heat stress) plays an important protective role in corals. More broadly, we conclude that CRISPR mutagenesis in corals should allow wide-ranging and rigorous tests of gene function in both larval and adult corals.

Multi-omics analysis reveals the functional transcription and potential translation of enhancers

Enhancer can transcribe RNAs, however, most of them were neglected in traditional RNA-seq analysis workflow. Here, we developed a Pipeline for Enhancer Transcription (PET, http://fun-science.club/PET) for quantifying enhancer RNAs (eRNAs) from RNA-seq. By applying this pipeline on lung cancer samples and cell lines, we showed that the transcribed enhancers are enriched with histone marks and transcription factor motifs (JUNB, Hand1-Tcf3 and GATA4). By training a machine learning model, we demonstrate that enhancers can predict prognosis better than their nearby genes. Integrating the Hi-C, ChIP-seq and RNA-seq data, we observe that transcribed enhancers associate with cancer hallmarks or oncogenes, among which LcsMYC-1 (Lung cancer-specific MYC eRNA-1) potentially supports MYC expression. Surprisingly, a significant proportion of transcribed enhancers contain small protein-coding open reading frames (sORFs) and can be translated into microproteins. Our study provides a computational method for eRNA quantification and deepens our understandings of the DNA, RNA and protein nature of enhancers.

MethHaplo: combining allele-specific DNA methylation and SNPs for haplotype region identification

Background

DNA methylation is an important epigenetic modification that plays a critical role in most eukaryotic organisms. Parental alleles in haploid genomes may exhibit different methylation patterns, which can lead to different phenotypes and even different therapeutic and drug responses to diseases. However, to our knowledge, no software is available for the identification of DNA methylation haplotype regions with combined allele-specific DNA methylation, single nucleotide polymorphisms (SNPs) and high-throughput chromosome conformation capture (Hi-C) data.

Results

In this paper, we developed a new method, MethHaplo, that identify DNA methylation haplotype regions with allele-specific DNA methylation and SNPs from whole-genome bisulfite sequencing (WGBS) data. Our results showed that methylation haplotype regions were ten times longer than haplotypes with SNPs only. When we integrate WGBS and Hi-C data, MethHaplo could call even longer haplotypes.

Conclusions

This study illustrates the usefulness of methylation haplotypes. By constructing methylation haplotypes for various cell lines, we provide a clearer picture of the effect of DNA methylation on gene expression, histone modification and three-dimensional chromosome structure at the haplotype level. Our method could benefit the study of parental inheritance-related disease and hybrid vigor in agriculture.

Improved haplotype inference by exploiting long-range linking and allelic imbalance in RNA-seq datasets

Haplotype reconstruction of distant genetic variants remains an unsolved problem due to the short-read length of common sequencing data. Here, we introduce HapTree-X, a probabilistic framework that utilizes latent long-range information to reconstruct unspecified haplotypes in diploid and polyploid organisms. It introduces the observation that differential allele-specific expression can link genetic variants from the same physical chromosome, thus even enabling using reads that cover only individual variants. We demonstrate HapTree-X's feasibility on in-house sequenced Genome in a Bottle RNA-seq and various whole exome, genome, and 10X Genomics datasets. HapTree-X produces more complete phases (up to 25%), even in clinically important genes, and phases more variants than other methods while maintaining similar or higher accuracy and being up to 10×  faster than other tools. The advantage of HapTree-X's ability to use multiple lines of evidence, as well as to phase polyploid genomes in a single integrative framework, substantially grows as the amount of diverse data increases.