# SNORA69 is up-regulated in the lateral habenula of individuals with major depressive disorder

Rixing Lin, Haruka Mitsuhashi, Laura M Fiori, Ryan Denniston, El Cherif

Ibrahim, Catherine Belzung, Naguib Mechawar, Gustavo Turecki

► To cite this version:

Rixing Lin, Haruka Mitsuhashi, Laura M Fiori, Ryan Denniston, El Cherif Ibrahim, et al.. SNORA69 is up-regulated in the lateral habenula of individuals with major depressive disorder. Scientific Reports, 2024, 14, pp.8258. 10.1038/s41598-024-58278-2 . hal-04672660

HAL Id: hal-04672660

https://hal.science/hal-04672660v1

Submitted on 19 Aug 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
<figure><img src="image_2.png" /></figure>

OPEN

# SNORA69 is up-regulated in the lateral habenula of individuals with major depressive disorder

Rixing Lin<fnref n="1" />, Haruka Mitsuhashi<fnref n="1" />, Laura M. Fiori<fnref n="1" />, Ryan Denniston<fnref n="1" />, El Cherif Ibrahim², Catherine Belzung³, Naguib Mechawar<fnref n="1" /> & Gustavo Turecki<fnref n="1" />✉

Major depressive disorder (MDD) is a complex and potentially debilitating illness whose etiology and pathology remains unclear. Non-coding RNAs have been implicated in MDD, where they display differential expression in the brain and the periphery. In this study, we quantified small nucleolar RNA (snoRNA) expression by small RNA sequencing in the lateral habenula (LHb) of individuals with MDD (n = 15) and psychiatrically-healthy controls (n = 15). We uncovered five snoRNAs that exhibited differential expression between MDD and controls (FDR < 0.01). Specifically, *SNORA69* showed increased expression in MDD and was technically validated via RT-qPCR. We further investigated the expression of *Snora69* in the LHb and peripheral blood of an unpredicted chronic mild stress (UCMS) mouse model of depression. *Snora69* was specifically up-regulated in mice that underwent the UCMS paradigm. *SNORA69* is known to guide pseudouridylation onto 5.8S and 18S rRNAs. We quantified the relative abundance of pseudouridines on 5.8S and 18S rRNA in human post-mortem LHb samples and found increased abundance of pseudouridines in the MDD group. Overall, our findings indicate the importance of brain snoRNAs in the pathology of MDD. Future studies characterizing *SNORA69*'s role in MDD pathology is warranted.

Major depressive disorder (MDD) is a prevalent illness and is the leading cause of disability and premature death worldwide<fnref n="1" />⁻³. Despite the prevalence and burden of MDD the underlying etiology and pathological processes that occur in the brain of depressed individuals are still not well characterized. The lateral habenula (LHb) is a key brain structure mediating the communication between the forebrain and midbrain/hindbrain, thus, regulating monoaminergic neurotransmission⁴⁻⁶. It is also involved in pathways related to reward processing, stress adaptation, sleep/circadian rhythm, and its dysfunction is linked to maladaptive processing of negative and positive valence; all of which are core behaviors dysregulated in MDD⁴⁻⁶. While many molecular factors have been associated with MDD, the primary focus of past molecular studies has been the investigation of gene expression and systems that regulate gene expression⁷⁻⁹. More recently, numerous studies have highlighted the role of various non-coding RNAs in MDD, or antidepressant response, including microRNAs, long non-coding RNAs and circular RNAs and more recently small nucleolar RNAs¹⁰⁻¹³.

The term snoRNAs was coined in the early 1980s; however, the first identification of snoRNAs can be traced back to the early 1960s¹⁴⁻¹⁷. RNA sequencing has allowed the identification and annotation of over 200 snoRNAs in humans; however, it is estimated that there are over 1000 unique snoRNAs¹⁸,¹⁹. SnoRNAs are essential, medium length (60–300 nucleotides long), non-coding RNAs present in all eukaryotes. Conserved sequence motifs (“boxes”) present on all snoRNAs further separate snoRNAs into two main distinct classes. Box C/D snoRNAs contain box C (5'-RUGAUGA-3') and box D (5'-CUGA-3') sequence motifs located near the 5' and 3' end respectively. Box H/ACA snoRNAs contain box H (5'-ANANNA-3') and a 5'-ACA-3' sequence present near the 3' end¹⁹,²⁰. Box C/D and box H/ACA snoRNAs recruit distinct sets of partner proteins, forming a ribonucleoprotein (snoRNP) complex¹⁹,²⁰. Classically, snoRNAs function by guiding chemical modifications onto ribosomal RNA (rRNA) or small nuclear RNAs (snRNAs), with box C/D snoRNAs guiding 2'-oxygen-ribose methylation and box H/ACA snoRNAs guiding pseudouridylation¹⁹,²⁰.

<footnote n="1">¹Department of Psychiatry, McGill Group for Suicide Studies, Douglas Mental Health University Institute, McGill University, Montreal, QC, Canada. ²CNRS, INT, Institute Neuroscience Timone, Aix-Marseille Université, Marseille, France. ³Imaging Brain and Neuropsychiatry iBraiN U1253, INSERM, Université de Tours, Tours, France. ✉email: gustavo.turecki@mcgill.ca</footnote>
In this study we used small RNA-sequencing, in human post-mortem LHB, to perform an un-biased profile of snoRNAs that are differentially expressed between individuals with depression and psychiatrically healthy controls.

# Materials and methods

## Human samples

The human post-mortem LHB samples included in this study partially overlap with those included in a previous study from our group¹⁰. LHB samples were obtained, in collaboration with the Quebec Coroner's Office, from the Douglas-Bell Canada Brain Bank (Douglas Mental Health University Institute, Montreal, Quebec, Canada). All individuals were of French-Canadian origin, a homogeneous population with a well-documented founder effect²¹. We analyzed a total of 30 LHB samples (15 individuals with MDD, and 15 psychiatrically healthy controls) (Table 1). One sample in the control group and seven samples in the MDD group had presence of psychotropic drugs at the time of death. Individuals with MDD all died by suicide (hanging, intoxication, or jumping) while control subjects died by natural causes or by accident. Groups were matched for age, sex, post-mortem interval (PMI), RNA integrity number (RIN), and pH (Table 1). Psychological autopsies were performed as described previously, based on DSM-IV criteria²². Written informed consent was obtained from next-of-kin. This study and the methods described were approved by the Douglas Hospital Research Centre institutional review board and performed in accordance with relevant guidelines and regulations.

## Small RNA sequencing and analysis

RNA was extracted from all brain samples using a combination of the miRNeasy Mini kit and the RNeasy MinElute Cleanup kit (Qiagen), with DNase treatment, and divided into small (< 200 nt) and large (> 200 nt) fractions¹⁰. RNA quality was assessed using the Agilent 2200 Tapestation. Small RNA libraries were prepared from the small RNA fraction, using the Illumina TruSeq Small RNA protocol following the manufacturer's instructions. Samples were sequenced at the McGill University and Genome Quebec Innovation Centre (Montreal, Canada) using the Illumina HiSeq2000 with 50 nt single-end reads¹⁰. Adapters and low-quality bases (Phred score < 15) were trimmed with cutadapt v2.10. The excerpt pipeline v4.6.3 with the options STARAlignEndsType=EndToEnd, STAR_outFilterMismatchNmax=0, and STAR_outFilterMatchNmin=16 was used to align trimmed reads to the genome (GRCh38) and quantify the snoRNAs into counts. Only snoRNAs with greater than or equal to 10 counts in 70% of the samples within either the control or MDD group were retained for later steps. The svaseq function from the sva v3.46.0 was used to identify surrogate variables representing known or unknown factors significantly influencing our dataset for inclusion in the model for differential expression. DESeq2 v1.38.1 was used to normalize the counts with the median of ratios method and to perform differential expression analysis on our dataset.

## Validation

RNA was reverse transcribed to cDNA using M-MLV reverse transcriptase (Gibco) following the manufacturer's protocol with random hexamers. Real-time quantitative PCR (RT-qPCR) reactions were run in triplicates using the QuantStudio 6 Flex Real-Time PCR System and data were collected using the QuantStudio Real-Time PCR Software v1.1 (Applied Biosystems). To measure RNA expression, we used custom-designed probes (Supplementary Table 1) with PowerUp SYBR Green Master Mix (ThermoFisher). Expression levels were calculated using the absolute quantitation standard curve method, using U6 as an endogenous control¹⁰. Student's t test and Pearson's correlation were performed in GraphPad Prism 9.

## Mouse samples

Habenula samples were collected from unpredictable chronic mild stress (UCMS) and control male BALB/c mice from Hervé et al.,²³. Mice underwent the UCMS exposure (UCMS; n = 8) or kept in standard housing conditions (control; n = 7) for eight weeks. Control mice were housed in groups of four in standard cages (21 × 38 cm) while UCMS mice were isolated in individual home cages (8.5 × 22 cm) with no physical contact with other mice. Stressors consisted of housing on damp bedding (~ 200 mL H₂O for 100 g bedding), placement in an empty cage (home cage without bedding), placement in an empty cage with H₂O (home cage without bedding filled with 1 cm H₂O at 21 °C), switching cages, cage tilting (45°), predator sounds, introduction of predator feces/fur, inversion of light/dark cycle, and confinement in small tubes (4 cm diameter × 5 cm length). Stressors were varied and applied in different sequences per week to avoid habituation. Weight and coat state were measured weekly,

<table><caption>**Table 1.** Human Lateral Habenula Cohort Demographics. (mean ± S.E.M.). ¹RIN: RNA integrity number, ²PMI: post-mortem interval (in hours).</caption><thead><tr><td></td><td>Control</td><td>MDD</td><td>p-value</td></tr></thead><tbody><tr><td>Sex</td><td>5F/10 M</td><td>5F/10 M</td><td>N/A</td></tr><tr><td>Age</td><td>55.33 ± 5.12</td><td>54.27 ± 3.80</td><td>0.87</td></tr><tr><td>pH</td><td>6.35 ± 0.08</td><td>6.49 ± 0.07</td><td>0.18</td></tr><tr><td>RIN¹</td><td>5.41 ± 0.44</td><td>5.98 ± 0.25</td><td>0.27</td></tr><tr><td>PMI²</td><td>62.45 ± 7.95</td><td>55.67 ± 6.51</td><td>0.51</td></tr></tbody></table>
as markers of UCMS-induced state of physiological stress, except for the last week before sacrifice. At the end of the eighth week, a complementary test of nest building, a common measure of well-being in mice, was performed just before sacrifice<sup>24,25</sup>. The nest building test is administered at the beginning of the activity period of the mice. The test mice are isolated in their home cages and pieces of pressed cotton are introduced in the cage. The quality of the nest is scored 5 h, and 24 h later according to the scale proposed by Deacon 2016<sup>26</sup>. Mice were euthanized by CO₂ inhalation. For more detailed description please see Hervé et al.,<sup>23</sup>. All experiments on mice were carried out according to policies on the care and use of laboratory animals of European Community legislation 2010/63/ EU and ARRIVE guidelines. The local Ethics Committee (Comité d'Ethique en Expérimentation Animale de Val de Loire (CEEAVdL) number 19) approved the protocols used in this study (protocol number 2011-06-10).

## Snora69 expression

Total RNA was extracted from blood using Mouse RiboPure-Blood RNA isolation kit (Life Technologies) and from habenula using mirVana miRNA isolation kit (Life Technologies) according to manufacturer's instructions. All RNA samples were subjected to DNase treatment. Concentration was assessed using nanodrop ND-1000 and RNA quality was assessed using Agilent 2100 Bioanalyzer. RNA was reverse-transcribed using M-MLV Reverse Transcriptase (200 U/µL) (ThermoFisher) with random hexamers. Expression of Snora69 was assessed via RT-qPCR, as detailed above, using custom probes (Supplementary Table 1). Expression levels were calculated using the absolute quantitation standard curve method, using U6 as an endogenous control. Student's t test and Pearson's correlation were performed in GraphPad Prism 9.

## Human neuronal culture

To test if antidepressant drugs have any effects on the expression of SNORA69 we utilized human NPC cultures were as described in Fiori et al.<sup>27</sup>. Human neuronal cultures were screened for cytotoxic effects using the MTT assay, and antidepressants were applied at nontoxic concentrations as described by Lopez et al.,<sup>28</sup>. Following 2 weeks of neuronal differentiation, culture media was supplemented with duloxetine (10 µM; Sigma-Aldrich), escitalopram (100 µM; Sigma-Aldrich), haloperidol (10 µM; Sigma-Aldrich), lithium (1 mM; Sigma-Aldrich), Aspirin (1 mM; Sigma-Aldrich), or left untreated (control). Cells for each drug treatment were incubated for 48 h before being harvested for RNA extractions. Each drug treatment was performed in triplicate. RNA was extracted using the Zymo DirectZol RNA Extraction kit and cDNA construction and RT-qPCR were as described above. One-way ANOVA was performed in IBM SPSS Statistics version 27.

## N-cyclohexyl-N'-β-(4-methylmorpholinium) ethylcarbodiimide (CMC) treatment

The CMC treatment protocol was adapted from Antonicka et al. and Nagasawa et al., with slight modifications, to measure pseudouridylation levels in association with SNORA69 expression<sup>29,30</sup>. CMC (Sigma-Aldrich) treatment was conducted on the above-mentioned RNA extracted from human post-mortem LHB samples. 1.5 ug of total RNA was used per sample and split into two separate aliquots (0.9 ug for CMC treatment and 0.6ug for negative control; each fraction was in a final volume of 20 uL). 2.9 uL of 40 mM EDTA (pH 8.0; Sigma-Aldrich) was pipette mixed, briefly centrifuged, placed into an 80 °C water bath for 3 min, and immediately placed on ice. Freshly made 0.4 M CMC was prepared in BEU buffer (50 mM bicine (pH 8.5; Sigma-Aldrich), 4 mM EDTA, and 7 M urea (Sigma-Aldrich)). 100 uL of CMC in BEU buffer or BEU buffer without CMC was added to each corresponding aliquot per sample (total of 20 RNA samples) pipette mixed and briefly centrifuged. Samples were incubated at 40 °C for 60 min at 1000 rpm with shaking. RNA was precipitated using 1000 uL 100% ethanol, 50 uL of 3 M sodium acetate (Sigma-Aldrich), and 1 uL glycogen (Invitrogen). Each sample was pipette mixed and briefly centrifuged before placing in -80 °C for 30 min. Afterwards, samples were centrifuged at 16,000 g for 30 min at 4 °C. The supernatant was carefully discarded, 500 uL of cold 70% ethanol was added, and samples were centrifuged at 16,000 g for 10 min at 4 °C (this step was performed twice). After discarding the supernatant from the last ethanol wash, the samples were left to air dry. 30 uL of 50 mM sodium carbonate (pH 10.4; BioShop) and 2 mM EDTA was added to each sample, pipette mixed, briefly centrifuged, and incubated at 50 °C for 2 h at 1000 rpm with shaking. RNA was precipitated using 90 uL 100% ethanol, 3.5 uL of 3 M sodium acetate, and 1 uL glycogen; briefly centrifuged and placed in -80 °C for 30 min and centrifuged at 16,000 g for 30 min at 4 °C. The supernatant was removed and 500uL of 70% ethanol was added and centrifuged at 16,000 g for 10 min at 4 °C (this step was performed twice). After discarding the supernatant from the last ethanol wash, the samples were left to air dry and then resuspended in 15 uL H2O.

## Pseudouridine (Ψ) quantification

RNA was reverse-transcribed as described above. We specifically quantified pseudouridylation at U<sup>36</sup> on 18S rRNA and U<sup>69</sup> on 5.8S rRNA using custom designed primers up- and down-stream of each pseudouridylation site (supplementary Fig. 1 & Supplementary Table 1). CMC treatment causes a stall in reverse transcription at pseudouridylated (Ψ) sites. This will cause a truncated fragment in CMC treated RNA compared to non-CMC treated RNA (Supplementary Fig. 1). The abundance of Ψ is quantified by the ratio (R) of fold-change of non-truncated fragments (long fragment) with or without CMC treatment to the fold-change of truncated fragments (short fragment) with or without CMC treatment using the following equation:

$$R = \frac{2^{\frac{\Delta \text{Ct long fragment (nonCMC mean} - \text{CMC mean)}}{\text{long fragment}}}}{2^{\frac{\Delta \text{Ct short fragment (nonCMC mean} - \text{CMC mean)}}{\text{short fragment}}}}$$
To better visualize the ratio, we plotted the inverse ratio (1/R) labeled as “relative Ψ”. Mann-Whitney U test and Pearson's correlation was performed in GraphPad Prism 9.

## Results

### SNORA69 is up-regulated in the lateral habenula of depressed humans and mice

We performed small RNA sequencing in 30 samples (15 individuals with MDD, and 15 psychiatrically healthy controls) from the lateral habenula (LHb). Groups were matched for age, sex, post-mortem interval (PMI), pH, and RIN (Table 1). We identified 184 unique snoRNAs that passed filtering (as described in the methods section), with 11 snoRNAs displaying significant (FDR < 0.05) differential expression between MDD and controls (Supplementary Table 2). To narrow down snoRNA candidates for technical validation we applied a more stringent FDR < 0.01 which left use with SNORA54, SNORA80E, SNORA26, and SNORA69 (Fig. 1A & Supplementary Table 2). Of the four candidate snoRNAs, SNORA69 and SNORA54 passed RT-qPCR validated (Fig. 1B-E).

Previous studies suggest that the LHb acts as a relay center for forebrain limbic regions and midbrain monoaminergic centers; playing a central role in reward or aversion related learning and memory associated to stressful and aversive stimuli³¹⁻³³. Thus, we investigated the expression of our candidate snoRNAs in an unpredictable chronic mild stress (UCMS) mice model, a common model of depression²³,³⁴,³⁵. Given that SNORD54 (Accession No. XR_004934678.1) does not have a well-defined ortholog in mice, we focused our attention exclusively on SNORA69 (Accession No. NR_002900.1) moving forward³⁶. UCMS mice showed a significant increase in weight gain compared to unstressed mice by the end of the UCMS protocol. Additionally, depressive-like physical alterations, evaluated through coat state score, distinguished unstressed from UCMS mice, while assessment via the nest-building test revealed reduced motivation and apathetic behavior in UCMS mice²³. Similar to what we observed in humans, Snora69 was significantly elevated in UCMS mice, in both the LHb and peripheral blood samples (Fig. 1F). Snora69 also displayed positive correlation between the LHb and peripheral blood suggesting Snora69 may be a potential biomarker candidate for MDD, although this did not reach statistical significance (Fig. 1F).

### SNORA69 expression is not altered by antidepressant drugs

We further investigated if antidepressant drugs have effects on SNORA69 expression by using human neuronal cultures treated with four different psychotherapeutic drugs (duloxetine, escitalopram, haloperidol, and lithium), a non-psychotropic drug control (aspirin) and no drug control. Interestingly, the expression of SNORA69 was not significantly altered by any of the abovementioned drugs (Supplementary Fig. 2).

### SNORA69 is associated with an increase in pseudouridylation on 5.8S rRNA and 18S rRNA

Pseudouridine is an abundant RNA modification and affects the structure of the modified RNA and their interactions with protein and/or RNA binding partners, making this modification especially important in rRNA functioning³⁷. SNORA69 is known to guide pseudouridylation onto 5.8S rRNA U⁶⁹ and 18S rRNA U³⁶ (Fig. 1G,H)³⁸⁻⁴¹. We investigated the amount of pseudouridines at these two locations on each corresponding rRNA in our human LHb samples. We observed a significant increase of pseudouridines on both 5.8S rRNA (U⁶⁹) and 18S rRNA (U³⁶) (Fig. 1G,H). However, correlation between SNORA69 expression and pseudouridylation was only observed for 18S rRNA (Supplementary Fig. 3A). There were no significant differences in the expression of either 5.8S rRNA nor 18S rRNA and no significant correlation between the expression of SNORA69 and expression of each corresponding rRNA (Supplementary Fig. 3B-C). This indicates that SNORA69 is not influencing the quantities of 5.8S rRNA and 18S rRNA, but more likely influencing the functioning of these rRNA through increasing the amount of pseudouridylation.

## Discussion

Despite the rapidly growing body of work highlighting the role of non-coding RNAs in MDD, as well as other psychiatric conditions, the involvement of snoRNAs in psychiatric disorders remains under investigated¹¹⁻¹³,²⁸,⁴². Our study suggests that snoRNA expression, in the LHb, are associated with MDD. From our discovery analysis we identified 184 unique snoRNAs with 22% falling into the H/ACA box snoRNA family and 74% falling into the C/D box snoRNA family. Of the snoRNAs identified, 11 snoRNAs were differentially expressed in individuals with MDD and died by suicide. It's interesting to note that although that majority of snoRNAs identified were C/D box snoRNAs those that were significantly differentially expressed in MDD were primarily H/ACA box snoRNAs. Of the top 11 snoRNAs seven were H/ACA box snoRNAs whereas only three were C/D box snoRNAs. Moreover, the top four snoRNAs that passed a more stringent FDR < 0.01 cutoff were all H/ACA box snoRNAs. The exact meaning and implication of these observations is unclear, but perhaps there might be some association between common mechanisms that elicit pseudouridylation onto target rRNAs with those involved in MDD. Upon technical validation and investigation in a UCMS mouse model, our study points to SNORA69 to be strongly associated with MDD. Furthermore, in our cell culture experiments, SNORA69 expression was not altered by psychotropic drugs including antidepressants.

Previous studies suggest that the dorsal raphe and LHb have reciprocal projections and that decreased serotonin transmission can disinhibit the LHb possibly mediating depressive symptoms⁴³⁻⁴⁵. Reciprocally the LHb has also been shown to modulate dorsal raphe nuclei activity⁴⁴,⁴⁵. This of particular interest since our observations for SNORA69 seem to be specific to MDD without any expressional change when exposed to antidepressant treatment, which primarily targets the serotonergic system.

While no previous studies have directly investigated snoRNAs in the context of MDD several studies have demonstrated the importance of snoRNAs in the brain⁴⁶,⁴⁷. Specifically, a cluster of brain enriched snoRNAs including SNORD115 (H/MBII-52) and SNORD116 (H/MBII-85) have received attention for their potential
<figure><img src="image_3.png" /><figcaption>**Figure 1.** SNORA69 is up-regulated in the LHB of individuals with MDD and regulates pseudouridine abundance. (A) Volcano plot of 184 snoRNA detected from small RNA sequencing in post-mortem human LHB tissue. We applied a cut-off of FDR < 0.01 resulting in four candidate snoRNAs: SNORA54, SNORA80E, SNORA26, and SNORA69. We also highlight SNORD115 and SNORD116 although they did not reach significance. (B-E) RT-qPCR validation of candidate snoRNAs. Only SNORA54 (B) and SNORA69 (C) were significantly validated. (F) SNORA69 expression was significantly up-regulated in the LHB (left) and peripheral blood (center) of mice that underwent unpredicted chronic mild stress (UCMS). SNORA69 expression was positively correlated between the LHB and peripheral blood (right). White circles are controls and pink circles are UCMS. (G, H) SNORA69 guides pseudouridylation onto 5.8S rRNA (G) and 18S rRNA (H). We observed an increased detection of pseudouridines on 5.8S rRNA (G) and 18S rRNA (H) in individuals with MDD compared to controls. All bar plots represent the mean with individual data points as dots. Ψ denotes pseudouridylation event and the superscript number indicates the nucleotide position on the SNORA69 sequence. Error bars represent S.E.M. Student's two-tailed t test (B-F) or Mann-Whitney U test (G-H) were used to assess group differences. Pearson correlation coefficient (r) was used for correlation analysis. \*p<0.05, \*\*p<0.01.</figcaption></figure>
involvement in Prader-Willi syndrome (PWS). PWS is a complex genetic disorder caused by the loss of function of paternally imprinted genes in the 15q11-q13 region, of which SNORD115 and SNORD116 clusters are located⁴⁸,⁴⁹. Several studies have created mice models where Snord116 is knocked out (KO) and mice exhibit phenotypic behaviors reminiscent of PWS⁵⁰⁻⁵⁴. Although it remains unclear how SNORD116 contributes to PWS, studies have shown that SNORD116 can influence gene expression and DNA methylation⁵⁵,⁵⁶. SNORD115 has been shown to regulate the editing or alternative splicing of serotonin receptor 2C⁵⁷⁻⁵⁹. SNORD115 KO studies in mice do not result in robust changes in behavior; however, SNORD115-KO mice present changes in monoaminergic circuits specifically increasing the firing rate of dopaminergic neurons in the ventral tegmental area and serotonergic neurons in the dorsal raphe nucleus⁶⁰,⁶¹. We note that, although we did identify SNORD116 and SNORD115 in our human post-mortem LHB samples, neither were significantly differentially expressed between MDD and controls (Fig. 1A).

Using small RNA sequencing we were able to associate SNORA69 with MDD. However, we recognize that our study is not without limitations. The sample size of our human and mouse LHB cohort is limited, but robust enough to identify and validate SNORA69. Additionally, our study focused on snoRNA expression in the LHB, precluding investigation in other brain areas implicated in MDD. While we did not observe any effect of SNORA69 expression with antidepressant treatment, our investigation was confined to human neuronal culture whereas different cell types may have different outcomes. Although we demonstrated that an up-regulation of SNORA69 is associated with an increase of pseudouridylation onto 5.8S and 18S rRNA, the direct relevance of this observation to the psychopathology of MDD remains uncertain. We speculate that alterations to rRNA pseudouridylation will influence binding to tRNA and mRNA, ultimately having implications on translation and proteomic profiles in MDD; however further studies are needed to confirm this⁶²,⁶³. Moreover, while the UCMS mouse model is recognized as a well-established mouse model of depression, we acknowledge the limitation of not conducting a comprehensive battery of behavioral tests on these animals. Lastly, due to the heterogenous nature of MDD and of the subjects included in the study, it remains unclear if SNORA69 has a general association with MDD or with a confounding factor that may be commonly associated with MDD. Further work is needed to better characterize the role of SNORA69 in MDD. Furthermore, our mouse data suggest that the expression of Snora69 is positively correlated between LHB and peripheral blood making it an interesting candidate for further investigations as a predictive biomarker for MDD or mediator marker contributing to MDD pathology.

## Data availability

Small RNA sequencing data used in this article can be obtained from NCBI GEO under accession number GSE244651.

Received: 19 September 2023; Accepted: 27 March 2024

Published online: 09 April 2024

## References

1. Collaborators C-MD. Global prevalence and burden of depressive and anxiety disorders in 204 countries and territories in 2020 due to the COVID-19 pandemic. Lancet. 398(10312): 1700-12 (2021).

2. Collaborators GBDMD. Global, regional, and national burden of 12 mental disorders in 204 countries and territories, 1990-2019: A systematic analysis for the Global Burden of Disease Study 2019. Lancet Psychiatry. 9(2): 137-50 (2022).

3. Diseases, G. B. D. & Injuries, C. Global burden of 369 diseases and injuries in 204 countries and territories, 1990-2019: A systematic analysis for the Global burden of disease study 2019. Lancet 396(10258), 1204-1222 (2020).

4. Browne, C. A., Hammack, R. & Lucki, I. Dysregulation of the lateral habenula in major depressive disorder. Front. Synaptic Neurosci. 10, 46 (2018).

5. Gold, P. W. & Kadriu, B. A major role for the lateral habenula in depressive illness: Physiologic and molecular mechanisms. Front. Psychiatry 10, 320 (2019).

6. Yang, Y., Wang, H., Hu, J. & Hu, H. Lateral habenula in the pathophysiology of depression. Curr. Opin. Neurobiol. 48, 90-96 (2018).

7. Aan Het Rot, M., Mathew, S. J. & Charney, D. S. Neurobiological mechanisms in major depressive disorder. CMAJ. 180(3), 305-313 (2009).

8. Turecki, G. The molecular bases of the suicidal brain. Nat. Rev. Neurosci. 15(12), 802-816 (2014).

9. Sun, H., Kennedy, P. J. & Nestler, E. J. Epigenetics of the depressed brain: Role of histone acetylation and methylation. Neuropsychopharmacology 38(1), 124-137 (2013).

10. Fiori, L. M. et al. miR-323a regulates ERBB4 and is involved in depression. Mol. Psychiatry 26(8), 4191-4204 (2021).

11. Issler, O. et al. Sex-specific role for the long non-coding RNA LINC00473 in depression. Neuron 106(6), 912-926 e5 (2020).

12. Zhang, Y. et al. CircDYM ameliorates depressive-like behavior by targeting miR-9 to regulate microglial activation via HSP90 ubiquitination. Mol. Psychiatry. 25(6), 1175-1190 (2020).

13. Lin, R., Kos, A., Lopez, J. P., Dine, J., Fiori, L. M., Yang, J., et al. SNORD90 induces glutamatergic signaling following treatment with monoaminergic antidepressants. Elife. 12 (2023).

14. Calvet, J. P. & Pederson, T. Base-pairing interactions between small nuclear RNAs and nuclear RNA precursors as revealed by psoralen cross-linking in vivo. Cell 26(3 Pt 1), 363-370 (1981).

15. Edstroem, J. E. & Gall, J. G. The base composition of ribonucleic acid in lampbrush chromosomes, nucleoli, nuclear sap, and cytoplasm of triturus oocytes. J. Cell Biol. 19, 279-284 (1963).

16. Stern, H., Johnston, F. B. & Setterfield, G. Some chemical properties of isolated pea nucleoli. J. Biophys. Biochem. Cytol. 6(1), 57-60 (1959).

17. Wise, J. A. & Weiner, A. M. Dictyostelium small nuclear RNA D2 is homologous to rat nucleolar RNA U3 and is encoded by a dispersed multigene family. Cell 22(1 Pt 1), 109-118 (1980).

18. Bachellerie, J. P., Cavaille, J. & Huttenhofer, A. The expanding snoRNA world. Biochimie 84(8), 775-790 (2002).

19. Kufel, J. & Grzechnik, P. Small nucleolar RNAs tell a different tale. Trends Genet. 35(2), 104-117 (2019).

20. Huang, Z. H., Du, Y. P., Wen, J. T., Lu, B. F. & Zhao, Y. snoRNAs: Functions and mechanisms in biological processes, and roles in tumor pathophysiology. Cell Death Discov. 8(1), 259 (2022).

21. Labuda, M. et al. Linkage disequilibrium analysis in young populations: Pseudo-vitamin D-deficiency rickets and the founder effect in French Canadians. Am. J. Hum. Genet. 59(3), 633-643 (1996).
22. Dumais, A. et al. Risk factors for suicide completion in major depression: A case-control study of impulsive and aggressive behaviors in men. *Am. J. Psychiatry* **162**(11), 2116–2124 (2005).

23. Herve, M. et al. Translational identification of transcriptional signatures of major depression and antidepressant response. *Front. Mol. Neurosci.* **10**, 248 (2017).

24. Jirkof, P. Burrowing and nest building behavior as indicators of well-being in mice. *J. Neurosci. Methods.* **234**, 139–146 (2014).

25. Kraeuter, A. K., Guest, P. C. & Sarnyai, Z. The nest building test in mice for assessment of general well-being. *Methods Mol. Biol.* **1916**, 87–91 (2019).

26. Deacon, R. M. Assessing nest building in mice. *Nat. Protoc.* **1**(3), 1117–1119 (2006).

27. Fiori, L. M. et al. Treatment-emergent and trajectory-based peripheral gene expression markers of antidepressant response. *Transl. Psychiatry.* **11**(1), 439 (2021).

28. Lopez, J. P. et al. miR-1202 is a primate-specific and brain-enriched microRNA involved in major depression and antidepressant treatment. *Nat. Med.* **20**(7), 764–768 (2014).

29. Antonicka, H. et al. A pseudouridine synthase module is essential for mitochondrial protein synthesis and cell viability. *EMBO Rep.* **18**(1), 28–38 (2017).

30. Nagasawa, C. K., Kibiryeva, N., Marshall, J., O'Brien, J. E. Jr. & Bittel, D. C. scaRNA1 levels alter pseudouridylation in spliceosomal RNA U2 affecting alternative mRNA splicing and embryonic development. *Pediatr. Cardiol.* **41**(2), 341–349 (2020).

31. Baker, P. M. et al. Lateral habenula beyond avoidance: Roles in stress, memory, and decision-making with implications for psychiatric disorders. *Front. Syst. Neurosci.* **16**, 826475 (2022).

32. Cohen, J. Y., Haesler, S., Vong, L., Lowell, B. B. & Uchida, N. Neuron-type-specific signals for reward and punishment in the ventral tegmental area. *Nature* **482**(7383), 85–88 (2012).

33. Proulx, C. D., Hikosaka, O. & Malinow, R. Reward processing by the lateral habenula in normal and depressive behaviors. *Nat. Neurosci.* **17**(9), 1146–1152 (2014).

34. Nollet, M. Models of depression: Unpredictable chronic mild stress in mice. *Curr. Protoc.* **1**(8), e208 (2021).

35. Burstein, O., Doron, R. The unpredictable chronic mild stress protocol for inducing anhedonia in mice. *J. Vis. Exp.* 2018(140).

36. Sayers, E. W. et al. Database resources of the national center for biotechnology information. *Nucleic Acids Res.* **50**(D1), D20–D26 (2022).

37. Borchardt, E. K., Martinez, N. M. & Gilbert, W. V. Regulation and function of RNA Pseudouridylation in human cells. *Annu. Rev. Genet.* **54**, 309–336 (2020).

38. Ganot, P., Bortolin, M. L. & Kiss, T. Site-specific pseudouridine formation in preribosomal RNA is guided by small nucleolar RNAs. *Cell* **89**(5), 799–809 (1997).

39. Ganot, P., Caizergues-Ferrer, M. & Kiss, T. The family of box ACA small nucleolar RNAs is defined by an evolutionarily conserved secondary structure and ubiquitous sequence elements essential for RNA accumulation. *Genes. Dev.* **11**(7), 941–956 (1997).

40. Lestrade, L. & Weber, M. J. snoRNA-LBME-db, a comprehensive database of human H/ACA and C/D box snoRNAs. *Nucleic Acids Res.* **34**, D158–D162 (2006).

41. Maden, E. H. & Wakeman, J. A. Pseudouridine distribution in mammalian 18 S ribosomal RNA. A major cluster in the central region of the molecule. *Biochem. J.* **249**(2), 459–464 (1988).

42. Issler, O. et al. The long noncoding RNA FEDORA is a cell type- and sex-specific regulator of depression. *Sci. Adv.* **8**(48), eabn9494 (2022).

43. Zhang, H. et al. Dorsal raphe projection inhibits the excitatory inputs on lateral habenula and alleviates depressive behaviors in rats. *Brain Struct. Funct.* **223**(5), 2243–2258 (2018).

44. Zhao, H., Zhang, B. L., Yang, S. J. & Rusak, B. The role of lateral habenula-dorsal raphe nucleus circuits in higher brain functions and psychiatric illness. *Behav. Brain Res.* **277**, 89–98 (2015).

45. Takahashi, A. et al. Lateral habenula glutamatergic neurons projecting to the dorsal raphe nucleus promote aggressive arousal in mice. *Nat. Commun.* **13**(1), 4039 (2022).

46. Rogelj, B. Brain-specific small nucleolar RNAs. *J. Mol. Neurosci.* **28**(2), 103–109 (2006).

47. Rogelj, B., Hartmann, C. E., Yeo, C. H., Hunt, S. P. & Giese, K. P. Contextual fear conditioning regulates the expression of brain-specific small nucleolar RNAs in hippocampus. *Eur. J. Neurosci.* **18**(11), 3089–3096 (2003).

48. Gallagher, R. C., Pils, B., Albalwi, M. & Francke, U. Evidence for the role of PWCR1/HBII-85 C/D box small nucleolar RNAs in Prader-Willi syndrome. *Am. J. Hum. Genet.* **71**(3), 669–678 (2002).

49. Runte, M., Varon, R., Horn, D., Horsthemke, B. & Buiting, K. Exclusion of the C/D box snoRNA gene cluster HBII-52 from a major role in Prader-Willi syndrome. *Hum. Genet.* **116**(3), 228–230 (2005).

50. Adhikari, A. et al. Cognitive deficits in the Snord116 deletion mouse model for Prader-Willi syndrome. *Neurobiol. Learn Mem.* **165**, 106874 (2019).

51. Polex-Wolf, J. et al. Hypothalamic loss of Snord116 recapitulates the hyperphagia of Prader-Willi syndrome. *J. Clin. Invest.* **128**(3), 960–969 (2018).

52. Qi, Y. et al. Snord116 is critical in the regulation of food intake and body weight. *Sci. Rep.* **6**, 18614 (2016).

53. Zieba, J. et al. Behavioural characteristics of the Prader-Willi syndrome related biallelic Snord116 mouse model. *Neuropeptides* **53**, 71–77 (2015).

54. Salles, J. et al. What can we learn from PWS and SNORD116 genes about the pathophysiology of addictive disorders?. *Mol Psychiatry* **26**(1), 51–59 (2021).

55. Coulson, R. L. et al. Snord116-dependent diurnal rhythm of DNA methylation in mouse cortex. *Nat. Commun.* **9**(1), 1616 (2018).

56. Falaleeva, M., Surface, J., Shen, M., de la Grange, P. & Stamm, S. SNORD116 and SNORD115 change expression of multiple genes and modify each other's activity. *Gene* **572**(2), 266–273 (2015).

57. Kishore, S. et al. The snoRNA MBII-52 (SNORD 115) is processed into smaller RNAs and regulates alternative splicing. *Hum. Mol. Genet.* **19**(7), 1153–1164 (2010).

58. Kishore, S. & Stamm, S. The snoRNA HBII-52 regulates alternative splicing of the serotonin receptor 2C. *Science* **311**(5758), 230–232 (2006).

59. Raabe, C. A. et al. Ectopic expression of Snord115 in choroid plexus interferes with editing but not splicing of 5-Ht2c receptor pre-mRNA in mice. *Sci. Rep.* **9**(1), 4300 (2019).

60. Hebras, J., Marty, V., Personnaz, J., Mercier, P., Krogh, N., Nielsen, H. et al. Reassessment of the involvement of Snord115 in the serotonin 2c receptor pathway in a genetically relevant mouse model. *Elife.* **9** (2020).

61. Marty, V., Butler, J., Basile, C., Chargui, O., Chagraoui, A., Guiard, B. P. et al. Deleting Snord115 genes in mice remodels monoaminergic systems activity in the brain toward cortico-subcortical imbalances. *Hum. Mol. Genet.* (2022).

62. Blaze, J. et al. Neuronal Nsun2 deficiency produces tRNA epitranscriptomic alterations and proteomic shifts impacting synaptic signaling and behavior. *Nat Commun.* **12**(1), 4913 (2021).

63. Penzo, M. & Montanaro, L. Turning uridines around: Role of rRNA pseudouridylation in ribosome biogenesis and ribosomal function. *Biomolecules* **8**(2), 38 (2018).
Acknowledgements

GT holds a Canada Research Chair (Tier 1) and is supported by Grants from the Canadian Institute of Health Research (CIHR) (FDN148374, EGM141899, ENP161427), and by the Fonds de recherche du Québec -Santé (FRQS) through the Quebec Network on Suicide, Mood Disorders, and Related Disorders.

Author contributions

R.L. and G.T. conceptualized the project. G.T. and N.M. contributed to the collection of human post-mortem brain samples. R.L. and L.M.F. performed RNA extraction, and small RNA sequencing from human post-mortem brain. R.D. performed bioinformatic analysis. E.C.I. and C.B. provided RNA samples from UCMS mouse model. R.L. performed cell culture experiments. R.L. and H.M. performed cDNA construction and RT-qPCR assays. R.L. and G.T. wrote the manuscript.

Competing interests

The authors declare no competing interests.

Additional information

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-024-58278-2.

Correspondence and requests for materials should be addressed to G.T.

Reprints and permissions information is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_4.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2024