Hostname: page-component-78c5997874-xbtfd Total loading time: 0 Render date: 2024-11-10T07:44:19.322Z Has data issue: false hasContentIssue false

RNA m6A methyltransferase activator affects anxiety-related behaviours, monoamines and striatal gene expression in the rat

Published online by Cambridge University Press:  09 October 2024

Margus Kanarik
Affiliation:
Division of Neuropsychopharmacology, Institute of Chemistry, University of Tartu, Tartu, Tartumaa, Estonia
Kristi Liiver
Affiliation:
School of Natural Sciences and Health, Tallinn University, Tallinn, Harjumaa, Estonia
Marianna Norden
Affiliation:
School of Natural Sciences and Health, Tallinn University, Tallinn, Harjumaa, Estonia Institute of Chemistry and Biotechnology, Tallinn University of Technology, Tallinn, Estonia
Indrek Teino
Affiliation:
Institute of Biotechnology, HiLIFE, University of Helsinki, Helsinki, Uusimaa, Finland Institute of Molecular and Cell Biology, University of Tartu, Tartu, Tartumaa, Estonia
Tõnis Org
Affiliation:
Institute of Genomics, University of Tartu, Tartu, Tartumaa, Estonia
Karita Laugus
Affiliation:
Division of Neuropsychopharmacology, Institute of Chemistry, University of Tartu, Tartu, Tartumaa, Estonia
Ruth Shimmo
Affiliation:
School of Natural Sciences and Health, Tallinn University, Tallinn, Harjumaa, Estonia
Mati Karelson
Affiliation:
Division of Molecular Technology, Institute of Chemistry, University of Tartu, Tartu, Estonia
Mart Saarma
Affiliation:
Institute of Biotechnology, HiLIFE, University of Helsinki, Helsinki, Uusimaa, Finland
Jaanus Harro*
Affiliation:
Division of Neuropsychopharmacology, Institute of Chemistry, University of Tartu, Tartu, Tartumaa, Estonia Institute of Biotechnology, HiLIFE, University of Helsinki, Helsinki, Uusimaa, Finland
*
Corresponding author: Jaanus Harro; Email: jaanus.harro@ut.ee
Rights & Permissions [Opens in a new window]

Abstract

Modification of mRNA by methylation is involved in post-transcriptional regulation of gene expression by affecting the splicing, transport, stability and translation of mRNA. Methylation of adenosine at N6 (m6A) is one of the most common and important cellular modification occurring in the mRNA of eukaryotes. Evidence that m6A mRNA methylation is involved in regulation of stress response and that its dysregulation may contribute to the pathogenesis of neuropsychiatric disorders is accumulating. We have examined the acute and subchronic (up to 18 days once per day intraperitoneally) effect of the first METTL3/METTL14 activator compound CHMA1004 (methyl-piperazine-2-carboxylate) at two doses (1 and 5 mg/kg) in male and female rats. CHMA1004 had a locomotor activating and anxiolytic-like profile in open field and elevated zero-maze tests. In female rats sucrose consumption and swimming in Porsolt’s test were increased. Nevertheless, CHMA1004 did not exhibit strong psychostimulant-like properties: CHMA1004 had no effect on 50-kHz ultrasonic vocalizations except that it reduced the baseline difference between male and female animals, and acute drug treatment had no effect on extracellular dopamine levels in striatum. Subchronic CHMA1004 altered ex vivo catecholamine levels in several brain regions. RNA sequencing of female rat striata after subchronic CHMA1004 treatment revealed changes in the expression of a number of genes linked to dopamine neuron viability, neurodegeneration, depression, anxiety and stress response. Conclusively, the first-in-class METTL3/METTL14 activator compound CHMA1004 increased locomotor activity and elicited anxiolytic-like effects after systemic administration, demonstrating that pharmacological activation of RNA m6A methylation has potential for neuropsychiatric drug development.

Type
Original Article
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.
Copyright
© The Author(s), 2024. Published by Cambridge University Press on behalf of Scandinavian College of Neuropsychopharmacology

Significant outcomes

  • Acute and subchronic systemic administration of the first-in-class METTL3/METTL14 activator drug CHMA1004 (methyl-piperazine-2-carboxylate) was found to exhibit activating and anxiotylic-like properties.

  • Behavioural activation by increased m6A mRNA methylation was not accompanied by behavioural or neurochemical effects typical for addictive psychostimulants.

  • Action on m6A mRNA methylation affected catecholamine levels in several brain regions and gene expression related to dopamine neuron viability and stress response.

Limitations

  • The number of male and female animals separately in these experiments was small.

  • Other types of behavioural tests will be required to fully characterise the psychopharmacological profile of METTL3/METTL14 activators.

  • While specific alterations in gene expression were demonstrated, it was not established whether these were directly owing to increased RNA methylation or downstream effects.

Introduction

Psychiatric disorders mostly develop in gene-environment interactions and have a chronic and relapsing course, likely mediated by multiple epigenetic and downstream regulatory mechanisms. Epigenetic DNA modifications, especially DNA methylation, have garnered much scrutiny in recent years (Juruena et al., Reference Juruena, Gadelrab, Cleare and Young2021; Panariello et al., Reference Panariello, Fanelli, Fabbri, Atti, De Ronchi and Serretti2022; Harvanek et al., Reference Harvanek, Boks, Vinkers and Higgins-Chen2023; Persaud and Cates, Reference Persaud and Cates2023; Pishva et al., Reference Pishva, van den Hove, Laroche, Lvovs, Roy, Ortega, Burrage, Veidebaum, Kanarik, Mill, Lesch and Harro2023). However, limited attention has been paid to modification of mRNA by methylation. Methylation of adenosine at N6 (m6A) is the most common and most important cellular modification occurring in the mRNA of eukaryotes. m6A methylation regulates splicing, transport, stability and translation of mRNA (Zaccara et al., Reference Zaccara, Ries and Jaffrey2019). m6A formation is catalysed by a methyltransferase complex containing methyltransferase-like 3 (METTL3), methyltransferase-like 14 (METTL14) and Wilms-tumor-1-associated protein (WTAP) and by a single subunit enzyme methyltransferase-like 16 (METTL16). RNA m6A demethylation is catalysed by the fat mass and obesity-associated protein FTO and the RNA demethylase enzyme ALKBH5. These enzymes modify usually only one but sometimes up to three adenosine residues in the mRNA molecules generally located in the 3´untranslated region. RNA sequences containing m6A residues are recognised by a dozen of RNA binding proteins, for example YTHDC1, YTHDC2, YTHDF1, YTHDF2 and YTHDF3 that regulate the fate of mRNA (Meyer and Jaffrey, Reference Meyer and Jaffrey2017). m6A modification occurs in all organs but is particularly abundant in mammalian brain. There is increasing evidence on how mRNA m6A methylation is involved in regulation of behaviour, suggestive that its dysregulation may contribute to neuropsychiatric disorders (Zaccara et al., Reference Zaccara, Ries and Jaffrey2019; Livneh et al., Reference Livneh, Moshitch-Moshkovitz, Amariglio, Rechavi and Dominissini2020; Selberg et al., Reference Selberg, Žusinaite, Herodes, Seli, Kankuri, Merits and Karelson2021b; Shafik et al., Reference Shafik, Allen and Jin2022; Su et al., Reference Su, Qu and Mu2023).

A pioneering study demonstrated that interference with RNA m6A methylation affects clock gene transcription and by this means the length of circadian period in mice (Fustin et al., Reference Fustin, Doi, Yamaguchi, Hida, Nishimura, Yoshida, Isagawa, Morioka, Kakeya, Manabe and Okamura2013). More recently, artificial lighting was shown to affect RNA m6A methylation and corresponding gene expression in the rat lateral habenula (Li Y et al., Reference Li, Ren, Zhang, Weng, Zhang, Zou, Wu and Hu2023), a brain region pivotal for fronto-limbic to brainstem connections and implicated in affective and substance use disorders (Ables et al., Reference Ables, Park and Ibañez-Tallon2023). Recent studies have suggested the involvement of cerebral RNA m6A methylation in stress response, whereas m6A and its regulatory gene expression levels can change in a regionally specific manner (Engel et al., Reference Engel, Eggert, Kaplick, Eder, Röh, Tietze, Namendorf, Arloth, Weber, Rex-Haffner, Geula, Jakovcevski, Hanna, Leshkowitz, Uhr, Wotjak, Schmidt, Deussing, Binder and Chen2018). Several studies have found changes in RNA m6A methylation after fear conditioning and effects of manipulation on m6A methylation enzymes on fear-conditioned behaviour (Widagdo et al., Reference Widagdo, Zhao, Kempen, Tan, Ratnu, Wei, Leighton, Spadaro, Edson, Anggono and Bredy2016; Walters et al., Reference Walters, Mercaldo, Gillon, Yip, Neve, Boyce, Frankland and Josselyn2017; Chang et al., Reference Chang, Zhang, Zhang, Zhu, Zheng, Qian, Wei, Liu, Liu, Ren and Han2022). In learned helpless rats, the expression of several plasticity genes in hippocampus was downregulated along with enrichment in their m6A methylation (Roy et al., Reference Roy, Ochi and Dwivedi2022). In summary, RNA m6A methylation serves as a crucial regulation mechanism at the epitranscriptomic level to shape synaptic plasticity and behaviour (Malovic and Pandey, Reference Malovic and Pandey2023) and involvement of RNA m6A methylation appears in multiple aversive behavioural conditions, but so far evidence for a role of m6A methylation in anxiety has remained correlational.

In silico-based approach that departed from the crystal structure of the METTL3/METTL14 complex with S-adenosyl-L-homocysteine (Wang P et al., Reference Wang, Doxtader and Nam2016) enabled some of us to identify the first small molecule ligands that bind to METTL3/METTL14/WTAP and function as RNA methylation activators (Selberg et al., Reference Selberg, Blokhina, Aatonen, Koivisto, Siltanen, Mervaala, Kankuri and Karelson2019). In particular, as measured using the surface plasmon resonance, the compounds increased 200-fold the binding efficiency of the RNA methylating substrate S-adenosylmethionine to the METTL3/METTL14/WTAP protein complex. The treatment of cells with these compounds activating the METTL3/METTL14/WTAP complex led to the increase of m6A-modified nucleotides in the substrate RNA (Selberg et al., Reference Selberg, Žusinaite, Herodes, Seli, Kankuri, Merits and Karelson2021a). The cytotoxicity screening on the number of endpoints did not show any aberrations up to the 200 µM concentration of compounds (acute oral toxicity LD50 > 2 000 mg/kg). The compounds do not inhibit hERG channel and the monoamine oxidase. The compounds have good water solubility (>100 µM), and the half-lives of the best compound CHMA1004 are t1/2 = 354 min in human plasma and t1/2 = 60.5 min in rat plasma, respectively (data available on request). The computational modelling on other methyltransferases (DNA and histone methyltransferases) did not show any strong binding of the m6A methylation activators to these enzymes, thus they are predicted to be very specific to METTL3/METTL14/WTAP protein complex. Two of these compounds were recently found to support the survival of mouse dopamine neurons in growth factor deprivation assay and to protect them from 6-OHDA toxicity. Importantly, CHMA1004 (methyl-piperazine-2-carboxylate) protected stem cell-derived human dopamine neurons against 6-OHDA-induced cell death. Furthermore, intrastriatal administration of CHMA1004 attenuated amphetamine-induced ipsilateral rotations after 6-OHDA treatment and increased dopamine neuron fibre density in the rat model of PD (Yu et al., Reference Yu, Selberg, Teino, Nam, Ivanova, Renzi, Seli, Kankuri, Voutilainen, Karelson and Saarma2023). Given the massive arborisation of dopaminergic axons (Matsuda et al., Reference Matsuda, Furuta, Nakamura, Hioki, Fujiyama, Arai and Kaneko2009), dopaminergic neural projections may be particularly sensitive to the epitranscriptomic regulation of mRNA translation that occurs in the axonal and dendritic ribosomes (Biever et al., Reference Biever, Glock, Tushev, Ciirdaeva, Dalmay, Langer and Schuman2020).

Thus, taking into account the indirect evidence reviewed above on the role of RNA m6A methylation in anxiety and depression, and the implication of dopaminergic neurotransmission in these disorders that often are co-morbid with addictions (Rehm et al., Reference Rehm, Anderson, Manthey, Shield, Struzzo, Wojnar and Gual2016), we aimed to examine the effect of the novel METTL3/METTL14 activator compound CHMA1004 in a battery of anxiolytic and antidepressant screening tests in the rat, its effect on emission of ultrasonic vocalisations indicative of positive and negative affect in the rat, and effects of acute amphetamine treatment. The effect of subchronic administration of CHMA1004 on monoamine neurochemistry in several brain regions and on gene expression in striatum, and its acute effect on dopamine overflow, were also examined.

Materials and methods

Animals

In the behavioural experiments, twenty-seven Wistar rats (15 female and 12 male; 3 months of age at the beginning of the experiment; body weight 290 ± 4 g and 499 ± 13 g, respectively; bred locally, parent rats obtained from Charles River) were given CHMA1004 (1 mg/kg or 5 mg/kg) dissolved in saline intraperitoneally (IP; 1 ml/kg) once in a day. Control group received vehicle (saline). Each group consisted of 5 female and 4 male rats. In the microdialysis experiment, 19 male and 16 female rats with the same background were used; their allocation to treatment groups is given in the Microdialysis section. Animals were housed 4–5 per cage in standard transparent polypropylene cages under controlled light cycle (lights on 08:00-20:00) and temperature (22°–24°C), with free access to water and food pellets (v1534–000 universal maintenance diet, ssniff Spezialdiäten GmbH, Soest, Germany). The studies were approved by the Ethical Committee for Animal Experiments of the Estonian Ministry of Rural Affairs (permission no 149).

General procedure

The procedure is illustrated in Figure 1. Thirty minutes after the first IP administration of CHMA1004 or vehicle, an open field test was conducted with recording of ultrasonic vocalisations (USVs). A series of behavioural tests were performed in the course of daily drug treatment, while all behavioural tests were started 30 min after drug administration. After 17 or 18 days of daily drug treatment, rats were sacrificed 30 min after the last injection by decapitation, the brains dissected as previously described (Liiver et al., 2023) and samples immediately frozen in liquid nitrogen. All procedures, except the sucrose preference test, occurred during the lights-on phase. Animals were transported 30 min before the start of testing to the room of behavioural experiments to habituate, and tested (within the sex group) in randomised order. All animal experiments were in compliance with the ARRIVE guidelines and were carried out in accordance with EU Directive 2010/63/EU for animal experiments; the experimental protocol was approved by the Animal Experimentation Committee at the Estonian Ministry of Rural Affairs.

Fig. 1. Experimental design of behavioural and ex vivo biochemical studies of the effect of acute and subchronic treatment with CHMA1004. Note: HPLC, high-performance liquid chromatography; IP, intraperitoneal; RNA-seq, RNA sequencing; USVs, recording of ultrasonic vocalisation.

Open field test

Open field test was carried out in a 0.5 × 1.0 m metal box with side walls 40 cm high (Tõnissaar et al., Reference Tõnissaar, Mällo, Eller, Häidkind, Kõiv and Harro2008). The area was divided into eight squares of equal size. Time spent moving around was recorded, and square crossings and rearing in 4 min were counted, as was the number of faecal boli left on the open field. The apparatus was cleaned with laboratory tissue dampened with water after each animal.

USVs

The microphone of an ultrasound recording device (Avisoft Ultra Sound Gate 116–200, Avisoft Bioacoustics, Berlin, Germany) was located about 30 cm from the floor of the apparatus, recording USVs with a sampling rate of 300 kHz in 16-bit format on a computer hard drive (Kõiv et al., Reference Kõiv, Matrov, Uusen and Harro2021). The files were recorded with Avisoft-RECORDER software and later analysed with Avisoft SASLab Pro (Avisoft Bioacoustics, Berlin, Germany) software. Spectrograms were created using the Fast Fourier Transform algorithm (1024 FFT length, 75% frame, Hamming window, and 75% time window overlap). USVs were manually marked and saved on the spectrogram by an observer blind to experimental conditions. Fifty kHz USVs with frequencies over 35 kHz, including flat, frequency-modulated and trill calls) and 22-kHz (<35 kHz) USVs were counted.

Elevated zero-maze

The elevated zero-maze test (Shepherd et al., Reference Shepherd, Grewal, Fletcher, Bill and Dourish1994) was conducted as previously described (Matto et al., Reference Matto, Harro and Allikmets1997). An elevated annular platform (outer Ø 105 cm, width 10 cm, 72 cm above the floor) was equally divided into two opposing enclosed quadrants (height of the walls 28 cm) that were connected by open quadrants (height of the edge 1 cm, 5–11 Lx in the open part). The open quadrants were also divided into three equidistant parts to quantify the locomotor activity. The rat was placed at the outset of one of the closed quadrants and was observed at a distance for 5 min. The measures taken were stretched-attend postures toward the open quadrants, latency to enter the open quadrant with all four paws (s), entries into the open quadrants, and time spent (s) in the open quadrants. The apparatus was cleaned with moist tissue after every rat.

Social interaction test

The social interaction test introduced by Sandra File (reviewed by File and Hyde, Reference File and Hyde1978) was conducted as described previously (Tõnissaar et al., Reference Tõnissaar, Philips, Eller and Harro2004) with the variance that animals were not single-housed. Two unfamiliar, weight-matched rats were placed in the opposite corners of a brightly lit chamber (30 × 30 × 60 cm) with floor covered with wood shavings. The behaviour of the animals was concurrently observed for 10 min and the total time spent in active social behaviour (allogrooming, sniffing the partner, crawling under and over, following) was recorded.

Forced swimming test

The forced swimming test (Porsolt et al., Reference Porsolt, Anton, Blavet and Jalfre1978) was carried out as described previously (Mällo et al., Reference Mällo, Matrov, Kõiv and Harro2009). On the first day of the experiment, the procedure lasted 15 min and the re-exposure 24 h later lasted 5 min. The measurements were based on the behavioural categories described by Armario et al. (Reference Armario, Gavaldà and Martí1988); the first 5 min period of the first day was counted. A rat was judged to be immobile when it remained floating in the water with all limbs motionless or made minimal movements in order to maintain its head above the water and to struggle whenever it made intense movements of all the four limbs with the two front paws breaking the surface of the water or touching the walls of the tank. The time spent in swimming represents the rat making active swimming motions, more than necessary to merely maintain its head above water.

Amphetamine treatment effect

D-Amphetamine (1 mg/kg IP) was dissolved in saline and administered 15 min prior to the animal being placed in a standard group-housing cage (55 × 33 cm) for assessment of stimulant-induced locomotor response and ultrasonic vocalisation (Kõiv et al., Reference Kõiv, Metelitsa, Vares, Tiitsaar, Raudkivi, Jaako, Vulla, Shimmo and Harro2016). Horizontal locomotion was assessed from a 12 min video by counting line crossings with four paws of the boundaries of six equal squares on the cage floor, rearings were counted when the animal stood on hind legs with the body at least at a 45° angle. Data were analysed by 2 min bins.

Locomotion and ultrasound vocalisations in a novel home-cage-like environment

Measurement of USVs and spontaneous horizontal and vertical locomotion was conducted in a standard novel home-cage-like box (for 8 min) with a similar protocol as previously described for amphetamine-induced behavioural measurements.

Measurement of sucrose intake and preference

The sucrose preference test was carried out during the dark phase with no previous food or water deprivation (Kõiv et al., Reference Kõiv, Vares, Kroon, Metelitsa, Tiitsaar, Laugus, Jaako and Harro2019). The animals were placed into single cages 1 h before the sucrose preference testing and group-housed again after that. Food was available during sucrose preference test. Two bottles were used, one filled with 1% sucrose solution and the other with water. Liquid consumption was measured by weighing previously weighed sucrose solution and water bottles 1h and 12 h after the beginning of the sucrose preference test.

High-performance liquid chromatography (HPLC) in ex vivo studies

Monoamines and their metabolites were assayed by HPLC with electrochemical (amperometric) detection as described previously (Kõiv et al., Reference Kõiv, Vares, Kroon, Metelitsa, Tiitsaar, Laugus, Jaako and Harro2019). Rat brain tissues were homogenised with an ultrasonic homogeniser (Bandelin Sonopuls, Bandelin Electronic, Berlin, Germany) in ice-cold solution of 0.1 M perchloric acid containing 5 mM of sodium bisulphite and 0.4 mM EDTA to avoid oxidation. The homogenate was then centrifuged at 14 000 r.p.m. for 10 min at 4°C. Aliquots (10 μl) of the supernatant obtained were chromatographed on a Luna C18(2) column (150 mm × 2 mm, 5 μm). The separation was done in isocratic elution mode at column temperature of 30°C, using the mobile phase containing 0.05 M sodium citrate buffer at pH 3.7, 0.02 mM EDTA, 1 mM KCl, 1 mM sodium octanesulphonate and 7.5% acetonitrile. The two used chromatography systems consisted of an isocratic pump (HP1100, Agilent, Waldbronn, Germany or Shimadzu LC-20AD, Japan), a temperature-regulated autosampler (HP1100 or Shimadzu SIL-20AC), a temperature-regulated column compartment and an HP 1049 electrochemical detector (Agilent, Waldbronn, Germany) with a glassy carbon electrode (HP flow cell or Antec VT-03 flow cell). The measurements were done at an electrode potential of + 0.7 V versus the Ag/AgCl reference electrode. The limits of detection at signal-to-noise ratio = 3 were as follows (expressed as pmol/mg tissue for each): 0.08 for dopamine (DA), 0.10 for homovanillic acid (HVA), 0.05 for 3,4-dihydroxyphenylacetic acid (DOPAC), 0.08 for serotonin (5-HT), 0.04 for 5-hydroxyindoleacetic acid (5-HIAA), 0.07 for noradrenaline (NA), 0.03 for normetanephrine (NMN), and 0.01 for 3-methoxytyramine (3-MT).

RNA extraction and mRNA sequencing

RNA was extracted from female rat striatal samples using FavorPrepTM Tissue Total RNA Mini Kit (FATRK 001-2, Favorgen, Taiwan) according to the manufacturer’s protocol. Briefly, the FARB lysis buffer, containing ß-mercaptoethanol, was added to the tissue sample, grinding was performed with a micropestle followed by homogenisation through a 20-G needle syringe. The lysate was cleared by using the filter column followed by addition of one volume of 70% ethanol. The sample was then transferred to the FARB Mini Column, washed and dried by centrifugation at 18,000 g for three minutes. Total RNA was eluted in 30 µl of nuclease-free water, concentration was measured by NanoDrop 1000 spectrophotometer (Thermo Fisher Scientific Inc., Waltham, MA, USA) and stored at −80°C.

Sequencing was performed by Novogene Co., Ltd (Beijing, China). Quality control and preprocessing of the fastq files, containing raw reads from Illumina 150bp paired end sequencing, was conducted with fastp v0.23.2 using default parameters (Chen S et al., Reference Chen, Hartman, Haavik, Harro, Klungsøyr, Hegvik, Wanders, Ottosen, Dalsgaard, Faraone and Larsson2018). Strand-specific alignment to the rat rn7 genome was carried out using HISAT2 v2.2.1. (paired end, Strand-FR). FeatureCounts v1.6.4 (Liao et al., Reference Liao, Smyth and Shi2014) was used for transcript counting using Rattus_norvegicus.mRatBN7.2.108.gtf annotation from Ensembl database as a reference (Martin et al., Reference Martin, Amode, Aneja, Austine-Orimoloye, Azov and Barnes2023). Differential gene expression between pairwise groups was performed with DeSeq2 v22.11.40.6 (Love et al., Reference Love, Huber and Anders2014) using poscounts method for size factor estimation. Genes with p-value smaller than 0.05 were considered differentially expressed. SynGO portal was used for ID conversion (Koopmans et al., Reference Koopmans, van Nierop, Andres-Alonso, Byrnes, Cijsouw and Coba2019). GO enrichment analysis was performed with DAVID 2021 (Sherman et al., Reference Sherman, Hao, Qiu, Jiao, Baseler, Lane, Imamichi and Chang2022). Two samples were excluded from the subsequent analysis because of low alignment rate (<5%) to the reference rn7 genome, and one sample was considered as an outlier by visual inspection of PCA plot and omitted from subsequent analysis; consequently, further analysis included 4 rats per group. Volcano plot was created using R ggplot2-based tool in Galaxy (Version 0.0.3) using output from DESeq2.

Microdialysis experiments

In vivo microdialysis was carried out to assess the effect of acute treatment with CHMA1004 on extracellular levels of dopamine in striatum, essentially as previously described (O’Leary et al., Reference O’Leary, Kõiv, Raudkivi and Harro2016). The animals were anaesthetised with ketamine/medetomide anaesthesia (60 mg/kg and 0.5 mg/kg IP, respectively) and mounted in a Stoelting stereotactic frame. A self-made concentric Y-shaped microdialysis probe with 7 mm shaft length and 3 mm active tip was implanted into the left dorsal striatum according to the following coordinates: AP + 0.7; ML + 3.0; DV –7.0, according to Paxinos and Watson (Reference Paxinos and Watson1986). The dialysis membrane used was polyacrylonitrile/sodium methyl sulphonate copolymer (Filtral 12; i.d.: 0.22 mm; o.d.: 0.31 mm; AN 69, Hospal, Bologna, Italy). Two stainless steel screws and dental cement were used to fix the probe to the scull. After the surgery, anaesthesia was reversed by administration of atipamezole (1 mg/kg, SC) and meloxicam (1–2 mg/kg, SC) was administered for post-operative analgesia. The rats were placed in 21 cm × 36 cm × 18 cm individual cages in which they remained throughout the microdialysis experiment. Rats were given about 24 h for recovery and microdialysis procedure was conducted in awake freely moving animals. The microdialysis probe was connected via FEP tubing (i.d. 0.12 mm, AgnTho’s AB, Lidingö, Sweden) to a syringe pump (SP101, World Precision Instruments, Inc., Sarasota, FL, USA) and to a refrigerated microsampler (Univentor 820, Univentor Limited, BLB029 Bulebel Industrial Estate, Zejtun ZTN 3000, MALTA) and perfused with Ringer solution (147 mM NaCl, 4 mM KCl, 1.2 mM CaCl2, 1.0 mM MgCl2, 1.0 mM Na2HPO4; pH 7.20–7.22) at a constant rate of 1.5 μl/min. After connecting the animal to the microdialysis system, the perfusate was discarded during the first 60 min to allow stabilisation. Six baseline samples were collected, followed by an i.p. injection of either saline (male rats; n = 2), CHMA1004 1 mg/kg (male, n = 6; female, n = 6), CHMA1004 5 mg/kg (male, n = 6; female, n = 6) or amphetamine 0.5 mg/kg (male, n = 4; female, n = 4), after which another 12 samples were collected. All samples were collected at 15 min intervals. Sample vials were prefilled with 7.5 μl of 0.02 M acetic acid to prevent oxidation of dopamine. The quantity of dopamine and serotonin in the microdialysis samples was determined by high-performance liquid chromatography with electrochemical detection. The chromatography system consisted of a Shimadzu LC-10AD pump and CBM-20A controller (Shimadzu Corporation, Kyoto, Japan), a Luna C18(2) 5 μm column (150 mm × 2 mm) kept at 30°C and Decade II digital electrochemical amperometric detector (Antec Leyden BV, the Netherlands) with electrochemical flow cell VT-03 (2 mm GC WE, ISAAC reference electrode, Antec Leyden BV, the Netherlands). The mobile phase consisted of 0.05 M sodium citrate buffered to pH 5.3, 2 mM KCl, 0.02 mM EDTA, 3.5 mM sodium octyl sulphonate and 14% acetonitrile. The mobile phase was pumped through the column at a rate of 0.2 ml/min. 5-HT and DA eluted from the column were measured with a glassy carbon working electrode maintained at a potential of + 0.4 V versus Ag/AgCl reference electrode. Data were acquired using a Shimadzu LC Solution system.

Data analysis

Behavioural and neurochemical data were first analysed by two-factor ANOVA (Sex and Treatment) with repeated measures (Test) if appropriate. Because the apparent dose dependency of the effect of CHMA1004 was in some tests different between sexes, this was followed by one-way ANOVA (CHMA1004) with male and female dosage group pooled or separate if the interaction term was closely missing p < 0.05. Owing to small group sizes and apparent difference in dose-response between male and female animals, subsequent pairwise comparisons were made by Fisher’s least significant difference test. Significance threshold was set at 0.05. Statistical analysis was made using the IBM SPSS Statistics 26 software.

Results

Body weight and weight of adrenals

Body weight gain was similar in all treatment groups. After the 18-day treatment with CHMA1004, no statistically significant effect was found on adrenal weight (data not shown).

Effect of acute and repeated administration of CHMA1004 on open field behaviour

Acute treatment with CHMA1004 increased time spent in active movement on the open field (Figure 2); in female rats, this effect was statistically significant at 1 mg/kg and still observable after repeated treatment (Supplementary Fig. 1) despite being smaller as the activity became lower (Treatment × Sex F2,21 = 3.49; p < 0.05; Test F2,42 = 33.5; p < 0.001). This increase of activity was statistically significant only in female rats while at the dose of 5 mg/kg a similar tendency appeared in males. Similar group differences were reflected in rearings (Treatment F2,21 = 4.27; p < 0.05; Sex F1,21 = 7.96; p < 0.05; Treatment × Sex F1,21 = 3,42; p = 0.052; Test F2,42 = 10.3; p < 0.001), and in the number of line crossings in which case the tendency of increase in males at the dose of 5 mg/kg was the closest to statistical significance (Treatment F2,21 = 3.37; p = 0.054; Sex F1,21 = 6.88; p < 0.05; Treatment × Sex F2,21 = 4.59; p < 0.05; Test F2,42 = 11.3; p < 0.001). The number of faecal boli left on the open field at the first test, a rather sensitive anxiety measure, was larger in male rats (Figure 2), and this was reduced by CHMA1004 treatment (Sex F1,21 = 21.4; p < 0.001; Treatment F1,21 = 3.46; p = 0.05; Treatment × Sex F2,42 = 4.78; p < 0.05). Thus, CHMA1004 treatment elicited anxiolytic- or psychostimulant-like locomotor activation, especially in female rats, and reduced emotionality in male rats as expressed in mitigation of defecation. Acute CHMA1004 treatment had no statistically significant effect on 50-kHz USVs (Figure 2). Only one rat produced 22-kHz USVs, suggestive that neither testing nor the drug was perceived as anxiogenic. Emission of 50-kHz USVs was also low, with a tendency of more vocalisation after CHMA1004 treatment. In conclusion, acute administration of CHMA1004 was without strong negative or positive emotional valence.

Fig. 2. Effect of acute administration of CHMA1004 on locomotion, defecation and emission of ultrasonic vocalisations in an open field test. Note: *p < 0.05; **p < 0.01 versus vehicle. Post hoc test on rearings in male rats CHMA1004 (5 mg/kg) p = 0.069.

Effect of repeated administration of CHMA1004 on behaviour in elevated zero-maze

Time spent in the open sections of the elevated zero-maze was longer on the second test (Test F1,20 = 16.7; p < 0.001), especially in females. CHMA1004 treatment had no statistically significant effect while the measure tended to be increased (Figure 3). Repeated CHMA1004 treatment increased the number of head-dips at the first zero-maze test (CHMA1004 F2,20 = 5.35; p < 0.05); this effect was more prominent in female rats. The number of entries into the open section was significantly higher at the second testing (Test F1,20 = 18.6; p < 0.001), and was significantly increased in male rats of the CHMA1004 5 mg/kg group as compared to the control at the second elevated zero-maze test (Sex × Test F1,20 = 5.77; p < 0.05). Stretched-attend postures were reduced at the second test (Test F1,20 = 19.3; p < 0.001), however the male CHMA1004 5 mg/kg treatment group had this measure reduced already at the first elevated zero-maze test (Treatment × Sex × Test F2,20 = 3.09; p = 0.067; CHMA1004 F2,20 = 8.03; p < 0.01). In sum, CHMA1004 had an anxiolytic-like effect after repeated administration.

Fig. 3. Effect of repeated administration of CHMA1004 in elevated zero-maze test. Note: SAPs, number of stretched-attend postures; *p < 0.05 versus vehicle.

Effect of repeated administration of CHMA1004 in social interaction test

In rats housed in groups, social interaction with an animal from another cage was very low (40–50s on average) in female rats, and was reduced in the male CHMA1004 5 mg/kg group to the level comparable to females (Sex F1,20 = 27.4; p < 0.001; Treatment × Sex F2,20 = 3.33; p = 0.057; data not shown).

Effect of repeated administration of CHMA1004 on novelty-induced USVs

Not a single rat produced any 22-kHz UVSs during this test, but 50-kHz USVs were emitted by all animals and frequently by most of them. Two-way ANOVA revealed no statistically significant Treatment effect but a significant Sex effect (F1,21 = 8.08; p < 0.01). Male control rats had higher 50-kHz emission than female controls, but CHMA1004 treatment tended to reduce this to the level of females (Supplementary Fig. 2).

Effect of repeated administration of CHMA1004 in the forced swimming test

In this screening test for antidepressant potential, CHMA1004 treatment had no significant effect on immobility (Figure 4), but increased swimming (Treatment F 2,21 = 3.75; p < 0.05; Test F1,21 = 40.6; p < 0.001; Sex × Test F1,21 = 16.5; p < 0.001) and reduced struggling (Test F1,21 = 13.8; p < 0.01; Sex × Test F1,21 = 11.9; p < 0.01) in female rats on both days of forced swimming, at both dose levels.

Effect of repeated administration of CHMA1004 on effects of amphetamine

This was studied in order to assess the potential of RNA m6A methylation to cause sensitisation to psychostimulant action. After acute administration of amphetamine (1 mg/kg) female rats made more rearings than males (Supplementary Fig. 3A), especially during the latter part of the test (Sex F1,21 = 27.8; p < 0.001; Test F5,105 = 26.4; p < 0.001; Sex × Test F1,21 = 2.88; p < 0.05), for horizontal activity the differences were not statistically significant. In male rats, locomotion tended to be higher in the CHMA1004 (5 mg/kg) group, but this effect was not statistically significant (Supplementary Fig. 3AB). After treatment with amphetamine, only three male rats and no female rat emitted 22-kHz USVs, but as expected the emission of 50-kHz USVs was robust (Supplementary Fig. 3C). Female rats produced more 50-kHz USVs than males, but this was dependent on CHMA1004 treatment (Sex F1,21 = 9.17; p < 0.01; Treatment × Sex F2,21 = 3.52; p < 0.05). Thus, in male rats subchronic CHMA1004 treatment tended to increase the emission of 50-kHz USVs, but in females this treatment reduced vocalisation, so that the large male-female difference in control groups was much smaller in CHMA1004 treatment groups.

Effect of repeated administration of CHMA1004 on sucrose intake and preference

CHMA1004 treatment did not statistically significantly affect sucrose intake or preference (data not shown), but if intake of 11 h was adjusted to body weight, an increase was found in female rats of the CHMA1004, significantly in the 1 mg/kg group, as compared to controls (Supplementary Fig. 4; Sex F1,21 = 31.2; p < 0.001).

Effect of repeated administration of CHMA1004 on monoamine neurochemistry

Subchronic treatment with most psychotropic medications causes adaptive changes in one or several monoamine systems, and indeed CHMA1004 affected the levels of monoamine neurotransmitters and their metabolites in several brain regions (Table 1). In striatum, CHMA1004 treatment increased the levels of NMN (Treatment F2,21 = 4.30; p < 0.05; Sex F1,21 = 3.41; p < 0.05) in females and dopamine turnover in males (Treatment × Sex F2,21 = 4.44; p < 0.05). Frontal cortical NMN levels were increased by subchronic treatment with CHMA1004 in males (Treatment F2,21 = 4.94; p < 0.05; Treatment × Sex F2,21 = 4.73; p < 0.05). Similar effect on noradrenaline levels closely missed the conventional level of significance (Treatment F2,21 = 3.40; p = 0.0525), and both NA and NMN tended to be higher also in females at the lower dose of CHMA1004. In locus coeruleus, CHMA1004 treatment reduced HVA levels in female rats, while the effect was opposite in males (Treatment × Sex F (2,21) = 4.66; p < 0.05). Ventral tegmental area was also dose-dependently affected by CHMA1004 treatment: Levels of noradrenaline and HVA were higher after 1 mg/kg (Treatment F2,21 = 5.22 and 4.21, respectively; p < 0.05; post hoc tests not significant owing to small groups by sex and large male-female differences), and 5-HIAA levels would be significantly affected in the same manner if one outlier in the vehicle group were excluded (Treatment F2,20 = 6.42; p < 0.01). In hippocampus and in the raphe area, two-way ANOVA did not reveal any statistically significant treatment effect.

Table 1. Monoamine levels after two and half weeks of administration of the METTL3/METTL14 activating compound CHMA1004 in male and female rats

Monoamine levels are expressed as pmol/mg tissue. *p < 0.05; **p < 0.01 versus vehicle; # p < 0.05 versus CHMA1004 (1 mg/kg); ¤ p = 0.055 versus vehicle.

Extracellular levels of dopamine in striatum

To assess whether METTL3/14 activation could be triggering dopamine release, we conducted a separate microdialysis experiment. As expected, acute administration of saline had no effect of striatal dopamine levels, while amphetamine (0.5 mg/kg) increased dopamine release in male (F51,187 = 6.88; p < 0.001) and female (F34,170 = 5.58; p < 0.001) rats. CHMA1004 at either dose did not elicit any significant effect on extracellular dopamine levels in striatum (Figure 5).

Fig. 4. Effect of repeated administration of CHMA1004 in the forced swimming test. Note: *p < 0.05 versus vehicle.

Fig. 5. Extracellular levels of dopamine after acute administration of amphetamine or CHMA1004. Note: A: female rats; B: male rats.

Effect of subchronic administration of CHMA1004 on striatal gene expression

Because it is not known which epitranscriptomic changes should be expected to be most responsive to METTL3/METTL14 activation, a hypothesis-free approach, RNA sequencing was performed. Differential expression analysis resulted in statistically significant change (Wald test p < 0.05) of 172 and 120 genes following daily 1 mg/kg and 5 mg/kg of CHMA1004, respectively (Supplementary Table 1). Of these genes, seven were overlapping (Galns, Hps3, Misp3, Nbl1, Samd4a, Slc28a2, Spon2), with similar change in expression. Among the genes that were affected with the highest statistical significance (Figure 6), most were of specific interest, for example Acsm5, Peli1, Pcdhb5, Psmb9, Fgr, Lrp10, Grb83, Hes1, Vwa5b2, Mx2, Serpinf2, Nbl1 and Mdga1; of other differentially expressed genes, Nxph3 would deserve particular attention (see Discussion). Pathway analyses (Supplementary Table 2) suggested the involvement of, for example, phospholipase C-activating G-protein coupled receptor signalling, positive regulation of phosphatidylinositol 3-kinase signalling metabolic pathways, cell adhesion, and mitochondrial matrix in the action of CHMA1004.

Fig. 6. Effect of repeated administration of CHMA1004 on striatal gene expression in female rats. Note: Volcano plots of the differentially expressed gene analysis with DESeq2 between control and CHMA1004 1 mg/kg (left) or 5 mg/kg (right). Significantly differentially expressed genes (p < 0.05) are labelled blue (downregulated genes) or red (upregulated genes). The names of the top 10 genes with the lowest p-values are shown.

Discussion

This study revealed that acute and subchronic administration of a novel METTL3/METTL14 activator compound CHMA1004 can elicit activating and anxiolytic-like effects in the rat. Subchronic treatment with CHMA1004 affected monoaminergic neurotransmission in the brainstem and forebrain areas. These effects were region- and dose-dependent and not all identical in male and female animals, but in general were suggestive of an increase in catecholamine function. Subchronic CHMA1004 treatment also elicited many alterations in striatal gene expression that are relevant to stress response.

Acute treatment effect of CHMA1004 was assessed in an open-field test. Increased locomotor activity that was observed can be indicative of anxiolytic action, but also of a psychostimulant effect. Similarly, lower levels of locomotion do not necessarily reflect anxiety, and in the rat, ultrasonic vocalisation provides the most reliable index of strong activation of emotion circuits (Panksepp, Reference Panksepp1998). Therefore, we simultaneously recorded USVs, as classic psychostimulants typically elicit 50-kHz USVs (Simola et al., Reference Simola, Ma and Schallert2010; Kõiv et al., Reference Kõiv, Vares, Kroon, Metelitsa, Tiitsaar, Laugus, Jaako and Harro2019; Wöhr, Reference Wöhr2022). No strong psychostimulant effect was observed, suggesting that the increased locomotor activity was reflective of either reduced anxiety or a non-dopaminergic psychostimulant effect. On the other hand, rats did not produce 22-kHz USVs. Previously, administration of neuropeptide S that also displays mixed anxiolytic/psychostimulant profile in exploration-based tests (Xu et al., Reference Xu, Reinscheid, Huitron-Resendiz, Clark, Wang, Lin, Brucher, Zeng, Ly, de Lecea and Civelli2004; Rizzi et al., Reference Rizzi, Vergura, Marzola, Ruzza, Guerrini, Salvadori, Regoli and Calo2008) was found to elicit substantial 22-kHz vocalisation revealing predominantly negative emotional response behind behavioural activation (Kõiv et al., Reference Kõiv, Matrov, Uusen and Harro2021). CHMA1004 did, in contrast, display an anxiolytic-like and stimulatory effect without any significant affective response in terms of USV production. That the effect of acute administration of CHMA1004 was anxiolytic, was further suggested by the reduction of defecations on open field in male rats. Previously, acute restraint stress was shown to rapidly downregulate Mettl3 expression in the prefrontal cortex and amygdala of mice (Engel et al., Reference Engel, Eggert, Kaplick, Eder, Röh, Tietze, Namendorf, Arloth, Weber, Rex-Haffner, Geula, Jakovcevski, Hanna, Leshkowitz, Uhr, Wotjak, Schmidt, Deussing, Binder and Chen2018). Altogether the facts that METTL3/METTL14 activating compound CHMA1004 can elicit an anxiolytic-like activating effect in conditions of novelty stress suggest that RNA m6A methylation is involved in positive stress response.

Repeated treatment with CHMA1004 did not appear to produce significant tolerance to its anxiolytic-like/activating effect. The anxiolytic effect of CHMA1004 was further substantiated by results of the elevated zero-maze test where stretched-attend postures were reduced, open quadrant entries and head-dips increased, which is not a profile of a psychostimulant effect. At the first glance, the tendency of decrease in social interaction time in male rats may appear as contradictory to an anxiolytic effect, but the test was not conducted as validated for this purpose: We could not use single housing in this experiment because this would sensitise animals to stressful stimulation. It is likely that at this locomotor activating dose CHMA1004 increased exploration of the novel physical environment on account of social behaviour. Similar is the effect of amphetamine in the social interaction test (File and Hyde, Reference File and Hyde1979; unpublished results from this laboratory).

The forced swimming test has been the preferred screening tool for novel drugs with antidepressant potential, while the interpretation of results with this test has been subject to multiple interpretations (e.g. Harro, Reference Harro2019). Subchronic CHMA1004 treatment did not affect immobility but increased swimming and reduced struggling in female rats. The characteristic response to the first generation antidepressant drugs such as tricyclics is reduction of immobility (Porsolt et al., Reference Porsolt, Anton, Blavet and Jalfre1978) while the selective serotonin reuptake inhibitors increase swimming (Armario et al., Reference Armario, Gavaldà and Martí1988). It should be noted that in the forced swimming test used as a screening tool mostly acute drug treatment after the first swimming session is applied. In the present study, the drug had been repeatedly administered before the test and immobility was relatively low in these repeatedly handled rats. As expected, swimming increased and struggling decreased when CHMA1004 treated animals were submitted to forced swimming the second time. Previous studies applying chronic stress or partial serotonergic denervation have suggested that a behavioural activation that is present already at the first swimming session is indicative of higher expression of impulsivity (Harro et al., Reference Harro, Tõnissaar, Eller, Kask and Oreland2001). Thus, the present results could best be interpreted as subchronic CHMA1004 treatment presenting an impulsivity-reducing effect. This possibility should be examined further in behavioural tests specifically targeting different aspects of impulsivity (Dalley and Robbins, Reference Dalley and Robbins2017). In another test of antidepressant prediction, the sucrose preference test, sucrose preference was not affected but consumption as adjusted to body weight was increased by CHMA1004 treatment in females. Again, the test is validated in conditions of chronic stress-induced reduction of sucrose intake, and as part of the present test, battery informs little about potential antidepressant effect but suggests that subchronic CHMA1004 treatment is not stressful/anhedonic and rather may increase the appetitive value of carbohydrates. In this context it is also relevant that long-term treatment with CHMA1004 did not increase adrenal weight, again suggesting that this treatment was not strongly stressful; indeed, it may be important to notice that at the doses that were behaviourally most effective (1 mg/kg in females and 5 mg/kg in males, respectively), a tendency of reduction of adrenal weight appeared.

Fifty kHz USVs are indicative of positive affect in the rat (Burgdorf et al., Reference Burgdorf, Knutson, Panksepp and Ikemoto2001; Brudzynski and Granon, Reference Brudzynski, Granon and Harro2023), and in adult animals, the most powerful stimulus to elicit this response is psychostimulant drugs such as amphetamine. Novel home-cage-like environment can also elicit the 50-kHz response. Subchronic CHMA1004 treatment did not modify the effect of either novel cage or amphetamine treatment significantly, but statistical power to detect its effect may have been too low, because substantial differences between male and female control groups appeared. A general limitation in this study is the number of animals of each sex was small. This was to reduce the overall number of rats used and simultaneously provide comparison of male and female animals. It was not foreseen that dose dependency varies between males and females, and the significantly larger body weight of males may have been a contributory factor to this. Nevertheless, it can be concluded that CHMA1004 did not elicit substantial psychostimulant-like effects or produced sensitisation on its own. CHMA1004 treatment enhanced locomotor activity in response to amphetamine challenge in male rats, but their activity remained below of females. Altogether the potential of CHMA1004 to increase amphetamine-like stimulant responses appeared low.

Recent studies have highlighted the fact that dysregulation of m6A mRNA methylation occurs in neurodegeneration, and dopamine neurons are dependent on proper regulation of m6A mRNA methylation (Chen et al., Reference Chen, Yu, Guo, Zheng, Ali, Huang, Zhang, Wang, Huang, Qie and Wang2019). Inhibition of m6A RNA demethylases Fto and Alkbh5, which results in increased m6A levels, can support the survival of dopamine neurons (Selberg et al., Reference Selberg, Žusinaite, Herodes, Seli, Kankuri, Merits and Karelson2021a). Furthermore, the METTL3/METTL14 activating compound CHMA1004 was also found to support survival of cultured dopamine neurons that were deprived of growth factors and to protect neurons from 6-OHDA toxicity (Yu et al., Reference Yu, Selberg, Teino, Nam, Ivanova, Renzi, Seli, Kankuri, Voutilainen, Karelson and Saarma2023). Importantly, intrastriatal administration of CHMA1004 improved motor behaviour and increased tyrosine hydroxylase positive fibre density in mice after 6-OHDA lesion; it also protected human cultured dopamine neurons in that study. Hence, we measured monoamine levels ex vivo in striatum after long-term administration of CHMA1004, and having observed alterations in catecholamines, performed further measurements in a few other brain regions. The most prominent effect of subchronic treatment with CHMA1004 was observed in VTA at the lower dose level. Dopaminergic measures tended to be altered also in the other brainstem regions and several alterations were observed in measures of catecholaminergic neurotransmission, mostly increasing with dose. These results suggest a complex reorganisation of monoamine systems by a METTL3/METTL14 activating drug, possibly primarily in the dopamine neurons. Acute CHMA1004 treatment did however not produce any increase in striatal extracellular dopamine levels, as might have been suggested by its behaviourally stimulatory effects, which is consistent with no increase in 50-kHz USVs after either acute or subchronic CHMA1004 treatment.

Most of proposed pathogenetic mechanisms of depression, relating to neurogenesis, activity of the HPA axis, glial function, neuroinflammation and synaptic transmission have all been associated with mRNA m6A methylation (Li C et al., Reference Li, Xu, Huang, Wang, Wu, Li, Peng, Wang, Li, He and Miao2023). One genetic association study could link a variance in the RNA m6A demethylase encoding gene ALKBH5 with major depression in a large sample (Du et al., Reference Du, Rao, Wu, Ye, Liu, Hu, Xiu, Shen and Xu2015) and another associated with depression in an FTO variant (Samaan et al., Reference Samaan, Anand, Zhang, Desai, Rivera, Pare, Thabane, Xie, Gerstein, Engert, Craig, Cohen-Woods, Mohan, Diaz, Wang, Liu, Corre, Preisig, Kutalik, Bergmann, Vollenweider, Waeber, Yusuf and Meyre2013). The biochemical pathways involved in the association of mRNA m6A methylation with anxiety or depression-related behavioural changes have however remained completely unknown. To extend knowledge beyond the monoamine systems and to analyse what mRNAs are regulated by METTL3/METTL14 activator CHMA1004 we also performed RNA sequencing in striatum, where CHMA1004 had been found acting on dopamine neurons in vivo (Yu et al., Reference Yu, Selberg, Teino, Nam, Ivanova, Renzi, Seli, Kankuri, Voutilainen, Karelson and Saarma2023) and where the authors had observed an increase in Mettl3 expression after an anxiety-reducing manipulation (unpublished). This was done in female rats. The low and high dose of CHMA1004 resulted in a different landscape of differential gene expression, suggesting that response to subchronic drug treatment involved dose-dependent adaptive changes that had also been apparent in behavioural and biochemical effects; however, expression of several genes was altered in a similar manner. Many effects of CHMA1004 on gene expression agree with what is known of m6A in neurodegeneration and affective disorders (see Introduction), and the capacity of the drug to protect dopamine neurons (Yu et al., Reference Yu, Selberg, Teino, Nam, Ivanova, Renzi, Seli, Kankuri, Voutilainen, Karelson and Saarma2023) as well as its behavioural effects in the present study. For example, Galns and Samd4a have been implicated in neuronal viability and plasticity; Galns that encodes a lysosomal exohydrolase has been considered a gene related to Parkinson’s disease (Wettergren et al., Reference Wettergren, Gussing, Quintino and Lundberg2012) and Samd4a is a RNA binding protein that can bind to a variety of target mRNAs through stem-loop structures, also known as Smaug recognition elements (SREs). Notably, Samd4a can regulate the mRNA stability, degradation and translation (Wang and Zhang, Reference Wang and Zhang2023). Fgr (Gardner-Rasheed feline sarcoma viral (v-fgr) oncogene homolog) and Psmb9 (proteasome subunit beta type-9), two innate immunity genes, both of which were affected by CHMA1004 treatment, each were recently associated with specific genetic forms of Parkinson’s disease (Shani et al., Reference Shani, Goldstein, Gana-Weisz, Bar-Shira, Thaler, Gurevich, Mirelman, Giladi, Alcalay and Orr-Urtreger2023). Another differentially expressed gene, Lrp10 (low-density lipoprotein receptor-related protein 10), has similarly been linked to Parkinson’s disease but also other neurodegenerative diseases, whereas its role has been attributed to astrocytic vesicle trafficking (Grochowska et al., Reference Grochowska, Carreras Mascaro, Boumeester, Natale, Breedveld, Geut, van Cappallen, Boon, Kievit, Sammler, Bank, Parchi, Cortelli, Alessi, van de Berg, Bonifati and Mandemakers2021), and its contribution to neurodegeneration suggested to be significant specifically in females (Guo et al., Reference Guo, Cao, Hou, Li, Huang, Zhu, Zhang, Lee, Lemos Duarte, Zhou, Wang, Liu, Martens, Chao, Goate, Bu, Haroutunian, Cai and Zhang2023). CHMA1004 treatment elicited upregulation of several cognition-related genes such as Acta2, Mx2, Spon2 and Peli1. Peli1, encoding Pellino E3 ubiquitin protein ligase 1, was recently shown to mediate methamphetamine-induced neuroinflammation whereas its downregulation was protective (Xu et al., Reference Xu, Yang, Lou, Chen, Wang, Hu, An, Gao, Wang and Chen2023); similarly, the reduction of Peli1 function reduced the lipopolysaccharide-induced inflammatory response in microglia (Dai et al., Reference Dai, Yuan, Wang, Xu, Mao and Xiao2019). Given that neuroinflammatory mechanisms are set in the brain by stress (Black, Reference Black2002), the reduction of Peli1 expression by CHMA1004 treatment may reflect a drug effect facilitating adaptive stress response to the experimental procedures. Mx2, an interferon-inducible GTP-binding protein, has been associated with the beneficial effect of physical exercise in a mouse model of Alzheimer’s disease (Widjaya et al., Reference Widjaya, Cheng, Kuo, Liu, Cheng and Lee2023), and Spon2 that can induce neurite outgrowth with learning effects in monkeys (Matsunaga et al., Reference Matsunaga, Nambu, Oka, Tanaka, Taoka and Iriki2015). Serpinf2 that encodes a serine protease inhibitor has a link to Alzheimer’s disease (De Jager et al., Reference De Jager, Srivastava, Lunnon, Burgess, Schalkwyk, Yu, Eaton, Keenan, Ernst, McCabe, Tang, Rai, Replogle, Brodeur, Gabriel, Chai, Younkin, Younkin, Zou, Szyf, Epstein, Schneider, Bernstein, Meissner, Ertekin-Taner, Chibnik, Kellis, Mill and Bennett2014), but was recently found as differentially expressed in a rat model of depression (Liao et al., Reference Liao, Mi, Zeng, Wei, Dai, Ye, Chen and Zhang2023). Depression-related is also Mdga1, encoding a cell adhesion molecule (Li et al., Reference Li, Kresock, Kuplicki, Savitz and McKinney2022), as well as a few other differentially expressed genes. Pcdhb5 or protocadherin beta-5 has been reported as a gene with altered expression in suicide victims (Thalmeier et al., Reference Thalmeier, Dickmann, Giegling, Schneider, Hartmann, Maurer, Schnabel, Kauert, Möller and Rujescu2008). It belongs to the group of cadherin-like cell adhesion molecules and knowledge on its function is limited but likely to include the establishment of specific neural connections. Nbl1, a candidate neuroblastoma suppressor gene has also been associated with suicide (Han et al., Reference Han, DiBlasi, Monson, Shabalin, Ferris, Chen, Fraser, Yu, Staley, Callor, Christensen, Crockett, Li, Willour, Bakian, Keeshin, Docherty, Eilbeck and Coon2023).

In turn, regulation of Vwa5b2 at the DNA methylation level has been related to lithium response in bipolar disorder (Marie-Claire et al., Reference Marie-Claire, Lejeune, Mundwiller, Ulveling, Moszer, Bellivier and Etain2020). Yet another prominent effect of CHMA1004 was downregulation of Acsm5 that encodes a mitochondrial enzyme that catalyses the activation of fatty acids by coenzyme A. While very little is known of the brain function of Acsm5, it has been identified as a top attention-deficit/hyperactivity disorder (ADHD)-related gene in a bioinformation transcriptome study (Qi et al., Reference Qi, Wang, Zhang, Liu, Wen, Ma, Cheng, Li, Cheng, Du, Liang, Zhao, Ding and Zhang2019), and ADHD is a vulnerability factor for several other psychiatric disorders including anxiety disorders, depression and bipolar disorder (Chen Q et al., Reference Chen, Zhou, Chen and Gu2018). A few other genes with substantially altered expression by CHMA1004 treatment have been linked to anxiety and stress response: expression of Hes1, an essential Notch signalling effector was associated with neuronal excitability in amygdala and with anxiety in the mouse (Matsuzaki et al., Reference Matsuzaki, Yoshihara, Ohtsuka and Kageyama2019), and knockdown of Gpr83, a glucocorticoid-regulated receptor has been associated with anxiety in female mice and a decrease in receptor expression to dexamethasone treatment, whereas anxiety was reduced in males (Fakira et al., Reference Fakira, Lueptow, Trimbake and Devi2021). Indeed, a thorough study on the role of m6A methylation in stress response regulation reported that the response of m6A levels upon glucocorticoid receptor stimulation was altered in patients with major depression (Engel et al., Reference Engel, Eggert, Kaplick, Eder, Röh, Tietze, Namendorf, Arloth, Weber, Rex-Haffner, Geula, Jakovcevski, Hanna, Leshkowitz, Uhr, Wotjak, Schmidt, Deussing, Binder and Chen2018). In that study, one of the most prominent alterations by knocking out the gene responsible for m6A methylation, Mettl3, was downregulation of Nxph3, encoding neurexophilin 3. The METTL3/METTL14 activator CHMA1004, in contrast, increased the expression of Nxph3 in the present study. Neurexophilin 3 is a ligand of synaptic α-neurexins thus contributing to the regulation of neurotransmitter release, and it has a role in sensorimotor gating and motor coordination (Beglopoulos et al., Reference Beglopoulos, Montag-Sallaz, Rohlmann, Piechotta, Ahmad, Montag and Missler2005).

The anxiolytic-like effects of a METTL3/METTL14 activating compound appear as inconsistent with findings that downregulation of Fto and increased m6A methylation can lead to anxiety- and depression-like behaviours in mice (Liu et al., Reference Liu, Xiu, Zhu, Meng, Li, Han, Du, Li, Xu, Liu, Zhu, Shen and Xu2021; Wang XL et al., Reference Wang, Wei, Yuan, Mao, Wang, Xing, Gu, Lin, Wang, Zhang and Xing2022). In these studies, the FTO-related anxious behaviours were associated with mRNA methylation and expression, respectively, of ß2-adrenoceptor in hippocampus or matrix metalloprotease 9 (MMP-9) in anterior cingulate cortex. In our study, neither Adrb2 nor MMP-9 expression was altered by CHMA1004 in striatum. While we did not examine the same brain regions, the contrasting behavioural findings suggest that METTL3/METTL14 activation by systemic administration of a small molecule compound affects RNA m6A methylation dissimilarly from downregulation of Fto by other type of interventions. Another study demonstrated that increased RNA m6A methylation of fatty acid amide hydrolase mRNA by reducing the function of ALKBH5 attenuated depressive-like behaviours (Huang et al., Reference Huang, Zhang, Bai, Han, Ju, Chen, Yang, Wang, Zhang, Zhang, Xie, Zhang and Yao2020), which is in line with our findings and further suggests that some of the multiple mechanisms recruited in modifying RNA m6A methylation have a potential in the management of psychiatric disorders. Furthermore, treadmill exercise was shown to counteract the anxiety-provoking effects of chronic stress in an RNA m6A methylation-dependent manner (Yan et al., Reference Yan, Wang, Yang, Wang, Wang, So and Zhang2023), providing support to the notion that METTL3/METTL14 activating drugs have anxiolytic and activating effects by promoting physiological mechanisms. As the next step, it is desirable to conduct transcriptomic analysis also in male rats, and RNA m6A sequencing to specify which transcriptomic changes correspond directly to drug effect at the epitranscriptomic levels.

Despite the study’s several limitations (small number of observations and limited statistical power due to different dose-response profile in males and females; potential carry-over effects in a test battery), the overall pattern of results suggested that treatment with CHMA1004 can elicit a variety of sizeable behavioural changes. Conclusively, the first-in-class METTL3/METTL14 activator compound methyl-piperazine-2-carboxylate (CHMA1004) elicited anxiolytic-like effects after acute systemic administration, and its anxiolytic-like and activating properties were consistently observed in a battery of behavioural tests during subchronic treatment. While behaviourally activating, CHMA1004 did not exhibit strong psychostimulant-like properties. Subchronic administration of CHMA1004 increased locomotor activity and led to alterations in catecholaminergic neurotransmission and changes in striatal gene expression linked to dopamine neuron viability, neurodegeneration, depression, anxiety and stress response. Thus, pharmacological activation of RNA m6A methylation has potential for neuropsychiatric drug development.

Supplementary material

The supplementary material for this article can be found at https://doi.org/10.1017/neu.2024.36.

Author contribution

MKk, MKn, MS and JH designed the study; MK, KL, MŠ, IT, KL and JH collected the data; MK, KL, IT, TO and KL analysed the data, all authors revised paper critically for important intellectual content.

Financial statement

This study was funded by the Estonian Research Council (PRG1213), the TLU Centre of Excellence in Behavioural and Neural Sciences and the School of Natural Sciences and Health Research Fund. JH was supported by Visiting Professorship of the Sigrid Jusélius Foundation. TO was supported by the European Union through Horizon 2020 research and innovation programme under grant number 810645. IT was funded by EU Marie Sklodowska Curie postdoctoral fellowship 101068830. MK and MS were supported from the Cure Parkinson’s Trust, UK grant (project code MS02).

Competing interests

MKn and JH are in the advisory board of Chemestmed Ltd that develops compounds that affect RNA m6A methylation. They are not paid for this role. Chemestmed did not provide support to the presented study except making the compound freely available.

JH has received lecturer fees from Lundbeck.

The authors declare no competing interests.

References

Ables, JL, Park, K and Ibañez-Tallon, I (2023) Understanding the habenula: a major node in circuits regulating emotion and motivation. Pharmacological Research 190, 106734. DOI: 10.1016/j.phrs.2023.106734.CrossRefGoogle Scholar
Armario, A, Gavaldà, A and Martí, O (1988) Forced swimming test in rats: effect of desipramine administration and the period of exposure to the test on struggling behavior, swimming, immobility and defecation rate. European Journal of Pharmacology 158, 207212. DOI: 10.1016/0014-2999(88)90068-4.CrossRefGoogle Scholar
Beglopoulos, V, Montag-Sallaz, M, Rohlmann, A, Piechotta, K, Ahmad, M, Montag, D and Missler, M (2005) Neurexophilin 3 is highly localized in cortical and cerebellar regions and is functionally important for sensorimotor gating and motor coordination. Molecular and Cellular Biology 25, 72787288. DOI: 10.1128/MCB.25.16.7278-7288.2005.CrossRefGoogle ScholarPubMed
Biever, A, Glock, C, Tushev, G, Ciirdaeva, E, Dalmay, T, Langer, JD and Schuman, EM (2020) Monosomes actively translate synaptic mRNAs in neuronal processes. Science 367, eaay4991. DOI: 10.1126/science.aay4991.CrossRefGoogle ScholarPubMed
Black, PH (2002) Stress and the inflammatory response: a review of neurogenic inflammation. Brain Behavior, and Immunity 16 622653. DOI: 10.1016/s0889-1591(02)00021-1.CrossRefGoogle ScholarPubMed
Brudzynski, SM and Granon, S (2023) Principles of recording rodent ultrasonic vocalizations,. In Harro, J (ed), Psychiatric Vulnerability, Mood, and Anxiety Disorders: Tests and Models in Mice and Rats, Neuromethods, 190, New York: Humana Press, pp. 329352.CrossRefGoogle Scholar
Burgdorf, J, Knutson, B, Panksepp, J and Ikemoto, S (2001) Nucleus accumbens amphetamine microinjections unconditionally elicit 50-kHz ultrasonic vocalizations in rats. Behavioral Neuroscience 115, 940944. DOI: 10.1037//0735-7044.115.4.940.CrossRefGoogle ScholarPubMed
Chang, Y, Zhang, J, Zhang, J, Zhu, W, Zheng, Q, Qian, Z, Wei, C, Liu, Y, Liu, Z, Ren, W and Han, J (2022) N6-methyladenosine RNA modification of glutamatergic neurons is associated with contextual fear discrimination. Physiology & Behavior 248, 113741. DOI: 10.1016/j.physbeh.2022.CrossRefGoogle ScholarPubMed
Chen, Q, Hartman, CA, Haavik, J, Harro, J, Klungsøyr, K, Hegvik, TA, Wanders, R, Ottosen, C, Dalsgaard, S, Faraone, SV and Larsson, H (2018) Common psychiatric and metabolic comorbidity of adult attention-deficit/hyperactivity disorder: a population-based cross-sectional study. Plos One 13, 204516. DOI: 10.1371/journal.pone 0. 204516.Google ScholarPubMed
Chen, S, Zhou, Y, Chen, Y and Gu, J (2018) Fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 34, i884890.CrossRefGoogle ScholarPubMed
Chen, X, Yu, C, Guo, M, Zheng, X, Ali, S, Huang, H, Zhang, L, Wang, S, Huang, Y, Qie, S and Wang, J (2019) Down-regulation of m6A mRNA methylation is involved in dopaminergic neuronal death. ACS Chemical Neuroscience 10, 23552363. DOI: 10.1021/acschemneuro.8b00657.CrossRefGoogle ScholarPubMed
Dai, D, Yuan, J, Wang, Y, Xu, J, Mao, C and Xiao, Y (2019) Peli1 controls the survival of dopaminergic neurons through modulating microglia-mediated neuroinflammation. Scientific Reports. 9, 8034. DOI: 10.1038/s41598-019-44573-w.CrossRefGoogle ScholarPubMed
Dalley, JW and Robbins, TW (2017) Fractionating impulsivity: neuropsychiatric implications. Nature Reviews Neuroscience 18, 158171. DOI: 10.1038/nrn.2017.8.CrossRefGoogle ScholarPubMed
De Jager, PL, Srivastava, G, Lunnon, K, Burgess, J, Schalkwyk, LC, Yu, L, Eaton, ML, Keenan, BT, Ernst, J, McCabe, C, Tang, A, Rai, T, Replogle, J, Brodeur, W, Gabriel, S, Chai, HS, Younkin, C, Younkin, SG, Zou, F, Szyf, M, Epstein, CB, Schneider, JA, Bernstein, BE, Meissner, A, Ertekin-Taner, N, Chibnik, LB, Kellis, M, Mill, J and Bennett, DA (2014) Alzheimer’s disease: early alterations in brain DNA methylation at ANK1, BIN1, RHBDF2 and other loci. Nature Neuroscience 17, 11561163. DOI: 10.1038/nn.3786.CrossRefGoogle ScholarPubMed
Du, T, Rao, S, Wu, L, Ye, N, Liu, Z, Hu, H, Xiu, J, Shen, Y and Xu, Q (2015) An association study of the m6A genes with major depressive disorder in chinese han population. Journal of Affective Disorders 183, 279286. DOI: 10.1016/j.jad.2015.05.025.CrossRefGoogle ScholarPubMed
Engel, M, Eggert, C, Kaplick, PM, Eder, M, Röh, S, Tietze, L, Namendorf, C, Arloth, J, Weber, P, Rex-Haffner, M, Geula, S, Jakovcevski, M, Hanna, JH, Leshkowitz, D, Uhr, M, Wotjak, CT, Schmidt, MV, Deussing, JM, Binder, EB and Chen, A (2018) The role of m6A/mRNA methylation in stress response regulation. Neuron 99, 389403.e9.e389.CrossRefGoogle Scholar
Fakira, AK, Lueptow, LM, Trimbake, NA and Devi, LA (2021) PEN receptor GPR83 in anxiety-like behaviors: differential regulation in global vs amygdalar knockdown. Frontiers in Neuroscience 15, 675769. DOI: 10.3389/fnins.2021.675769.CrossRefGoogle ScholarPubMed
File, SE and Hyde, JR (1978) Can social interaction be used to measure anxiety? British Journal of Pharmacology 62, 1924.CrossRefGoogle ScholarPubMed
File, SE and Hyde, JR (1979) A test of anxiety that distinguishes between the actions of benzodiazepines and those of other minor tranquillizers and of stimulants. Pharmacology Biochemistry & Behavior 11, 6569. DOI: 10.1016/0091-3057(79)90298-3.CrossRefGoogle ScholarPubMed
Fustin, JM, Doi, M, Yamaguchi, Y, Hida, H, Nishimura, S, Yoshida, M, Isagawa, T, Morioka, MS, Kakeya, H, Manabe, I and Okamura, H (2013) RNA-methylation-dependent RNA processing controls the speed of the circadian clock.. Cell 155:793806 DOI: 10.1016/j.cell.2013.10.026.CrossRefGoogle Scholar
Grochowska, MM, Carreras Mascaro, A, Boumeester, V, Natale, D, Breedveld, GJ, Geut, H, van Cappallen, WA, Boon, AJW, Kievit, AJA, Sammler, E, Bank, Netherlands Brain, Parchi, P, Cortelli, P, Alessi, DR, van de Berg, WDJ, Bonifati, V and Mandemakers, W (2021) LRP10 interacts with SORL1 in the intracellular vesicle trafficking pathway in non-neuronal brain cells and localises to Lewy bodies in parkinson’s disease and dementia with Lewy bodies. Acta Neuropathologica 142, 117137. DOI: 10.1007/s00401-021-02313-3.CrossRefGoogle ScholarPubMed
Guo, L, Cao, J, Hou, J, Li, Y, Huang, M, Zhu, L, Zhang, L, Lee, Y, Lemos Duarte, M, Zhou, X, Wang, M, Liu, CC, Martens, Y, Chao, M, Goate, A, Bu, G, Haroutunian, V, Cai, D and Zhang, B (2023) Sex specific molecular networks and key drivers of alzheimer’s disease. Molecular Neurodegeneration 18, 39. DOI: 10.1186/s13024-023-00624-5.CrossRefGoogle ScholarPubMed
Han, S, DiBlasi, E, Monson, ET, Shabalin, A, Ferris, E, Chen, D, Fraser, A, Yu, Z, Staley, M, Callor, WB, Christensen, EB, Crockett, DK, Li, QS, Willour, V, Bakian, AV, Keeshin, B, Docherty, AR, Eilbeck, K and Coon, H (2023) Whole-genome sequencing analysis of suicide deaths integrating brain-regulatory eQTLs data to identify risk loci and genes. Molecular Psychiatry 28,39093919. DOI: 10.1038/s41380-023-02282-x.CrossRefGoogle ScholarPubMed
Harro, J (2019) Animal models of depression: pros and cons. Cell and Tissue Research 377, 520. DOI: 10.1007/s00441-018-2973-0.CrossRefGoogle ScholarPubMed
Harro, J, Tõnissaar, M, Eller, M, Kask, A and Oreland, L (2001) Chronic variable stress and partial 5-HT denervation by parachloroamphetamine treatment in the rat: effects on behavior and monoamine neurochemistry. Brain Research 899, 227239. DOI: 10.1016/s0006-8993(01)02256-9.CrossRefGoogle ScholarPubMed
Harvanek, ZM, Boks, MP, Vinkers, CH and Higgins-Chen, AT (2023) The cutting edge of epigenetic clocks: In search of mechanisms linking aging and mental health. Biological Psychiatry 94, 694705. DOI: 10.1016/j.biopsych.2023.02.001.CrossRefGoogle ScholarPubMed
Huang, R, Zhang, Y, Bai, Y, Han, B, Ju, M, Chen, B, Yang, L, Wang, Y, Zhang, H, Zhang, H, Xie, C, Zhang, Z and Yao, H (2020) N6-methyladenosine modification of fatty amide hydrolase messenger RNA in circular RNA STAG1-regulated astrocyte dysfunction and depressive-like behaviors. Biological Psychiatry 88, 392404. DOI: 10.1016/j.biopsych.2020.02.018.CrossRefGoogle ScholarPubMed
Juruena, MF, Gadelrab, R, Cleare, AJ and Young, AH (2021) Epigenetics: a missing link between early life stress and depression. Progress in Neuro-Psychopharmacology & Biological Psychiatry 109, 110231. DOI: 10.1016/j.pnpbp.2020.CrossRefGoogle ScholarPubMed
Kõiv, K, Matrov, D, Uusen, T and Harro, J (2021) Effect of neuropeptide S administration on ultrasonic vocalizations and behaviour in rats with low vs. high exploratory activity. Pharmaceuticals (Basel) 14, 524. DOI: 10.3390/ph14060524.CrossRefGoogle ScholarPubMed
Kõiv, K, Metelitsa, M, Vares, M, Tiitsaar, K, Raudkivi, K, Jaako, K, Vulla, K, Shimmo, R and Harro, J (2016) Chronic variable stress prevents amphetamine-elicited 50-kHz calls in rats with low positive affectivity. European Neuropsychopharmacology 26(4), 631643. DOI: 10.1016/j.euroneuro.2016.02.011.CrossRefGoogle ScholarPubMed
Kõiv, K, Vares, M, Kroon, C, Metelitsa, M, Tiitsaar, K, Laugus, K, Jaako, K and Harro, J (2019) Effect of chronic variable stress on sensitization to amphetamine in high and low sucrose-consuming rats. Journal of Psychopharmacology 33, 15121523. DOI: 10.1177/0269881119856000.CrossRefGoogle ScholarPubMed
Koopmans, F, van Nierop, P, Andres-Alonso, M, Byrnes, A, Cijsouw, T, Coba, MP and et al. (2019) SynGO: An evidence-based, expert-curated knowledge base for the synapse. Neuron 103(2), 217234.e4.CrossRefGoogle ScholarPubMed
Li, C, Xu, P, Huang, Y, Wang, Y, Wu, Y, Li, H, Peng, Y, Wang, X, Li, G, He, J and Miao, C (2023) RNA methylations in depression, from pathological mechanisms to therapeutic potential. Biochemical Pharmacology 215, 115750.CrossRefGoogle ScholarPubMed
Li, Y, Ren, J, Zhang, Z, Weng, Y, Zhang, J, Zou, X, Wu, S and Hu, H (2023) Modification and expression of mRNA m6A in the lateral habenula of rats after long-term exposure to blue light during the sleep period. Genes (Basel) 14, 143. DOI: 10.3390/genes14010143.CrossRefGoogle ScholarPubMed
Li, YJ, Kresock, E, Kuplicki, R, Savitz, J and McKinney, BA (2022) Differential expression of MDGA1 in major depressive disorder. Brain, Behavior, & Immunity - Health 26, 100534. DOI: 10.1016/j.bbih.2022.100534.CrossRefGoogle ScholarPubMed
Liao, J, Mi, X, Zeng, G, Wei, Y, Dai, X, Ye, Q, Chen, X and Zhang, J (2023) Circuit-wide proteomics profiling reveals brain region-specific protein signatures in the male WKY rats with endogeneous depression. Journal of Affective Disorders 320, 98107. DOI: 10.1016/j.jad.2022.09.086.CrossRefGoogle ScholarPubMed
Liao, Y, Smyth, GK and Shi, W (2014) FeatureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 30, 923930.CrossRefGoogle ScholarPubMed
Liiver K, Imbeault S, Školnaja M, Kaart T, Kanarik M, Laugus K, De Wettinck J, Pulver A, Shimmo R and Harro J (2023) Active vs passive novelty-related strategies: sex differences in exploratory behaviour and monoaminergic systems. Behavioural Brain Research 441, 1 14297. DOI: 10.1016/j.bbr.2023.114297.Google Scholar
Liu, S, Xiu, J, Zhu, C, Meng, K, Li, C, Han, R, Du, T, Li, L, Xu, L, Liu, R, Zhu, W, Shen, Y and Xu, Q (2021) Fat mass and obesity-associated protein regulates RNA methylation associated with depression-like behavior in mice. Nature Communications. 12, 6937. DOI: 10.1038/s41467-021-27044-7.CrossRefGoogle ScholarPubMed
Livneh, I, Moshitch-Moshkovitz, S, Amariglio, N, Rechavi, G and Dominissini, D (2020) The m6a epitranscriptome: transcriptome plasticity in brain development and function. Nature Reviews Neuroscience 21, 3651. DOI: 10.1038/s41583-019-0244-z.CrossRefGoogle Scholar
Love, MI, Huber, W and Anders, S (2014) Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biology 15, 550.CrossRefGoogle ScholarPubMed
Mällo, T, Matrov, D, Kõiv, K and Harro, J (2009) Effect of chronic stress on behavior and cerebral oxidative metabolism in rats with high or low positive affect. Neuroscience 164, 963974. DOI: 10.1016/j.neuroscience.2009.08.041.CrossRefGoogle ScholarPubMed
Malovic, E and Pandey, SC (2023) N6-methyadenosine (m6A) epitranscriptomics in synaptic plasticity and behaviors. Neuropsychopharmacology 48, 221222.CrossRefGoogle Scholar
Marie-Claire, C, Lejeune, FX, Mundwiller, E, Ulveling, D, Moszer, I, Bellivier, F and Etain, B (2020) A DNA methylation signature discriminates between excellent and non-response to lithium in patients with bipolar disorder type 1. Scientific Reports. 10, 12239. DOI: 10.1038/s41598-020-69073-0.CrossRefGoogle ScholarPubMed
Martin, FJ, Amode, MR, Aneja, A, Austine-Orimoloye, O, Azov, AG, Barnes, I, et al. (2023) Ensembl 2023. Nucleic Acids Research. 51, D933D941.CrossRefGoogle ScholarPubMed
Matsuda, W, Furuta, T, Nakamura, KC, Hioki, H, Fujiyama, F, Arai, R and Kaneko, T (2009) Single nigrostriatal dopaminergic neurons form widely spread and highly dense axonal arborizations in the neostriatum. Journal of Neuroscience 29, 444453. DOI: 10.1523/JNEUROSCI.4029-08.2009.CrossRefGoogle ScholarPubMed
Matsunaga, E, Nambu, S, Oka, M, Tanaka, M, Taoka, M and Iriki, A (2015) Identification of tool use acquisition-associated genes in the primate neocortex. Development, Growth & Differentiation 57, 484495. DOI: 10.1111/dgd.CrossRefGoogle ScholarPubMed
Matsuzaki, T, Yoshihara, T, Ohtsuka, T and Kageyama, R (2019) Hes1 expression in mature neurons in the adult mouse is required for normal behaviors. Scientific Reports 9, 8251. DOI: 10.1038/s41598-019-44698-y.CrossRefGoogle ScholarPubMed
Matto, V, Harro, J and Allikmets, L (1997) The effects of cholecystokinin a and B receptor antagonists on exploratory behaviour in the elevated zero-maze in rat. Neuropharmacology 36, 389396. DOI: 10.1016/S0028-3908(97)00011-7.CrossRefGoogle Scholar
Meyer, KD and Jaffrey, SR (2017) Rethinking m6a readers, writers, and erasers. Annual Review of Cell and Developmental Biology 33, 319342. DOI: 10.1146/annurev-cellbio-100616-060758.CrossRefGoogle ScholarPubMed
O’Leary, A, Kõiv, K, Raudkivi, K and Harro, J (2016) Antidepressants differentially affect striatal amphetamine-stimulated dopamine and serotonin release in rats with high and low novelty-oriented behaviour. Pharmacological Research 113, 739746.CrossRefGoogle ScholarPubMed
Panariello, F, Fanelli, G, Fabbri, C, Atti, AR, De Ronchi, D and Serretti, A (2022) Epigenetic basis of psychiatric disorders: a narrative review. CNS & Neurological Disorders - Drug Targets 21, 302315. DOI: 10.2174/1871527320666210825101915.CrossRefGoogle ScholarPubMed
Panksepp, J (1998) Affective Neuroscience: The Foundations of Human and Animal Emotions. Oxford: Oxford University Press.CrossRefGoogle Scholar
Paxinos, G and Watson, C (1986) The Rat Brain Stereotaxic Coordinates. San Diego: Academic Press.Google Scholar
Persaud, NS and Cates, HM (2023) The epigenetics of anxiety pathiphysiology: a DNA methylation and histone modification focused review. eNeuro 10, ENEURO.0109–21.2021. DOI: 10.1523/ENEURO.0109-21.2021.CrossRefGoogle ScholarPubMed
Pishva, E, van den Hove, DLA, Laroche, V, Lvovs, A, Roy, A, Ortega, G, Burrage, J, Veidebaum, T, Kanarik, M, Mill, J, Lesch, KP and Harro, J (2023) Genome-wide DNA methylation analysis of aggressive behaviour: a longitudinal population-based study. Journal of Child Psychology and Psychiatry 64, 9981006. DOI: 10.1111/jcpp.13782.CrossRefGoogle ScholarPubMed
Porsolt, RD, Anton, G, Blavet, N and Jalfre, M (1978) Behavioural despair in rats: a new model sensitive to antidepressant treatments. European Journal of Pharmacology 47, 379391. DOI: 10.1016/0014-2999(78)90118-8.CrossRefGoogle ScholarPubMed
Qi, X, Wang, S, Zhang, L, Liu, L, Wen, Y, Ma, M, Cheng, S, Li, P, Cheng, B, Du, Y, Liang, X, Zhao, Y, Ding, M and Zhang, F (2019) An integrative analysis of transcriptome-wide association study and mRNA expression profile identified candidate genes for attention-deficit/hyperactivity disorder. Psychiatry Research 282, 112639.CrossRefGoogle ScholarPubMed
Rehm, J, Anderson, P, Manthey, J, Shield, KD, Struzzo, P, Wojnar, M and Gual, A (2016) Alcohol use disorders in primary health care: what do we know and where do we go? Alcohol and Alcoholism 51, 422427. DOI: 10.1093/alcalc/agv127.CrossRefGoogle ScholarPubMed
Rizzi, A, Vergura, R, Marzola, G, Ruzza, C, Guerrini, R, Salvadori, S, Regoli, D and Calo, G (2008) Neuropeptide S is a stimulatory anxiolytic agent: a behavioural study in mice. British Journal of Pharmacology 154, 471479. DOI: 10.1038/bjp.2008.96.CrossRefGoogle ScholarPubMed
Roy, B, Ochi, S and Dwivedi, Y (2022) M6a RNA methylation-based epitranscriptomic modifications in plasticity-related genes via miR-124-C/EBPα-FTO-transcriptional axis in the hippocampus of learned helplessness rats. International Journal of Neuropsychopharmacology 25, 10371049.CrossRefGoogle ScholarPubMed
Samaan, Z, Anand, SS, Zhang, X, Desai, D, Rivera, M, Pare, G, Thabane, L, Xie, C, Gerstein, H, Engert, JC, Craig, I, Cohen-Woods, S, Mohan, V, Diaz, R, Wang, X, Liu, L, Corre, T, Preisig, M, Kutalik, Z, Bergmann, S, Vollenweider, P, Waeber, G, Yusuf, S and Meyre, D (2013) The protective effect of the obesity-associated rs9939609 a variant in fat mass- and obesity-related gene on depression. Molecular Psychiatry 18, 12811286. DOI: 10.1038/mp.2012.160.CrossRefGoogle Scholar
Selberg, S, Blokhina, D, Aatonen, M, Koivisto, P, Siltanen, A, Mervaala, E, Kankuri, E and Karelson, M (2019) Discovery of small molecules that activate RNA methylation through cooperative binding to the METTL3-14-WTAP complex active site. Cell Reports 26, 37623771.e5. DOI: 10.1016/j.celrep.2019.02.100.CrossRefGoogle Scholar
Selberg, S, Yu, LY, Bondarenko, O, Kankuri, E, Seli, N, Kovaleva, V, Herodes, K, Saarma, M and Karelson, M (2021b) Small-molecule inhibitors of the RNA m6A demethylases FTO potently support the survival of dopamine neurons. International Journal of Molecular Science 22, 4537. DOI: 10.3390/ijms22094537.CrossRefGoogle ScholarPubMed
Selberg, S, Žusinaite, E, Herodes, K, Seli, N, Kankuri, E, Merits, A and Karelson, M (2021a) HIV replication is increased by RNA methylation METTL3/METTL14/WTAP compex activators. ACS Omega 6, 1595715963. DOI: 10.1021/acsomega 1626.CrossRefGoogle Scholar
Shafik, AM, Allen, EG and Jin, P (2022) Epitranscriptomic dynamics in brain development and disease. Molecular Psychiatry 27, 36333646. DOI: 10.1038/s41380-022-01570-2.CrossRefGoogle ScholarPubMed
Shani, S, Goldstein, O, Gana-Weisz, M, Bar-Shira, A, Thaler, A, Gurevich, T, Mirelman, A, Giladi, N, Alcalay, RN and Orr-Urtreger, A (2023) Variants in PSMB9 and FGR differentially affect parkinson’s disease risk in GBA and LRRK2 mutation carriers. Parkinsonism and Related Disorders 111, 105398. DOI: 10.1016/j.parkreldis.2023.105398.CrossRefGoogle ScholarPubMed
Shepherd, JK, Grewal, SS, Fletcher, A, Bill, DJ and Dourish, CT (1994) Behavioural and pharmacological characterisation of the elevated, zero-maze, as an animal model of anxiety. Psychopharmacology (Berl) 116, 5664. DOI: 10.1007/BF02244871.CrossRefGoogle ScholarPubMed
Sherman, BT, Hao, M, Qiu, J, Jiao, X, Baseler, MW, Lane, HC, Imamichi, T and Chang, W (2022) DAVID: a web server for functional enrichment analysis and functional annotation of gene lists (2021 update). Nucleic Acids Research 50, W216W221.CrossRefGoogle Scholar
Simola, N, Ma, ST and Schallert, T (2010) Influence of acute caffeine on 50-kHz ultrasonic vocalizations in male adult rats and relevance to caffeine-mediated psychopharmacological effects. International Journal of Neuropsychopharmacology 13, 123132. DOI: 10.1017/S1461145709990113.CrossRefGoogle ScholarPubMed
Su, X, Qu, Y and Mu, D (2023) The regulatory network of METTL3 in the nervous system: diagnostic biomarkers and therapeutic targets. Biomolecules 13,664. DOI: 10.3390/biom13040664.CrossRefGoogle ScholarPubMed
Thalmeier, A, Dickmann, M, Giegling, I, Schneider, BM, Hartmann, A, Maurer, K, Schnabel, A, Kauert, G, Möller, HJ and Rujescu, D (2008) Gene expression profiling of post-mortem orbitofrontal cortex in violent suicide victims. International Journal of Neuropsychopharmacology 11, 217228. DOI: 10.1017/S1461145707007894.CrossRefGoogle ScholarPubMed
Tõnissaar, M, Mällo, T, Eller, M, Häidkind, R, Kõiv, K and Harro, J (2008) Rat behavior after chronic variable stress and partial lesioning of 5-HT-ergic neurotransmission: effects of citalopram. Progress in Neuro-Psychopharmacology & Biological Psychiatry 32, 164177. DOI: 10.1016/j.pnpbp.2007.08.001.CrossRefGoogle ScholarPubMed
Tõnissaar, M, Philips, MA, Eller, M and Harro, J (2004) Sociability trait and serotonin metabolism in the rat social interaction test. Neuroscience Letters 367, 309312. DOI: 10.1016/j.neulet.2004.06.023.CrossRefGoogle ScholarPubMed
Walters, BJ, Mercaldo, V, Gillon, CJ, Yip, M, Neve, RL, Boyce, FM, Frankland, PW and Josselyn, SA (2017) The role of the RNA demethylase FTO (fat mass and obesity-associated) and mRNA methylation in hippocampal memory formation. Neuropsychopharmacology 42, 15021510. DOI: 10.1038/npp.2017.31.CrossRefGoogle ScholarPubMed
Wang, P, Doxtader, KA and Nam, Y (2016) Structural basis for cooperative function of Mettl3 and Mettl14 methyltransferases. Molecular Cell 63, 306317.CrossRefGoogle ScholarPubMed
Wang, XL, Wei, X, Yuan, JJ, Mao, YY, Wang, ZY, Xing, N, Gu, HW, Lin, CH, Wang, WT, Zhang, W and Xing, F (2022) Downregulation of fat mass and obesity-related protein in the anterior cingulate cortex participates in anxiety- and depression-like behaviors induced by neuropathic pain. Frontiers in Cellular Neuroscience 16, 884296. DOI: 10.3389/fncel.2022.884296.Google ScholarPubMed
Wang, XY and Zhang, LN (2023) RNA binding protein SAMD4: current knowledge and future perspectives. BMC Cell & Bioscience 13, 21.CrossRefGoogle ScholarPubMed
Wettergren, EE, Gussing, F, Quintino, L and Lundberg, C (2012) Novel disease-specific promoters for use in gene therapy for parkinson’s disease. Neuroscience Letters 530, 2934. DOI: 10.1016/j.neulet.2012.09.059.CrossRefGoogle ScholarPubMed
Widagdo, J, Zhao, Q, Kempen, MJ, Tan, MC, Ratnu, VS, Wei, W, Leighton, L, Spadaro, PA, Edson, J, Anggono, V and Bredy, TW (2016) Experience-dependent accumulation of N6-methyadenosine in the prefrontal cortex is associated with memory processes in mice. Journal of Neuroscience 36, 67716777. DOI: 10.1523/JNEUROSCI.4053-15.2016.CrossRefGoogle Scholar
Widjaya, MA, Cheng, YJ, Kuo, YM, Liu, CH, Cheng, WC and Lee, SD (2023) Transcriptomic analysis of exercise training in alzheimer’s disease cerebral cortex. Journal of Alzheimers Disease 93, 349363. DOI: 10.3233/JAD-221139.CrossRefGoogle ScholarPubMed
Wöhr, M (2022) Measuring mania-like elevated mood through amphetamine-induce 50-kHz ultrasonic vocalizations in rats. British Journal of Pharmacology 179, 42014219. DOI: 10.1111/bph.15487.CrossRefGoogle ScholarPubMed
Xu, W, Yang, T, Lou, X, Chen, J, Wang, X, Hu, M, An, D, Gao, R, Wang, J and Chen, X (2023) Role of the Peli1-RIPK1 signaling axis in methamphetamine-induced neuroinflammation. ACS Chemical Neuroscience 14, 864874. DOI: 10.1021/acschemneuro.2c00623.CrossRefGoogle ScholarPubMed
Xu, YLL, Reinscheid, RK, Huitron-Resendiz, S, Clark, SD, Wang, Z, Lin, SH, Brucher, FA, Zeng, J, Ly, NK, de Lecea, LHenriksen SJ and Civelli, O (2004) Neuropeptide S: a neuropeptide promoting arousal and anxiolytic-like effects. Neuron 43, 487497. DOI: 10.1016/j.neuron.2004.08.005.CrossRefGoogle Scholar
Yan, L, Wang, M, Yang, F, Wang, Y, Wang, S, So, KF and Zhang, L (2023) Physical exercise mediates a cortical FMRP-mTOR pathway to improve resilience against chronic stress in adolescent mice. Translational Psychiatry 13, 16. DOI: 10.1038/s41398-023-02311-x.CrossRefGoogle ScholarPubMed
Yu, LY, Selberg, S, Teino, I, Nam, J, Ivanova, I, Renzi, R, Seli, N, Kankuri, E, Voutilainen, MH, Karelson, M and Saarma, M (2023) Small molecule activation of m6A mRNA methylation as a novel approach for neuroprotection. bioRxiv. DOI: 10.1101/2023.07.05.547860.Google Scholar
Zaccara, S, Ries, RJ and Jaffrey, SR (2019) Reading, writing and erasing mRNA methylation. Nature Reviews Molecular Cell Biology 20, 608624. DOI: 10.1038/s41580-019-0168-5.CrossRefGoogle ScholarPubMed
Figure 0

Fig. 1. Experimental design of behavioural and ex vivo biochemical studies of the effect of acute and subchronic treatment with CHMA1004. Note: HPLC, high-performance liquid chromatography; IP, intraperitoneal; RNA-seq, RNA sequencing; USVs, recording of ultrasonic vocalisation.

Figure 1

Fig. 2. Effect of acute administration of CHMA1004 on locomotion, defecation and emission of ultrasonic vocalisations in an open field test. Note: *p < 0.05; **p < 0.01 versus vehicle. Post hoc test on rearings in male rats CHMA1004 (5 mg/kg) p = 0.069.

Figure 2

Fig. 3. Effect of repeated administration of CHMA1004 in elevated zero-maze test. Note: SAPs, number of stretched-attend postures; *p < 0.05 versus vehicle.

Figure 3

Table 1. Monoamine levels after two and half weeks of administration of the METTL3/METTL14 activating compound CHMA1004 in male and female rats

Figure 4

Fig. 4. Effect of repeated administration of CHMA1004 in the forced swimming test. Note: *p < 0.05 versus vehicle.

Figure 5

Fig. 5. Extracellular levels of dopamine after acute administration of amphetamine or CHMA1004. Note: A: female rats; B: male rats.

Figure 6

Fig. 6. Effect of repeated administration of CHMA1004 on striatal gene expression in female rats. Note: Volcano plots of the differentially expressed gene analysis with DESeq2 between control and CHMA1004 1 mg/kg (left) or 5 mg/kg (right). Significantly differentially expressed genes (p < 0.05) are labelled blue (downregulated genes) or red (upregulated genes). The names of the top 10 genes with the lowest p-values are shown.

Supplementary material: File

Kanarik et al. supplementary material 1

Kanarik et al. supplementary material
Download Kanarik et al. supplementary material 1(File)
File 240.4 KB
Supplementary material: File

Kanarik et al. supplementary material 2

Kanarik et al. supplementary material
Download Kanarik et al. supplementary material 2(File)
File 62.8 KB
Supplementary material: File

Kanarik et al. supplementary material 3

Kanarik et al. supplementary material
Download Kanarik et al. supplementary material 3(File)
File 252 KB
Supplementary material: File

Kanarik et al. supplementary material 4

Kanarik et al. supplementary material
Download Kanarik et al. supplementary material 4(File)
File 64.9 KB
Supplementary material: File

Kanarik et al. supplementary material 5

Kanarik et al. supplementary material
Download Kanarik et al. supplementary material 5(File)
File 63.5 KB
Supplementary material: File

Kanarik et al. supplementary material 6

Kanarik et al. supplementary material
Download Kanarik et al. supplementary material 6(File)
File 19.8 KB