Abstract
Allogeneic hematopoietic stem cell transplantation (HSCT) using αβ T-/B-cell–depleted grafts recently emerged as a transplant strategy and highlighted the potential role of γδ T cells on HSCT outcomes. Our aim was to scrutinize available evidence of γδ T-cell impact on relapse, infections, survival, and acute graft-versus-host disease (aGVHD). We performed a systematic review and meta-analysis of studies assessing γδ T cells in HSCT. We searched PubMed, Web of Science, Scopus, and conference abstracts from inception to March 2019 for relevant studies. We included all studies that assessed γδ T cells associated with HSCT. Data were extracted independently by 2 investigators based on strict selection criteria. A random-effects model was used to pool outcomes across studies. Primary outcome was disease relapse. We also assessed infections, survival, and aGVHD incidence. The review was registered with PROSPERO (CRD42019133344). Our search returned 2412 studies, of which 11 (919 patients) were eligible for meta-analysis. Median follow-up was 30 months (interquartile range, 22-32). High γδ T-cell values after HSCT were associated with less disease relapse (risk ratio [RR], 0.58; 95% confidence interval [95% CI], 0.40-0.84; P = .004; I2 = 0%), fewer viral infections (RR, 0.59; 95% CI, 0.43-0.82; P = .002; I2 = 0%) and higher overall (HR, 0.28; 95% CI, 0.18-0.44; P < .00001; I2 = 0%) and disease-free survivals (HR 0.29; 95% CI, 0.18-0.48; P < .00001; I2 = 0%). We found no association between high γδ T-cell values and aGVHD incidence (RR, 0.72; 95% CI, 0.41-1.27; P = .26; I2 = 0%). In conclusion, high γδ T cells after HSCT is associated with a favorable clinical outcome but not with aGVHD development, suggesting that γδ T cells have a significant effect on the success of HSCT. This study was registered with PROSPERO as #CRD42019133344.
Introduction
Allogeneic hematopoietic stem cell transplantation (HSCT) offers a potential cure for a variety of life-threatening hematological diseases. The most common causes of posttransplantation mortality are underlying disease relapse, infections, and graft-versus-host disease (GVHD). In this context, graft manipulation strategies have developed and been refined over the past 20 years to improve HSCT outcomes, moving from CD34 selection to CD3 depletion and, more recently, to αβ T-/B-cell depletion.1 These strategies aim to reduce GVHD while retaining cells that can mediate graft-versus-leukemia (GVL) effect and control infections. Several concerns exist with these strategies, primarily whether remaining cells will be sufficient to prevent disease recurrence and infections without causing GVHD.
γδ T cells are a unique population of lymphocytes that mediate innate immunity against a wide variety of infections and exert effective antitumor activity.2 Understanding the role of γδ T cells in HSCT has been the subject of numerous studies in the past decade.1,3 Early studies observed that long-term disease-free survival (DFS) of leukemia patients who received αβ T-cell–depleted (TCD) partially mismatched related donor (PMRD) transplants had high numbers of circulating γδ T cells after HSCT.4,5 More recent studies with non-TCD PMRD6 and autologous HSCT7 further corroborated these findings, describing an improved survival in patients with higher γδ T-cell counts posttransplantation. These studies suggested that the recovery of γδ T cells after HSCT is critical for an efficient GVL effect3 and possibly to control infections,8 helping to pave the way to a more broad application of graft manipulation strategies. As result, haploidentical αβ T-/B-cell–depleted HSCT is currently used with great success to treat children with malignant9,10 and nonmalignant disorders.11 These works argue that the clinical improvement after HSCT relies on the spared γδ T cells in the grafts, which exhibit strong antileukemia potential without mediating GVHD,9-11 and can participate in the control of opportunistic infections.3 Indeed, recent studies indicate that γδ T cells are capable of adaptive responses and undergo clonal expansion after cytomegalovirus (CMV) reactivation,12 but it is unclear whether this decreases infection frequency after HSCT.
Preclinical studies from the late 1990s has shown a beneficial role of γδ T cells on improving allogeneic engraftment without causing lethal GVHD,13,14 indicating potential benefits of this T-cell subpopulation on HSCT setting. Conversely, several mice studies suggest that γδ T cells contribute toward GVHD development,15,16 and a key study described that the graft content of donor γδ T cells predicted the risk of acute GVHD (aGVHD).17 Altogether, these works propose that γδ T cells might be beneficial in HSCT but with the cost of higher incidence of GVHD. Considering that αβ T-/B-cell–depleted HSCT is increasingly used in patients with no matched donor9-11 and that γδ T cell–enriched donor lymphocyte products are being under current investigation,18 consistent evidence is needed that supports beneficial effects of γδ T cell enrichment with no detrimental effect to the patients. We aimed to determine whether the concentration of γδ T cells in the graft or during immune reconstitution influenced the clinical outcome following HSCT and aGVHD incidence.
Methods
Search strategy and selection criteria
This systematic review and meta-analysis was performed following the PRISMA guidelines for conducting and reporting systematic reviews.19 In March 2019, we conducted a comprehensive literature search for potentially relevant studies of γδ T cells and their effect on HSCT outcomes (relapse, infections, overall and DFS, and aGVHD) with no publication time limits. Inclusion criteria were: (1) original studies (randomized, cohort, case-control, prospective and retrospective observational studies) enrolling adult or pediatric patients who had undergone peripheral blood, umbilical cord blood, or bone marrow HSCT (allogeneic and autologous) as therapy for any condition (malignant or not); and (2) have had γδ T cells (absolute count, percentage, or subsets) measured in the graft before transplantation or during immune reconstitution process.
Studies that reported dichotomous data (high vs low γδ T-cell groups) together with number of events (relapse, infections, or aGVHD) and time-to-event data (overall or DFS) during the follow-up period were included in the meta-analysis, whereas those who reported associative data were included in the qualitative synthesis. We excluded reviews, studies that reported insufficient data, experimental in vitro studies of expanded or modified γδ T cells, animal models, or studies in which presenting outcomes were not relevant to the searching protocol. Articles not in English were excluded if translations of the abstracts were not available.
Two authors (L.C.M.A. and A.G.) searched PubMed/Medline, Web of Science, and Scopus, in duplicate and independently of each other with the following search terms: “(hematopoietic stem cell transplantation OR HSCT) AND (γ δ OR γδ)” and “(bone marrow transplantation OR BMT) AND (gamma delta OR γδ).” We also checked for unpublished relevant studies and conference abstracts of the American Society of Blood and Marrow Transplantation, the European Group for Blood and Marrow Transplantation, and the American Society of Hematology. We checked references of selected publications for additional potentially relevant studies.
Two investigators (L.C.M.A. and A.G.) assessed all studies’ eligibility based on title and abstract. Potentially eligible studies were retrieved, and the full study report evaluated. We resolved disagreements by consensus or discussion, or they were adjudicated by a third reviewer (M.U.). We contacted study authors to clarify results details when necessary; if no answer was obtained or these data were not available, the record was excluded.
Data analysis
Data were extracted onto data extraction forms by 2 reviewers independently (L.C.M.A. and A.G.) and included author name, year of publication, journal, number of subjects, donor and patient age/sex, underlying disease, donor type, graft source, conditioning regiment, GVHD prophylaxis, median and range of follow-up, timepoint, γδ T-cell phenotype, and outcomes. Groups definition was based on γδ T-cell content reported by the works: high or low γδ T-cell numbers. We extracted the number of events in each group by annotating the number of relapses, infection reactivation events after HSCT, number of patients alive (overall and disease-free) and aGVHD events in each group.20 The number of subjects in each group is denoted as “total” throughout the figures. Duplicate publications and meeting abstracts were not included in the final selection. We did subgroup analyses on the basis of sample origin: blood, when γδ T-cell numbers were obtained in patients’ blood after HSCT; and grafts, when the intragraft γδ T-cell content was assessed and associated with the outcome. For each study included in the meta-analysis, we extracted the number of events in each group and the total number of subjects. In case of time-to-event analysis, we extracted the hazard ratios (HR) of high vs low γδ T-cell counts effect on overall survival (OS) and DFS. If not available, HR was estimated from published summary statistics by using the spreadsheet provided by Tierney et al.21 Disagreements on data extraction were resolved by consensus with the supervision of a third reviewer (M.U.).
To judge study quality, we used for the current systematic review a modified Newcastle-Ottawa Scale and Research Triangle Institute Item Bank to assess the risk of bias and confounding factors in observational studies,22,23 as recommended by the Cochrane Collaboration (supplemental Appendices).19,20 Risk of bias was assessed by 2 reviewers (L.C.M.A. and A.G.) independently.
Underlying disease relapse was defined as the primary outcome of this review. Secondary outcomes were infections, OS, DFS, and aGVHD incidence. Meta-analyses were performed using Review Manager (RevMan),24 version 5.3, by Mantel-Haenszel25 method for dichotomous data, or generic Inverse Variance26 for time-to-event outcome. We used random effects to calculate pooled risks on the basis of assumption that the true effect size of γδ T cells effect would vary between studies.27 Results are shown as risk ratio (RR) or HR with 95% confidence interval (CI), as calculated following the Cochrane handbook.20 The RR and HR for immune reconstitution studies and graft subgroup evaluations were calculated separately. Then, we grouped all studies to calculate a pooled risk effect for each outcome. The overall effect significance was calculated by the z test,20 with P < .05 set as significant.
We assessed heterogeneity in the meta-analysis with the I2 statistics.28,29 The test use χ2 and degrees of freedom to describe the percentage of the variability in effect estimates that is due to heterogeneity rather than sampling error (chance).20 I2 reflects the percentage of total variation across studies, and values greater than 25%, 50%, or 75% were considered to respectively indicate low, moderate, or high heterogeneity.29 If P < .05, the pooled analysis was considered significantly heterogeneous.20 We also used τ2 to estimate the dispersion of true effect sizes between studies, with low values meaning low dispersion and consequently low heterogeneity.20
Results
The initial literature search found 2412 potentially eligible records. After removing duplicates and screening titles/abstracts, we fully reviewed 78 reports, of which 43 were excluded (Figure 1). We included 24 studies (21 full papers and 3 meeting abstracts) in qualitative synthesis,12,30-52 summarized in Table 1. Eleven studies were used in meta-analysis, enrolling 919 patients.4-8,17,53-57 From those reports, 8 evaluated γδ T-cell reconstitution after HSCT and 3 evaluated γδ T-cell content in the graft (Table 2).
All meta-analyzed studies were single cohorts that followed the patients for a median follow-up of 30 months. Most enrolled adult patients (median age, 32 years), males (median, 60%), with acute leukemias that received matched allogeneic bone marrow stem cells after a myeloablative conditioning (Table 2). GVHD prophylaxis was consistent among the studies and was based on cyclosporine and methotrexate. aGVHD incidence was reported in 7 studies,4,6,17,53-55,57 and the number of relapses was reported in 7 studies.4,5,7,53-56 Viral infections after HSCT were reported in 6 studies; CMV reactivation was reported in 5,4,6,53,55,56 whereas Epstein-Barr virus (EBV) reactivation in 1 study.8 OS was reported in 5 studies,4,7,53,55,56 whereas DFS was reported in 4.4,5,54,55 γδ T-cell content was defined primarily by anti-pan γδ T-cell receptor marker, despite some studies using anti-Vδ28 or anti-Vδ1.6,57 Six studies stratified patients using γδ T-cell percentage,5,6,53,54,56,57 whereas 5 used absolute counts.4,7,8,17,55 Three studies checked the intragraft γδ T-cell content,17,53,57 whereas 8 studies evaluated its immune reconstitution after HSCT and used 100 days posttransplantation as the median sampling timepoint for group definition.4-8,54-56
The risk of bias and confounding assessment depicted that all studies reported well-documented patient’s baseline characteristics, were selected appropriately, and had outcome measures consistently defined across all participants (supplemental Table 1). γδ T-cell stratification was consistently performed based on median values, with the exception of 2 studies that arbitrarily defined the threshold5,54 and 2 other studies that used health donors’ median γδ T-cell distribution as the cutoff.4,55 Potential confounders were taken into account in most of studies, including disease risk category as a competing risk (supplemental Table 1).5,53,54 No study reported blinding medical practitioners to γδ T-cell data.
Patients with high γδ T-cells during immune reconstitution after HSCT were more likely to present less relapse than patients with low γδ T-cell values (RR, 0.58; 95% CI, 0.40-0.84; P = .004; Figure 2). All studies reported a positive association between high γδ T cells and less incidence of relapse (I2 = 0%, P = .54). If the autologous HSCT is removed, there is still statistical significance and low heterogeneity across studies (RR, 0.50; 95% CI, 0.28-0.89; P = .002; I2 = 0%; supplemental Figure 1). Two qualitative studies also described that γδ T-cell clonotypes were associated with less relapse after HSCT (Table 1).45,47 Only 1 study assessed the γδ T-cell graft composition and observed no effect on relapse incidence (RR, 0.99; 95% CI, 0.76-1.29; P = .95).53 The pooled risk effect of both immune reconstitution and graft content further confirmed improved outcome for patients with high γδ T cells (RR, 0.65; 95% CI, 0.42-1.29; P = .05), although there was evidence of subgroup heterogeneity (I2 = 81.2%, P = .02).
Higher γδ T-cell values after HSCT were also associated with lower incidence of viral infections (RR, 0.59; 95% CI, 0.43-0.82; P = .002; Figure 3). Statistical analysis revealed homogeneity of the data (I2 = 0%, P = .56). The sole study on grafts observed no correlation between γδ T-cell graft content and CMV reactivation (RR, 1.05; 95% CI, 0.78-1.42; P = .74).53 The pooled risk effect also indicates lower incidence of infections in patients with high γδ T cells after HSCT, although this was not significant (RR, 0.68; 95% CI, 0.45-1.02; P = .06). The studies included in the qualitative synthesis highlighted that the Vδ1 subtype mediates the antiviral effect31,32,38,48,51 and that distinct γδ T-cell clones are important in control of viral infection (Table 1).12,45
The OS and DFS follow-up period were not consistently reported among the studies, ranging from 2,7 2.5,54 3,5,56 4,53 and 5,55 up to 7 years.4 Only 1 study reported HR between high vs low γδ T-cell groups4 ; for all the others, we estimated the HR following standard guidelines.21 Patients presenting a higher count of γδ T cells after HSCT tended to experience higher OS (HR, 0.28; 95% CI, 0.18-0.44; P < .00001; Figure 4; Table 2) and DFS (HR, 0.29; 95% CI, 0.18-0.48; P < .00001; Figure 5). The heterogeneity was absent for both outcomes (I2 = 0%, P > .05). If the autologous HSCT is removed, there is still statistical significance and low heterogeneity across studies (HR, 0.23; 95% CI, 0.13-0.41; P < .00001; I2 = 0%; supplemental Figure 2). The intragraft γδ T-cell evaluation did not show a significant effect on OS (HR, 1.34; 95% CI, 0.59-3.05; P = .49), but depicted a high heterogeneity across studies (I2 = 62%, P = .02), although the overall effect remained statistically significant (HR, 0.36; 95% CI, 0.18-0.70; P = .003).
The effect of high γδ T-cell numbers after HSCT on aGVHD incidence was not significant (RR, 0.72; 95% CI, 0.41-1.27; P = .26; Figure 6), with no evidence of interstudy heterogeneity (I2 = 0%, P = .73). Intragraft γδ T-cell content studies reported discrepant results and presented high heterogeneity (I2 = 82%, P = .004). One study reported a higher incidence of aGVHD in the high γδ T-cell patient subgroup (RR, 1.70; 95% CI, 1.02-2.82),17 whereas more recent studies describe the protective role of these cells (RR, 0.30; 95% CI, 0.13-0.72)57 or that they do not affect aGVHD development (RR, 1.13; 95% CI, 0.68-1.88).53 The pooled risk effect further confirmed that high γδ T-cell content is not associated with aGVHD development (RR, 0.82; 95% CI, 0.50-1.35; P = .44), with high between-subgroup homogeneity (I2 = 0%, P = .66). Studies included on qualitative synthesis also indicated the lack of association between γδ T cells and aGVHD development (Table 1).30,35,36,40,41,43,45,46,48,49,52
Discussion
Our systematic review and meta-analysis show that higher numbers of γδ T cells in peripheral blood after HSCT is associated with less risk of relapse, fewer infection events, and higher survival, with no risk association with GVHD development. γδ T cells are a unique and conserved population of innate immunity lymphocytes that play key roles in immune surveillance and tissue homeostasis.58 They represent just a small fraction of circulating T cells, but display the ability to expand in response to infections12 and exert antitumor effect.3 In contrast to αβ T cells, γδ T cells are mainly CD4−/CD8− and are not HLA-restricted. Because donor-derived γδ T cells may exert GVL effect without causing GVHD, large-scale methods to enrich, isolate, expand, and manipulate these cells for HSCT application are in progress and will clarify their full function in this setting.3 This meta-analysis is the first to suggest that the use γδ T cells in HSCT might be beneficial.
Observations that T cells are the key mediators of GVHD development led clinicians to ex vivo deplete T cells through CD34+ cell selection or removal of CD3+ T cells.1 These approaches result in loss of certain cell subsets that may play a beneficial role in the recipient. CD34+ selection was associated with slow immune recovery and high infection rate,59 whereas the removal of CD3+ T cells presented high infection and relapse rates.60 In fact, although αβ T cells mediate GVHD development, γδ T cells have lower alloreactivity and contribute to important anti-infectious activity,58 in addition to a possible antileukemia role.2,3 Despite reports of a positive association between γδ T cells and less disease relapse,4,5,54 no strong and unbiased evidence exists that these cells are indeed key players for successful HSCT. Our systematic review and meta-analysis suggest that in TCD PMRD,4,5,54 matched unrelated donor (MUD),55 autologous,7 and TCD haplo-HSCT,56 increased reconstitution of γδ T cells are positively associated with a significantly decreased risk of relapse. These results are sustained by specific γδ T-cell clonotypes from the graft donor45 that can expand after HSCT and exert antileukemic effect,47 supporting the notion that these cells could be broadly used in HSCT through enrichment methods or by post-HSCT infusions. Our results support the initial observations that the prompt reconstitution of γδ T cells after αβ T-/B-cell–depleted HSCT might be associated with the transplantation efficacy.11 More recently, an Italian multicenter study showed that αβ T-/B-cell depletion presented better relapse-free survival than MRD and MUD HSCT, highlighting the role of γδ T cells in protecting the host against leukemia relapse.10 The utilization of γδ T cells to mitigate the risk of relapse and to enhance immune reconstitution after HSCT continues to be under investigation. The successful use of αβ T-cell–depleted cell products as stem cell boosters after HSCT in patients with poor graft function, primary graft failure, and/or infectious complications was recently reported.18
CMV viremia is one of the most life-threatening infections after HSCT, being associated with high nonrelapse mortality in the first months posttransplant. Several works have reported a possible protective role of γδ T cells on infections.37 The first report described that high γδ T cells were not only protective of leukemia relapse but also in the control of CMV.4 Latter studies have shown that the Vδ1 subpopulation expands after CMV reactivation, playing a role in infection control.31,33,34,38,51,52 Recently, it was shown that CMV reactivation leads to clonal proliferation of individual virus-reactive γδ T-cell receptor (TCR) sequences, suggesting an adaptive antiviral γδ T-cell immune response.12 The same was observed in grafts from CMV+ donors in which several CMV-associated γδ T-cell clones had clonally expanded.45 Here we show that, although intragraft γδ T-cell counts were not associated with fewer infections, they play a role in host infection control and increased numbers of these cells after HSCT might be warranted to protect patients from opportunistic infections. Recently, αβ T cell–depleted donor lymphocyte infusions to treat patients with poor immune reconstitution and infections were reported with encouraging results.18 Eleven of 12 patients presented favorable responses with no increased risk of GVHD development, indicating that γδ T-cell–enriched donor lymphocyte infusion products are a viable option to improve HSCT efficacy and safety after transplantation.
The first reports evaluating the role of γδ T cells in GVHD pathogenesis are from the early 1990s.30,41 These studies provided no evidence that γδ T cells could mediate the pathogenesis of gut, liver nor epithelial lesions associated with GVHD,30 that were, in fact, mediated by αβ T cells.41 In the 2000s, several complementary studies used flow cytometry and molecular techniques (TCR spectratype) to corroborate that they are not mediating GVHD pathogenesis,46,48 with some controversies.39,47,49 Recently, long-term observation of 80 children given αβ T-/B-cell–depleted grafts showed no severe (grade 3-4) aGVHD, further indicating that γδ T cells might not be involved in GVHD pathogenesis.9 Thirty percent of the patients presented skin-only grade 1-2 aGVHD, and no extensive chronic GVHD was reported, despite not receiving any GVHD prophylaxis. Additionally, in an Italian multicenter study, 98 children receiving the same HSCT protocol showed significantly less GVHD incidence than MUD- and mismatched unrelated donor–treated patients.10 These studies advocate that the low GVHD incidence is due to αβ T-cell graft depletion together with γδ T-cell enrichment, which quickly reconstitute posttransplantation and gives support to hematopoiesis without triggering GVHD. Pabst et al reported an increased risk of aGVHD development in patient with enriched γδ T cells in the graft,17 whereas other reports showed no relation between γδ T cells and GVHD.40,42,43,53,57 Here, we show that there is no association between γδ T-cell reconstitution and GVHD development.
This systematic review depicts that only a few studies have assessed γδ T cells in graft composition and their effect on clinical outcome. More studies are necessary to better describe the role of these cells. Across the 3 studies reported, there were divergent results regarding aGVHD incidence. Pabst et al reported an increased risk of aGVHD development in patient with enriched γδ T cells in the graft,17 whereas other reports showed no relation between γδ T cells and GVHD.40,42,43,53,57
γδ T cells quickly reconstitute after HSCT9,11,12,34 and, at 1 month after HSCT, can expand and compose around 80% of total T cells.34 These cells are cytotoxic and effectively kill leukemia.5,34 Ravens et al reported a heterogeneous overlap of γδ T-cell sequences between donor and recipients repertoires posttransplantation, indicating that donor-derived γδ T cells are able to contribute to recipient’s pool, together with the de novo generation of thymic-derived cells.12 Long-term maintenance of γδ T cells might be key to achieve the beneficial outcomes of less relapse and infection events here cited, and it is conceivable that grafts enriched in γδ T-cell content will result in higher counts after HSCT. Indeed, sparing γδ T cells from the graft through αβ T-cell depletion resulted in significantly higher reconstitution of these cells when compared with pan T-cell depletion.5 In contrast, total T-cell depletion resulted in impaired γδ T-cell reconstitution and less beneficial effects.5 The survival advantage associated with high circulating numbers of γδ T cells is durable over 7 years following HSCT.4
Our study has limitations. The number of patients is low and the threshold definitions of high and low γδ T cells were divergent across reports and dependent on the timepoint of assessment. If assessed on day 30 after HSCT, patients with increased γδ T-cell numbers were those with >70%.56 But, if assessed at 100 days after HSCT, >10% is considered a high value.5,54 We propose that absolute counts are better and that the threshold of 1.75 × 105 γδ T cells/mL is consistent across studies to define the high γδ T-cell content.4,55 Another limitation is the inclusion of heterogeneous transplantation methods, donors, and patients, which can jeopardize the effect of our findings. Of the 11 studies and 919 patients included in this meta-analysis, 378 (41%) patients were TCD transplants using PMRD, 172 (19%) other PMRD transplants, 251 (27%) MUD or MRD, and 101 (11%) autologous. The beneficial effects of γδ T cells were mainly from the 8 studies (721 patients) examining γδ T-cell reconstitution. Of this, only 70 patients (10%) were MUD/MRD, with the remainder being TCR, PMRD, or autografts. Additionally, of the 3 studies (198 patients) examining graft content, all were MUD/MRD and with essentially no impact of γδ T cells on infection or relapse. This indicates that, although γδ T cells are clearly a population of interest in HSCT, it is quite likely that their effect is context dependent, with more impact in the TCD setting. More studies are warranted to full address the role of these cells in HSCT.
In summary, our findings indicate that γδ T cells may play an important role in HSCT efficacy and safety, participating in both leukemia and infection control and resulting in higher survival of the patients, with no association with GVHD development.
Acknowledgment
The authors thank Olle Ringdén (Department of Laboratory Medicine, Karolinska Institute, Stockholm, Sweden) for the critical review of this manuscript.
Authorship
Contribution: L.C.M.A. searched the published work, produced the figures, collected, analyzed, and interpreted data, and wrote the report; A.G. searched the published work, collected and analyzed data, and wrote the report; M.U. analyzed and interpreted data and wrote the report; and all authors reviewed and revised the manuscript and approved the final version.
Conflict-of-interest disclosure: The authors declare no competing financial interests.
Correspondence: Lucas C. M. Arruda, Department of CLINTEC, Karolinska Institutet, ANA Futura, Alfred Nobel Allé 8, SE-141 52 Stockholm, Sweden; e-mail: lucas.arruda@ki.se.
References
Author notes
The full-text version of this article contains a data supplement.