Introduction
Over 100 million people have been infected with severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), and over 1 million have died from coronavirus disease 2019 (COVID-19) worldwide (COVID-19 Map, 2021). Several nationwide studies have suggested that individuals with severe mental disorders may have worse clinical outcomes from COVID-19 (Egede et al., Reference Egede, Campbell, Walker, Garacci, Dawson and Egede2021; Fond et al., Reference Fond, Pauly, Leone, Llorca, Orleans, Loundou and Boyer2020a, Reference Fond, Pauly, Orleans, Antonini, Fabre, Sanz and Boyer2020b; Lee et al., Reference Lee, Yang, Moon, Yoo, Ha, Kim and Yon2020; Li, Li, Fortunati, & Krystal, Reference Li, Li, Fortunati and Krystal2020; Nemani et al., Reference Nemani, Li, Olfson, Blessing, Razavian, Chen and Goff2021; Poblador-Plou et al., Reference Poblador-Plou, Carmona-Pírez, Ioakeim-Skoufa, Poncel-Falcó, Bliek-Bueno and Cano-Del Pozo2020; Tzur Bitan et al., Reference Tzur Bitan, Krieger, Kridin, Komantscher, Scheinman, Weinstein and Feingold2021; Wang, Xu, & Volkow, Reference Wang, Xu and Volkow2020). The US and Korean studies reported a global increased risk of COVID-19 mortality in patients with various mental disorders [e.g. bipolar disorders (BDs), schizophrenia, depression, attention deficit hyperactivity disorder, eating disorders, and personality disorders] (Lee et al., Reference Lee, Yang, Moon, Yoo, Ha, Kim and Yon2020; Wang et al., Reference Wang, Xu and Volkow2020). The French study targeted COVID-19 patients with schizophrenia and reported that increased mortality was not homogeneous according to age or other clinical characteristics, suggesting that COVID-19 health care strategies should be personalized according to vulnerability levels (Fond et al., Reference Fond, Pauly, Leone, Llorca, Orleans, Loundou and Boyer2020a, Reference Fond, Pauly, Orleans, Antonini, Fabre, Sanz and Boyer2020b). A better understanding of the vulnerability to severe COVID-19 in mental disorders other than schizophrenia is an urgent research priority for proposing personalized interventions for COVID-19. To date, no study has targeted COVID-19 patients with BDs, even though BD is one of the most frequent and severe mental disorders, affecting more than 1% of the world's population (McIntyre et al., Reference McIntyre, Berk, Brietzke, Goldstein, López-Jaramillo, Kessing and Mansur2020). BD is associated with a 10–20-year loss of life expectancy, mostly due to cardiovascular diseases and multiple metabolic and addiction comorbidities aggravated by poor compliance with medical follow-up (McIntyre et al., Reference McIntyre, Berk, Brietzke, Goldstein, López-Jaramillo, Kessing and Mansur2020). BD is also associated with socioeconomic deprivation (Schoeyen et al., Reference Schoeyen, Birkenaes, Vaaler, Auestad, Malt, Andreassen and Morken2011). All these factors were reported to be associated with an increased risk of severe COVID-19 (Myers, Parodi, Escobar, & Liu, Reference Myers, Parodi, Escobar and Liu2020; Schoeyen et al., Reference Schoeyen, Birkenaes, Vaaler, Auestad, Malt, Andreassen and Morken2011; Yancy, Reference Yancy2020; Zhou et al., Reference Zhou, Yu, Du, Fan, Liu, Liu and Cao2020), suggesting that patients with BD may be at an increased risk of COVID-19-related death. In contrast with this hypothesis, a recent large US study suggested that diagnoses of mood disorders were not associated with increased mortality (Nemani et al., Reference Nemani, Li, Olfson, Blessing, Razavian, Chen and Goff2021). It remains unknown whether COVID-19 patients with BD are at an increased risk of mortality and should be targeted as a high-risk population for severe forms of COVID-19. More information is thus needed about the clinical characteristics, care and outcomes of COVID-19 patients with BD.
In this nationwide study, we aimed to determine whether clinical characteristics, outcomes, and care differed between patients with BD and patients without a diagnosis of severe mental illness. The primary objective was to compare in-hospital mortality between patients with BD and patients without a diagnosis of severe mental illness after adjustment for the main clinical risk factors for COVID-19. The secondary objective was to compare intensive care unit (ICU) admissions between patients with BD and patients without a diagnosis of severe mental illness.
Methods
Study design and data sources
In this population-based cohort study, we used data from the Programme de Médicalisation des Systèmes d'Information (PMSI database), the French national hospital database in which administrative and medical data are systematically collected for acute and psychiatric care. The PMSI database is based on diagnosis-related groups, with all diagnoses coded according to the 10th revision of the International Classification of Diseases (ICD-10) and procedural codes from the Classification Commune des Actes Médicaux. In our study, we included all hospitalized patients in acute care settings between 1 February 2020 and 9 June 2020, aged 15 years or older, with identified COVID-19 (ICD-10 codes: U07.10 or U07.12 or U07.14) and respiratory symptoms (ICD-10 = U07.10 or U07.11), and a length of hospital stay >24 h (to exclude COVID-19 forms that did not require hospitalization), except if the patients died within 24 h. The threshold of 15 years was chosen because it is the usual threshold to distinguish adult psychiatry from pedopsychiatry in France. We excluded patients with a severe mental illness diagnosis other than BD, such as schizophrenia or recurrent major depression (ICD-10 codes = F20*, F22*, F25*, or F33*).
The PMSI database is used to determine financial resources and is frequently and thoroughly verified by both its producer and the paying party, with possible financial and legal consequences for inaccuracies (Boudemaghe & Belhadj, Reference Boudemaghe and Belhadj2017). Data from the PMSI database are deidentified and can be reused for research purposes (Fond et al., Reference Fond, Salas, Pauly, Baumstarck, Bernard, Orleans and Boyer2019b; Revon-Rivière et al., Reference Revon-Rivière, Pauly, Baumstarck, Bernard, André, Gentet and Boyer2019). Due to its suitable accuracy and exhaustive data collection, no patients were lost to follow-up during the study period.
Ethical considerations
Because this study was strictly observational and based on anonymous data, written informed consent from the participants or the authorization from an ethical committee for dealing with human issues was not required in accordance with the French laws.
Procedures
We defined four groups following two steps.
Step 1. Cases were patients who had a diagnosis of BD according to specific ICD-10 codes [i.e. F30* (manic episode) or F31* (bipolar affective disorders)] in either the acute care or psychiatric PMSI database. Controls were patients who did not have a diagnosis of severe mental illness according to specific ICD-10 codes in the acute care PMSI database and who were not listed in the PMSI psychiatry database.
Step 2. To obtain comparability between the groups defined in step 1, we conducted propensity score (PS) matching. Two matched cohorts of patients were then formed: matched cases and matched controls.
Outcomes and collected data
The primary outcome was in-hospital mortality during the index stay for COVID-19. The secondary outcome was ICU admission. We gathered patient sociodemographic data [age classes: <55, 55–65, 65–80 and >80 years; sex: male, female; social deprivation: socially favored, socially deprived (Rey, Jougla, Fouillet, & Hémon, Reference Rey, Jougla, Fouillet and Hémon2009)], clinical data at baseline using ICD-10 codes [smoking addiction: yes, no; weight: normal weight, overweight and obesity; Charlson comorbidity index score categorized as follows: 0, 1–2, 3 and more (Bannay et al., Reference Bannay, Chaignot, Blotière, Basson, Weill, Ricordeau and Alla2016) and main comorbidities: yes, no], stay data (origin of patients: from home or from hospital institution; length of ICU and hospital stay in days), and hospital data (hospital category: public, university and private; number of hospital stays for COVID-19). Social deprivation involves four socioeconomic ecological variables: the proportion of residents who graduated from high school, median household income, the percentage of residents who are blue-collar workers, and the unemployment rate (Rey et al., Reference Rey, Jougla, Fouillet and Hémon2009). It has been categorized into two categories favored v. deprived using zero as the threshold.
Statistical analysis
The baseline characteristics of our study subjects are represented as the counts (proportions) and medians (interquartile ranges) for categorical and continuous variables, respectively.
Analyses were then carried out in two steps. First (step 1), we used either the χ2 test or Student's t test to compare sociodemographic, clinical, and hospital data between cases and unmatched controls. Second (step 2), PSs (Austin, Reference Austin2011; Marrie, Dawson, & Garland, Reference Marrie, Dawson and Garland2009) were used for matching, considering the following characteristics (p values <0.2 from step 1): age classes (<55, 55–65, 65–75, and ⩾80 years), sex (male, female), social deprivation (favored, deprived), smoking addiction (yes, no), weight (normal weight, overweight, and obesity), and origin of patients (from home or from hospital institution). In addition, all Charlson comorbidities (yes, no) were included regardless of the p-value. We matched patients with the closest PS inside the hospital to control for confounders at the hospital level using a greedy 1:20 max algorithm without replacement and requiring that the logit of the PS of a patient with BD and one without a diagnosis of severe mental illness be within 0.20 standard deviations of one another. Age, sex, and dementia were exactly matched due to persistent heterogeneity between groups, even when these variables were involved in the PS score. Standardized differences were used to assess the balance between the matched groups for baseline characteristics using weights to normalize the distribution of patients with BD to equal the number of patients without a diagnosis of severe mental illness. An absolute standardized difference (SD) of ⩽0.20 was chosen to indicate a negligible difference in the mean or prevalence of a variable between groups. To study the association between BD and the two outcomes, the odds ratio (OR) and 95% confidence interval (95% CI) were estimated using a multivariable generalized linear model (logit function), with a robust variance estimator to account for clustering within matched pairs. Several models were performed for each outcome: model 1 without adjustment; model 2: adding to model 1 the variables with the most important differences after PS matching (SD > 0.1); model 3 (only for mortality): adding ICU admission to model 1; and model 4: adding to model 2 three interaction terms (BD × age, BD × sex, and BD × social deprivation) to check whether the association between BD and the two outcomes was homogenous according to age, sex, and social deprivation. This interaction was determined based on a previous study reporting the influence of age on the COVID-19 prognosis in schizophrenia (Fond et al., Reference Fond, Pauly, Leone, Llorca, Orleans, Loundou and Boyer2020a, Reference Fond, Pauly, Orleans, Antonini, Fabre, Sanz and Boyer2020b; Lakbar, Luque-Paz, Mege, Einav, & Leone, Reference Lakbar, Luque-Paz, Mege, Einav and Leone2020; Lewis et al., Reference Lewis, Friedrichs, Wagstaff, Sage, LaCross, Bui and Dunn2020).
To test the robustness of our results, three series of sensitivity analyses were performed. We repeated the analyses for the two outcomes: (i) on the two subgroups of patients: F30* (manic episode) or F31* (bipolar affective disorders); (ii) by completing the matching approach described above with exact matching for tobacco smoking and a caliper = 0.10; and (iii) by matching on all the variables with SD > 0.05 using a narrower caliper = 0.05.
A significance threshold of p < 0.05 was used. All analyses were performed using SAS (version 9.4).
Results
Overall, 95 662 patients with COVID-19 were hospitalized between 1 February 2020 and 9 June 2020. After excluding patients aged <15 years, those without identified virus and respiratory symptoms, those with a length of stay <24 h (except if they died at the hospital) and those with a diagnosis of schizophrenia or recurrent major depressive disorder, 50 407 patients were included in the current study (flow chart, Fig. 1). The overall in-hospital mortality rate was 21.8%.
In the first step, 480 patients with BD (449 with F31 code and 31 with F30 code) were compared to 49 927 controls without a diagnosis of severe mental illness (Table 1). Compared to controls, patients with BD had an increased in-hospital mortality rate (26.0% v. 21.7%, p = 0.0341), with no difference in the ICU admission rate (27.9% v. 28.4%, p = 0.7985). In addition, patients with BD were slightly older [mean 70.92 standard deviation (13.48) v. 68.72 (17.07), p < 0.0001], more frequently women (62.1% v. 43.1%, p < 0.0001), more socially favored (57.9% v. 48.6%, p = 0.0037), particularly in patients with BD older than 55 years, more frequently diagnosed with smoking addiction (9.0% v. 4.2%, p = 0.0003), higher Charlson comorbidity score (p = 0.0002), more frequent chronic renal disease (16.5% v. 12.1%, p = 0.0106), dementia (21.3% v. 9.6%, p < 0.0001), cerebrovascular disease (9.8% v. 6.3%, p = 0.0084), and chronic obstructive pulmonary disease (18.5% v. 12.6%, p = 0.0021). They were less frequently admitted from home (84.8% v. 90.3%, p = 0.0020) and had a shorter length of hospital stay [median 9-day interquartile range (6–19) v. 12 days (5–16)]. The distribution of hospitalizations by month was similar between patients with BD and controls (online Supplementary Fig. S1).
N: effective; %: percentage; IQR: interquartile range; BD: bipolar disorder.
p-value in bold: statistical significance.
a Body mass index: <25: normal weight; 25–30: overweight; >30: obesity.
In the second step, 452 patients with BD were matched with 6171 controls without a diagnosis of severe mental illness (Table 2). The patients with BD were found to have an increased in-hospital mortality rate compared to their matched controls (26.6% v. 21.4%) regardless of the tested models, while no significant difference was observed for ICU admission (27.7% v. 26.2%) (Figs 2 and 3; Table 3). Significant interactions between BD and age and between BD and social deprivation for in-hospital mortality were found, highlighting that the most important inequalities in mortality were observed in the youngest (OR, 95% CI 2.28 (1.18–4.41), p = 0.0015) and most deprived patients with BD (OR, 95% CI 1.60 (1.33–1.92), p < 0.001). The youngest patients with BD had more direct ICU admission with a more severe state of health than controls. They were also more likely to require invasive mechanical ventilation (online Supplementary File 1).
N, effective; %, percentage; IQR, interquartile range; BD, bipolar disorder; SD, standardized difference.
p-Value in bold: statistical significance.
SD > 0.2 indicated an important imbalance, and SD > 0.1 was included in the multivariate models.
a Body mass index: <25: normal weight; 25–30: overweight; >30: obesity.
ICU, intensive care unit.
Model 1: unadjusted model.
Model 2: model 1 plus adjustment for variables with SD > 0.1.
Model 3: model 1 plus adjustment for ICU admission.
Model 4: model 2 plus the interaction terms. Reference: matched controls.
OR (95% CI): odds ratio (95% confidence interval).
p-Value in bold: statistical significance.
In the sensitivity analyses, the results were maintained when excluding the 31 patients with F30* code (manic episode) [adjusted OR for mortality: 1.30 (1.16–1.45), p < 0.0001, p > 0.05 for ICU admission, analyses not shown], with the exact matching for tobacco smoking and caliper = 0.10 [adjusted OR for mortality: 1.30 (1.13–1.45), p < 0.0001, p > 0.05 for ICU admission, see online Supplementary File 2] and with the matching on all the variables with SD > 0.05 using caliper = 0.05 [adjusted OR for mortality: 1.27 (1.04–1.54), p = 0.017, p > 0.05 for ICU admission, analyses not shown].
Discussion
To our knowledge, this is the largest series of patients with BD with COVID-19 reported to date, including 480 patients with BD hospitalized for COVID-19 infection with a nationwide geographical distribution. COVID-19 patients with BD were at an increased risk of mortality after adjustment for clinical risk factors for COVID-19, which was exacerbated in the youngest and most deprived patients with BD. In contrast to patients with schizophrenia, no difference was reported for ICU admission (Fond et al., Reference Fond, Pauly, Leone, Llorca, Orleans, Loundou and Boyer2020a, Reference Fond, Pauly, Orleans, Antonini, Fabre, Sanz and Boyer2020b).
Our prematching findings inform us about BD patient characteristics, highlighting several risk factors for COVID-19 aggravation (i.e. obesity, renal disease, dementia and cerebrovascular disease, and chronic obstructive pulmonary disease). Obesity has been reported as a risk factor for COVID-19 hospital admission and worse outcomes (Lighter et al., Reference Lighter, Phillips, Hochman, Sterling, Johnson, Francois and Stachel2020). In a nationwide study, hospitalized COVID-19 patients with chronic renal disease had significantly higher mortality than those without renal disease (Ozturk et al., Reference Ozturk, Turgutalp, Arici, Odabas, Altiparmak, Aydin and Ates2020). A preexisting diagnosis of dementia was found to be an independent risk factor for COVID-19 hospitalization and mortality (Atkins et al., Reference Atkins, Masoli, Delgado, Pilling, Kuo, Kuchel and Melzer2020). Chronic obstructive pulmonary disease has also been reported to worsen the progression and prognosis of COVID-19 (Zhao et al., Reference Zhao, Meng, Kumar, Wu, Huang, Lian and Lin2020). Beyond the strict question of COVID-19, these data are a strong reminder of the need to pursue efforts to prevent and manage somatic diseases in BD (Heiberg et al., Reference Heiberg, Jacobsen, Balteskard, Bramness, Naess, Ystrom and Høye2019; Nielsen, Kugathasan, Straszek, Jensen, & Licht, Reference Nielsen, Kugathasan, Straszek, Jensen and Licht2019). In contrast, one BD patient characteristic was less expected. Patients with BD were largely more socially favored than controls, although BD is known to induce impaired professional and social functioning and poor socioeconomic outcomes (Hakulinen, Elovainio, et al., Reference Hakulinen, Elovainio, Arffman, Lumme, Pirkola, Keskimäki and Böckerman2019; Hakulinen, Musliner, & Agerbo, Reference Hakulinen, Musliner and Agerbo2019). Social deprivation is an important cause of premature mortality among patients with BD (Fond et al., Reference Fond, Pauly, Bege, Orleans, Braunstein, Leone and Boyer2019a; Martin et al., Reference Martin, McLean, Park, Martin, Connolly, Mercer and Smith2014), which can thus explain a broader representation of socially favored patients with BD, particularly in the oldest patients with BD. The overrepresentation of socially favored patients with BD may also suggest lower access to hospital care in deprived patients with BD than in controls. Previous studies have suggested that, among deprived individuals, patients with severe mental illnesses, such as BD, are some of the most underserved and undertreated populations (Nielsen et al., Reference Nielsen, Kugathasan, Straszek, Jensen and Licht2019; Tinland et al., Reference Tinland, Loubière, Boucekine, Boyer, Fond, Girard and Auquier2020). In a recent study on homeless individuals, the main negative predictors for hospital use were a lower social functioning score and BDs (Loubière et al., Reference Loubière, Tinland, Taylor, Loundou, Girard, Boyer and Auquier2020). Further studies are warranted to explore disease management during the period between the onset of symptoms and hospitalization in patients with BD.
In contrast with a recent US study reporting that mood disorders were not associated with increased mortality (Nemani et al., Reference Nemani, Li, Olfson, Blessing, Razavian, Chen and Goff2021), we found that COVID-19 patients with BD were at an increased risk of mortality after adjustment for the main clinical risk factors for COVID-19 (i.e. age, obesity, renal disease, dementia and cerebrovascular disease, and chronic obstructive pulmonary disease). This discrepancy suggests that mortality is not homogeneous within affective disorders and that patients with BD should be targeted as a high-risk population for severe forms of COVID-19. Several factors can be advanced to explain this health inequity between BD and the general population. Abnormal mood may alter decision-making or treatment adherence, which may impact prognosis. A delay in obtaining medical attention and in access to hospital care may explain the higher severity of patients with BD, as suggested in the youngest patients with BD who had more frequent direct ICU admissions. The barriers to access to somatic care for patients with BD have been described in previous studies (Hernández-Gómez, Andrade-González, Lahera, & Vieta, Reference Hernández-Gómez, Andrade-González, Lahera and Vieta2020). Immunological disturbances have also been suggested to play an important role in the mortality of young patients with BD without classic comorbidities (Lu et al., Reference Lu, Huang, Wang, Tang, Wang, Gao and Zhao2021). Human leukocyte antigen predominantly regulates viral infection, especially COVID-19 (Tamouza, Krishnamoorthy, & Leboyer, Reference Tamouza, Krishnamoorthy and Leboyer2021), and has recently been suggested to play a role in the etiology of BD (Tamouza et al., Reference Tamouza, Oliveira, Etain, Bengoufa, Hamdani, Manier and Leboyer2018). Patients with BD may also be at higher risk of hypovitaminosis D, contributing to poor COVID-19 prognosis (Di Nicola et al., Reference Di Nicola, Dattoli, Moccia, Pepe, Janiri, Fiorillo and Sani2020; Hastie et al., Reference Hastie, Mackay, Ho, Celis-Morales, Katikireddi, Niedzwiedz and Pell2020). Psychotropic drugs are also modulators of the immune system. Lithium, widely used in treating BD, also exhibits antiviral activity that has not been shown to protect patients with BD (Ishii, Terao, & Hirakawa, Reference Ishii, Terao and Hirakawa2021). Antipsychotics are used in both schizophrenia and BD treatment and modify immune dysfunction (May, Slitzky, Rostama, Barlow, & Houseknecht, Reference May, Slitzky, Rostama, Barlow and Houseknecht2020), which may explain the shared increased mortality between BD and schizophrenia patients (26% in both populations) (Fond et al., Reference Fond, Pauly, Leone, Llorca, Orleans, Loundou and Boyer2020a, Reference Fond, Pauly, Orleans, Antonini, Fabre, Sanz and Boyer2020b). The background regimen of patients with BD may also have drug interactions with anti-COVID-19 therapies, such as antimicrobials or immunomodulators. Patients with BD are often treated with long-term benzodiazepines. They may therefore require higher doses of sedatives due to tolerance issues, which may lead to increased complications. Thus, in the ICU, these patients may require muscle relaxant agents more frequently, which may lead to increased complications upon mechanical ventilation withdrawal. Further studies are warranted to better understand and to address modifiable factors that can prevent this increased mortality in COVID-19 patients with BD. Our study also underscores the exacerbation of health disparities in the youngest and most deprived patients with BD. Targeted interventions are needed for the youngest and most deprived patients with BD to reduce these stark disparities.
Strengths
The current study is the first national-based study of COVID-19 mortality in patients with BD. Many confounding factors were captured in the PS matching, which ensured the robustness of our results. The PMSI has a system of coding at discharge with strict variable definitions (primary and secondary diagnoses using the ICD-10; as well as procedure codes associated with the care provided) and a subset of records audited on a regular basis to avoid excessively high rates of coding errors.
Limits
Only the in-hospital mortality rate was explored in the current study. Ambulatory variables (treatment observance, medical follow-up, and social isolation) and ethnicity were not captured in the database and have therefore not been explored and may have impacted the COVID-19 prognosis. Some data are known to be insufficiently coded in medico-administrative databases (e.g. tobacco smoking and obesity) (Fond et al., Reference Fond, Pauly, Leone, Llorca, Orleans, Loundou and Boyer2020a, Reference Fond, Pauly, Orleans, Antonini, Fabre, Sanz and Boyer2020b) and lifestyle factors were not captured. The mood state of the patients was not available in the PMSI database, and it is not known whether an acute manic or depressed episode may have impacted the COVID-19 prognosis of patients with BD. As the precise cause of death was not captured in the database, we cannot definitely conclude that COVID-19 was the cause of death. However, the patients died in median during the 11 days following hospital admission for COVID-19 infection, we can therefore reasonably assume that COVID-19 infection played an important role in their death. As mentioned in our rationale, our purpose was to explore specifically BDs, the most severe adult psychiatric disorder with schizophrenia that was already explored in our previous paper (Fond et al., Reference Fond, Pauly, Leone, Llorca, Orleans, Loundou and Boyer2020a, Reference Fond, Pauly, Orleans, Antonini, Fabre, Sanz and Boyer2020b). These diagnoses are the most well identified and coded psychiatric diagnoses in the PMSI database, contrary to other mood or personality disorders, especially if patients have not been hospitalized or are followed-up in the private sector. Future studies should specifically explore mood and personality disorders considering these limitations.
Conclusion
COVID-19 patients with BD were at an increased risk of mortality after adjustment for the main clinical risk factors for COVID-19 (i.e. age, obesity, renal disease, dementia and cerebrovascular disease, and chronic obstructive pulmonary disease), suggesting the existence of other factors that lead to this health inequity (i.e. delay in access to care, immunological disturbances, and the effects of psychotropic drugs). This increased mortality was exacerbated in the youngest and most deprived patients with BD. COVID-19 patients with BD should thus be targeted as a high-risk population for severe forms of COVID-19, requiring enhanced preventive care and disease management strategies.
Supplementary material
The supplementary material for this article can be found at https://doi.org/10.1017/S0033291721003676
Acknowledgements
The authors thank David Braunstein, Myriam Dubuc, Cyprien Fabre, Marie-Thérèse Jimeno, Vincent Pradel, Anne Remacle, Fanny Romain, Catherine Seyler, Françoise Volot, and all the other members of the Department of Medical Information.
Author contributions
Veronica Orleans and Vanessa Pauly had full access to all of the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis. Concept and design: Laurent Boyer and Guillaume Fond. Acquisition, analysis, or interpretation of data: Pascal Auquier, Karine Baumstarck, Laurent Boyer, Christophe Lancon, Pierre-Michel Llorca, Guillaume Fond, Marc Leone, Vanessa Pauly, and Veronica Orleans. Drafting of the manuscript: Laurent Boyer and Guillaume Fond. Critical revision of the manuscript for important intellectual content: All the authors. Statistical analysis: Vanessa Pauly. Administrative, technical, or material support: Veronica Orleans. Supervision: Laurent Boyer.
Conflict of interest
The authors declare no competing interests.