Comprehensive characterization of somatic variants associated with intronic polyadenylation in human cancers.

Somatic single nucleotide variants (SNVs) in cancer genome affect gene expression through various mechanisms depending on their genomic location. While somatic SNVs near canonical splice sites have been reported to cause abnormal splicing of cancer-related genes, whether these SNVs can affect gene expression through other mechanisms remains an open question. Here, we analyzed RNA sequencing and exome data from 4,998 cancer patients covering ten cancer types and identified 152 somatic SNVs near splice sites that were associated with abnormal intronic polyadenylation (IPA). IPA-associated somatic variants favored the localization near the donor splice sites compared to the acceptor splice sites. A proportion of SNV-associated IPA events overlapped with premature cleavage and polyadenylation events triggered by U1 small nuclear ribonucleoproteins (snRNP) inhibition. GC content, intron length and polyadenylation signal were three genomic features that differentiated between SNV-associated IPA and intron retention. Notably, IPA-associated SNVs were enriched in tumor suppressor genes (TSGs), including the well-known TSGs such as PTEN and CDH1 with recurrent SNV-associated IPA events. Minigene assay confirmed that SNVs from PTEN, CDH1, VEGFA, GRHL2, CUL3 and WWC2 could lead to IPA. This work reveals that IPA acts as a novel mechanism explaining the functional consequence of somatic SNVs in human cancer.

Pan-Cancer Profiling of Intron Retention and Its Clinical Significance in Diagnosis and Prognosis.

Alternative splicing can produce transcripts that affect cancer development and thus shows potential for cancer diagnosis and treatment. However, intron retention (IR), a type of alternative splicing, has been studied less in cancer biology research. Here, we generated a pan-cancer IR landscape for more than 10,000 samples across 33 cancer types from The Cancer Genome Atlas (TCGA). We characterized differentially retained introns between tumor and normal samples and identified retained introns associated with survival. We discovered 988 differentially retained introns in 14 cancers, some of which demonstrated diagnostic potential in multiple cancer types. We also inferred a large number of prognosis-related introns in 33 cancer types, and the associated genes included well-known cancer hallmarks such as angiogenesis, metastasis, and DNA mutations. Notably, we discovered a novel intron retention inside the 5'UTR of STN1 that is associated with the survival of lung cancer patients. The retained intron reduces translation efficiency by producing upstream open reading frames (uORFs) and thereby inhibits colony formation and cell migration of lung cancer cells. Besides, the IR-based prognostic model achieved good stratification in certain cancers, as illustrated in acute myeloid leukemia. Taken together, we performed a comprehensive IR survey at a pan-cancer level, and the results implied that IR has the potential to be diagnostic and prognostic cancer biomarkers, as well as new drug targets.

The landscape of cryptic antisense transcription in human cancers reveals an oncogenic noncoding RNA in lung cancer.

Cryptic transcription initiation has been previously linked to activation of oncogenic transcripts. However, the prevalence and impact of cryptic antisense transcription from the opposite strand of protein-coding genes were mostly unknown in cancer. Applying a robust computational pipeline to publicly available transcriptome and epigenome datasets, we identified hundreds of previously unannotated cryptic antisense polyadenylated transcripts (CAPTs) that were enriched in tumor samples. We showed that the activation of cryptic antisense transcription was associated with increased chromatin accessibility and active histone marks. Accordingly, we found that many of the antisense transcripts were inducible by treatment of epigenetic drugs. Moreover, CRISPR-mediated epigenetic editing assays revealed that transcription of a noncoding RNA LRRK1-CAPT promoted LUSC cell proliferation, suggesting its oncogenic role. Our findings largely expand our understanding of cancer-associated transcription events, which may facilitate the development of novel strategies for cancer diagnosis and treatment.

Low level of ARID1A contributes to adaptive immune resistance and sensitizes triple-negative breast cancer to immune checkpoint inhibitors.

BACKGROUND: Immune checkpoint inhibitors (ICIs) shed new light on triple-negative breast cancer (TNBC), but only a minority of patients demonstrate response. Therefore, adaptive immune resistance (AIR) needs to be further defined to guide the development of ICI regimens. METHODS: Databases, including The Cancer Genome Atlas, Gene Ontology Resource, University of California Santa Cruz Genome Browser, and Pubmed, were used to screen epigenetic modulators, regulators for CD8+ T cells, and transcriptional regulators of programmed cell death-ligand 1 (PD-L1). Human peripheral blood mononuclear cell (Hu-PBMC) reconstruction mice were adopted for xenograft transplantation. Tumor specimens from a TNBC cohort and the clinical trial CTR20191353 were retrospectively analyzed. RNA-sequencing, Western blotting, qPCR and immunohistochemistry were used to assess gene expression. Coculture assays were performed to evaluate the regulation of TNBC cells on T cells. Chromatin immunoprecipitation and transposase-accessible chromatin sequencing were used to determine chromatin-binding and accessibility. RESULTS: The epigenetic modulator AT-rich interaction domain 1A (ARID1A) gene demonstrated the highest expression association with AIR relative to other epigenetic modulators in TNBC patients. Low ARID1A expression in TNBC, causing an immunosuppressive microenvironment, promoted AIR and inhibited CD8+ T cell infiltration and activity through upregulating PD-L1. However, ARID1A did not directly regulate PD-L1 expression. We found that ARID1A directly bound the promoter of nucleophosmin 1 (NPM1) and that low ARID1A expression increased NPM1 chromatin accessibility as well as gene expression, further activating PD-L1 transcription. In Hu-PBMC mice, atezolizumab demonstrated the potential to reverse ARID1A deficiency-induced AIR in TNBC by reducing tumor malignancy and activating anti-tumor immunity. In CTR20191353, ARID1A-low patients derived more benefit from pucotenlimab compared to ARID1A-high patients. CONCLUSIONS: In AIR epigenetics, low ARID1A expression in TNBC contributed to AIR via the ARID1A/NPM1/PD-L1 axis, leading to poor outcome but sensitivity to ICI treatment.

Lung cancer cells expressing a shortened CDK16 3'UTR escape senescence through impaired miR-485-5p targeting.

Inducing senescence in cancer cells is an emerging strategy for cancer therapy. The dysregulation and mutation of genes encoding cyclin-dependent kinases (CDKs) have been implicated in various human cancers. However, whether CDK can induce cancer cell senescence remains poorly understood. We observed that CDK16 expression was high in multiple cancer types, including lung cancer, whereas various replicative senescence models displayed low CDK16 expression. CDK16 knockdown caused senescence-associated phenotypes in lung cancer cell lines. Interestingly, the CDK16 3' UTR was shortened in cancer and lengthened in senescence models, which was regulated by alternative polyadenylation (APA). The longer 3'UTR [using the distal polyA (pA) site] generated less protein than the shorter one (using the proximal pA site). Since microRNAs (miRNAs) usually bind to the 3'UTR of target genes to suppress their expression, we investigated whether miRNAs targeting the region between the shortened and longer 3'UTR are responsible for the reduced expression. We found that miR-485-5p targeted the 3'UTR between the distal and proximal pA site and caused senescence-associated phenotypes by reducing protein production from the longer CDK16 transcript. Of note, CDK16 knockdown led to a reduced expression of MYC proto-oncogene, bHLH transcription factor (MYC) and CD274 molecule (PD-L1), which in turn enhanced the tumor-suppressive effects of senescent cancer cells. The present study discovered that CDK16, whose expression is under the regulation of APA and miR-485-5p, is a potential target for prosenescence therapy for lung cancer.

Integrative Analysis of Somatic Mutations in Non-coding Regions Altering RNA Secondary Structures in Cancer Genomes.

RNA secondary structure may influence many cellular processes, including RNA processing, stability, localization, and translation. Single-nucleotide variations (SNVs) that alter RNA secondary structure, referred to as riboSNitches, are potentially causative of human diseases, especially in untranslated regions (UTRs) and noncoding RNAs (ncRNAs). The functions of somatic mutations that act as riboSNitches in cancer development remain poorly understood. In this study, we developed a computational pipeline called SNIPER (riboSNitch-enriched or depleted elements in cancer genomes), which employs MeanDiff and EucDiff to detect riboSNitches and then identifies riboSNitch-enriched or riboSNitch-depleted non-coding elements across tumors. SNIPER is available at github: https://github.com/suzhixi/SNIPER/ . We found that riboSNitches were more likely to be pathogenic. Moreover, we predicted several UTRs and lncRNAs (long non-coding RNA) that significantly enriched or depleted riboSNitches in cancer genomes, indicative of potential cancer driver or essential noncoding elements. Our study highlights the possibly neglected importance of RNA secondary structure in cancer genomes and provides a new strategy to identify new cancer-associated genes.

Structural genomic changes underlie alternative reproductive strategies in the ruff (Philomachus pugnax).

The ruff is a Palearctic wader with a spectacular lekking behavior where highly ornamented males compete for females. This bird has one of the most remarkable mating systems in the animal kingdom, comprising three different male morphs (independents, satellites and faeders) that differ in behavior, plumage color and body size. Remarkably, the satellite and faeder morphs are controlled by dominant alleles. Here we have used whole-genome sequencing and resolved the enigma of how such complex phenotypic differences can have a simple genetic basis. The Satellite and Faeder alleles are both associated with a 4.5-Mb inversion that occurred about 3.8 million years ago. We propose an evolutionary scenario where the Satellite chromosome arose by a rare recombination event about 500,000 years ago. The ruff mating system is the result of an evolutionary process in which multiple genetic changes contributing to phenotypic differences between morphs have accumulated within the inverted region.

Genome sequence of cultivated Upland cotton (Gossypium hirsutum TM-1) provides insights into genome evolution.

Gossypium hirsutum has proven difficult to sequence owing to its complex allotetraploid (AtDt) genome. Here we produce a draft genome using 181-fold paired-end sequences assisted by fivefold BAC-to-BAC sequences and a high-resolution genetic map. In our assembly 88.5% of the 2,173-Mb scaffolds, which cover 89.6%∼96.7% of the AtDt genome, are anchored and oriented to 26 pseudochromosomes. Comparison of this G. hirsutum AtDt genome with the already sequenced diploid Gossypium arboreum (AA) and Gossypium raimondii (DD) genomes revealed conserved gene order. Repeated sequences account for 67.2% of the AtDt genome, and transposable elements (TEs) originating from Dt seem more active than from At. Reduction in the AtDt genome size occurred after allopolyploidization. The A or At genome may have undergone positive selection for fiber traits. Concerted evolution of different regulatory mechanisms for Cellulose synthase (CesA) and 1-Aminocyclopropane-1-carboxylic acid oxidase1 and 3 (ACO1,3) may be important for enhanced fiber production in G. hirsutum.