- 1First Clinical College, Shandong University of Traditional Chinese Medicine, Jinan, Shandong, China
- 2Department of Oncology, Affiliated Hospital of Shandong University of Traditional Chinese Medicine, Jinan, Shandong, China
- 3Department of Cardiology, Longhua Hospital, Shanghai University of Traditional Chinese Medicine, Shanghai, China
- 4Department of Oncology, Longhua Hospital Affiliated to Shanghai University of Traditional Chinese Medicine, Shanghai, China
Background: The causal relationship between lipid-lowering drug (LLD) use and lung cancer risk is controversial, and the role of sphingolipid metabolism in this effect remains unclear.
Methods: Genome-wide association study data on low-density lipoprotein (LDL), apolipoprotein B (ApoB), and triglycerides (TG) were used to develop genetic instrumental variables (IVs) for LLDs. Two-step Mendelian randomization analyses were performed to examine the causal relationship between LLDs and lung cancer risk. The effects of ceramide, sphingosine-1-phosphate (S1P), and ceramidases on lung cancer risk were explored, and the proportions of the effects of LLDs on lung cancer risk mediated by sphingolipid metabolism were calculated.
Results: APOB inhibition decreased the lung cancer risk in ever-smokers via ApoB (odds ratio [OR] 0.81, 95% confidence interval [CI] 0.70–0.92, p = 0.010), LDL (OR 0.82, 95% CI 0.71–0.96, p = 0.040), and TG (OR 0.63, 95% CI 0.46–0.83, p = 0.015) reduction by 1 standard deviation (SD), decreased small-cell lung cancer (SCLC) risk via LDL reduction by 1 SD (OR 0.71, 95% CI 0.56–0.90, p = 0.016), and decreased the plasma ceramide level and increased the neutral ceramidase level. APOC3 inhibition decreased the lung adenocarcinoma (LUAD) risk (OR 0.60, 95% CI 0.43–0.84, p = 0.039) but increased SCLC risk (OR 2.18, 95% CI 1.17–4.09, p = 0.029) via ApoB reduction by 1 SD. HMGCR inhibition increased SCLC risk via ApoB reduction by 1 SD (OR 3.04, 95% CI 1.38–6.70, p = 0.014). The LPL agonist decreased SCLC risk via ApoB (OR 0.20, 95% CI 0.07–0.58, p = 0.012) and TG reduction (OR 0.58, 95% CI 0.43–0.77, p = 0.003) while increased the plasma S1P level. PCSK9 inhibition decreased the ceramide level. Neutral ceramidase mediated 8.1% and 9.5% of the reduced lung cancer risk in ever-smokers via ApoB and TG reduction by APOB inhibition, respectively, and mediated 8.7% of the reduced LUAD risk via ApoB reduction by APOC3 inhibition.
Conclusion: We elucidated the intricate interplay between LLDs, sphingolipid metabolites, and lung cancer risk. Associations of APOB, APOC3, and HMGCR inhibition and LPL agonist with distinct lung cancer risks underscore the multifaceted nature of these relationships. The observed mediation effects highlight the considerable influence of neutral ceramidase on the lung cancer risk reduction achieved by APOB and APOC3 inhibition.
1 Introduction
The incidence and mortality rates of lung cancer have consistently increased in accordance with global patterns of malignant tumor (Sharma, 2022). Despite continuous advancements and refinements in lung cancer treatment, the pursuit of effective prevention methods to mitigate the risk of developing this disease remains paramount (Kordiak et al., 2022). Observational studies have suggested a potential association between statin use and a reduced risk of lung cancer (Khurana et al., 2007; Yang et al., 2015; Kwon et al., 2019). However, a case–control study conducted within the UK population reported an increased risk of lung cancer associated with long-term statin exposure (Vinogradova et al., 2011). Furthermore, recent Mendelian randomization (MR) studies have challenged the notion of a causal relationship (Carter et al., 2020; Min et al., 2023). Investigations of lipidomic profiles revealed that the use of statins, PCSK9 inhibitors, and fenofibrate, instead of NPC1L1 inhibitors, results in a reduction in plasma ceramide and sphingomyelin levels (Ng et al., 2014; Tarasov et al., 2014; Ng et al., 2015; Croyal et al., 2018). Nonetheless, findings regarding the impact of statins on plasma sphingosine-1-phosphate (S1P) levels have exhibited inconsistency (Egom et al., 2013; Therond and Chapman, 2022). Preclinical investigations have elucidated the significant role of sphingolipid metabolism in the pathogenesis of lung cancer (Ogretmen, 2018; Meng et al., 2021; Tang et al., 2023), although the precise relationship remains unclear (Lin et al., 2022). Among the sphingolipids, ceramide and S1P have garnered particular attention within lung cancer studies (Goldkorn et al., 2013). Notably, ceramidases play a pivotal role in the conversion of ceramide to S1P (Pyne and Pyne, 2010). While a case–control study has uncovered associations between higher concentrations of plasma sphingosine-1-phosphate and ceramides and the increasing risk of lung cancer, these findings contradict the preclinical results that suggest that ceramides promote apoptosis in lung cancer cells (Ogretmen and Hannun, 2004; Carpinteiro et al., 2008; Alberg et al., 2013). Evidence on how lipid-lowering drugs (LLDs) may influence the risk of lung cancer by modulating sphingolipid metabolites remains elusive, along with the potential underlying mechanisms.
To explore the intricate relationships among LLDs, sphingolipid metabolites (including plasma ceramide, ceramidase, and S1P), and the risk of lung cancer, we employed a drug-target MR analysis (Holmes et al., 2021) to address a critical knowledge gap and provide valuable insights into potential chemoprevention in lung cancer.
2 Materials and methods
2.1 Study design
The results obtained from MR analysis, using data from genome-wide association studies (GWAS) to analyze the effects of exposure factors, closely resemble the findings of randomized controlled trials. This resemblance can be attributed to the random allocation of genetic variants during meiosis (Yarmolinsky et al., 2022). The expression and functionality of drug targets can be profoundly affected by genetic variations. Moreover, the impact of drugs can be anticipated through the genetic variability present in the genes encoding their protein targets (Chauquet et al., 2021). Similar to the previous study, instrumental variables (IVs) for LLDs were extracted from GWAS summary statistics related to low-density lipoprotein (LDL), apolipoprotein B (ApoB), and triglycerides (TG), enabling an analysis of the causal relationship between LLD exposure and the risk of developing lung cancer (Williams et al., 2020; Xiao et al., 2023). To clarify the mediating role of sphingolipid metabolites, we employed a two-step MR analysis. This involves utilizing the outcome variable from the initial MR analysis as the exposure variable for the subsequent MR analysis (Adams and Boutwell, 2021; Xu et al., 2022). Our analysis investigates the impact of LLDs on plasma sphingolipid metabolite (ceramide, S1P, and ceramidase) levels, alongside establishing a causal link between these sphingolipid metabolites and lung cancer risk. Additionally, we quantify the mediating influence of sphingolipid metabolites on the association between LLDs and lung cancer risk.
2.2 Genetic variant selection
Information about chromosomal and gene loci for the eight LLD targets (Williams et al., 2020), namely, ANGPTL3, APOB, APOC3, HMGCR, LPL, NPC1L1, PCSK9, and PPARA, was obtained from the National Center for Biotechnology Information gene database (https://www.ncbi.nlm.nih.gov/gene) (see Table 1 and Supplementary Table S1 for details).
We accessed the Medical Research Council Integrative Epidemiology Unit (IEU) OpenGWAS database (https://gwas.mrcieu.ac.uk/) and obtained GWAS summary statistics for LDL (ieu-a-300, n = 173,082), TG (ieu-a-302, n = 177,861), and ApoB (ieu-b-108, n = 439,214). These lipid traits were selected as they represent downstream substances influenced by the effectiveness of LLDs (Williams et al., 2020).
To generate IVs for LLDs, we used the LLD target genes, involving four targets of drugs for decreasing LDL (APOB, HMGCR, NPC1L1, and PCSK9), six targets of drugs for decreasing ApoB (ANGPTL3, APOB, APOC3, LPL, PPARA, and PCSK9), and five targets of drugs for decreasing TG (ANGPTL3, APOB, APOC3, PPARA, and LPL) (Xiao et al., 2023). We confined the single-nucleotide polymorphisms (SNPs) to those located within a 100-kb range surrounding the LLD target genes that exhibited genome-wide significance in their association with the lipid traits (p < 5.0 × 10−8). To maximize the IV strength for each LLD target, SNPs were allowed to be in weak linkage disequilibrium (r2 < 0.30, window size = 10,000 kb) with each other (Xie et al., 2023). Given the absence of any genetic variation in PPARA during the selection process, it has been excluded from subsequent evaluations (Li et al., 2023).
IV strength was assessed using the F-statistic, with values > 10 indicating non-weak IVs (Duan et al., 2021). Additionally, for a positive control analysis, GWAS summary statistics for coronary heart disease (CHD) obtained from CARDIoGRAM (ieu-a-8, n = 86,995) were utilized as the outcome measure (Huang et al., 2021).
2.3 Sources of GWAS summary statistics for ceramidases and lung cancer
GWAS summary statistics for plasma ceramide (GCST90025189, n = 6,057) and S1P levels (GCST90199657, n = 8,246) were obtained from GWAS Catalog (Cadby et al., 2022; Chen et al., 2023). Ceramidases are classified into acid, neutral, and alkaline ceramidases based on their pH characteristics. However, in the available datasets, there were only GWAS summary statistics for acid and neutral ceramidases (prot-a-178 and prot-a-179, n = 3,301) from the IEU OpenGWAS database (Sun et al., 2018). The GWAS summary statistics for lung cancer were sourced from the publicly available dataset published by McKay et al. (2017), which included lung cancer cases (n = 85,716) as well as subgroups of ever-smokers (n = 40,187), never-smokers (n = 9,859), lung squamous cell carcinoma (LUSC) (n = 63,053), lung adenocarcinoma (LUAD) (n = 66,756), and small-cell lung cancer (SCLC) (n = 24,108). Table 2 lists the demographic characteristics of detailed lung cancer GWAS summary statistics.
2.4 Statistical analysis
MR analysis for the effects of LLDs on the risk of lung cancer and sphingolipid metabolite levels was calculated. The clumping procedure (r2 < 0.001, window size = 10,000 kb) was performed to remove linkage disequilibrium, and SNPs that exhibited genome-wide significance (p < 5.0 × 10−8) were selected as IVs for ceramide, S1P, and acid/neutral ceramidases (Fang et al., 2023) (Supplementary Table S3).
When only one SNP was available as the IV, we utilized the Wald ratio method. Otherwise, we employed five methods: inverse variance-weighted (IVW), MR-Egger, simple mode, weighted median, and weighted mode. A reliable causal relationship was defined as follows: estimates from the IVW method or the Wald ratio method were statistically significant (p < 0.05) and were consistent with the direction of causal estimates computed by the MR-Egger method (Wang et al.).
To explore the mediating effect of sphingolipid metabolites regarding the effect of LLDs on lung cancer risk, a two-step MR analysis was performed. The proportion mediated by sphingolipid metabolites was calculated using the following formula (Chen et al., 2022; Zhao et al., 2022):
To assess the sensitivity of the MR results, we conducted MR-Egger intercept tests and calculated Cochran’s Q statistic to evaluate pleiotropy and heterogeneity. A significance level of p < 0.05 was applied to determine statistical significance. All MR analyses were performed using the “TwoSampleMR” package (version 0.5.6) in R software (version 4.1.1), and forest plots were generated using the “forestploter” package.
3 Results
3.1 Positive control analysis
We investigated the causal associations between the IVs for inhibiting the LLD targets and CHD. Notably, the effect of ANGPTL3 inhibition, whether via ApoB or TG reduction, did not show statistically significant results in the positive control analysis. Consequently, ANGPTL3 inhibition was excluded from the subsequent MR analysis. The inhibition of the remaining six LLD targets exhibited a significant causal association, with a reduced risk of CHD. Figure 1 and Supplementary Table S4 show detailed results, illustrating the magnitude and significance of these associations.
FIGURE 1. Results of positive control analysis. LDL, low-density lipoprotein; ApoB, apolipoprotein B; TG, triglycerides; SNPs, single-nucleotide polymorphisms; FDR, false discovery rate.
3.2 The effects of LLDs on lung cancer risk
Figures 2, 3 and Supplementary Tables S5–S10 show the results of our MR analysis, estimating the effects of LLDs on the risk of lung cancer (overall) and in five subgroups (lung cancer in ever-smokers, lung cancer in never-smokers, LUSC, LUAD, and SCLC).
FIGURE 3. Estimated effects of LLD on the risk of (A) lung cancer (overall); (B) lung cancer in ever-smokers; (C) lung cancer in never-smokers; (D) LUSC; (E) LUAD; and (F) SCLC. FDR-adjusted p-values with statistical significance (<0.05) are shown in bold. LDL, low-density lipoprotein; ApoB, apolipoprotein B; TG, triglycerides; LUAD, lung adenocarcinoma; LUSC, lung squamous cell carcinoma; SCLC, small-cell lung cancer; FDR, false discovery rate.
APOB inhibition significantly decreased the risk of lung cancer in ever-smokers [odds ratio (OR) 0.81, 95% confidence interval (CI) 0.70–0.92, p = 0.01] and SCLC (OR 0.70, 95% CI 0.55–0.88, p = 0.014) via ApoB reduction by 1 standard deviation (SD). Additionally, APOB inhibition decreased the risk of lung cancer in ever-smokers via LDL (OR 0.82, 95% CI 0.71–0.96, p = 0.040) and TG (OR 0.62, 95% CI 0.46–0.83, p = 0.015) reduction by 1 SD. Furthermore, APOB inhibition decreased the risk of SCLC via LDL reduction (OR 0.71, 95% CI 0.56–0.90, p = 0.016).
APOC3 inhibition decreased the risk of LUAD via ApoB reduction by 1 SD (OR 0.60, 95% CI 0.43–0.84, p = 0.039). However, it increased the risk of SCLC via ApoB reduction (OR 2.18, 95% CI 1.17–4.09, p = 0.029).
Interestingly, HMGCR inhibition increased the risk of SCLC via ApoB reduction by 1 SD (OR 3.04, 95% CI 1.38–6.70, p = 0.014). On the other hand, the LPL agonist decreased the risk of SCLC via ApoB (OR 0.20, 95% CI 0.07–0.58, p = 0.012) and TG (OR 0.58, 95% CI 0.43–0.77, p = 0.003) reduction by 1 SD. The MR analysis conducted exhibited neither pleiotropy nor heterogeneity.
3.3 The effects of LLDs on ceramide, S1P, and ceramidases
Figure 4 and Supplementary Tables S11–S14 show the results of MR analysis, estimating the effects of LLDs on plasma ceramide, S1P, and ceramidases levels. APOB inhibition decreased the plasma ceramide level via ApoB (OR 0.67, 95% CI 0.55–0.82, p = 5.21 × 10−4), LDL (OR 0.60, 95% CI 0.48–0.74, p = 3.76 × 10−5), and TG (OR 0.47, 95% CI 0.29–0.77, p = 0.007) reduction by 1 SD. Similarly, PCSK9 inhibition decreased the plasma ceramide level via ApoB (OR 0.59, 95% CI 0.43–0.80, p = 0.003) and LDL (OR 0.76, 95% CI 0.60–0.98, p = 0.061) reduction by 1 SD. Conversely, the LPL agonist increased the plasma ceramide level via ApoB (OR 3.45, 95% CI 1.23–9.66, p = 0.044) reduction by 1 SD and increased the plasma S1P level via ApoB (OR 2.71, 95% CI 1.63–4.51, p = 7.03 × 10−4) and TG (OR 1.34, 95% CI 1.19–1.52, p = 2.88 × 10−5) reduction by 1 SD. Interestingly, APOB inhibition increased the plasma neutral ceramidase level via ApoB (OR 1.27, 95% CI 1.07–1.51, p = 0.03) and TG (OR 1.90, 95% CI 1.17–3.09, p = 0.034) reduction by 1 SD. Although APOB inhibition increased the plasma neutral ceramidase level via LDL reduction, it did not reach statistical significance after correction (OR 1.25, 95% CI 1.03–1.51, p = 0.072). However, LLDs did not significantly impact plasma acid ceramidase levels. The MR analysis exhibited neither pleiotropy nor heterogeneity.
FIGURE 4. Estimated effects of LLD on plasma (A) ceramide; (B) S1P; (C) neutral ceramidase; and (D) acid ceramidase levels. FDR-adjusted p-values with statistical significance (<0.05) are shown in bold. LDL, low-density lipoprotein; ApoB, apolipoprotein B; TG, triglycerides; FDR, false discovery rate; S1P, sphingosine-1-phosphate.
3.4 The effects of sphingolipid metabolites on lung cancer risk
Our MR analysis revealed that an increase in the plasma neutral ceramidase level by 1 SD decreased the risk of overall lung cancer (OR 0.95, 95% CI 0.91–0.99, p = 0.032), lung cancer in ever-smokers (OR 0.92, 95% CI 0.88–0.98, p = 0.024), LUAD (OR 0.93, 95% CI 0.88–0.99, p = 0.043), and LUSC (OR 0.93, 95% CI 0.86–0.99, p = 0.041) (Figure 5A; Supplementary Table S15). However, there was no evidence of a causal relationship between the plasma acid ceramidase level and risk of lung cancer (Figure 5B; Supplementary Table S16). The plasma ceramide levels exhibited a positive correlation with heightened lung cancer risk in never-smokers; however, this relationship did not achieve statistical significance following correction (OR 1.49, 95% CI 1.10–2.01, p = 0.054) (Figure 5C; Supplementary Table S17). As the IVs of S1P did not intersect with the SNPs present in the GWAS data for lung cancer, the execution of this MR analysis was rendered unfeasible. The MR analysis exhibited neither pleiotropy nor heterogeneity.
FIGURE 5. Estimated effects of plasma (A) neutral ceramidase; (B) acid ceramidase; and (C) ceramide levels on lung cancer risk. FDR-adjusted p-values with statistical significance (<0.05) are shown in bold. LUAD, lung adenocarcinoma; LUSC, lung squamous cell carcinoma; SCLC, small-cell lung cancer; FDR, false discovery rate.
3.5 Analysis of the mediating effects of sphingolipid metabolites
The aforementioned results suggested that neutral ceramidase plays a mediating role in the effect of LLDs on lung cancer risk. Table 3 shows the proportions mediated by neutral ceramidase regarding the negative effect of APOB inhibition on lung cancer risk in smokers via ApoB and TG reduction (8.1% and 9.5%, respectively) and the negative effect of APOC3 inhibition on LUAD risk via ApoB reduction (8.7%). Figure 6 illustrates schematic diagrams of the direct and indirect effects of LLDs on lung cancer and the mediation by neutral ceramidase.
TABLE 3. Proportions mediated by neutral ceramidase regarding the effects of LLD on lung cancer risk.
FIGURE 6. Overview of the direct and indirect effects of LLD on lung cancer risk, and the effects of neutral ceramidase. The effects of LLD on lung cancer risk involving mediating effects of neutral ceramidase are highlighted in red. LDL, low-density lipoprotein; ApoB, apolipoprotein B; TG, triglycerides.
4 Discussion
In this drug-target MR analysis, we assessed the effects of LLDs on lung cancer risk via six LLD targets (APOB, APOC3, HMGCR, LPL, NPC1L1, and PCSK9), comprising 12 lipid-lowering pathways. Our findings revealed that APOB inhibition decreased the lung cancer risk in ever-smokers and SCLC. APOC3 inhibition decreased the LUAD risk but increased the SCLC risk. Additionally, HMGCR inhibition increased the risk of SCLC, whereas the LPL agonist decreased the risk of SCLC.
We delved deeper into the mediating mechanisms involving sphingolipid metabolites within these processes. Our investigation revealed that inhibiting APOB, APOC3, and HMGCR increased plasma neutral ceramidase levels. Notably, we observed that neutral ceramidase played a protective role in various forms of lung cancer, including overall lung cancer, lung cancer in ever-smokers, LUAD, and LUSC. These findings highlight the potential therapeutic significance of neutral ceramidase as a mitigating factor in the development and progression of various types of lung cancer. In addition, we found that APOB and PCSK9 inhibition led to a reduction in plasma ceramide levels, while the LPL agonist increased plasma ceramide and S1P levels. Nonetheless, no substantiated evidence has emerged to establish a causal link between ceramides and the risk of lung cancer.
Several MR studies have explored the potential preventive effects of LLDs on malignancies (Carter et al., 2020; Liu et al., 2021a; Min et al., 2023); however, these investigations have primarily focused on statins and have not included stratified analysis of different subtypes of lung cancer. In contrast to previous MR studies, our study expands the range of LLD targets under consideration. In our positive control analysis, we observed that ANGPTL3 inhibition did not achieve statistical significance. This finding aligns with the results reported by Wang et al. (2021). We made the decision not to conduct further analysis for ANGPTL3 inhibition in subsequent MR analysis.
In contrast to previous MR analysis of the effects of statins on lung cancer risk, our study reveals novel findings that HMGCR inhibition increases the risk of SCLC, as does APOC3 inhibition. This aligns with the findings reported by Vinogradova et al., where the link between statins and heightened lung cancer risk persisted even after adjusting for cardiovascular factors. However, their investigation lacked a subgroup analysis focusing on lung cancer types (Vinogradova et al., 2011). Moreover, evidence hints at the potential for long-term statin use to heighten cancer risk in women (Friedman et al., 2008). The mechanistic underpinnings behind the potential of statins in increasing SCLC risk remain elusive. We posited that this effect might stem from off-target consequences associated with prolonged statin use (Liu et al., 2021b; Jiang et al., 2023).
Immunohistochemical assessments of lung cancer tissue samples have unveiled lower protein expression of APOC3 in SCLC compared to normal lung tissue. This discrepancy might propose a pathway through which APOC3 inhibitors heighten the vulnerability to SCLC (Shi et al., 2016). Furthermore, an investigation under the UK’s Early Access to Medicines Scheme reported instances of lung cancer metastasis after the prolonged utilization of APOC3 inhibitors (volanesorsen) (Jones et al., 2023). Our findings strongly advocate for heightened vigilance concerning the potential for APOC3 inhibitors to increase lung cancer risk.
While alipogene tiparvovec acts as an LPL agonist, its limited market presence and high cost (Senior, 2017) have left the landscape of observational inquiries into its long-term cancer risk ambiguous. Nonetheless, our study propounds a shielding effect of LPL agonists against SCLC susceptibility. Moreover, Cerne et al. (2007) reported diminished LPL gene expression within lung cancer tissues as opposed to normal tissues, potentially reinforcing the notion of LPL agonists mitigating lung cancer risk.
The observed causal relationships between APOB inhibition and reduced lung cancer risk in ever-smokers and LUAD risk are particularly noteworthy, as neutral ceramidase plays a significant mediating role in these relationships. In recent years, there has been increasing research into the role and mechanisms of neutral ceramidase in malignancies (Coant and Hannun, 2019). Preclinical studies have demonstrated that inhibiting neutral ceramidase can prevent the occurrence and progression of colon cancer (Garcia-Barros et al., 2016), whereas inhibiting acid ceramidase can inhibit the proliferation of non-small-cell lung cancer and enhance its sensitivity to cisplatin (Yildiz-Ozer et al., 2018; White-Gilbertson et al., 2019). However, observational studies have presented contrasting findings, with high levels of acid ceramidase being associated with improved prognosis in breast and ovarian cancers (Ruckhaberle et al., 2009; Hanker et al., 2013; Sanger et al., 2015). The conflicting results from clinical observational studies and preclinical studies suggest that ceramidases may have pleiotropic effects on the occurrence, treatment, and prognosis of cancers. These findings highlight the complex and multifaceted nature of ceramidases in cancer-related processes. In our MR analysis, we identified a causal relationship between increased neutral ceramidase and a reduced risk of lung cancer in ever-smokers. This suggests that neutral ceramidase may play a critical role in the development of smoking-induced lung cancer, with APOB inhibition potentially acting as a protective factor. Further research is warranted to elucidate the precise mechanisms underlying their effects and to reconcile the disparities observed among different types of cancer.
Ceramidases play essential regulatory roles in the onset and advancement of various cancer forms (Sanger et al., 2015; Garcia-Barros et al., 2016) and are pivotal in sphingolipid metabolism for converting ceramides into sphingosine, S1P, and fatty acids (Parveen et al., 2019). The sphingolipid metabolism pathway proves indispensable in the metastatic process of lung cancer (Pyne et al., 2018; Coant and Hannun, 2019); its potential as a novel therapeutic target holds significant promise for augmenting the efficacy of tumor treatments (Vijayan et al., 2019), and animal experiments have demonstrated that NPC1L1 inhibition can modulate sphingolipid metabolism (Yamanashi et al., 2020). Although observational studies have indicated an association between statin use and reduced plasma ceramide levels, there is limited information on the impact on ceramidase levels.
In alignment with the findings of a previous study, our study identifies increased plasma ceramide levels as a contributing factor to heightened lung cancer risk among non-smokers, albeit with nominal statistical significance (Alberg et al., 2013). Preclinical studies have shown that exposure to cigarette smoke leads to increased ceramide levels in the lung tissue of mice (Filosto et al., 2011; Goldkorn et al., 2014; Lavrynenko et al., 2020); the contribution of ceramides to the heightened risk of smoking-related lung cancer necessitates prospective investigations for conclusive evidence. Consistent with previous observational studies, we found that PCSK9 inhibition reduced plasma ceramide levels (Tarasov et al., 2014), and APOB inhibition had the same effect. In addition, our results revealed that the inhibition of LPL could increase plasma S1P levels. However, owing to the non-convergence of the IVs for S1P in the lung cancer GWAS dataset, an analysis of the association between S1P and lung cancer risk was precluded. Notably, Alberg et al. (2013) observed an increased risk of lung cancer with increased S1P levels.
The strength of our study resides in its augmentation of LLD targets and lipid-lowering pathways. It amplifies the scrutiny of distinct lung cancer subgroups. More importantly, we delve into the mediating effect of sphingolipid metabolites, thus elucidating the intricate mechanism underpinning how LLD influences lung cancer risk. However, it is important to acknowledge the limitations of our study. First, MR analyses could not be performed on two of the target genes of LLD: ANGPTL3 and PPARA. This limitation arises from the inability of IVs to establish positive control associations or the unavailability of suitable genetic variants. Furthermore, given the intricate pharmacological mechanisms underlying these drugs, our MR analysis was unable to assess their potential off-target effects. Second, IVs employed as proxies for lifetime exposure to LLD may be limited in their ability to capture the effects of short-term drug exposure. Additionally, our study is constrained to examining causal directions and does not enable a precise estimation of dosage effects or the cumulative impact of multiple medications. Third, although the F-statistic for the genetic variants exceeds 10, indicating a low probability of weak instrument bias, statistical power is constrained due to the limited availability of only two SNPs as IVs for acidic and neutral ceramidases (Larsson et al., 2019). Moreover, the relatively small sample size utilized for the assessment of sphingolipids raises concerns regarding potential selection bias. As a result, the interpretation of our findings should be exercised with caution. Fourth, despite sensitivity analyses not yielding statistically significant outcomes, it remains essential to acknowledge the potential impact of confounding variables and horizontal pleiotropy, given the possibility of SNPs residing in a state of weak linkage disequilibrium (r2 < 0.30). Fifth, the absence of GWAS summary statistics for sphingolipids and lung cancer across diverse populations precludes validation. MR studies cannot substitute randomized controlled trials, and further clinical and pharmacoepidemiological studies are essential for triangulating. Sixth, the study population predominantly comprised individuals of European ancestry, which may limit the generalizability of our findings to other racial and ethnic groups.
5 Conclusion
We elucidated the intricate interplay between LLDs, sphingolipid metabolites, and lung cancer risk. Associations of APOB, APOC3, and HMGCR inhibition and the LPL agonist with distinct lung cancer risks underscore the multifaceted nature of these relationships. The observed mediation effects highlight the considerable influence of neutral ceramidase on lung cancer risk reduction achieved by APOB and APOC3 inhibition.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics statement
The studies involving humans did not require additional ethical approval as the publicly available GWAS summary datasets were collected in studies that had already obtained the necessary ethical approvals. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required from the participants or the participants’ legal guardians/next of kin.
Author contributions
HoL: conceptualization, data curation, methodology, software, visualization, writing–original draft, and writing–review and editing. LZ: conceptualization, data curation, methodology, software, visualization, writing–original draft, and writing–review and editing. FY: data curation, software, visualization, and writing–review and editing. XF: data curation and writing–review and editing. RF: data curation, software, and writing–review and editing. RZ: software and writing–review and editing. XL: writing–review and editing. HuL: conceptualization, data curation, formal analysis, investigation, supervision, visualization, and writing–review and editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. The study was funded by the Natural Science Foundation of Shandong Province (ZR2021LZY029).
Conflict of interest
The 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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2023.1269291/full#supplementary-material
References
Adams, C. D., and Boutwell, B. B. (2021). Using multiple mendelian randomization approaches and genetic correlations to understand obesity, urate, and gout. Sci. Rep. 11 (1), 17799. doi:10.1038/s41598-021-97410-4
Alberg, A. J., Armeson, K., Pierce, J. S., Bielawski, J., Bielawska, A., Visvanathan, K., et al. (2013). Plasma sphingolipids and lung cancer: a population-based, nested case-control study. Cancer Epidemiol. Biomarkers Prev. 22 (8), 1374–1382. doi:10.1158/1055-9965.EPI-12-1424
Cadby, G., Giles, C., Melton, P. E., Huynh, K., Mellett, N. A., Duong, T., et al. (2022). Comprehensive genetic analysis of the human lipidome identifies loci associated with lipid homeostasis with links to coronary artery disease. Nat. Commun. 13 s(1), 3124. doi:10.1038/s41467-022-30875-7
Carpinteiro, A., Dumitru, C., Schenck, M., and Gulbins, E. (2008). Ceramide-induced cell death in malignant cells. Cancer Lett. 264 (1), 1–10. doi:10.1016/j.canlet.2008.02.020
Carter, P., Vithayathil, M., Kar, S., Potluri, R., Mason, A. M., Larsson, S. C., et al. (2020). Predicting the effect of statins on cancer risk using genetic variants from a mendelian randomization study in the UK biobank. Elife 9, e57191. doi:10.7554/eLife.57191
Cerne, D., Melkic, E., Trost, Z., Sok, M., and Marc, J. (2007). Lipoprotein lipase activity and gene expression in lung cancer and in adjacent noncancer lung tissue. Exp. Lung Res. 33 (5), 217–225. doi:10.1080/01902140701481054
Chauquet, S., Zhu, Z., O'Donovan, M. C., Walters, J. T. R., Wray, N. R., and Shah, S. (2021). Association of antihypertensive drug target genes with psychiatric disorders: a mendelian randomization study. JAMA Psychiatry 78 (6), 623–631. doi:10.1001/jamapsychiatry.2021.0005
Chen, J., Zhou, H., Jin, H., and Liu, K. (2022). Role of inflammatory factors in mediating the effect of lipids on nonalcoholic fatty liver disease: a two-step, multivariable mendelian randomization study. Nutrients 14 (20), 4434. doi:10.3390/nu14204434
Chen, Y., Lu, T., Pettersson-Kymmer, U., Stewart, I. D., Butler-Laporte, G., Nakanishi, T., et al. (2023). Genomic atlas of the plasma metabolome prioritizes metabolites implicated in human diseases. Nat. Genet. 55 (1), 44–53. doi:10.1038/s41588-022-01270-1
Coant, N., and Hannun, Y. A. (2019). Neutral ceramidase: advances in mechanisms, cell regulation, and roles in cancer. Adv. Biol. Regul. 71, 141–146. doi:10.1016/j.jbior.2018.10.005
Croyal, M., Kaabia, Z., Leon, L., Ramin-Mangata, S., Baty, T., Fall, F., et al. (2018). Fenofibrate decreases plasma ceramide in type 2 diabetes patients: a novel marker of cvd? Diabetes Metab. 44 (2), 143–149. doi:10.1016/j.diabet.2017.04.003
Duan, C., Shi, J., Yuan, G., Shou, X., Chen, T., Zhu, X., et al. (2021). Causal association between heart failure and alzheimer's disease: a two-sample bidirectional mendelian randomization study. Front. Genet. 12, 772343. doi:10.3389/fgene.2021.772343
Egom, E. E., Rose, R. A., Neyses, L., Soran, H., Cleland, J. G., and Mamas, M. A. (2013). Activation of sphingosine-1-phosphate signalling as a potential underlying mechanism of the pleiotropic effects of statin therapy. Crit. Rev. Clin. Lab. Sci. 50 (3), 79–89. doi:10.3109/10408363.2013.813013
Fang, S., Yarmolinsky, J., Gill, D., Bull, C. J., Perks, C. M., Consortium, P., et al. (2023). Association between genetically proxied Pcsk9 inhibition and prostate cancer risk: a mendelian randomisation study. PLoS Med. 20 (1), e1003988. doi:10.1371/journal.pmed.1003988
Filosto, S., Castillo, S., Danielson, A., Franzi, L., Khan, E., Kenyon, N., et al. (2011). Neutral sphingomyelinase 2: a novel target in cigarette smoke-induced apoptosis and lung injury. Am. J. Respir. Cell Mol. Biol. 44 (3), 350–360. doi:10.1165/rcmb.2009-0422OC
Friedman, G. D., Flick, E. D., Udaltsova, N., Chan, J., Quesenberry, C. P., and Habel, L. A. (2008). Screening statins for possible carcinogenic risk: up to 9 Years of follow-up of 361,859 recipients. Pharmacoepidemiol Drug Saf. 17 (1), 27–36. doi:10.1002/pds.1507
Garcia-Barros, M., Coant, N., Kawamori, T., Wada, M., Snider, A. J., Truman, J. P., et al. (2016). Role of neutral ceramidase in colon cancer. FASEB J. 30 (12), 4159–4171. doi:10.1096/fj.201600611R
Goldkorn, T., Chung, S., and Filosto, S. (2013). Lung cancer and lung injury: the dual role of ceramide. Handb. Exp. Pharmacol. (216), 93–113. doi:10.1007/978-3-7091-1511-4_5
Goldkorn, T., Filosto, S., and Chung, S. (2014). Lung injury and lung cancer caused by cigarette smoke-induced oxidative stress: molecular mechanisms and therapeutic opportunities involving the ceramide-generating machinery and epidermal growth factor receptor. Antioxid. Redox Signal 21 (15), 2149–2174. doi:10.1089/ars.2013.5469
Hanker, L. C., Karn, T., Holtrich, U., Gatje, R., Rody, A., Heinrich, T., et al. (2013). Acid ceramidase (Ac)--a key enzyme of sphingolipid metabolism--correlates with better prognosis in epithelial ovarian cancer. Int. J. Gynecol. Pathol. 32 (3), 249–257. doi:10.1097/PGP.0b013e3182673982
Holmes, M. V., Richardson, T. G., Ference, B. A., Davies, N. M., and Davey Smith, G. (2021). Integrating Genomics with biomarkers and therapeutic targets to invigorate cardiovascular drug development. Nat. Rev. Cardiol. 18 (6), 435–453. doi:10.1038/s41569-020-00493-1
Huang, W., Xiao, J., Ji, J., and Chen, L. (2021). Association of lipid-lowering drugs with covid-19 outcomes from a mendelian randomization study. Elife 10, e73873. doi:10.7554/eLife.73873
Jiang, J. C., Hu, C., McIntosh, A. M., and Shah, S. (2023). Investigating the potential anti-depressive mechanisms of statins: a transcriptomic and mendelian randomization analysis. Transl. Psychiatry 13 (1), 110. doi:10.1038/s41398-023-02403-8
Jones, A., Peers, K., Wierzbicki, A. S., Ramachandran, R., Mansfield, M., Dawson, C., et al. (2023). Long-term effects of volanesorsen on triglycerides and pancreatitis in patients with familial chylomicronaemia syndrome (fcs) in the UK early access to Medicines Scheme (eams). Atherosclerosis 375, 67–74. doi:10.1016/j.atherosclerosis.2023.05.008
Khurana, V., Bejjanki, H. R., Caldito, G., and Owens, M. W. (2007). Statins reduce the risk of lung cancer in humans: a large case-control study of us veterans. Chest 131 (5), 1282–1288. doi:10.1378/chest.06-0931
Kordiak, J., Bielec, F., Jablonski, S., and Pastuszak-Lewandoska, D. (2022). Role of beta-carotene in lung cancer primary chemoprevention: a systematic review with meta-analysis and meta-regression. Nutrients 14 (7), 1361. doi:10.3390/nu14071361
Kwon, Y. J., You, N. Y., Lee, J. W., Kim, J., and Kang, H. T. (2019). High receipt of statins reduces the risk of lung cancer in current smokers with hypercholesterolemia: the national health insurance service-health screening cohort. Clin. Lung Cancer 20 (2), e177–e185. doi:10.1016/j.cllc.2018.11.009
Larsson, S. C., Traylor, M., and Markus, H. S. (2019). Homocysteine and small vessel stroke: a mendelian randomization analysis. Ann. Neurol. 85 (4), 495–501. doi:10.1002/ana.25440
Lavrynenko, O., Titz, B., Dijon, S., Santos, D. D., Nury, C., Schneider, T., et al. (2020). Ceramide ratios are affected by cigarette smoke but not heat-not-burn or E-vapor aerosols across four independent mouse studies. Life Sci. 263, 118753. doi:10.1016/j.lfs.2020.118753
Li, Z., Zhang, B., Liu, Q., Tao, Z., Ding, L., Guo, B., et al. (2023). Genetic association of lipids and lipid-lowering drug target genes with non-alcoholic fatty liver disease. EBioMedicine 90, 104543. doi:10.1016/j.ebiom.2023.104543
Lin, M., Li, Y., Wang, S., Cao, B., Li, C., and Li, G. (2022). Sphingolipid metabolism and signaling in lung cancer: a potential therapeutic target. J. Oncol. 2022, 9099612. doi:10.1155/2022/9099612
Liu, G., Shi, M., Mosley, J. D., Weng, C., Zhang, Y., Lee, M. T. M., et al. (2021b). A mendelian randomization approach using 3-Hmg-Coenzyme-a reductase gene variation to evaluate the association of statin-induced low-density lipoprotein cholesterol lowering with noncardiovascular disease phenotypes. JAMA Netw. Open 4 (6), e2112820. doi:10.1001/jamanetworkopen.2021.12820
Liu, L., Sheng, C., Lyu, Z., Dai, H., and Chen, K. (2021a). Association between genetically proxied lipid-lowering drug targets and renal cell carcinoma: a mendelian randomization study. Front. Nutr. 8, 755834. doi:10.3389/fnut.2021.755834
Liu, N., Wang, G., Liu, C., Liu, J., Huang, S., Zhou, Y., et al. (2023). Non-alcoholic fatty liver disease and complications in type 1 and type 2 diabetes: a mendelian randomization study. Diabetes Obes. Metab. 25 (2), 365–376. doi:10.1111/dom.14877
McKay, J. D., Hung, R. J., Han, Y., Zong, X., Carreras-Torres, R., Christiani, D. C., et al. (2017). Large-scale association analysis identifies new lung cancer susceptibility loci and heterogeneity in genetic susceptibility across histological subtypes. Nat. Genet. 49 (7), 1126–1132. doi:10.1038/ng.3892
Meng, Q., Hu, X., Zhao, X., Kong, X., Meng, Y. M., Chen, Y., et al. (2021). A circular network of coregulated sphingolipids dictates lung cancer growth and progression. EBioMedicine 66, 103301. doi:10.1016/j.ebiom.2021.103301
Min, Y., Wei, X., Liu, Z., Wei, Z., Pei, Y., Li, R., et al. (2023). Assessing the role of lipid-lowering therapy on multi-cancer prevention: a mendelian randomization study. Front. Pharmacol. 14, 1109580. doi:10.3389/fphar.2023.1109580
Ng, T. W., Ooi, E. M., Watts, G. F., Chan, D. C., Meikle, P. J., and Barrett, P. H. (2015). Association of plasma ceramides and sphingomyelin with vldl apob-100 fractional catabolic rate before and after rosuvastatin treatment. J. Clin. Endocrinol. Metab. 100 (6), 2497–2501. doi:10.1210/jc.2014-4348
Ng, T. W., Ooi, E. M., Watts, G. F., Chan, D. C., Weir, J. M., Meikle, P. J., et al. (2014). Dose-dependent effects of rosuvastatin on the plasma sphingolipidome and phospholipidome in the metabolic syndrome. J. Clin. Endocrinol. Metab. 99 (11), E2335–E2340. doi:10.1210/jc.2014-1665
Ogretmen, B. (2018). Sphingolipid metabolism in cancer signalling and therapy. Nat. Rev. Cancer 18 (1), 33–50. doi:10.1038/nrc.2017.96
Ogretmen, B., and Hannun, Y. A. (2004). Biologically active sphingolipids in cancer pathogenesis and treatment. Nat. Rev. Cancer 4 (8), 604–616. doi:10.1038/nrc1411
Parveen, F., Bender, D., Law, S. H., Mishra, V. K., Chen, C. C., and Ke, L. Y. (2019). Role of ceramidases in sphingolipid metabolism and human diseases. Cells 8 (12), 1573. doi:10.3390/cells8121573
Pyne, N. J., El Buri, A., Adams, D. R., and Pyne, S. (2018). Sphingosine 1-phosphate and cancer. Adv. Biol. Regul. 68, 97–106. doi:10.1016/j.jbior.2017.09.006
Pyne, N. J., and Pyne, S. (2010). Sphingosine 1-phosphate and cancer. Nat. Rev. Cancer 10 (7), 489–503. doi:10.1038/nrc2875
Ruckhaberle, E., Holtrich, U., Engels, K., Hanker, L., Gatje, R., Metzler, D., et al. (2009). Acid ceramidase 1 expression correlates with a better prognosis in Er-positive breast cancer. Climacteric 12 (6), 502–513. doi:10.3109/13697130902939913
Sanger, N., Ruckhaberle, E., Gyorffy, B., Engels, K., Heinrich, T., Fehm, T., et al. (2015). Acid ceramidase is associated with an improved prognosis in both dcis and invasive breast cancer. Mol. Oncol. 9 (1), 58–67. doi:10.1016/j.molonc.2014.07.016
Senior, M. (2017). After glybera's withdrawal, what's next for gene therapy? Nat. Biotechnol. 35 (6), 491–492. doi:10.1038/nbt0617-491
Sharma, R. (2022). Mapping of global, regional and national incidence, mortality and mortality-to-incidence ratio of lung cancer in 2020 and 2050. Int. J. Clin. Oncol. 27 (4), 665–675. doi:10.1007/s10147-021-02108-2
Shi, J., Yang, H., Duan, X., Li, L., Sun, L., Li, Q., et al. (2016). Apolipoproteins as differentiating and predictive markers for assessing clinical outcomes in patients with small cell lung cancer. Yonsei Med. J. 57 (3), 549–556. doi:10.3349/ymj.2016.57.3.549
Sun, B. B., Maranville, J. C., Peters, J. E., Stacey, D., Staley, J. R., Blackshaw, J., et al. (2018). Genomic atlas of the human plasma proteome. Nature 558 (7708), 73–79. doi:10.1038/s41586-018-0175-2
Tang, Z., Liang, D., Deubler, E. L., Sarnat, J. A., Chow, S. S., Diver, R., et al. (2023). Lung cancer metabolomics: a pooled analysis in the cancer prevention studies. Alexandria, Virginia: American Society of Clinical Oncology.
Tarasov, K., Ekroos, K., Suoniemi, M., Kauhanen, D., Sylvanne, T., Hurme, R., et al. (2014). Molecular lipids identify cardiovascular risk and are efficiently lowered by simvastatin and Pcsk9 deficiency. J. Clin. Endocrinol. Metab. 99 (1), E45–E52. doi:10.1210/jc.2013-2559
Therond, P., and Chapman, M. J. (2022). Sphingosine-1-Phosphate: metabolism, transport, atheroprotection and effect of statin treatment. Curr. Opin. Lipidol. 33 (3), 199–207. doi:10.1097/MOL.0000000000000825
Vijayan, Y., Lankadasari, M. B., and Harikumar, K. B. (2019). Acid ceramidase: a novel therapeutic target in cancer. Curr. Top. Med. Chem. 19 (17), 1512–1520. doi:10.2174/1568026619666190227222930
Vinogradova, Y., Coupland, C., and Hippisley-Cox, J. (2011). Exposure to statins and risk of common cancers: a series of nested case-control studies. BMC Cancer 11, 409. doi:10.1186/1471-2407-11-409
Wang, Q., Oliver-Williams, C., Raitakari, O. T., Viikari, J., Lehtimaki, T., Kahonen, M., et al. (2021). Metabolic profiling of angiopoietin-like protein 3 and 4 inhibition: a drug-target mendelian randomization analysis. Eur. Heart J. 42 (12), 1160–1169. doi:10.1093/eurheartj/ehaa972
Wang, S., Zhu, H., Pan, L., Zhang, M., Wan, X., Xu, H., et al. (2023). Systemic inflammatory regulators and risk of acute-on-chronic liver failure: a bidirectional mendelian-randomization study. Front. Cell Dev. Biol. 11, 1125233. doi:10.3389/fcell.2023.1125233
White-Gilbertson, S., Lu, P., Norris, J. S., and Voelkel-Johnson, C. (2019). Genetic and pharmacological inhibition of acid ceramidase prevents asymmetric cell division by neosis. J. Lipid Res. 60 (7), 1225–1235. doi:10.1194/jlr.M092247
Williams, D. M., Bandres-Ciga, S., Heilbron, K., Hinds, D., Noyce, A. J., andMe Research, T., et al. (2020). Evaluating lipid-lowering drug targets for Parkinson's disease prevention with mendelian randomization. Ann. Neurol. 88 (5), 1043–1047. doi:10.1002/ana.25880
Xiao, J., Ji, J., Zhang, N., Yang, X., Chen, K., Chen, L., et al. (2023). Association of genetically predicted lipid traits and lipid-modifying targets with heart failure. Eur. J. Prev. Cardiol. 30 (4), 358–366. doi:10.1093/eurjpc/zwac290
Xie, W., Li, J., Du, H., and Xia, J. (2023). Causal relationship between Pcsk9 inhibitor and autoimmune diseases: a drug target mendelian randomization study. Arthritis Res. Ther. 25 (1), 148. doi:10.1186/s13075-023-03122-7
Xu, M., Zheng, J., Hou, T., Lin, H., Wang, T., Wang, S., et al. (2022). Sglt2 inhibition, choline metabolites, and cardiometabolic diseases: a mediation mendelian randomization study. Diabetes Care 45 (11), 2718–2728. doi:10.2337/dc22-0323
Yamanashi, Y., Takada, T., Yamamoto, H., and Suzuki, H. (2020). Npc1l1 facilitates sphingomyelin absorption and regulates diet-induced production of vldl/ldl-associated S1p. Nutrients 12 (9), 2641. doi:10.3390/nu12092641
Yang, T. Y., Lin, W. M., Lin, C. L., Sung, F. C., and Kao, C. H. (2015). Correlation between use of simvastatin and lovastatin and female lung cancer risk: a nationwide case-control study. Int. J. Clin. Pract. 69 (5), 571–576. doi:10.1111/ijcp.12598
Yarmolinsky, J., Diez-Obrero, V., Richardson, T. G., Pigeyre, M., Sjaarda, J., Pare, G., et al. (2022). Genetically proxied therapeutic inhibition of antihypertensive drug targets and risk of common cancers: a mendelian randomization analysis. PLoS Med. 19 (2), e1003897. doi:10.1371/journal.pmed.1003897
Yildiz-Ozer, M., Oztopcu-Vatan, P., and Kus, G. (2018). The investigation of ceranib-2 on apoptosis and drug interaction with carboplatin in human non small cell lung cancer cells in vitro. Cytotechnology 70 (1), 387–396. doi:10.1007/s10616-017-0154-8
Keywords: lipid-lowering drug, ceramide, ceramidase, sphingosine-1-phosphate, sphingolipid metabolism, lung cancer, Mendelian randomization study
Citation: Li H, Zhang L, Yang F, Feng X, Fu R, Zhao R, Li X and Li H (2023) Lipid-lowering drugs affect lung cancer risk via sphingolipid metabolism: a drug-target Mendelian randomization study. Front. Genet. 14:1269291. doi: 10.3389/fgene.2023.1269291
Received: 01 August 2023; Accepted: 02 November 2023;
Published: 16 November 2023.
Edited by:
Jared C. Roach, Institute for Systems Biology (ISB), United StatesReviewed by:
Aleksandra Czumaj, Medical University of Gdansk, PolandYeun-po Chiang, Emory University, United States
Copyright © 2023 Li, Zhang, Yang, Feng, Fu, Zhao, Li 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: Huijie Li, MjAwOGxpaHVpamllQDE2My5jb20=
†These authors have contributed equally to this work