Cis-regulatory effect of HPV integration is constrained by host chromatin architecture in cervical cancers.

Human papillomavirus (HPV) infections are the primary drivers of cervical cancers, and often HPV DNA gets integrated into the host genome. Although the oncogenic impact of HPV encoded genes is relatively well known, the cis-regulatory effect of integrated HPV DNA on host chromatin structure and gene regulation remains less understood. We investigated genome-wide patterns of HPV integrations and associated host gene expression changes in the context of host chromatin states and topologically associating domains (TADs). HPV integrations were significantly enriched in active chromatin regions and depleted in inactive ones. Interestingly, regardless of chromatin state, genomic regions flanking HPV integrations showed transcriptional upregulation. Nevertheless, upregulation (both local and long-range) was mostly confined to TADs with integration, but not affecting adjacent TADs. Few TADs showed recurrent integrations associated with overexpression of oncogenes within them (e.g. MYC, PVT1, TP63 and ERBB2) regardless of proximity. Hi-C and 4C-seq analyses in cervical cancer cell line (HeLa) demonstrated chromatin looping interactions between integrated HPV and MYC/PVT1 regions (~ 500 kb apart), leading to allele-specific overexpression. Based on these, we propose HPV integrations can trigger multimodal oncogenic activation to promote cancer progression.

An androgen receptor regulated gene score is associated with epithelial to mesenchymal transition features in triple negative breast cancers.

BACKGROUND: Androgen receptor (AR) is considered a marker of better prognosis in hormone receptor positive breast cancers (BC), however, its role in triple negative breast cancer (TNBC) is controversial. This may be attributed to intrinsic molecular differences or scoring methods for AR positivity. We derived AR regulated gene score and examined its utility in BC subtypes. METHODS: AR regulated genes were derived by applying a bioinformatic pipeline on publicly available microarray data sets of AR+ BC cell lines and gene score was calculated as average expression of six AR regulated genes. Tumors were divided into AR high and low based on gene score and associations with clinical parameters, circulating androgens, survival and epithelial to mesenchymal transition (EMT) markers were examined, further evaluated in invitro models and public datasets. RESULTS: 53% (133/249) tumors were classified as AR gene score high and were associated with significantly better clinical parameters, disease-free survival (86.13 vs 72.69 months, log rank p = 0.032) when compared to AR low tumors. 36% of TNBC (N = 66) were AR gene score high with higher expression of EMT markers (p = 0.024) and had high intratumoral levels of 5α-reductase, enzyme involved in intracrine androgen metabolism. In MDA-MB-453 treated with dihydrotestosterone, SLUG expression increased, E-cadherin decreased with increase in migration and these changes were reversed with bicalutamide. Similar results were obtained in public datasets. CONCLUSION: Deciphering the role of AR in BC is difficult based on AR protein levels alone. Our results support the context dependent function of AR in driving better prognosis in ER positive tumors and EMT features in TNBC tumors.

Prognostic association of immunoproteasome expression in solid tumours is governed by the immediate immune environment.

Induction of immunoproteasome (IP) expression in tumour cells can enhance antigen presentation and immunogenicity. Recently, the overexpression of IP genes has been associated with better prognosis and response to immune checkpoint blockade (ICB) therapies in melanoma. However, the extent of this association in other solid tumours and how that is influenced by tumour cell-intrinsic and cell-extrinsic factors remain unclear. Here, we address this by exploring the gene expression patterns from available bulk and single-cell transcriptomic data of primary tumours. We find that tumours with high-IP expression exhibit cytotoxic immune cell infiltration and upregulation of IFN-γ and TNF-α pathways in tumour cells. However, the association of IP expression with overall survival (TCGA cohort) and response to ICB therapy (non-TCGA cohorts) is tumour-type specific (better in non-small-cell lung, breast, bladder and thymus; and worse in glioma and renal) and is greatly influenced by pro- or antitumourigenic immune cell infiltration patterns. This emphasises the need for considering immune cell infiltration patterns, along with IP expression, as a prognostic biomarker to predict overall survival or response to ICB therapies in solid tumours, besides melanoma.

Functional analysis of structural variants in single cells using Strand-seq.

Somatic structural variants (SVs) are widespread in cancer, but their impact on disease evolution is understudied due to a lack of methods to directly characterize their functional consequences. We present a computational method, scNOVA, which uses Strand-seq to perform haplotype-aware integration of SV discovery and molecular phenotyping in single cells by using nucleosome occupancy to infer gene expression as a readout. Application to leukemias and cell lines identifies local effects of copy-balanced rearrangements on gene deregulation, and consequences of SVs on aberrant signaling pathways in subclones. We discovered distinct SV subclones with dysregulated Wnt signaling in a chronic lymphocytic leukemia patient. We further uncovered the consequences of subclonal chromothripsis in T cell acute lymphoblastic leukemia, which revealed c-Myb activation, enrichment of a primitive cell state and informed successful targeting of the subclone in cell culture, using a Notch inhibitor. By directly linking SVs to their functional effects, scNOVA enables systematic single-cell multiomic studies of structural variation in heterogeneous cell populations.

Differential chromatin accessibility landscape of gain-of-function mutant p53 tumours.

BACKGROUND: Mutations in TP53 not only affect its tumour suppressor activity but also exerts oncogenic gain-of-function activity. While the genome-wide mutant p53 binding sites have been identified in cancer cell lines, the chromatin accessibility landscape driven by mutant p53 in primary tumours is unknown. Here, we leveraged the chromatin accessibility data of primary tumours from The Cancer Genome Atlas (TCGA) to identify differentially accessible regions in mutant p53 tumours compared to wild-type p53 tumours, especially in breast and colon cancers. RESULTS: We identified 1587 lost and 984 gained accessible chromatin regions in breast, and 1143 lost and 640 gained regions in colon cancers. However, only less than half of those regions in both cancer types contain sequence motifs for wild-type or mutant p53 binding. Whereas, the remaining showed enrichment for master transcriptional regulators, such as FOX-Family TFs and NF-kB in lost and SMAD and KLF TFs in gained regions of breast. In colon, ATF3 and FOS/JUN TFs were enriched in lost, and CDX family TFs and HNF4A in gained regions. By integrating the gene expression data, we identified known and novel target genes regulated by the mutant p53. CONCLUSION: This study reveals the direct and indirect mechanisms by which gain-of-function mutant p53 targets the chromatin and subsequent gene expression patterns in a tumour-type specific manner. This furthers our understanding of the impact of mutant p53 in cancer development.

Variable interplay of UV-induced DNA damage and repair at transcription factor binding sites.

An abnormally high rate of UV-light related mutations appears at transcription factor binding sites (TFBS) across melanomas. The binding of transcription factors (TFs) to the DNA impairs the repair of UV-induced lesions and certain TFs have been shown to increase the rate of generation of these lesions at their binding sites. However, the precise contribution of these two elements to the increase in mutation rate at TFBS in these malignant cells is not understood. Here, exploiting nucleotide-resolution data, we computed the rate of formation and repair of UV-lesions within the binding sites of TFs of different families. We observed, at certain dipyrimidine positions within the binding site of TFs in the Tryptophan Cluster family, an increased rate of formation of UV-induced lesions, corroborating previous studies. Nevertheless, across most families of TFs, the observed increased mutation rate within the entire DNA region covered by the protein results from the decreased repair efficiency. While the rate of mutations across all TFBS does not agree with the amount of UV-induced lesions observed immediately after UV exposure, it strongly agrees with that observed after 48 h. This corroborates the determinant role of the impaired repair in the observed increase of mutation rate.

Exome sequencing identifies ARID2 as a novel tumor suppressor in early-onset sporadic rectal cancer.

Early-onset sporadic rectal cancer (EOSRC) is a unique and predominant colorectal cancer (CRC) subtype in India. In order to understand the tumorigenic process in EOSRC, we performed whole-exome sequencing of 47 microsatellite stable EOSRC samples. Signature 1 was the predominant mutational signature in EOSRC, as previously shown in other CRC exome studies. More importantly, we identified TP53, KRAS, APC, PIK3R1, SMAD4 and ZNF880 as significantly mutated (q < 0.1) and ARID1A and ARID2 as near-significantly mutated (restricted hypothesis testing; q < 0.1) candidate drivers. Unlike the other candidates, the tumorigenic potential of ARID2, encoding a component of the SWI/SNF chromatin remodeling complex, is largely unexplored in CRC. shRNA-mediated ARID2 knockdown performed in different CRC cell lines resulted in significant alterations in transcript levels of cancer-related target genes. More importantly, ARID2 knockdown promoted several tumorigenic features including cell viability, proliferation, ability to override contact inhibition of growth, and migration besides significantly increasing tumor formation ability in nude mice. The observed gain in tumorigenic features was rescued upon ectopic expression of wild type but not mutant ARID2. Analyses of the TCGA pan-cancer dataset revealed several modes of ARID2 inactivation and of the CRC dataset revealed poorer survival in patients with ARID2 alterations. We therefore propose ARID2 as a novel tumor suppressor in CRC.

A rare variant of African ancestry activates 8q24 lncRNA hub by modulating cancer associated enhancer.

Genetic variation at the 8q24 locus is linked with the greater susceptibility to prostate cancer in men of African ancestry. One such African ancestry specific rare variant, rs72725854 (A>G/T) (~6% allele frequency) has been associated with a ~2-fold increase in prostate cancer risk. However, the functional relevance of this variant is unknown. Here we show that the variant rs72725854 is present in a prostate cancer-specific enhancer at 8q24 locus. Chromatin-conformation capture and dCas9 mediated enhancer blocking establish a direct regulatory link between this enhancer and lncRNAs PCAT1, PRNCR1 and PVT1. The risk allele ('T') is associated with higher expression of PCAT1, PVT1 and c-myc in prostate tumors. Further, enhancer with the risk allele gains response to androgen stimulation by recruiting the transcription factor SPDEF whereas, non-risk alleles remain non-responsive. Elevated expression of these lncRNAs and c-myc in risk allele carriers may explain their greater susceptibility to prostate cancer.

Twist1 induces chromosomal instability (CIN) in colorectal cancer cells.

Twist1 is a basic helix-loop-helix transcription factor, essential during early development in mammals. While Twist1 induces epithelial-to-mesenchymal transition (EMT), here we show that Twist1 overexpression enhances nuclear and mitotic aberrations. This is accompanied by an increase in whole chromosomal copy number gains and losses, underscoring the role of Twist1 in inducing chromosomal instability (CIN) in colorectal cancer cells. Array comparative genomic hybridization (array CGH) analysis further shows sub-chromosomal deletions, consistent with an increased frequency of DNA double strand breaks (DSBs). Remarkably, Twist1 overexpression downmodulates key cell cycle checkpoint factors-Bub1, BubR1, Mad1 and Mad2-that regulate CIN. Mathematical simulations using the RACIPE tool show a negative correlation of Twist1 with E-cadherin and BubR1. Data analyses of gene expression profiles of patient samples from The Cancer Genome Atlas (TCGA) reveal a positive correlation between Twist1 and mesenchymal genes across cancers, whereas the correlation of TWIST1 with CIN and DSB genes is cancer subtype-specific. Taken together, these studies highlight the mechanistic involvement of Twist1 in the deregulation of factors that maintain genome stability during EMT in colorectal cancer cells. Twist1 overexpression enhances genome instability in the context of EMT that further contributes to cellular heterogeneity. In addition, these studies imply that Twist1 downmodulates nuclear lamins that further alter spatiotemporal organization of the cancer genome and epigenome. Notwithstanding their genetic background, colorectal cancer cells nevertheless maintain their overall ploidy, while the downstream effects of Twist1 enhance CIN and DNA damage enriching for sub-populations of aggressive cancer cells.

Analyses of non-coding somatic drivers in 2,658 cancer whole genomes.

The discovery of drivers of cancer has traditionally focused on protein-coding genes1-4. Here we present analyses of driver point mutations and structural variants in non-coding regions across 2,658 genomes from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium5 of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). For point mutations, we developed a statistically rigorous strategy for combining significance levels from multiple methods of driver discovery that overcomes the limitations of individual methods. For structural variants, we present two methods of driver discovery, and identify regions that are significantly affected by recurrent breakpoints and recurrent somatic juxtapositions. Our analyses confirm previously reported drivers6,7, raise doubts about others and identify novel candidates, including point mutations in the 5' region of TP53, in the 3' untranslated regions of NFKBIZ and TOB1, focal deletions in BRD4 and rearrangements in the loci of AKR1C genes. We show that although point mutations and structural variants that drive cancer are less frequent in non-coding genes and regulatory sequences than in protein-coding genes, additional examples of these drivers will be found as more cancer genomes become available.

The repertoire of mutational signatures in human cancer.

Somatic mutations in cancer genomes are caused by multiple mutational processes, each of which generates a characteristic mutational signature1. Here, as part of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium2 of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA), we characterized mutational signatures using 84,729,690 somatic mutations from 4,645 whole-genome and 19,184 exome sequences that encompass most types of cancer. We identified 49 single-base-substitution, 11 doublet-base-substitution, 4 clustered-base-substitution and 17 small insertion-and-deletion signatures. The substantial size of our dataset, compared with previous analyses3-15, enabled the discovery of new signatures, the separation of overlapping signatures and the decomposition of signatures into components that may represent associated-but distinct-DNA damage, repair and/or replication mechanisms. By estimating the contribution of each signature to the mutational catalogues of individual cancer genomes, we revealed associations of signatures to exogenous or endogenous exposures, as well as to defective DNA-maintenance processes. However, many signatures are of unknown cause. This analysis provides a systematic perspective on the repertoire of mutational processes that contribute to the development of human cancer.

Local Determinants of the Mutational Landscape of the Human Genome.

Large-scale chromatin features, such as replication time and accessibility influence the rate of somatic and germline mutations at the megabase scale. This article reviews how local chromatin structures -e.g., DNA wrapped around nucleosomes, transcription factors bound to DNA- affect the mutation rate at a local scale. It dissects how the interaction of some mutagenic agents and/or DNA repair systems with these local structures influence the generation of mutations. We discuss how this local mutation rate variability affects our understanding of the evolution of the genomic sequence, and the study of the evolution of organisms and tumors.

Somatic and Germline Mutation Periodicity Follow the Orientation of the DNA Minor Groove around Nucleosomes.

Mutation rates along the genome are highly variable and influenced by several chromatin features. Here, we addressed how nucleosomes, the most pervasive chromatin structure in eukaryotes, affect the generation of mutations. We discovered that within nucleosomes, the somatic mutation rate across several tumor cohorts exhibits a strong 10 base pair (bp) periodicity. This periodic pattern tracks the alternation of the DNA minor groove facing toward and away from the histones. The strength and phase of the mutation rate periodicity are determined by the mutational processes active in tumors. We uncovered similar periodic patterns in the genetic variation among human and Arabidopsis populations, also detectable in their divergence from close species, indicating that the same principles underlie germline and somatic mutation rates. We propose that differential DNA damage and repair processes dependent on the minor groove orientation in nucleosome-bound DNA contribute to the 10-bp periodicity in AT/CG content in eukaryotic genomes.

A Pan-cancer Landscape of Interactions between Solid Tumors and Infiltrating Immune Cell Populations.

Purpose: Throughout their development, tumors are challenged by the immune system, and they acquire features to evade its surveillance. A systematic view of these traits, which shed light on how tumors respond to immunotherapies, is still lacking.Experimental Design: Here, we computed the relative abundance of an array of immune cell populations to measure the immune infiltration pattern of 9,174 tumors of 29 solid cancers. We then clustered tumors with similar infiltration pattern to define immunophenotypes. Finally, we identified genomic and transcriptomic traits associated to these immunophenotypes across cancer types.Results: In highly cytotoxic immunophenotypes, we found tumors with low clonal heterogeneity enriched for alterations of genes involved in epigenetic regulation, ubiquitin-mediated proteolysis, antigen presentation, and cell-cell communication, which may drive resistance in combination with the ectopic expression of negative immune checkpoints. Tumors with immunophenotypes of intermediate cytotoxicity are characterized by an upregulation of processes involved in neighboring tissue invasion and remodeling that may foster the recruitment of immunosuppressive cells. Tumors with poorly cytotoxic immunophenotype tend to be of more advanced stages and bear a greater burden of copy number alterations and frequent alterations of cell cycle, hedgehog, ß-catenin, and TGFß pathways, which may cause immune depletion.Conclusions: We provide a comprehensive landscape of the characteristics of solid tumors that may influence (or be influenced by) the characteristics of their immune infiltrate. These results may help interpret the response of solid tumors to immunotherapies and guide the development of novel drug combination strategies. Clin Cancer Res; 24(15); 3717-28. ©2018 AACR.

Reduced mutation rate in exons due to differential mismatch repair.

While recent studies have identified higher than anticipated heterogeneity of mutation rate across genomic regions, mutations in exons and introns are assumed to be generated at the same rate. Here we find fewer somatic mutations in exons than expected from their sequence content and demonstrate that this is not due to purifying selection. Instead, we show that it is caused by higher mismatch-repair activity in exonic than in intronic regions. Our findings have important implications for understanding of mutational and DNA repair processes and knowledge of the evolution of eukaryotic genes, and they have practical ramifications for the study of evolution of both tumors and species.

Whole-genome landscapes of major melanoma subtypes.

Melanoma of the skin is a common cancer only in Europeans, whereas it arises in internal body surfaces (mucosal sites) and on the hands and feet (acral sites) in people throughout the world. Here we report analysis of whole-genome sequences from cutaneous, acral and mucosal subtypes of melanoma. The heavily mutated landscape of coding and non-coding mutations in cutaneous melanoma resolved novel signatures of mutagenesis attributable to ultraviolet radiation. However, acral and mucosal melanomas were dominated by structural changes and mutation signatures of unknown aetiology, not previously identified in melanoma. The number of genes affected by recurrent mutations disrupting non-coding sequences was similar to that affected by recurrent mutations to coding sequences. Significantly mutated genes included BRAF, CDKN2A, NRAS and TP53 in cutaneous melanoma, BRAF, NRAS and NF1 in acral melanoma and SF3B1 in mucosal melanoma. Mutations affecting the TERT promoter were the most frequent of all; however, neither they nor ATRX mutations, which correlate with alternative telomere lengthening, were associated with greater telomere length. Most melanomas had potentially actionable mutations, most in components of the mitogen-activated protein kinase and phosphoinositol kinase pathways. The whole-genome mutation landscape of melanoma reveals diverse carcinogenic processes across its subtypes, some unrelated to sun exposure, and extends potential involvement of the non-coding genome in its pathogenesis.

OncodriveFML: a general framework to identify coding and non-coding regions with cancer driver mutations.

Distinguishing the driver mutations from somatic mutations in a tumor genome is one of the major challenges of cancer research. This challenge is more acute and far from solved for non-coding mutations. Here we present OncodriveFML, a method designed to analyze the pattern of somatic mutations across tumors in both coding and non-coding genomic regions to identify signals of positive selection, and therefore, their involvement in tumorigenesis. We describe the method and illustrate its usefulness to identify protein-coding genes, promoters, untranslated regions, intronic splice regions, and lncRNAs-containing driver mutations in several malignancies.