TO THE EDITOR:

Multiple myeloma (MM) remains an incurable plasma cell (PC) malignancy.1,2 MM cells residing in the bone marrow (BM) engage in extensive crosstalk with the BM microenvironment. These interactions are leveraged by immunotherapies, but many patients acquire resistance, suffer dose-limiting toxicities, or remain unresponsive, highlighting an unmet need.

Natural killer (NK) cells demonstrate strong cytotoxicity against stressed cells, such as tumor cells and malignant PCs.3 The efficacy of some existing MM therapies, including immunomodulatory drugs and monoclonal antibodies (mAbs), rely on antibody-dependent cellular cytotoxicity (ADCC) driven by endogenous NK cells.4 Yet, MM often escapes NK cell–mediated cytotoxicity and induces NK cell dysfunction.4,5 

Dysfunctional NK cells are a hallmark of MM disease6,7 with decreased cytotoxicity, reduced activating receptors expression (NKG2D,8,9 CD16, 2B4,10 and DNAM-111), and increased inhibitory receptors expression (KIR2DL112 and PD113). To overcome this, the adoptive transfer of ex vivo–activated or engineered autologous or allogeneic NK cells has been used.14 Cytokine-induced memory-like NK cells are of particular interest because they acquire memory-like features, including elevated activity upon restimulation, greater interferon gamma (IFN-γ) responses, and antitumor cytotoxicity.15,16 Memory-like NK cells safely induce clinical remission in adult and pediatric patients with acute myeloid leukemia17-19 but have not been studied in MM.

The successful adoptive transfer of memory-like NK cells in MM depends on the ability of the cells to recognize malignant PCs. NKG2D, a major NK cell–activating receptor, is upregulated in memory-like NK cells compared with that in conventional NK cells (cNK cells).17 NKG2D recognizes malignant PCs via its ligands major histocompatibility complex class I chain–related protein A (MICA) and major histocompatibility complex class I chain–related protein B (MICB).11,20,21 MM cells evade recognition by proteolytically shedding surface MICA/B, which leads to NKG2D internalization.9 

The novel mAb 7C6 binds to the proteolytic MICA α3 domain, stabilizing surface MICA/B and enhancing NK cell–mediated ADCC.22,23 We hypothesized that 7C6 would simultaneously blockade of MICA/B shedding and increase ADCC by memory-like NK cells, resulting in enhanced MM cytotoxicity.

Previous studies have shown that surface MICA/B decreases and soluble MICA/B increases stepwise from monoclonal gammopathy of undetermined significance to overt MM due to shedding induced by the coordinated action of disulfide isomerase ERp5 and proteases, including a disintegrin and metalloproteinase domain-containing protein 10 (ADAM10) and ADAM17.9,24 We investigated the extent of surface MICA/B shedding in MM cells. MICA/B shedding involves ERp5 disulfide isomerase, which removes a disulfide bond from the α3 domain of MICA/B,25 allowing subsequent cleavage by proteases from the ADAM family.26 Therefore, MICA, MICB, ADAM10, and ERp5 transcript levels in human MM cell lines (HMCLs) were analyzed in the Cancer Cell Line Encyclopedia (CCLE) (21Q4).27 MICB and PDIA6 (encoding ERp5) transcription were the highest in HMCLs among the 1389 tumor cell lines (Figure 1A); this was not the case for MICA (Figure 1A). In public single-cell sequencing data from tumors of patients with newly diagnosed MM,28 PDIA6 transcription appeared higher in primary PCs from MM BM than in healthy control PCs (supplemental Figure 1A). Although MICA/B surface protein levels appeared to be significantly higher in primary PCs from patients with newly diagnosed multiple myeloma than in patients with monoclonal gammopathy of undetermined significance, it should be noted that it is difficult to maintain good viability of primary PCs ex vivo, which likely affects MICA/B shedding and thus surface MICA/B levels. Furthermore, the MICA and MICB transcript and protein levels in primary PCs were low, limiting our ability to interpret changes with disease progression (supplemental Figure 1B-C). Overall, prerequisites for proteolytic MICA/B cleavage were present in MM.

Figure 1.

Inhibition of MICA/B shedding enhances memory-like NK cell–mediated cytotoxicity against MM. (A) MICA, MICB, PDIA6 (ERp5), and ADAM10 gene expression levels in 1389 tumor cell lines from the Cancer Cell Line Encyclopedia database. Gene expression levels are presented as log2(transcripts per million [TPM] + 1). (B) MICA/B surface expression in 10 different HMCLs, as measured by flow cytometry. (C) MICA/B gene expression levels in 9 different HMCLs from the publicly available Broad CCLE data set. (D) Dose-dependent increase in MICA/B surface expression in HMCL MM.1S in response to 7C6 treatment. Significance was calculated using the t test. (E) Shed MICB was determined using sandwich enzyme-linked immunosorbent assay. A decrease in soluble MICB was seen in response to the 7C6 mAb treatment. For functional assays, effector cells (conventional or memory-like NK cells) were cocultured with MM.1S at a 5:1 effector-to-target ratio for 6 hours. (F) Increased CD107a degranulation (n = 7), IFN-γ production (n = 10), and tumor necrosis factor-α (TNF-α) production (n = 7) measured by flow cytometry in response to 7C6 treatment of cNK cells. (G) Summary of cytotoxicity (n = 9) of cNK and memory-like NK cells upon 7C6 mAb treatment of MM.1S cells. Data are depicted as the mean ± standard deviation. Significance was calculated using the Wilcoxon test for comparison of IC and 7C6 treatment conditions, and the Mann-Whitney U test was used for the remainder of the comparisons; ns, P > .05; ∗P ≤ .05; ∗∗P ≤ .01; ∗∗∗P ≤ .001; ∗∗∗∗P ≤ .0001. E:T, effector:target; IC, isotype control; ns, not significant.

Figure 1.

Inhibition of MICA/B shedding enhances memory-like NK cell–mediated cytotoxicity against MM. (A) MICA, MICB, PDIA6 (ERp5), and ADAM10 gene expression levels in 1389 tumor cell lines from the Cancer Cell Line Encyclopedia database. Gene expression levels are presented as log2(transcripts per million [TPM] + 1). (B) MICA/B surface expression in 10 different HMCLs, as measured by flow cytometry. (C) MICA/B gene expression levels in 9 different HMCLs from the publicly available Broad CCLE data set. (D) Dose-dependent increase in MICA/B surface expression in HMCL MM.1S in response to 7C6 treatment. Significance was calculated using the t test. (E) Shed MICB was determined using sandwich enzyme-linked immunosorbent assay. A decrease in soluble MICB was seen in response to the 7C6 mAb treatment. For functional assays, effector cells (conventional or memory-like NK cells) were cocultured with MM.1S at a 5:1 effector-to-target ratio for 6 hours. (F) Increased CD107a degranulation (n = 7), IFN-γ production (n = 10), and tumor necrosis factor-α (TNF-α) production (n = 7) measured by flow cytometry in response to 7C6 treatment of cNK cells. (G) Summary of cytotoxicity (n = 9) of cNK and memory-like NK cells upon 7C6 mAb treatment of MM.1S cells. Data are depicted as the mean ± standard deviation. Significance was calculated using the Wilcoxon test for comparison of IC and 7C6 treatment conditions, and the Mann-Whitney U test was used for the remainder of the comparisons; ns, P > .05; ∗P ≤ .05; ∗∗P ≤ .01; ∗∗∗P ≤ .001; ∗∗∗∗P ≤ .0001. E:T, effector:target; IC, isotype control; ns, not significant.

Close modal

If proteolytic shedding indeed occurs in MM, it would result in a discordance between MICA/B transcription and surface protein expression. We measured surface MICA/B levels by flow cytometry in 10 HMCLs and observed heterogeneous expression (Figure 1C). We next collected MICA/B transcript levels for these HMCLs from CCLE (Figure 1D). There was no correlation between transcript levels and protein expression (supplemental Figure 1D), suggesting that posttranslational processes determine MM MICA/B surface expression. This agrees with studies showing NKG2D ligand expression is finely regulated by transcriptional, posttranscriptional, and posttranslational mechanisms.29 Finally, to determine whether MICA and MICB transcript levels could be of prognostic value in MM, we analyzed the Multiple Myeloma Research Foundation CoMMpass data. High MICB (but not MICA) expression was significantly associated with inferior progression-free survival and overall survival (supplemental Figure 2A-B).

Antibodies targeting the MICA/B proteolytic cleavage site, the α3 domain,22 prevent shedding in solid malignancies. We treated HMCL with the novel anti-α3 domain mAb, 7C6. Flow cytometry of HMCL MM.1S revealed that 7C6 induced a dose-dependent increase in surface MICA/B (Figure 1D) and a corresponding decrease in shed MICA/B (Figure 1F). These results suggest that proteolytic cleavage of MICA/B occurs in MM and is effectively inhibited by 7C6 mAb. To assess the functional consequences of MICA/B stabilization, we examined NK cell activation in coculture assays of NK cells with KMS-12-BM, MM.1S, and RPMI-8226 pretreated with 7C6 or isotype control. We hypothesized that MICA/B stabilization would increase NK cell activation against HMCLs with sufficient MICA/B surface levels (MM.1S and RPMI-8226). Indeed, MICA/B stabilization led to increased NK cell activation, with increased cytotoxicity, degranulation (CD107a), and IFN-γ production against MM.1S (Figure 1G) and RPMI-8226 (supplemental Figure 3). This demonstrates that a MICA/B-stabilizing antibody can enhance human NK cell activation by MM cells.

The therapeutic potential of MICA/B stabilization in MM depends on the intrinsic functionality of NK cells in patients with MM. However, dysfunctional NK cells with reduced NKG2D expression have been described in MM8,9 and could hamper the efficacy of MICA/B stabilization, whereas adoptive transfer of memory-like NK cells could mitigate this limitation. We measured the reactivity of memory-like NK cells to MM target cells directly by coculturing memory-like NK cells generated from healthy human peripheral blood mononuclear cells with MM.1S. Memory-like NK cells exhibited enhanced cytotoxicity against MM target cells compared with cNK cells from the same donors (Figure 1I; supplemental Figure 4A). We repeated the coculture experiments by pretreating MM.1S with 7C6 or isotype control. Memory-like NK cells demonstrated significantly higher cytotoxicity against 7C6-pretreated MM.1S (Figure 1H; supplemental Figure 4B), likely due to increased MICA/B and NKG2D expression (supplemental Figure 5). Increased IFN-γ and tumor necrosis factor-α secretion was observed upon combining memory-like NK cells with 7C6 treatment, whereas CD107a degranulation was only significantly increased when compared with the control and 7C6 treatment (supplemental Figure 4C-F). These data demonstrated that human memory-like NK cells are superior to cNK cells in response to MM cells. MICA/B shedding further increased anti-MM memory-like NK cell activity.

To test the combined effect of MICA/B stabilization and memory-like NK cell adoptive transfer, we used a xenograft MM model17,22,30 (supplemental Methods; Figure 2A). Animal experiments were in accordance with and approved by the Institutional Animal Care and Use Committee. Luciferase-expressing MM.1S cells were engrafted into irradiated NOD scid gamma (NSG) mice, and memory-like NK cells were adoptively transferred on day 3. NK cells were isolated from 7 independent donors and then differentiated into memory-like NK cells by activating briefly (12-16 hours) with interleukin-12 (IL-12), IL-15, and IL-18. Memory-like NK cells were adoptively transferred into 7 independent cohorts of tumor-bearing NSG mice. IL-2 was injected intraperitoneal every other day for the duration of the study to support NK cell survival. Starting 1 day after NK cell transfer, 7C6-hIgG1 or isotype control antibodies were administered once per week, and tumor progression was tracked via bioluminescence imaging. Memory-like NK cells more effectively controlled early MM tumor cell outgrowth in mice treated with 7C6-hIgG1 mAb than in those treated with control antibodies (Figure 2B-C). However, early MM control did not result in improved survival of these mice (supplemental Figure 6A; P = .89). This was partly explained by donor-dependent responses to 7C6 therapy (supplemental Figure 6B-C). Of note, when the survival of “responding” donors was analyzed separately, a statistical difference in survival was detected (data not shown) but was limited by the small number of tested mice. Whether the donor-dependent responses result from MICA/B surface variation, in vivo persistence, exhaustion of adoptively transferred NK cells, or technical variance remains unknown.

Figure 2.

Combined MICA/B stabilization and memory-like NK cell transfer leads to early tumor control in vivo. (A) Experimental design. (B) Representative bioluminescence imaging (BLI) of recipient mice 9 and 12 days after engraftment with MM.1S-luc. Note the scale bar with lower intensities for the BLI signal in the 7C6-treated mice. (C) Enumeration of serial BLI measurements that show reduced tumor burden in mice receiving 7C6 mAb treatment compared with the isotype control. Differences were determined using the Mann-Whitney test; ns, P > .05; ∗P ≤ .05; ∗∗P ≤ .01; ∗∗∗P ≤ .001; ∗∗∗∗P ≤ .0001. IC, isotype control; ns, not significant.

Figure 2.

Combined MICA/B stabilization and memory-like NK cell transfer leads to early tumor control in vivo. (A) Experimental design. (B) Representative bioluminescence imaging (BLI) of recipient mice 9 and 12 days after engraftment with MM.1S-luc. Note the scale bar with lower intensities for the BLI signal in the 7C6-treated mice. (C) Enumeration of serial BLI measurements that show reduced tumor burden in mice receiving 7C6 mAb treatment compared with the isotype control. Differences were determined using the Mann-Whitney test; ns, P > .05; ∗P ≤ .05; ∗∗P ≤ .01; ∗∗∗P ≤ .001; ∗∗∗∗P ≤ .0001. IC, isotype control; ns, not significant.

Close modal

Here, we demonstrated that MICA/B shedding occurs in MM, and inhibition by the novel 7C6 mAb enhances memory-like NK cell–mediated responses against MM cells. Our data support the hypothesis that this combination of targeting MICA/B on MM cells and introducing memory-like NK cells results in improved NK cell–mediated immune surveillance in MM, making it a strong case for evaluation as a potential immunotherapeutic approach in MM.

The review boards of all participating centers approved the study, and the study was performed in accordance with the Declaration of Helsinki.

Acknowledgments: Anna V. Justis, a medical writer employed by the Dana-Farber Cancer Institute, provided writing and editorial support to the authors.

This work was supported by National Institutes of Health grant R01CA238039 (K.W.W.). R.R. is a member of the Parker Institute for Cancer Immunotherapy at the Dana-Farber Cancer Institute, and his work is supported, in part, by the Parker Institute for Cancer Immunotherapy. R.R. is a recipient of the Career Development Award from the Leukemia and Lymphoma Society. I.M.G. was supported by funding from the Dr. Miriam and Sheldon G. Adelson Medical Research Foundation. L.F.d.A. was funded by Cancer Research Institute award CRI3483, Leukemia and Lymphoma Society award 6647-23, US Department of Defense project number CA210940, the Tisch Cancer Institute Development Funds Award (National Cancer Institute grant P30CA196521), and the Elsa U. Pardee Foundation.

Contribution: S.T., K.W.W., R.R., and I.M.G. conceptualized the study; S.T., L.F.d.A., K.W.W., R.R., T.C., G.G., and I.M.G. developed the methodology; S.T., N.K.S., S.P., L.L., H.D., J.V.C., N.P., M.R., Y.S., L.B.F., A.C., J.-B.A., R.B., and E.L. conducted the investigations; S.T. wrote the original draft of the manuscript; K.W.W. and I.M.G. acquired funding; T.C., R.R., and I.M.G. supervised the study; and all authors reviewed and edited the final manuscript.

Conflict-of-interest disclosure: I.M.G. has acted as a consultant for Bristol Myers Squibb, Novartis, Amgen, Takeda, Noxxon Pharmaceuticals, Celgene, Sanofi, Genentech, GlaxoSmithKline, Gene Network Sciences Healthcare, Karyopharm Therapeutics, Adaptive Biotechnologies, Janssen, Medscape, and AbbVie, and has received honoraria from Celgene, Bristol Myers Squibb, Takeda, Amgen, Janssen, Karyopharm Therapeutics, Cellectar Biosciences, Adaptive Biotechnologies, Sanofi, and Medscape. P.S. reports consulting fees and research support from Celgene, Janssen, Karyopharm Therapeutics, and SkylineDx. K.W.W. serves on the scientific advisory boards of SQZ Biotech, Nextech Invest, Bisou Bioscience, and TScan Therapeutics; receives sponsored research funding from Novartis; and is a scientific cofounder of Immunitas Therapeutics. R.R. receives sponsored research funding from CRISPR Therapeutics and Skyline Therapeutics; owns equity in InnDura Therapeutics; and serves on the scientific advisory board of Glycostem Therapeutics and XNK Therapeutics. L.F.d.A. is a coinventor on an issued patent for an α3 domain-specific antibody and has served as a consultant for Cullinan Oncology. The remaining authors declare no competing financial interests.

Correspondence: Rizwan Romee, Division of Cellular Therapy and Stem Cell Transplantation, Dana-Farber Cancer Institute, 450 Brookline Ave, Boston, MA, 02215; email: rizwan_romee@dfci.harvard.edu; and Irene M. Ghobrial, Department of Medical Oncology, Dana-Farber Cancer Institute, 450 Brookline Ave, Boston, MA, 02215; email: irene_ghobrial@dfci.harvard.edu.

1.
Kumar
SK
,
Rajkumar
V
,
Kyle
RA
, et al
.
Multiple myeloma
.
Nat Rev Dis Primers
.
2017
;
3
:
17046
.
2.
Fonseca
R
,
Abouzaid
S
,
Bonafede
M
, et al
.
Trends in overall survival and costs of multiple myeloma, 2000-2014
.
Leukemia
.
2017
;
31
(
9
):
1915
-
1921
.
3.
Frohn
C
,
Hoppner
M
,
Schlenke
P
,
Kirchner
H
,
Koritke
P
,
Luhm
J
.
Anti-myeloma activity of natural killer lymphocytes
.
Br J Haematol
.
2002
;
119
(
3
):
660
-
664
.
4.
Clara
JA
,
Childs
RW
.
Harnessing natural killer cells for the treatment of multiple myeloma
.
Semin Oncol
.
2022
;
49
(
1
):
69
-
85
.
5.
Venglar
O
,
Bago
JR
,
Motais
B
,
Hajek
R
,
Jelinek
T
.
Natural killer cells in the malignant niche of multiple myeloma
.
Front Immunol
.
2021
;
12
:
816499
.
6.
Jurisic
V
,
Srdic
T
,
Konjevic
G
,
Markovic
O
,
Colovic
M
.
Clinical stage-depending decrease of NK cell activity in multiple myeloma patients
.
Med Oncol
.
2007
;
24
(
3
):
312
-
317
.
7.
Famularo
G
,
D'Ambrosio
A
,
Quintieri
F
, et al
.
Natural killer cell frequency and function in patients with monoclonal gammopathies
.
J Clin Lab Immunol
.
1992
;
37
(
3
):
99
-
109
.
8.
Costello
RT
,
Boehrer
A
,
Sanchez
C
, et al
.
Differential expression of natural killer cell activating receptors in blood versus bone marrow in patients with monoclonal gammopathy
.
Immunology
.
2013
;
139
(
3
):
338
-
341
.
9.
Jinushi
M
,
Vanneman
M
,
Munshi
NC
, et al
.
MHC class I chain-related protein A antibodies and shedding are associated with the progression of multiple myeloma
.
Proc Natl Acad Sci U S A
.
2008
;
105
(
4
):
1285
-
1290
.
10.
Fauriat
C
,
Mallet
F
,
Olive
D
,
Costello
RT
.
Impaired activating receptor expression pattern in natural killer cells from patients with multiple myeloma
.
Leukemia
.
2006
;
20
(
4
):
732
-
733
.
11.
El-Sherbiny
YM
,
Meade
JL
,
Holmes
TD
, et al
.
The requirement for DNAM-1, NKG2D, and NKp46 in the natural killer cell-mediated killing of myeloma cells
.
Cancer Res
.
2007
;
67
(
18
):
8444
-
8449
.
12.
Konjevic
G
,
Vuletic
A
,
Mirjacic Martinovic
K
,
Colovic
N
,
Colovic
M
,
Jurisic
V
.
Decreased CD161 activating and increased CD158a inhibitory receptor expression on NK cells underlies impaired NK cell cytotoxicity in patients with multiple myeloma
.
J Clin Pathol
.
2016
;
69
(
11
):
1009
-
1016
.
13.
Benson
DM
,
Bakan
CE
,
Mishra
A
, et al
.
The PD-1/PD-L1 axis modulates the natural killer cell versus multiple myeloma effect: a therapeutic target for CT-011, a novel monoclonal anti-PD-1 antibody
.
Blood
.
2010
;
116
(
13
):
2286
-
2294
.
14.
Shimasaki
N
,
Jain
A
,
Campana
D
.
NK cells for cancer immunotherapy
.
Nat Rev Drug Discov
.
2020
;
19
(
3
):
200
-
218
.
15.
O'Sullivan
TE
,
Sun
JC
,
Lanier
LL
.
Natural killer cell memory
.
Immunity
.
2015
;
43
(
4
):
634
-
645
.
16.
Romee
R
,
Schneider
SE
,
Leong
JW
, et al
.
Cytokine activation induces human memory-like NK cells
.
Blood
.
2012
;
120
(
24
):
4751
-
4760
.
17.
Romee
R
,
Rosario
M
,
Berrien-Elliott
MM
, et al
.
Cytokine-induced memory-like natural killer cells exhibit enhanced responses against myeloid leukemia
.
Sci Transl Med
.
2016
;
8
(
357
):
357ra123
.
18.
Bednarski
JJ
,
Zimmerman
C
,
Berrien-Elliott
MM
, et al
.
Donor memory-like NK cells persist and induce remissions in pediatric patients with relapsed AML after transplant
.
Blood
.
2022
;
139
(
11
):
1670
-
1683
.
19.
Shapiro
RM
,
Birch
GC
,
Hu
G
, et al
.
Expansion, persistence, and efficacy of donor memory-like NK cells infused for post-transplant relapse
.
J Clin Invest
.
2022
;
132
(
11
):
e154334
.
20.
Carbone
E
,
Neri
P
,
Mesuraca
M
, et al
.
HLA class I, NKG2D, and natural cytotoxicity receptors regulate multiple myeloma cell recognition by natural killer cells
.
Blood
.
2005
;
105
(
1
):
251
-
258
.
21.
Lu
X
,
Zhu
A
,
Cai
X
, et al
.
Role of NKG2D in cytokine-induced killer cells against multiple myeloma cells
.
Cancer Biol Ther
.
2012
;
13
(
8
):
623
-
629
.
22.
Ferrari de Andrade
L
,
Tay
RE
,
Pan
D
, et al
.
Antibody-mediated inhibition of MICA and MICB shedding promotes NK cell-driven tumor immunity
.
Science
.
2018
;
359
(
6383
):
1537
-
1542
.
23.
Badrinath
S
,
Dellacherie
MO
,
Li
A
, et al
.
A vaccine targeting resistant tumours by dual T cell plus NK cell attack
.
Nature
.
2022
;
606
(
7916
):
992
-
998
.
24.
Zingoni
A
,
Vulpis
E
,
Cecere
F
, et al
.
MICA-129 dimorphism and soluble MICA are associated with the progression of multiple myeloma
.
Front Immunol
.
2018
;
9
:
926
.
25.
Kaiser
BK
,
Yim
D
,
Chow
IT
, et al
.
Disulphide-isomerase-enabled shedding of tumour-associated NKG2D ligands
.
Nature
.
2007
;
447
(
7143
):
482
-
486
.
26.
Waldhauer
I
,
Goehlsdorf
D
,
Gieseke
F
, et al
.
Tumor-associated MICA is shed by ADAM proteases
.
Cancer Res
.
2008
;
68
(
15
):
6368
-
6376
.
27.
Barretina
J
,
Caponigro
G
,
Stransky
N
, et al
.
The Cancer Cell Line Encyclopedia enables predictive modelling of anticancer drug sensitivity
.
Nature
.
2012
;
483
(
7391
):
603
-
607
.
28.
de Jong
MME
,
Kellermayer
Z
,
Papazian
N
, et al
.
The multiple myeloma microenvironment is defined by an inflammatory stromal cell landscape
.
Nat Immunol
.
2021
;
22
(
6
):
769
-
780
.
29.
Raulet
DH
,
Gasser
S
,
Gowen
BG
,
Deng
W
,
Jung
H
.
Regulation of ligands for the NKG2D activating receptor
.
Annu Rev Immunol
.
2013
;
31
:
413
-
441
.
30.
Kurdi
AT
,
Glavey
SV
,
Bezman
NA
, et al
.
Antibody-dependent cellular phagocytosis by macrophages is a novel mechanism of action of elotuzumab
.
Mol Cancer Ther
.
2018
;
17
(
7
):
1454
-
1463
.

Author notes

S.T. and S.P. are joint first authors.

I.M.G. and R.R. are joint senior authors.

Only publicly available data sets have been used in this manuscript and referenced accordingly.

The full-text version of this article contains a data supplement.

Supplemental data