- 1The Clinical Medicine Postdoctoral Research Station, Department of Hematology, First Affiliated Hospital; Jinan University, Guangzhou, China
- 2Key Laboratory for Regenerative Medicine of Ministry of Education, Institute of Hematology, School of Medicine; Jinan University, Guangzhou, China
- 3Laboratory Center, Tianhe Nuoya Bio-Engineering Co. Ltd, Guangzhou, China
- 4Department of Clinical Laboratory, First Affiliated Hospital, Jinan University, Guangzhou, China Guangzhou, China
- 5Department of Hematology, First Affiliated Hospital, Guangzhou Medical University, Guangzhou, China, China
- 6Department of Hematology, Guangdong Provincial People’s Hospital, Guangdong Academy of Medical Sciences, Guangzhou, China
Despite the great success of immune-checkpoint inhibitor (ICI) treatment for multiple cancers, evidence for the clinical use of ICIs in acute myeloid leukemia (AML) remains inadequate. Further exploration of the causes of immune evasion in the bone marrow (BM) environment, the primary leukemia site, and peripheral blood (PB) and understanding how T cells are affected by AML induction chemotherapy or the influence of age may help to select patients who may benefit from ICI treatment. In this study, we comprehensively compared the distribution of PD-1 and TIGIT, two of the most well-studied IC proteins, in PB and BM T cells from AML patients at the stages of initial diagnosis, complete remission (CR), and relapse-refractory (R/R) disease after chemotherapy. Our results show that PD-1 was generally expressed higher in PB and BM T cells from de novo (DN) and R/R patients, while it was partially recovered in CR patients. The expression of TIGIT was increased in the BM of CD8+ T cells from DN and R/R patients, but it did not recover with CR. In addition, according to age correlation analysis, we found that elderly AML patients possess an even higher percentage of PD-1 and TIGIT single-positive CD8+ T cells in PB and BM, which indicate greater impairment of T cell function in elderly patients. In addition, we found that both DN and R/R patients accumulate a higher frequency of PD-1+ and TIGIT+ CD8+ T cells in BM than in corresponding PB, indicating that a more immunosuppressive microenvironment in leukemia BM may promote disease progression. Collectively, our study may help guide the combined use of anti-PD-1 and anti-TIGIT antibodies for treating elderly AML patients and pave the way for the exploration of strategies for reviving the immunosuppressive BM microenvironment to improve the survival of AML patients.
Introduction
In the last few decades, immunotherapy has emerged as the fourth pillar following surgery, radiation/chemotherapy and targeted therapy for solid tumor and leukemia therapy. One of the most effective immunotherapies includes resolving T cell dysfunction with immune checkpoint (IC) inhibitors (ICIs), such as programmed death receptor-1 (PD-1) and its ligand (PD-L1), cytotoxic T lymphocyte-associated antigen-4 (CTLA-4), and T cell immunoglobulin mucin-3 (Tim-3) (1–5). IC proteins are co-inhibitory receptors that distributed on the surface of several immune cells. After ligand binding, these regulators can transducing inhibitory signals to avoid immune injury under physiological conditions. However, under pathological conditions, such as chronic inflamation and cancer, continuous antigen stimulation and several immunosupressive factors could drive T cells to develop into a “unfitness” state called T cell exhaustion. One of the features of T cell exhaustion is the occurrence of elevated expession of several IC proteins. When these proteins bind their ligands, they reduce the anti-virus or anti-tumor effects of T cells, ultimately resulting in the immune escape of pathogens or tumors (3, 6–8).
Acute myeloid leukemia (AML) is a genetically, epigenetically, and clinically heterogeneous disease characterized by a failure in normal hematopoiesis and the accumulation of immature myeloid cells in the bone marrow (BM) and peripheral blood (PB) (9). For all AML patients (except for those with acute promyelocytic leukemia) who are medically able to tolerate chemotherapy, treatment with a 7-day infusion of cytarabine and a 3-day infusion of daunorubicin (“7+3”), has not changed in the past 50 years. Additionally, post-remission therapy aimed to prevent relapse of the disease. Additional cytotoxic chemotherapies (such as high or intermediate dose cytarabine), with or without targeted therapies, and allogenic hematopoietic stem cell transplantation (Allo-SCT) are the two commonly employed strategies. Several targeted therapies have been developed and approved for AML patients with special molecular and cytogenetic alterations, including venetoclax to target B-cell lymphoma 2, midostaurin, and gilteritinib to target FLT3, and ivosidenib and enasidenib to target mutant isocitrate dehydrogenase 1 and 2 (IDH1 and IDH2). The outcomes of AML patients who have a favorable or intermediate risk prognostic classification have improved with the addition of various targeted drugs; however, longer-term overall survival (OS) beyond 3–5 years remains low for adverse risk patients, particularly older AML patients (9–14). In recent studies, ICIs have demonstrated encouraging results when combined with hypomethylating agents (HMAs) in relapsed/refractory AML patients, particularly for those who were not exposed to HMA treatment, but most of the other clinical trials did not achieve a satisfactory response when using ICIs alone (15–19). One of the reasons for this discrepancy may be that the patients involved in those trials often experienced failure with multiple lines of conventional therapy, which could have a long-term impact on immunological fitness and clinical responses (20, 21). Therefore, comparing the immune status of T cells in AML patients before and after chemotherapy may help determine which type of patient could benefit from novel immune therapies. Additionally, cancer cells can orchestrate surrounding cells, such as regulatory T cells (Tregs), myeloid-derived suppressor cells (MDSCs), and plasmacytoid dendritic cells (pDCs), to construct an immunosuppressive tumor microenvironment (TME), enabling immune evasion and immunotherapy resistance (22, 23). In the BM, which is the nest of the leukemia progenitors generated and output in AML, the environment is quite complicated and resembles that of the TME in solid tumors to some degree. For example, the BM of AML patients also accumulates Treg, MDSC, and pDC cells that could inhibit the anti-leukemia immune response of T cells (24–27). Increasing evidence has shown that the BM immune environment of AML patients is profoundly altered, contributing to the severity of the disease; however, there have been limited studies comparing differences in T cell dysfunction between PB and BM. Thus, further elucidating the characteristics of the immune microenvironment in the BM of AML patients may help guide the use of immunotherapy drugs and facilitate the exploration of new immune targets.
PD-1, CTLA-4, and Tim-3 are the most well-studied IC proteins, and multiple other IC proteins are also targeted by immune checkpoint blockades or agonists in clinical research, such as lymphocyte activation Gene-3 (LAG-3), T cell immunoreceptor with immunoglobulin, and ITIM domain (TIGIT), and tumor necrosis factor receptor (OX40) (28–31). TIGIT, which is express on activated T cells, Treg cells, and NK cells, has been identified as a promising new target for cancer immunotherapy in recent years (2, 3, 7, 28, 32–34). Previously, high expression of PD-1 and TIGIT on PB T cells from AML patients has been reported (33, 35, 36), and higher PD-1 expression also detected on BM T cells in our previous study (37); however, how the expression of TIGIT and its co-expression with PD-1 changes in the BM of de novo (DN) AML patients remains unknown. Thus, we comprehensively compared the single and dual expression of PD-1 and TIGIT in both PB and BM T cells from AML patients at initial diagnosis. In addition, since ICI treatment typically used for patients with relapsed-refractory (R/R) disease, it is necessary to check the fitness of T cells from patients who have received chemotherapy. Thus, we also analyzed the expression of PD-1 and TIGIT in the PB and BM of AML patients at completed remission (CR) and disease relapse after chemotherapy.
Materials and Methods
Sample Collection and Preparation
A total of 38 PB samples and 32 BM samples from de novo AML patients [(median age, range), (PB: 54.5, 11-88; BM: 57.5, 14-80)] were collected, including 26 paired samples. A total of 17 PB samples and 21 BM samples from AML patients who achieved CR (PB: 46, 13-80; BM: 47, 13-73) were collected, including 16 paired samples. A total of 10 PB samples and 8 BM samples from AML patients with relapsed-refractory disease (PB: 51, 23-80; BM: 49, 23-67) were collected, including 6 paired samples. All of the patients included in our CR and R/R cohorts are transplant-naïve patients. In addition, 36 PB samples from healthy donors (53, 13-65), and 14 unpaired BM samples from hematopoietic stem cell transplantation donors or iron deficiency anemia patients (47, 17-68) underwent bone marrow aspiration served as the control group. Sample details shown in Table 1. All samples collected with informed consent. Ethical approval for the study was obtained from the Ethics Committee of the First Affiliated Hospital of Jinan University. All procedures were conducted according to the guidelines of the Medical Ethics Committees of the Health Bureau of the Guangdong Province in China. All samples were freshly obtained and subjected to immediate preparation.
Flow Cytometry Analysis
Cell surface staining for flow cytometry was performed using the following antibodies: CD45-BUV395 (clone HI30, BD), CD3-AF700 (clone UCHT1, BD), CD4-APC-H7 (clone RPA-T4, BD), CD8-APC-H7 (clone SK1, BD), TIGIT-PE/Dazzle™ 594 (clone A15153G, Biolegend), TIGIT-BV421 (clone A15153G, Biolegend), and PD1-BV421 (clone EH12.2H7, Biolegend). These antibodies were used to analyze surface receptors in two different panels. Isotype-matched antibodies, labeled with the proper fluorochromes, were used as negative controls. Extracellular staining was performed according to the manufacturer’s instructions. Briefly, an antibody cocktail was added to whole blood or BM samples, which were then incubated at room temperature for 15 minutes in the dark, followed by red cell lysis with lysis buffer (BD; Cat: 555899). The lysed cells were washed twice with phosphate buffer saline (PBS), and then 350 μl stain buffer was added for further flow cytometer analysis. Twenty microliters of absolute count microsphere (Thermo; Cat: C36950) were added to the samples for absolute number analysis. Cells were analyzed using a BD Fortessa flow cytometer (BD Biosciences), and data analysis performed with Flowjo 10.6 software.
Statistical Analysis
Data were analyzed by correlation, linear regression, independent t-test, or paired t-test depending on the experimental design. Calculations were performed using GraphPad Prism, version 8.02 software. Significance is indicated as *p < 0.05, **p < 0.01, ***p < 0.001, ****p < 0.0001. Values of p < 0.05 were considered significant.
Results
The Distribution of PD-1 and TIGIT in the PB and BM of AML Patients During Diagnosis, Complete Remission, and in Relapsed-Refractory Patients
We first compared the PD-1 and TIGIT expression frequency on CptD4+ and CD8+ T cells from the PB and BM of AML patients in the DN, CR, and R/R stages. The results demonstrated that PD-1 expression increased on both PB CD4+ (20.63% vs 15.86%, p=0.029) and CD8+ (24.77% vs 17.49%, p=0.011) T cells in DN patients when compared with healthy individuls (HIs), while the same result was found when comparing the BM CD4+ (26.17% vs 18.48%, p=0.010) and CD8+ (41.17% vs 21.68%, p=0.0003) T cell populations. When patients achieved CR, we found that PD-1 expression recovered in most T cell populations except for the PB CD8+ population (27.00% vs 17.49%, p=0.014). In addition, we also detected higher PD-1 expression on PB CD4+ (23.22% vs 15.86%, p=0.023) and CD8+ (26.84% vs 17.49%, p=0.021) T cells as well as BM CD8+ (37.14% vs 21.68%, p=0.018) T cells from R/R patients when compared with HIs (Figure 1A). Unlike the high expression of PD-1 on both PB and BM T cells in DN AML patients, TIGIT expression was increased only in the BM CD8+ (59.18% vs 38.69%, p=0.003) T cell population. Interestingly, the expression of TIGIT was also higher in the BM CD8+ population from CR (55.08% vs 38.69%, p=0.002) and R/R (65.94% vs 38.69%, p<0.0001) patients compared with HIs (Figure 1B). The above results indicated that higher expression of PD-1 in PB and BM T cells is a characteristic of patients with DN AML and might be a sign of disease relapse, while PD-1 expression recovery might be a favorable signal of disease remisson. Whether PD-1 expression could predict the prognosis or survival of AML patients needs to be confirmed in a future study.
Figure 1 T cells from the PB and BM of DN AML and R/R AML patients generally have high expression of PD-1, and BM CD8+ T cells from DN, CR, and R/R patients have increased TIGIT expression. (A) The top figure shows the gating strategy for PD-1 in the CD4+ and CD8+ populations by flow cytometry. The frequency of PD-1 by subset in the PB of HIs (CD4, n=33; CD8, n=36) and patients in the AML-DN (CD4, n=34; CD8, n=38), AML-CR (n=17), and AML-R/R (n=10) cohorts (Upper); The frequency of PD-1 by subset in the BM of HIs (n=14), and patients in the AML-DN (CD4, n=29; CD8, n=32), AML-CR (n=21), and AML-R/R (CD4, n=7; CD8, n=8) cohorts (Lower). The PD-1 expression frequency generally increased in T cells from DN and R/R patients. (B) The top figure shows the gating strategy for TIGIT in the CD4+ and CD8+ populations by flow cytometry. The frequency of TIGIT by subset in the PB (Upper) and BM of HIs, and patients in the AML-DN, AML-CR, and AML-R/R cohorts (Lower); The TIGIT expression frequency increased in BM CD8+ T cells from DN, CR, and R/R patients compared with HIs. The p values shown are from independent t-tests between groups. *p < 0.05, **p < 0.01, ***p < 0.001, ns denotes not significant.
Age Correlation of PD-1 and TIGIT Expression on PB and BM T Cells
TIGIT was previously reported to have increased expression in blood T cells of elderly healthy individuals (38); thus, we analyzed the correlation between age and the expression frequency of PD-1 and TIGIT on PB and BM CD4+ and CD8+ T cells in each group. We found that PD-1 had no age correlation with any T cell subset in the HI group; however, in the DN group, the expression of PD-1 on PB and BM CD8+ T cells was positively correlated with age. In addition, the expression of PD-1 on BM CD4+ T cells also demonstrated a trend toward a positive correlation (Figure 2A). These results suggest that the expression of PD-1 in elderly AML patients may be even higher than that in young patients, which may weaken the anti-leukemia T cell response of elderly patients. Elderly AML generally defined as AML in a patient who is greater than 60 years of age. Hence, we further compared the expression of PD-1 on CD4+ and CD8+ T cells in the PB and BM of young (< 60 years) and elderly (>= 60 years) AML patients. The results demonstrated that PD-1 is expressed higher on all T cell subsets except for the CD4+ subset in the PB of elderly AML patients (Figure 2B). As shown in Figure 2A, there was also no correlation between age and PD-1 expression for any of the T cell subsets in the CR and R/R groups.
Figure 2 PD-1+ and TIGIT+ CD8+ T cells increase in the BM of older DN-AML patients. (A) The upper panel shows the correlation between the frequency of PD-1 (red dots) or TIGIT (blue triangles) and age in the PB CD4+ and CD8+ T cell populations in HIs (CD4, n=33; CD8, n=36) and patients in the AML-DN (CD4, n=34; CD8, n=38), AML-CR (n=17), and AML-R/R (n=10) cohorts; the lower panel shows the correlation between the frequency of PD-1 (the red dots) or TIGIT (blue triangles) and age in the BM CD4+ and CD8+ T cell populations in HIs (n=14) and patients in the AML-DN (CD4, n=29; CD8, n=32), AML-CR (n=21), and AML-R/R (CD4, n=7; CD8, n=8) cohorts. Spearman’s nonparametric test used to test for correlations between each group. (B) Comparison of the expression of PD-1 on CD4+ and CD8+ T cells from the PB of DNy (CD4, n=19; CD8, n=21) and DNo (CD4, n=15; CD8, n=17) (younger and older than 60 years) patients (Left) from the BM of DNy (CD4, n=16; CD8, n=17) and DNo (CD4, n=13; CD8, n=15) patients (Right). (C) Comparison of the expression of TIGIT on CD4+ and CD8+ T cell subsets from the PB and BM of DNy and DNo AML patients. The sample cases used for TIGIT are the same as those used for the PD-1 comparison. The p values shown in B and C are from the independent student’s t-test. *p < 0.05, ns denotes not significant.
Regarding the expression of TIGIT, we found the exact age correlation for its expression on CD4+ and CD8+ T cells in HI PB, which was dramatically evident in the CD8+ population. Additionally, we found that the positive correlation between TIGIT and age in the PB CD4+ T cell population lost in the DN, CR, and R/R groups. In contrast, in the PB CD8+ T cell population, it was diminished in the DN group but recovered in the CR group. Considering the strong age correlation of TIGIT in the PB of HIs, we further dissected the HI, DN, and CR groups into young and elderly cohorts and compared the expression of TIGIT on each subset. The results demonstrated that in the cohorts under age 60, the frequency of TIGIT on the CD4+ T cells of DN AML patients was higher than that on HIs (26.31% vs 20.14%, p=0.038); however, there was no difference in other comparisons between AML and HIs in the same age cohorts (Supplementary Figure 1). Regarding the expression of TIGIT on HI BM T cells, there was no age correlation found in either the CD4+ or CD8+ population, but a significant age correlation was found for the DN AML group (Figure 2A). Further analysis revealed that the frequency of TIGIT in the CD8+ population in BM was higher in elderly vs. younger DN AML patients (Figure 2C).
In summary, we found that the frequency of PD-1 expression on T cells in HIs had no age correlation in PB or BM; however, strikingly, there was higher PD-1 expression in elderly AML patients. Moreover, the above results confirmed that the expression of TIGIT on T cells is closely correlated with age in PB but not BM. However, higher TIGIT expression is still detected in the CD4+ population of DN patients under 60 years. In addition, relatively higher TIGIT expression was also detected in the BM CD8+ T cell subset of elderly DN patients.
Higher PD-1 and TIGIT Expression Detected in the BM of DN and R/R AML Patients
Previous studies have found that the leukemia BM micro- environment possesses several immune suppressive factors that may protect malignant hematopoietic stem cells from immunological surveillance. We and others have also reported that PD-1 increased in BM CD8+ T cells from DN AML patients, which may promote immune evasion for leukemia cells in the BM (36, 37, 39, 40). In this study, we further compared the expression of PD-1 and TIGIT in paired PB and BM from DN, CR, and R/R patients. The results demonstrated that both the CD4+ and CD8+ T cell subsets in the BM of DN AML patients possessed a higher percentage of PD-1 (Figure 3A) and TIGIT (Figure 3B) when compared with corresponding PB, and the increase in the CD8+ population was most apparent. In addition, the above differences did not exist for CR patients and remained in the CD8+ population for R/R patients. These results indicate that the BM microenviroment of AML patients possesses an increased number of immunosupressive factors compared with PB, and reviving the immunosupressive microenviroment might be a key factor for treating AML and preventing relapse.
Figure 3 BM CD4+ and CD8+ T cells from DN-AML and R/R-AML patients have a higher percentage of PD-1- and TIGIT-expressing T cells than matched PB. (A) The flow-cytometry analysis detected an increase in the frequency of PD-1-expressing CD4+ and CD8+ T cells in the BM of DN-AML patients compared with matched PB (CD4, n=24; CD8, n=27); the expression of PD-1 on CD4+ and CD8+ T cells between PB and BM was not different in the CR-AML cohort (CD4, n=15; CD8, n=16). Increased PD-1 expression was also detected in BM CD8+ T cells from R/R patients (n=6). (B) Flow-cytometry analysis detected an increased frequency of TIGIT-expressing CD4+ and CD8+ T cells in the BM of DN-AML patients compared with matched PB (CD4, n=24; CD8, n=27); the expression of TIGIT on CD4+ and CD8+ T cells between PB and BM was not different in the CR-AML cohort (CD4, n=15; CD8, n=16); increased TIGIT expression was also detected in the BM CD8+ T cells of R/R patients (n=6). The p values in (A, B) were calculated by the paired student’s t-test. *p < 0.05, **p < 0.01, ***p < 0.001, ns denotes not significant.
Higher Co-Expression of PD-1 and TIGIT was Detected in the DN and R/R AML Cohorts, and DN-AML Patients Accumulated More PD-1+TIGIT+ CD8+ T Cells in BM Than PB
Co-expression of several IC proteins is a characteristic of exhausted T cells. A high percentage of PD-1+TIGIT+CD8+ T cells was found in the PB of DN AML patients with the exhaustion phenotype; however, it remains unknown whether they change in the BM and in patients who achieve CR or those with R/R disease. Thus, we further examined the co-expression of PD-1 and TIGIT in the PB and BM of each cohort. As shown in Figure 4A, higher co-expression of PD-1 and TIGIT only detected on PB CD4+ T cells from DN patients, while, as shown in Figure 4B, higher co-expression detected on PB CD8+ T cells from DN and R/R AML patients but not on the cells in the CR group. However, increased co-expression was also found in the BM CD8+ T cell population of AML patients in all disease conditions. However, statistical significance was found only for DN patients based on the limited CR and R/R patient samples analyzed (Figure 4B). We further compared the co-expression of PD-1 and TIGIT between PB and BM in AML patients at different disease stages. We found no significant co-expression difference in the CD4+ population in any of the AML sub-groups (Figure 4C, left). In terms of the CD8+ population, higher co-expression was detected in the DN and CR cohorts (Figure 4C, right). These results further suggest a more immunosupressive BM microenvironment in AML patients at initial diagnosis and relapse.
Figure 4 Co-expression of PD-1 and TIGIT is increased in the PB and BM T cells of DN-AML patients. (A) Shown is the co-expression of PD-1 and TIGIT in PB and BM CD4+ T cells from HIs (PB, n=12; BM, n=11) and patients in the DN-AML (PB, n=16; BM, n=15), CR-AML (PB, n=11; BM, n=14), and R/R-AML (n=5) cohorts. (B) Co-expression of PD-1 and TIGIT in PB and BM CD8+ T cells from HIs and patients in the DN-AML, CR-AML, and R/R-AML cohorts. The cases in the figure are the same as those in the CD4+ population. The p values in (A, B) were calculated by the independent student’s t-test. (C) Comparison of PD-1 and TIGIT co-expression between the paired PB and BM CD4+ and CD8+ subsets. Higher co-expression existed in the BM CD8+ T cell from most of the DN-AML and CR-AML patients. The p values shown in C were calculated by the paired student’s t-test. *p < 0.05, ns denotes not significant.
Discussion
Previous studies have reported the expression of PD-1 and TIGIT in AML patients at diagnosis and after hematopoietic stem cell transplantation (33, 35–37, 41, 42); however, most of these studies were only concerned about T cells in the PB. Some of these studies ignored the impact of patient age; thus, controversial results exist. In this study, we comprehensively compared the distribution of PD-1 and TIGIT, two of the most well-studied IC proteins in PB and BM T cells from AML patients at different disease stages. The results demonstrated that PD-1 is highly expressed in almost all CD4+ and CD8+ T cells in the PB and BM of DN AML patients and demonstrates higher expression in the elderly than younger patients. Expression in the BM was higher than that in PB. When patients achieve CR after induction chemotherapy, most of T cells can restore the average expression level of PD-1 except for the CD8+ T cells in PB. In addition, patients who were not response to or relapsed from chemotherapy maintained a high expression of PD-1 in most of the T cell subsets except for CD4+ T cells in BM. These results indicated that the function of T cells in elderly patients at first diagnosis might be worse than that in younger patients, which may be a reason for their chemotherapy intolerance, but this may be an advantage for receiving ICI therapy. Although most of the clinical trials using PD-1/PD-L1 inhibitors have been administered to R/R AML patients who experienced intensive induction chemotherapy at least twice (16, 32, 43, 44), a multi-center phase 2 study (NCT02845297) of pembrolizumab (PD-1–blocking antibody) and azacitidine (AZA) in patients with R/R AML and a small cohort of newly diagnosed AML patients (median age and range: 75, 67-83 years) suggested a favorable overall response rate (ORR; 71%) and survival for newly diagnosed and unfit older AML patients.
Moreover, this study also found that the AZA/pembro combination is safe, feasible, and well tolerated by these elderly patients (45). Moreover, a higher survival rate when using drugs targeting PD-1 in metastatic melanoma patients between the ages of 70 and 80 was reported in a large population-based Danish patient cohort (1,082 cases) (46). In addition to age, cytogenetic abnormalities are another critical prognostic factor for AML. Thus, analyzing changes in T cell status in patients with different karyotypes and gene mutations is also essential. Previously Williams et al. reported that PD-L1 is higher expressed in BM blasts from patients with TP53 mutation and also have a higher expression trend in patients with adverse cytogenetics; however, they did not find any correlation between cytogenetic abnormalities and the expression of IC proteins in T cells (36). We also analyzed the expression of PD-1 and TIGIT in patients with or without FLT3-ITD mutation and complex karyotype. However, there was no difference to be found in either group (data not shown). To explore the relationship of PD1 and TIGIT expression on T cells with other gene mutations, such as TP-53 and IDH1/2, would require a larger cohort study in the future.
The high expression of PD-1/PD-L1 is an essential indicator for effective ICI therapy; however, whether sufficient tumor-infiltrating T cells in the tumor microenvironment is also a prerequisite (32). We also analyzed the absolute numbers of T cells in the PB and BM samples included in this study. We found increased numbers of infiltrating T cells in the BM of DN patients, and CR and R/R patients had a lower T cell count in the PB and comparable numbers of T cells in the BM compared with HIs (Supplementary Figure 2). These results suggest that the T cells were activated and responded to leukemia antigen stimulation in DN patients. However, most of them may be functionally impaired by expressing IC proteins and other immunosuppressive factors. When receiving induction chemotherapy, most of the T cells were killed simultaneously by the unspecific cytotoxicity of the chemotherapy drugs; thus, there might not be sufficient numbers of T cells to respond to ICI treatment. This reason may be one of the causes for anti-PD-1 treatment without satisfactory effects in R/R patients; however, whether a single anti-PD-1 agent would affect the overall survival (OS) of untreated younger AML patients remains unknown. Recently, a phase 2 study (NCT02464657) using nivolumab combined with induction cytarabine plus idarubicin produced an ORR of 80% and a median OS of 18.54 months in untreated, younger AML patients. Interestingly, responding patients who continued idarubicin, cytarabine, and nivolumab beyond remission had similar overall survival and event-free survival compared with those bridged to allogeneic stem cell transplantation, suggesting the possibility of nivolumab efficacy in restoring host antitumor immune surveillance. Moreover, the authors also found that a higher prevalence of CD3-positive T cells in the bone marrow before the therapy administration appeared to predict the response (47). This study, combined with our findings here, may shed light on using anti-PD-1 as a single agent or in combination with other drugs for treating DN AML patients, particularly elderly DN AML patients. Previously, two studies also compared the absolute number of T cells in AML patients; one used flow-cytometry and found a higher CD3+ T cell number in the PB of DN AML patients compared with age-matched HIs. The different results reported here may be because we used anti-CD45 to exclude the impact of many leukemia cells in the samples (48). In another study, the authors used immunochemistry to compare the CD3+ T cell number in the BM and found a comparable amount of CD3+ T cells between DN patients and HIs (36). Considering the high sensitivity of flow cytometry in minimal residual disease detection and the greater number of samples analyzed in our study (31 vs 14), we believe that our results have higher accuracy. Whether the increased absolute count of T cells in the PB and BM of DN AML patients would predict the outcome and survival of DN AML patients is worth further exploration with a large cohort.
With regards to the expression of TIGIT, we confirmed Yangzi Song’s result that the TIGIT+ cell frequencies among PB CD8+ T cells were significant age correlated, whereas CD4+ T cells exhibited exhibited a weak correlation (38). Based on this age-induced change in PB T cells, the different expression of TIGIT in PB T cells between AMLy and AMLo patients is most likely to be a distinction that also exists in the HIy and HIo groups. In 2016, Yaxian Kong had reported that TIGIT expression elevated in the PB CD8+ T cells of AML patients, and high TIGIT was associated with poor clinical outcome in AML (49); however, we only found higher expression of TIGIT on PB CD4+ T cells of DN AML patients compared with HIs in the younger age group and on PB CD8+ T cells in both young and older patients. This inconsistency may arise from the average age of their healthy cohort appears to be younger than that of the AML cohort (51 vs. 60). Interestingly, the expression of TIGIT on BM T cells demonstrated no age correlation for HIs. Thus, we speculate that TIGIT function in BM T cells might be different from that in PB. Thus, the elevated TIGIT expression on the BM CD8+ T cells of DN, CR, and R/R patients may be mainly impacted by the leukemia microenvironment itself, while the BM microenvironment of elderly patients may further exacerbate the high expression of TIGIT. Although the immunosupressive functions of TIGIT have been demonstrated many times, whether its high expression on the BM CD8+ T cells of AML patients correlates with prognosis and survival remains to be seen in the future. In addition, further studies should elucidate the functional heterogeneity of TIGIT+ T cells co-expressed with other IC proteins and with its competitive activation receptor CD226, considering the complex mechanisms of TIGIT in suppressing T cell function.
Co-expression of several IC proteins is usually thought to be one of the characteristics of T cell exhaustion. Two previous studies have reported higher dual expression of PD-1 and TIGIT in PB T cells from DN AML patients, and functional experiments supported that PD-1+TIGIT+ CD8+ T cells were functionally exhausted (35, 50). In this study, we further found that the percentage of PD-1 and TIGIT double-positive T cells not only increased in PB but also in BM CD8+ T cells from DN AML patients. In addition, it was shown that there was an increase in PB and a trend toward an increase in BM CD8+ T cells in R/R AML patients. All these results indicated that a more immunosuppressive microenvironment might exist in the BM of AML patients. Previously Jia’s study found that intracellular expression of Granzyme B in BM CD8+ T cells was significantly lower compared with that of PB, which supports our results (51); however, a study from Lamble et al. reported that there was no significant proliferation difference between BM and PB CD8+ T cells from AML patients (52), which may be due to the high patient-to-patient variability when looking into their data. Such individual differences also existed in our study; thus, to precisely understand the dysfunctional status of BM T cells in each AML patient, a comprehensive assessment of the immune checkpoints expression, cytokine secretion, and proliferation capacity of paired PB and BM samples is needed. In addition, whether these results support the combined use of anti-PD-1 and anti-TIGIT ICIs for treating DN and R/R patients remains to be seen.
Taken together, in addition to the previously reported high expression of PD-1 and co-expression of PD-1 and TIGIT in the PB of DN AML patients and increased PD-1 expression in the BM of DN AML patients, we further identified the accumulation of TIGIT+ CD8+ T cells in the BM of DN, CR, and R/R patients. Moreover, we found a higher percentage of PD-1 and TIGIT co-expressing CD8+ T cells in the BM of DN AML patients. These results support the idea that the BM microenvironment of AML patients possesses many more immunosupressive factors than PB. Future treatment strategies focused on reviving the immunosuppressive BM microenvironment may improve AML patients’ survival. Moreover, we also found that increased PD-1 and TIGIT single-positive T cells exist in the PB and BM of older AML patients compared with younger patients, suggesting that T cell function in elderly patients might be even worse than that in younger patients. This finding may partially explain the chemotherapy intolerance of older AML patients, which brings hope for using ICIs to treat DN elderly AML patients who may have sufficient T cells to respond.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics Statement
The studies involving human participants were reviewed and approved by Ethics Committee of the First Affiliated Hospital of Jinan University. Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin.
Author Contributions
YQL and XH contributed to the concept development and study design. LX, LL, DLY, XBZ, YKZ, and SHS performed the laboratory studies. JL, ZYJ, XFZ, RHZ, YHL, MML, and XH collected the clinical samples. LL collected the clinical information of patients. LX and YQL drafted the manuscript. SHC managed the laboratory reagents and financial affairs. LX and YQL helped modify the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This study was supported by grants from the National Natural Science Foundation of China (Nos. 82000108, 82070152, 91642111, 81890991, and 81800143), China postdoctoral science foundation (2018M640884), Guangdong Provincial Outstanding Young Medical Talents Supporting Research Foundation (KJ012019459), Guangdong Provincial Natural Science Foundation (2020A1515110310), Guangzhou Science and Technology Project (201807010004 and 201803040017), and Key Laboratory for Regenerative Medicine of Ministry of Education Project (ZSYXM202001).
Conflict of Interest
Author YKZ was employed by company Tianhe Nuoya Bio-Engineerinig Co. Ltd.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We would like to acknowledge the Institute of Tumour Pharmacology, College of Pharmacology, College of Pharmacy of Jinan University for providing the flow cytometer for our project. We also would like to thank the healthy volunteers who donated blood or bone marrow for this project.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2021.686156/full#supplementary-material
Supplementary Figure 1 | Expression frequency of TIGIT on the PB CD4+ and CD8+ subsets of the HI, DN-AML, and CR-AML cohorts in the young and old age groups. Flow-cytometry analysis shows that the expression frequency of TIGIT increased in the PB CD4+ population of DNy patients when compared with HIy.
Supplementary Figure 2 | The absolute number of CD3+, CD4+, and CD8+ T cells in the PB and BM of the HI and AML cohorts. The total CD3+ T cell number in PB heavily decreased in CR and R/R patients who received induction chemotherapy, and total CD3+ T cells in the BM increased in DN-AML patients when compared with HIs.
References
1. Gong J, Chehrazi-Raffle A, Reddi S, Salgia R. Development of PD-1 and PD-L1 Inhibitors as a Form of Cancer Immunotherapy: A Comprehensive Review of Registration Trials and Future Considerations. J Immunother Cancer (2018) 6(1):8. doi: 10.1186/s40425-018-0316-z
2. Rotte A. Combination of CTLA-4 and PD-1 Blockers for Treatment of Cancer. J Exp Clin Cancer Res (2019) 38(1):255. doi: 10.1186/s13046-019-1259-z
3. Li X, Shao C, Shi Y, Han W. Lessons Learned From the Blockade of Immune Checkpoints in Cancer Immunotherapy. J Hematol Oncol (2018) 11(1):31. doi: 10.1186/s13045-018-0578-4
4. Wang JC, Chen C, Kundra A, Kodali S, Pandey A, Wong C, et al. Programmed Cell Death Receptor (PD-1) Ligand (PD-L1) Expression in Philadelphia Chromosome-Negative Myeloproliferative Neoplasms. Leuk Res (2019) 79:52–9. doi: 10.1016/j.leukres.2019.02.010
5. Kohler N, Ruess DA, Kesselring R, Zeiser R. The Role of Immune Checkpoint Molecules for Relapse After Allogeneic Hematopoietic Cell Transplantation. Front Immunol (2021) 12:634435. doi: 10.3389/fimmu.2021.634435
6. Kim N, Kim HS. Targeting Checkpoint Receptors and Molecules for Therapeutic Modulation of Natural Killer Cells. Front Immunol (2018) 9:2041. doi: 10.3389/fimmu.2018.02041
7. Yu X, Gao R, Li Y, Zeng C. Regulation of PD-1 in T Cells for Cancer Immunotherapy. Eur J Pharmacol (2020) 881:173240. doi: 10.1016/j.ejphar.2020.173240
8. Tan J, Chen S, Lu Y, Yao D, Xu L, Zhang Y, et al. Higher PD-1 Expression Concurrent With Exhausted CD8+ T Cells in Patients With De Novo Acute Myeloid Leukemia. Chin J Cancer Res (2017) 29(5):463–70. doi: 10.1016/j.exphem.2017.06.190
9. Khanal N, Upadhyay Banskota S, Bhatt VR. Novel Treatment Paradigms in Acute Myeloid Leukemia. Clin Pharmacol Ther (2020) 108(3):506–14. doi: 10.1002/cpt.1962
10. Xuan L, Wang Y, Huang F, Fan Z, Xu Y, Sun J, et al. Sorafenib Maintenance in Patients With FLT3-ITD Acute Myeloid Leukaemia Undergoing Allogeneic Haematopoietic Stem-Cell Transplantation: An Open-Label, Multicentre, Randomised Phase 3 Trial. Lancet Oncol (2020) 21(9):1201–12. doi: 10.1016/S1470-2045(20)30455-1
11. Cerrano M, Itzykson R. New Treatment Options for Acute Myeloid Leukemia in 2019. Curr Oncol Rep (2019) 21(2):16. doi: 10.1007/s11912-019-0764-8
12. Kayser S, Levis MJ. Advances in Targeted Therapy for Acute Myeloid Leukaemia. Br J Haematol (2018) 180(4):484–500. doi: 10.1111/bjh.15032
13. Perl AE. The Role of Targeted Therapy in the Management of Patients With AML. Blood Adv (2017) 1(24):2281–94. doi: 10.1182/bloodadvances.2017009829
14. Bazinet A, Assouline S. A Review of FDA-Approved Acute Myeloid Leukemia Therapies Beyond '7+3'. Expert Rev Hematol (2021) 14(2):185–97. doi: 10.1080/17474086.2021.1875814
15. Giannopoulos K. Targeting Immune Signaling Checkpoints in Acute Myeloid Leukemia. J Clin Med (2019) 8(2):236. doi: 10.3390/jcm8020236
16. Daver N, Garcia-Manero G, Basu S, Boddu PC, Alfayez M, Cortes JE, et al. Efficacy, Safety, and Biomarkers of Response to Azacitidine and Nivolumab in Relapsed/Refractory Acute Myeloid Leukemia: A Nonrandomized, Open-Label, Phase II Study. Cancer Discov (2019) 9(3):370–83. doi: 10.1158/2159-8290.CD-18-0774
17. Stahl M, Goldberg AD. Immune Checkpoint Inhibitors in Acute Myeloid Leukemia: Novel Combinations and Therapeutic Targets. Curr Oncol Rep (2019) 21(4):37. doi: 10.1007/s11912-019-0781-7
18. Przespolewski A, Szeles A, Wang ES. Advances in Immunotherapy for Acute Myeloid Leukemia. Future Oncol (2018) 14(10):963–78. doi: 10.2217/fon-2017-0459
19. Gale RP. Will Immune Therapy Cure Acute Myeloid Leukemia? Blood Sci (2019) 1: (1):2–3. doi: 10.1097/BS9.0000000000000024
20. Davis JE, Handunnetti SM, Ludford-Menting M, Sharpe C, Blombery P, Anderson MA, et al. Immune Recovery in Patients With Mantle Cell Lymphoma Receiving Long-Term Ibrutinib and Venetoclax Combination Therapy. Blood Adv (2020) 4(19):4849–59. doi: 10.1182/bloodadvances.2020002810
21. Cooke RE, Quinn KM, Quach H, Harrison S, Prince HM, Koldej R, et al. Conventional Treatment for Multiple Myeloma Drives Premature Aging Phenotypes and Metabolic Dysfunction in T Cells. Front Immunol (2020) 11:2153. doi: 10.3389/fimmu.2020.02153
22. Hinshaw DC, Shevde LA. The Tumor Microenvironment Innately Modulates Cancer Progression. Cancer Res (2019) 79(18):4557–66. doi: 10.1158/0008-5472.CAN-18-3962
23. Anderson NM, Simon MC. The Tumor Microenvironment. Curr Biol (2020) 30(16):R921–5. doi: 10.1016/j.cub.2020.06.081
24. Mendez-Ferrer S, Bonnet D, Steensma DP, Hasserjian RP, Ghobrial IM, Gribben JG, et al. Bone Marrow Niches in Haematological Malignancies. Nat Rev Cancer (2020) 20(5):285–98. doi: 10.1038/s41568-020-0245-2
25. Swatler J, Turos-Korgul L, Kozlowska E, Piwocka K. Immunosuppressive Cell Subsets and Factors in Myeloid Leukemias. Cancers (Basel) (2021) 13(6):1203. doi: 10.3390/cancers13061203
26. Mu LL, Ke F, Guo XL, Cai JJ, Hong DL. Neoplasms in the Bone Marrow Niches: Disturbance of the Microecosystem. Int J Hematol (2017) 105(5):558–65. doi: 10.1007/s12185-017-2193-5
27. Zhu L, Wang P, Zhang W, Li Q, Xiong J, Li J, et al. Plasmacytoid Dendritic Cell Infiltration in Acute Myeloid Leukemia. Cancer Manag Res (2020) 12:11411–9. doi: 10.2147/CMAR.S260825
28. Rodriguez-Abreu D, Johnson ML, Hussein MA, Cobo M, Patel AJ, Secen NM, et al. Primary Analysis of a Randomized, Double-Blind, Phase II Study of the Anti-TIGIT Antibody Tiragolumab (Tira) Plus Atezolizumab (Atezo) Versus Placebo Plus Atezo as First-Line (1L) Treatment in Patients With PD-L1-Selected NSCLC (CITYSCAPE). J Clin Oncol (2020) 38: (15_suppl):9503–3. doi: 10.1200/JCO.2020.38.15_suppl.9503
29. Maruhashi T, Sugiura D, Okazaki IM, Okazaki T. LAG-3: From Molecular Functions to Clinical Applications. J Immunother Cancer (2020) 8(2):e001014. doi: 10.1136/jitc-2020-001014
30. Curti BD, Kovacsovics-Bankowski M, Morris N, Walker E, Chisholm L. Floyd K Et Al: OX40 Is a Potent Immune-Stimulating Target in Late-Stage Cancer Patients. Cancer Res (2013) 73(24):7189–98. doi: 10.1158/0008-5472.CAN-12-4174
31. Wang Z, Chen J, Wang M, Zhang L, Yu L. One Stone, Two Birds: The Roles of Tim-3 in Acute Myeloid Leukemia. Front Immunol (2021) 12:618710. doi: 10.3389/fimmu.2021.618710
32. Vago L, Gojo I. Immune Escape and Immunotherapy of Acute Myeloid Leukemia. J Clin Invest (2020) 130(4):1552–64. doi: 10.1172/JCI129204
33. Zhu L, Kong Y, Zhang J, Claxton DF, Ehmann WC, Rybka WB, et al. Blimp-1 Impairs T Cell Function via Upregulation of TIGIT and PD-1 in Patients With Acute Myeloid Leukemia. J Hematol Oncol (2017) 10(1):124. doi: 10.1186/s13045-017-0486-z
34. Catakovic K, Gassner FJ, Ratswohl C, Zaborsky N, Rebhandl S, Schubert M, et al. TIGIT Expressing CD4+T Cells Represent a Tumor-Supportive T Cell Subset in Chronic Lymphocytic Leukemia. Oncoimmunology (2017) 7(1):e1371399. doi: 10.1080/2162402X.2017.1371399
35. Wang M, Bu J, Zhou M, Sido J, Lin Y, Liu G, et al. Cd8(+)T Cells Expressing Both PD-1 and TIGIT But Not CD226 Are Dysfunctional in Acute Myeloid Leukemia (AML) Patients. Clin Immunol (2018) 190:64–73. doi: 10.1016/j.clim.2017.08.021
36. Williams P, Basu S, Garcia-Manero G, Hourigan CS, Oetjen KA, Cortes JE, et al. The Distribution of T-Cell Subsets and the Expression of Immune Checkpoint Receptors and Ligands in Patients With Newly Diagnosed and Relapsed Acute Myeloid Leukemia. Cancer (2019) 125(9):1470–81. doi: 10.1002/cncr.31896
37. Tan J, Yu Z, Huang J, Chen Y, Huang S, Yao D, et al. Increased PD-1+Tim-3+ Exhausted T Cells in Bone Marrow may Influence the Clinical Outcome of Patients With AML. Biomark Res (2020) 8:6. doi: 10.1186/s40364-020-0185-8
38. Song Y, Wang B, Song R, Hao Y, Wang D, Li Y, et al. T-Cell Immunoglobulin and ITIM Domain Contributes to CD8(+) T-Cell Immunosenescence. Aging Cell (2018) 17(2):e12716. doi: 10.1111/acel.12716
39. Zeidner JF, Knaus HA, Zeidan AM, Blackford AL, Montiel-Esparza R, Hackl H, et al. Immunomodulation With Pomalidomide at Early Lymphocyte Recovery After Induction Chemotherapy in Newly Diagnosed AML and High-Risk MDS. Leukemia (2020) 34(6):1563–76. doi: 10.1038/s41375-019-0693-4
40. Zhao E, Xu H, Wang L, Kryczek I, Wu K, Hu Y, et al. Bone Marrow and the Control of Immunity. Cell Mol Immunol (2012) 9(1):11–9. doi: 10.1038/cmi.2011.47
41. Knaus HA, Berglund S, Hackl H, Blackford AL, Zeidner JF, Montiel-Esparza R, et al. Signatures of CD8+ T Cell Dysfunction in AML Patients and Their Reversibility With Response to Chemotherapy. JCI Insight (2018) 3(21):e120974. doi: 10.1172/jci.insight.120974
42. Schnorfeil FM, Lichtenegger FS, Emmerig K, Schlueter M, Neitz JS, Draenert R, et al. T Cells Are Functionally Not Impaired in AML: Increased PD-1 Expression Is Only Seen at Time of Relapse and Correlates With a Shift Towards the Memory T Cell Compartment. J Hematol Oncol (2015) 8:93. doi: 10.1186/s13045-015-0189-2
43. Davids MS, Kim HT, Bachireddy P, Costello C, Liguori R, Savell A, et al. Ipilimumab for Patients With Relapse After Allogeneic Transplantation. New Engl J Med (2016) 375(2):143–53. doi: 10.1056/NEJMoa1601202
44. Berger R, Rotem-Yehudar R, Slama G, Landes S, Kneller A, Leiba M, et al. Phase I Safety and Pharmacokinetic Study of CT-011, a Humanized Antibody Interacting With PD-1, in Patients With Advanced Hematologic Malignancies. Clin Cancer Res Off J Am Assoc Cancer Res (2008) 14(10):3044–51. doi: 10.1158/1078-0432.CCR-07-4079
45. Gojo I, Stuart RK, Webster J, Blackford A, Varela JC, Morrow J, et al. Multi-Center Phase 2 Study of Pembroluzimab (Pembro) and Azacitidine (AZA) in Patients With Relapsed/Refractory Acute Myeloid Leukemia (AML) and in Newly Diagnosed (≥65 Years) AML Patients. Blood (2019) 134(Supplement_1):832–2. doi: 10.1182/blood-2019-127345
46. Bastholt L, Schmidt H, Bjerregaard JK, Herrstedt J, Svane IM. Age Favoured Overall Survival in a Large Population-Based Danish Patient Cohort Treated With Anti-PD1 Immune Checkpoint Inhibitor for Metastatic Melanoma. Eur J Cancer (2019) 119:122–31. doi: 10.1016/j.ejca.2019.06.022
47. Ravandi F, Assi R, Daver N, Benton CB, Kadia T, Thompson PA, et al. Idarubicin, Cytarabine, and Nivolumab in Patients With Newly Diagnosed Acute Myeloid Leukaemia or High-Risk Myelodysplastic Syndrome: A Single-Arm, Phase 2 Study. Lancet Haematol (2019) 6(9):e480–8. doi: 10.1016/S2352-3026(19)30114-0
48. Le Dieu R, Taussig DC, Ramsay AG, Mitter R, Miraki-Moud F, Fatah R, et al. Peripheral Blood T Cells in Acute Myeloid Leukemia (AML) Patients at Diagnosis Have Abnormal Phenotype and Genotype and Form Defective Immune Synapses With AML Blasts. Blood (2009) 114(18):3909–16. doi: 10.1182/blood-2009-02-206946
49. Kong Y, Zhu L, Schell TD, Zhang J, Claxton DF, Ehmann WC, et al. T-Cell Immunoglobulin and ITIM Domain (TIGIT) Associates With CD8+ T-Cell Exhaustion and Poor Clinical Outcome in AML Patients. Clin Cancer Res Off J Am Assoc Cancer Res (2016) 22(12):3057–66. doi: 10.1158/1078-0432.CCR-15-2626
50. Zhang RJ, Zhai JH, Zhang ZJ, Yang LH, Wang MF, Dong CX. Hypomethylating Agents for Elderly Patients With Acute Myeloid Leukemia: A PRISMA Systematic Review and Meta-Analysis. Eur Rev Med Pharmaco (2021) 25(6):2577–90. doi: 10.26355/eurrev_202103_25421
51. Jia B, Wang L, Claxton DF, Ehmann WC, Rybka WB, Mineishi S, et al. Bone Marrow CD8 T Cells Express High Frequency of PD-1 and Exhibit Reduced Anti-Leukemia Response in Newly Diagnosed AML Patients. Blood Cancer J (2018) 8(3):34. doi: 10.1038/s41408-018-0069-4
Keywords: PD-1, TIGIT, acute myeloid leukemia, immune-checkpoints, T cell, bone marrow
Citation: Xu L, Liu L, Yao D, Zeng X, Zhang Y, Lai J, Zhong J, Zha X, Zheng R, Lu Y, Li M, Jin Z, Hebbar Subramanyam S, Chen S, Huang X and Li Y (2021) PD-1 and TIGIT Are Highly Co-Expressed on CD8+ T Cells in AML Patient Bone Marrow. Front. Oncol. 11:686156. doi: 10.3389/fonc.2021.686156
Received: 31 March 2021; Accepted: 30 July 2021;
Published: 18 August 2021.
Edited by:
Antonio Curti, University of Bologna, ItalyReviewed by:
Evan F. Lind, Oregon Health and Science University, United StatesRory Shallis, Yale University, United States
Copyright © 2021 Xu, Liu, Yao, Zeng, Zhang, Lai, Zhong, Zha, Zheng, Lu, Li, Jin, Hebbar Subramanyam, Chen, Huang and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xin Huang, huangmingxin2001@sina.com; Yangqiu Li, yangqiuli@hotmail.com
†ORCID: Ling Xu, orcid.org/0000-0002-7044-7663
Danlin Yao, orcid.org/0000-0002-0963-8700
Xianfeng Zha, orcid.org/0000-0002-4970-1305
Shaohua Chen, orcid.org/0000-0003-4945-5914
Yangqiu Li, orcid.org/0000-0002-0974-4036
‡These authors have contributed equally to this work