CIGAR‐seq,a CRISPR/Cas‐based method for unbiased screening of novel mRNA modification regulators

Cellular RNA is decorated with over 170 types of chemical modifications. Many modifications in mRNA, including m⁶A and m⁵C, have been associated with critical cellular functions under physiological and/or pathological conditions. To understand the biological functions of these modifications, it is vital to identify the regulators that modulate the modification rate. However, a high‐throughput method for unbiased screening of these regulators is so far lacking. Here, we report such a method combining pooled CRISPR screen and reporters with RNA modification readout, termed CRISPR integrated gRNA and reporter sequencing (CIGAR‐seq). Using CIGAR‐seq, we discovered NSUN6 as a novel mRNA m⁵C methyltransferase. Subsequent mRNA bisulfite sequencing in HAP1 cells without or with NSUN6 and/or NSUN2 knockout showed that NSUN6 and NSUN2 worked on non‐overlapping subsets of mRNA m⁵C sites and together contributed to almost all the m⁵C modification in mRNA. Finally, using m¹A as an example, we demonstrated that CIGAR‐seq can be easily adapted for identifying regulators of other mRNA modification.     

NSUN4 Is a Dual Function Mitochondrial Protein Required for Both Methylation of 12S rRNA and Coordination of Mitoribosomal Assembly

Biogenesis of mammalian mitochondrial ribosomes requires a concerted maturation of both the small (SSU) and large subunit (LSU). We demonstrate here that the m⁵C methyltransferase NSUN4, which forms a complex with MTERF4, is essential in mitochondrial ribosomal biogenesis as mitochondrial translation is abolished in conditional Nsun4 mouse knockouts. Deep sequencing of bisulfite-treated RNA shows that NSUN4 methylates cytosine 911 in 12S rRNA (m5C911) of the SSU. Surprisingly, NSUN4 does not need MTERF4 to generate this modification. Instead, the NSUN4/MTERF4 complex is required to assemble the SSU and LSU to form a monosome. NSUN4 is thus a dual function protein, which on the one hand is needed for 12S rRNA methylation and, on the other hand interacts with MTERF4 to facilitate monosome assembly. The presented data suggest that NSUN4 has a key role in controlling a final step in ribosome biogenesis to ensure that only the mature SSU and LSU are assembled.     

NSUN6 is a human RNA methyltransferase that catalyzes formation of m5C72 in specific tRNAs

Many cellular RNAs require modification of specific residues for their biogenesis, structure, and function. 5-methylcytosine (m⁵C) is a common chemical modification in DNA and RNA but in contrast to the DNA modifying enzymes, only little is known about the methyltransferases that establish m⁵C modifications in RNA. The putative RNA methyltransferase NSUN6 belongs to the family of Nol1/Nop2/SUN domain (NSUN) proteins, but so far its cellular function has remained unknown. To reveal the target spectrum of human NSUN6, we applied UV crosslinking and analysis of cDNA (CRAC) as well as chemical crosslinking with 5-azacytidine. We found that human NSUN6 is associated with tRNAs and acts as a tRNA methyltransferase. Furthermore, we uncovered tRNA^Cys and tRNA^Thr as RNA substrates of NSUN6 and identified the cytosine C72 at the 3′ end of the tRNA acceptor stem as the target nucleoside. Interestingly, target recognition in vitro depends on the presence of the 3′-CCA tail. Together with the finding that NSUN6 localizes to the cytoplasm and largely colocalizes with marker proteins for the Golgi apparatus and pericentriolar matrix, our data suggest that NSUN6 modifies tRNAs in a late step in their biogenesis.     

NSUN3 and ABH1 modify the wobble position of mt‐tRNAMet to expand codon recognition in mitochondrial translation

Mitochondrial gene expression uses a non‐universal genetic code in mammals. Besides reading the conventional AUG codon, mitochondrial (mt‐)tRNA^Met mediates incorporation of methionine on AUA and AUU codons during translation initiation and on AUA codons during elongation. We show that the RNA methyltransferase NSUN3 localises to mitochondria and interacts with mt‐tRNA^Met to methylate cytosine 34 (C34) at the wobble position. NSUN3 specifically recognises the anticodon stem loop (ASL) of the tRNA, explaining why a mutation that compromises ASL basepairing leads to disease. We further identify ALKBH1/ABH1 as the dioxygenase responsible for oxidising m⁵C34 of mt‐tRNA^Met to generate an f⁵C34 modification. In vitro codon recognition studies with mitochondrial translation factors reveal preferential utilisation of m⁵C34 mt‐tRNA^Met in initiation. Depletion of either NSUN3 or ABH1 strongly affects mitochondrial translation in human cells, implying that modifications generated by both enzymes are necessary for mt‐tRNA^Met function. Together, our data reveal how modifications in mt‐tRNA^Met are generated by the sequential action of NSUN3 and ABH1, allowing the single mitochondrial tRNA^Met to recognise the different codons encoding methionine.     

A new m6A methylation-related gene signature for prognostic value in patient with urothelial carcinoma of the bladder

Background: Bladder cancer (BC) is one of the most common malignant urological cancer in the world. Because of its characteristic of easy-recurrence and muscle-invasive, advances in our genetic understanding of bladder cancer should be translated into prognostic indicators.Methods: We investigated 16 m⁶A RNA methylation regulators from The Cancer Genome Atlas (TCGA) database and The Human Protein Atlas (HPA) database. The expression profile, clinical application as well as prognostic value of these genes in UC were investigated. Moreover, we further explored the correlation between RNA methylation genes and biological functions, pathways and immune status.Results: Five m⁶A-related genes (HNRNPC, YTHDF2, YTHDF1, HNRNPA2B1, METTL3) up-regulated in UC tissues, while three regulators (ZC3H13, METTL16, FTO) down-regulated in UC. FTO and YTHDF2 show biomarker potential for the prognosis of UC patients. In addition, these identified genes may related with essential functions and core molecular pathways.Conclusions: Our research shows that two m⁶A RNA methylation regulators can serve as reliable prognostic biomarkers of UC, which might be exerted as potential targets of therapeutic strategies.     

HNRNPC可能成为肺腺癌独立的预后因子

为探讨RNA m⁶A甲基化调节因子在肺腺癌中的作用,从TCGA数据库下载肺腺癌患者的RNA表达数据和临床数据。通过limma软件包分析12种m⁶A调节剂的表达情况。使用Pheatmap、vioplot和corrplot软件包生成热图、小提琴图和表达相关图。采用Kaplan-Meier方法分别计算肺腺癌中12种RNA m⁶A调节因子的生存曲线。使用Cox回归和Kaplan-Meier方法分析TCGA肺腺癌患者的总体存活相关的临床病理学特征。最后用Kruskal(KS)检验和logistic回归分析临床病理学特征与HNRNPC表达的关系。 在肺腺癌的TCGA队列中,发现HNRNPC、WTAP、YTHDF3、FTO、ZC3H13、METTL14、METTL3、YTHDF1、YTHDF2这些基因是差异表达的。Kaplan-Meier生存分析显示,在这些差异表达的基因中仅仅HNRNPC和YTHDF2的表达与生存显著相关。然后,通过多因素Cox回归结果表明HNRNPC的表达在肺腺癌TCGA队列中是个独立危险因素。最后,HNRNPC在肺腺癌中的表达与临床分期(IV vs I, OR=3.692 308)和组织浸润(T2 vs T1, OR=1.776 471;T4 vs T1, OR=6.303 03)显著相关(所有p<0.05)。 结论认为HNRNPC可能作为肺腺癌的独立的预后因子。     

Prognostic signature and immune efficacy of m1A-, m5C- and m6A-related regulators in cutaneous melanoma

Cutaneous melanoma (CM) is an aggressive cancer; given that initial and specific signs are lacking, diagnosis is often late and the prognosis is poor. RNA modification has been widely studied in tumour progression. Nevertheless, little progress has been made in the signature of N¹-methyladenosine (m¹A), 5-methylcytosine (m⁵C), N⁶-methyladenosine (m⁶A)-related regulators and the tumour microenvironment (TME) cell infiltration in CM. Our study identified the characteristics of m¹A-, m⁵C- and m⁶A-related regulators based on 468 CM samples from the public database. Using univariate, multivariate and LASSO Cox regression analysis, a risk model of regulators was established and validated by a nomogram on independent prognostic factors. The gene set variation analysis (GSVA) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) clarified the involved functional pathways. A combined single-sample gene set enrichment analysis (ssGSEA) and CIBERSORT approach revealed TME of regulator-related prognostic signature. The nine-gene signature stratified the patients into distinct risk subgroups for personalized prognostic assessment. Additionally, functional enrichment, immune infiltration and immunotherapy response analysis indicated that the high-risk group was correlated with T-cell suppression, while the low-risk group was more sensitive to immunotherapy. The findings presented here contribute to our understanding of the TME molecular heterogeneity in CM. Nine m¹A-, m⁵C- and m⁶A-related regulators may also be promising biomarkers for future research.     

Identification of a novel 5-methylcytosine-related signature for prognostic prediction of kidney renal papillary cell carcinoma and a Putative target for drug repurposing

BackgroundMany studies have demonstrated the crucial roles of 5-methylcytosine (m5C) RNA methylation in cancer pathogenesis.MethodsTwo datasets, including TCGA-KIRP and ICGC, and related clinical information were downloaded, where the expression of 13 m5C regulators was examined. We applied LASSO regression to construct a multi-m5C-regulator-based signature in the TCGA cohort, which was further validated using the ICGC cohort. Univariate and multivariate Cox regressions were applied to evaluate the independent prognostic value of our model. The differences in biological functions and immune characterizations between high and low-risk groups divided based on the risk scores were also investigated via multiple approaches, such as enrichment analyses, mutation mining, and immune scoring. Finally, the sensitivities of commonly used targeted drugs were tested, and the connectivity MAP (cMAP) was utilized to screen potentially effective molecules for patients in the high-risk group. Experimental validation was done following qPCR tests in Caki-2 and HK-2 cell lines.Results3 m5C regulators, including ALYREF, DNMT3B and YBX1, were involved in our model. Survival analysis revealed a worse prognosis for patients in the high-risk group. Cox regression results indicated our model's superior predictive performance compared to single-factor prognostic evaluation. Functional enrichment analyses indicated a higher mutation frequency and poorer tumor microenvironment of patients in the high-risk group. qPCR-based results revealed that ALYREF, DNMT3B, and YBX1 were significantly up-regulated in Caki-2 cell lines compared with HK-2 cell lines. Molecules like BRD-K72451865, Levosimendan, and BRD-K03515135 were advised by cMAP for patients in the high-risk group.ConclusionOur study presented a novel predictive model for KIRP prognosis. Furthermore, the results of our analysis provide new insights for investigating m5C events in KIRP pathogenesis.     

High Bone Sialoprotein (BSP) Expression Correlates with Increased Tumor Grade and Predicts a Poorer Prognosis of High-Grade Glioma Patients

Objectives

To investigate the expression and prognostic value of bone sialoprotein (BSP) in glioma patients.

Methods

We determined the expression of BSP using real-time RT-PCR and immunohistochemistry in tissue microarrays containing 15 normal brain and 270 glioma samples. Cumulative survival was calculated by the Kaplan-Meier method and analyzed by the log-rank test. Univariate and multivariate analyses were performed by the stepwise forward Cox regression model.

Results

Both BSP mRNA and protein levels were significantly elevated in high-grade glioma tissues compared with those of normal brain and low-grade glioma tissues, and BSP expression positively correlated with tumor grade (P<0.001). Univariate and multivariate analysis showed high BSP expression was an independent prognostic factor for a shorter progression-free survival (PFS) and overall survival (OS) in both grade III and grade IV glioma patients [hazard ratio (HR) = 2.549 and 3.154 for grade III glioma, and HR = 1.637 and 1.574 for grade IV glioma, respectively]. Patients with low BSP expression had a significantly longer median OS and PFS than those with high BSP expression. Small extent of resection and lineage of astrocyte served as independent risk factors of both shorter PFS and OS in grade III glioma patients; GBM patients without O⁶-methylguanine (O⁶-meG) DNA methyltransferase (MGMT) methylation and Karnofsky performance score (KPS) less than 70 points were related to poor prognosis. Lack of radiotherapy related to shorter OS but not affect PFS in both grade III and grade IV glioma patients.

Conclusion

High BSP expression occurs in a significant subset of high-grade glioma patients and predicts a poorer outcome. The study identifies a potentially useful molecular marker for the categorization and targeted therapy of gliomas.