High-Resolution Genomic Profiling of Liver Cancer Links Etiology With Mutation and Epigenetic Signatures.

BACKGROUND & AIMS: Hepatocellular carcinoma (HCC) is a model of a diverse spectrum of cancers because it is induced by well-known etiologies, mainly hepatitis C virus (HCV) and hepatitis B virus. Here, we aimed to identify HCV-specific mutational signatures and explored the link between the HCV-related regional variation in mutations rates and HCV-induced alterations in genome-wide chromatin organization. METHODS: To identify an HCV-specific mutational signature in HCC, we performed high-resolution targeted sequencing to detect passenger mutations on 64 HCC samples from 3 etiology groups: hepatitis B virus, HCV, or other. To explore the link between the genomic signature and genome-wide chromatin organization we performed chromatin immunoprecipitation sequencing for the transcriptionally permissive H3K4Me3, H3K9Ac, and suppressive H3K9Me3 modifications after HCV infection. RESULTS: Regional variation in mutation rate analysis showed significant etiology-dependent regional mutation rates in 12 genes: LRP2, KRT84, TMEM132B, DOCK2, DMD, INADL, JAK2, DNAH6, MTMR9, ATM, SLX4, and ARSD. We found an enrichment of C->T transversion mutations in the HCV-associated HCC cases. Furthermore, these cases showed regional variation in mutation rates associated with genomic intervals in which HCV infection dictated epigenetic alterations. This signature may be related to the HCV-induced decreased expression of genes encoding key enzymes in the base excision repair pathway. CONCLUSIONS: We identified novel distinct HCV etiology-dependent mutation signatures in HCC associated with HCV-induced alterations in histone modification. This study presents a link between cancer-causing mutagenesis and the increased predisposition to liver cancer in chronic HCV-infected individuals, and unveils novel etiology-specific mechanisms leading to HCC and cancer in general.

HCV Infection Increases the Expression of ACE2 Receptor, Leading to Enhanced Entry of Both HCV and SARS-CoV-2 into Hepatocytes and a Coinfection State.

Recent studies suggest the enhancement of liver injury in COVID-19 patients infected with Hepatitis C virus (HCV). Hepatocytes express low levels of angiotensin-converting enzyme 2 (ACE2), the SARS-CoV-2 entry receptor, raising the possibility of HCV-SARS-CoV-2 coinfection in the liver. This work aimed to explore whether HCV and SARS-CoV-2 coinfect hepatocytes and the interplay between these viruses. We demonstrate that SARS-CoV-2 coinfects HCV-infected Huh7.5 (Huh7.5HCV) cells. Both viruses replicated efficiently in the coinfected cells, with HCV replication enhanced in coinfected compared to HCV-mono-infected cells. Strikingly, Huh7.5HCV cells were eight fold more susceptible to SARS-CoV-2 pseudoviruses than naive Huh7.5 cells, suggesting enhanced SARS-CoV-2 entry into HCV-preinfected hepatocytes. In addition, we observed increased binding of spike receptor-binding domain (RBD) protein to Huh7.5HCV cells, as well as enhanced cell-to-cell fusion of Huh7.5HCV cells with spike-expressing Huh7.5 cells. We explored the mechanism of enhanced SARS-CoV-2 entry and identified an increased ACE2 mRNA and protein levels in Huh7.5HCV cells, primary hepatocytes, and in data from infected liver biopsies obtained from database. Importantly, higher expression of ACE2 increased HCV infection by enhancing its binding to the host cell, underscoring its role in the HCV life cycle as well. Transcriptome analysis revealed that shared host signaling pathways were induced in HCV-SARS-CoV-2 coinfection. This study revealed complex interactions between HCV and SARS-CoV-2 infections in hepatocytes, which may lead to the increased liver damage recently reported in HCV-positive COVID-19 patients. IMPORTANCE Here, we provide the first experimental evidence for the coexistence of SARS-CoV-2 infection with HCV, and the interplay between them. The study revealed a complex relationship of enhancement between the two viruses, where HCV infection increased the expression of the SARS-CoV-2 entry receptor ACE2, thus facilitating SARS-CoV-2 entry, and potentially, also HCV entry. Thereafter, SARS-CoV-2 infection enhanced HCV replication in hepatocytes. This study may explain the aggravation of liver damage that was recently reported in COVID-19 patients with HCV coinfection and suggests preinfection with HCV as a risk factor for severe COVID-19. Moreover, it highlights the possible importance of HCV treatment for coinfected patients. In a broader view, these findings emphasize the importance of identifying coinfecting pathogens that increase the risk of SARS-CoV-2 infection and that may accelerate COVID-19-related co-morbidities.

Ig Enhancers Increase RNA Polymerase II Stalling at Somatic Hypermutation Target Sequences.

Somatic hypermutation (SHM) drives the genetic diversity of Ig genes in activated B cells and supports the generation of Abs with increased affinity for Ag. SHM is targeted to Ig genes by their enhancers (diversification activators [DIVACs]), but how the enhancers mediate this activity is unknown. We show using chicken DT40 B cells that highly active DIVACs increase the phosphorylation of RNA polymerase II (Pol II) and Pol II occupancy in the mutating gene with little or no accompanying increase in elongation-competent Pol II or production of full-length transcripts, indicating accumulation of stalled Pol II. DIVAC has similar effect also in human Ramos Burkitt lymphoma cells. The DIVAC-induced stalling is weakly associated with an increase in the detection of ssDNA bubbles in the mutating target gene. We did not find evidence for antisense transcription, or that DIVAC functions by altering levels of H3K27ac or the histone variant H3.3 in the mutating gene. These findings argue for a connection between Pol II stalling and cis-acting targeting elements in the context of SHM and thus define a mechanistic basis for locus-specific targeting of SHM in the genome. Our results suggest that DIVAC elements render the target gene a suitable platform for AID-mediated mutation without a requirement for increasing transcriptional output.

Image-based pooled whole-genome CRISPRi screening for subcellular phenotypes.

Genome-wide CRISPR screens have transformed our ability to systematically interrogate human gene function, but are currently limited to a subset of cellular phenotypes. We report a novel pooled screening approach for a wider range of cellular and subtle subcellular phenotypes. Machine learning and convolutional neural network models are trained on the subcellular phenotype to be queried. Genome-wide screening then utilizes cells stably expressing dCas9-KRAB (CRISPRi), photoactivatable fluorescent protein (PA-mCherry), and a lentiviral guide RNA (gRNA) pool. Cells are screened by using microscopy and classified by artificial intelligence (AI) algorithms, which precisely identify the genetically altered phenotype. Cells with the phenotype of interest are photoactivated and isolated via flow cytometry, and the gRNAs are identified by sequencing. A proof-of-concept screen accurately identified PINK1 as essential for Parkin recruitment to mitochondria. A genome-wide screen identified factors mediating TFEB relocation from the nucleus to the cytosol upon prolonged starvation. Twenty-one of the 64 hits called by the neural network model were independently validated, revealing new effectors of TFEB subcellular localization. This approach, AI-photoswitchable screening (AI-PS), offers a novel screening platform capable of classifying a broad range of mammalian subcellular morphologies, an approach largely unattainable with current methodologies at genome-wide scale.

ATM and PRDM9 regulate SPO11-bound recombination intermediates during meiosis.

Meiotic recombination is initiated by SPO11-induced double-strand breaks (DSBs). In most mammals, the methyltransferase PRDM9 guides SPO11 targeting, and the ATM kinase controls meiotic DSB numbers. Following MRE11 nuclease removal of SPO11, the DSB is resected and loaded with DMC1 filaments for homolog invasion. Here, we demonstrate the direct detection of meiotic DSBs and resection using END-seq on mouse spermatocytes with low sample input. We find that DMC1 limits both minimum and maximum resection lengths, whereas 53BP1, BRCA1 and EXO1 play surprisingly minimal roles. Through enzymatic modifications to END-seq, we identify a SPO11-bound meiotic recombination intermediate (SPO11-RI) present at all hotspots. We propose that SPO11-RI forms because chromatin-bound PRDM9 asymmetrically blocks MRE11 from releasing SPO11. In Atm-/- spermatocytes, trapped SPO11 cleavage complexes accumulate due to defective MRE11 initiation of resection. Thus, in addition to governing SPO11 breakage, ATM and PRDM9 are critical local regulators of mammalian SPO11 processing.

53BP1 Enforces Distinct Pre- and Post-resection Blocks on Homologous Recombination.

53BP1 activity drives genome instability and lethality in BRCA1-deficient mice by inhibiting homologous recombination (HR). The anti-recombinogenic functions of 53BP1 require phosphorylation-dependent interactions with PTIP and RIF1/shieldin effector complexes. While RIF1/shieldin blocks 5'-3' nucleolytic processing of DNA ends, it remains unclear how PTIP antagonizes HR. Here, we show that mutation of the PTIP interaction site in 53BP1 (S25A) allows sufficient DNA2-dependent end resection to rescue the lethality of BRCA1Δ11 mice, despite increasing RIF1 "end-blocking" at DNA damage sites. However, double-mutant cells fail to complete HR, as excessive shieldin activity also inhibits RNF168-mediated loading of PALB2/RAD51. As a result, BRCA1Δ1153BP1S25A mice exhibit hallmark features of HR insufficiency, including premature aging and hypersensitivity to PARPi. Disruption of shieldin or forced targeting of PALB2 to ssDNA in BRCA1D1153BP1S25A cells restores RNF168 recruitment, RAD51 nucleofilament formation, and PARPi resistance. Our study therefore reveals a critical function of shieldin post-resection that limits the loading of RAD51.

Topologically Associated Domains Delineate Susceptibility to Somatic Hypermutation.

Somatic hypermutation (SHM) introduces point mutations into immunoglobulin (Ig) genes but also causes mutations in other parts of the genome. We have used lentiviral SHM reporter vectors to identify regions of the genome that are susceptible ("hot") and resistant ("cold") to SHM, revealing that SHM susceptibility and resistance are often properties of entire topologically associated domains (TADs). Comparison of hot and cold TADs reveals that while levels of transcription are equivalent, hot TADs are enriched for the cohesin loader NIPBL, super-enhancers, markers of paused/stalled RNA polymerase 2, and multiple important B cell transcription factors. We demonstrate that at least some hot TADs contain enhancers that possess SHM targeting activity and that insertion of a strong Ig SHM-targeting element into a cold TAD renders it hot. Our findings lead to a model for SHM susceptibility involving the cooperative action of cis-acting SHM targeting elements and the dynamic and architectural properties of TADs.

Intra-Vκ Cluster Recombination Shapes the Ig Kappa Locus Repertoire.

During V(D)J recombination, RAG proteins introduce DNA double-strand breaks (DSBs) at recombination signal sequences (RSSs) that contain either 12- or 23-nt spacer regions. Coordinated 12/23 cleavage predicts that DSBs at variable (V) gene segments should equal the level of breakage at joining (J) segments. Contrary to this, here we report abundant RAG-dependent DSBs at multiple Vκ gene segments independent of V-J rearrangement. We find that a large fraction of Vκ gene segments are flanked not only by a bone-fide 12 spacer but also an overlapping, 23-spacer flipped RSS. These compatible pairs of RSSs mediate recombination and deletion inside the Vκ cluster even in the complete absence of Jκ gene segments and support a V(D)J recombination center (RC) independent of the conventional Jκ-centered RC. We propose an improved model of Vκ-Jκ repertoire formation by incorporating these surprisingly frequent, evolutionarily conserved intra-Vκ cluster recombination events.

Topoisomerase II-Induced Chromosome Breakage and Translocation Is Determined by Chromosome Architecture and Transcriptional Activity.

Topoisomerase II (TOP2) relieves torsional stress by forming transient cleavage complex intermediates (TOP2ccs) that contain TOP2-linked DNA breaks (DSBs). While TOP2ccs are normally reversible, they can be "trapped" by chemotherapeutic drugs such as etoposide and subsequently converted into irreversible TOP2-linked DSBs. Here, we have quantified etoposide-induced trapping of TOP2ccs, their conversion into irreversible TOP2-linked DSBs, and their processing during DNA repair genome-wide, as a function of time. We find that while TOP2 chromatin localization and trapping is independent of transcription, it requires pre-existing binding of cohesin to DNA. In contrast, the conversion of trapped TOP2ccs to irreversible DSBs during DNA repair is accelerated 2-fold at transcribed loci relative to non-transcribed loci. This conversion is dependent on proteasomal degradation and TDP2 phosphodiesterase activity. Quantitative modeling shows that only two features of pre-existing chromatin structure-namely, cohesin binding and transcriptional activity-can be used to predict the kinetics of TOP2-induced DSBs.

Dual Roles of Poly(dA:dT) Tracts in Replication Initiation and Fork Collapse.

Replication origins, fragile sites, and rDNA have been implicated as sources of chromosomal instability. However, the defining genomic features of replication origins and fragile sites are among the least understood elements of eukaryote genomes. Here, we map sites of replication initiation and breakage in primary cells at high resolution. We find that replication initiates between transcribed genes within nucleosome-depleted structures established by long asymmetrical poly(dA:dT) tracts flanking the initiation site. Paradoxically, long (>20 bp) (dA:dT) tracts are also preferential sites of polar replication fork stalling and collapse within early-replicating fragile sites (ERFSs) and late-replicating common fragile sites (CFSs) and at the rDNA replication fork barrier. Poly(dA:dT) sequences are fragile because long single-strand poly(dA) stretches at the replication fork are unprotected by the replication protein A (RPA). We propose that the evolutionary expansion of poly(dA:dT) tracts in eukaryotic genomes promotes replication initiation, but at the cost of chromosome fragility.

The Energetics and Physiological Impact of Cohesin Extrusion.

Cohesin extrusion is thought to play a central role in establishing the architecture of mammalian genomes. However, extrusion has not been visualized in vivo, and thus, its functional impact and energetics are unknown. Using ultra-deep Hi-C, we show that loop domains form by a process that requires cohesin ATPases. Once formed, however, loops and compartments are maintained for hours without energy input. Strikingly, without ATP, we observe the emergence of hundreds of CTCF-independent loops that link regulatory DNA. We also identify architectural "stripes," where a loop anchor interacts with entire domains at high frequency. Stripes often tether super-enhancers to cognate promoters, and in B cells, they facilitate Igh transcription and recombination. Stripe anchors represent major hotspots for topoisomerase-mediated lesions, which promote chromosomal translocations and cancer. In plasmacytomas, stripes can deregulate Igh-translocated oncogenes. We propose that higher organisms have coopted cohesin extrusion to enhance transcription and recombination, with implications for tumor development.

Genome Organization Drives Chromosome Fragility.

In this study, we show that evolutionarily conserved chromosome loop anchors bound by CCCTC-binding factor (CTCF) and cohesin are vulnerable to DNA double strand breaks (DSBs) mediated by topoisomerase 2B (TOP2B). Polymorphisms in the genome that redistribute CTCF/cohesin occupancy rewire DNA cleavage sites to novel loop anchors. While transcription- and replication-coupled genomic rearrangements have been well documented, we demonstrate that DSBs formed at loop anchors are largely transcription-, replication-, and cell-type-independent. DSBs are continuously formed throughout interphase, are enriched on both sides of strong topological domain borders, and frequently occur at breakpoint clusters commonly translocated in cancer. Thus, loop anchors serve as fragile sites that generate DSBs and chromosomal rearrangements. VIDEO ABSTRACT.

RAG1 targeting in the genome is dominated by chromatin interactions mediated by the non-core regions of RAG1 and RAG2.

The RAG1/RAG2 endonuclease initiates V(D)J recombination at antigen receptor loci but also binds to thousands of places outside of these loci. RAG2 localizes directly to lysine 4 trimethylated histone 3 (H3K4me3) through a plant homeodomain (PHD) finger. The relative contribution of RAG2-dependent and RAG1-intrinsic mechanisms in determining RAG1 binding patterns is not known. Through analysis of deep RAG1 ChIP-seq data, we provide a quantitative description of the forces underlying genome-wide targeting of RAG1. Surprisingly, sequence-specific DNA binding contributes minimally to RAG1 targeting outside of antigen receptor loci. Instead, RAG1 binding is driven by two distinct modes of interaction with chromatin: the first is driven by H3K4me3, promoter-focused and dependent on the RAG2 PHD, and the second is defined by H3K27Ac, enhancer-focused and dependent on 'non-core' portions of RAG1. Based on this and additional chromatin and genomic features, we formulated a predictive model of RAG1 targeting to the genome. RAG1 binding sites predicted by our model correlate well with observed patterns of RAG1-mediated breaks in human pro-B acute lymphoblastic leukemia. Overall, this study provides an integrative model for RAG1 genome-wide binding and off-target activity and reveals a novel role for the RAG1 non-core region in RAG1 targeting.

Bcl6 Is Required for Somatic Hypermutation and Gene Conversion in Chicken DT40 Cells.

The activation induced cytosine deaminase (AID) mediates diversification of B cell immunoglobulin genes by the three distinct yet related processes of somatic hypermutation (SHM), class switch recombination (CSR), and gene conversion (GCV). SHM occurs in germinal center B cells, and the transcription factor Bcl6 is a key regulator of the germinal center B cell gene expression program, including expression of AID. To test the hypothesis that Bcl6 function is important for the process of SHM, we compared WT chicken DT40 B cells, which constitutively perform SHM/GCV, to their Bcl6-deficient counterparts. We found that Bcl6-deficient DT40 cells were unable to perform SHM and GCV despite enforced high level expression of AID and substantial levels of AID in the nucleus of the cells. To gain mechanistic insight into the GCV/SHM dependency on Bcl6, transcriptional features of a highly expressed SHM target gene were analyzed in Bcl6-sufficient and -deficient DT40 cells. No defect was observed in the accumulation of single stranded DNA in the target gene as a result of Bcl6 deficiency. In contrast, association of Spt5, an RNA polymerase II (Pol II) and AID binding factor, was strongly reduced at the target gene body relative to the transcription start site in Bcl6-deficient cells as compared to WT cells. However, partial reconstitution of Bcl6 function substantially reconstituted Spt5 association with the target gene body but did not restore detectable SHM. Our observations suggest that in the absence of Bcl6, Spt5 fails to associate efficiently with Pol II at SHM targets, perhaps precluding robust AID action on the SHM target DNA. Our data also suggest, however, that Spt5 binding is not sufficient for SHM of a target gene even in DT40 cells with strong expression of AID.

Estimate of within population incremental selection through branch imbalance in lineage trees.

Incremental selection within a population, defined as limited fitness changes following mutation, is an important aspect of many evolutionary processes. Strongly advantageous or deleterious mutations are detected using the synonymous to non-synonymous mutations ratio. However, there are currently no precise methods to estimate incremental selection. We here provide for the first time such a detailed method and show its precision in multiple cases of micro-evolution. The proposed method is a novel mixed lineage tree/sequence based method to detect within population selection as defined by the effect of mutations on the average number of offspring. Specifically, we propose to measure the log of the ratio between the number of leaves in lineage trees branches following synonymous and non-synonymous mutations. The method requires a high enough number of sequences, and a large enough number of independent mutations. It assumes that all mutations are independent events. It does not require of a baseline model and is practically not affected by sampling biases. We show the method's wide applicability by testing it on multiple cases of micro-evolution. We show that it can detect genes and inter-genic regions using the selection rate and detect selection pressures in viral proteins and in the immune response to pathogens.

RAG Represents a Widespread Threat to the Lymphocyte Genome.

The RAG1 endonuclease, together with its cofactor RAG2, is essential for V(D)J recombination but is a potent threat to genome stability. The sources of RAG1 mis-targeting and the mechanisms that have evolved to suppress it are poorly understood. Here, we report that RAG1 associates with chromatin at thousands of active promoters and enhancers in the genome of developing lymphocytes. The mouse and human genomes appear to have responded by reducing the abundance of "cryptic" recombination signals near RAG1 binding sites. This depletion operates specifically on the RSS heptamer, whereas nonamers are enriched at RAG1 binding sites. Reversing this RAG-driven depletion of cleavage sites by insertion of strong recombination signals creates an ectopic hub of RAG-mediated V(D)J recombination and chromosomal translocations. Our findings delineate rules governing RAG binding in the genome, identify areas at risk of RAG-mediated damage, and highlight the evolutionary struggle to accommodate programmed DNA damage in developing lymphocytes.

Genomic landscape of cutaneous T cell lymphoma.

Cutaneous T cell lymphoma (CTCL) is a non-Hodgkin lymphoma of skin-homing T lymphocytes. We performed exome and whole-genome DNA sequencing and RNA sequencing on purified CTCL and matched normal cells. The results implicate mutations in 17 genes in CTCL pathogenesis, including genes involved in T cell activation and apoptosis, NF-κB signaling, chromatin remodeling and DNA damage response. CTCL is distinctive in that somatic copy number variants (SCNVs) comprise 92% of all driver mutations (mean of 11.8 pathogenic SCNVs versus 1.0 somatic single-nucleotide variant per CTCL). These findings have implications for new therapeutics.

Viral CD8 T cell epitope nucleotide composition shows evidence of short- and long-term evolutionary strategies.

Viral epitopes have a distinct codon usage that reflects their dual role in infection and immunity. On the one hand, epitopes are part of proteins important to viral function; on the other hand, they are targets of the immune response. Studies of selection are most commonly based on changes of amino acid and seen through the accumulation of non-synonymous mutations. An independent measure of selection is the codon usage and underlying changeability of the nucleotide sequences. We here use multiple tools and a large-scale analysis of viral genomes to demonstrate that viral epitopes have a distinct codon usage and that this codon usage reflects distinct short- and long-term types of selection during viral evolution. We show that CD8(+) T cell epitopes are encoded by codons more distant from stop codons and more changeable than codons outside epitopes. This biased codon usage reflects the viral population toggling back and forth from a wild-type sequence to an escape mode, which enable them to avoid immune detection when needed, and go back to the functionally favorable form when the threat is removed (i.e., in a new host).

Viral proteome size and CD8+ T cell epitope density are correlated: the effect of complexity on selection.

The relation between the complexity of organisms and proteins and their evolution rates has been discussed in the context of multiple generic models. The main robust claim from most such models is the negative relation between complexity and the accumulation rate of mutations. Viruses accumulate escape mutations in their epitopes to avoid detection and destruction of their host cell by CD8+ T cells. The extreme regime of immune escape, namely, strong selection and high mutation rate, provide an opportunity to extend and validate the existing models of relation between complexity and evolution rate as proposed by Fisher and Kimura. Using epitope prediction algorithms to compute the epitopes presented on the most frequent human HLA alleles in over 100 fully sequenced human viruses, and over 900 non-human viruses, we here study the correlation between viruses/proteins complexity (as measured by the number of proteins in the virus and the length of each protein, respectively) and the rate of accumulation of escape mutation. The latter is evaluated by measuring the normalized epitope density of viral proteins. If the virus/protein complexity prevents the accumulation of escape mutations, the epitope density is expected to be positively correlated with both the number of proteins in the virus and the length of proteins. We show that such correlations are indeed observed for most human viruses. For non-human viruses the correlations were much less significant, indicating that the correlation is indeed induced by human HLA molecules.