Mutation in WDR4 impairs tRNA m(7)G46 methylation and causes a distinct form of microcephalic primordial dwarfism

Construction of prognostic signature of breast cancer based on N7-Methylguanosine-Related LncRNAs and prediction of immune response

Background: Long non-coding RNA (LncRNA) is a prognostic factor for malignancies, and N7-Methylguanosine (m7G) is crucial in the occurrence and progression of tumors. However, it has not been documented how well m7G-related LncRNAs predict the development of breast cancer (BC). This study aims to develop a predictive signature based on long non-coding RNAs (LncRNAs) associated with m7G to predict the prognosis of breast cancer patients.

Methods: The Cancer Genome Atlas (TCGA) database provided us with the RNA-seq data and matching clinical information of individuals with breast cancer. To identify the signature of N7-Methylguanosine-Related LncRNAs and create a prognostic model, we employed co-expression network analysis, least absolute shrinkage selection operator (LASSO) regression analysis, univariate Cox regression analysis, and multivariate Cox regression analysis. The signature was assessed using the Kaplan-Meier analysis and Receiver Operating Characteristic (ROC) curve. A nomogram and principal component analysis (PCA) were employed to confirm the predictive signature’s usefulness. Then, we examined the drug sensitivity between the two risk groups and utilized single-sample gene set enrichment analysis (ssGSEA) to investigate the association between predictive factors and the tumor immune microenvironment in high-risk and low-risk groups.

Results: Nine m7G-related LncRNAs (LINC01871, AP003469.4, Z68871.1, AC245297.3, EGOT, TFAP2A-AS1, AL136531.1, SEMA3B-AS1, AL606834.2) that are independently associated with the overall survival time (OS) of BC patients make up the signature we developed. For predicting 1-, 3-, and 5-year survival rates, the areas under the ROC curve (AUC) were 0.715, 0.724, and 0.726, respectively. The Kaplan-Meier analysis revealed that the prognosis of BC patients in the high-risk group was worse than that of those in the low-risk group. When compared to clinicopathological variables, multiple regression analysis demonstrated that risk score was a significant independent predictive factor for BC patients. The results of the ssGSEA study revealed a substantial correlation between the predictive traits and the BC patients’ immunological status, low-risk BC patients had more active immune systems, and they responded better to PD1/L1 immunotherapy.

Conclusion: The prognostic signature, which is based on m7G-related LncRNAs, can be utilized to inform patients’ customized treatment plans by independently predicting their prognosis and how well they would respond to immunotherapy.

A novel prognostic signature based on N7-methylguanosine-related long non-coding RNAs in breast cancer

Long non-coding RNA (lncRNA) are closely associated with the occurrence and progression of tumors. However, research on N7-methylguanosine (m7G)-related lncRNA in breast cancer is lacking. Therefore, the present study explored the prognostic value, gene expression characteristics, and effects of m7G-related lncRNA on tumor immune cell infiltration and tumor mutational burden (TMB) in breast cancer. lncRNA expression matrices and clinical follow-up data of patients with breast cancer were obtained from The Cancer Genome Atlas, revealing eight significantly differentially expressed and prognostically relevant m7G-related lncRNAs in breast cancer tissues: BAIAP2-DT, COL4A2-AS1, FARP1-AS1, RERE-AS1, NDUFA6-DT, TFAP2A-AS1, LINC00115, and MIR302CHG. A breast cancer prognostic signature was created based on these m7G-related lncRNAs according to least absolute shrinkage and selection operator Cox regression. The prognostic signature combined with potential prognostic factors showed independent prognostic value, reliability, and specificity. Meanwhile, we constructed a risk score-based nomogram to assist clinical decision-making. Gene set enrichment analysis revealed that low- and high-risk group were associated with metabolism-related pathways. Our study demonstrated the association between tumor immune cell infiltration based on analyses with the CIBERSORT algorithm and prognostic signature. We also assessed the correlation between prognostic signature and TMB. Lastly, quantitative real-time polymerase chain reaction analysis was performed to validate differentially expressed lncRNAs. The effective prognostic signature based on m7G-related lncRNAs has the potential to predict the survival prognosis of patients with breast cancer. The eight m7G-related lncRNAs identified in this study might represent potential biomarkers and therapeutic targets of breast cancer.

Novel roles of METTL1/WDR4 in tumor via m7G methylation

As one of the prevalent posttranscriptional modifications of RNA, N⁷-methylguanosine (m⁷G) plays essential roles in RNA processing, metabolism, and function, mainly regulated by the methyltransferase-like 1 (METTL1) and WD repeat domain 4 (WDR4) complex. Emerging evidence suggests that the METTL1/WDR4 complex promoted or inhibited the processes of many tumors, including head and neck, lung, liver, colon, bladder cancer, and teratoma, dependent on close m⁷G methylation modification of tRNA or microRNA (miRNA). Therefore, METTL1 and m⁷G modification can be used as biomarkers or potential intervention targets, providing new possibilities for early diagnosis and treatment of tumors. This review will mainly focus on the mechanisms of METTL1/WDR4 via m⁷G in tumorigenesis and the corresponding detection methods.

METTL1 promotes neuroblastoma development through m7G tRNA modification and selective oncogenic gene translation

Background

Neuroblastoma (NBL) is the most common extra-cranial solid tumour in childhood, with prognosis ranging from spontaneous remission to high risk for rapid and fatal progression. Despite existing therapy approaches, the 5-year event-free survival (EFS) for patients with advanced NBL remains below 30%, emphasizing urgent necessary for novel therapeutic strategies. Studies have shown that epigenetic disorders play an essential role in the pathogenesis of NBL. However, the function and mechanism of N7-methylguanosine (m⁷G) methyltransferase in NBL remains unknown.

Methods

The expression levels of m⁷G tRNA methyltransferase Methyltransferase-like 1 (METTL1) were analyzed by querying the Gene Expression Omnibus (GEO) database and further confirmed by immunohistochemistry (IHC) assay. Kaplan-Meier, univariate and multivariate cox hazard analysis were performed to reveal the prognostic role of METTL1. Cell function assays were performed to evaluate how METTL1 works in proliferation, apoptosis and migration in cell lines and xenograft mouse models. The role of METTL1 on mRNA translation activity of NBL cells was measured using puromycin intake assay and polysome profiling assay. The m⁷G modified tRNAs were identified by tRNA reduction and cleavage sequencing (TRAC-seq). Ribosome nascent-chain complex-bound mRNA sequencing (RNC-seq) was utilized to identify the variation of gene translation efficiency (TE). Analyzed the codon frequency decoded by m⁷G tRNA to clarify the translation regulation and mechanism of m⁷G modification in NBL.

Results

This study found that METTL1 were significantly up-regulated in advanced NBL, which acted as an independent risk factor and predicted poor prognosis. Further in NBL cell lines and BALB/c-nu female mice, we found METTL1 played a crucial role in promoting NBL progression. Furthermore, m⁷G profiling and translation analysis revealed downregulation of METTL1 would inhibit puromycin intake efficiency of NBL cells, indicating that METTL1 did count crucially in regulation of NBL cell translation. With all tRNAs with m⁷G modification identified in NBL cells, knockdown of METTL1 would significantly reduce the levels of both m⁷G modification and m⁷G tRNAs expressions. Result of RNC-seq shew there were 339 overlapped genes with impaired translation in NBL cells upon METTL1 knockdown. Further analysis revealed these genes contained higher frequency of codons decoded by m⁷G-modified tRNAs and were enriched in oncogenic pathways.

Conclusion

This study revealed the critical role and mechanism of METTL1-mediated tRNA m⁷G modification in regulating NBL progression, providing new insights for developing therapeutic approaches for NBL patients.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40364-022-00414-z.

The role of RNA modification in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is a highly mortal type of primary liver cancer. Abnormal epigenetic modifications are present in HCC, and RNA modification is dynamic and reversible and is a key post-transcriptional regulator. With the in-depth study of post-transcriptional modifications, RNA modifications are aberrantly expressed in human cancers. Moreover, the regulators of RNA modifications can be used as potential targets for cancer therapy. In RNA modifications, N6-methyladenosine (m6A), N7-methylguanosine (m7G), and 5-methylcytosine (m5C) and their regulators have important regulatory roles in HCC progression and represent potential novel biomarkers for the confirmation of diagnosis and treatment of HCC. This review focuses on RNA modifications in HCC and the roles and mechanisms of m6A, m7G, m5C, N1-methyladenosine (m1A), N3-methylcytosine (m3C), and pseudouridine (ψ) on its development and maintenance. The potential therapeutic strategies of RNA modifications are elaborated for HCC.

The pattern of expression and prognostic value of key regulators for m7G RNA methylation in hepatocellular carcinoma

N7-methylguanosine (m7G) modification on internal RNA positions plays a vital role in several biological processes. Recent research shows m7G modification is associated with multiple cancers. However, in hepatocellular carcinoma (HCC), its implications remain to be determined. In this place, we need to interrogate the mRNA patterns for 29 key regulators of m7G RNA modification and assess their prognostic value in HCC. Initial, the details from The Cancer Genome Atlas (TCGA) database concerning transcribed gene data and clinical information of HCC patients were inspected systematically. Second, according to the mRNA profiles of 29 m7G RNA methylation regulators, two clusters (named 1 and 2, respectively) were identified by consensus clustering. Furthermore, robust risk signature for seven m7G RNA modification regulators was constructed. Last, we used the Gene Expression Omnibus (GEO) dataset to validate the prognostic associations of the seven-gene risk signature. We figured out that 24/29 key regulators of m7G RNA modification varied remarkably in their grades of expression between the HCC and the adjacent tumor control tissues. Cluster one compared with cluster two had a substandard prognosis and was also positively correlated with T classification (T), pathological stage, and vital status (fustat) significantly. Consensus clustering results suggested the expression pattern of m7G RNA modification regulators was correlated with the malignancy of HCC strongly. In addition, cluster one was extensively enriched in metabolic-related pathways. Seven optimal genes (METTL1, WDR4, NSUN2, EIF4E, EIF4E2, NCBP1, and NCBP2) were selected to establish the risk model for HCC. Indicating by further analyses and validation, the prognostic model has fine anticipating command and this probability signature might be a self supporting presage factor for HCC. Finally, a new prognostic nomogram based on age, gender, pathological stage, histological grade, and prospects were established to forecast the prognosis of HCC patients accurately. In essence, we detected association of HCC severity and expression levels of m7G RNA modification regulators, and developed a risk score model for predicting prognosis of HCC patients’ progression.

Molecular characterization, clinical relevance and immune feature of m7G regulator genes across 33 cancer types

Over 170 RNA modifications have been identified after transcriptions, involving in regulation of RNA splicing, processing, translation and decay. Growing evidence has unmasked the crucial role of N6-methyladenosine (m6A) in cancer development and progression, while, as a relative newly found RNA modification, N7-methylguanosine (m7G) is also certified to participate in tumorigenesis via different catalytic machinery from that of m6A. However, system analysis on m7G RNA modification-related regulator genes is lack. In this study, we first investigated the genetic alteration of m7G related regulator genes in 33 cancers, and found mRNA expression levels of most regulator genes were positively correlated with copy number variation (CNV) and negatively correlated with methylation in most cancers. We built a m7G RNA modification model based on the enrichment of the regulator gene scores to evaluate the m7G modification levels in 33 cancers, and investigated the connections of m7G scores to clinical outcomes. Furthermore, we paid close attention to the role of m7G in immunology due to the widely used immune checkpoint blockade therapy. Our results showed the higher m7G scores related to immunosuppression of tumor cells. Further confirmation with phase 3 clinical data with application of anti-PDL1/PDL indicated the impact of m7G modification level on immunotherapy effect. Relevance of m7G regulator genes and drug sensitivity was also evaluated to provide a better treatment choice when treating cancers. In summary, our study uncovered the profile of m7G RNA modification through various cancers, and figured out the connection of m7G modification levels with therapeutical outcomes, providing potential better options of cancer treatment.

The Epitranscriptome in miRNAs: Crosstalk, Detection, and Function in Cancer

The epitranscriptome encompasses all post-transcriptional modifications that occur on RNAs. These modifications can alter the function and regulation of their RNA targets, which, if dysregulated, result in various diseases and cancers. As with other RNAs, miRNAs are highly modified by epitranscriptomic modifications such as m6A methylation, 2'-O-methylation, m5C methylation, m7G methylation, polyuridine, and A-to-I editing. miRNAs are a class of small non-coding RNAs that regulates gene expression at the post-transcriptional level. miRNAs have gathered high clinical interest due to their role in disease, development, and cancer progression. Epitranscriptomic modifications alter the targeting, regulation, and biogenesis of miRNAs, increasing the complexity of miRNA regulation. In addition, emerging studies have revealed crosstalk between these modifications. In this review, we will summarize the epitranscriptomic modifications-focusing on those relevant to miRNAs-examine the recent crosstalk between these modifications, and give a perspective on how this crosstalk expands the complexity of miRNA biology.

N7-Methylguanosine Genes Related Prognostic Biomarker in Hepatocellular Carcinoma

Background: About 90% of liver cancer-related deaths are caused by hepatocellular carcinoma (HCC). N7-methylguanosine (m7G) modification is associated with the biological process and regulation of various diseases. To the best of our knowledge, its role in the pathogenesis and prognosis of HCC has not been thoroughly investigated. Aim: To identify N7-methylguanosine (m7G) related prognostic biomarkers in HCC. Furthermore, we also studied the association of m7G-related prognostic gene signature with immune infiltration in HCC. Methods: The TCGA datasets were used as a training and GEO dataset "GSE76427" for validation of the results. Statistical analyses were performed using the R statistical software version 4.1.2. Results: Functional enrichment analysis identified some pathogenesis related to HCC. We identified 3 m7G-related genes (CDK1, ANO1, and PDGFRA) as prognostic biomarkers for HCC. A risk score was calculated from these 3 prognostic m7G-related genes which showed the high-risk group had a significantly poorer prognosis than the low-risk group in both training and validation datasets. The 3- and 5-years overall survival was predicted better with the risk score than the ideal model in the entire cohort in the predictive nomogram. Furthermore, immune checkpoint genes like CTLA4, HAVCR2, LAG3, and TIGT were expressed significantly higher in the high-risk group and the chemotherapy sensitivity analysis showed that the high-risk groups were responsive to sorafenib treatment. Conclusion: These 3 m7G genes related signature model can be used as prognostic biomarkers in HCC and a guide for immunotherapy and chemotherapy response. Future clinical study on this biomarker model is required to verify its clinical implications.

Identification and Verification of m7G Modification Patterns and Characterization of Tumor Microenvironment Infiltration via Multi-Omics Analysis in Clear Cell Renal Cell Carcinoma

The epigenetic modification of tumorigenesis and progression in neoplasm has been demonstrated in recent studies. Nevertheless, the underlying association of N7-methylguanosine (m⁷G) regulation with molecular heterogeneity and tumor microenvironment (TME) in clear cell renal cell carcinoma (ccRCC) remains unknown. We explored the expression profiles and genetic variation features of m⁷G regulators and identified their correlations with patient outcomes in pan-cancer. Three distinct m⁷G modification patterns, including MGCS1, MGCS2, and MGCS3, were further determined and systematically characterized via multi-omics data in ccRCC. Compared with the other two subtypes, patients in MGCS3 exhibited a lower clinical stage/grade and better prognosis. MGCS1 showed the lowest enrichment of metabolic activities. MGCS2 was characterized by the suppression of immunity. We then established and validated a scoring tool named m7Sig, which could predict the prognosis of ccRCC patients. This study revealed that m⁷G modification played a vital role in the formation of the tumor microenvironment in ccRCC. Evaluating the m⁷G modification landscape helps us to raise awareness and strengthen the understanding of ccRCC’s characterization and, furthermore, to guide future clinical decision making.

The potential role of N7-methylguanosine (m7G) in cancer

N⁷-methylguanosine (m7G), one of the most prevalent RNA modifications, has recently attracted significant attention. The m7G modification actively participates in biological and pathological functions by affecting the metabolism of various RNA molecules, including messenger RNA, ribosomal RNA, microRNA, and transfer RNA. Increasing evidence indicates a critical role for m7G in human disease development, especially cancer, and aberrant m7G levels are closely associated with tumorigenesis and progression via regulation of the expression of multiple oncogenes and tumor suppressor genes. Currently, the underlying molecular mechanisms of m7G modification in cancer are not comprehensively understood. Here, we review the current knowledge regarding the potential function of m7G modifications in cancer and discuss future m7G-related diagnostic and therapeutic strategies.

Case Report: Biallelic Variant in the tRNA Methyltransferase Domain of the AlkB Homolog 8 Causes Syndromic Intellectual Disability

Intellectual disability (ID) has become very common and is an extremely heterogeneous disorder, where the patients face many challenges with deficits in intellectual functioning and adaptive behaviors. A single affected family revealed severe disease phenotypes such as ID, developmental delay, dysmorphic facial features, postaxial polydactyly type B, and speech impairment. DNA of a single affected individual was directly subjected to whole exome sequencing (WES), followed by Sanger sequencing. Data analysis revealed a novel biallelic missense variant (c.1511G>C; p.(Trp504Ser)) in the ALKBH8 gene, which plays a significant role in tRNA modifications. Our finding adds another variant to the growing list of ALKBH8-associated tRNA modifications causing ID and additional phenotypic manifestations. The present study depicts the key role of the genes associated with tRNA modifications, such as ALKBH8, in the development and pathophysiology of the human brain.

A N7-Methylguanine-Related Gene Signature Applicable for the Prognosis and Microenvironment of Prostate Cancer

Background

Despite the constant iteration of small-molecule inhibitors and immune checkpoint inhibitors, PRAD (prostate adenocarcinoma) patients with distant metastases and biochemical recurrence maintain a poor survival outcome along with an increasing morbidity in recent years. N⁷-Methylguanine, a new-found type of RNA modification, has demonstrated an essential role in tumor progression but has hardly been studied for its effect on prostate carcinoma. The current study aimed to seek m⁷G (N7-methylguanosine) related prognostic biomarkers and potential targets for PRAD treatment.

Methods

42 genes related to m⁷G were collected from former literatures and GSEA (Gene Set Enrichment Analysis) website. Then, RNA-seq (RNA sequencing) and clinical data from TCGA-PRAD (The Cancer Genome Atlas-Prostate) cohort were retrieved to screen the differentially expressed m⁷G genes to further construct a multivariate Cox prognostic model for PRAD. Next, GSE116918, a prostate cancer cohort acquired from GEO (Gene Expression Omnibus) database, was analyzed for the external validation group to assess the ability to predict BFFS (biochemical failure-free survival) of our m⁷G prognostic signature. Kaplan-Meier, ROC (receiver operator characteristic), AUC (areas under ROC curve), and calibration curves were adopted to display the performance of this prognostic signature. In addition, immune infiltration analysis was implemented to evaluate the effect of these m⁷G genes on immunoinfiltrating cells. Correlation with drug susceptibility of the m⁷G signature was also analyzed by matching drug information in CellMiner database.

Results

The m⁷G-related prognostic signature, including three genes (EIF3D, EIF4A1, LARP1) illustrated superior prognostic ability for PRAD in both training and validation cohorts. The 5-year AUC were 0.768 for TCGA-PRAD and 0.608 for GSE116918. It can well distinguish patients into different risk groups of biochemical recurrence (p =1e-04 for TCGA-PRAD and p =0.0186 for GSE116918). Immune infiltration analysis suggested potential regulation of m⁷G genes on neutrophils and dendritic cells in PRAD.

Conclusions

A m⁷G-related prognostic signature was constructed and validated in the current study, giving new sights of m⁷G methylation in predicting the prognostic and improving the treatment of PRAD.

Identification of a Trm732 Motif Required for 2'-O-methylation of the tRNA Anticodon Loop by Trm7

Posttranscriptional tRNA modifications are essential for proper gene expression, and defects in the enzymes that perform tRNA modifications are associated with numerous human disorders. Throughout eukaryotes, 2'-O-methylation of residues 32 and 34 of the anticodon loop of tRNA is important for proper translation, and in humans, a lack of these modifications results in non-syndromic X-linked intellectual disability. In yeast, the methyltransferase Trm7 forms a complex with Trm732 to 2'-O-methylate tRNA residue 32 and with Trm734 to 2'-O-methylate tRNA residue 34. Trm732 and Trm734 are required for the methylation activity of Trm7, but the role of these auxiliary proteins is not clear. Additionally, Trm732 and Trm734 homologs are implicated in biological processes not directly related to translation, suggesting that these proteins may have additional cellular functions. To identify critical amino acids in Trm732, we generated variants and tested their ability to function in yeast cells. We identified a conserved RRSAGLP motif in the conserved DUF2428 domain of Trm732 that is required for tRNA modification activity by both yeast Trm732 and its human homolog, THADA. The identification of Trm732 variants that lack tRNA modification activity will help to determine if other biological functions ascribed to Trm732 and THADA are directly due to tRNA modification or to secondary effects due to other functions of these proteins.

THUMPD1 bi-allelic variants cause loss of tRNA acetylation and a syndromic neurodevelopmental disorder

Covalent tRNA modifications play multi-faceted roles in tRNA stability, folding, and recognition, as well as the rate and fidelity of translation, and other cellular processes such as growth, development, and stress responses. Mutations in genes that are known to regulate tRNA modifications lead to a wide array of phenotypes and diseases including numerous cognitive and neurodevelopmental disorders, highlighting the critical role of tRNA modification in human disease. One such gene, THUMPD1, is involved in regulating tRNA N4-acetylcytidine modification (ac4C), and recently was proposed as a candidate gene for autosomal-recessive intellectual disability. Here, we present 13 individuals from 8 families who harbor rare loss-of-function variants in THUMPD1. Common phenotypic findings included global developmental delay, speech delay, moderate to severe intellectual deficiency, behavioral abnormalities such as angry outbursts, facial dysmorphism, and ophthalmological abnormalities. We demonstrate that the bi-allelic variants identified cause loss of function of THUMPD1 and that this defect results in a loss of ac4C modification in small RNAs, and of individually purified tRNA-Ser-CGA. We further corroborate this effect by showing a loss of tRNA acetylation in two CRISPR-Cas9-generated THUMPD1 KO cell lines. In addition, we also show the resultant amino acid substitution that occurs in a missense THUMPD1 allele identified in an individual with compound heterozygous variants results in a marked decrease in THUMPD1 stability and RNA-binding capacity. Taken together, these results suggest that the lack of tRNA acetylation due to THUMPD1 loss of function results in a syndromic form of intellectual disability associated with developmental delay, behavioral abnormalities, hearing loss, and facial dysmorphism.

N7-methylguanosine tRNA modification promotes tumorigenesis and chemoresistance through WNT/β-catenin pathway in nasopharyngeal carcinoma

Treatment selections are very limited for patients with advanced nasopharyngeal carcinoma (NPC) experiencing disease progression. Uncovering mechanisms underlying NPC progression is crucial for the development of novel treatments. Here we show that N7-methylguanosine (m7G) tRNA modification enzyme METTL1 and its partner WDR4 are significantly elevated in NPC and are associated with poor prognosis. Loss-of-function and gain-of-function assays demonstrated that METTL1/WDR4 promotes NPC growth and metastasis in vitro and in vivo. Mechanistically, ARNT was identified as an upstream transcription factor regulating METTL1 expression in NPC. METTL1 depletion resulted in decreased m7G tRNA modification and expression, which led to impaired codon recognition during mRNA translation, therefore reducing the translation efficiencies of mRNAs with higher m7G codons. METTL1 upregulated the WNT/β-catenin signaling pathway and promoted NPC cell epithelial-mesenchymal transition (EMT) and chemoresistance to cisplatin and docetaxel in vitro and in vivo. Overexpression of WNT3A bypassed the requirement of METTL1 for EMT and chemoresistance. This work uncovers novel insights into tRNA modification-mediated mRNA translation regulation and highlights the critical function of tRNA modification in cancer progression.

N7-methylguanosine tRNA modification promotes esophageal squamous cell carcinoma tumorigenesis via the RPTOR/ULK1/autophagy axis

Mis-regulated RNA modifications promote the processing and translation of oncogenic mRNAs to facilitate cancer progression, while the molecular mechanisms remain unclear. Here we reveal that tRNA m⁷G methyltransferase complex proteins METTL1 and WDR4 are significantly up-regulated in esophageal squamous cell carcinoma (ESCC) tissues and associated with poor ESCC prognosis. In addition, METTL1 and WDR4 promote ESCC progression via the tRNA m⁷G methyltransferase activity in vitro and in vivo. Mechanistically, METTL1 or WDR4 knockdown leads to decreased expression of m⁷G-modified tRNAs and reduces the translation of a subset of oncogenic transcripts enriched in RPTOR/ULK1/autophagy pathway. Furthermore, ESCC models using Mettl1 conditional knockout and knockin mice uncover the essential function of METTL1 in promoting ESCC tumorigenesis in vivo. Our study demonstrates the important oncogenic function of mis-regulated tRNA m⁷G modification in ESCC, and suggest that targeting METTL1 and its downstream signaling axis could be a promising therapeutic target for ESCC treatment.

Deregulation of METTL1-mediated N7- methylguanosine tRNA modification can promote oncogenesis. Here, the authors report that this modification regulates the translation of proteins in both the mTOR and negative regulators of autophagy pathways, resulting in the progression of esophageal squamous cell carcinoma.

RNA nucleotide methylation: 2021 update

Among RNA modifications, transfer of methylgroups from the typical cofactor S-adenosyl-l-methionine by methyltransferases (MTases) to RNA is by far the most common reaction. Since our last review about a decade ago, the field has witnessed the re-emergence of mRNA methylation as an important mechanism in gene regulation. Attention has then spread to many other RNA species; all being included into the newly coined concept of the "epitranscriptome." The focus moved from prokaryotes and single cell eukaryotes as model organisms to higher eukaryotes, in particular to mammals. The perception of the field has dramatically changed over the past decade. A previous lack of phenotypes in knockouts in single cell organisms has been replaced by the apparition of MTases in numerous disease models and clinical investigations. Major driving forces of the field include methylation mapping techniques, as well as the characterization of the various MTases, termed "writers." The latter term has spilled over from DNA modification in the neighboring epigenetics field, along with the designations "readers," applied to mediators of biological effects upon specific binding to a methylated RNA. Furthermore "eraser" enzymes effect the newly discovered oxidative removal of methylgroups. A sense of reversibility and dynamics has replaced the older perception of RNA modification as a concrete-cast, irreversible part of RNA maturation. A related concept concerns incompletely methylated residues, which, through permutation of each site, lead to inhomogeneous populations of numerous modivariants. This review recapitulates the major developments of the past decade outlined above, and attempts a prediction of upcoming trends. This article is categorized under: RNA Processing > RNA Editing and Modification.

METTL1 promotes hepatocarcinogenesis via m7 G tRNA modification-dependent translation control

BACKGROUND

N7 -methylguanosine (m7 G) modification is one of the most common transfer RNA (tRNA) modifications in humans. The precise function and molecular mechanism of m7 G tRNA modification in hepatocellular carcinoma (HCC) remain poorly understood.

METHODS

The prognostic value and expression level of m7 G tRNA methyltransferase complex components methyltransferase-like protein-1 (METTL1) and WD repeat domain 4 (WDR4) in HCC were evaluated using clinical samples and TCGA data. The biological functions and mechanisms of m7 G tRNA modification in HCC progression were studied in vitro and in vivo using cell culture, xenograft model, knockin and knockout mouse models. The m7 G reduction and cleavage sequencing (TRAC-seq), polysome profiling and polyribosome-associated mRNA sequencing methods were used to study the levels of m7 G tRNA modification, tRNA expression and mRNA translation efficiency.

RESULTS

The levels of METTL1 and WDR4 are elevated in HCC and associated with advanced tumour stages and poor patient survival. Functionally, silencing METTL1 or WDR4 inhibits HCC cell proliferation, migration and invasion, while forced expression of wild-type METTL1 but not its catalytic dead mutant promotes HCC progression. Knockdown of METTL1 reduces m7 G tRNA modification and decreases m7 G-modified tRNA expression in HCC cells. Mechanistically, METTL1-mediated tRNA m7 G modification promotes the translation of target mRNAs with higher frequencies of m7 G-related codons. Furthermore, in vivo studies with Mettl1 knockin and conditional knockout mice reveal the essential physiological function of Mettl1 in hepatocarcinogenesis using hydrodynamics transfection HCC model.

CONCLUSIONS

Our work reveals new insights into the role of the misregulated tRNA modifications in liver cancer and provides molecular basis for HCC diagnosis and treatment.

METTL1-m7 G-EGFR/EFEMP1 axis promotes the bladder cancer development

BACKGROUND

The posttranscriptional modifications of transfer RNA (tRNA) are critical for all aspects of the tRNA function and have been implicated in the tumourigenesis and progression of many human cancers. By contrast, the biological functions of methyltransferase-like 1 (METTL1)-regulated m7 G tRNA modification in bladder cancer (BC) remain obscure.

RESULTS

In this research, we show that METTL1 was highly expressed in BC, and its level was correlated with poor patient prognosis. Silencing METTL1 suppresses the proliferation, migration and invasion of BC cells in vitro and in vivo. Multi-omics analysis reveals that METTL1-mediated m7 G tRNA modification altered expression of certain target genes, including EGFR/EFEMP1. Mechanistically, METTL1 regulates the translation of EGFR/EFEMP1 via modifying certain tRNAs. Furthermore, forced expression of EGFR/EFEMP1 partially rescues the effect of METTL1 deletion on BC cells.

CONCLUSIONS

Our findings demonstrate the oncogenic role of METTL1 and the pathological significance of the METTL1-m7 G-EGFR/EFEMP1 axis in the BC development, thus providing potential therapeutic targets for the BC treatment.

Insight into ALKBH8-related intellectual developmental disability based on the first pathogenic missense variant

ALKBH8 is a methyltransferase that modifies tRNAs by methylating the anticodon wobble uridine residue. The syndrome of ALKBH8-related intellectual developmental disability (MRT71) has thus far been reported solely in the context of homozygous truncating variants that cluster in the last exon. This raises interesting questions about the disease mechanism, because these variants are predicted to escape nonsense mediated decay and yet they appear to be loss of function. Furthermore, the limited class of reported variants complicates the future interpretation of missense variants in ALKBH8. Here, we report a consanguineous family in which two children with MRT71-compatible phenotype are homozygous for a novel missense variant in the methyltransferase domain. We confirm the pathogenicity of this variant by demonstrating complete absence of ALKBH8-dependent modifications in patient cells. Targeted proteomics analysis of ALKBH8 indicates that the variant does not lead to loss of ALKBH8 protein expression. This report adds to the clinical delineation of MRT71, confirms loss of function of ALKBH8 as the disease mechanism and expands the repertoire of its molecular lesions.

Advances in RNA Epigenetic Modifications in Hepatocellular Carcinoma and Potential Targeted Intervention Strategies

The current interventions for hepatocellular carcinoma (HCC) are not satisfactory, and more precise targets and promising strategies need to be explored. Recent research has demonstrated the non-negligible roles of RNA epigenetic modifications such as N6-methyladenosine (m6A) and 5-methylcytosine (m5C) in various cancers, including HCC. However, the specific targeting mechanisms are not well elucidated. In this review, we focus on the occurrence and detailed physiopathological roles of multiple RNA modifications on diverse RNAs closely related to the HCC process. In particular, we highlight fresh insights into the impact mechanisms of these posttranscriptional modifications on the whole progression of HCC. Furthermore, we analyzed the possibilities and significance of these modifications and regulators as potential therapeutic targets in HCC treatment, which provides the foundation for exploring targeted intervention strategies. This review will propel the identification of promising therapeutic targets and novel strategies that can be translated into clinical applications for HCC treatment.

METTL1/WDR4 mediated m7G tRNA modifications and m7G codon usage promote mRNA translation and lung cancer progression

Mis-regulated epigenetic modifications in RNAs are associated with human cancers. The transfer RNAs (tRNAs) are the most heavily modified RNA species in cells, however, little is known about the functions of tRNA modifications in cancers. In this study, we uncovered that the expression levels of tRNA N⁷-methylguanosine (m⁷G) methyltransferase complex components METTL1 and WDR4 are significantly elevated in human lung cancer samples and negatively associated with patient prognosis. Impaired m⁷G tRNA modification upon METTL1/WDR4 depletion resulted in decreased cell proliferation, colony formation, cell invasion and impaired tumorigenic capacities of lung cancer cells in vitro and in vivo. Moreover, gain-of-function and mutagenesis experiments revealed that METTL1 promoted lung cancer growth and invasion through regulation of m⁷G tRNA modifications. Profiling of tRNA methylation and mRNA translation revealed that highly translated mRNAs have higher frequencies of m⁷G tRNA decoded codons and knockdown of METTL1 resulted in decreased translation of mRNAs with higher frequencies of m⁷G tRNA codons, suggesting that tRNA modifications and codon usage play essential function in mRNA translation regulation. Our data uncovered novel insights on mRNA translation regulation through tRNA modifications and the corresponding mRNA codon compositions in lung cancer, providing new molecular basis underlying lung cancer progression.

Genetic Association of the Functional WDR4 Gene in Male Fertility

Infertility is one of the important problems in the modern world. Male infertility is characterized by several clinical manifestations, including low sperm production (oligozoospermia), reduced sperm motility (asthenozoospermia), and abnormal sperm morphology (teratozoospermia). WDR4, known as Wuho, controls fertility in Drosophila. However, it is unclear whether WDR4 is associated with clinical manifestations of male fertility in human. Here, we attempted to determine the physiological functions of WDR4 gene. Two cohorts were applied to address this question. The first cohort was the general population from Taiwan Biobank. Genomic profiles from 68,948 individuals and 87 common physiological traits were applied for phenome-wide association studies (PheWAS). The second cohort comprised patients with male infertility from Wan Fang Hospital, Taipei Medical University. In total, 81 male participants were recruited for the genetic association study. Clinical records including gender, age, total testosterone, follicle-stimulating hormone (FSH), luteinizing hormone (LH), total sperm number, sperm motility, and sperm morphology were collected. In the first cohort, results from PheWAS exhibited no associations between WDR4 genetic variants and 87 common physiological traits. In the second cohort, a total of four tagging single-nucleotide polymorphisms (tSNPs) from WDR4 gene (rs2298666, rs465663, rs2248490, and rs3746939) were selected for genotyping. We found that SNP rs465663 solely associated with asthenozoospermia. Functional annotations through the GTEx portal revealed the correlation between TT or TC genotype and low expression of WDR4. Furthermore, we used mouse embryonic fibroblasts cells from mwdr4 heterozygous (+/‒) mice for functional validation by western blotting. Indeed, low expression of WDR4 contributed to ROS-induced DNA fragmentation. In conclusion, our results suggest a critical role of WDR4 gene variant as well as protein expression in asthenozoospermia.

N7-Methylguanosine tRNA modification enhances oncogenic mRNA translation and promotes intrahepatic cholangiocarcinoma progression

Cancer cells selectively promote translation of specific oncogenic transcripts to facilitate cancer survival and progression, but the underlying mechanisms are poorly understood. Here, we find that N⁷-methylguanosine (m⁷G) tRNA modification and its methyltransferase complex components, METTL1 and WDR4, are significantly upregulated in intrahepatic cholangiocarcinoma (ICC) and associated with poor prognosis. We further reveal the critical role of METTL1/WDR4 in promoting ICC cell survival and progression using loss- and gain-of-function assays in vitro and in vivo. Mechanistically, m⁷G tRNA modification selectively regulates the translation of oncogenic transcripts, including cell-cycle and epidermal growth factor receptor (EGFR) pathway genes, in m⁷G-tRNA-decoded codon-frequency-dependent mechanisms. Moreover, using overexpression and knockout mouse models, we demonstrate the crucial oncogenic function of Mettl1-mediated m⁷G tRNA modification in promoting ICC tumorigenesis and progression in vivo. Our study uncovers the important physiological function and mechanism of METTL1-mediated m⁷G tRNA modification in the regulation of oncogenic mRNA translation and cancer progression.

Multi-omics integration of methyltransferase-like protein family reveals clinical outcomes and functional signatures in human cancer

Human methyltransferase-like (METTL) proteins transfer methyl groups to nucleic acids, proteins, lipids, and other small molecules, subsequently playing important roles in various cellular processes. In this study, we performed integrated genomic, transcriptomic, proteomic, and clinicopathological analyses of 34 METTLs in a large cohort of primary tumor and cell line data. We identified a subset of METTL genes, notably METTL1, METTL7B, and NTMT1, with high frequencies of genomic amplification and/or up-regulation at both the mRNA and protein levels in a spectrum of human cancers. Higher METTL1 expression was associated with high-grade tumors and poor disease prognosis. Loss-of-function analysis in tumor cell lines indicated the biological importance of METTL1, an m⁷G methyltransferase, in cancer cell growth and survival. Furthermore, functional annotation and pathway analysis of METTL1-associated proteins revealed that, in addition to the METTL1 cofactor WDR4, RNA regulators and DNA packaging complexes may be functionally interconnected with METTL1 in human cancer. Finally, we generated a crystal structure model of the METTL1–WDR4 heterodimeric complex that might aid in understanding the key functional residues. Our results provide new information for further functional study of some METTL alterations in human cancer and might lead to the development of small inhibitors that target cancer-promoting METTLs.

tRNA-derived fragments: A new class of non-coding RNA with key roles in nervous system function and dysfunction

tRNA-derived small RNAs (tsRNA) are a recently identified family of non-coding RNA that have been associated with a variety of cellular functions including the regulation of protein translation and gene expression. Recent sequencing and bioinformatic studies have identified the broad spectrum of tsRNA in the nervous system and demonstrated that this new class of non-coding RNA is produced from tRNA by specific cleavage events catalysed by ribonucleases such as angiogenin and dicer. Evidence is also accumulating that production of tsRNA is increased during disease processes where they regulate stress responses, proteostasis, and neuronal survival. Mutations to tRNA cleaving and modifying enzymes have been implicated in several neurodegenerative disorders, and tsRNA levels in the blood are advancing as biomarkers for neurological disease. In this review we summarize the physiological importance of tsRNA in the central nervous system and their relevance to neurological disease.

Division of labor in epitranscriptomics: What have we learnt from the structures of eukaryotic and viral multimeric RNA methyltransferases?

The translation of an mRNA template into the corresponding protein is a highly complex and regulated choreography performed by ribosomes, tRNAs, and translation factors. Most RNAs involved in this process are decorated by multiple chemical modifications (known as epitranscriptomic marks) contributing to the efficiency, the fidelity, and the regulation of the mRNA translation process. Many of these epitranscriptomic marks are written by holoenzymes made of a catalytic subunit associated with an activating subunit. These holoenzymes play critical roles in cell development. Indeed, several mutations being identified in the genes encoding for those proteins are linked to human pathologies such as cancers and intellectual disorders for instance. This review describes the structural and functional properties of RNA methyltransferase holoenzymes, which when mutated often result in brain development pathologies. It illustrates how structurally different activating subunits contribute to the catalytic activity of these holoenzymes through common mechanistic trends that most likely apply to other classes of holoenzymes. This article is categorized under: RNA Processing > RNA Editing and Modification RNA Processing > Capping and 5' End Modifications.

Structural model of the M7G46 Methyltransferase TrmB in complex with tRNA

TrmB belongs to the class I S-adenosylmethionine (SAM)-dependent methyltransferases (MTases) and introduces a methyl group to guanine at position 7 (m⁷G) in tRNA. In tRNAs m⁷G is most frequently found at position 46 in the variable loop and forms a tertiary base pair with C13 and U22, introducing a positive charge at G46. The TrmB/Trm8 enzyme family is structurally diverse, as TrmB proteins exist in a monomeric, homodimeric, and heterodimeric form. So far, the exact enzymatic mechanism, as well as the tRNA-TrmB crystal structure is not known. Here we present the first crystal structures of B. subtilis TrmB in complex with SAM and SAH. The crystal structures of TrmB apo and in complex with SAM and SAH have been determined by X-ray crystallography to 1.9 Å (apo), 2.5 Å (SAM), and 3.1 Å (SAH). The obtained crystal structures revealed Tyr193 to be important during SAM binding and MTase activity. Applying fluorescence polarization, the dissociation constant K_d of TrmB and tRNA^Phe was determined to be 0.12 µM ± 0.002 µM. Luminescence-based methyltransferase activity assays revealed cooperative effects during TrmB catalysis with half-of-the-site reactivity at physiological SAM concentrations. Structural data retrieved from small-angle x-ray scattering (SAXS), mass-spectrometry of cross-linked complexes, and molecular docking experiments led to the determination of the TrmB-tRNA^Phe complex structure.

m7GDisAI: N7-methylguanosine (m7G) sites and diseases associations inference based on heterogeneous network

Background

Recent studies have confirmed that N7-methylguanosine (m⁷G) modification plays an important role in regulating various biological processes and has associations with multiple diseases. Wet-lab experiments are cost and time ineffective for the identification of disease-associated m⁷G sites. To date, tens of thousands of m⁷G sites have been identified by high-throughput sequencing approaches and the information is publicly available in bioinformatics databases, which can be leveraged to predict potential disease-associated m⁷G sites using a computational perspective. Thus, computational methods for m⁷G-disease association prediction are urgently needed, but none are currently available at present.

Results

To fill this gap, we collected association information between m⁷G sites and diseases, genomic information of m⁷G sites, and phenotypic information of diseases from different databases to build an m⁷G-disease association dataset. To infer potential disease-associated m⁷G sites, we then proposed a heterogeneous network-based model, m⁷G Sites and Diseases Associations Inference (m⁷GDisAI) model. m⁷GDisAI predicts the potential disease-associated m⁷G sites by applying a matrix decomposition method on heterogeneous networks which integrate comprehensive similarity information of m⁷G sites and diseases. To evaluate the prediction performance, 10 runs of tenfold cross validation were first conducted, and m⁷GDisAI got the highest AUC of 0.740(± 0.0024). Then global and local leave-one-out cross validation (LOOCV) experiments were implemented to evaluate the model’s accuracy in global and local situations respectively. AUC of 0.769 was achieved in global LOOCV, while 0.635 in local LOOCV. A case study was finally conducted to identify the most promising ovarian cancer-related m⁷G sites for further functional analysis. Gene Ontology (GO) enrichment analysis was performed to explore the complex associations between host gene of m⁷G sites and GO terms. The results showed that m⁷GDisAI identified disease-associated m⁷G sites and their host genes are consistently related to the pathogenesis of ovarian cancer, which may provide some clues for pathogenesis of diseases.

Conclusion

The m⁷GDisAI web server can be accessed at http://180.208.58.66/m7GDisAI/, which provides a user-friendly interface to query disease associated m⁷G. The list of top 20 m⁷G sites predicted to be associted with 177 diseases can be achieved. Furthermore, detailed information about specific m⁷G sites and diseases are also shown.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04007-9.

ASHG 2020 Curt Stern Award introduction: Fowzan Sami Alkuraya

This article is based on the address given by the author at the 2020 virtual meeting of the American Society of Human Genetics (ASHG) on October 26, 2020. The video of the original address can be found at the ASHG website.

HN-CNN: A Heterogeneous Network Based on Convolutional Neural Network for m7 G Site Disease Association Prediction

N⁷-methylguanosine (m⁷G) is a typical positively charged RNA modification, playing a vital role in transcriptional regulation. m⁷G can affect the biological processes of mRNA and tRNA and has associations with multiple diseases including cancers. Wet-lab experiments are cost and time ineffective for the identification of disease-related m⁷G sites. Thus, a heterogeneous network method based on Convolutional Neural Networks (HN-CNN) has been proposed to predict unknown associations between m⁷G sites and diseases. HN-CNN constructs a heterogeneous network with m⁷G site similarity, disease similarity, and disease-associated m⁷G sites to formulate features for m⁷G site-disease pairs. Next, a convolutional neural network (CNN) obtains multidimensional and irrelevant features prominently. Finally, XGBoost is adopted to predict the association between m⁷G sites and diseases. The performance of HN-CNN is compared with Naive Bayes (NB), Random Forest (RF), Support Vector Machine (SVM), as well as Gradient Boosting Decision Tree (GBDT) through 10-fold cross-validation. The average AUC of HN-CNN is 0.827, which is superior to others.

The expanding world of tRNA modifications and their disease relevance

Transfer RNA (tRNA) is an adapter molecule that links a specific codon in mRNA with its corresponding amino acid during protein synthesis. tRNAs are enzymatically modified post-transcriptionally. A wide variety of tRNA modifications are found in the tRNA anticodon, which are crucial for precise codon recognition and reading frame maintenance, thereby ensuring accurate and efficient protein synthesis. In addition, tRNA-body regions are also frequently modified and thus stabilized in the cell. Over the past two decades, 16 novel tRNA modifications were discovered in various organisms, and the chemical space of tRNA modification continues to expand. Recent studies have revealed that tRNA modifications can be dynamically altered in response to levels of cellular metabolites and environmental stresses. Importantly, we now understand that deficiencies in tRNA modification can have pathological consequences, which are termed 'RNA modopathies'. Dysregulation of tRNA modification is involved in mitochondrial diseases, neurological disorders and cancer.

Human transfer RNA modopathies: diseases caused by aberrations in transfer RNA modifications

tRNA molecules are post-transcriptionally modified by tRNA modification enzymes. Although composed of different chemistries, more than 40 types of human tRNA modifications play pivotal roles in protein synthesis by regulating tRNA structure and stability as well as decoding genetic information on mRNA. Many tRNA modifications are conserved among all three kingdoms of life, and aberrations in various human tRNA modification enzymes cause life-threatening diseases. Here, we describe the class of diseases and disorders caused by aberrations in tRNA modifications as 'tRNA modopathies'. Aberrations in over 50 tRNA modification enzymes are associated with tRNA modopathies, which most frequently manifest as dysfunctions of the brain and/or kidney, mitochondrial diseases, and cancer. However, the molecular mechanisms that link aberrant tRNA modifications to human diseases are largely unknown. In this review, we provide a comprehensive compilation of human tRNA modification functions, tRNA modification enzyme genes, and tRNA modopathies, and we summarize the elucidated pathogenic mechanisms underlying several tRNA modopathies. We will also discuss important questions that need to be addressed in order to understand the molecular pathogenesis of tRNA modopathies.

Transcriptome-Wide Detection of Internal N7-Methylguanosine

m⁷G-seq detects internal 7-methylguanosine (m⁷G) sites within mRNAs and noncoding RNAs by misincorporation signatures. A chemical-assisted sequencing approach selectively converts internal m⁷G sites into abasic sites, triggering misincorporation at these sites in the presence of a specific reverse transcriptase. The further enrichment of m⁷G-induced abasic sites by biotin pull-down reveals hundreds of internal m⁷G sites in human mRNA. The misincorporation ratio before pull-down enrichment can be used for estimating the methylation fraction of some highly methylated m⁷G sites.

A novel assay provides insight into tRNAPhe retrograde nuclear import and re-export in S. cerevisiae

In eukaryotes, tRNAs are transcribed in the nucleus and subsequently exported to the cytoplasm where they serve as essential adaptor molecules in translation. However, tRNAs can be returned to the nucleus by the evolutionarily conserved process called tRNA retrograde nuclear import, before relocalization back to the cytoplasm via a nuclear re-export step. Several important functions of these latter two trafficking events have been identified, yet the pathways are largely unknown. Therefore, we developed an assay in Saccharomyces cerevisiae to identify proteins mediating tRNA retrograde nuclear import and re-export using the unique wybutosine modification of mature tRNAPhe. Our hydrochloric acid/aniline assay revealed that the karyopherin Mtr10 mediates retrograde import of tRNAPhe, constitutively and in response to amino acid deprivation, whereas the Hsp70 protein Ssa2 mediates import specifically in the latter. Furthermore, tRNAPhe is re-exported by Crm1 and Mex67, but not by the canonical tRNA exporters Los1 or Msn5. These findings indicate that the re-export process occurs in a tRNA family-specific manner. Together, this assay provides insights into the pathways for tRNAPhe retrograde import and re-export and is a tool that can be used on a genome-wide level to identify additional gene products involved in these tRNA trafficking events.

Next generation sequencing reveals novel homozygous frameshift in PUS7 and splice acceptor variants in AASS gene leading to intellectual disability, developmental delay, dysmorphic feature and microcephaly

Intellectual developmental disorder with abnormal behavior, microcephaly and short stature (IDDABS), (OMIM# 618342) is an autosomal recessive condition described as developmental delay, poor or absent speech, intellectual disability, short stature, mild to progressive microcephaly, delayed psychomotor development, hyperactivity, seizure, along with mild to swear aggressive behavior. Homozygous frameshift mutation in Pseudouridine Synthase 7, Putative; (PUS7) OMIM# 616,261 NM_019042.3 and splice acceptor variants in Alpha-Aminoadipic Semialdehyde Synthase; (AASS) OMIM# 605,113 NM_005763.3 was funded. Whole exome sequencing (WES) technique was used as tool to identify the molecular diagnostic test. Different bioinformatics analysis done for WES data and we identified two novel mutations one as frameshift mutation c.606_607delGA, p.Ser282CysfsTer9 in the PUS7 gene and splice acceptor variants c.1767–1 G > A in the AASS gene has been reported. The pattern of family segregation maintained the pathogenicity of this variation associated with abnormal behavior, intellectual developmental disorder, microcephaly along with short stature IDDABS. Further, the WES data was validated in the family having other affected individuals and healthy controls (n = 100) was done using Sanger sequencing. Finally, our results further explained the role of WES in the disease diagnosis and elucidated that the mutation in PUS7 and AASS genes may lead an important role for the development of IDDABS in Saudi family.

WD40-Repeat Proteins in Ciliopathies and Congenital Disorders of Endocrine System

WD40-repeat (WDR)-containing proteins constitute an evolutionarily conserved large protein family with a broad range of biological functions. In human proteome, WDR makes up one of the most abundant protein-protein interaction domains. Members of the WDR protein family play important roles in nearly all major cellular signalling pathways. Mutations of WDR proteins have been associated with various human pathologies including neurological disorders, cancer, obesity, ciliopathies and endocrine disorders. This review provides an updated overview of the biological functions of WDR proteins and their mutations found in congenital disorders. We also highlight the significant role of WDR proteins in ciliopathies and endocrine disorders. The new insights may help develop therapeutic approaches targeting WDR motifs.

Hypomodified tRNA in evolutionarily distant yeasts can trigger rapid tRNA decay to activate the general amino acid control response, but with different consequences

All tRNAs are extensively modified, and modification deficiency often results in growth defects in the budding yeast Saccharomyces cerevisiae and neurological or other disorders in humans. In S. cerevisiae, lack of any of several tRNA body modifications results in rapid tRNA decay (RTD) of certain mature tRNAs by the 5'-3' exonucleases Rat1 and Xrn1. As tRNA quality control decay mechanisms are not extensively studied in other eukaryotes, we studied trm8Δ mutants in the evolutionarily distant fission yeast Schizosaccharomyces pombe, which lack 7-methylguanosine at G46 (m7G46) of their tRNAs. We report here that S. pombe trm8Δ mutants are temperature sensitive primarily due to decay of tRNATyr(GUA) and that spontaneous mutations in the RAT1 ortholog dhp1+ restored temperature resistance and prevented tRNA decay, demonstrating conservation of the RTD pathway. We also report for the first time evidence linking the RTD and the general amino acid control (GAAC) pathways, which we show in both S. pombe and S. cerevisiae. In S. pombe trm8Δ mutants, spontaneous GAAC mutations restored temperature resistance and tRNA levels, and the trm8Δ temperature sensitivity was precisely linked to GAAC activation due to tRNATyr(GUA) decay. Similarly, in the well-studied S. cerevisiae trm8Δ trm4Δ RTD mutant, temperature sensitivity was closely linked to GAAC activation due to tRNAVal(AAC) decay; however, in S. cerevisiae, GAAC mutations increased tRNA loss and exacerbated temperature sensitivity. A similar exacerbated growth defect occurred upon GAAC mutation in S. cerevisiae trm8Δ and other single modification mutants that triggered RTD. Thus, these results demonstrate a conserved GAAC activation coincident with RTD in S. pombe and S. cerevisiae, but an opposite impact of the GAAC response in the two organisms. We speculate that the RTD pathway and its regulation of the GAAC pathway is widely conserved in eukaryotes, extending to other mutants affecting tRNA body modifications.

METTL1-mediated m7G methylation maintains pluripotency in human stem cells and limits mesoderm differentiation and vascular development

Background

7-Methylguanosine (m⁷G) is one of the most conserved modifications in nucleosides within tRNAs and rRNAs. It plays essential roles in the regulation of mRNA export, splicing, and translation. Recent studies highlighted the importance of METTL1-mediated m⁷G tRNA methylome in the self-renewal of mouse embryonic stem cells (mESCs) through its ability to regulate mRNA translation. However, the exact mechanisms by which METTL1 regulates pluripotency and differentiation in human induced pluripotent stem cells (hiPSCs) remain unknown. In this study, we evaluated the functions and underlying molecular mechanisms of METTL1 in regulating hiPSC self-renewal and differentiation in vivo and in vitro.

Methods

By establishing METTL1 knockdown (KD) hiPSCs, gene expression profiling was performed by RNA sequencing followed by pathway analyses. Anti-m⁷G northwestern assay was used to identify m⁷G modifications in tRNAs and mRNAs. Polysome profiling was used to assess the translation efficiency of the major pluripotent transcription factors. Moreover, the in vitro and in vivo differentiation capacities of METTL1-KD hiPSCs were assessed in embryoid body (EB) formation and teratoma formation assays.

Results

METTL1 silencing resulted in alterations in the global m⁷G profile in hiPSCs and reduced the translational efficiency of stem cell marker genes. METTL1-KD hiPSCs exhibited reduced pluripotency with slower cell cycling. Moreover, METTL1 silencing accelerates hiPSC differentiation into EBs and promotes the expression of mesoderm-related genes. Similarly, METTL1 knockdown enhances teratoma formation and mesoderm differentiation in vivo by promoting cell proliferation and angiogenesis in nude mice.

Conclusion

Our findings provided novel insight into the critical role of METTL1-mediated m⁷G modification in the regulation of hiPSC pluripotency and differentiation, as well as its potential roles in vascular development and the treatment of vascular diseases.

m7GHub: deciphering the location, regulation and pathogenesis of internal mRNA N7-methylguanosine (m7G) sites in human

MOTIVATION

Recent progress in N7-methylguanosine (m7G) RNA methylation studies has focused on its internal (rather than capped) presence within mRNAs. Tens of thousands of internal mRNA m7G sites have been identified within mammalian transcriptomes, and a single resource to best share, annotate and analyze the massive m7G data generated recently are sorely needed.

RESULTS

We report here m7GHub, a comprehensive online platform for deciphering the location, regulation and pathogenesis of internal mRNA m7G. The m7GHub consists of four main components, including: the first internal mRNA m7G database containing 44 058 experimentally validated internal mRNA m7G sites, a sequence-based high-accuracy predictor, the first web server for assessing the impact of mutations on m7G status, and the first database recording 1218 disease-associated genetic mutations that may function through regulation of m7G methylation. Together, m7GHub will serve as a useful resource for research on internal mRNA m7G modification.

AVAILABILITY AND IMPLEMENTATION

m7GHub is freely accessible online at www.xjtlu.edu.cn/biologicalsciences/m7ghub.

CONTACT

kunqi.chen@liverpool.ac.uk.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Anticodon Wobble Uridine Modification by Elongator at the Crossroad of Cell Signaling, Differentiation, and Diseases

First identified 20 years ago as an RNA polymerase II-associated putative histone acetyltransferase, the conserved Elongator complex has since been recognized as the central player of a complex, regulated, and biologically relevant epitranscriptomic pathway targeting the wobble uridine of some tRNAs. Numerous studies have contributed to three emerging concepts resulting from anticodon modification by Elongator: the codon-specific control of translation, the ability of reprogramming translation in various physiological or pathological contexts, and the maintenance of proteome integrity by counteracting protein aggregation. These three aspects of tRNA modification by Elongator constitute a new layer of regulation that fundamentally contributes to gene expression and are now recognized as being critically involved in various human diseases.

Structure of tRNA methyltransferase complex of Trm7 and Trm734 reveals a novel binding interface for tRNA recognition

The complex between Trm7 and Trm734 (Trm7–Trm734) from Saccharomyces cerevisiae catalyzes 2′-O-methylation at position 34 in tRNA. We report biochemical and structural studies of the Trm7–Trm734 complex. Purified recombinant Trm7–Trm734 preferentially methylates tRNA^Phe transcript variants possessing two of three factors (Cm32, m¹G37 and pyrimidine34). Therefore, tRNA^Phe, tRNA^Trp and tRNA^Leu are specifically methylated by Trm7–Trm734. We have solved the crystal structures of the apo and S-adenosyl-L-methionine bound forms of Trm7–Trm734. Small angle X-ray scattering reveals that Trm7–Trm734 exists as a hetero-dimer in solution. Trm7 possesses a Rossmann-fold catalytic domain, while Trm734 consists of three WD40 β-propeller domains (termed BPA, BPB and BPC). BPA and BPC form a unique V-shaped cleft, which docks to Trm7. The C-terminal region of Trm7 is required for binding to Trm734. The D-arm of substrate tRNA is required for methylation by Trm7–Trm734. If the D-arm in tRNA^Phe is docked onto the positively charged area of BPB in Trm734, the anticodon-loop is located near the catalytic pocket of Trm7. This model suggests that Trm734 is required for correct positioning of tRNA for methylation. Additionally, a point-mutation in Trm7, which is observed in FTSJ1 (human Trm7 ortholog) of nosyndromic X-linked intellectual disability patients, decreases the methylation activity.

tRNA 2'-O-methylation by a duo of TRM7/FTSJ1 proteins modulates small RNA silencing in Drosophila

2′-O-Methylation (Nm) represents one of the most common RNA modifications. Nm affects RNA structure and function with crucial roles in various RNA-mediated processes ranging from RNA silencing, translation, self versus non-self recognition to viral defense mechanisms. Here, we identify two Nm methyltransferases (Nm-MTases) in Drosophila melanogaster (CG7009 and CG5220) as functional orthologs of yeast TRM7 and human FTSJ1. Genetic knockout studies together with MALDI-TOF mass spectrometry and RiboMethSeq mapping revealed that CG7009 is responsible for methylating the wobble position in tRNA^Phe, tRNA^Trp and tRNA^Leu, while CG5220 methylates position C32 in the same tRNAs and also targets additional tRNAs. CG7009 or CG5220 mutant animals were viable and fertile but exhibited various phenotypes such as lifespan reduction, small RNA pathways dysfunction and increased sensitivity to RNA virus infections. Our results provide the first detailed characterization of two TRM7 family members in Drosophila and uncover a molecular link between enzymes catalyzing Nm at specific tRNAs and small RNA-induced gene silencing pathways.

An intellectual disability-associated missense variant in TRMT1 impairs tRNA modification and reconstitution of enzymatic activity

The human TRMT1 gene encodes an RNA methyltransferase enzyme responsible for catalyzing dimethylguanosine (m2,2G) formation in transfer RNAs (tRNAs). Frameshift mutations in TRMT1 have been shown to cause autosomal-recessive intellectual disability (ID) in the human population but additional TRMT1 variants remain to be characterized. Here, we describe a homozygous TRMT1 missense variant in a patient displaying developmental delay, ID, and epilepsy. The missense variant changes an arginine residue to a cysteine (R323C) within the methyltransferase domain and is expected to perturb protein folding. Patient cells expressing TRMT1-R323C exhibit a deficiency in m2,2G modifications within tRNAs, indicating that the mutation causes loss of function. Notably, the TRMT1 R323C mutant retains tRNA binding but is unable to rescue m2,2G formation in TRMT1-deficient human cells. Our results identify a pathogenic point mutation in TRMT1 that perturbs tRNA modification activity and demonstrate that m2,2G modifications are disrupted in the cells of patients with TRMT1-associated ID disorders.

Put the Pedal to the METTL1: Adding Internal m7G Increases mRNA Translation Efficiency and Augments miRNA Processing

Complementary papers by Zhang, Liu, and colleagues (Zhang et al., 2019) and Pandolfini, Barbieri, and colleagues (Pandolfini et al., 2019) develop new sequencing techniques that reveal that METTL1 N7-methylates internal guanosines in mRNAs and miRNAs to increase translation efficiency and miRNA processing, respectively.

Ontogenic mRNA expression of RNA modification writers, erasers, and readers in mouse liver

RNA modifications are recently emerged epigenetic modifications. These diverse RNA modifications have been shown to regulate multiple biological processes, including development. RNA modifications are dynamically controlled by the “writers, erasers, and readers”, where RNA modifying proteins are able to add, remove, and recognize specific chemical modification groups on RNAs. However, little is known about the ontogenic expression of these RNA modifying proteins in various organs, such as liver. In the present study, the hepatic mRNA expression of selected RNA modifying proteins involve in m⁶A, m¹A, m⁵C, hm⁵C, m⁷G, and Ψ modifications was analyzed using the RNA-seq technique. Liver samples were collected from male C57BL/6 mice at several ages from prenatal through neonatal, infant, child to young adult. Results showed that most of the RNA modifying proteins were highly expressed in prenatal mouse liver with a dramatic drop at birth. After birth, most of the RNA modifying proteins showed a downregulation trend during liver maturation. Moreover, the RNA modifying proteins that belong to the same enzyme family were expressed at different abundances at the same ages in mouse liver. In conclusion, this study unveils that the mRNA expression of RNA modifying proteins follows specific ontogenic expression patterns in mice liver during maturation. These data indicated that the changes in expression of RNA modifying proteins might have a potential role to regulate gene expression in liver through alteration of RNA modification status.

Small noncoding RNA discovery and profiling with sRNAtools based on high-throughput sequencing

Small noncoding RNAs (sRNA/sncRNAs) are generated from different genomic loci and play important roles in biological processes, such as cell proliferation and the regulation of gene expression. Next-generation sequencing (NGS) has provided an unprecedented opportunity to discover and quantify diverse kinds of sncRNA, such as tRFs (tRNA-derived small RNA fragments), phasiRNAs (phased, secondary, small-interfering RNAs), Piwi-interacting RNA (piRNAs) and plant-specific 24-nt short interfering RNAs (siRNAs). However, currently available web-based tools do not provide approaches to comprehensively analyze all of these diverse sncRNAs. This study presents a novel integrated platform, sRNAtools (https://bioinformatics.caf.ac.cn/sRNAtools), that can be used in conjunction with high-throughput sequencing to identify and functionally annotate sncRNAs, including profiling microRNAss, piRNAs, tRNAs, small nuclear RNAs, small nucleolar RNAs and rRNAs and discovering isomiRs, tRFs, phasiRNAs and plant-specific 24-nt siRNAs for up to 21 model organisms. Different modules, including single case, batch case, group case and target case, are developed to provide users with flexible ways of studying sncRNA. In addition, sRNAtools supports different ways of uploading small RNA sequencing data in a very interactive queue system, while local versions based on the program package/Docker/virtureBox are also available. We believe that sRNAtools will greatly benefit the scientific community as an integrated tool for studying sncRNAs.

TrmB, a tRNA m7G46 methyltransferase, plays a role in hydrogen peroxide resistance and positively modulates the translation of katA and katB mRNAs in Pseudomonas aeruginosa

Cellular response to oxidative stress is a crucial mechanism that promotes the survival of Pseudomonas aeruginosa during infection. However, the translational regulation of oxidative stress response remains largely unknown. Here, we reveal a tRNA modification-mediated translational response to H₂O₂ in P. aeruginosa. We demonstrated that the P. aeruginosa trmB gene encodes a tRNA guanine (46)-N7-methyltransferase that catalyzes the formation of m⁷G46 in the tRNA variable loop. Twenty-three tRNA substrates of TrmB with a guanosine residue at position 46 were identified, including 11 novel tRNA substrates. We showed that loss of trmB had a strong negative effect on the translation of Phe- and Asp-enriched mRNAs. The trmB-mediated m⁷G modification modulated the expression of the catalase genes katA and katB, which are enriched with Phe/Asp codons at the translational level. In response to H₂O₂ exposure, the level of m⁷G modification increased, consistent with the increased translation efficiency of Phe- and Asp-enriched mRNAs. Inactivation of trmB led to decreased KatA and KatB protein abundance and decreased catalase activity, resulting in H₂O₂-sensitive phenotype. Taken together, our observations reveal a novel role of m⁷G46 tRNA modification in oxidative stress response through translational regulation of Phe- and Asp-enriched genes, such as katA and katB.

Nucleotide resolution profiling of m7G tRNA modification by TRAC-Seq

Precise identification of sites of RNA modification is key to studying the functional role of such modifications in the regulation of gene expression and for elucidating relevance to diverse physiological processes. tRNA reduction and cleavage sequencing (TRAC-Seq) is a chemically based approach for the unbiased global mapping of 7-methylguansine (m⁷G) modification of tRNAs at single-nucleotide resolution throughout the tRNA transcriptome. m⁷G TRAC-Seq involves the treatment of size-selected (<200 nt) RNAs with the demethylase AlkB to remove major tRNA modifications, followed by sodium borohydride (NaBH₄) reduction of m⁷G sites and subsequent aniline-mediated cleavage of the RNA chain at the resulting abasic sites. The cleaved sites are subsequently ligated with adaptors for the construction of libraries for high-throughput sequencing. The m⁷G modification sites are identified using a bioinformatic pipeline that calculates the cleavage scores at individual sites on all tRNAs. Unlike antibody-based methods, such as methylated RNA immunoprecipitation and sequencing (meRIP-Seq) for enrichment of methylated RNA sequences, chemically based approaches, including TRAC-Seq, can provide nucleotide-level resolution of modification sites. Compared to the related method AlkAniline-Seq (alkaline hydrolysis and aniline cleavage sequencing), TRAC-Seq incorporates small RNA selection, AlkB demethylation, and sodium borohydride reduction steps to achieve specific and efficient single-nucleotide resolution profiling of m⁷G sites in tRNAs. The m⁷G TRAC-Seq protocol could be adapted to chemical cleavage-mediated detection of other RNA modifications. The protocol can be completed within ~9 d for four biological replicates of input and treated samples.