key: cord-0021159-ttbai8h2 authors: Houist, Anne-Laure; Bondeelle, Louise; Salmona, Maud; LeGoff, Jérôme; de Latour, Régis Peffault; Rivière, Frédéric; Soler, Charles; Houdouin, Véronique; Dalle, Jean-Hugues; Robin, Christine; Fourati, Slim; Griscelli, Franck; Coman, Tereza; Chevret, Sylvie; Bergeron, Anne title: Evaluation of prognostic scores for respiratory syncytial virus infection in a French multicentre cohort of allogeneic haematopoietic stem cell transplantation recipients date: 2021-09-21 journal: Bone Marrow Transplant DOI: 10.1038/s41409-021-01462-z sha: 5cde244b875ff34731099c1e9052f6d75224ac14 doc_id: 21159 cord_uid: ttbai8h2 Haematopoietic stem cell transplantation (HSCT) recipients are at risk for severe respiratory syncytial virus (RSV) infection. Two prognostic scores have been proposed to predict the risk of progression from upper respiratory tract infection (URTI) to lower respiratory tract infection (LRTI) and death. This was a multicentre study of allogeneic HSCT recipients diagnosed with an RSV infection between 2010 and 2019 who were retrospectively stratified by the immunodeficiency scoring index (ISI) and the severe immunodeficiency (SID) score. Endpoints were overall survival, RSV-attributable mortality and progression to LRTI after URTI. Prognostic analyses were performed using Cox regression models. We included 147 consecutive patients, including 94 (63.9%) initially diagnosed with URTI and 53 (36.1%) with LRTI. At 90 days, 14 patients had died (survival rate, 90.5%; 95% CI: 85.9–95.3), and nine deaths were attributable to RSV (attributable mortality rate, 5.4%; 95% CI: 2.5–10.0). The cumulative 90-day incidence of LRTI after URTI was 13.8% (95% CI: 7.8–21.6). Neither score showed prognostic value for mortality, while the ISI allowed the prediction of progression to LRTI (p = 0.0008). Our results do not fully replicate the results previously reported in cohorts of HSCT recipients. This may reflect the recent epidemiology of RSV infections in this HSCT cohort. In the course of allogeneic haematopoietic stem cell transplantation (HSCT), respiratory syncytial virus (RSV) infection can lead to severe infection with a poor prognosis [1] [2] [3] [4] [5] . RSV infection is seasonal, with recurrent winter outbreaks with various annual incidences ranging from 2 to 49% [6, 7] in allogeneic HSCT recipients. Upper respiratory tract infection (URTI) is usually minor, but the risk of progression to lower respiratory tract infection (LRTI) is high, especially in the early phase of HSCT, and it has been associated with mortality [7] [8] [9] [10] . Although the widespread use of molecular diagnostics allows earlier diagnosis and HSCT postponement in the case of a pretransplant RSV-positive test [11] , the mortality rate attributable to RSV remains significant, being approximately 10% in adults [9, 12] and <5% in children [13] in recent studies. The identification of risk factors for progression to LRTI and death in allogeneic HSCT patients has been the subject of several publications [4, 14, 15] , with the aim of defining the population of patients who would benefit from antiviral treatment. Two severity scores have been developed to classify patients according to the risk for progression to LRTI and the risk of RSVattributable mortality [11, 16] . The "Immunodeficiency Scoring Index" (ISI) from the MD Anderson Cancer Center was established in 2014 [16, 17] , and the "Severe Immuno-Deficiency" (SID) score was proposed by a single-institutional study from the University of Basel in 2008 [11] . These two grading systems aim to classify patients into 3 risk categories (low, moderate and high for the ISI and moderate, severe, and very severe for the SID score), relying on both clinical and biological criteria (Supplementary Table 1) . Current treatments for RSV are limited, but the guidelines of the 8th European Conference on Infections in Leukaemia [18] recommend systemic or aerosolized ribavirin, a broad-spectrum molecule used against many RNA viruses, with or without the associated use of intravenous polyclonal immunoglobulins (IVIGs) in both cases of LRTI and URTI with a high risk of progression to LRTI. Indeed, the pooling of cohort studies assessing the efficacy of ribavirin has suggested a reduction in both LRTI progression and mortality [6, 17, [19] [20] [21] [22] . However, meta-analyses have included studies that are difficult to compare, especially studies involving different RSV diagnostic techniques and heterogeneous management strategies [11] . Only one randomized controlled study was conducted to investigate the efficacy of aerosolized ribavirin, but difficulty in enrolling patients led to early termination; it was, therefore, impossible to draw a firm conclusion [23] . Other molecules are being developed, in particular presatovir, an RSV fusion inhibitor. A randomized placebo-controlled study failed to demonstrate a significant antiviral effect and a reduction in LRTI progression in HSCT patients with RSV-URTI, although post hoc analysis found an effect specifically in patients with lymphopenia [24] . Similarly, presatovir did not improve virologic or clinical outcomes compared with placebo in HSCT patients with RSV-LRTI [25] . At a time when new molecules are being evaluated, similar future trials should judiciously select suitable at-risk patients to maximize the potential benefits. Because LRTI increases the risk of mortality, the prompt diagnosis of and early intervention for RSV-URTI in high-risk patients with effective antiviral agents are imperative to improve the clinical outcomes. The aim of our study was to assess the prognostic value of both the ISI and the SID score using a recent multicentre cohort of allogeneic HSCT recipients. Both paediatric and adult allogeneic HSCT patients diagnosed with RSV infection by polymerase chain reaction (PCR) at the following five French hospitals were retrospectively included in the study: Saint Louis Hospital (Paris), Robert Debré Hospital (Paris), Gustave Roussy Institute (Villejuif), Percy Hospital (Clamart), and Henri Mondor Hospital (Créteil). All RSV-positive respiratory samples collected in the virology laboratory of each site were screened from the date of implementation of PCR (depending on each site) until December 2019. PCR for RSV was being used routinely first at Saint Louis Hospital and Robert Debré Hospital in December 2010 and January 2011, respectively; Percy Hospital, Henri Mondor Hospital, and Gustave Roussy Institute have been using it since December 2011, January 2015, and January 2016, respectively. Clinical data were collected from patient's medical charts. All chest imaging data were reviewed (ALH, LB, AB). Patients were included if PCR for RSV yielded a positive result from 10 days prior to HSCT to any time after transplantation. This study was approved by the Institutional Review Board of the French Learned Society for Respiratory Medicine, Société de Pneumologie de Langue Française (CEPRO 2019-040), and the data collected was anonymized prior to processing (Redcap software). The primary endpoint was overall survival (OS). Secondary endpoints were early mortality (any death within 90 days after RSV infection), RSV-attributable mortality, and, for patients with URTI, progression to LRTI. The diagnosis of URTI was defined by detection of the virus in the upper respiratory tract with symptoms of URTI (rhinitis, pharyngitis, laryngitis, conjunctivitis, sinusitis) without hypoxemia or infiltration on chest imaging. RSV infections were classified according to the ECIL8 definitions [18] . Briefly, the diagnosis of LRTI was defined by the association of respiratory symptoms with (i) pulmonary infiltration on chest imaging (within five days before and after infection diagnosis) or hypoxemia and (ii) the identification of RSV in the upper or lower respiratory tract. LRTI was classified as "probable" or "presumptive" LRTI when the virus was detected in the upper airway (i.e. by nasal swab, nasopharyngeal lavage, or nasopharyngeal aspiration) or in the lower respiratory tract (i.e. by endotracheal aspiration or bronchoalveolar lavage), respectively. Progression of URTI to LRTI was defined by new-onset lung infiltration on radiology, hypoxemia, or lower respiratory tract virus detection in a patient who initially presented with URTI. Coinfection was defined by the identification of another pathogen (bacterial, viral, or fungal) within 3 days before or after a positive RSV-PCR result. Data regarding coinfection at the time of progression to LRTI were also collected. An episode was considered multiple episodes when there was a four-week time interval between two positive PCR results for RSV or a negative PCR result between two episodes. Nosocomial RSV was defined by the occurrence of RSV infection at least 72 h after the patient's hospital admission, provided that signs and symptoms related to the infection were not present at the time of admission. Lymphoproliferative diseases included lymphoma, chronic lymphoid leukaemia, and myeloma. All-cause mortality was defined as any death within the 90 days following the diagnosis of RSV infection. RSV-attributable mortality was defined as death due to respiratory failure with no cause other than identified RSV. Acute or chronic graft-versus-host disease (GVHD) was evaluated at the time of RSV diagnosis [26] . ISI and SID score All patients were classified according to the criteria of both the ISI and the SID score at the time of RSV diagnosis (Supplementary Table 1 ). The ISI ranged from 0 to 12, calculated as follows [16, 17] : three points for either neutropenia ≤0.5 × 10 9 /l or lymphopenia ≤0.2 × 10 9 /l; two points for age ≥40 years; and one point each for myeloablative conditioning regimen, GVHD, corticosteroid therapy within 30 days before RSV infection, recent HSCT (less than 1 month prior) or no engraftment. Allogeneic HSCT recipients were then stratified into risk levels for severe RSV infection: low risk (ISI, 0-2); moderate risk (ISI, 3-6); and high risk (ISI, 7-12). The SID score was also used to classify patients into three categories: moderate immunodeficiency (MID); severe immunodeficiency (SID); and very severe immunodeficiency (verySID). MID was defined as follows: HSCT > 6 months prior; B or T lymphocyte depletion >3 months; neutropenia >0.5 × 10 9 /l and lymphocyte count >0.1 × 10 9 /l; GVHD < grade 2 or untreated; and hypogammaglobulinemia >6.5 g/L. In contrast, SID was defined by one and verySID was defined by two or more of the following risk factors: HSCT ≤ 6 months prior; B or T lymphocyte depletion ≤3 months; neutropenia ≤0.5 × 10 9 /l; lymphopenia ≤0.1 × 10 9 /l; GVHD grade ≥2 or treated; and hypogammaglobulinemia ≤6.5 g/L [11, 27] . Descriptive statistics, namely, the median with the interquartile range (IQR), or percentages are reported. Only the first episode of infection was considered for survival analyses, starting from the date of diagnosis to the date of death or the last follow-up. The OS distribution was estimated by the Kaplan-Meier method, while mortality attributable to RSV was estimated by cumulative incidence curves in a competing-risks setting where death from other causes was defined as a competing risk event. For the cumulative incidence of LRTI, only subjects with URTI were considered, with the death without progression defined as a competing risk outcome. Prediction models were cause-specific Cox models, with the cause-specific hazard ratio (HR) as the measure of the effect reported with the 95% confidence interval (CI). Univariable models were fitted, and multivariable models were used, including factors associated with the outcome at the 10% level. A stepwise selection procedure was used based on the Akaike criterion. Log-linear relationships were checked using generalized additive models with splines, while proportional hazards assumptions were checked by the Grambsch and Therneau statistical test [28] . Last, we reran the analyses specifically on the adult patients (age ≥ 18 yrs), given the SID score has been only developed on such a population [11] while the cohort used for the ISI also included children [16, 17] . All the analyses were carried out in R 3.6.2 (http://www.Rproject.org). Over a nine-year period, 251 samples from allogeneic HSCT patients tested positive for RSV by PCR, including samples from 147 patients who developed a total of 154 RSV episodes: 140 patients developed one episode, and 7 patients developed 2 episodes. The study population included adults and children (age < 18, n = 29), with a median age at the time of HSCT of 40.2 years (IQR: 22-56.6). Baseline characteristics are summarized in Table 1 . Of these 147 patients, 94 (63.9%) presented with URTI and 53 with LRTI (36.1%) at the time of diagnosis with RSV infection. There were 11 (7.5%) cases of nosocomial infection (4 LRTI and 7 URTI cases). The characteristics of the first episode of RSV infection are shown in Table 2 . The median follow-up was 23 months (IQR: 9.1-36.6). A total of 53 deaths occurred, including 14 within the first 90 days after RSV diagnosis, with an estimated 90-day survival of 90.5% (95% CI: 85.9-95.3%) (Fig. 1 ). Nine deaths were attributable to RSV, all of which occurred within the first 90 days. The cumulative RSVattributable mortality rate at 90 days was 5.4% (95% CI: 2.5-10.0) (Fig. 2) . There was no significant difference in survival between LRTI and URTI (Fig. 1) . The stratification of disease severity according to either the ISI or SID score did not demonstrate any difference in either overall mortality or RSV-attributable mortality both for the entire cohort and those restricted to the adults ( The results of the prognostic analyses are reported in Table 3 . In the univariate analyses of 90-day survival, the factors selected at the 5% level were the same, with the addition of hypogammaglobulinemia (HR = 3.91; 95% CI: 1.09-14.01; p = 0.036). In the multivariable analysis, only age was added to the other prognostic information, encompassing the whole follow-up period ( Table 3 , Supplementary Fig. 1 ). When restricting on adult patients, only leukocytes count below 2 Giga/L (HR = 2.56; 95%CI: 1.20-5.56; p = 0.015),and use of corticosteroids within 30 days of RSV infection (HR = 2.16; 95%CI: 1.22-3.84; p = 0.009) were associated with the outcome. Of the 94 patients with URTI, 13 progressed to LRTI, with a cumulative incidence of progression of 13.8% at 90 days (95% CI: 7.8-21.6) (Fig. 3) . The ISI identified three groups of patients with a specific risk for progression to LRTI (p = 0.0008). In contrast, the SID score could not predict progression to LRTI (p = 0.43) (Fig. 3) . The results of the univariable analyses identifying risk factors for progression are reported in Figure 4) . Our results show that the ISI but not the SID score was predictive of progression from RSV-URTI to LRTI in a multicentre retrospective cohort of allogeneic HSCT recipients. However, we found that neither of these scores predicted either overall mortality or RSV-attributable mortality. RSV-associated mortality was the primary endpoint of previous studies, including overall RSV-associated mortality for the SID score and 90-day RSV-associated mortality for the ISI. Because it is usually difficult to establish a single cause of death in allogeneic HSCT patients, in whom mortality is often multifactorial and involves numerous associated complications [21, 29] , classification errors are probable in retrospective cohorts. Therefore, as carried out in previous studies [9, [30] [31] [32] , we primarily analysed OS and truncated survival at 90 days as prognostic factors. We subsequently extended the study to include RSV-associated mortality for comparison with previous studies. Our study highlights the recent epidemiology of RSV infection in allogeneic HSCT patients. We found an overall mortality rate of 9.5% at 90 days and a cumulative RSV-attributable mortality rate of 5.4%. Furthermore, few patients required transfer to the intensive care unit (ICU). In studies conducted in the 1990s-2000s, the RSV-attributable mortality rate was between 30% and 80% [10, 33, 34] , whereas more recent studies found an RSVattributable mortality rate of 6-10% [9, 12, 30] . In the specific cohorts of the ISI and the SID score, the RSV-attributable mortality rate was 8% and 18%, respectively. The differences between our cohort and both the SID score and ISI cohorts may explain the lack of full validation of both scores. Both the ISI and the SID score are based on older cohorts of All patient received rituximab (10 patients within the three months prior to RSV infection and four patients between four months and six months prior to RSV infection). The diagnosis of RSV infection was based exclusively on viral culture or direct immunofluorescence antibody testing (DIFT) in the ISI cohort, whereas only RSV-PCR, which is actually used in clinical practice, was used for diagnosis in our cohort. More than 90% of the patients from the SID score cohort underwent RSV-PCR, but the study included only 34 patients, with 27 HSCT recipients [11] . PCR, used in the most recent studies, is much more sensitive than viral culture or DIFT for the identification of any community-acquired respiratory virus (CARV) and allows the detection of lower viral loads [11] . Therefore, it enables the diagnosis of milder RSV infections. Furthermore, testing patients before they undergo HSCT during RSV epidemics allows the delay of HSCT if the test is positive, thereby preventing severe RSV infection occurring early post-HSCT [35] . Thus, more recent studies include a higher proportion of URTI cases than older studies, including more RSV-LRTI cases, with higher mortality rates [36, 37] . We included only allogeneic HSCT patients, including both children and adults, in contrast to the SID score cohort, which considered only adults treated with both autologous and allogeneic HSCT. The prognosis of RSV infection in autologous HSCT patients compared to allogeneic HSCT patients is discordant in the literature [14] ; some studies found no difference [15, 21, 31] , while others showed increased severity in allogeneic HSCT patients [38, 39] . Data in paediatric allogeneic HSCT patients are limited. However, while children most often receive myeloablative conditioning, the RSV-attributable mortality rate seems to be lower in children than in adults (<5%) [40] . The characteristics of our patients also differed somewhat from those in the ISI cohort. First, the median time from HSCT to RSV infection was longer in our study (196 days) than in the ISI study (155 days). Second, our patients were younger (40.2 vs 47 years). Finally, we included patients with LRTI, whereas the ISI cohort focused on patients with URTI. Other studies have assessed these scores in cohorts of allogeneic HSCT recipients infected with various CARVs, including variable proportions of patients infected with RSV [8] . The results regarding the validation of these scores are discordant among the different studies regarding both mortality and progression to LRTI [9, 12, 30-32, 37, 41, 42] . Only one study focused on 45 RSV infections in adult allogeneic HCT recipients exclusively diagnosed by PCR [42] . The authors could validate the ISI for mortality but not for LRTI progression [42] . We found that only advanced age was an independent risk factor for overall and 90-day mortality in our entire cohort. Age was taken into consideration only in the ISI and not in the SID score [11, 16] . We identified other risk factors for mortality in univariable analyses that are consistent with the variables included in either the ISI or SID score, i.e. GVHD and corticosteroid therapy. We also found that hypogammaglobulinemia ≤6.5 g/l was associated with 90-day mortality, whereas this variable is only part of the SID score and not the ISI. Surprisingly, we did not identify low lymphocyte or neutrophil counts as prognostic factors for either overall or 90-day mortality, whereas they were associated with progression to LRTI. Similarly, we did not find any difference in prognostic impact between URTI and LRTI. The lack of a difference in mortality between LRTI and URTI can be explained by our low rate of presumptive LRTI among all LRTI cases (5.7% vs 58.3% in the SID score cohort). Waghmare et al. showed that possible LRTI can be misleading and that the prognosis was worse in presumptive LRTI [43] . In contrast to other studies, we found that myeloablative conditioning was associated with better survival. Our discordant result could be attributable to the widespread use of nonmyeloablative conditioning regimens in adults, representing 60% of our entire cohort and 73% of our adult population. Only 58% of patients in the ISI cohort and 33% in the SID score cohort underwent nonmyeloablative conditioning. Notably, 46% of our patients who underwent myeloablative conditioning were less than 18 years old, and no deaths were observed in this young population. Multivariate analysis confirmed that both moderate-and highrisk patients according to the ISI were at risk for progression to LRTI. The prediction of LRTI progression was not significant for the SID score, although no patients progressed to LRTI in the lowerrisk group. Because RSV infections in pediatric allogeneic HCT recipients may differ from those occurring in adult HCT recipients, we reran the analyses by removing the 29 children from our cohort. Of note, the ISI cohort included children but not the SID score cohort. Results did not change regarding both predictions for mortality and LRTI. However, a low leukocytes count was identified as another factor associated with mortality. Based on a recent, large multicentric retrospective cohort of allogeneic HSCT recipients with PCR-confirmed RSV infection, we confirmed a reduction in morbidity and mortality compared with historical cohorts. This result could explain why we could not fully validate either the ISI or the SID score, specifically regarding the prediction of mortality. However, we confirmed that immunosuppressive factors were determinants for both progression from URTI to LRTI and mortality. The impact of RSV treatment was not found to be a significant prognostic factor, but our study was not appropriately designed to assess this outcome. Our study provides recent information on cases of RSV infection diagnosed by PCR and should help to better design future trials. All patients received rituximab (within the three months prior to RSV infection for x patients and between four months and six months prior to RSV infection for y patients). c Ribavirine and/or intravenous polyclonal immunoglobulins. A.-L. Houist et al. Survival in neutropenic patients with severe sepsis or septic shock* Respiratory syncytial virus: infection, detection, and new options for prevention and treatment Respiratory virus infection in immunocompromised patients Timeline, epidemiology, and risk factors for bacterial, fungal, and viral infections in children and adolescents after allogeneic hematopoietic stem cell transplantation Virus infections post transplant: risk and immunity Management of respiratory viral infections in hematopoietic cell transplant recipients and patients with hematologic malignancies Management of respiratory viral infections in hematopoietic cell transplant recipients Community-acquired respiratory viruses in transplant patients: diversity, impact, unmet clinical needs Community-acquired respiratory paramyxovirus infection after allogeneic hematopoietic cell transplantation: a single-center experience Prospective study of the incidence, clinical features, and outcome of symptomatic upper and lower respiratory tract infections by respiratory viruses in adult recipients of hematopoietic stem cell transplants for hematologic malignancies Respiratory syncytial virus infection in patients with hematological diseases: single-center study and review of the literature Respiratory virus infection after allogeneic hematopoietic stem cell transplant in a tropical center: Predictive value of the immunodeficiency scoring index A multicenter investigation of respiratory syncytial viral infection in children with hematopoietic cell transplantation Respiratory syncytial virus in hematopoietic cell transplant recipients: factors determining progression to lower respiratory tract disease Respiratory syncytial virus in hematopoietic cell transplant recipients: factors determining progression to lower respiratory tract disease Immunodeficiency scoring index to predict poor outcomes in hematopoietic cell transplant recipients with RSV infections Impact of aerosolized ribavirin on mortality in 280 allogeneic haematopoietic stem cell transplant recipients with respiratory syncytial virus infections Update on Community-acquired respiratory viruses in hematology patients -Final Slide Set Management of RSV infections in adult recipients of hematopoietic stem cell transplantation Oral ribavirin for treatment of respiratory syncitial virus and parainfluenza 3 virus infections post allogeneic haematopoietic stem cell transplantation Risk factors and containment of respiratory syncytial virus outbreak in a hematology and transplant unit Oral ribavirin is a highly effective treatment for lower respiratory tract infections due to respiratory syncytial virus or parainfluenza after allogeneic stem cell transplantation Randomized controlled multicenter trial of aerosolized ribavirin for respiratory syncytial virus upper respiratory tract infection in hematopoietic cell transplant recipients A phase 2, randomized, double-blind, placebo-controlled trial of presatovir for the treatment of respiratory syncytial virus upper respiratory tract infection in hematopoietic-cell transplant recipients A phase 2b, randomized, double-blind, placebo-controlled multicenter study evaluating antiviral effects, pharmacokinetics, safety, and tolerability of presatovir in hematopoietic cell transplant recipients with respiratory syncytial virus (RSV) infection of the lower respiratory tract National institutes of health consensus development project on criteria for clinical trials in chronic graft-versus-host disease: I. the 2014 diagnosis and staging working group report Outcome of influenza infections in outpatients after allogeneic hematopoietic stem cell transplantation Proportional hazards tests and diagnostics based on weighted residuals Burden of severe RSV disease among immunocompromised children and adults: a 10 year retrospective study The effect of timing on community acquired respiratory virus infection mortality during the first year after allogeneic hematopoietic stem cell transplantation: a prospective epidemiological survey Risk factors for mortality after respiratory syncytial virus lower respiratory tract infection in adults with hematologic malignancies Respiratory syncytial virus and human metapneumovirus after allogeneic hematopoietic stem cell transplantation: Impact of the immunodeficiency scoring index, viral load, and ribavirin treatment on the outcomes Community respiratory virus infections in bone marrow transplant recipients: the M.D. Anderson Cancer Center experience An outbreak of respiratory syncytial virus in a bone marrow transplant center Pretransplantation respiratory syncytial virus infection: impact of a strategy to delay transplantation Respiratory syncytial virus in hematopoietic cell transplant recipients and patients with hematologic malignancies Community acquired respiratory virus infections in adult patients undergoing umbilical cord blood transplantation Respiratory virus infections after stem cell transplantation: a prospective study from the Infectious Diseases Working Party of the European Group for Blood and Marrow Transplantation Respiratory viral infections after bone marrow/peripheral stem-cell transplantation: the Christie hospital experience 1 A multicenter consortium to define the epidemiology and outcomes of inpatient respiratory viral infections in pediatric hematopoietic stem cell transplant recipients Community-acquired respiratory virus lower respiratory tract disease in allogeneic stem cell transplantation recipient: Risk factors and mortality from pulmonary virus-bacterial mixed infections Corticosteroid use as adjunct therapy for respiratory syncytial virus infection in adult allogeneic stem cell transplant recipients Supplemental oxygen-free days in hematopoietic cell transplant recipients with respiratory syncytial virus The authors thank Stéphane Cassonnet for his technical assistance. The authors thank EGMOS, a patient association, for their support. ALH, SC and AB conceived the study. ALH, LB, MS, JL, SC and AB designed the work that led to the submission. ALH, LB, MS, JL, RPL, FR, CS, VH, JHD, CR, SF, FG, TC, SC and AB acquired data, and played an important role in interpreting the results. ALH, SC and AB drafted the manuscript. All authors revised the manuscript, approved the final version and agreed to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. The authors declare no competing interests. Supplementary information The online version contains supplementary material available at https://doi.org/10.1038/s41409-021-01462-z.Correspondence and requests for materials should be addressed to Anne Bergeron.Reprints and permission information is available at http://www.nature.com/ reprintsPublisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.