- 1Department of Research and Development, Kibio Inc., Houston, TX, United States
- 2Department of Biomedical Data Science, Stanford University, Stanford, CA, United States
- 3Division of Nephrology, Critical Care and Sleep Medicine, Department of Internal Medicine, University of California, Davis, Davis, CA, United States
- 4Division of Pulmonary, Critical Care and Sleep Medicine, Department of Internal Medicine, University of California, Davis, Davis, CA, United States
- 5Comprehensive Cancer Center, University of California, Davis, Davis, CA, United States
- 6Department of Urology, Chang Gung Memorial Hospital, Linko, Taiwan
- 7Department of Translational Molecular Pathology, The University of Texas MD Anderson Cancer Center, Houston, TX, United States
Our study aimed to elucidate the role of Galectin-1 (Gal-1) role in the immunosuppressive tumor microenvironment (TME) of prostate cancer (PCa). Our previous findings demonstrated a correlation between elevated Gal-1 expression and advanced PCa stages. In this study, we also observed that Gal-1 is expressed around the tumor stroma and its expression level is associated with PCa progression. We identified that Gal-1 could be secreted by PCa cells, and secreted Gal-1 has the potential to induce T cell apoptosis. Gal-1 knockdown or inhibition of Gal-1 function by LLS30 suppresses T cell apoptosis resulting in increased intratumoral T cell infiltration. Importantly, LLS30 treatment significantly improved the antitumor efficacy of anti-PD-1 in vivo. Mechanistically, LLS30 binds to the carbohydrate recognition domain (CRD) of Gal-1, disrupting its binding to CD45 leading to the suppression of T cell apoptosis. In addition, RNA-seq analysis revealed a novel mechanism of action for LLS30, linking its tumor-intrinsic oncogenic effects to anti-tumor immunity. These findings suggested that tumor-derived Gal-1 contributes to the immunosuppressive TME in PCa by inducing apoptosis in effector T cells. Targeting Gal-1 with LLS30 may offer a strategy to enhance anti-tumor immunity and improve immunotherapy.
1 Introduction
Prostate cancer (PCa) stands as the prevalent malignancy and ranks as the second-leading cause of cancer-related mortality in men. Initial treatment for prostate cancer involves androgen deprivation therapy (1). However, a significant challenge arises for patients facing metastatic castration-resistant prostate cancer (mCRPC). mCRPC poses a critical issue, marked by a poor prognosis, and a majority of patients typically undergo disease recurrence within a span of 16 to 18 months (2). Most conventional anticancer therapies and immunotherapy are designed to specifically target PCa cells. However, the immunosuppressive tumor microenvironment (TME) in PCa plays a crucial role in promoting resistance to these therapies (3). For example, T-cell exclusion is prominent in PCa, limiting direct engagement with cancer cells. Immune cells are confined to adjacent stroma and benign areas, often exhibiting immunosuppressive traits (4, 5). This microenvironment is driven by elevated levels of indoleamine 2, 3-dioxygenase (IDO), nitric oxide (NO), interleukin 10 (IL10), transforming growth factor-beta (TGF-β), arginase and adenosine (5–8). Therefore, the development of strategies that target the TME emerges as an attractive approach for the effective treatment of PCa.
The 135 kDa Galectin-1 (Gal-1) protein is encoded by the gene LGALS1 at 22q13.1 (9). Gal-1 is a family of carbohydrate-binding lectins that bind to β-galactoside-containing glycoconjugates (10). Gal-1 is upregulated in cancers (11–15), located both extracellularly and intracellularly, and contributes to many cancer related events including cell proliferation (16), T cell apoptosis (17, 18), angiogenesis (19), and metastatic spread of cancer (20–22). Several Gal-1-binding proteins have been discovered, including H-Ras, integrins, laminins, fibronectin, vitronectin, osteopontin, neuropilin-1, CD45, CD146, and CD326 (23, 24). For example, extracellular Gal-1 recognizes terminal galactose residues β-1,4-linked to LacNAc, which is present in CD45 T cell receptors. Through the binding of LacNAc, Gal-1 can stimulate apoptosis of effector leukocytes (17). Extracellular Gal-1-binding usually occurs through protein-carbohydrate interactions which could be inhibited by lactose (25). In our previous study, we demonstrated that Gal-1 is upregulated in PCa patients and is highly expressed in CRPC cells, but not in androgen-sensitive cells (26). In addition, elevated levels of Gal-1 were observed in the extracellular matrix, and increased stromal Gal-1 expression is associated with poorer outcomes (27). These findings illustrate the importance of Gal-1 in PCa TME. However, whether Gal-1 secreted by PCa cells and its impact on T cell death remain unexplored.
Given that Gal-1 expression in PCa is associated with adverse clinical outcomes, the inhibition of Gal-1 should be considered as a potential treatment approach for CRPC. To target Gal-1, we developed a novel small molecule inhibitor of Gal-1 named LLS30 (26). LLS30 is a benzimidazole-based small molecule with the ability to inhibit the growth and metastasis of PCa xenografts in athymic mice (26, 28, 29). LLS30 functions as a Gal-1 allosteric inhibitor and leads to a decrease in the binding affinity of Gal-1 to its partners. Notably, LLS30 potentiates the antitumor effect of docetaxel and leads to a complete regression of CRPC tumors cells in vivo (26). Furthermore, LLS30 exhibits the capability to inhibit tumor invasion and metastasis in vivo (26).
In this study, our research results demonstrated the secretion of Gal-1 by PCa cells and its potential to induce T cell apoptosis. Furthermore, Gal-1 inhibitor LLS30 was found to significantly suppress the Gal-1 induced T cells apoptosis in vitro and in vivo. Importantly, LLS30 can significantly enhance the antitumor efficacy of anti-PD-1 in vivo. These findings illustrate a potential mechanism of immusuppressive TME in PCa and provide a strategy to improve the therapy outcomes for immunotherapy resistant PCa.
2 Materials and methods
2.1 Cell lines
The cell lines include human AR-positive 22RV1 and AR-negative PC3, as well as mouse AR-positive Myc-CaP PCa cells. These cells were cultured in RPMI1640 (22RV1 and PC3) or DMEM (Myc-CaP), supplemented with 10% fetal bovine serum and 1% penicillin/streptomycin, and maintained in a 5% CO2 atmosphere at 37°C. Routine monthly testing for Mycoplasma contamination was conducted.
2.2 Patients samples
Patient samples were obtained from US Biomax (PR803b, 66cases) and the University of California Davis Comprehensive Cancer Center (83cases) under the protocol titled ‘De-identification and Usage of Existing Pathology Specimens for Research Purposes by the UC Davis Pathology Biorepository,’ with the Institutional Review Board number 293828. Gleason grading was done on tissue core based on 2005 ISUP modified system (30): Gleason score<6 (low grade), = 7 (intermediate grade), and >8 (high grade). Human tissues covered 38 cases of benign hyperplasia (BHP), 32 cases of low grade, 24 cases of intermediate grade and 55 cases of high grade PCa.
2.3 Co-cultures of PCa cells with peripheral blood mononuclear cells
To initiate the experiment, 5 × 103 22RV1 cells were seeded in each well of a 48-well plate containing 100 μL of complete culture medium (RPMI1640 supplemented with 10% fetal bovine serum and 1% penicillin/streptomycin). The plates were then left to incubate overnight, allowing the cells to adhere, a process taking approximately 24 hours. By the following day, it was estimated that the number of cancer cells had doubled to 1 × 104 cells per well, forming the basis for calculating Effector-to-Target (E:T) ratios. For co-culture experiments, PBMCs serving as effector cells, were mixed with 22Rv1 cells at a ratio of 3:1. PBMCs were cultured in complete medium and activated with IL-2 (5 ng/mL). After 24 hours of activation, PBMCs were preincubated with anti-PD-1 (Nivolumab, selleckchem) at concentrations of 10 μg/ml, 20 μg/ml, and 40 μg/ml for 30 minutes. Subsequently, they were added to the cancer cells in the presence or absence of 1 μM or 2 μM LLS30 for an additional 24 hours. The culture with PBMCs and 22RV1 cells without nivolumab and LLS30 addition served as a control cultures. Following the 24-hour incubation period, the culture medium was aspirated, and the cells were washed twice with PBS to remove any suspended PBMCs and dead 22RV1 cells. The viability of the adherent 22RV1 cells was assessed using the MTT assay.
2.4 Gal-1 knockdown cells establishment
PCa cells were seeded on the 6 well plates and transfected with 50 nM negative control mimic or mixed siRNA against Gal-1 or control siRNA (Qiagen)), using Lipofectamine 2000 transfection reagent (Life Technologies) according to the manufacturer’s instructions. PCa cells were infected with control or Gal-1 shRNA lentiviral particles (Santa Cruz) at an MOI of 10 for 24 h in the presence of 8 μg/ml polybrene. Immunoblotting was employed to assess Gal-1 expression following transfection with siGal-1 or shGal-1.
2.5 Immunoblotting analysis
Immunoblotting was done as described previously (26). Antibodies used in immunoblotting were Gal-1 (Abcam), beta-actin (Cell Signaling) and CD45 (Cell signaling). Immunoblotting bands were quantified by Image J software.
2.6 Quantification of secreted Gal-1 in conditioned medium from PCa cells
2 × 105 22RV1, PC3 and Myc-CaP cells were seeded per well into 6-well plates and allowed to attach for 24 hours. Following the 24-hour incubation period, cell supernatants were collected as conditioned medium. The concentration of secreted Gal-1 in the conditioned medium was measured using a Gal-1 ELISA kit (R&D Systems) following the manufacturer’s protocol.
2.7 Gal-1 mediated T cell apoptosis assay
Cells transfected with sicontrol or siGal-1 in culture medium were seeded onto a 6-well tissue culture plate at a density of 5×105 cells per well. Following an overnight incubation, conditioned medium was obtained by centrifugation of the supernatants. To pharmacologically inhibit secreted Gal-1, LLS30 2µM was added to the conditioned medium from wild-type PCa cells. CD8+ T cells were isolated from human and mouse PBMCs using the Dynabeads CD8 Positive Isolation Kit as per the manufacturer’s instructions (Invitrogen). Then these cells were cultured in conditioned medium for 24 hours, followed by an apoptosis assay to detect Caspase-3/7 activity using a luminescent Caspase-Glo 3/7 assay kit (Promega).
2.8 Gal-1 binding to T cells
Fluorescein isothiocyanate (FITC) from Sigma was utilized to label recombinant Gal-1 according to the manufacturer’s instructions. This labeling occurs through the reaction of the isothiocyanate group with amino terminal and primary amines in proteins. Gal-1 at a concentration of 2 µM was treated with LLS30 at 2 µM or DMSO 0.02% for 4 hours at 4°C. T cells were collected by centrifugation and fixed in 4% paraformaldehyde. Non-specific protein binding was blocked by adding 5% BSA. T Cells were then incubated with Gal-1/LLS30 or Gal-1/DMSO. Cytospin 200 μL of the stained cell suspension onto a glass slide. The cells were mounted and photographed using a fluorescence microscope (Zeiss).
2.9 Co-immunoprecipitation assay
We first treated recombinant Gal-1 2 µM with LLS30 2 µM or DMSO 0.02% for 4 hours at 4°C. In addition, 2 x 106 T cells were cultured for membrane protein extraction using the ProteoExtract® Native Membrane Protein Extraction Kit (Millipore). 30 μg of membrane protein was added in Gal-1/LLS30 or Gal-1/DMSO for 4 hours at 4°C. In the immunoprecipitation step, 10 µL of anti-Gal-1 (Abcam) was gently mixed with 100 µL of a Protein A/G Sepharose slurry (Abcam) and incubated for 4 hours at 4°C. Subsequently, the Sepharose slurry was added to protein mixtures with rotary agitation overnight at 4°C. Sepharose was washed with PBS, and then elution of bound proteins was performed with 150 mM glycine, pH 2.5, for 10 minutes. The eluent was then neutralized by adding 10 mL neutralization buffer (Tris, pH 8.0) and subject to immunoblotting to assess the precipitation of proteins.
2.10 Immunohistochemistry
Tissue sections were de-waxed using xylene twice and rehydration with 100% ethanol for 5 minutes each, followed by 95% and 80% ethanol for 5 minutes each, and then rinsed in PBS. Antigen retrieval was carried out in a 10 mM, pH 6.0 sodium citrate buffer at 95–100°C for 20 minutes. After cooling to room temperature, the sections were rinsed with PBS, and endogenous peroxidase was blocked with 1% H2O2. Non-specific binding sites were blocked with Power Block (BioGenex) for 5 minutes at room temperature. Next, the tissue sections were incubated overnight with specific antibodies against Gal-1 (Abcam), Caspase-3 (Cell Signaling), or CD8 (Cell Signaling). Following rinsing with PBS, the sections were incubated with a biotin-conjugated goat anti-rabbit IgG (BioGenex) as the second antibody. Then the sections were incubated with streptavidin-conjugated HRP (BioGenex) for 20 minutes at room temperature. HRP activity was detected using diaminobenzidine tetrahydrochloride (DAB) as the substrate (BioGenex). Nuclei were counterstained with hematoxylin (Cell Signaling). CD8-positive cells were counted in three randomly chosen areas.
2.11 In vivo xenograft tumor
LLS30 stock solutions (6X) were prepared in a mixture of 50% absolute alcohol and 50% Tween-80 to achieve a concentration of 6 mg/ml for Myc-CaP xenografts. Male FVB/N mice were obtained from the Jackson Laboratory, and 1 x 106 Myc-CaP cells were subcutaneously injected into the right side of the mouse dorsal flank. The tumors were allowed to grow to approximately 100 mm3. Four groups of mice bearing Myc-CaP xenografts were designed, with 6 mice per group. Mice bearing Myc-CaP xenografts received intravenous administration for 2 weeks with one of the following treatments (1): vehicle control (8.7% alcohol and 8.7% Tween-80 in PBS), (2) 10 mg/kg LLS30 daily for 5 successive days, followed by a two-day break, and then for another 5 successive days, (3) 10 mg/kg anti-PD-1 given every other day for 4 shots, or (4) a combination of LLS30 and anti-PD-1 via intravenous injection. Tumor size and body weight were measured twice a week, and tumor volumes were measured using the formula (length x width2)/2.
2.12 Transcriptome sequencing and enrichment analysis
The RNA-seq data analysis involved extracting total RNA from both control and LLS30-treated 22RV1 cells using the PureLink RNA Mini Kit (Invitrogen) following the manufacturer’s instructions. RNA quality was evaluated using the Agilent 2100 Bioanalyzer system (Agilent Technologies). The mRNA sequencing library was prepared, and paired-end sequencing was conducted using the Illumina HiSeq 4000 Sequencing System. Differentially expressed genes (DEGs) were identified with a 1.5-fold change and with a significance level of P < 0.05. Volcano plots were generated using the Sigmaplot software. GO and Reactome analyses were performed on the Enrichr website (https://maayanlab.cloud/Enrichr/) using a list of gene symbols.
2.13 Statistical analysis
Gal-1 expression levels in tumor stroma were assessed using the following scoring system: 0 for negative, 1 for low intensity, 2 for moderate intensity, 3 for high intensity, and 4 for very high intensity. The CompuSyn software was utilized to calculate the combination index, with values below 1 indicating synergy. In vitro experiments were conducted in triplicate across two independent experiments, and the results are presented as the mean ± SD. The student’s t-test (two-tailed) was employed for comparing datasets between two groups with similar variance. P value < 0.05 was considered indicative of a statistically significant difference. Statistical differences, when compared with controls, are denoted as * (P < 0.05), ** (P < 0.01), or *** (P < 0.001).
3 Results
3.1 Gal-1 overexpression in the PCa TME
We previously demonstrated a correlation between elevated Gal-1 expression and advanced stages of PCa. However, our investigation did not encompass the assessment of Gal-1 expression in the PCa stroma. To address this gap, we reexamined our earlier TMA IHC analysis of human PCa to investigate the presence of Gal-1 in tumor stroma. The analysis results showed that, in addition to its expression in glandular epithelial cells, overexpressed Gal-1 is also observed in the stroma (Figure 1A). Of note, the stromal accumulation of Gal-1 exhibited a significantly progressive upregulation, transitioning from low to intermediate grades (p < 0.05 for low vs. intermediate-grade PCa) and further escalating in high-grade PCa (p < 0.001 for intermediate- vs. high-grade PCa) (Figure 1B). Furthermore, these findings are consistent with Castronovo’s study, wherein their dataset indicated that Gal-1 expression in cancer-associated stroma is associated with malignant tissue and tumor aggressiveness (27). These clinical observations underscore the association between increased Gal-1 expression in the stroma and the progression of PCa. This suggests that secreted Gal-1 may play a role in modulating the immune response, potentially fostering tumor progression.
Figure 1 Expression of Gal-1 in the human PCa microenvironment. (A) The expression levels of Gal-1were detected by IHC. IHC showed that Gal-1 was localized in tumor and tumor stroma. (B) Tissue Gal-1 expression in 38, 32, 24 and 55 samples of BHP, low-, intermediate- and high-grade PCa tissues, respectively. Sampling distribution of Gal-1 expression was displayed by Box-Plot (dash line: mean; lines above and below the dash line, third quartile to the first quartile; lines above and below the box, maximum and minimum; dots, outlier values). *P < 0.05, ***P < 0.001; Student’s t test.
3.2 Conditioned medium from Gal-1 expressing PCa cells induces T cell apoptosis
It is known that tumor-secreted Gal-1 can bind to glycosylated receptors on immune cells and trigger the apoptosis of T cells in the TME (31). However, in the context of PCa, it remains unclear whether Gal-1 can disseminate from cancer cells and induce T cell apoptosis. To investigate this, we first elucidate whether Gal-1 is secreted from cancer cells. We performed immunoblot analysis to verify the endogenous Gal-1 expression in Gal-1 expressing PCa cells including human 22RV1 (AR positive), PC3 (AR negative) and mouse Myc-CaP (AR positive) (Figure 2A). Following this, we cultured these PCa cells in RPMI1640 (for 22RV1 and PC3) or DMEM (for Myc-CaP), supplemented with 10% fetal bovine serum and 1% penicillin/streptomycin for 48 hours to generate conditioned medium. We then used an ELISA assay to detect the secreted Gal-1 in the conditioned medium. The ELISA results confirmed that Gal-1 was indeed secreted from PCa cells into the culture medium (Figure 2B). Next, conditioned medium from Gal-1 expressing cells were collected and treated with T cells to assessed for its apoptotic inducing effect on T cells. As shown in Figure 2C, the conditioned medium from PCa cells exerts the ability to induce apoptosis. We next determine the effects of tumor secreted Gal-1 on induced T-cell apoptosis in vivo. To achieve this, we subcutaneously injected the Myc-CaP cells into syngeneic FVB/N mice. Tumors were collected upon reaching a diameter of 1 cm and prepared for clinical pathology examination and IHC analysis. HE stained showed that lymphocytes have dark nuclei and a slender periphery of pale cytoplasm. Of note, we observed apoptotic cell debris with DNA fragmentation around the tumor stroma expressing Gal-1 (Figure 2D, basophilic nuclear remnants; arrows). Given that Caspase-3 is essential for DNA fragmentation during apoptosis, we detected caspase-3 in these tissues. IHC results confirmed the presence of caspase-3 in these tissues, supporting the apoptotic processes observed. These findings highlight the role of secreted Gal-1 in inducing T cell apoptosis in PCa.
Figure 2 Conditioned medium from Gal-1 expressing PCa cells induces T cell apoptosis. (A) Immunoblotting illustrated endogenous Gal-1 expression in PCa cells. Conditioned medium (CM) from PCa cells were collected for ELISA and apoptosis assays, as depicted in (B, C), respectively. (B) ELISA results showed the detection of secreted Gal-1 in CM from PCa cells. (C) Apoptosis assay for CD8+ T cells. CM from PCa cells induced T cell apoptosis compared to the medium-only control. (D) H&E and enlarged views (upper), Gal-1 staining (lower left), and Caspase-3 staining (lower right). Arrows in the image indicate the location of apoptotic cell debris with DNA fragmentation. *P < 0.05, ***P < 0.001; Student’s t test, n = 3. Data shown are mean ± s.d.
3.3 Gal-1 knockdown PCa cells prevented T cell apoptosis
To investigate the specific effects of secreted Gal-1 on T cell apoptosis, Gal-1 knockdown was conducted in Gal-1-expressing PCa cell lines. Immunoblotting confirmed the efficiency of Gal-1 knockdown (Figure 3A) and reduced secretion of Gal-1 as detected by ELISA (Figure 3B). Conditioned medium from sicontrol or siGal-1 PCa cells was then used to treat T cells. The CM from PCa Gal-1 knockdown cells showed a reduced ability to induce apoptosis compared to Gal-1 sicontrol cells (Figure 3C). In in vivo experiments, histologic analysis of the resected tumor revealed fewer apoptotic T cells around the shGal-1 tumor, indicating that Gal-1 knockdown prevents T cell apoptosis. Notably, CD8 positive cells were upregulated in the Gal-1 knockdown tumor (Figure 3D). These results collectively suggested that the expression of Gal-1 in the TME may induce T cell apoptosis, potentially contributing to the absence of tumor-infiltrating lymphocytes (TILs).
Figure 3 Gal-1 knockdown in PCa cells prevents T cell apoptosis. (A) Endogenous Gal-1 expression in PCa cells was effectively suppressed by two different siRNAs at 72 hours post-transfection. The conditioned medium (CM) from PCa were then collected for ELISA and apoptosis assays. (B) ELISA results illustrated the detection of secreted Gal-1 in CM from PCa sicontrol and PCa siGal-1 cells. (C) In the apoptosis assay for CD8+ T cells, CM from PCa sicontrol cells induced T cell apoptosis, but not PCa siGal-1 cells. (D) H&E-stained images (upper) and CD8+ stained images (lower). ***P < 0.001; Student’s t test, n = 3. Data shown are mean ± s.d.
3.4 Pharmacological inhibition of secreted Gal-1 function by LLS30 suppressed T cell apoptosis
Extracellular Gal-1 carbohydrate recognition domain (CRD) recognized terminal galactose residues β-1,4-linked to LacNAc as presented in CD45 T cell receptors resulting in the apoptosis of effector leukocytes (17). Our previous studies showed that LLS30 binds the CRD of Gal-1 as an allosteric inhibitor leading to a decrease in binding affinity of Gal-1 to its binding partners. We further conducted a cell-based surface glycan binding assay to confirm the inhibitory effect of LLS30 on the adhesion of Gal-1 to glycosylated receptors on T cells. FITC was utilized for the labeling of Gal-1. We observed a significant decrease in fluorescence signals in the T cells in the presence of LLS30 compared to DMSO control cells (Figure 4A), indicating that LLS30 blocks FITC-labeled Gal-1 binding to glycosylated receptors on T cell surfaces. Given that CD45 on T cells is a primary binding partner for Gal-1 in apoptosis induction, we used Co-immunoprecipitation (co-IP) to examine whether LLS30 reduces the binding of CD45 to Gal-1. The co-IP results confirmed this reduction in CD45 binding to Gal-1 in the presence of LLS30 (Figure 4B). Next, we examined whether interfering Gal-1/CD45 interaction by LLS30 would reduce Gal-1-mediated apoptosis in T cells. We pretreated conditioned medium from Gal-1 wild-type PCa cells with or without LLS30 (2 µM), and co-culture with T cells to assessed the T cell apoptosis. The results showed that Gal-1 triggered T cells apoptosis can be inhibited with LLS30 at 2 µM which is not toxic for human PBMCs cells (Figure 4C). Taken together, these studies indicated that LLS30 binds to the Gal-1 CRD and interfere Gal-1 binding to CD45, resulting in a suppression of T cell apoptosis.
Figure 4 LLS30 disrupts Gal-1-carbohydrate interactions. (A) Cell surface glycan binding assay; 2 μM of LLS30 inhibited the adhesion of FITC labeled Gal-1 to cell surface glycan. (B) Left: A schematic illustrates the workflow for co-immunoprecipitation experiments analyzing the interaction between Gal-1 and CD45 after treatment with LLS30 or vehicle DMSO. Right: Co-immunoprecipitation of Gal-1 with CD45. Membrane protein extracts from T cells was treated with Gal-1/LLS30 or Gal-1/DMSO and immunoprecipitated with anti-Gal-1 antibody. Gal-1/LLS30 or Gal-1/DMSO treated T cell membrane proteins without immunoprecipitated with anti-Gal-1 antibody serve as inputs. Both inputs and the resultant immunoprecipitates were examined via immunoblotted with CD45 antibodies. (C) Apoptosis assay for CD8+ T cells. LLS30 treated CM from PCa cells suppressed induced T cell apoptosis compared to non LLS30-treated CM. **P < 0.01, ***P < 0.001; Student’s t test, n = 3. Data shown are mean ± s.d.
3.5 LLS30 potentiates antitumor activity of anti-PD1 in immunotherapy resistant PCa
To assess the impact of Gal-1 inhibition by LLS30 on the TME, we subcutaneously injected the Myc-CaP cells into syngeneic FVB/N mice. We observed enhanced infiltration of CD8+ T cells in the tumors of mice treated with LLS30 (Figure 5A). These results suggest that LLS30 treatment prevents T cell apoptosis, resulting in the increased presence of TILs. PCa is typically considered immunologically cold tumors with minimal T cell infiltrates and limited responsiveness to checkpoint anti-PD-1 therapy. Given the increase in TILs with LLS30 treatment, we investigated whether LLS30 treatment could enhance the outcomes of anti-PD-1 therapy. We conducted in-vitro T cell killing assays by co-culturing 22RV1 cells with peripheral blood mononuclear cells (PBMCs). The results revealed that the combination of anti-PD-1 with LLS30 synergistically induced cell death, as evidenced by a combination index ranging from 0.2 to 0.4 (Figure 5B). We further assessed the combination effects on tumor growth in vivo. Considering the observed deficiency of T cells in Myc-CaP tumors (Figure 5A), we hypothesized that anti-PD1 therapy might not effectively suppress Myc-CaP tumor growth, and anticipated that inhibiting Gal-1 with LLS30 could enhance the efficacy of anti-PD1 therapy by improving T cell infiltration. In vivo study showed that anti-PD1 antibody alone had no effect on Myc-CaP tumor growth (Figure 5C). However, LLS30 at 10 mg/kg dose significantly suppressed the growth of Myc-CaP tumor, and combination of LLS30 with anti-PD-1 caused better tumor regression (Figure 5C). Moreover, CD-8+ cells were observed in LLS30-treated and combination LLS30/anti-PD1–treated tumor (Figure 5D).
Figure 5 LLS30 can potentiate the anticancer effects of anti-PD1 in vivo. (A) CD8+ stained images after treatment with 10 µM LLS30. (B) Quantification of 22RV1 cells in mixed PBMCs-cancer co-cultures was conducted after 24 hours of treatment with anti-PD-1 and LLS30 at indicated concentrations. Combination index (C.I.) below 1 indicating synergy. (C) Tumor growth curves (n=6). Effect of LLS30, anti-PD1 and combination on growth of Myc-CaP tumor xenografts in vivo. (D) H&E and CD8+ stained images. Data shown are mean ± s.d. **P < 0.01, ***P< 0.001.
3.6 LLS30 modulates cancer cell-intrinsic pathways impacting anti-tumor immunity
Increased evidence indicates that tumor-intrinsic signaling can modulate the immune response to the tumor (32). Therefore, in addition to establishing that one of the functions of LLS30 is to prevent Gal-1-mediated T cell apoptosis, we conducted a detailed exploration of the gene regulatory processes responsible for inducing cell death through the effects of LLS30. Our analysis involved examining RNA-Seq datasets to identify functional enrichments of differentially expressed genes in 22RV1 cells. RNA-Seq analysis unveiled the upregulation of 968 genes and down-regulation of 408 genes (with a 1.5-fold change and P < 0.05) following LLS30 treatment (Figure 6A). Functional pathway analysis performed in the differentially expressed genes revealed that LLS30 had effects on multiple important cellular pathways. In the Reactome database, the top five upregulated pathways include the unfolded protein response (UPR), PERK regulates gene expression, Asparagine N-linked glycosylation, ATF4 activates genes, and RAF-independent MAPK1/3 activation (Figure 6B). UPR, PERK and ATF4 are essential for induction of ER stress, suggesting LLS30 has a notable impact on induction of ER stress. Further exploration of significant pathway using the Gene Ontology (GO) biological processes database revealed that among the top five upregulated pathways, LLS30 significantly influenced ER stress (Figure 6B). In addition, the Reactome analysis highlighted down-regulated pathways were implicated in cellular communication and the regulation of protein folding. The GO database indicated that the down-regulated genes were associated with the maintenance of cell survival and proliferation (Figure 6C). Importantly, the activation of ER stress signaling has been reported to initiate a cascade of events for immunogenic apoptosis in tumor cells and subsequent recruitment and activation of T cells (33). This suggests the potential of LLS30 to induce immunogenic cell death, providing insights into its role in enhancing anti-tumor immunity and improving immunotherapy.
Figure 6 Transcriptomic analysis for LLS30 effects in 22RV1 PCa cells. (A) Volcano plot. The log2 fold change and log10 p-value indicated the expression level and significance for each gene. Each dot represents one gene. Black dots represent no significant differential expressed genes between control and LLS30 treatment, the green dots represent downregulated genes (FC<1.5, p<0.05) and red dots represent upregulated genes (FC>1.5, p<0.05). Reactome and GO biological process analysis of the (B) upregulated genes and (C) downregulated genes between control and LLS30 treatment. The top 10 most enriched pathways in each category were listed. Pathways significantly enriched for upregulated genes are represented by red bars, while those enriched for downregulated genes are denoted by green bars.
4 Discussion
Within the TME, Gal-1 significantly contributes to tumor immune evasion. In head and neck cancer, Chawla et al. discovered that 58.8% of the patient cohort exhibited moderate to marked lymphocyte infiltrates, consisting of T cells, B cells, and FoxP3-expressing T cells, while the presence of Gal-1 staining within lymphocyte areas of the tumor was significantly correlated with a poorer patient progonosis (34). In colorectal cancer, Rabinovich et al. demonstrated Gal-1 derived from the tumor promotes immunosuppression in a syngeneic colorectal cancer model by inducing CD8+ regulatory T cells (35). Our research findings contribute to the existing knowledge by providing evidence that Gal-1 plays a role in orchestrating an immunosuppressive microenvironment in PCa through the induction of apoptosis in effector T cells.
Tumor with poor T cell trafficking is one of the challenges for checkpoint blockade in cancer immunotherapy (36). Here we demonstrated that the inhibition of Gal-1 by LLS30 suppressed T cell apoptosis, contributing the presence of TILs and consequently enhancing the effects of anti-PD1 immunotherapy to inhibit tumor growth. Our research highlights the potential synergy between Gal-1 inhibition and an-PD-1 immunotherapy. In a recent investigation conducted by He et al., their IHC studies revealed that high Gal-1 expression was associated with a reduced number of TILs and PD-1 expression in lung, esophageal, and colorectal cancers (37). To assess the impact of combination therapy, they employed anti-PD-1 treatment in conjunction with the Gal-1 inhibitor OTX008, examining its effects on tumor growth in LLC and B16-F10 tumor models. They found that the efficacy of the combination of anti-PD-1 and Gal-1 inhibitor was significantly better than that of therapy alone. Moreover, Nambiar et al., demonstrated that expression of Gal-1 is associated with poor overall survival in head and neck cancer patients treated with immune checkpoint inhibitors (38). In vivo experiments also showed that Gal-1 deletion or blockade with an antibody increased intratumorally T-cell infiltration and improve the response to anti-PD1 therapy with or without radiotherapy (38). Mechanistically, Gal-1 causes immune evasion by preventing T-cell migration into tumors, by reprogramming the tumor endothelium to upregulate cell-surface PD-L1 and Gal-9 (38). Collectively, the results from these studies, indicating that Gal-1 inhibition enhances anti-PD-1 therapy, align with our findings and substantiate our notions.
Secreted Gal-1 in the TME may also function on other immune cells. Verschuere et al. demonstrated that the suppression of Gal-1 derived from glioma significantly reduced the presence of brain-infiltrating macrophages and myeloid-derived suppressor cells in an orthotopic GL261 mouse glioma model (39). Wandall et al. demonstrated that Gal-1 induces a tumor-associated macrophage-like phenotype in monocyte-derived macrophages, leading to the induced expression of PD-L1 and IDO1 in M2-like macrophages, which are key factors in immune regulation by suppressing immune responses within TME (40).
Considerable research has been conducted regarding inhibiting Gal-1 in malignancy in experimental models. Thiodigalactoside, a disaccharides, has proven effective in reducing the progression of breast cancer when administered in conjunction with vaccine immunotherapy (41–44). OTX008, a small molecule Gal-1 inhibitor, exhibits better stability compared to saccharides. Mechanistically, OTX008 inhibits ERK and causes G2/M cell cycle arrest in a panel of human cancer cell lines (28). Experimental data suggest that OTX008 alone or in combination with other standards of care may be an effective treatment for solid tumor (45, 46). In 2012, the first phase of clinical trial aimed at evaluating the effects of OTX008 for the treatment of advanced solid tumors was conducted (ClinicalTrials.gov: NCT01724320). So far, the fate of this clinical trial remains unknown because no reporting of outcome in this clinical has been released. Notably, the existing repertoire of Gal-1 inhibitors is limited in effectiveness, and none have progressed for human use. This signifies the urgency to develop a Gal-1 inhibitor that is both effective and safe for human application. To our knowledge, our group stands as the singular entity delving into the study of Gal-1 inhibition by LLS30 in PCa. We have demonstrated that LLS30 can inhibit the growth and metastasis of PCa xenografts (26). Moreover, LLS30 has the ability to potentiate the antitumor effects of docetaxel chemotherapy and anti-PD-1 immunotherapy in vivo. LLS30 is a promising small molecule compound that warrants further development for human clinical trials.
In conclusion, the high expression of Gal-1 in PCa TME contributes to the induction of T cell apoptosis. Importantly, by inhibiting Gal-1 function with LLS30, we successfully suppressed T cell apoptosis, resulting in increased TILs and enhanced the antitumor efficacy of anti-PD-1. Furthermore, the RNA-seq analysis results unveiled a novel mechanism of action for LLS30, highlighting a crucial link between its tumor-intrinsic oncogenic effects and anti-tumor immunity. Targeting Gal-1 with LLS30 represents a promising and novel therapeutic strategy to enhance immunotherapeutic approaches.
Data availability statement
The data presented in the study are deposited in the BioProject database, accession number PRJNA1120033.
Ethics statement
Human tissues were obtained from US Biomax (PR803b, 66cases) and the University of California Davis Comprehensive Cancer Center (83cases) under the protocol titled 'De-identification and Usage of Existing Pathology Specimens for Research Purposes by the UC Davis Pathology Biorepository,' with the Institutional Review Board number 293828. The animal study was approved by University of California, Davis IACUC protocols. The study was conducted in accordance with the local legislation and institutional requirements.
Author contributions
H-CW: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Resources, Validation, Writing – original draft, Writing – review & editing. RX: Conceptualization, Data curation, Formal analysis, Investigation, Resources, Software, Validation, Visualization, Writing – review & editing. W-HC: Data curation, Formal analysis, Investigation, Software, Writing – review & editing. S-WH: Data curation, Formal analysis, Investigation, Methodology, Software, Writing – review & editing. C-TW: Data curation, Formal analysis, Writing – review & editing. C-HC: Investigation, Methodology, Resources, Supervision, Writing – review & editing. T-CS: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study is supported by grants funded by the U.S. Department of Defense Prostate Cancer Research Program (W81XWH-19-1-0532) to T-CS.
Acknowledgments
We appreciate Dr. Allen C. Gao and Dr. Emanual Maverakis for their valuable contributions to the scientific discussions.
Conflict of interest
Author H-CW was employed by Kibio Inc.
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.
References
1. Huggins C HC. The effect of castration, of estrogen and of androgen injection on serum phosphatases in metastatic carcinoma of the prostate. Cancer Res. (1941) 1:293–7. doi: 10.3322/canjclin.22.4.232
2. Marques RB, Dits NF, Erkens-Schulze S, van Weerden WM, Jenster G. Bypass mechanisms of the androgen receptor pathway in therapy-resistant prostate cancer cell models. PloS One. (2010) 5:e13500. doi: 10.1371/journal.pone.0013500
3. Ge R, Wang Z, Cheng L. Tumor microenvironment heterogeneity an important mediator of prostate cancer progression and therapeutic resistance. NPJ Precis Oncol. (2022) 6:31. doi: 10.1038/s41698-022-00272-w
4. Bilusic M, Madan RA, Gulley JL. Immunotherapy of prostate cancer: facts and hopes. Clin Cancer Res an Off J Am Assoc Cancer Res. (2017) 23:6764–70. doi: 10.1158/1078-0432.CCR-17-0019
5. Molon B, Ugel S, Del Pozzo F, Soldani C, Zilio S, Avella D, et al. Chemokine nitration prevents intratumoral infiltration of antigen-specific T cells. J Exp Med. (2011) 208:1949–62. doi: 10.1084/jem.20101956
6. Munn DH, Bronte V. Immune suppressive mechanisms in the tumor microenvironment. Curr Opin Immunol. (2016) 39:1–6. doi: 10.1016/j.coi.2015.10.009
7. Joyce JA, Fearon DT. T cell exclusion, immune privilege, and the tumor microenvironment. Sci (New York NY). (2015) 348:74–80. doi: 10.1126/science.aaa6204
8. Wu AA, Drake V, Huang HS, Chiu S, Zheng L. Reprogramming the tumor microenvironment: tumor-induced immunosuppressive factors paralyze T cells. Oncoimmunology. (2015) 4:e1016700. doi: 10.1080/2162402X.2015.1016700
9. Barondes SH, Castronovo V, Cooper DN, Cummings RD, Drickamer K, Feizi T, et al. Galectins: a family of animal beta-galactoside-binding lectins. Cell. (1994) 76:597–8. doi: 10.1016/0092-8674(94)90498-7
10. Earl LA, Bi S, Baum LG. Galectin multimerization and lattice formation are regulated by linker region structure. Glycobiology. (2011) 21:6–12. doi: 10.1093/glycob/cwq144
11. Zhang P, Zhang P, Shi B, Zhou M, Jiang H, Zhang H, et al. Galectin-1 overexpression promotes progression and chemoresistance to cisplatin in epithelial ovarian cancer. Cell Death disease. (2014) 5:e991. doi: 10.1038/cddis.2013.526
12. Carlini MJ, Roitman P, Nuñez M, Pallotta MG, Boggio G, Smith D, et al. Clinical relevance of galectin-1 expression in non-small cell lung cancer patients. Lung Cancer (Amsterdam Netherlands). (2014) 84:73–8. doi: 10.1016/j.lungcan.2014.01.016
13. Jung EJ, Moon HG, Cho BI, Jeong CY, Joo YT, Lee YJ, et al. Galectin-1 expression in cancer-associated stromal cells correlates tumor invasiveness and tumor progression in breast cancer. Int J Cancer. (2007) 120:2331–8. doi: 10.1002/ijc.22434
14. White NM, Masui O, Newsted D, Scorilas A, Romaschin AD, Bjarnason GA, et al. Galectin-1 has potential prognostic significance and is implicated in clear cell renal cell carcinoma progression through the HIF/mTOR signaling axis. Br J Cancer. (2014) 110:1250–9. doi: 10.1038/bjc.2013.828
15. Martinez-Bosch N, Fernandez-Barrena MG, Moreno M, Ortiz-Zapater E, Munne-Collado J, Iglesias M, et al. Galectin-1 drives pancreatic carcinogenesis through stroma remodeling and hedgehog signaling activation. Cancer Res. (2014). doi: 10.1158/0008-5472.Can-13-3013
16. Liu FT, Rabinovich GA. Galectins as modulators of tumour progression. Nat Rev Cancer. (2005) 5:29–41. doi: 10.1038/nrc1527
17. Perillo NL, Pace KE, Seilhamer JJ, Baum LG. Apoptosis of T cells mediated by galectin-1. Nature. (1995) 378:736–9. doi: 10.1038/378736a0
18. Rubinstein N, Alvarez M, Zwirner NW, Toscano MA, Ilarregui JM, Bravo A, et al. Targeted inhibition of galectin-1 gene expression in tumor cells results in heightened T cell-mediated rejection. Cancer Cell. (2004) 5:241–51. doi: 10.1016/S1535-6108(04)00024-8
19. Stanley P. Galectin-1 pulls the strings on VEGFR2. Cell. (2014) 156:625–6. doi: 10.1016/j.cell.2014.01.059
20. Hsu YL, Wu CY, Hung JY, Lin YS, Huang MS, Kuo PL. Galectin-1 promotes lung cancer tumor metastasis by potentiating integrin alpha6beta4 and Notch1/Jagged2 signaling pathway. Carcinogenesis. (2013) 34:1370–81. doi: 10.1093/carcin/bgt040
21. Qian D, Lu Z, Xu Q, Wu P, Tian L, Zhao L, et al. Galectin-1-driven upregulation of SDF-1 in pancreatic stellate cells promotes pancreatic cancer metastasis. Cancer letters. (2017) 397:43–51. doi: 10.1016/j.canlet.2017.03.024
22. Wu M-H, Hong T-M, Cheng H-W, Pan S-H, Liang Y-R, Hong H-C, et al. Galectin-1-mediated tumor invasion and metastasis, up-regulated matrix metalloproteinase expression, and reorganized actin cytoskeletons. Mol Cancer Res. (2009) 7:311–8. doi: 10.1158/1541-7786.MCR-08-0297
23. Thijssen VL, Heusschen R, Caers J, Griffioen AW. Galectin expression in cancer diagnosis and prognosis: A systematic review. Biochim Biophys Acta. (2015) 1855:235–47. doi: 10.1016/j.bbcan.2015.03.003
24. Gordon-Alonso M, Bruger AM, van der Bruggen P. Extracellular galectins as controllers of cytokines in hematological cancer. Blood. (2018) 132:484–91. doi: 10.1182/blood-2018-04-846014
25. Camby I, Le Mercier M, Lefranc F, Kiss R. Galectin-1: a small protein with major functions. Glycobiology. (2006) 16:137r–57r. doi: 10.1093/glycob/cwl025
26. Shih TC, Liu R, Wu CT, Li X, Xiao W, Deng X, et al. Targeting galectin-1 impairs castration-resistant prostate cancer progression and invasion. Clin Cancer Res. (2018) 24:4319–31. doi: 10.1158/1078-0432.CCR-18-0157
27. van den Brule AF, Castronovo V, Waltregny D. Increased expression of galectin-1 in carcinoma-associated stroma predicts poor outcome in prostate carcinoma patients. J Pathol. (2001) 193:80–7. doi: 10.1002/(ISSN)1096-9896
28. Astorgues-Xerri L, Riveiro ME, Tijeras-Raballand A, Serova M, Rabinovich GA, Bieche I, et al. OTX008, a selective small-molecule inhibitor of galectin-1, downregulates cancer cell proliferation, invasion and tumour angiogenesis. Eur J Cancer. (2014) 50:2463–77. doi: 10.1016/j.ejca.2014.06.015
29. Dings RPM, Miller MC, Nesmelova I, Astorgues-Xerri L, Kumar N, Serova M, et al. Antitumor agent calixarene 0118 targets human galectin-1 as an allosteric inhibitor of carbohydrate binding. J Med Chem. (2012) 55:5121–9. doi: 10.1021/jm300014q
30. Epstein JI, Allsbrook WCJ, Amin MB, Egevad LL, Committee ISUP Grading. The 2005 international society of urological pathology (ISUP) consensus conference on gleason grading of prostatic carcinoma. Am J Surg Pathology. (2005) 29:1228–42. doi: 10.1097/01.pas.0000173646.99337.b1
31. Huang Y, Wang HC, Zhao J, Wu MH, Shih TC. Immunosuppressive roles of galectin-1 in the tumor microenvironment. Biomolecules. (2021) 11. doi: 10.3390/biom11101398
32. Nguyen KB, Spranger S. Modulation of the immune microenvironment by tumor-intrinsic oncogenic signaling. J Cell Biol. (2020) 219. doi: 10.1083/jcb.201908224
33. Troitskaya OS, Novak DD, Richter VA, Koval OA. Immunogenic cell death in cancer therapy. Acta naturae. (2022) 14:40–53. doi: 10.32607/actanaturae.11523
34. Chawla S, Warren TA, Wockner LF, Lambie DL, Brown IS, Martin TP, et al. Galectin-1 is associated with poor prognosis in patients with cutaneous head and neck cancer with perineural spread. Cancer immunology immunotherapy CII. (2016) 65:213–22. doi: 10.1007/s00262-015-1788-z
35. Cagnoni AJ, Giribaldi ML, Blidner AG, Cutine AM, Gatto SG, Morales RM, et al. Galectin-1 fosters an immunosuppressive microenvironment in colorectal cancer by reprogramming CD8(+) regulatory T cells. Proc Natl Acad Sci United States America. (2021) 118. doi: 10.1073/pnas.2102950118
36. Fares CM, Allen EMV, Drake CG, Allison JP, Hu-Lieskovan S. Mechanisms of resistance to immune checkpoint blockade: why does checkpoint inhibitor immunotherapy not work for all patients? Am Soc Clin Oncol Educ book. (2019) 39):147–64. doi: 10.1200/EDBK_240837
37. Zheng Y, Zhang H, Xiao C, Deng Z, Fan T, Zheng B, et al. KLF12 overcomes anti-PD-1 resistance by reducing galectin-1 in cancer cells. J immunotherapy Cancer. (2023) 11. doi: 10.1136/jitc-2023-007286
38. Nambiar DK, Aguilera T, Cao H, Kwok S, Kong C, Bloomstein J, et al. Galectin-1-driven T cell exclusion in the tumor endothelium promotes immunotherapy resistance. J Clin Invest. (2019) 129:5553–67. doi: 10.1172/JCI129025
39. Verschuere T, Toelen J, Maes W, Poirier F, Boon L, Tousseyn T, et al. Glioma-derived galectin-1 regulates innate and adaptive antitumor immunity. Int J Cancer. (2014) 134:873–84. doi: 10.1002/ijc.28426
40. Rudjord-Levann AM, Ye Z, Hafkenscheid L, Horn S, Wiegertjes R, Nielsen MAI, et al. Galectin-1 induces a tumor-associated macrophage phenotype and upregulates indoleamine 2,3-dioxygenase-1. iScience. (2023) 26:106984. doi: 10.1016/j.isci.2023.106984
41. Stannard KA, Collins PM, Ito K, Sullivan EM, Scott SA, Gabutero E, et al. Galectin inhibitory disaccharides promote tumour immunity in a breast cancer model. Cancer letters. (2010) 299:95–110. doi: 10.1016/j.canlet.2010.08.005
42. Blanchard H, Bum-Erdene K, Bohari MH, Yu X. Galectin-1 inhibitors and their potential therapeutic applications: a patent review. Expert Opin Ther Patents. (2016) 26:537–54. doi: 10.1517/13543776.2016.1163338
43. Sethi A, Sanam S, Alvala R, Alvala M. An updated patent review of galectin-1 and galectin-3 inhibitors and their potential therapeutic applications (2016-present). Expert Opin Ther Pat. (2021) 31:709–21. doi: 10.1080/13543776.2021.1903430
44. Laaf D, Bojarová P, Elling L, Křen V. Galectin-carbohydrate interactions in biomedicine and biotechnology. Trends Biotechnol. (2019) 37:402–15. doi: 10.1016/j.tibtech.2018.10.001
45. Koonce NA, Griffin RJ, Dings RPM. Galectin-1 inhibitor OTX008 induces tumor vessel normalization and tumor growth inhibition in human head and neck squamous cell carcinoma models. Int J Mol Sci. (2017) 18. doi: 10.3390/ijms18122671
Keywords: galectin-1, TME, immunotherapy, LLS30, prostate cancer
Citation: Wang H-C, Xia R, Chang W-H, Hsu S-W, Wu C-T, Chen C-H and Shih T-C (2024) Improving cancer immunotherapy in prostate cancer by modulating T cell function through targeting the galectin-1. Front. Immunol. 15:1372956. doi: 10.3389/fimmu.2024.1372956
Received: 18 January 2024; Accepted: 27 May 2024;
Published: 17 June 2024.
Edited by:
Zahid Pranjol, University of Sussex, United KingdomReviewed by:
Nadjia Bekkari, University of Science and Technology Houari Boumediene, AlgeriaJale Yuzugulen, Eastern Mediterranean University, Türkiye
Copyright © 2024 Wang, Xia, Chang, Hsu, Wu, Chen and Shih. 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: Tsung-Chieh Shih, dHNoaWgyQG1kYW5kZXJzb24ub3Jn
†These authors have contributed equally to this work and share first authorship