# Decreased telomere length in a subgroup of young individuals with bipolar disorders: replication in the FACE-BD cohort and association with the shelterin component POT1

Luana Spano, Cynthia Marie-Claire, Ophélia Godin, Apolline Lebras, Cindie Courtin, Jean-Louis Laplanche, Marion Leboyer, Bruno Aouizerate, Antoine Lefrere, Raoul Belzeaux, et al.

► To cite this version:

Luana Spano, Cynthia Marie-Claire, Ophélia Godin, Apolline Lebras, Cindie Courtin, et al. Decreased telomere length in a subgroup of young individuals with bipolar disorders: replication in the FACE-BD cohort and association with the shelterin component POT1. Translational Psychiatry, 2024, 14 (1), pp.131. 10.1038/s41398-024-02824-z. hal-04492634

HAL Id: hal-04492634

https://hal.science/hal-04492634

Submitted on 21 May 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>

Distributed under a Creative Commons Attribution 4.0 International License
# Decreased telomere length in a subgroup of young individuals with bipolar disorders: replication in the FACE-BD cohort and association with the shelterin component POT1

Luana Spano, Cynthia Marie-Claire, Ophélia Godin, Apolline Lebras, Cindie Courtin, Jean-Louis Laplanche, Marion Leboyer, Bruno Aouizerate, Antoine Lefrere, Raoul Belzeaux, et al.

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

## ► To cite this version:

Luana Spano, Cynthia Marie-Claire, Ophélia Godin, Apolline Lebras, Cindie Courtin, et al. Decreased telomere length in a subgroup of young individuals with bipolar disorders: replication in the FACE-BD cohort and association with the shelterin component POT1. Translational Psychiatry, 2024, 14 (1), pp.131. 10.1038/s41398-024-02824-z. hal-04492634

HAL Id: hal-04492634

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

Submitted on 21 May 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_2.png" /></figure>
**ARTICLE** **OPEN**

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

Check for updates

# Decreased telomere length in a subgroup of young individuals with bipolar disorders: replication in the FACE-BD cohort and association with the shelterin component POT1

Luana Spano<fnref n="1" />, Cynthia Marie-Claire ID<fnref n="1" />✉, Ophélia Godin ID²,³, Apolline Lebras<fnref n="1" />, Cindie Courtin<fnref n="1" />, Jean-Louis Laplanche<fnref n="1" />,⁴, Marion Leboyer ID²,³,⁵, Bruno Auizerate²,⁶, Antoine Lefrere²,⁷, Raoul Belzeaux²,⁸,⁹, Philippe Courtet²,¹⁰, Emilie Olié²,¹⁰, Caroline Dubertret²,¹¹,¹², Raymund Schwan²,¹³, Valérie Aubin²,¹⁴, Paul Roux²,¹⁵,¹⁶, Mircea Polosan²,¹⁷, Ludovic Samalin ID²,¹⁸, Emmanuel Haffen²,¹⁹, Fondamental Advanced Centers Of Expertise In Bipolar Disorders (Face-Bd) Collaborators*, Frank Bellivier<fnref n="1" />,²,²⁰ and Bruno Etain ID<fnref n="1" />,²,²⁰

© The Author(s) 2024

Bipolar disorder (BD) has been associated with premature cellular aging with shortened telomere length (TL) as compared to the general population. We recently identified a subgroup of young individuals with prematurely shortened TL. The aims of the present study were to replicate this observation in a larger sample and analyze the expression levels of genes associated with age or TL in a subsample of these individuals. TL was measured on peripheral blood DNA using quantitative polymerase chain reaction in a sample of 542 individuals with BD and clustering analyses were performed. Gene expression level of 29 genes, associated with aging or with telomere maintenance, was analyzed in RNA samples from a subsample of 129 individuals. Clustering analyses identified a group of young individuals (mean age 29.64 years), with shorter TL. None of the tested clinical variables were significantly associated with this subgroup. Gene expression level analyses showed significant downregulation of *MYC*, *POT1*, and *CD27* in the prematurely aged young individuals compared to the young individuals with longer TL. After adjustment only *POT1* remained significantly differentially expressed between the two groups of young individuals. This study confirms the existence of a subgroup of young individuals with BD with shortened TL. The observed decrease of *POT1* expression level suggests a newly described cellular mechanism in individuals with BD, that may contribute to telomere shortening.

*Translational Psychiatry* (2024)14:131; https://doi.org/10.1038/s41398-024-02824-z

## INTRODUCTION

A 10–15 years decrease in life expectancy has been repeatedly observed in individuals with bipolar disorder (BD) as compared to the general population [1, 2]. Individuals with BD experience higher rates of comorbidities associated with aging, such as type II diabetes, metabolic syndrome, cancer, immune dysregulation, cardiovascular and cerebrovascular disorders [3–5]. These observations suggest that BD may be associated with a premature cellular aging [6], which can be estimated by various molecular markers [7].

Telomere length (TL) is the most well-described biomarker of cellular aging. Telomeres are repeated sequences of non-coding nucleotides (TTAGGG). Since they are located at the end of chromosomes, they protect chromosomes from their degradation and fusion [8]. However, at each cell division, the replication of DNA is incomplete, leading to telomere shortening with aging [9]. Several mechanisms are at stake to preserve TL, such as the “shelterin complex” which involves the interaction between six proteins and the telomerase, an enzyme that adds TTAGGG

<footnote n="1">¹Université Paris Cité, INSERM UMR-S 1144, Optimisation Thérapeutique en Neuropsychopharmacologie OTeN, Paris, France. ²Fondation FondaMental, Créteil, France. ³Université Paris Est Créteil, INSERM U955, IMRB, Translational NeuroPsychiatry Laboratory, Créteil, France. ⁴Département de Biochimie et Biologie Moléculaire, DMU BioGeM, Hôpitaux Lariboisière-Fernand Widal, GHU APHP.Nord – Université de Paris, Paris, France. ⁵AP-HP, Hôpitaux Universitaires Henri Mondor, Département Médico-Universitaire de Psychiatrie et d'Addictologie (DMUIMPACT), Fédération Hospitalo-Universitaire de Médecine de Précision en Psychiatrie (FHU ADAPT), Créteil, France. ⁶Centre Hospitalier Charles Perrens, Laboratoire NutriNeuro (UMR INRA 1286), Université de Bordeaux, Bordeaux, France. ⁷Pôle de Psychiatrie, Assistance Publique Hôpitaux de Marseille, Marseille, France. ⁸INT-UMR7289, CNRS Aix-Marseille Université, Marseille, France. ⁹Université Montpellier, Montpellier, France. ¹⁰Department of Emergency Psychiatry and Acute Care, CHU Montpellier, IGF, Univ. Montpellier, CNRS, INSERM, Montpellier, France. ¹¹AP-HP, Groupe Hospitalo-Universitaire AP-HP Nord, DMU ESPRIT, Service de Psychiatrie et Addictologie, Hôpital Louis Mourier, Colombes, France. ¹²Université de Paris, Inserm UMR1266, Sorbonne Paris Cité, Faculté de Médecine, Paris, France. ¹³Université de Lorraine, Centre Psychothérapique de Nancy, Inserm U1254, Nancy, France. ¹⁴Pôle de Psychiatrie, Centre Hospitalier Princesse Grace, Monaco, Monaco. ¹⁵Centre Hospitalier de Versailles, Service Universitaire de Psychiatrie d'Adulte et d'Addictologie, Le Chesnay, France. ¹⁶Équipe DisAP-PsyDev, CESP, Université Versailles Saint- Quentin-en-Yvelines – Paris-Saclay, Inserm, Villejuif, France. ¹⁷Université Grenoble Alpes, Inserm, U1216, CHU Grenoble Alpes, Grenoble Institut Neurosciences, Grenoble, France. ¹⁸Centre Hospitalier et Universitaire, Département de Psychiatrie, Université Clermont Auvergne, CNRS, Clermont Auvergne INP, Institut Pascal (UMR 6602), Clermont-Ferrand, France. ¹⁹Service de Psychiatrie de l'Adulte, CIC-1431 INSERM, CHU de Besançon, Laboratoire de Neurosciences, UFC, UBFC, Besançon, France. ²⁰AP-HP, Groupe Hospitalo-Universitaire AP-HP Nord, DMU Neurosciences, Hôpital Fernand Widal, Département de Psychiatrie et de Médecine Addictologique, Paris, France. *A list of authors and their affiliations appears at the end of the paper.
✉email: cynthia.marie-claire@u-paris.fr</footnote>
**ARTICLE** **OPEN**

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

Check for updates

# Decreased telomere length in a subgroup of young individuals with bipolar disorders: replication in the FACE-BD cohort and association with the shelterin component POT1

Luana Spano<fnref n="1" />, Cynthia Marie-Claire ID <fnref n="1" />✉, Ophélia Godin ID ²,³, Apolline Lebras<fnref n="1" />, Cindie Courtin<fnref n="1" />, Jean-Louis Laplanche<fnref n="1" />,⁴, Marion Leboyer ID ²,³,⁵, Bruno Auizerate²,⁶, Antoine Lefrere²,⁷, Raoul Belzeaux²,⁸,⁹, Philippe Courtet²,¹⁰, Emilie Olié²,¹⁰, Caroline Dubertret²,¹¹,¹², Raymund Schwan²,¹³, Valérie Aubin²,¹⁴, Paul Roux²,¹⁵,¹⁶, Mircea Polosan²,¹⁷, Ludovic Samalin ID ²,¹⁸, Emmanuel Haffen²,¹⁹, Fondamental Advanced Centers Of Expertise In Bipolar Disorders (Face-Bd) Collaborators*, Frank Bellivier<fnref n="1" />,²,²⁰ and Bruno Etain ID <fnref n="1" />,²,²⁰

© The Author(s) 2024

Bipolar disorder (BD) has been associated with premature cellular aging with shortened telomere length (TL) as compared to the general population. We recently identified a subgroup of young individuals with prematurely shortened TL. The aims of the present study were to replicate this observation in a larger sample and analyze the expression levels of genes associated with age or TL in a subsample of these individuals. TL was measured on peripheral blood DNA using quantitative polymerase chain reaction in a sample of 542 individuals with BD and clustering analyses were performed. Gene expression level of 29 genes, associated with aging or with telomere maintenance, was analyzed in RNA samples from a subsample of 129 individuals. Clustering analyses identified a group of young individuals (mean age 29.64 years), with shorter TL. None of the tested clinical variables were significantly associated with this subgroup. Gene expression level analyses showed significant downregulation of *MYC*, *POT1*, and *CD27* in the prematurely aged young individuals compared to the young individuals with longer TL. After adjustment only *POT1* remained significantly differentially expressed between the two groups of young individuals. This study confirms the existence of a subgroup of young individuals with BD with shortened TL. The observed decrease of *POT1* expression level suggests a newly described cellular mechanism in individuals with BD, that may contribute to telomere shortening.

*Translational Psychiatry* (2024)14:131; https://doi.org/10.1038/s41398-024-02824-z

## INTRODUCTION

A 10–15 years decrease in life expectancy has been repeatedly observed in individuals with bipolar disorder (BD) as compared to the general population [1, 2]. Individuals with BD experience higher rates of comorbidities associated with aging, such as type II diabetes, metabolic syndrome, cancer, immune dysregulation, cardiovascular and cerebrovascular disorders [3–5]. These observations suggest that BD may be associated with a premature cellular aging [6], which can be estimated by various molecular markers [7].

Telomere length (TL) is the most well-described biomarker of cellular aging. Telomeres are repeated sequences of non-coding nucleotides (TTAGGG). Since they are located at the end of chromosomes, they protect chromosomes from their degradation and fusion [8]. However, at each cell division, the replication of DNA is incomplete, leading to telomere shortening with aging [9]. Several mechanisms are at stake to preserve TL, such as the “shelterin complex” which involves the interaction between six proteins and the telomerase, an enzyme that adds TTAGGG

<footnote n="1">¹Université Paris Cité, INSERM UMR-S 1144, Optimisation Thérapeutique en Neuropsychopharmacologie OTeN, Paris, France. ²Fondation FondaMental, Créteil, France. ³Université Paris Est Créteil, INSERM U955, IMRB, Translational NeuroPsychiatry Laboratory, Créteil, France. ⁴Département de Biochimie et Biologie Moléculaire, DMU BioGeM, Hôpitaux Lariboisière-Fernand Widal, GHU APHP.Nord – Université de Paris, Paris, France. ⁵AP-HP, Hôpitaux Universitaires Henri Mondor, Département Médico-Universitaire de Psychiatrie et d'Addictologie (DMUIMPACT), Fédération Hospitalo-Universitaire de Médecine de Précision en Psychiatrie (FHU ADAPT), Créteil, France. ⁶Centre Hospitalier Charles Perrens, Laboratoire NutriNeuro (UMR INRA 1286), Université de Bordeaux, Bordeaux, France. ⁷Pôle de Psychiatrie, Assistance Publique Hôpitaux de Marseille, Marseille, France. ⁸INT-UMR7289, CNRS Aix-Marseille Université, Marseille, France. ⁹Université Montpellier, Montpellier, France. ¹⁰Department of Emergency Psychiatry and Acute Care, CHU Montpellier, IGF, Univ. Montpellier, CNRS, INSERM, Montpellier, France. ¹¹AP-HP, Groupe Hospitalo-Universitaire AP-HP Nord, DMU ESPRIT, Service de Psychiatrie et Addictologie, Hôpital Louis Mourier, Colombes, France. ¹²Université de Paris, Inserm UMR1266, Sorbonne Paris Cité, Faculté de Médecine, Paris, France. ¹³Université de Lorraine, Centre Psychothérapique de Nancy, Inserm U1254, Nancy, France. ¹⁴Pôle de Psychiatrie, Centre Hospitalier Princesse Grace, Monaco, Monaco. ¹⁵Centre Hospitalier de Versailles, Service Universitaire de Psychiatrie d'Adulte et d'Addictologie, Le Chesnay, France. ¹⁶Équipe DisAP-PsyDev, CESP, Université Versailles Saint- Quentin-en-Yvelines – Paris-Saclay, Inserm, Villejuif, France. ¹⁷Université Grenoble Alpes, Inserm, U1216, CHU Grenoble Alpes, Grenoble Institut Neurosciences, Grenoble, France. ¹⁸Centre Hospitalier et Universitaire, Département de Psychiatrie, Université Clermont Auvergne, CNRS, Clermont Auvergne INP, Institut Pascal (UMR 6602), Clermont-Ferrand, France. ¹⁹Service de Psychiatrie de l'Adulte, CIC-1431 INSERM, CHU de Besançon, Laboratoire de Neurosciences, UFC, UBFC, Besançon, France. ²⁰AP-HP, Groupe Hospitalo-Universitaire AP-HP Nord, DMU Neurosciences, Hôpital Fernand Widal, Département de Psychiatrie et de Médecine Addictologique, Paris, France. *A list of authors and their affiliations appears at the end of the paper.
✉email: cynthia.marie-claire@u-paris.fr</footnote>
sequences to the telomeres [9]. TL, as a marker of aging, has been shown to be decreased in individuals with BD compared to the unaffected or healthy population in a meta-analysis [10]. Moreover, in BD, we have recently identified a subgroup of individuals in their thirties but with a TL that was as short as the one of those individuals in their sixties [11]. However, it remains unclear whether this decrease of TL observed in individuals with BD is associated with specific clinical characteristic, even if some potential determinants have been proposed [12]. For example, long-term lithium treatment would have a positive effect on telomere attrition [13] while childhood maltreatment [14], as well as a higher number of mood episodes [15], would have a negative impact on TL. Some alterations of circadian rhythms (i.e. circadian misalignment) might also contribute to shorten telomeres, as we suggested in a recent article [16].

The identification of genes whose expression levels change with age and, independently of chronological age, predict accelerated decline of health is an expanding area of research [17]. A transcriptomic age predictor has been proposed using a whole-blood gene expression meta-analysis in 14,983 individuals of European ancestry where 1,497 genes were identified as differently expressed with chronological age [18].

Among these age-related genes, those involved in telomere structure and maintenance have been also recently reviewed [19] and their expression levels studied in the general population. For example, a recent study of the data from a whole genome sequencing of 250 Dutch family trios analyzed TL inheritance patterns and associated genetic factors. *RRM1*, a gene playing a role in cellular proliferation, was suggested to be involved in TL regulation [20]. To our knowledge, no previous study has analyzed simultaneously the expression of age and TL-related genes in individuals with BD. Such a study in individuals with BD would help identifying potential molecular determinants of telomere shortening in this population.

In this study, we first analyzed TL in individuals with BD to replicate, in a larger cohort and with a larger set of clinical variables, the identification of a sub-group of young individuals with telomere shortening. Second, we analyzed expression levels of genes associated with age or TL in a subsample of these individuals.

# MATERIALS AND METHODS

## Participants

542 individuals from the multicenter cross-sectional study, FACE-BD (FondaMental Advanced Center of Expertise for Bipolar Disorders), were included in this study. All participants were aged 16 years or older and were diagnosed with BD type I, type II or NOS (not otherwise specified) according to DSM-IV criteria (Diagnostic and Statistical Manual of Mental Disorders - fourth edition) [21]. This cohort has been previously described in details [22]. Individual were selected for having the most completed data collected on age, sex, Montgomery Asberg Depression Scale (MADRS) [23], Young Mania Rating Scale (YMRS) [24], number of mood episodes, age of onset, Body mass Index (BMI), tobacco and childhood traumatic events (assessed using the Childhood Trauma Questionnaire, CTQ) [25]. They also completed questionnaires to assess their sleep quality and rhythm with the Pittsburgh Sleep Quality Index (PSQI) [26] and the Composite Scale of Morningness (CSM) [27]. The assessment protocol was approved by the institutional review board (Comité de Protection des Personnes Ile de France IX; January 18, 2010), in accordance with the French laws for non-interventional studies and requires only an informational letter.

## qPCR for TL measurements

A blood sample was collected at inclusion for some of the individuals who participated to the BioPsy or PsyCoh study (ANR-11-IDEX-0004-02 and ANR-10-COHO-10-01). A specific informed consent form was signed by all participants to have an additional blood sample for research purpose. DNA was extracted from peripheral blood samples of 542 patients. TL was analyzed using real-time quantitative polymerase chain reaction (qPCR) carried out on a CFX384 Touch Real-Time PCR Detection System (Biorad),

as previously described [11]. To summarize, Beta-hemoglobin was used as a single-copy standard for a relative quantification of TL and all samples were assessed in triplicate. qPCRs were performed in 10 µL of reaction containing 25 ng of DNA, 300 nmol/L of each primer and 5 µL of 2x SsoAdvanced Universal SYBR Green Supermix (Biorad). T/S ratio was obtained using the delta delta Cq calculation method.

# Gene expression measurements

To measure gene expression levels, we performed RT-qPCR (reverse transcriptase - qPCR) in 129 individuals. These samples were selected on the basis of the availability of RNA Sample. Total RNA was isolated using a Maxwell 16 LEV simply RNA Purification Kit (Promega, France). For the 129 individuals, 200 ng of RNA were reverse transcribed in a final volume of 60 µL using the iScript Reverse Transcription Supermix (Bio-Rad Laboratories, France). Following the manufacturer's protocol, the reaction mix was incubated 5 min at 25 °C, 20 min at 46 °C and 1 min at 95 °C. At the end of the reaction, cDNA samples were stored at -20 °C. Based on their expression levels in blood from the GTEX portal (https://gtexportal.org/home/), 40 genes of interest implicated in the homeostasis of TL [19, 20] or whose expression levels have been previously associated with age [18, 28, 29] were selected (Supplementary 1). Six reference genes B2M, *HPRT1*, HSP90AB1, RPL30, RPS18, and *SDHA* were chosen based on previous data from our group [30] (Supplementary 1). Custom designed 384 wells Prime PCR plates were preplated (Bio-Rad laboratories, France). Quantitative PCR were performed following the manufacturer's protocol for PrimePCR Assays (Bio-Rad laboratories, France) in a final volume of 10 µL with 2X of SsoAdvanced universal Supermix (Bio-Rad laboratories, France) and 0.5 µL of cDNA sample. The thermal cycling conditions were as follow: polymerase activation step at 95 °C for 2 min followed by 40 cycles of 95 °C for 5 sec and 60 °C for 30 s for denaturation and annealing/extension step. At the end of amplification, a melting curve was used to verify the specificity of primers. These amplifications were carried out on a CFX384 Touch Real-Time PCR Detection System (Biorad laboratories, France). Delta delta Cq calculation method was used to obtain relative level of gene expression. All samples were assessed in duplicate. For eleven genes (CCDC88C, CHAF1A, CHAF1B, CRP, DCLRE1B, EXO1, EZH2, TERF2, TERT, TP53, and WRN), the RT-qPCR results did not meet the quality control criteria of the MIQE guidelines (Cq > 30 or multiple peaks melting curve) [31]. Therefore, a total of 29 genes could be analyzed. The built-in GeNorm algorithm, in CFX Maestro Software (Bio-Rad laboratories, France) was used to calculate normalization factor from the six reference genes.

# Statistical analysis

Statistical analysis were performed on JASP (Version 0.14.1), JAMOVI (Version 1.6.23.0), Prism (Version 9), and Rstudio 2 (Version 1.2.1335) softwares. When *p*-values were < 0.05, results were considered significant.

A k-means clustering analysis was used to classify the 542 subjects using age and TL (standardized values) as classification variables with the following options : maximal number of clusters set at 10, means used as center type and Lloyd-Forgy algorithm. The elbow method was used to determine the optimal number of cluster. Kruskal-Wallis followed by pairwise comparisons Dunn's tests for continuous variables (clinical variables and levels of expression) and Chi2 tests for categorical variables were performed to compare the clusters. A logistic regression was used to test for the association between the genes found significantly different between the two youngest clusters with adjustments for age, sex, BMI and tobacco. Spearman's correlations were used to investigate correlation between genes and TL in between the two youngest clusters. In this exploratory study of gene expression no correction for multiple comparison was made.

# RESULTS

## Sample characteristics

Socio-demographic and clinical characteristics of the sample (N = 542) are presented in Supplementary 2. In summary, the median age of the sample was 33 years old and there were more women than men (64.4%). BD type 1 was diagnosed among 46.4% of the participants. Almost half of the individuals (47.44%) were current tobacco smokers and a quarter (25.74%) had a past lifetime alcohol misuse. More than a third of the participants were currently treated with lithium (39.44%) and antidepressants (ATD)
sequences to the telomeres [9]. TL, as a marker of aging, has been shown to be decreased in individuals with BD compared to the unaffected or healthy population in a meta-analysis [10]. Moreover, in BD, we have recently identified a subgroup of individuals in their thirties but with a TL that was as short as the one of those individuals in their sixties [11]. However, it remains unclear whether this decrease of TL observed in individuals with BD is associated with specific clinical characteristic, even if some potential determinants have been proposed [12]. For example, long-term lithium treatment would have a positive effect on telomere attrition [13] while childhood maltreatment [14], as well as a higher number of mood episodes [15], would have a negative impact on TL. Some alterations of circadian rhythms (i.e. circadian misalignment) might also contribute to shorten telomeres, as we suggested in a recent article [16].

The identification of genes whose expression levels change with age and, independently of chronological age, predict accelerated decline of health is an expanding area of research [17]. A transcriptomic age predictor has been proposed using a whole-blood gene expression meta-analysis in 14,983 individuals of European ancestry where 1,497 genes were identified as differently expressed with chronological age [18].

Among these age-related genes, those involved in telomere structure and maintenance have been also recently reviewed [19] and their expression levels studied in the general population. For example, a recent study of the data from a whole genome sequencing of 250 Dutch family trios analyzed TL inheritance patterns and associated genetic factors. *RRM1*, a gene playing a role in cellular proliferation, was suggested to be involved in TL regulation [20]. To our knowledge, no previous study has analyzed simultaneously the expression of age and TL-related genes in individuals with BD. Such a study in individuals with BD would help identifying potential molecular determinants of telomere shortening in this population.

In this study, we first analyzed TL in individuals with BD to replicate, in a larger cohort and with a larger set of clinical variables, the identification of a sub-group of young individuals with telomere shortening. Second, we analyzed expression levels of genes associated with age or TL in a subsample of these individuals.

# MATERIALS AND METHODS

## Participants

542 individuals from the multicenter cross-sectional study, FACE-BD (FondaMental Advanced Center of Expertise for Bipolar Disorders), were included in this study. All participants were aged 16 years or older and were diagnosed with BD type I, type II or NOS (not otherwise specified) according to DSM-IV criteria (Diagnostic and Statistical Manual of Mental Disorders - fourth edition) [21]. This cohort has been previously described in details [22]. Individual were selected for having the most completed data collected on age, sex, Montgomery Asberg Depression Scale (MADRS) [23], Young Mania Rating Scale (YMRS) [24], number of mood episodes, age of onset, Body mass Index (BMI), tobacco and childhood traumatic events (assessed using the Childhood Trauma Questionnaire, CTQ) [25]. They also completed questionnaires to assess their sleep quality and rhythm with the Pittsburgh Sleep Quality Index (PSQI) [26] and the Composite Scale of Morningness (CSM) [27]. The assessment protocol was approved by the institutional review board (Comité de Protection des Personnes Ile de France IX; January 18, 2010), in accordance with the French laws for non-interventional studies and requires only an informational letter.

## qPCR for TL measurements

A blood sample was collected at inclusion for some of the individuals who participated to the BioPsy or PsyCoh study (ANR-11-IDEX-0004-02 and ANR-10-COHO-10-01). A specific informed consent form was signed by all participants to have an additional blood sample for research purpose. DNA was extracted from peripheral blood samples of 542 patients. TL was analyzed using real-time quantitative polymerase chain reaction (qPCR) carried out on a CFX384 Touch Real-Time PCR Detection System (Biorad),

as previously described [11]. To summarize, Beta-hemoglobin was used as a single-copy standard for a relative quantification of TL and all samples were assessed in triplicate. qPCRs were performed in 10 µL of reaction containing 25 ng of DNA, 300 nmol/L of each primer and 5 µL of 2x SsoAdvanced Universal SYBR Green Supermix (Biorad). T/S ratio was obtained using the delta delta Cq calculation method.

# Gene expression measurements

To measure gene expression levels, we performed RT-qPCR (reverse transcriptase - qPCR) in 129 individuals. These samples were selected on the basis of the availability of RNA Sample. Total RNA was isolated using a Maxwell 16 LEV simply RNA Purification Kit (Promega, France). For the 129 individuals, 200 ng of RNA were reverse transcribed in a final volume of 60 µL using the iScript Reverse Transcription Supermix (Bio-Rad Laboratories, France). Following the manufacturer's protocol, the reaction mix was incubated 5 min at 25 °C, 20 min at 46 °C and 1 min at 95 °C. At the end of the reaction, cDNA samples were stored at -20 °C. Based on their expression levels in blood from the GTEX portal (https://gtexportal.org/home/), 40 genes of interest implicated in the homeostasis of TL [19, 20] or whose expression levels have been previously associated with age [18, 28, 29] were selected (Supplementary 1). Six reference genes *B2M*, *HPRT1*, *HSP90AB1*, *RPL30*, *RPS18*, and *SDHA* were chosen based on previous data from our group [30] (Supplementary 1). Custom designed 384 wells Prime PCR plates were preplated (Bio-Rad laboratories, France). Quantitative PCR were performed following the manufacturer's protocol for PrimePCR Assays (Bio-Rad laboratories, France) in a final volume of 10 µL with 2X of SsoAdvanced universal Supermix (Bio-Rad laboratories, France) and 0.5 µL of cDNA sample. The thermal cycling conditions were as follows: polymerase activation step at 95 °C for 2 min followed by 40 cycles of 95 °C for 5 sec and 60 °C for 30 s for denaturation and annealing/extension step. At the end of amplification, a melting curve was used to verify the specificity of primers. These amplifications were carried out on a CFX384 Touch Real-Time PCR Detection System (Biorad laboratories, France). Delta delta Cq calculation method was used to obtain relative level of gene expression. All samples were assessed in duplicate. For eleven genes (CCDC88C, CHAF1A, CHAF1B, CRP, DCLRE1B, EXO1, EZH2, TERF2, TERT, TP53, and WRN), the RT-qPCR results did not meet the quality control criteria of the MIQE guidelines (Cq > 30 or multiple peaks melting curve) [31]. Therefore, a total of 29 genes could be analyzed. The built-in GeNorm algorithm, in CFX Maestro Software (Bio-Rad laboratories, France) was used to calculate normalization factor from the six reference genes.

# Statistical analysis

Statistical analysis were performed on JASP (Version 0.14.1), JAMOVI (Version 1.6.23.0), Prism (Version 9), and Rstudio 2 (Version 1.2.1335) softwares. When *p*-values were < 0.05, results were considered significant.

A k-means clustering analysis was used to classify the 542 subjects using age and TL (standardized values) as classification variables with the following options : maximal number of clusters set at 10, means used as center type and Lloyd-Forgy algorithm. The elbow method was used to determine the optimal number of cluster. Kruskal-Wallis followed by pairwise comparisons Dunn's tests for continuous variables (clinical variables and levels of expression) and Chi2 tests for categorical variables were performed to compare the clusters. A logistic regression was used to test for the association between the genes found significantly different between the two youngest clusters with adjustments for age, sex, BMI and tobacco. Spearman's correlations were used to investigate correlation between genes and TL in between the two youngest clusters. In this exploratory study of gene expression no correction for multiple comparison was made.

# RESULTS

## Sample characteristics

Socio-demographic and clinical characteristics of the sample (N = 542) are presented in Supplementary 2. In summary, the median age of the sample was 33 years old and there were more women than men (64.4%). BD type 1 was diagnosed among 46.4% of the participants. Almost half of the individuals (47.44%) were current tobacco smokers and a quarter (25.74%) had a past lifetime alcohol misuse. More than a third of the participants were currently treated with lithium (39.44%) and antidepressants (ATD)
<figure><img src="image_3.png" /><figcaption>**Fig. 1** Clustering of 542 BD patients with age and TL as classification variables. Standardized values of age and TL were used to create clusters.</figcaption></figure>

(34.06%). Almost half of them were currently treated with atypical antipsychotic agents (APA) (44.69%) or with anticonvulsant agents (ACA) (50%). The median childhood trauma score of this sample was 38 and the median relative TL was 2.36. TL was negatively correlated with age in the sample (Spearman's rho = -0.236, *p* < 0.001).

## Clustering analysis with age and TL

The elbow method identified three subgroups of individuals (Supplementary 3a). As shown in Fig. 1, the first cluster named "Elderly" was a subgroup with an older age and a decreased TL (N = 138; mean age = 51.73; mean TL = 2). The second named "Young Aged" consisted of individuals with a young age but with a similar TL as the oldest cluster (N = 265; mean of age = 29.64; mean of TL = 1.96). The last subgroup named "Young" consisted in young individuals with preserved TL (N = 139; mean of age = 29.02; mean of TL = 4.36). A graph with the clusters of points resulting from this k-means clustering analysis is shown in Supplementary 3b

These three clusters were compared for all available clinical variables. Results are shown in Table 1. Most differences are found between the "Elderly" cluster and the two other clusters (i.e. youngest individuals). Individuals in "Elderly" cluster had higher

<table><caption>**Table 1.** Comparisons between the three clusters for clinical variables.</caption><thead><tr><th>Variables</th><th>C1 Elderly (N = 138)</th><th>C2 Young Aged (N = 265)</th><th>C3 Young (N = 139)</th><th>p value global</th><th>Pairwise comparisons</th></tr></thead><tbody><tr><td>Age (years), mean ± SD</td><td>51.73 ± 7.97</td><td>29.64 ± 6.11</td><td>29.02 ± 7.09</td><td>&lt;0.001<sup>a</sup></td><td>1 &gt; 2 = 3 (p &lt; 0.001<sup>c</sup>)</td></tr><tr><td>TL, mean ± SD</td><td>2.00 ± 0.86</td><td>1.96 ± 0.71</td><td>4.36 ± 1.19</td><td>&lt;0.001<sup>a</sup></td><td>1 = 2 &lt; 3 (p &lt; 0.001<sup>c</sup>)</td></tr><tr><td>Sex (Female), %</td><td>88 (63.8%)</td><td>171 (64.5%)</td><td>90 (64.7%)</td><td>0.984<sup>b</sup></td><td></td></tr><tr><td>BMI (kg/m<sup>2</sup>), mean ± SD</td><td>26.58 ± 5.78</td><td>25.36 ± 5.12</td><td>24.92 ± 5.19</td><td>0.039<sup>a</sup></td><td>1 &gt; 3 (p = 0.014<sup>c</sup>)<br>1 = 2 (p = 0.050<sup>c</sup>)<br>2 = 3 (p = 0.384<sup>c</sup>)</td></tr><tr><td>Tobacco current smoker (yes), %</td><td>40 (29.0%)</td><td>127 (47.9%)</td><td>74 (53.2%)</td><td>&lt;0.001<sup>b</sup></td><td>1 &lt; 2 = 3 (p &lt; 0.001<sup>b</sup>)</td></tr><tr><td>Past lifetime alcohol misuse (yes), %</td><td>26 (18.8%)</td><td>66 (24.9%)</td><td>29 (20.9%)</td><td>0.384<sup>b</sup></td><td></td></tr><tr><td>Past lifetime cannabis misuse (yes), %</td><td>14 (10.1%)</td><td>72 (27.2%)</td><td>32 (23.0%)</td><td>&lt;0.001<sup>b</sup></td><td>1 &lt; 2 (p &lt; 0.001<sup>b</sup>)<br>1 &lt; 3 (p = 0.002<sup>b</sup>)<br>2 = 3 (p = 0.493<sup>b</sup>)</td></tr><tr><td>MADRS, mean ± SD</td><td>9.59 ± 8.23</td><td>10.73 ± 9.51</td><td>10.64 ± 9.50</td><td>0.767<sup>a</sup></td><td></td></tr><tr><td>YMRS, mean ± SD</td><td>2.53 ± 3.64</td><td>2.22 ± 3.51</td><td>1.93 ± 2.66</td><td>0.543<sup>a</sup></td><td></td></tr><tr><td>BD type 1 (yes), %</td><td>62 (44.9%)</td><td>115 (43.4%)</td><td>74 (53.2%)</td><td>0.158<sup>b</sup></td><td></td></tr><tr><td>Age of onset (years), mean ± SD</td><td>27.77 ± 11.12</td><td>19.90 ± 5.75</td><td>19.48 ± 4.94</td><td>&lt;0.001<sup>a</sup></td><td>1 &gt; 2 = 3 (p &lt; 0.001<sup>c</sup>)</td></tr><tr><td>Duration of illness (years), mean ± SD</td><td>23.97 ± 11.09</td><td>9.77 ± 5.99</td><td>9.56 ± 6.48</td><td>&lt;0.001<sup>a</sup></td><td>1 &gt; 2 = 3 (p &lt; 0.001<sup>c</sup>)</td></tr><tr><td>Number of mood episodes, mean ± SD</td><td>10.72 ± 9.90</td><td>7.31 ± 7.17</td><td>7.12 ± 7.22</td><td>&lt;0.001<sup>a</sup></td><td>1 &gt; 2 = 3 (p &lt; 0.001<sup>c</sup>)</td></tr><tr><td>CTQ total score, mean ± SD</td><td>42.98 ± 14.72</td><td>40.72 ± 13.25</td><td>41.53 ± 14.73</td><td>0.237<sup>a</sup></td><td></td></tr><tr><td>Lithium (yes), %</td><td>56 (45.2%)</td><td>67 (34.2%)</td><td>47 (42.3%)</td><td>0.113<sup>b</sup></td><td></td></tr><tr><td>ACA (yes), %</td><td>68 (55.3%)</td><td>101 (50.8%)</td><td>48 (42.9%)</td><td>0.157<sup>b</sup></td><td></td></tr><tr><td>APA (yes), %</td><td>50 (36.2%)</td><td>85 (32.1%)</td><td>50 (36.0%)</td><td>0.716<sup>b</sup></td><td></td></tr><tr><td>ATD (yes), %</td><td>45 (32.6%)</td><td>63 (23.8%)</td><td>33 (23.7%)</td><td>0.566<sup>b</sup></td><td></td></tr><tr><td>Total number of medications, mean ± SD</td><td>2.16 ± 1.14</td><td>1.87 ± 1.03</td><td>1.82 ± 0.96</td><td>0.013<sup>a</sup></td><td>1 &gt; 2 = 3 (p = 0.01<sup>c</sup>)</td></tr><tr><td>PSQI, mean ± SD</td><td>6.75 ± 3.90</td><td>7.00 ± 3.45</td><td>6.61 ± 3.38</td><td>0.398<sup>a</sup></td><td></td></tr><tr><td>CSM, mean ± SD</td><td>37.94 ± 7.47</td><td>32.94 ± 7.41</td><td>32.60 ± 7.10</td><td>&lt;0.001<sup>a</sup></td><td>1 &gt; 2 = 3 (p &lt; 0.001<sup>c</sup>)</td></tr></tbody></table>

TL Telomere Length, BMI Body Mass Index, MADRS Montgomery Asberg Depression Rating Scale, YMRS Young Mania Rating Scale, BD Bipolar Disorder, CTQ Childhood Trauma Questionnaire, AC Anticonvulsant agent, APA Atypical Antipsychotic agent, ATD Antidepressant, PSQI Pittsburgh Sleep Quality Index, CSM Composite Scale of Morningness.

<footnote><sup>a</sup>Kruskal-Wallis test; <sup>b</sup>chi2 test; <sup>c</sup>Dunn's test.</footnote>
<figure><img src="image_4.png" /><figcaption>**Fig. 1** Clustering of 542 BD patients with age and TL as classification variables. Standardized values of age and TL were used to create clusters.</figcaption></figure>

(34.06%). Almost half of them were currently treated with atypical antipsychotic agents (APA) (44.69%) or with anticonvulsant agents (ACA) (50%). The median childhood trauma score of this sample was 38 and the median relative TL was 2.36. TL was negatively correlated with age in the sample (Spearman's rho = -0.236, *p* < 0.001).

## Clustering analysis with age and TL

The elbow method identified three subgroups of individuals (Supplementary 3a). As shown in Fig. 1, the first cluster named "Elderly" was a subgroup with an older age and a decreased TL (N = 138; mean age = 51.73; mean TL = 2). The second named "Young Aged" consisted of individuals with a young age but with a similar TL as the oldest cluster (N = 265; mean of age = 29.64; mean of TL = 1.96). The last subgroup named "Young" consisted in young individuals with preserved TL (N = 139; mean of age = 29.02; mean of TL = 4.36). A graph with the clusters of points resulting from this k-means clustering analysis is shown in Supplementary 3b

These three clusters were compared for all available clinical variables. Results are shown in Table 1. Most differences are found between the "Elderly" cluster and the two other clusters (i.e. youngest individuals). Individuals in "Elderly" cluster had higher

<table><caption>**Table 1.** Comparisons between the three clusters for clinical variables.</caption><thead><tr><th>Variables</th><th>C1 Elderly (N = 138)</th><th>C2 Young Aged (N = 265)</th><th>C3 Young (N = 139)</th><th>p value global</th><th>Pairwise comparisons</th></tr></thead><tbody><tr><td>Age (years), mean ± SD</td><td>51.73 ± 7.97</td><td>29.64 ± 6.11</td><td>29.02 ± 7.09</td><td>&lt;0.001<sup>a</sup></td><td>1 &gt; 2 = 3 (p &lt; 0.001<sup>c</sup>)</td></tr><tr><td>TL, mean ± SD</td><td>2.00 ± 0.86</td><td>1.96 ± 0.71</td><td>4.36 ± 1.19</td><td>&lt;0.001<sup>a</sup></td><td>1 = 2 &lt; 3 (p &lt; 0.001<sup>c</sup>)</td></tr><tr><td>Sex (Female), %</td><td>88 (63.8%)</td><td>171 (64.5%)</td><td>90 (64.7%)</td><td>0.984<sup>b</sup></td><td></td></tr><tr><td>BMI (kg/m<sup>2</sup>), mean ± SD</td><td>26.58 ± 5.78</td><td>25.36 ± 5.12</td><td>24.92 ± 5.19</td><td>0.039<sup>a</sup></td><td>1 &gt; 3 (p = 0.014<sup>c</sup>)<br>1 = 2 (p = 0.050<sup>c</sup>)<br>2 = 3 (p = 0.384<sup>c</sup>)</td></tr><tr><td>Tobacco current smoker (yes), %</td><td>40 (29.0%)</td><td>127 (47.9%)</td><td>74 (53.2%)</td><td>&lt;0.001<sup>b</sup></td><td>1 &lt; 2 = 3 (p &lt; 0.001<sup>b</sup>)</td></tr><tr><td>Past lifetime alcohol misuse (yes), %</td><td>26 (18.8%)</td><td>66 (24.9%)</td><td>29 (20.9%)</td><td>0.384<sup>b</sup></td><td></td></tr><tr><td>Past lifetime cannabis misuse (yes), %</td><td>14 (10.1%)</td><td>72 (27.2%)</td><td>32 (23.0%)</td><td>&lt;0.001<sup>b</sup></td><td>1 &lt; 2 (p &lt; 0.001<sup>b</sup>)<br>1 &lt; 3 (p = 0.002<sup>b</sup>)<br>2 = 3 (p = 0.493<sup>b</sup>)</td></tr><tr><td>MADRS, mean ± SD</td><td>9.59 ± 8.23</td><td>10.73 ± 9.51</td><td>10.64 ± 9.50</td><td>0.767<sup>a</sup></td><td></td></tr><tr><td>YMRS, mean ± SD</td><td>2.53 ± 3.64</td><td>2.22 ± 3.51</td><td>1.93 ± 2.66</td><td>0.543<sup>a</sup></td><td></td></tr><tr><td>BD type 1 (yes), %</td><td>62 (44.9%)</td><td>115 (43.4%)</td><td>74 (53.2%)</td><td>0.158<sup>b</sup></td><td></td></tr><tr><td>Age of onset (years), mean ± SD</td><td>27.77 ± 11.12</td><td>19.90 ± 5.75</td><td>19.48 ± 4.94</td><td>&lt;0.001<sup>a</sup></td><td>1 &gt; 2 = 3 (p &lt; 0.001<sup>c</sup>)</td></tr><tr><td>Duration of illness (years), mean ± SD</td><td>23.97 ± 11.09</td><td>9.77 ± 5.99</td><td>9.56 ± 6.48</td><td>&lt;0.001<sup>a</sup></td><td>1 &gt; 2 = 3 (p &lt; 0.001<sup>c</sup>)</td></tr><tr><td>Number of mood episodes, mean ± SD</td><td>10.72 ± 9.90</td><td>7.31 ± 7.17</td><td>7.12 ± 7.22</td><td>&lt;0.001<sup>a</sup></td><td>1 &gt; 2 = 3 (p &lt; 0.001<sup>c</sup>)</td></tr><tr><td>CTQ total score, mean ± SD</td><td>42.98 ± 14.72</td><td>40.72 ± 13.25</td><td>41.53 ± 14.73</td><td>0.237<sup>a</sup></td><td></td></tr><tr><td>Lithium (yes), %</td><td>56 (45.2%)</td><td>67 (34.2%)</td><td>47 (42.3%)</td><td>0.113<sup>b</sup></td><td></td></tr><tr><td>ACA (yes), %</td><td>68 (55.3%)</td><td>101 (50.8%)</td><td>48 (42.9%)</td><td>0.157<sup>b</sup></td><td></td></tr><tr><td>APA (yes), %</td><td>50 (36.2%)</td><td>85 (32.1%)</td><td>50 (36.0%)</td><td>0.716<sup>b</sup></td><td></td></tr><tr><td>ATD (yes), %</td><td>45 (32.6%)</td><td>63 (23.8%)</td><td>33 (23.7%)</td><td>0.566<sup>b</sup></td><td></td></tr><tr><td>Total number of medications, mean ± SD</td><td>2.16 ± 1.14</td><td>1.87 ± 1.03</td><td>1.82 ± 0.96</td><td>0.013<sup>a</sup></td><td>1 &gt; 2 = 3 (p = 0.01<sup>c</sup>)</td></tr><tr><td>PSQI, mean ± SD</td><td>6.75 ± 3.90</td><td>7.00 ± 3.45</td><td>6.61 ± 3.38</td><td>0.398<sup>a</sup></td><td></td></tr><tr><td>CSM, mean ± SD</td><td>37.94 ± 7.47</td><td>32.94 ± 7.41</td><td>32.60 ± 7.10</td><td>&lt;0.001<sup>a</sup></td><td>1 &gt; 2 = 3 (p &lt; 0.001<sup>c</sup>)</td></tr></tbody></table>

TL Telomere Length, BMI Body Mass Index, MADRS Montgomery Asberg Depression Rating Scale, YMRS Young Mania Rating Scale, BD Bipolar Disorder, CTQ Childhood Trauma Questionnaire, AC Anticonvulsant agent, APA Atypical Antipsychotic agent, ATD Antidepressant, PSQI Pittsburgh Sleep Quality Index, CSM Composite Scale of Morningness.

<footnote><sup>a</sup>Kruskal-Wallis test; <sup>b</sup>chi2 test; <sup>c</sup>Dunn's test.</footnote>
<table><caption>**Table 2.** Comparisons of genes expression between the three clusters.</caption><thead><tr><th>Genes</th><th>C1 Elderly mean ± SD (n = 16)</th><th>C2 Young aged mean ± SD (n = 68)</th><th>C3 Young mean ± SD (n = 45)</th><th>Statistic</th><th>p value<sup>a</sup></th><th>Pairwise comparisons<sup>b</sup></th></tr></thead><tbody><tr><td>CD27</td><td>0.19 ± 0.07</td><td>0.34 ± 0.29</td><td>0.43 ± 0.46</td><td>11.74</td><td>0.003</td><td>1 &lt; 2 (p = 0.04)<br>1 &lt; 3 (p &lt; 0.001)<br>2 &lt; 3 (p = 0.038)</td></tr><tr><td>LDHB</td><td>1.28 ± 0.71</td><td>1.63 ± 0.64</td><td>1.80 ± 0.68</td><td>7.11</td><td>0.029</td><td>1 = 2 (p = 0.063)<br>1 &lt; 3 (p = 0.008)<br>2 = 3 (p = 0.19)</td></tr><tr><td>LEF1</td><td>0.31 ± 0.16</td><td>0.51 ± 0.48</td><td>0.60 ± 0.48</td><td>7.63</td><td>0.022</td><td>1 = 2 (p = 0.118)<br>1 &lt; 3 (p = 0.008)<br>2 = 3 (p = 0.079)</td></tr><tr><td>LTB</td><td>0.52 ± 0.18</td><td>0.76 ± 0.42</td><td>0.74 ± 0.33</td><td>8.63</td><td>0.013</td><td>1 &lt; 2 (p = 0.005)<br>1 &lt; 3 (p = 0.006)<br>2 = 3 (p = 0.92)</td></tr><tr><td>MYC</td><td>0.11 ± 0.05</td><td>0.20 ± 0.24</td><td>0.30 ± 0.43</td><td>10.62</td><td>0.005</td><td>1 = 2 (p = 0.091)<br>1 &lt; 3 (p = 0.002)<br>2 &lt; 3 (p = 0.027)</td></tr><tr><td>PIK3IP1</td><td>0.20 ± 0.10</td><td>0.30 ± 0.27</td><td>0.38 ± 0.39</td><td>7.01</td><td>0.030</td><td>1 = 2 (p = 0.230)<br>1 &lt; 3 (p = 0.015)<br>2 = 3 (p = 0.053)</td></tr><tr><td>POT1</td><td>1.26 ± 0.73</td><td>1.40 ± 0.66</td><td>1.67 ± 0.68</td><td>6.43</td><td>0.040</td><td>1 = 2 (p = 0.57)<br>1 &lt; 3 (p = 0.045)<br>2 &lt; 3 (p = 0.027)</td></tr><tr><td>RPS6</td><td>1.20 ± 0.45</td><td>1.48 ± 0.44</td><td>1.53 ± 0.41</td><td>6.60</td><td>0.037</td><td>1 &lt; 2 (p = 0.046)<br>1 &lt; 3 (p = 0.010)<br>2 = 3 (p = 0.316)</td></tr></tbody></table>

BMI (p = 0.039), age of onset (p < 0.001), CSM score (p < 0.001) and total number of medications (p = 0.013), as compared to the two youngest clusters. This "Elderly" cluster also had fewest current smokers (p < 0.001) and cannabis users (p < 0.001), as compared to the two youngest clusters.

Individuals in "Young Aged" and "Young" clusters were similar in age, but differed only for their TL. Therefore, we did not identify any clinical variable that differentiated these 2 clusters.

### Associations between age-related genes and clusters

Gene expression levels of the 29 age-related genes were compared using Kruskal-Wallis analysis between the 3 clusters (Supplementary 4) and only genes that differed at a *p*-value threshold of 0.05 are presented in Table 2. Expression levels of 8 genes were found to be different for this threshold between the three clusters: CD27, LDHB, LEF1, LTB, MYC, PIK3IP1, POT1, and RPS6. For seven of the 8 genes (except LTB), "Young Aged" cluster had mean values of expression that were numerically between values observed for the "Young" cluster and for the "Elderly" cluster.

Pairwise comparisons tests showed that only three genes CD27 (p = 0.038), MYC (p = 0.027), and POT1 (p = 0.027) were significantly differentially expressed between the two clusters of young individuals (29.09 years and 29.31 years respectively) but with a different level of TL (2.12 and 4.35 respectively). For these three genes, a lower expression level was observed in the cluster with shorter TL ("Young Aged" cluster) as compared to the "Young" cluster (Table 2).

In order to compare these 2 youngest clusters and adjust for covariates such as Age, Sex, BMI, and tobacco, a logistic regression was performed for the 3 genes identified previously. As seen in Table 3, only POT1 (p = 0.037) and BMI (p = 0.019) significantly differed between the two youngest clusters.

To investigate if these three genes were correlated to telomere length in each cluster of the youngest populations, Spearman's correlations were performed. POT1 was the only gene to correlate

<table><caption>**Table 3.** Genes differentially expressed between Young Aged and Young clusters after adjustment for covariates.</caption><thead><tr><th>Variables</th><th>Beta</th><th>SE</th><th>p value*</th></tr></thead><tbody><tr><td>Intercept</td><td>0.645</td><td>1.581</td><td>0.639</td></tr><tr><td>Age</td><td>0.044</td><td>0.038</td><td>0.250</td></tr><tr><td>Sex</td><td>-0.374</td><td>0.472</td><td>0.428</td></tr><tr><td>BMI</td><td>-0.141</td><td>0.060</td><td><strong>0.019</strong></td></tr><tr><td>Current smoker</td><td>0.107</td><td>0.446</td><td>0.811</td></tr><tr><td>CD27</td><td>-0.310</td><td>1.810</td><td>0.864</td></tr><tr><td>MYC</td><td>1.226</td><td>2.069</td><td>0.553</td></tr><tr><td>POT1</td><td>0.676</td><td>0.325</td><td><strong>0.037</strong></td></tr></tbody></table>

SE Standard Error, BMI Body Mass Index.

*Logistic regression.

with telomere length significantly and this correlation was observed only in the "Young Aged" cluster (Table 4).

### DISCUSSION

The aim of this study was to analyze TL in a sample of 542 individuals with BD and to investigate gene expression levels of 29 age/TL-associated genes in a subgroup of 129 individuals. We first replicated the identification of a subgroup of young individuals, but with already degraded TL in an independent and larger cohort. In addition, we identified a decrease in POT1 gene expression in the "Young Aged" group and showed that POT1 expression level correlated with TL only in this specific subgroup.

TL is the most studied marker of aging in psychiatric diseases including BD [12]. TL is shortened in individuals with BD compared to the general population in most previous studies,

<footnote><sup>a</sup>Kruskal-Wallis test; <sup>b</sup>Dunn's test.</footnote>
<table><caption>**Table 2.** Comparisons of genes expression between the three clusters.</caption><thead><tr><th>Genes</th><th>C1 Elderly mean ± SD (n = 16)</th><th>C2 Young aged mean ± SD (n = 68)</th><th>C3 Young mean ± SD (n = 45)</th><th>Statistic</th><th>p value<sup>a</sup></th><th>Pairwise comparisons<sup>b</sup></th></tr></thead><tbody><tr><td>CD27</td><td>0.19 ± 0.07</td><td>0.34 ± 0.29</td><td>0.43 ± 0.46</td><td>11.74</td><td>0.003</td><td>1 &lt; 2 (p = 0.04)<br>1 &lt; 3 (p &lt; 0.001)<br><b>2 &lt; 3 (p = 0.038)</b></td></tr><tr><td>LDHB</td><td>1.28 ± 0.71</td><td>1.63 ± 0.64</td><td>1.80 ± 0.68</td><td>7.11</td><td>0.029</td><td>1 = 2 (p = 0.063)<br>1 &lt; 3 (p = 0.008)<br>2 = 3 (p = 0.19)</td></tr><tr><td>LEF1</td><td>0.31 ± 0.16</td><td>0.51 ± 0.48</td><td>0.60 ± 0.48</td><td>7.63</td><td>0.022</td><td>1 = 2 (p = 0.118)<br>1 &lt; 3 (p = 0.008)<br>2 = 3 (p = 0.079)</td></tr><tr><td>LTB</td><td>0.52 ± 0.18</td><td>0.76 ± 0.42</td><td>0.74 ± 0.33</td><td>8.63</td><td>0.013</td><td>1 &lt; 2 (p = 0.005)<br>1 &lt; 3 (p = 0.006)<br>2 = 3 (p = 0.92)</td></tr><tr><td>MYC</td><td>0.11 ± 0.05</td><td>0.20 ± 0.24</td><td>0.30 ± 0.43</td><td>10.62</td><td>0.005</td><td>1 = 2 (p = 0.091)<br>1 &lt; 3 (p = 0.002)<br><b>2 &lt; 3 (p = 0.027)</b></td></tr><tr><td>PIK3IP1</td><td>0.20 ± 0.10</td><td>0.30 ± 0.27</td><td>0.38 ± 0.39</td><td>7.01</td><td>0.030</td><td>1 = 2 (p = 0.230)<br>1 &lt; 3 (p = 0.015)<br>2 = 3 (p = 0.053)</td></tr><tr><td>POT1</td><td>1.26 ± 0.73</td><td>1.40 ± 0.66</td><td>1.67 ± 0.68</td><td>6.43</td><td>0.040</td><td>1 = 2 (p = 0.57)<br>1 &lt; 3 (p = 0.045)<br><b>2 &lt; 3 (p = 0.027)</b></td></tr><tr><td>RPS6</td><td>1.20 ± 0.45</td><td>1.48 ± 0.44</td><td>1.53 ± 0.41</td><td>6.60</td><td>0.037</td><td>1 &lt; 2 (p = 0.046)<br>1 &lt; 3 (p = 0.010)<br>2 = 3 (p = 0.316)</td></tr></tbody></table>

BMI (p = 0.039), age of onset (p < 0.001), CSM score (p < 0.001) and total number of medications (p = 0.013), as compared to the two youngest clusters. This "Elderly" cluster also had fewest current smokers (p < 0.001) and cannabis users (p < 0.001), as compared to the two youngest clusters.

Individuals in "Young Aged" and "Young" clusters were similar in age, but differed only for their TL. Therefore, we did not identify any clinical variable that differentiated these 2 clusters.

### Associations between age-related genes and clusters

Gene expression levels of the 29 age-related genes were compared using Kruskal-Wallis analysis between the 3 clusters (Supplementary 4) and only genes that differed at a *p*-value threshold of 0.05 are presented in Table 2. Expression levels of 8 genes were found to be different for this threshold between the three clusters: CD27, LDHB, LEF1, LTB, MYC, PIK3IP1, POT1, and RPS6. For seven of the 8 genes (except LTB), "Young Aged" cluster had mean values of expression that were numerically between values observed for the "Young" cluster and for the "Elderly" cluster.

Pairwise comparisons tests showed that only three genes CD27 (p = 0.038), MYC (p = 0.027), and POT1 (p = 0.027) were significantly differentially expressed between the two clusters of young individuals (29.09 years and 29.31 years respectively) but with a different level of TL (2.12 and 4.35 respectively). For these three genes, a lower expression level was observed in the cluster with shorter TL ("Young Aged" cluster) as compared to the "Young" cluster (Table 2).

In order to compare these 2 youngest clusters and adjust for covariates such as Age, Sex, BMI, and tobacco, a logistic regression was performed for the 3 genes identified previously. As seen in Table 3, only POT1 (p = 0.037) and BMI (p = 0.019) significantly differed between the two youngest clusters.

To investigate if these three genes were correlated to telomere length in each cluster of the youngest populations, Spearman's correlations were performed. POT1 was the only gene to correlate

<table><caption>**Table 3.** Genes differentially expressed between Young Aged and Young clusters after adjustment for covariates.</caption><thead><tr><th>Variables</th><th>Beta</th><th>SE</th><th>p value*</th></tr></thead><tbody><tr><td>Intercept</td><td>0.645</td><td>1.581</td><td>0.639</td></tr><tr><td>Age</td><td>0.044</td><td>0.038</td><td>0.250</td></tr><tr><td>Sex</td><td>-0.374</td><td>0.472</td><td>0.428</td></tr><tr><td>BMI</td><td>-0.141</td><td>0.060</td><td><b>0.019</b></td></tr><tr><td>Current smoker</td><td>0.107</td><td>0.446</td><td>0.811</td></tr><tr><td>CD27</td><td>-0.310</td><td>1.810</td><td>0.864</td></tr><tr><td>MYC</td><td>1.226</td><td>2.069</td><td>0.553</td></tr><tr><td>POT1</td><td>0.676</td><td>0.325</td><td><b>0.037</b></td></tr></tbody></table>

SE Standard Error, BMI Body Mass Index.

*Logistic regression.

with telomere length significantly and this correlation was observed only in the "Young Aged" cluster (Table 4).

### DISCUSSION

The aim of this study was to analyze TL in a sample of 542 individuals with BD and to investigate gene expression levels of 29 age/TL-associated genes in a subgroup of 129 individuals. We first replicated the identification of a subgroup of young individuals, but with already degraded TL in an independent and larger cohort. In addition, we identified a decrease in POT1 gene expression in the "Young Aged" group and showed that POT1 expression level correlated with TL only in this specific subgroup.

TL is the most studied marker of aging in psychiatric diseases including BD [12]. TL is shortened in individuals with BD compared to the general population in most previous studies,

<footnote><sup>a</sup>Kruskal-Wallis test; <sup>b</sup>Dunn's test.</footnote>
<table><caption>**Table 4.** Correlation between TL and genes (*CD27*, *MYC* and *POT1*) in the two clusters of young individuals.</caption><thead><tr><th rowspan="2">Genes</th><th colspan="2">Young (n = 45)</th><th colspan="2">Young Aged (n = 68)</th></tr><tr><th>Spearman's rho</th><th>p value</th><th>Spearman's rho</th><th>p value</th></tr></thead><tbody><tr><td>CD27</td><td>-0.108</td><td>0.478</td><td>0.018</td><td>0.881</td></tr><tr><td>MYC</td><td>-0.194</td><td>0.200</td><td>0.056</td><td>0.651</td></tr><tr><td>POT1</td><td>-0.046</td><td>0.764</td><td>0.357</td><td><strong>0.003</strong></td></tr></tbody></table>

but without any consensus about the potential clinical determinants [10, 12]. As expected, TL was negatively correlated with age in our sample [8].

First, we confirmed our previous finding of a subgroup of young individuals with an already pronounced telomere shortening [11] in a larger and independent sample of 542 individuals. As previously found, the two clusters of young patients—but with different TL—do not differ for all tested clinical variables. These negative results contrast with previous published ones. Shortened TL has been suggested to be associated with factors such as a high early life trauma score in both general population [32, 33] in individuals with psychiatric disorders including BD [14, 34], but also with a higher number of mood episodes [15]. However, as in our previously published study [11], we did not find an association between TL shortening and these variables in our sample. Another correlation has been suggested between long-term lithium treatment and longer TL [13, 35-37], however, the duration of lithium treatment was available for only 13% of the individuals in our database. Therefore we could not test this hypothesis. In addition, the previously observed higher anticonvulsant current use in the cluster of young patients with low TL [11] was not found in this sample. We therefore cannot conclude about the clinical determinants of accelerated cellular aging in the subgroup of young individuals with shortened TL.

Second, we analyzed the expression levels of 29 genes involved in the maintenance of TL or whose expression has been associated with age [18-20, 28, 29]. However, mRNA were available only for a subset of the individuals in the cohort, thus leading to some potential power issues. Nevertheless, expression levels of three genes (*CD27*, *MYC* and *POT1*) were significantly different between the two clusters of young individuals. After adjustments with age, sex, BMI and current smoker status only BMI and *POT1* remained significantly associated with the "Young Aged" group. The observed negative association of BMI with TL in the "Young Aged" subgroup is consistent with in a meta-analysis of 87 studies in the general population [38].

*POT1* (Protection Of Telomeres 1) encodes a protein of the "shelterin complex". This complex includes six proteins POT1, TRF1, TRF2, TPP1 (ACD), TIN2 and RAP1 and participates in telomere maintenance [39]. POT1 directly interacts with telomeric DNA with the single-stranded G-rich tip where it prevents telomeres from degradation [39]. In addition, its affinity with telomeric DNA is enhanced by TPP1 binding [40] and the POT1-TPP1 complex has an essential role in telomere extension by recruiting telomerase to allow DNA synthesis of telomeres [41]. In this study, we found that *POT1* expression was decreased in the "Young Aged" cluster (shorter TL) as compared to the "Young" cluster suggesting that the formation of the POT1-TPP1 complex might be impaired in the "Young Aged" group. We hypothesize that this may result in less protection of the telomeric ends and a less effective recruitment of the telomerase. This is consistent with the increase telomere instability induced by *POT1* depletion in mice models [42, 43]. Furthermore, depletion of *POT1* has been shown to promote telomere fragility in several human cell lines [44-46] In addition, decreased *POT1* expression level has been associated with decreased TL in individuals with severe aplastic anemia [47], which is consistent with the correlation we observed specifically in the "Young Aged" cluster. To our knowledge, this is

the first evidence of *POT1* decrease of expression in BD or any other psychiatric disorders. Moreover, we found that *POT1* expression level correlated with TL only in the "Young Aged" cluster. This suggests that, in this cluster, *POT1* expression level has a high impact on TL.

This study has several strengths and limitations. First, a large cohort of 542 individuals with BD was used for TL analyzes with a large spectrum of clinical variables available. However, some variables could not be analyzed due to missing data. This was the case especially for all somatic comorbidities variables due to the limited number of patients in each group and missing data (for example 11 subjects with diabetes versus 509 without diabetes). Indeed, many somatic comorbidities associated with an accelerated aging [48, 49] are also associated with BD. These are therefore important additional factors to analyze, but comorbidities data are often not collected properly and few studies can analyze them. Second, DNA and RNA samples were available in the same cohort, allowing us to analyze both TL and gene expression levels in the same individuals. However, RNA was not available for the whole sample and this may have led to a decrease in the power of gene expression analyses. Here, we cannot exclude some false negative results. It is also important to note that this part of the study was exploratory and carried out on a subgroup of our sample. Therefore, no correction for multiple testing was made. The aim of this study was to investigate the association between TL and expression levels of age-related genes with clinical variables in individuals with BD. However, it would have been interesting to compare *POT1* mRNA levels in healthy controls samples and individuals with BD. Finally, although we have included telomerase (TERT) in the tested genes, its expression level was too low and could not be analyzed. Furthermore, as blood cells were not available in this sample, telomerase activity could not be analyzed in this study. Therefore, we can't rule out a role for telomerase (in addition to the decreased *POT1* gene expression) in the observed decrease in telomere length in the "Young Aged" group compared to the "Young" group.

In conclusion, this study confirmed the existence of a subgroup of young individuals with BD presenting with shorter TL in a larger and independent sample. None of the available clinical variables in this sample was associated with this prematurely aged phenotype. For the first time, we studied 29 genes whose expression levels are associated with age or TL in a sample of individuals with BD. We showed that a component of the shelterin complex (*POT1*) is under-expressed in the subgroup of prematurely aged patients compared to the subgroup of similar age. This downregulation is in line with previously described mechanisms of telomere instability. Therefore, our results suggest a newly described cellular mechanism in individuals with BD, that may, at least in part, contribute to telomere shortening. Further studies in independent and larger samples are required to replicate the downregulation of *POT1* and provide a better understanding of the mechanisms at stake in cellular aging in individuals with BD.

DATA AVAILABILITY

The RT-qPCR data that support the findings of this study are available from the corresponding author upon reasonable request.
<table><caption>**Table 4.** Correlation between TL and genes (*CD27*, *MYC* and *POT1*) in the two clusters of young individuals.</caption><thead><tr><th rowspan="2">Genes</th><th colspan="2">Young (n = 45)</th><th colspan="2">Young Aged (n = 68)</th></tr><tr><th>Spearman's rho</th><th>p value</th><th>Spearman's rho</th><th>p value</th></tr></thead><tbody><tr><td>CD27</td><td>-0.108</td><td>0.478</td><td>0.018</td><td>0.881</td></tr><tr><td>MYC</td><td>-0.194</td><td>0.200</td><td>0.056</td><td>0.651</td></tr><tr><td>POT1</td><td>-0.046</td><td>0.764</td><td>0.357</td><td><strong>0.003</strong></td></tr></tbody></table>

but without any consensus about the potential clinical determinants [10, 12]. As expected, TL was negatively correlated with age in our sample [8].

First, we confirmed our previous finding of a subgroup of young individuals with an already pronounced telomere shortening [11] in a larger and independent sample of 542 individuals. As previously found, the two clusters of young patients—but with different TL—do not differ for all tested clinical variables. These negative results contrast with previous published ones. Shortened TL has been suggested to be associated with factors such as a high early life trauma score in both general population [32, 33] in individuals with psychiatric disorders including BD [14, 34], but also with a higher number of mood episodes [15]. However, as in our previously published study [11], we did not find an association between TL shortening and these variables in our sample. Another correlation has been suggested between long-term lithium treatment and longer TL [13, 35-37], however, the duration of lithium treatment was available for only 13% of the individuals in our database. Therefore we could not test this hypothesis. In addition, the previously observed higher anticonvulsant current use in the cluster of young patients with low TL [11] was not found in this sample. We therefore cannot conclude about the clinical determinants of accelerated cellular aging in the subgroup of young individuals with shortened TL.

Second, we analyzed the expression levels of 29 genes involved in the maintenance of TL or whose expression has been associated with age [18-20, 28, 29]. However, mRNA were available only for a subset of the individuals in the cohort, thus leading to some potential power issues. Nevertheless, expression levels of three genes (*CD27*, *MYC* and *POT1*) were significantly different between the two clusters of young individuals. After adjustments with age, sex, BMI and current smoker status only BMI and *POT1* remained significantly associated with the "Young Aged" group. The observed negative association of BMI with TL in the "Young Aged" subgroup is consistent with in a meta-analysis of 87 studies in the general population [38].

*POT1* (Protection Of Telomeres 1) encodes a protein of the "shelterin complex". This complex includes six proteins POT1, TRF1, TRF2, TPP1 (ACD), TIN2 and RAP1 and participates in telomere maintenance [39]. POT1 directly interacts with telomeric DNA with the single-stranded G-rich tip where it prevents telomeres from degradation [39]. In addition, its affinity with telomeric DNA is enhanced by TPP1 binding [40] and the POT1-TPP1 complex has an essential role in telomere extension by recruiting telomerase to allow DNA synthesis of telomeres [41]. In this study, we found that *POT1* expression was decreased in the "Young Aged" cluster (shorter TL) as compared to the "Young" cluster suggesting that the formation of the POT1-TPP1 complex might be impaired in the "Young Aged" group. We hypothesize that this may result in less protection of the telomeric ends and a less effective recruitment of the telomerase. This is consistent with the increase telomere instability induced by *POT1* depletion in mice models [42, 43]. Furthermore, depletion of *POT1* has been shown to promote telomere fragility in several human cell lines [44-46] In addition, decreased *POT1* expression level has been associated with decreased TL in individuals with severe aplastic anemia [47], which is consistent with the correlation we observed specifically in the "Young Aged" cluster. To our knowledge, this is

the first evidence of *POT1* decrease of expression in BD or any other psychiatric disorders. Moreover, we found that *POT1* expression level correlated with TL only in the "Young Aged" cluster. This suggests that, in this cluster, *POT1* expression level has a high impact on TL.

This study has several strengths and limitations. First, a large cohort of 542 individuals with BD was used for TL analyzes with a large spectrum of clinical variables available. However, some variables could not be analyzed due to missing data. This was the case especially for all somatic comorbidities variables due to the limited number of patients in each group and missing data (for example 11 subjects with diabetes versus 509 without diabetes). Indeed, many somatic comorbidities associated with an accelerated aging [48, 49] are also associated with BD. These are therefore important additional factors to analyze, but comorbidities data are often not collected properly and few studies can analyze them. Second, DNA and RNA samples were available in the same cohort, allowing us to analyze both TL and gene expression levels in the same individuals. However, RNA was not available for the whole sample and this may have led to a decrease in the power of gene expression analyses. Here, we cannot exclude some false negative results. It is also important to note that this part of the study was exploratory and carried out on a subgroup of our sample. Therefore, no correction for multiple testing was made. The aim of this study was to investigate the association between TL and expression levels of age-related genes with clinical variables in individuals with BD. However, it would have been interesting to compare *POT1* mRNA levels in healthy controls samples and individuals with BD. Finally, although we have included telomerase (TERT) in the tested genes, its expression level was too low and could not be analyzed. Furthermore, as blood cells were not available in this sample, telomerase activity could not be analyzed in this study. Therefore, we can't rule out a role for telomerase (in addition to the decreased *POT1* gene expression) in the observed decrease in telomere length in the "Young Aged" group compared to the "Young" group.

In conclusion, this study confirmed the existence of a subgroup of young individuals with BD presenting with shorter TL in a larger and independent sample. None of the available clinical variables in this sample was associated with this prematurely aged phenotype. For the first time, we studied 29 genes whose expression levels are associated with age or TL in a sample of individuals with BD. We showed that a component of the shelterin complex (*POT1*) is under-expressed in the subgroup of prematurely aged patients compared to the subgroup of similar age. This downregulation is in line with previously described mechanisms of telomere instability. Therefore, our results suggest a newly described cellular mechanism in individuals with BD, that may, at least in part, contribute to telomere shortening. Further studies in independent and larger samples are required to replicate the downregulation of *POT1* and provide a better understanding of the mechanisms at stake in cellular aging in individuals with BD.

DATA AVAILABILITY

The RT-qPCR data that support the findings of this study are available from the corresponding author upon reasonable request.
REFERENCES

1. Kessing LV, Vradi E, Andersen PK. Life expectancy in bipolar disorder. Bipolar Disord. 2015;17:543–8.

2. Chan JKN, Tong CHY, Wong CSM, Chen EYH, Chang WC. Life expectancy and years of potential life lost in bipolar disorder: systematic review and meta-analysis. Br J Psychiatry: J Ment Sci. 2022;221:567–76.

3. Kessing LV, Vradi E, McIntyre RS, Andersen PK. Causes of decreased life expectancy over the life span in bipolar disorder. J Affect Disord. 2015;180:142–7.

4. McIntyre RS, Berk M, Brietzke E, Goldstein BI, López-Jaramillo C, Kessing LV, et al. Bipolar disorders. Lancet. 2020;396:1841–56.

5. Yuan M, Xiao Z-L, Zhou H-Y, Rao W, Huang G, Nie H-B, et al. Bipolar disorder and the risk for stroke incidence and mortality: a meta-analysis. Neurological Sci. 2022;43:467–76.

6. Rizzo LB. The theory of bipolar disorder as an illness of accelerated aging: Implications for clinical care and research. Neuroscience Biobehav Rev. 2014;42:157–69.

7. Jylhävä J, Pedersen NL, Hägg S. Biological age predictors. EBioMedicine. 2017;21:29–36.

8. Aubert G, Lansdorp PM. Telomeres and aging. Physiol Rev. 2008;88:23.

9. Blackburn EH, Epel ES, Lin J. Human telomere biology: a contributory and interactive factor in aging, disease risks, and protection. Science. 2015;350:1193–8.

10. Huang Y-C, Wang L-J, Tseng P-T, Hung C-F, Lin P-Y. Leukocyte telomere length in patients with bipolar disorder: An updated meta-analysis and subgroup analysis by mood status. Psychiatry Res. 2018;270:41–9.

11. Spano L, Etain B, Meyrel M, Hennion V, Gross G, Laplanche J-L, et al. Telomere length and mitochondrial DNA copy number in bipolar disorder: identification of a subgroup of young individuals with accelerated cellular aging. Transl Psychiatry. 2022;12:135.

12. Fries GR, Zamzow MJ, Andrews T, Pink O, Scaini G, Quevedo J. Accelerated aging in bipolar disorder: A comprehensive review of molecular findings and their clinical implications. Neurosci Biobehav Rev. 2020;112:107–16.

13. Martinsson L, Wei Y, Xu D, Melas PA, Mathé AA, Schalling M, et al. Long-term lithium treatment in bipolar disorder is associated with longer leukocyte telomeres. Transl Psychiatry. 2013;3:e261–e261.

14. Aas M, Elvsåshagen T, Westlye LT, Kaufmann T, Athanasiu L, Djurovic S, et al. Telomere length is associated with childhood trauma in patients with severe mental disorders. Transl Psychiatry. 2019;9:97.

15. Elvsåshagen T, Vera E, Boen E, Bratlie J, Andreassen OA, Josefsen D, et al. The load of short telomeres is increased and associated with lifetime number of depressive episodes in bipolar II disorder. J Affect Disord. 2011;135:43–50.

16. Spano L, Hennion V, Marie-Claire C, Bellivier F, Scott J, Etain B. Associations between circadian misalignment and telomere length in BD: an actigraphy study. Int J Bipolar Disord. 2022;10:14.

17. Cardoso AL, Fernandes A, Aguilar-Pimentel JA, de Angelis MH, Guedes JR, Brito MA, et al. Towards frailty biomarkers: Candidates from genes and pathways regulated in aging and age-related diseases. Ageing Res Rev. 2018;47:214–77.

18. Peters MJ, Joehanes R, Pilling LC, Schurmann C, Conneely KN, Powell J, et al. The transcriptional landscape of age in human peripheral blood. Nat Commun. 2015;6:8570.

19. Aksenova AY, Mirkin SM. At the Beginning of the End and in the Middle of the Beginning: Structure and Maintenance of Telomeric DNA Repeats and Interstitial Telomeric Sequences. Genes. 2019;10:118.

20. Nersisyan L, Nikoghosyan M, The Genome of the Netherlands consortium, Francioli LC, Menelaou A, Pulit SL, et al. WGS-based telomere length analysis in Dutch family trios implicates stronger maternal inheritance and a role for RRM1 gene. Sci Rep. 2019;9:18758.

21. Bell CC. DSM-IV: Diagnostic and Statistical Manual of Mental Disorders. JAMA: J Am Med Assoc. 1994;272:828.

22. Grillault Laroche D, Godin O, Belzeaux M, R'Bailara K, Loftus J, Courtet P, et al. Association between childhood maltreatment and the clinical course of bipolar disorders: A survival analysis of mood recurrences. Acta Psychiatr Scandinavica. 2022;145:373–83.

23. Montgomery SA, Åsberg M. A new depression scale designed to be sensitive to change. Br J Psychiatry. 1979;134:382–9.

24. Bech P, Rafaelsen OJ, Kramp P, Bolwig TG. The mania rating scale: scale construction and inter-observer agreement. Neuropharmacology. 1978;17:430–1.

25. Bernstein DP, Fink L. Childhood Trauma Questionnaire: A Retrospective Self-Report: Manual. Orlando: Psychological Corporation; 1998.

26. Buysse DJ, Reynolds CF, Monk TH, Berman SR, Kupfer DJ. The Pittsburgh Sleep Quality Index: a new instrument for psychiatric practice and research. Psychiatry Res. 1989;28:193–213.

27. Smith CS, Reilly C, Midkiff K. Evaluation of three circadian rhythm questionnaires with suggestions for an improved measure of morningness. J Appl Psychol. 1989;74:728–38.

28. Lin C-W, Chang L-C, Ma T, Oh H, French B, Puralewski R, et al. Older molecular brain age in severe mental illness. Mol Psychiatry. 2021;26:3646–56.

29. Craig T, Smelick C, Tacutu R, Wuttke D, Wood SH, Stanley H, et al. The Digital Ageing Atlas: integrating the diversity of age-related changes into a unified resource. Nucleic Acids Res. 2015;43:D873–878.

30. Curis E, Nepost C, Grillault Laroche D, Courtin C, Laplanche J-L, Etain B, et al. Selecting reference genes in RT-qPCR based on equivalence tests: a network based approach. Sci Rep. 2019;9:16231.

31. Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, Kubista M, et al. The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin Chem. 2009;55:611–22.

32. Li Z, He Y, Wang D, Tang J, Chen X. Association between childhood trauma and accelerated telomere erosion in adulthood: a meta-analytic study. J Psychiatr Res. 2017;93:64–71.

33. Woo JMP, Parks CG, Hyde EE, Auer PL, Simanek AM, Konkel RH, et al. Early life trauma and adult leucocyte telomere length. Psychoneuroendocrinology. 2022;144:105876.

34. Riley G, Perrin M, Vaez-Azizi LM, Ruby E, Goetz RR, Draxler R, et al. Telomere length and early trauma in schizophrenia. Schizophr Res. 2018;199:426–30.

35. Squassina A, Pisanu C, Congiu D, Caria P, Frau D, Niola P, et al. Leukocyte telomere length positively correlates with duration of lithium treatment in bipolar disorder patients. Eur Neuropsychopharmacol. 2016;26:1241–7.

36. Pisanu C, Congiu D, Manchia M, Caria P, Cocco C, Dettori T, et al. Differences in telomere length between patients with bipolar disorder and controls are influenced by lithium treatment. Pharmacogenomics. 2020;21:533–40.

37. Coutts F, Palmos AB, Duarte RRR, de Jong S, Lewis CM, Dima D, et al. The polygenic nature of telomere length and the anti-ageing properties of lithium. Neuropsychopharmacol. 2019;44:757–65.

38. Gielen M, Hageman GJ, Antoniou EE, Nordfjall K, Mangino M, Balasubramanyam M, et al. Body mass index is negatively associated with telomere length: a collaborative cross-sectional meta-analysis of 87 observational studies. Am J Clin Nutr. 2018;108:453–75.

39. Smith EM, Pendlebury DF, Nandakumar J. Structural biology of telomeres and telomerase. Cell Mol Life Sci. 2020;77:61–79.

40. Aramburu T, Plucinsky S, Skordalakes E. POT1-TPP1 telomere length regulation and disease. Comput Struct Biotechnol J. 2020;18:1939–46.

41. Wang F, Podell ER, Zaug AJ, Yang Y, Baciu P, Cech TR, et al. The POT1-TPP1 telomere complex is a telomerase processivity factor. Nature. 2007;445:506–10.

42. He H, Wang Y, Guo X, Ramchandani S, Ma J, Shen M-F, et al. Pot1b deletion and telomerase haploinsufficiency in mice initiate an ATR-dependent DNA damage response and elicit phenotypes resembling dyskeratosis congenita. Mol Cell Biol. 2009;29:229–40.

43. Wu L, Multani AS, He H, Cosme-Blanco W, Deng Y, Deng JM, et al. Pot1 deficiency initiates DNA damage checkpoint activation and aberrant homologous recombination at telomeres. Cell. 2006;126:49–62.

44. Wang H, Ni J, Guo X, Zhou T, Ma X, Xue J, et al. Shelterin differentially respond to oxidative stress induced by TiO2-NPs and regulate telomere length in human hepatocytes and hepatocarcinoma cells in vitro. Biochem Biophys Res Commun. 2018;503:697–702.

45. Chuenwisad K, More-Krong P, Tubsaeng P, Chotechuang N, Srisa-Art M, Storer RJ, et al. Premature Senescence and Telomere Shortening Induced by Oxidative Stress From Oxalate, Calcium Oxalate Monohydrate, and Urine From Patients With Calcium Oxalate Nephrolithiasis. Front Immunol. 2021;12:696486.

46. Glousker G, Briod A-S, Quadroni M, Lingner J. Human shelterin protein POT1 prevents severe telomere instability induced by homology-directed DNA repair. EMBO J. 2020;39:e104500.

47. Wang T, Mei S, Fu R, Wang H, Shao Z. Expression of shelterin component POT1 is associated with decreased telomere length and immunity condition in humans with severe aplastic anemia. J Immunol Res. 2014;2014:439530.

48. Shalev I, Entringer S, Wadhwa PD, Wolkowitz OM, Puterman E, Lin J, et al. Stress and telomere biology: a lifespan perspective. Psychoneuroendocrinology. 2013;38:1835–42.

49. Patel J, Baptiste BA, Kim E, Hussain M, Croteau DL, Bohr VA. DNA damage and mitochondria in cancer and aging. Carcinogenesis. 2020;41:1625–34.

ACKNOWLEDGEMENTS

The authors thank the FondaMental Foundation (www-fondation-fondamental.org) for its support for the infrastructure of Bipolar Expert Centers. We thank the clinicians and nursing staff of the Center for Expertise in BD for the recruitment and clinical characterization of individuals with BD. The biobank is hosted at the Plateforme de Ressources Biologique, Hôpital Henri Mondor, Créteil (Dr C. Barau). Dr. C. Marie-Claire is supported by the Centre National de la Recherche Scientifique.
REFERENCES

1. Kessing LV, Vradi E, Andersen PK. Life expectancy in bipolar disorder. Bipolar Disord. 2015;17:543–8.

2. Chan JKN, Tong CHY, Wong CSM, Chen EYH, Chang WC. Life expectancy and years of potential life lost in bipolar disorder: systematic review and meta-analysis. Br J Psychiatry: J Ment Sci. 2022;221:567–76.

3. Kessing LV, Vradi E, McIntyre RS, Andersen PK. Causes of decreased life expectancy over the life span in bipolar disorder. J Affect Disord. 2015;180:142–7.

4. McIntyre RS, Berk M, Brietzke E, Goldstein BI, López-Jaramillo C, Kessing LV, et al. Bipolar disorders. Lancet. 2020;396:1841–56.

5. Yuan M, Xiao Z-L, Zhou H-Y, Rao W, Huang G, Nie H-B, et al. Bipolar disorder and the risk for stroke incidence and mortality: a meta-analysis. Neurological Sci. 2022;43:467–76.

6. Rizzo LB. The theory of bipolar disorder as an illness of accelerated aging: Implications for clinical care and research. Neuroscience Biobehav Rev. 2014;42:157–69.

7. Jylhävä J, Pedersen NL, Hägg S. Biological age predictors. EBioMedicine. 2017;21:29–36.

8. Aubert G, Lansdorp PM. Telomeres and aging. Physiol Rev. 2008;88:23.

9. Blackburn EH, Epel ES, Lin J. Human telomere biology: a contributory and interactive factor in aging, disease risks, and protection. Science. 2015;350:1193–8.

10. Huang Y-C, Wang L-J, Tseng P-T, Hung C-F, Lin P-Y. Leukocyte telomere length in patients with bipolar disorder: An updated meta-analysis and subgroup analysis by mood status. Psychiatry Res. 2018;270:41–9.

11. Spano L, Etain B, Meyrel M, Hennion V, Gross G, Laplanche J-L, et al. Telomere length and mitochondrial DNA copy number in bipolar disorder: identification of a subgroup of young individuals with accelerated cellular aging. Transl Psychiatry. 2022;12:135.

12. Fries GR, Zamzow MJ, Andrews T, Pink O, Scaini G, Quevedo J. Accelerated aging in bipolar disorder: A comprehensive review of molecular findings and their clinical implications. Neurosci Biobehav Rev. 2020;112:107–16.

13. Martinsson L, Wei Y, Xu D, Melas PA, Mathé AA, Schalling M, et al. Long-term lithium treatment in bipolar disorder is associated with longer leukocyte telomeres. Transl Psychiatry. 2013;3:e261–e261.

14. Aas M, Elvsåshagen T, Westlye LT, Kaufmann T, Athanasiu L, Djurovic S, et al. Telomere length is associated with childhood trauma in patients with severe mental disorders. Transl Psychiatry. 2019;9:97.

15. Elvsåshagen T, Vera E, Boen E, Bratlie J, Andreassen OA, Josefsen D, et al. The load of short telomeres is increased and associated with lifetime number of depressive episodes in bipolar II disorder. J Affect Disord. 2011;135:43–50.

16. Spano L, Hennion V, Marie-Claire C, Bellivier F, Scott J, Etain B. Associations between circadian misalignment and telomere length in BD: an actigraphy study. Int J Bipolar Disord. 2022;10:14.

17. Cardoso AL, Fernandes A, Aguilar-Pimentel JA, de Angelis MH, Guedes JR, Brito MA, et al. Towards frailty biomarkers: Candidates from genes and pathways regulated in aging and age-related diseases. Ageing Res Rev. 2018;47:214–77.

18. Peters MJ, Joehanes R, Pilling LC, Schurmann C, Conneely KN, Powell J, et al. The transcriptional landscape of age in human peripheral blood. Nat Commun. 2015;6:8570.

19. Aksenova AY, Mirkin SM. At the Beginning of the End and in the Middle of the Beginning: Structure and Maintenance of Telomeric DNA Repeats and Interstitial Telomeric Sequences. Genes. 2019;10:118.

20. Nersisyan L, Nikoghosyan M, The Genome of the Netherlands consortium, Francioli LC, Menelaou A, Pulit SL, et al. WGS-based telomere length analysis in Dutch family trios implicates stronger maternal inheritance and a role for RRM1 gene. Sci Rep. 2019;9:18758.

21. Bell CC. DSM-IV: Diagnostic and Statistical Manual of Mental Disorders. JAMA: J Am Med Assoc. 1994;272:828.

22. Grillault Laroche D, Godin O, Belzeaux M, R'Bailara K, Loftus J, Courtet P, et al. Association between childhood maltreatment and the clinical course of bipolar disorders: A survival analysis of mood recurrences. Acta Psychiatr Scandinavica. 2022;145:373–83.

23. Montgomery SA, Åsberg M. A new depression scale designed to be sensitive to change. Br J Psychiatry. 1979;134:382–9.

24. Bech P, Rafaelsen OJ, Kramp P, Bolwig TG. The mania rating scale: scale construction and inter-observer agreement. Neuropharmacology. 1978;17:430–1.

25. Bernstein DP, Fink L. Childhood Trauma Questionnaire: A Retrospective Self-Report: Manual. Orlando: Psychological Corporation; 1998.

26. Buysse DJ, Reynolds CF, Monk TH, Berman SR, Kupfer DJ. The Pittsburgh Sleep Quality Index: a new instrument for psychiatric practice and research. Psychiatry Res. 1989;28:193–213.

27. Smith CS, Reilly C, Midkiff K. Evaluation of three circadian rhythm questionnaires with suggestions for an improved measure of morningness. J Appl Psychol. 1989;74:728–38.

28. Lin C-W, Chang L-C, Ma T, Oh H, French B, Puralewski R, et al. Older molecular brain age in severe mental illness. Mol Psychiatry. 2021;26:3646–56.

29. Craig T, Smelick C, Tacutu R, Wuttke D, Wood SH, Stanley H, et al. The Digital Ageing Atlas: integrating the diversity of age-related changes into a unified resource. Nucleic Acids Res. 2015;43:D873–878.

30. Curis E, Nepost C, Grillault Laroche D, Courtin C, Laplanche J-L, Etain B, et al. Selecting reference genes in RT-qPCR based on equivalence tests: a network based approach. Sci Rep. 2019;9:16231.

31. Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, Kubista M, et al. The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin Chem. 2009;55:611–22.

32. Li Z, He Y, Wang D, Tang J, Chen X. Association between childhood trauma and accelerated telomere erosion in adulthood: a meta-analytic study. J Psychiatr Res. 2017;93:64–71.

33. Woo JMP, Parks CG, Hyde EE, Auer PL, Simanek AM, Konkel RH, et al. Early life trauma and adult leucocyte telomere length. Psychoneuroendocrinology. 2022;144:105876.

34. Riley G, Perrin M, Vaez-Azizi LM, Ruby E, Goetz RR, Draxler R, et al. Telomere length and early trauma in schizophrenia. Schizophr Res. 2018;199:426–30.

35. Squassina A, Pisanu C, Congiu D, Caria P, Frau D, Niola P, et al. Leukocyte telomere length positively correlates with duration of lithium treatment in bipolar disorder patients. Eur Neuropsychopharmacol. 2016;26:1241–7.

36. Pisanu C, Congiu D, Manchia M, Caria P, Cocco C, Dettori T, et al. Differences in telomere length between patients with bipolar disorder and controls are influenced by lithium treatment. Pharmacogenomics. 2020;21:533–40.

37. Coutts F, Palmos AB, Duarte RRR, de Jong S, Lewis CM, Dima D, et al. The polygenic nature of telomere length and the anti-ageing properties of lithium. Neuropsychopharmacol. 2019;44:757–65.

38. Gielen M, Hageman GJ, Antoniou EE, Nordfjall K, Mangino M, Balasubramanyam M, et al. Body mass index is negatively associated with telomere length: a collaborative cross-sectional meta-analysis of 87 observational studies. Am J Clin Nutr. 2018;108:453–75.

39. Smith EM, Pendlebury DF, Nandakumar J. Structural biology of telomeres and telomerase. Cell Mol Life Sci. 2020;77:61–79.

40. Aramburu T, Plucinsky S, Skordalakes E. POT1-TPP1 telomere length regulation and disease. Comput Struct Biotechnol J. 2020;18:1939–46.

41. Wang F, Podell ER, Zaug AJ, Yang Y, Baciu P, Cech TR, et al. The POT1-TPP1 telomere complex is a telomerase processivity factor. Nature. 2007;445:506–10.

42. He H, Wang Y, Guo X, Ramchandani S, Ma J, Shen M-F, et al. Pot1b deletion and telomerase haploinsufficiency in mice initiate an ATR-dependent DNA damage response and elicit phenotypes resembling dyskeratosis congenita. Mol Cell Biol. 2009;29:229–40.

43. Wu L, Multani AS, He H, Cosme-Blanco W, Deng Y, Deng JM, et al. Pot1 deficiency initiates DNA damage checkpoint activation and aberrant homologous recombination at telomeres. Cell. 2006;126:49–62.

44. Wang H, Ni J, Guo X, Zhou T, Ma X, Xue J, et al. Shelterin differentially respond to oxidative stress induced by TiO2-NPs and regulate telomere length in human hepatocytes and hepatocarcinoma cells in vitro. Biochem Biophys Res Commun. 2018;503:697–702.

45. Chuenwisad K, More-Krong P, Tubsaeng P, Chotechuang N, Srisa-Art M, Storer RJ, et al. Premature Senescence and Telomere Shortening Induced by Oxidative Stress From Oxalate, Calcium Oxalate Monohydrate, and Urine From Patients With Calcium Oxalate Nephrolithiasis. Front Immunol. 2021;12:696486.

46. Glousker G, Briod A-S, Quadroni M, Lingner J. Human shelterin protein POT1 prevents severe telomere instability induced by homology-directed DNA repair. EMBO J. 2020;39:e104500.

47. Wang T, Mei S, Fu R, Wang H, Shao Z. Expression of shelterin component POT1 is associated with decreased telomere length and immunity condition in humans with severe aplastic anemia. J Immunol Res. 2014;2014:439530.

48. Shalev I, Entringer S, Wadhwa PD, Wolkowitz OM, Puterman E, Lin J, et al. Stress and telomere biology: a lifespan perspective. Psychoneuroendocrinology. 2013;38:1835–42.

49. Patel J, Baptiste BA, Kim E, Hussain M, Croteau DL, Bohr VA. DNA damage and mitochondria in cancer and aging. Carcinogenesis. 2020;41:1625–34.

ACKNOWLEDGEMENTS

The authors thank the FondaMental Foundation (www-fondation-fondamental.org) for its support for the infrastructure of Bipolar Expert Centers. We thank the clinicians and nursing staff of the Center for Expertise in BD for the recruitment and clinical characterization of individuals with BD. The biobank is hosted at the Plateforme de Ressources Biologique, Hôpital Henri Mondor, Créteil (Dr C. Barau). Dr. C. Marie-Claire is supported by the Centre National de la Recherche Scientifique.
## AUTHOR CONTRIBUTIONS

BE and CMC designed the study. LS and CC developed the experiments. LS and AL performed the experiments. LS performed data analysis and interpretation with inputs from BE and CMC. LS drafted and revised the manuscript. BE and CMC reviewed and edited the manuscript. All authors contributed to the supervision of the project and revised the manuscript. All authors approved the submitted version of the manuscript.

## FUNDING

This work was supported (in part) by the Fondation FondaMental, Créteil, France and by the Investissements d'Avenir programs managed by the ANR under references ANR-11-IDEX-0004-02 and ANR-10-COH0-10-01. This project was also supported by a grant from Fondation FondaMental (Prix Marcel Dassault 2020 to B. Etain), a grant from Fondation de France (2018 to B. Etain). L. Spano is supported by a grant from Fondation FondaMental (Prix FACE 2021).

## 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/s41398-024-02824-z.

**Correspondence** and requests for materials should be addressed to Cynthia Marie-Claire.

**Reprints and permission information** is available at http://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

## FONDAMENTAL ADVANCED CENTERS OF EXPERTISE IN BIPOLAR DISORDERS (FACE-BD) COLLABORATORS

B. Etain<fnref n="21" />, E. Olié<fnref n="21" />, M. Leboyer<fnref n="21" />, E. Haffen<fnref n="21" />, P. M. Llorca<fnref n="21" />, V. Barteau²², S. Bensalem²², O. Godin²², H. Laouamri²², K. Souryis²², S. Hotier²³,²⁴, A. Pelletier²³,²⁴, N. Drancourt²³,²⁴, J. P. Sanchez²³,²⁴, E. Saliou²³,²⁴, C. Hebbache²³,²⁴, J. Petrucci²³,²⁴, L. Willaume²³,²⁴, E. Bourdin²³,²⁴, F. Bellivier²⁵, B. Etain²⁵, V. Hennion²⁵, E. Marlinge²⁵, P. Lebard²⁵, B. Antoniol²⁶, A. Desage²⁶, S. Gard²⁶, A. Jutant²⁶, K. Mbailara²⁶, I. Minois²⁶, L. Zanouy²⁶, C. Abettan²⁷, L. Bardin²⁷, A. Cazals²⁷, P. Courtet²⁷, B. Deffinis²⁷, D. Ducasse²⁷, M. Gachet²⁷, A. Henrion²⁷, E. Martinerie²⁷, F. Molière²⁷, B. Noisette²⁷, E. Olié²⁷, G. Tarquini²⁷, J. M. Azorin²⁸, R. Belzeaux²⁸, N. Correard²⁸, J. L. Consoloni²⁸, F. Groppi²⁸, L. Lescalier²⁸, J. Montant²⁸, M. Rebattu²⁸, N. Viglianese²⁸, R. Cohen²⁹, J. P. Kahn²⁹, M. Milazzo²⁹, O. Wajsbrot-Elgrabli²⁹, T. Bougerol³⁰, B. Fredembach³⁰, Q. Denoual³⁰, A. Bertrand³⁰, A. Pouchon³⁰, M. Polosan³⁰, L. Brehon³¹, G. Bony³¹, L. Durand³¹, V. Feuga³¹, N. Kayser³¹, C. Passerieux³¹, P. Roux³¹, V. Aubin³², I. Cussac³², M. A. Dupont³², J. Loftus³², I. Medecin³², C. Dubertret³³, N. Mazer³³, C. Portalier³³, C. Scognamiglio³³, A. Bing³³, P. Laurent³³, C. Beal³⁴, O. Blanc³⁴, T. Bonnet³⁴, D. Lacelle³⁴, P. M. Llorca³⁴, M. Mennetrier³⁴, L. Samalin³⁴ and M. Vayssié³⁴

<footnote n="21">²¹FACE-BD Clinical Coordinating Center (Fondation FondaMental), Créteil, France. ²²FACE-BD Data Coordinating Center (Fondation FondaMental), Créteil, France. ²³FACE-BD Clinical Sites and Principal Collaborators in France, Montpellier, France. ²⁴AP-HP, DHU PePSY, Pôle de Psychiatrie et d'Addictologie des Hôpitaux Universitaires H Mondor, Créteil, France. ²⁵AP-HP, GH Saint-Louis-Lariboisière-Fernand Widal, Pôle Neurosciences, Paris, France. ²⁶Hôpital C. Perrens, Centre Expert Trouble Bipolaire, Service de Psychiatrie Adulte, Pôle 3-4-7, Bordeaux, France. ²⁷Département d'Urgence et Post Urgence Psychiatrique, CHRU Montpellier, Montpellier, France. ²⁸Pôle de Psychiatrie, addictologie et pédopsychiatrie, Hôpital Sainte Marguerite, Marseille, France. ²⁹Service de Psychiatrie et Psychologie Clinique, CHU de Nancy, Hôpitaux de Brabois, Vandoeuvre Les Nancy, France. ³⁰Service Universitaire de Psychiatrie, CHU de Grenoble et des Alpes, Grenoble, France. ³¹Centre Hospitalier de Versailles, Service Universitaire de Psychiatrie d'adultes, Le Chesnay, France. ³²Centre Hospitalier Princesse Grace, Monaco, Monaco. ³³APHP, Groupe Hospitalo-universitaire AP-HP Nord, DMU ESPRIT, Service de Psychiatrie et Addictologie, Hôpital Louis Mourier, Colombes, France. ³⁴Service de Psychiatrie d'adultes B, Centre Expert Trouble Bipolaire, CHU de Clermont-Ferrand, Clermont-Ferrand, France.</footnote>
## AUTHOR CONTRIBUTIONS

BE and CMC designed the study. LS and CC developed the experiments. LS and AL performed the experiments. LS performed data analysis and interpretation with inputs from BE and CMC. LS drafted and revised the manuscript. BE and CMC reviewed and edited the manuscript. All authors contributed to the supervision of the project and revised the manuscript. All authors approved the submitted version of the manuscript.

## FUNDING

This work was supported (in part) by the Fondation FondaMental, Créteil, France and by the Investissements d'Avenir programs managed by the ANR under references ANR-11-IDEX-0004-02 and ANR-10-COH0-10-01. This project was also supported by a grant from Fondation FondaMental (Prix Marcel Dassault 2020 to B. Etain), a grant from Fondation de France (2018 to B. Etain). L. Spano is supported by a grant from Fondation FondaMental (Prix FACE 2021).

## 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/s41398-024-02824-z.

**Correspondence** and requests for materials should be addressed to Cynthia Marie-Claire.

**Reprints and permission information** is available at http://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_5.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

## FONDAMENTAL ADVANCED CENTERS OF EXPERTISE IN BIPOLAR DISORDERS (FACE-BD) COLLABORATORS

B. Etain<fnref n="21" />, E. Olié<fnref n="21" />, M. Leboyer<fnref n="21" />, E. Haffen<fnref n="21" />, P. M. Llorca<fnref n="21" />, V. Barteau²², S. Bensalem²², O. Godin²², H. Laouamri²², K. Souryis²², S. Hotier²³,²⁴, A. Pelletier²³,²⁴, N. Drancourt²³,²⁴, J. P. Sanchez²³,²⁴, E. Saliou²³,²⁴, C. Hebbache²³,²⁴, J. Petrucci²³,²⁴, L. Willaume²³,²⁴, E. Bourdin²³,²⁴, F. Bellivier²⁵, B. Etain²⁵, V. Hennion²⁵, E. Marlinge²⁵, P. Lebard²⁵, B. Antoniol²⁶, A. Desage²⁶, S. Gard²⁶, A. Jutant²⁶, K. Mbailara²⁶, I. Minois²⁶, L. Zanouy²⁶, C. Abettan²⁷, L. Bardin²⁷, A. Cazals²⁷, P. Courtet²⁷, B. Deffinis²⁷, D. Ducasse²⁷, M. Gachet²⁷, A. Henrion²⁷, E. Martinerie²⁷, F. Molière²⁷, B. Noisette²⁷, E. Olié²⁷, G. Tarquini²⁷, J. M. Azorin²⁸, R. Belzeaux²⁸, N. Correard²⁸, J. L. Consoloni²⁸, F. Groppi²⁸, L. Lescalier²⁸, J. Montant²⁸, M. Rebattu²⁸, N. Viglianese²⁸, R. Cohen²⁹, J. P. Kahn²⁹, M. Milazzo²⁹, O. Wajsbrot-Elgrabli²⁹, T. Bougerol³⁰, B. Fredembach³⁰, Q. Denoual³⁰, A. Bertrand³⁰, A. Pouchon³⁰, M. Polosan³⁰, L. Brehon³¹, G. Bony³¹, L. Durand³¹, V. Feuga³¹, N. Kayser³¹, C. Passerieux³¹, P. Roux³¹, V. Aubin³², I. Cussac³², M. A. Dupont³², J. Loftus³², I. Medecin³², C. Dubertret³³, N. Mazer³³, C. Portalier³³, C. Scognamiglio³³, A. Bing³³, P. Laurent³³, C. Beal³⁴, O. Blanc³⁴, T. Bonnet³⁴, D. Lacelle³⁴, P. M. Llorca³⁴, M. Mennetrier³⁴, L. Samalin³⁴ and M. Vayssié³⁴

<footnote n="21">²¹FACE-BD Clinical Coordinating Center (Fondation FondaMental), Créteil, France. ²²FACE-BD Data Coordinating Center (Fondation FondaMental), Créteil, France. ²³FACE-BD Clinical Sites and Principal Collaborators in France, Montpellier, France. ²⁴AP-HP, DHU PePSY, Pôle de Psychiatrie et d'Addictologie des Hôpitaux Universitaires H Mondor, Créteil, France. ²⁵AP-HP, GH Saint-Louis-Lariboisière-Fernand Widal, Pôle Neurosciences, Paris, France. ²⁶Hôpital C. Perrens, Centre Expert Trouble Bipolaire, Service de Psychiatrie Adulte, Pôle 3-4-7, Bordeaux, France. ²⁷Département d'Urgence et Post Urgence Psychiatrique, CHRU Montpellier, Montpellier, France. ²⁸Pôle de Psychiatrie, addictologie et pédopsychiatrie, Hôpital Sainte Marguerite, Marseille, France. ²⁹Service de Psychiatrie et Psychologie Clinique, CHU de Nancy, Hôpitaux de Brabois, Vandoeuvre Les Nancy, France. ³⁰Service Universitaire de Psychiatrie, CHU de Grenoble et des Alpes, Grenoble, France. ³¹Centre Hospitalier de Versailles, Service Universitaire de Psychiatrie d'adultes, Le Chesnay, France. ³²Centre Hospitalier Princesse Grace, Monaco, Monaco. ³³APHP, Groupe Hospitalo-universitaire AP-HP Nord, DMU ESPRIT, Service de Psychiatrie et Addictologie, Hôpital Louis Mourier, Colombes, France. ³⁴Service de Psychiatrie d'adultes B, Centre Expert Trouble Bipolaire, CHU de Clermont-Ferrand, Clermont-Ferrand, France.</footnote>