RBM45 is an m6A-binding protein that affects neuronal differentiation and the splicing of a subset of mRNAs

Single-cell m6A profiling in the mouse brain uncovers cell type-specific RNA methylomes and age-dependent differential methylation

N6-methyladenosine (m6A) is an abundant mRNA modification in the brain that has important roles in neurodevelopment and brain function. However, because of technical limitations, global profiling of m6A sites within the individual cell types that make up the brain has not been possible. Here, we develop a mouse model that enables transcriptome-wide m6A detection in any tissue of interest at single-cell resolution. We use these mice to map m6A across different brain regions and within single cells of the mouse cortex and discover a high degree of shared methylation across brain regions and cell types. However, we also identify a small number of differentially methylated mRNAs in neurons that encode important regulators of neuronal signaling, and we discover that microglia have lower levels of m6A than other cell types. Finally, we perform single-cell m6A mapping in aged mice and identify many transcripts with age-dependent changes in m6A.

F11R RNA trinucleotide over-edited by ADAR in gastric and colorectal cancers: Cross-cohort validation, gene expression regulation, and diagnostic significance

The F11 receptor (F11R) gene encoding junctional adhesion molecule A has been associated with gastric cancer (GC) and colorectal cancer (CRC), in which its role and regulation remain to be further elucidated. Recently F11R was also identified as a potential target of adenosine-to-inosine (A-to-I) mediated by the adenosine deaminases acting on RNA (ADARs). Herein, using RNA-Seq and experimental validation, our current study revealed an F11R RNA trinucleotide over-edited by ADAR, with its regulation of gene expression and clinical significance in four GC and three CRC cohorts. Our results found an over-edited AAA trinucleotide in an AluSg located in the F11R 3'-untranslated region (3'-UTR), which showed editing levels correlated with elevated ADAR expression across all GC and CRC cohorts in our study. Overexpression and knockdown of ADAR in GC and CRC cells, followed by RNA-Seq and Sanger sequencing, confirmed the ADAR-mediated F11R 3'-UTR trinucleotide editing, which potentially disrupted an RBM45 binding site identified by crosslinking immunoprecipitation sequencing (CLIP-seq) and regulated F11R expression in luciferase reporter assays. Moreover, the F11R trinucleotide editing showed promising predictive performance for diagnosing GC and CRC across GC and CRC cohorts. Our findings thus highlight both the potential biological and clinical significance of an ADAR-edited F11R trinucleotide in GC and CRC, providing new insights into its application as a novel diagnostic biomarker for both cancers.

Structural basis for RNA recognition by the C-terminal RRM domain of human RBM45

RBM45 is an RNA-binding protein with roles in neural development by regulating RNA splicing. Its dysfunction and aggregation are associated with neurodegenerative diseases such as amyotrophic lateral sclerosis (ALS) and frontotemporal lobar dementia (FTLD). RBM45 harbors three RRM domains that potentially bind RNA. While the recognitions of RNA by its N-terminal tandem RRM domains (RRM1 and RRM2) have been well understood, the RNA-binding property of its C-terminal RRM (RRM3) remains unclear. In this work, we identified that the RRM3 of the RBM45 sequence specifically binds RNA with a GACG sequence, similar but not identical to those recognized by the RRM1 and RRM2. Further, we determined the crystal structure of RBM45RRM3 in complex with a GACG sequence-containing single-stranded DNA. Our structural results, together with the RNA-binding assays of mutants at key amino acid residues, revealed the molecular mechanism by which RBM45RRM3 recognizes an RNA sequence. Our finding on the RNA-binding property of the individual RRM module of RBM45 provides the foundation for unraveling the RNA-binding characteristics of full-length RBM45 and for understanding the biological functions of RBM45.

Loss of RBM45 inhibits breast cancer progression by reducing the SUMOylation of IRF7 to promote IFNB1 transcription

Type I interferons exhibit anti-proliferative and anti-cancer activities, but their detailed regulatory mechanisms in cancer have not been fully elucidated yet. RNA binding proteins are master orchestrators of gene regulation, which are closely related to tumor progression. Here we show that the upregulated RNA binding protein RBM45 correlates with poor prognosis in breast cancer. Depletion of RBM45 suppresses breast cancer progression both in cultured cells and xenograft mouse models. Mechanistically, RBM45 ablation inhibits breast cancer progression through regulating type I interferon signaling, particularly by elevating IFN-β production. Importantly, RBM45 recruits TRIM28 to IRF7 and stimulates its SUMOylation, thereby repressing IFNB1 transcription. Loss of RBM45 reduced the SUMOylation of IRF7 by reducing the interaction between TRIM28 and IRF7 to promote IFNB1 transcription, leading to the inhibition of breast cancer progression. Taken together, our finding uncovers a vital role of RBM45 in modulating type I interferon signaling and cancer aggressive progression, implicating RBM45 as a potential therapeutic target in breast cancer.

Mutations causing premature termination codons discriminate and generate cellular and clinical variability in HHT

ABSTRACT

For monogenic diseases caused by pathogenic loss-of-function DNA variants, attention focuses on dysregulated gene-specific pathways, usually considering molecular subtypes together within causal genes. To better understand phenotypic variability in hereditary hemorrhagic telangiectasia (HHT), we subcategorized pathogenic DNA variants in ENG/endoglin, ACVRL1/ALK1, and SMAD4 if they generated premature termination codons (PTCs) subject to nonsense-mediated decay. In 3 patient cohorts, a PTC-based classification system explained some previously puzzling hemorrhage variability. In blood outgrowth endothelial cells (BOECs) derived from patients with ACVRL1+/PTC, ENG+/PTC, and SMAD4+/PTC genotypes, PTC-containing RNA transcripts persisted at low levels (8%-23% expected, varying between replicate cultures); genes differentially expressed to Bonferroni P < .05 in HHT+/PTC BOECs clustered significantly only to generic protein terms (isopeptide-bond/ubiquitin-like conjugation) and pulse-chase experiments detected subtle protein maturation differences but no evidence for PTC-truncated protein. BOECs displaying highest PTC persistence were discriminated in unsupervised hierarchical clustering of near-invariant housekeeper genes, with patterns compatible with higher cellular stress in BOECs with >11% PTC persistence. To test directionality, we used a HeLa reporter system to detect induction of activating transcription factor 4 (ATF4), which controls expression of stress-adaptive genes, and showed that ENG Q436X but not ENG R93X directly induced ATF4. AlphaFold accurately modeled relevant ENG domains, with AlphaMissense suggesting that readthrough substitutions would be benign for ENG R93X and other less rare ENG nonsense variants but more damaging for Q436X. We conclude that PTCs should be distinguished from other loss-of-function variants, PTC transcript levels increase in stressed cells, and readthrough proteins and mechanisms provide promising research avenues.

RNA m6A modification, signals for degradation or stabilisation?

The RNA modification N6-methyladenosine (m6A) is conserved across eukaryotes, and profoundly influences RNA metabolism, including regulating RNA stability. METTL3 and METTL14, together with several accessory components, form a 'writer' complex catalysing m6A modification. Conversely, FTO and ALKBH5 function as demethylases, rendering m6A dynamic. Key to understanding the functional significance of m6A is its 'reader' proteins, exemplified by YTH-domain-containing proteins (YTHDFs) canonical reader and insulin-like growth factor 2 mRNA-binding proteins (IGF2BPs) non-canonical reader. These proteins play a crucial role in determining RNA stability: YTHDFs mainly promote mRNA degradation through different cytoplasmic pathways, whereas IGF2BPs function to maintain mRNA stability. Additionally, YTHDC1 functions within the nucleus to degrade or protect certain m6A-containing RNAs, and other non-canonical readers also contribute to RNA stability regulation. Notably, m6A regulates retrotransposon LINE1 RNA stability and/or transcription via multiple mechanisms. However, conflicting observations underscore the complexities underlying m6A's regulation of RNA stability depending upon the RNA sequence/structure context, developmental stage, and/or cellular environment. Understanding the interplay between m6A and other RNA regulatory elements is pivotal in deciphering the multifaceted roles m6A plays in RNA stability regulation and broader cellular biology.

Studying m6A in the brain: a perspective on current methods, challenges, and future directions

A major mechanism of post-transcriptional RNA regulation in cells is the addition of chemical modifications to RNA nucleosides, which contributes to nearly every aspect of the RNA life cycle. N6-methyladenosine (m6A) is a highly prevalent modification in cellular mRNAs and non-coding RNAs, and it plays important roles in the control of gene expression and cellular function. Within the brain, proper regulation of m6A is critical for neurodevelopment, learning and memory, and the response to injury, and m6A dysregulation has been implicated in a variety of neurological disorders. Thus, understanding m6A and how it is regulated in the brain is important for uncovering its roles in brain function and potentially identifying novel therapeutic pathways for human disease. Much of our knowledge of m6A has been driven by technical advances in the ability to map and quantify m6A sites. Here, we review current technologies for characterizing m6A and highlight emerging methods. We discuss the advantages and limitations of current tools as well as major challenges going forward, and we provide our perspective on how continued developments in this area can propel our understanding of m6A in the brain and its role in brain disease.

HNRNPA2B1 and HNRNPR stabilize ASCL1 in an m6A-dependent manner to promote neuroblastoma progression

HNRNPA2B1 and HNRNPR stabilize ASCL1 mRNA in neuroblastoma, but whether their regulatory effects depend on m6A modification and whether their function involves ASCL1 remain unknown. This study investigated the m6A-dependent binding of HNRNPA2B1 and HNRNPR to ASCL1 and subsequent regulation, as well as the expression, clinical significance, and function of HNRNPA2B1 and HNRNPR in neuroblastoma. We revealed that METTL14 mediated ASCL1 m6A modification to stabilize ASCL1. HNRNPA2B1 and HNRNPR significantly enriched ASCL1 mRNA by binding to the 5' and 3' untranslated regions, respectively, and METTL14 knockdown reduced this enrichment. Mutations in m6A sites in the untranslated regions of ASCL1 mRNA considerably decreased probe capacity to engage HNRNPA2B1 and HNRNPR. HNRNPR interacts with IGF2BP1, and knocking down either impaired binding to ASCL1 mRNA. HNRNPA2B1 and HNRNPR knockdown suppressed neuroblastoma cell growth and invasion, while ASCL1 overexpression restored these effects. The high HNRNPA2B1 and HNRNPR expression in neuroblastoma correlated with ASCL1 expression. Thus, HNRNPA2B1 and HNRNPR bind and stabilize ASCL1 mRNA in an m6A-dependent manner to promote neuroblastoma progression. This study not only discovered a new mechanism underlying the high ASCL1 expression in neuroblastoma but also identified the HNRNPA2B1/HNRNPR/ASCL1 axis as a promising target for inhibiting neuroblastoma progression.

Modifying the antiviral innate immune response by selective writing, erasing, and reading of m6A on viral and cellular RNA

Viral infection and the antiviral innate immune response are regulated by the RNA modification m6A. m6A directs nearly all aspects of RNA metabolism by recruiting RNA-binding proteins that mediate the fate of m6A-containing RNA. m6A controls the antiviral innate immune response in diverse ways, including shielding viral RNA from detection by antiviral sensors and influencing the expression of cellular mRNAs encoding antiviral signaling proteins, cytokines, and effector proteins. While m6A and the m6A machinery are important for the antiviral response, the precise mechanisms that determine how the m6A machinery selects specific viral or cellular RNA molecules for modification during infection are not fully understood. In this review, we highlight recent findings that shed light on how viral infection redirects the m6A machinery during the antiviral response. A better understanding of m6A targeting during viral infection could lead to new immunomodulatory and therapeutic strategies against viral infection.

RBM45 reprograms lipid metabolism promoting hepatocellular carcinoma via Rictor and ACSL1/ACSL4

Reprogramming of lipid metabolism during hepatocarcinogenesis is not well elucidated. Here, we aimed to explore pivotal RNA-binding motif proteins (RBMs) in lipid metabolism and their therapeutic potential in hepatocellular carcinoma (HCC). Through bioinformatic analysis, we identified RBM45 as a critical gene of interest among differentially expressed RBMs in HCC, with significant prognostic relevance. RBM45 influenced the malignant biological phenotype and lipid metabolism of HCC cells. Mechanically, RBM45 promotes de novo lipogenesis in HCC by directly targeting two key enzymes involved in long-chain fatty acid synthesis, ACSL1 and ACSL4. RBM45 also targets Rictor, which has been demonstrated to modulate lipid metabolism profoundly. RBM45 also aided lipid degradation through activating a key fatty acid β oxidation enzyme, CPT1A. Thus, RBM45 boosted lipid synthesis and decomposition, indicating an enhanced utility of lipid fuels in HCC. Clinically, body mass index was positively correlated with RBM45 in human HCCs. The combination of a PI3K/AKT/mTOR pathway inhibitor in vitro or Sorafenib in orthotopic liver cancer mouse models with shRBM45 has a more significant therapeutic effect on liver cancer than the drug alone. In summary, our findings highlight the versatile roles of RBM45 in lipid metabolism reprogramming and its therapeutic potential in HCC. Lipids induced RBM45 expression. In turn, RBM45 promoted the utility of lipid in HCCs through accelerating both de novo lipogenesis and fatty acid β oxidation, which required the participation of Rictor, a core component of mTORC2 that has been demonstrated to modulate lipid metabolism potently, as well as ACSL1/ACSL4, two key enzymes of long-chain fatty acid synthesis. When the first-line chemotherapy drug sorafenib is combined with a PI3K/AKT/mTOR pathway inhibitor (MK2206 is an AKT inhibitor, rapamycin is a mTOR inhibitor, and inhibiting RBM45 can significantly inhibit Rictor), cell cycle, proliferation, lipid metabolism reprogramming, and hepatocarcinogenesis can be significantly inhibited, while apoptosis can be significantly enhanced.

Captive ERVWE1 triggers impairment of 5-HT neuronal plasticity in the first-episode schizophrenia by post-transcriptional activation of HTR1B in ALKBH5-m6A dependent epigenetic mechanisms

BACKGROUND

Abnormalities in the 5-HT system and synaptic plasticity are hallmark features of schizophrenia. Previous studies suggest that the human endogenous retrovirus W family envelope (ERVWE1) is an influential risk factor for schizophrenia and inversely correlates with 5-HT4 receptor in schizophrenia. To our knowledge, no data describes the effect of ERVWE1 on 5-HT neuronal plasticity. N6-methyladenosine (m6A) regulates gene expression and impacts synaptic plasticity. Our research aims to systematically investigate the effects of ERVWE1 on 5-HT neuronal plasticity through m6A modification in schizophrenia.

RESULTS

HTR1B, ALKBH5, and Arc exhibited higher levels in individuals with first-episode schizophrenia compared to the controls and showed a strong positive correlation with ERVWE1. Interestingly, HTR1B was also correlated with ALKBH5 and Arc. Further analyses confirmed that ALKBH5 may be an independent risk factor for schizophrenia. In vitro studies, we discovered that ERVWE1 enhanced HTR1B expression, thereby activating the ERK-ELK1-Arc pathway and reducing the complexity and spine density of 5-HT neurons. Furthermore, ERVWE1 reduced m6A levels through ALKBH5 demethylation. ERVWE1 induced HTR1B upregulation by improving its mRNA stability in ALKBH5-m6A-dependent epigenetic mechanisms. Importantly, ALKBH5 mediated the observed alterations in 5-HT neuronal plasticity induced by ERVWE1.

CONCLUSIONS

Overall, HTR1B, Arc, and ALKBH5 levels were increased in schizophrenia and positively associated with ERVWE1. Moreover, ALKBH5 was a novel risk gene for schizophrenia. ERVWE1 impaired 5-HT neuronal plasticity in ALKBH5-m6A dependent mechanism by the HTR1B-ERK-ELK1-Arc pathway, which may be an important contributor to aberrant synaptic plasticity in schizophrenia.

Climate Adaptation and Drift Shape the Genomes of Two Eel-Goby Sister Species Endemic to Contrasting Latitude

Deciphering the role of climate adaptation in generating genetic divergence and hence speciation is a central question in evolution. Comparisons of genomes of closely related species spanning selective climate gradients are particularly informative in discerning the signatures of selection and thereby providing valuable information concerning the role of climate adaptation in speciation. Here we re-sequenced 99 genomes of the two sister eel-goby species Odontamblyopus lacepedii and O. rebecca, which are endemic to tidal mudflats spanning contrasting latitude gradients, to estimate the influence of divergent climate selection on shaping genome-wide patterns of divergence. The results indicated that genome-wide differentiation between the two species was evident (genome-wide FST = 0.313). Against a background of high baseline genomic divergence, 588 and 1202 elevated divergent loci were detected to be widespread throughout their genomes, as opposed to focused within small islands of genomic regions. These patterns of divergence may arise from divergent climate selection in addition to genetic drift acting through past glacial segregation (1.46 million years ago). We identified several candidate genes that exhibited elevated divergence between the two species, including genes associated with substance metabolism, energy production, and response to environmental cues, all putative candidates closely linked to thermal adaptation expected from the latitude gradient. Interestingly, several candidates related to gamete recognition and time of puberty, and also exhibited elevated divergence, indicating their possible role in pre-zygote isolation and speciation of the two species. Our results would expand our knowledge on the roles of latitude climate adaptation and genetic drift in generating and maintaining biodiversity in marine teleosts.

RNA binding motif protein 45-mediated phosphorylation enhances protein stability of ASCT2 to promote hepatocellular carcinoma progression

Targeting metabolic remodeling represents a potentially promising strategy for hepatocellular carcinoma (HCC) therapy. In-depth understanding on the regulation of the glutamine transporter alanine-serine-cysteine transporter 2 (ASCT2) contributes to the development of novel promising therapeutics. As a developmentally regulated RNA binding protein, RBM45 is capable to shuttle between nucleus and cytoplasm, and directly interacts with proteins. By bioinformatics analysis, we screened out that RBM45 was elevated in the HCC patient specimens and positively correlated with poor prognosis. RBM45 promoted cell proliferation, boosted xenograft tumorigenicity and accelerated HCC progression. Using untargeted metabolomics, it was found that RBM45 interfered with glutamine metabolism. Further results demonstrated that RBM45 positively associated with ASCT2 in human and mouse specimens. Moreover, RBM45 enhanced ASCT2 protein stability by counteracting autophagy-independent lysosomal degradation. Significantly, wild-type ASCT2, instead of phospho-defective mutants, rescued siRBM45-suppressed HCC cell proliferation. Using molecular docking approaches, we found AG-221, a mutant isocitrate dehydrogenase 2 (mIDH2) inhibitor for acute myeloid leukemia therapy, pharmacologically perturbed RBM45-ASCT2 interaction, decreased ASCT2 stability and suppressed HCC progression. These findings provide evidence that RBM45 plays a crucial role in HCC progression via interacting with and counteracting the degradation of ASCT2. Our findings suggest a novel alternative structural sites for the design of ASCT2 inhibitors and the agents interfering with RBM45-ASCT2 interaction may be a potential direction for HCC drug development.

Emerging Roles for DNA 6mA and RNA m6A Methylation in Mammalian Genome

Epigenetic methylation has been shown to play an important role in transcriptional regulation and disease pathogenesis. Recent advancements in detection techniques have identified DNA N6-methyldeoxyadenosine (6mA) and RNA N6-methyladenosine (m6A) as methylation modifications at the sixth position of adenine in DNA and RNA, respectively. While the distributions and functions of 6mA and m6A have been extensively studied in prokaryotes, their roles in the mammalian brain, where they are enriched, are still not fully understood. In this review, we provide a comprehensive summary of the current research progress on 6mA and m6A, as well as their associated writers, erasers, and readers at both DNA and RNA levels. Specifically, we focus on the potential roles of 6mA and m6A in the fundamental biological pathways of the mammalian genome and highlight the significant regulatory functions of 6mA in neurodegenerative diseases.

Transcriptome analysis identification of A-to-I RNA editing in granulosa cells associated with PCOS

Background

Polycystic ovary syndrome (PCOS) is a complex, multifactor disorder in women of reproductive age worldwide. Although RNA editing may contribute to a variety of diseases, its role in PCOS remains unclear.

Methods

A discovery RNA-Seq dataset was obtained from the NCBI Gene Expression Omnibus database of granulosa cells from women with PCOS and women without PCOS (controls). A validation RNA-Seq dataset downloaded from the European Nucleotide Archive Databank was used to validate differential editing. Transcriptome-wide investigation was conducted to analyze adenosine-to-inosine (A-to-I) RNA editing in PCOS and control samples.

Results

A total of 17,395 high-confidence A-to-I RNA editing sites were identified in 3,644 genes in all GC samples. As for differential RNA editing, there were 545 differential RNA editing (DRE) sites in 259 genes with Nucleoporin 43 (NUP43), Retinoblastoma Binding Protein 4 (RBBP4), and leckstrin homology-like domain family A member 1 (PHLDA) showing the most significant three 3'-untranslated region (3'UTR) editing. Furthermore, we identified 20 DRE sites that demonstrated a significant correlation between editing levels and gene expression levels. Notably, MIR193b-365a Host Gene (MIR193BHG) and Hook Microtubule Tethering Protein 3 (HOOK3) exhibited significant differential expression between PCOS and controls. Functional enrichment analysis showed that these 259 differentially edited genes were mainly related to apoptosis and necroptosis pathways. RNA binding protein (RBP) analysis revealed that RNA Binding Motif Protein 45 (RBM45) was predicted as the most frequent RBP binding with RNA editing sites. Additionally, we observed a correlation between editing levels of differential editing sites and the expression level of the RNA editing enzyme Adenosine Deaminase RNA Specific B1 (ADARB1). Moreover, the existence of 55 common differentially edited genes and nine differential editing sites were confirmed in the validation dataset.

Conclusion

Our current study highlighted the potential role of RNA editing in the pathophysiology of PCOS as an epigenetic process. These findings could provide valuable insights into the development of more targeted and effective treatment options for PCOS.

The Proteins of mRNA Modification: Writers, Readers, and Erasers

Over the past decade, mRNA modifications have emerged as important regulators of gene expression control in cells. Fueled in large part by the development of tools for detecting RNA modifications transcriptome wide, researchers have uncovered a diverse epitranscriptome that serves as an additional layer of gene regulation beyond simple RNA sequence. Here, we review the proteins that write, read, and erase these marks, with a particular focus on the most abundant internal modification, N6-methyladenosine (m6A). We first describe the discovery of the key enzymes that deposit and remove m6A and other modifications and discuss how our understanding of these proteins has shaped our views of modification dynamics. We then review current models for the function of m6A reader proteins and how our knowledge of these proteins has evolved. Finally, we highlight important future directions for the field and discuss key questions that remain unanswered.

N6-methyladenosine-mediated gene regulation and therapeutic implications

N6-methyladenosine (m⁶A) RNA methylation is the most abundant form of mRNA modification in eukaryotes and is at the front line of biological and biomedical research. This dynamic and reversible m⁶A RNA modification determines the fates of modified RNA molecules at the post-transcriptional level, affecting almost all important biological processes. Notably, m⁶A is also involved in chromatin and transcriptional regulation, while m⁶A dysregulation is implicated in various diseases. Here, we review current knowledge of post-transcriptional and transcriptional regulatory mechanisms involving m⁶A modification. We also discuss their involvement in the occurrence and development of diseases, including cancer, as well as potential theranostic targets, in hope of facilitating the translation of preclinical findings to the clinic.

Integrated analysis of the voltage-gated potassium channel-associated gene KCNH2 across cancers

KCNH2 encodes the human ether-a-go-go-related gene (hERG) potassium channel and is an important repolarization reserve for regulating cardiac electrical activity. Increasing evidence suggests that it is involved in the development of various tumours, yet a thorough analysis of the underlying process has not been performed. Here, we have comprehensively examined the role of KCNH2 in multiple cancers by assessing KCNH2 gene expression, diagnostic and prognostic value, genetic alterations, immune infiltration correlations, RNA modifications, mutations, clinical correlations, interacting proteins, and associated signalling pathways. KCNH2 is differentially expressed in over 30 cancers and has a high diagnostic value for 10 tumours. Survival analysis showed that high expression of KCNH2 was associated with a poor prognosis in glioblastoma multiforme (GBM) and hepatocellular carcinoma (LIHC). Mutations and RNA methylation modifications (especially m6A) of KCNH2 are associated with its expression in multiple tumours. KCNH2 expression is correlated with tumour mutation burden, microsatellite instability, neoantigen load, and mutant-allele tumour heterogeneity. In addition, KCNH2 expression is associated with the tumour immune microenvironment and its immunosuppressive phenotype. KEGG signalling pathway enrichment analysis revealed that KCNH2 and its interacting molecules are involved in a variety of pathways related to carcinogenesis and signal regulation, such as the PI3K/Akt and focal adhesion pathways. Overall, we found that KCNH2 and its interaction molecular are expected to be immune-related biomarkers for cancer diagnosis and prognosis evaluation, and are potential regulatory targets of singalling pathways for tumour development due to their significant role in cancers.