A Cluster of Evolutionarily Recent KRAB Zinc Finger Proteins Protects Cancer Cells from Replicative Stress-Induced Inflammation

Heterochromatin loss and genetic instability enhance cancer progression by favoring clonal diversity, yet uncontrolled replicative stress leads to mitotic catastrophe and inflammatory responses that promote immune rejection. KRAB domain-containing zinc finger proteins (KZFP) contribute to heterochromatin maintenance at transposable elements (TE). Here, we identified an association of upregulation of a cluster of primate-specific KZFPs with poor prognosis, increased copy-number alterations, and changes in the tumor microenvironment in diffuse large B-cell lymphoma (DLBCL). Depleting two of these KZFPs targeting evolutionarily recent TEs, ZNF587 and ZNF417, impaired the proliferation of cells derived from DLBCL and several other tumor types. ZNF587 and ZNF417 depletion led to heterochromatin redistribution, replicative stress, and cGAS-STING-mediated induction of an interferon/inflammatory response, which enhanced susceptibility to macrophage-mediated phagocytosis and increased surface expression of HLA-I, together with presentation of a neoimmunopeptidome. Thus, cancer cells can exploit KZFPs to dampen TE-originating surveillance mechanisms, which likely facilitates clonal expansion, diversification, and immune evasion.

SIGNIFICANCE

Upregulation of a cluster of primate-specific KRAB zinc finger proteins in cancer cells prevents replicative stress and inflammation by regulating heterochromatin maintenance, which could facilitate the development of improved biomarkers and treatments.

SARS-CoV-2 hijacks a cell damage response, which induces transcription of a more efficient Spike S-acyltransferase

SARS-CoV-2 infection requires Spike protein-mediated fusion between the viral and cellular membranes. The fusogenic activity of Spike depends on its post-translational lipid modification by host S-acyltransferases, predominantly ZDHHC20. Previous observations indicate that SARS-CoV-2 infection augments the S-acylation of Spike when compared to mere Spike transfection. Here, we find that SARS-CoV-2 infection triggers a change in the transcriptional start site of the zdhhc20 gene, both in cells and in an in vivo infection model, resulting in a 67-amino-acid-long N-terminally extended protein with approx. 40 times higher Spike acylating activity, resulting in enhanced fusion of viruses with host cells. Furthermore, we observed the same induced transcriptional change in response to other challenges, such as chemically induced colitis and pore-forming toxins, indicating that SARS-CoV-2 hijacks an existing cell damage response pathway to optimize it fusion glycoprotein.

LnCompare: gene set feature analysis for human long non-coding RNAs

Interest in the biological roles of long noncoding RNAs (lncRNAs) has resulted in growing numbers of studies that produce large sets of candidate genes, for example, differentially expressed between two conditions. For sets of protein-coding genes, ontology and pathway analyses are powerful tools for generating new insights from statistical enrichment of gene features. Here we present the LnCompare web server, an equivalent resource for studying the properties of lncRNA gene sets. The Gene Set Feature Comparison mode tests for enrichment amongst a panel of quantitative and categorical features, spanning gene structure, evolutionary conservation, expression, subcellular localization, repetitive sequences and disease association. Moreover, in Similar Gene Identification mode, users may identify other lncRNAs by similarity across a defined range of features. Comprehensive results may be downloaded in tabular and graphical formats, in addition to the entire feature resource. LnCompare will empower researchers to extract useful hypotheses and candidates from lncRNA gene sets.

Pan-cancer analysis of whole genomes

Cancer is driven by genetic change, and the advent of massively parallel sequencing has enabled systematic documentation of this variation at the whole-genome scale1-3. Here we report the integrative analysis of 2,658 whole-cancer genomes and their matching normal tissues across 38 tumour types from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). We describe the generation of the PCAWG resource, facilitated by international data sharing using compute clouds. On average, cancer genomes contained 4-5 driver mutations when combining coding and non-coding genomic elements; however, in around 5% of cases no drivers were identified, suggesting that cancer driver discovery is not yet complete. Chromothripsis, in which many clustered structural variants arise in a single catastrophic event, is frequently an early event in tumour evolution; in acral melanoma, for example, these events precede most somatic point mutations and affect several cancer-associated genes simultaneously. Cancers with abnormal telomere maintenance often originate from tissues with low replicative activity and show several mechanisms of preventing telomere attrition to critical levels. Common and rare germline variants affect patterns of somatic mutation, including point mutations, structural variants and somatic retrotransposition. A collection of papers from the PCAWG Consortium describes non-coding mutations that drive cancer beyond those in the TERT promoter4; identifies new signatures of mutational processes that cause base substitutions, small insertions and deletions and structural variation5,6; analyses timings and patterns of tumour evolution7; describes the diverse transcriptional consequences of somatic mutation on splicing, expression levels, fusion genes and promoter activity8,9; and evaluates a range of more-specialized features of cancer genomes8,10-18.

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.

Global Positioning System: Understanding Long Noncoding RNAs through Subcellular Localization

The localization of long noncoding RNAs (lncRNAs) within the cell is the primary determinant of their molecular functions. LncRNAs are often thought of as chromatin-restricted regulators of gene transcription and chromatin structure. However, a rich population of cytoplasmic lncRNAs has come to light, with diverse roles including translational regulation, signaling, and respiration. RNA maps of increasing resolution and scope are revealing a subcellular world of highly specific localization patterns and hint at sequence-based address codes specifying lncRNA fates. We propose a new framework for analyzing sequencing-based data, which suggests that numbers of cytoplasmic lncRNA molecules rival those in the nucleus. New techniques promise to create high-resolution, transcriptome-wide maps associated with all organelles of the mammalian cell. Given its intimate link to molecular roles, subcellular localization provides a means of unlocking the mystery of lncRNA functions.

Ancient exapted transposable elements promote nuclear enrichment of human long noncoding RNAs

The sequence domains underlying long noncoding RNA (lncRNA) activities, including their characteristic nuclear enrichment, remain largely unknown. It has been proposed that these domains can originate from neofunctionalized fragments of transposable elements (TEs), otherwise known as RIDLs (repeat insertion domains of lncRNA), although just a handful have been identified. It is challenging to distinguish functional RIDL instances against a numerous genomic background of neutrally evolving TEs. We here show evidence that a subset of TE types experience evolutionary selection in the context of lncRNA exons. Together these comprise an enrichment group of 5374 TE fragments in 3566 loci. Their host lncRNAs tend to be functionally validated and associated with disease. This RIDL group was used to explore the relationship between TEs and lncRNA subcellular localization. By using global localization data from 10 human cell lines, we uncover a dose-dependent relationship between nuclear/cytoplasmic distribution and evolutionarily conserved L2b, MIRb, and MIRc elements. This is observed in multiple cell types and is unaffected by confounders of transcript length or expression. Experimental validation with engineered transgenes shows that these TEs drive nuclear enrichment in a natural sequence context. Together these data reveal a role for TEs in regulating the subcellular localization of lncRNAs.

Ancient exapted transposable elements promote nuclear enrichment of human long noncoding RNAs

The sequence domains underlying long noncoding RNA (lncRNA) activities, including their characteristic nuclear enrichment, remain largely unknown. It has been proposed that these domains can originate from neofunctionalized fragments of transposable elements (TEs), otherwise known as RIDLs (repeat insertion domains of lncRNA), although just a handful have been identified. It is challenging to distinguish functional RIDL instances against a numerous genomic background of neutrally evolving TEs. We here show evidence that a subset of TE types experience evolutionary selection in the context of lncRNA exons. Together these comprise an enrichment group of 5374 TE fragments in 3566 loci. Their host lncRNAs tend to be functionally validated and associated with disease. This RIDL group was used to explore the relationship between TEs and lncRNA subcellular localization. By using global localization data from 10 human cell lines, we uncover a dose-dependent relationship between nuclear/cytoplasmic distribution and evolutionarily conserved L2b, MIRb, and MIRc elements. This is observed in multiple cell types and is unaffected by confounders of transcript length or expression. Experimental validation with engineered transgenes shows that these TEs drive nuclear enrichment in a natural sequence context. Together these data reveal a role for TEs in regulating the subcellular localization of lncRNAs.

LncATLAS database for subcellular localization of long noncoding RNAs

The subcellular localization of long noncoding RNAs (lncRNAs) holds valuable clues to their molecular function. However, measuring localization of newly discovered lncRNAs involves time-consuming and costly experimental methods. We have created "lncATLAS," a comprehensive resource of lncRNA localization in human cells based on RNA-sequencing data sets. Altogether, 6768 GENCODE-annotated lncRNAs are represented across various compartments of 15 cell lines. We introduce relative concentration index (RCI) as a useful measure of localization derived from ensemble RNA-seq measurements. LncATLAS is accessible through an intuitive and informative webserver, from which lncRNAs of interest are accessed using identifiers or names. Localization is presented across cell types and organelles, and may be compared to the distribution of all other genes. Publication-quality figures and raw data tables are automatically generated with each query, and the entire data set is also available to download. LncATLAS makes lncRNA subcellular localization data available to the widest possible number of researchers. It is available at lncatlas.crg.eu.

Cytoplasmic long noncoding RNAs are frequently bound to and degraded at ribosomes in human cells

Recent footprinting studies have made the surprising observation that long noncoding RNAs (lncRNAs) physically interact with ribosomes. However, these findings remain controversial, and the overall proportion of cytoplasmic lncRNAs involved is unknown. Here we make a global, absolute estimate of the cytoplasmic and ribosome-associated population of stringently filtered lncRNAs in a human cell line using polysome profiling coupled to spike-in normalized microarray analysis. Fifty-four percent of expressed lncRNAs are detected in the cytoplasm. The majority of these (70%) have >50% of their cytoplasmic copies associated with polysomal fractions. These interactions are lost upon disruption of ribosomes by puromycin. Polysomal lncRNAs are distinguished by a number of 5' mRNA-like features, including capping and 5'UTR length. On the other hand, nonpolysomal "free cytoplasmic" lncRNAs have more conserved promoters and a wider range of expression across cell types. Exons of polysomal lncRNAs are depleted of endogenous retroviral insertions, suggesting a role for repetitive elements in lncRNA localization. Finally, we show that blocking of ribosomal elongation results in stabilization of many associated lncRNAs. Together these findings suggest that the ribosome is the default destination for the majority of cytoplasmic long noncoding RNAs and may play a role in their degradation.