RNA editing, epitranscriptomics, and processing in cancer progression

Ocular A-to-I RNA editing signatures associated with SARS-CoV-2 infection

Ophthalmic manifestations have recently been observed in acute and post-acute complications of COVID-19 caused by SARS-CoV-2 infection. Our precious study has shown that host RNA editing is linked to RNA viral infection, yet ocular adenosine to inosine (A-to-I) RNA editing during SARS-CoV-2 infection remains uninvestigated in COVID-19. Herein we used an epitranscriptomic pipeline to analyze 37 samples and investigate A-to-I editing associated with SARS-CoV-2 infection, in five ocular tissue types including the conjunctiva, limbus, cornea, sclera, and retinal organoids. Our results revealed dramatically altered A-to-I RNA editing across the five ocular tissues. Notably, the transcriptome-wide average level of RNA editing was increased in the cornea but generally decreased in the other four ocular tissues. Functional enrichment analysis showed that differential RNA editing (DRE) was mainly in genes related to ubiquitin-dependent protein catabolic process, transcriptional regulation, and RNA splicing. In addition to tissue-specific RNA editing found in each tissue, common RNA editing was observed across different tissues, especially in the innate antiviral immune gene MAVS and the E3 ubiquitin-protein ligase MDM2. Analysis in retinal organoids further revealed highly dynamic RNA editing alterations over time during SARS-CoV-2 infection. Our study thus suggested the potential role played by RNA editing in ophthalmic manifestations of COVID-19, and highlighted its potential transcriptome impact, especially on innate immunity.

N7-methylguanosine methylation-related regulator genes as biological markers in predicting prognosis for melanoma

The aim of this study is to find those N7-methylguanosine (m7G) methylation-related regulator genes (m7GMRRGs) which were associated with melanoma prognosis and use them to develop a prognostic prediction model. Clinical information was retrieved online from The Cancer Gene Atlas (TCGA) and the Gene Expression Omnibus (GEO). R software was used to extract m⁷GMRRGs by differential expression analysis. To create a prognostic risk model, univariate and multivariate Cox regression analyses were employed for the evaluation of the prognostic significance of m⁷G methylation modifiers. Internal validation using cohort from TCGA (training set) and external validation using cohort from GEO (validation set) of the model were carried out. The model's predictive performance was confirmed by using the Kaplan-Meier, univariate, and multivariate Cox regression, and receiver operating characteristic curve (ROC) by constructing column line plots incorporating clinical factor characteristics. Immune infiltration analyses were performed to assess the immune function of m⁷GMRRGs. Drug sensitivity analysis was conducted to study chemotherapeutic drug treatment cues. Prognostic models using four m⁷GMRRGs (EIF4E3, LARP1, NCBP3, and IFIT5) showed good prognostic power in training and validation sets. The area under the curve (AUC) at 1, 3, and 5 years for GEO-melanoma were 0.689, 0.704, and 0.726, respectively. The prediction model could distinctly classify patients with melanoma into different risk subgroups (P < 0.001 for TCGA-melanoma and P < 0.05 for GEO-melanoma). Clinical characteristics were taken into account in Cox regression and AUC analysis, which highlighted that the risk score served as an independent risk factor determining the prognosis of patients with melanoma. Immuno-infiltration analysis showed that m⁷GMRRGs could potentially regulate CD8⁺ T cells as well as regulatory T cells (Treg cells). Results of our study indicate a association between m⁷GMRRGs and melanoma prognosis, and the prognostic prediction model using m⁷GMRRGs may predict the prognosis of patients with melanoma well. Nevertheless, these results may provide a clue for potential better options of melanoma treatment but need further validation in futural studies.

Identification and Validation of a Prognostic Prediction Model of m6A Regulator-Related LncRNAs in Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) is a highly invasive malignancy prone to recurrence, and patients with HCC have a low 5-year survival rate. Long non-coding RNAs (lncRNAs) play a vital role in the occurrence and development of HCC. N6-methyladenosine methylation (m6A) is the most common modification influencing cancer development. Here, we used the transcriptome of m6A regulators and lncRNAs, along with the complete corresponding clinical HCC patient information obtained from The Cancer Genome Atlas (TCGA), to explore the role of m6A regulator-related lncRNA (m6ARlnc) as a prognostic biomarker in patients with HCC. The prognostic m6ARlnc was selected using Pearson correlation and univariate Cox regression analyses. Moreover, three clusters were obtained via consensus clustering analysis and further investigated for differences in immune infiltration, immune microenvironment, and prognosis. Subsequently, nine m6ARlncs were identified with Lasso-Cox regression analysis to construct the prognostic signature m6A-9LPS for patients with HCC in the training cohort (n = 226). Based on m6A-9LPS, the risk score for each case was calculated. Patients were then divided into high- and low-risk subgroups based on the cutoff value set by the X-tile software. m6A-9LPS showed a strong prognosis prediction ability in the validation cohort (n = 116), the whole cohort (n = 342), and even clinicopathological stratified survival analysis. Combining the risk score and clinical characteristics, we established a nomogram for predicting the overall survival (OS) of patients. To further understand the mechanism underlying the m6A-9LPS-based classification of prognosis differences, KEGG and GO enrichment analyses, competitive endogenous RNA (ceRNA) network, chemotherapeutic agent sensibility, and immune checkpoint expression level were assessed. Taken together, m6A-9LPS could be used as a precise prediction model for the prognosis of patients with HCC, which will help in individualized treatment of HCC.

Finding the hot spot: identifying immune sensitive gastrointestinal tumors

Although researchers have been trying to harness the immune system for over 100 years, the advent of immune checkpoint blockers (ICB) marks an era of significant clinical outcomes in various metastatic solid tumors, characterized by complete and durable responses. ICBs are monoclonal antibodies that target either of a pair of transmembrane molecules in tumors or T-cells involved in immune evasion. Currently 2 ICBs targeting the checkpoint program death 1 (PD-1), nivolumab and pembrolizumab, and one cytotoxic lymphocyte antigen-4 (CTLA-4) inhibitor (ipilimumab) are approved in gastrointestinal malignancies. We review herein the current evidence on predictive biomarkers for ICB response in gastrointestinal tumors. A review of literature based on the National Cancer Institute list of FDA-approved drugs for neoplasms and FDA-approved therapies at the FDA website was performed. An initial literature review was based on the American Association for Clinical Research meeting 2019, the American Society of Clinical Oncology meeting 2019 and the European Society of Medical Oncology 2019 proceedings. A systematic search of PubMed was performed involving MeSH browser terms such as biomarkers, immunotherapy, gastrointestinal diseases and neoplasms. When appropriate, American and British terms were used in the search. The most relevant predictor of response to ICBs is microsatellite instability (MSI) and the data is strongest for colorectal cancer. At least 3 prospective trials show evidence of PD-L1 as a predictive biomarker for ICB response in gastroesophageal malignancies. At least one prospective trial has described tumor mutational burden high (TMB-H), independent of MSI, as predictive of response in anal and biliary tract carcinomas. DNA Polymerase Epsilon (POLE) or delta (POL-D) mutations have been implicated in a subset of MSS colorectal cancer with TMB-H but this biomarker requires prospective validation. There is evolving data based on retrospective observations that gene alterations predicting acquired resistance and hyper-progression. Ongoing clinical research is assessing the role of the human microbiome and RNA-editing complex mutations as predictive biomarkers of response to ICBs. MSI has the strongest predictive power among current biomarkers for ICB response in gastrointestinal cancers. Data continue to accumulate from ongoing clinical trials and new biomarkers are emerging from pre-clinical studies, suggesting that drug combinations targeting pathways complimentary to the PD-1/PD-L1 axis inhibition will define a robust field of clinical research.

Host-Directed Antiviral Therapy

Antiviral drugs have traditionally been developed by directly targeting essential viral components. However, this strategy often fails due to the rapid generation of drug-resistant viruses. Recent genome-wide approaches, such as those employing small interfering RNA (siRNA) or clustered regularly interspaced short palindromic repeats (CRISPR) or those using small molecule chemical inhibitors targeting the cellular "kinome," have been used successfully to identify cellular factors that can support virus replication. Since some of these cellular factors are critical for virus replication, but are dispensable for the host, they can serve as novel targets for antiviral drug development. In addition, potentiation of immune responses, regulation of cytokine storms, and modulation of epigenetic changes upon virus infections are also feasible approaches to control infections. Because it is less likely that viruses will mutate to replace missing cellular functions, the chance of generating drug-resistant mutants with host-targeted inhibitor approaches is minimized. However, drug resistance against some host-directed agents can, in fact, occur under certain circumstances, such as long-term selection pressure of a host-directed antiviral agent that can allow the virus the opportunity to adapt to use an alternate host factor or to alter its affinity toward the target that confers resistance. This review describes novel approaches for antiviral drug development with a focus on host-directed therapies and the potential mechanisms that may account for the acquisition of antiviral drug resistance against host-directed agents.

Significant epitranscriptomes in heterogeneous cancer

Precision medicine places significant emphasis on techniques for the identification of DNA mutations and gene expression by deep sequencing of gene panels to obtain medical data. However, other diverse information that is not easily readable using bioinformatics, including RNA modifications, has emerged as a novel diagnostic and innovative therapy owing to its multifunctional aspects. It is suggested that this breakthrough innovation might open new avenues for the elucidation of uncharacterized cancer cellular functions to develop more precise medical applications. The functional characteristics and regulatory mechanisms of RNA modifications, ie, the epitranscriptome (ETR), which reflects RNA metabolism, remains unclear, mainly due to detection methods being limited. Recent studies have revealed that N6‐methyl adenosine, the most common modification in mRNA in eukaryotes, is affected in various types of cancer and, in some cases, cancer stem cells, but also affects cellular responses to viral infections. The ETR can control cancer cell fate through mRNA splicing, stability, nuclear export, and translation. Here we report on the recent progress of ETR detection methods, and biological findings regarding the significance of ETR in cancer precision medicine.

Epigenetic and epitranscriptomic changes in colorectal cancer: Diagnostic, prognostic, and treatment implications

A cancer cell is the final product of a complex mixture of genetic, epigenetic and epitranscriptomic alterations, whose final interplay contribute to cancer onset and progression. This is specifically true for colorectal cancer, a tumor with a strong epigenetic component, which acts earlier than any other genetic alteration in promoting cancer cell malignant transformation. The pattern of progressive, and usually subtype-specific, DNA and histone modifications that occur in colorectal cancer has been extensively studied in the last decade, providing plenty of data to explore. For this tumor, it became recently evident that also RNA modifications play a relevant role in the activation of oncogenes or repression of tumor suppressor genes. In this review we provide a brief overview of all epigenetic and epitranscriptomic changes that have been found associated to colorectal cancer till now. We explore the impact of these alterations in cancer prognosis and response to treatment and discuss their potential use as cancer biomarkers.

Global analysis of H3K27me3 as an epigenetic marker in prostate cancer progression

BACKGROUND

H3K27me3 histone marks shape the inhibition of gene transcription. In prostate cancer, the deregulation of H3K27me3 marks might play a role in prostate tumor progression.

METHODS

We investigated genome-wide H3K27me3 histone methylation profile using chromatin immunoprecipitation (ChIP) and 2X400K promoter microarrays to identify differentially-enriched regions in biopsy samples from prostate cancer patients. H3K27me3 marks were assessed in 34 prostate tumors: 11 with Gleason score > 7 (GS > 7), 10 with Gleason score ≤ 7 (GS ≤ 7), and 13 morphologically normal prostate samples.

RESULTS

Here, H3K27me3 profiling identified an average of 386 enriched-genes on promoter regions in healthy control group versus 545 genes in GS ≤ 7 and 748 genes in GS > 7 group. We then ran a factorial discriminant analysis (FDA) and compared the enriched genes in prostate-tumor biopsies and normal biopsies using ANOVA to identify significantly differentially-enriched genes. The analysis identified ALG5, EXOSC8, CBX1, GRID2, GRIN3B, ING3, MYO1D, NPHP3-AS1, MSH6, FBXO11, SND1, SPATS2, TENM4 and TRA2A genes. These genes are possibly associated with prostate cancer. Notably, the H3K27me3 histone mark emerged as a novel regulatory mechanism in poor-prognosis prostate cancer.

CONCLUSIONS

Our findings point to epigenetic mark H3K27me3 as an important event in prostate carcinogenesis and progression. The results reported here provide new molecular insights into the pathogenesis of prostate cancer.

RNA binding to APOBEC deaminases; Not simply a substrate for C to U editing

Apolipoprotein B mRNA Editing Catalytic Polypeptide-like 1 or APOBEC1 was discovered in 1993 as the zinc-dependent cytidine deaminase responsible for the production of an in frame stop codon in apoB mRNA through modification of cytidine at nucleotide position 6666 to uridine. At the time of this discovery there was much speculation concerning the mechanism of base modification RNA editing which has been rekindled by the discovery of multiple C to U RNA editing events in the 3′ UTRs of mRNAs and the finding that other members of the APOBEC family while able to bind RNA, have the biological function of being DNA mutating enzymes. Current research is addressing the mechanism for these nucleotide modification events that appear not to adhere to the mooring sequence-dependent model for APOBEC1 involving the assembly of a multi protein containing editosome. This review will summarize our current understanding of the structure and function of APOBEC proteins and examine how RNA binding to them may be a regulatory mechanism.

Diastereoselectivity of 5-Methyluridine Osmylation Is Inverted inside an RNA Chain

In this study, we investigated the reaction of the osmium tetroxide-bipyridine complex with pyrimidines in RNA. This reagent, which reacts with the diastereotopic 5-6 double bond, thus leading to the formation of two diastereomers, was used in the past to label thymidine and 5-methylcytosine in DNA. In light of the growing interest in post-transcriptional RNA modifications, we addressed the question of whether this reagent could be used for labeling of the naturally occurring RNA modifications 5-methylcytosine and 5-methyluridine. On nucleoside level, 5-methylcytosine and 5-methyluridine revealed a 5- and 12-fold preference, respectively, over their nonmethylated equivalents. Performing the reaction on an RNA level, we could show that the steric environment of a pentanucleotide has a major detrimental impact on the reaction rate of osmylation. Interestingly, this drop in reactivity was due to a dramatic change in diastereoselectivity, which in turn resulted from impediment of the preferred attack via the si side. Thus, while on the nucleoside level, the absolute configuration of the major product of osmylation of 5-methyluridine was (5R,6S)-5-methyluridine glycol-dioxoosmium-bipyridine, reaction with an RNA pentanucleotide afforded the corresponding (5S,6R)-diastereomer as the major product. The change in diastereoselectivity lead to an almost complete loss of selectivity toward 5-methylcytosine in a pentanucleotide context, while 5-methyluridine remained about 8 times more reactive than the canonical pyrimidines. On the basis of these findings, we evaluate the usefulness of osmium tetroxide-bipyridine as a potential label for the 5-methyluridine modification in transcriptome-wide studies.

Principles Governing A-to-I RNA Editing in the Breast Cancer Transcriptome

Little is known about how RNA editing operates in cancer. Transcriptome analysis of 68 normal and cancerous breast tissues revealed that the editing enzyme ADAR acts uniformly, on the same loci, across tissues. In controlled ADAR expression experiments, the editing frequency increased at all loci with ADAR expression levels according to the logistic model. Loci-specific “editabilities,” i.e., propensities to be edited by ADAR, were quantifiable by fitting the logistic function to dose-response data. The editing frequency was increased in tumor cells in comparison to normal controls. Type I interferon response and ADAR DNA copy number together explained 53% of ADAR expression variance in breast cancers. ADAR silencing using small hairpin RNA lentivirus transduction in breast cancer cell lines led to less cell proliferation and more apoptosis. A-to-I editing is a pervasive, yet reproducible, source of variation that is globally controlled by 1q amplification and inflammation, both of which are highly prevalent among human cancers.

•

A-to-I editing is a major source of mRNA variability in breast and other cancers

•

RNA editing is globally controlled by tumor interferon and ADAR copy number

•

Both these factors are highly prevalent among human cancers

•

RNA editing sites might represent a new class of therapeutic targets

Genomic variations of the mevalonate pathway in porokeratosis

Porokeratosis (PK) is a heterogeneous group of keratinization disorders. No causal genes except MVK have been identified, even though the disease was linked to several genomic loci. Here, we performed massively parallel sequencing and exonic CNV screening of 12 isoprenoid genes in 134 index PK patients (61 familial and 73 sporadic) and identified causal mutations in three novel genes (PMVK, MVD, and FDPS) in addition to MVK in the mevalonate pathway. Allelic expression imbalance (AEI) assays were performed in 13 lesional tissues. At least one mutation in one of the four genes in the mevalonate pathway was found in 60 (98%) familial and 53 (73%) sporadic patients, which suggests that isoprenoid biosynthesis via the mevalonate pathway may play a role in the pathogenesis of PK. Significantly reduced expression of the wild allele was common in lesional tissues due to gene conversion or some other unknown mechanism. A G-to-A RNA editing was observed in one lesional tissue without AEI. In addition, we observed correlations between the mutations in the four mevalonate pathway genes and clinical manifestations in the PK patients, which might support a new and simplified classification of PK under the guidance of genetic testing.

DOI:http://dx.doi.org/10.7554/eLife.06322.001

Porokeratosis refers to a group of around twenty skin conditions that involve a build-up of a protein called keratin in skin cells. Keratin forms the tough fibres that give strength to hair and nails, and people suffering from porokeratosis develop hardened skin lesions. Porokeratosis is an uncommon condition; most cases are inherited and have been linked to exposure to ultraviolet light and having a weakened immune system.

Mutations in one gene called MVK are known to cause two forms of the disorder, but it is suspected that other genetic causes of porokeratosis will also be identified. The MVK gene encodes an enzyme that is involved in making chemicals called isoprenoids. This large and diverse class of chemicals provides the building blocks for making many other important molecules in all living species. Zhang, Li et al. have now analysed genetic material from 134 different porokeratosis patients to search for mutations in other genes involved in the production of isoprenoids. The patients examined include 61 people with a family history of the disorder, and 73 cases in which the condition seems to be a one-off occurrence.

This search identified mutations in three additional genes (called PMVK, MVD and FDPS) that are all linked to porokeratosis. Further analysis of these three genes and MVK revealed that about half of the patients with mutations in the MVK gene developed large lesions (that were over 5 centimetres in diameter). However, those with mutations in the other three genes did not develop such large lesions. Mutations in some of the newly identified genes were instead linked to porokeratosis affecting specific areas of the body; for example, PMVK and MVD mutations are linked to porokeratosis localized to the genitals and around the eyes, respectively. This means that, in the future, doctors might be able to simplify the diagnosis of the different varieties of porokeratosis based on information gained via genetic tests.

DOI:http://dx.doi.org/10.7554/eLife.06322.002