Abstract

Warfarin, a vitamin K antagonist (VKA), is known to promote arterial calcification (AC). In the present study, we conducted a case-cohort study within the Multi-Ethnic Study of Atherosclerosis (MESA); 6655 participants were included. From MESA data, we found that AC was related to both age and vitamin K; furthermore, the score of AC increased with SASP marker including interlukin-6 (IL-6) and tumor necrosis factor alpha (TNF-α) rising. Next, a total of 79 warfarin users in our center developed significantly more calcified coronary plaques as compared to non-VKA users. We investigated the role of warfarin in phosphate-induced AC in different ages by in vitro experimental study. Furthermore, dose-time-response of warfarin was positively correlated with AC score distribution and plasma levels of the SASP maker IL-6 among years, but not among years. In addition, in vitro research suggested that warfarin treatment tended to deteriorate calcification in young VSMC at the early stage of calcification. Our results suggested that aging and warfarin-treatment were independently related to increased AC. Younger patients were more sensitive to warfarin-related AC than older patients, which was possibly due to accumulated warfarin-induced cellular senescence.

1. Introduction

Aortic calcification (AC) is a pathological condition with increasing prevalence of morbidity and mortality. AC is a process of osteoblast-like cell accumulation in the muscular layer of arteries. Enhanced stiffness of the arteries in AC might lead to severe vascular complications in the brain, heart, and kidneys. AC is a strong, independent predictor of cardiovascular disease (CVD) and cardiac adverse events. Risk factors of AC include advanced age, cigarette smoking, diabetes mellitus, hypertension, and kidney disease. In addition, cytokines and growth factors that also play a role in AC include proinflammatory cytokines (interleukin 6 (IL-6) and tumor necrosis factor-a (TNF-α)), osteoprotegerin, sclerostin, and matrix gamma carboxy glutamic acid-rich protein (MGP). Apart from the abovementioned traditional risk factors, warfarin, a vitamin K antagonist and currently the most frequently used anticoagulant in atrial fibrillation (AF) treatment, has been suggested to be associated with increased risk of AC [1]. Both in vitro and in vivo studies have demonstrated that warfarin could inactivate MGP, a vitamin K-dependent protein, and eventually lead to AC [25].

Aging is related to increased incidence, prevalence, and mortality associated with cardiovascular diseases. Accumulation of senescent cells within vascular tissues can potentially lead to aging-related biological dysfunction. Numerous studies have proved that aging is an independent risk factor of AC [6]. Besides, cellular senescence plays a pivotal role in the developmental pathophysiology of AC. It has been suggested that senescent vascular smooth muscle cells (VSMC) contribute to vascular dysfunction through induction to AC. However, the long-term effects of warfarin therapy and VSMC senescence have not been studied in detail.

It has been observed that there is a systemic increase level of secreted proteins with aging, which contain several proinflammatory cytokines, chemokines, tissue-damaging proteases, and growth factors. These secreted factors affect the microenvironment of tissue, in which they could propagate the stress response and regulate neighboring cells. These phenotypes are termed the senescence-associated secretory phenotype (SASP) [7]. The SASP is a critical intrinsic characteristic of cellular senescence resulting in a chronic low-grade state of inflammation that has been implicated in the development of several chronic diseases of aging including AC. Among SASP, cytokines and growth factors are important in the differentiation of senescent VSMC into calcified cells, and the effect of IL-6 on osteoblast differentiation and AC has been discussed previously [8]. So far, the long-term effects of warfarin treatment in AC and cellular senescence in different ages have not been studied in detail. However, recent studies showed that the association between warfarin treatment and AC might be influenced by the level of MGP which was not independently correlated with the mechanism of AC [4]. In this study, we conducted a retrospective analysis to investigate the contribution to warfarin use on AC and cellular senescence. Besides, we investigated the association between vitamin K and AC with data from the Multi-Ethnic Study of Atherosclerosis (MESA) cohort [6].

2. Materials and Methods

2.1. Study Population and Data Collection

The MESA is a multicenter, prospective study of risk factors affecting cardiovascular disease (CVD) progression. The 6655 MESA participants were enrolled from a multiethnical background, consisting of 6655 MESA participants that were enrolled from a multiethnical background, consisting of 2624 Caucasians, 1849 African Americans, 747 Chinese Americans, and 1435 Hispanics. The MESA recruited overlapping ethnic groups among field centers to minimize confounding by ethnicity by site. The design and study protocol of the MESA have been described elsewhere [9].

At the baseline examination, participants’ dietary intake was assessed using a modified 120-item Block-style food-frequency questionnaire (FFQ) which includes Chinese and Hispanic foods to accommodate the MESA population. The FFQ inquired about serving size, frequency, and duration of intake for selected foods and beverages, allowing quantification of nutrient intake from supplements. Daily vitamin K intake from foods was estimated by multiplying the reported amount consumed by its nutrient content. Although FFQ cannot provide estimates of absolute vitamin K intake, it correctly ranked people within the population according to relative vitamin K intake. Adjustment for total dietary intake improves accuracy and reduces extraneous between-person variation in vitamin K intake, which might increase precision of the estimates due to cancellation of correlated measurement errors for total energy and food of interest.

This study cohort included 79 patients (40 males and 39 females) who attended Sun Yat-sen Memorial Hospital between January 2013 and November 2018. We used the following inclusion and exclusion criteria to determine the study cohort. The inclusion criteria included patients diagnosed with nonvalvular AF and were at risk for stroke (CHA2DS2-VASc ), aged from 18 to 74 years. Patients with a history of primary parathyroid disease, chronic kidney disease (CKD), heart failure, cirrhotic hepatic disease, acute or chronic inflammatory disease, malignancy, and diabetes mellitus (DM) were excluded from the study. Patients with hypertension were classified based upon clearly documented medical history with systolic blood  mmHg or diastolic blood  mmHg. Patients’ baseline characteristics and the medication history of warfarin were extracted from the hospital information system. Data collection included demographic data, smoking history, medical history, and anthropometric measurements. The study protocol conformed to the ethical guidelines on the 1975 Declaration of Helsinki as reflected in a prior approval by the Ethics Committee of Sun Yat-sen University.

Dose-time-response models are used for the prediction of the response to conventional small drug molecules where the drug effects happen on much longer timescales, that is, months or years. In this study, dose-time-response of warfarin included two quantitative factors: dosage of warfarin (×1, mg) and postexposure time (×2, months). The dose-time-response value was calculated as ().

2.2. Collection of Laboratory Parameters

All the patients were required for at least 10-hour fasting before the blood samples were drawn. Laboratory parameters which included measurements of serum calcium, phosphorus, alkaline phosphatase (ALP), total cholesterol (TC), triglyceride (TG), high-density lipoprotein cholesterol (HDL-c), low-density lipoprotein cholesterol (LDL-c), international normalized ratio (INR), fasting glucose, and high-sensitivity C-reactive protein (hsCRP) were analyzed using a standardized and certified TBA-120 autoanalyzer (Toshiba Medical Systems, Japan) in the institutional central laboratory of the Sun Yat-sen Memorial Hospital [10, 11].

2.3. Multidetector Row Computerized Tomographic Analysis

All the patients underwent a nonenhanced thoracoabdominal CT scan to evaluate AC. Calcification was identified as a plaque of ≥1 mm2 with a density of >130 Hounsfield units and quantified by the Agatston scoring method as described previously [1214]. All imaging procedures were done on the same equipment using the same parameters. To measure AC, a single CT scan was performed using a scan collimation of 0.6 mm, slice thickness of 1 mm, reconstruction using 1 mm slice with 35 cm field of view, and normal kernel. Blood pressure and heart rate were monitored before the examination. Scans were read centrally by the SIEMENS Syngo CT Workplace, and calcification in the abdominal aorta was used for analysis. All CT scan data were analyzed at random order by two blinded and experienced investigators. According to the aortic calcium score, patients were grouped as having no detectable artery calcification (Agatston ) and mild (1-100), moderate (101-400), and severe (>400) AC [5, 15].

2.4. Echocardiography

All patients underwent standard two-dimensional transthoracic echocardiography which was carried out on Doppler echocardiography (Vivid 3®, GE Vingmed Ultrasound; Haifa, Israel) using a 2.5 MHz transducer. Echocardiographic parameters included left ventricular ejection fraction, left atrium, and left ventricular posterior wall depth [16].

2.5. ELISA

Plasma samples from AF patients were collected for IL-6 laboratory testing. Peripheral blood samples were collected between 8 and 10 am using vials with ethylene diamine tetraacetic anticoagulant (EDTA) and stored at 4°C for a maximum of 2 h. Plasma was sampled after centrifugation at 3,000 rpm for 10 min at 4°C and stored at −80°C until analysis. Commercial enzyme-linked immunosorbent assay (ELISA) kits purchased from Meimian Biotechnology (Yancheng, Jiangsu, China) were used to measure human IL-6 level in human plasma samples according to the manufacturer’s recommendations. Absorbances were measured in a Tecan Spark 10 M multimode microplate reader at 450 nm.

2.6. Cell Culture and Calcification Induction

VSMC were isolated from male Sprague–Dawley rats (3 weeks old, 40–60 g). Briefly, rats were euthanized using an intraperitoneal injection of sodium pentobarbital (60 mg/kg). The thoracic aorta was obtained, and the surrounding fat and connective tissues were discarded. It was slit longitudinally, and its lumen surface was scraped with a razor blade to remove the intima and then cut into 3–5 mm2 pieces. It was explanted lumen side down on collagen-coated culture dishes. After 7 days, tissue fragments were discarded and sprouted VSMC were collected (referred to as P0). Primary VSMC were cultured in Dulbecco’s modified Eagle’s medium (DMEM, Thermo Fisher Scientific) with 10% fetal bovine serum (FBS, Gibco) and 1% antibiotics (Gibco), at 37°C in 95% humidified air with 5% CO2. Calcification was induced by a solution to inorganic phosphate (Pi) (Na2HPO4 and NaH2PO4, pH 7.4, Sigma) added to serum supplemented-DMEM at concentrations of 3.0 mM (calcification medium).

2.7. Arterial Ring Calcification

Aortas (from the thoracic to the iliac arteries) were removed in a sterile manner from male Sprague–Dawley rats (150–200 g), and the adventitia and endothelium were carefully removed. The vessels were cut into 2–3 mm rings and placed in either calcification medium or normal culture medium at 37°C under 5% CO2 for 7 days, with medium changes once every 2 days. The arterial ring organ culture experiments were independently performed 3 times.

2.8. Alizarin Red S Staining

To observe calcium deposition, fixed VSMC or tissue sections were stained with 2% Alizarin Red S solution (Leagene) for 15 min at room temperature. Then, the samples were washed in PBS for 3 times. Samples were rinsed, and the stained calcium deposits were photographed.

2.9. Statistical Analysis

Data of the continuous variables were presented in mean and standard deviation (SD), skewed data in median and interquartile range, and categorical variables in absolute numbers and percentages. The baseline data between patients with AC and patients without AC were compared with the methods of the unpaired -test, Mann-Whitney test, chi-squared test, and Fisher’s exact test for distributed variables. Spearman’s rank correlation coefficient was used for exploring the correlation between the Agatston score and the baseline characteristics. Multivariable logistic regression analysis was used to explore the factors affecting calcification in patients with AF which was expressed in terms of odds ratio (OR) and 95% confidence interval (CI).

All statistical analyses were performed by using the software SPSS version 20.0 (SPSS, Inc., Chicago, Illinois, USA). And statistical significance was set to [17, 18]. For all statistical tests, two-sided was chosen to be the threshold of the statistical significance.

3. Results

3.1. Study Characteristics

At baseline, 51.8% of the population was male and the mean age was 61.90 years. Among the study population, 3356 (50.4%) participants had hypertension and 2058 (30.9%) had DM history. In addition, 3405 (51.1%) individuals were free from AC at baseline and AC was found in 3250 (48.9%) individuals.

6655 cases selected from the MESA were divided into two groups based on whether AC has been detected. The baseline characteristics across these 2 groups have been shown in Table 1. In the study population, participants with AC were older and had higher body mass index (BMI) and higher serum homocysteine and creatinine level than those without AC. Besides, participants with AC were more likely to be diabetic and had higher serum levels of proinflammatory markers and less accumulated vitamin K intake.

Similarly, baseline characteristics of the AF patients on warfarin are shown in Table 2. A total of 51.0% population was male, and the mean age was 64.0 years. Among the 79 patients, 40 (51.0%) participants had hypertension and 17 (21.0%) underwent statin therapy. In addition, 27 (34.2%) individuals were free from AC at baseline and AC was found in 52 (65.8%) individuals. The 79 patients were divided into the AC group and no-AC group. Participants with AC showed similar characteristics as the MESA cohort.

3.2. The Relationship among Vitamin K Intake, AC, and SASP

To establish the relationship among warfarin, senescence, and AC, 6655 cases from the MESA were analyzed. Since dose information of warfarin was not available in the MESA, accumulated vitamin K intake was analyzed to mimic the effect of warfarin. As shown in Figures 1(a) and 1(b), SASP marker IL-6 (, ) and TNF-α (, ) were both positively correlated with AC (Figures 1(a)1(b)) according to a linear regression analysis. Further, we divided participants into 4 groups depending on vitamin K intake. Statistical results showed that intake of vitamin K was negatively related to Agatston score (). Interestingly, individuals with higher vitamin K intake were more likely to show lower IL-6 () and TNF-α () level in the circulation (Figures 1(c) and 1(d)).

3.3. Comparison of Demographic and Clinical Characteristics between Absence and Presence of AC in AF Patients

Clinical data from 79 patients with NVAF were obtained for this analysis. These patients were divided as the noncalcification group () and calcification group () depending on the Agatston score of AC. Comparisons between demographic, biochemical parameters, and medication history across these two groups are shown in Table 2. The mean age of the calcification group was significantly higher than that of the noncalcification group (67 vs. 59 years, ). The calcification group also tended to have significantly higher levels of ALP and hsCRP in comparison to the noncalcification group (ALP, 77 U/L (57, 93) versus 62 U/L (47, 74), ; hsCRP, 5.26 mg/dL (1.67, 17.91) versus 1.49 mg/dL (0.88, 3.32), , Table 2). Interestingly, the dose and the duration of the warfarin treatment and INR control level showed no significant differences between the two groups (, Table 2). Moreover, no statistical significance was found in BMI, smoking, systolic blood pressure (SBP) and diastolic blood pressure (DBP), total cholesterol (TC), LDL-c, HDL-c, triglycerides, fasting glucose, and serum calcium and phosphorus (, Table 2).

3.4. Correlation between AC and Clinical, Biochemical, and Echocardiographic Parameters and Medication History in AF Patients

To better understand the relationship between clinical variables and AC in AF patients, correlation analyses were performed. The results are shown in Table 3. It was found that age, ALP, and hsCRP were positively correlated with AC (age, , ; ALP, , ; hsCRP, , ). However, no significant correlation was observed among AC and the dose of warfarin treatment, the duration of warfarin treatment, and the value of INR (dose of warfarin, , ; duration of warfarin, , ; and INR, , , Table 3).

3.5. Risk Factor Analysis of AC in AF Patients

To elucidate the independent predictors of AC in AF patients, multivariable logistic regression analysis was performed. As expected, age and ALP were independent risk factors of AC (, , Table 4). Of note, neither the dose nor the duration of the warfarin treatment was shown to be an independent risk factor of AC in AF patients (dose of warfarin, ; duration of warfarin, , Table 4).

3.6. Correlation between Warfarin Use and AC in Different Age Groups of AF Patients

To exclude the influence of age on the association between use of warfarin and AC, participants were divided into 2 groups. We classified the degree of calcification according to the Agatston score. The results showed that within the two age groups, no significant differences of AC were found in doses and durations of warfarin treatment (Figures 2(a) and 2(b)). To further explain the relationship between warfarin treatment and AC in AF patients, the dose-time-response of warfarin was introduced. Dose-time-response of warfarin was defined as the product of dose and duration of warfarin. As shown in Figure 2(c), there was no significant correlation between dose-time-response of warfarin and the Agatston score in all AF patients (, ). After dividing the patients into two age groups and ranged calcification groups according to Agatston score, there was still no significant in years (, Figure 2(e)). However, the dose-time-response of warfarin was significantly associated with increased Agatston score in the mild calcification group (Agatston -100) in years (, Figure 2(d)), which suggested that accumulating use of warfarin may promote low-grade calcification in young patients.

3.7. Correlation between Dose-Time-Response of Warfarin and SASP Marker Expression in Different Ages of AF Patients

The proinflammatory cytokine IL-6 is a proven marker of SASP. In order to explore the relationship between warfarin treatment and IL-6 concentration, other 17 plasma samples from AF patients were collected. IL-6 level was detected by an ELISA kit. As shown in Figure 3, the dose-time-response of warfarin was positively correlated with plasma IL-6 level in years but not in years. The correlation coefficient was 0.0075 and 0.0055, respectively. It is noteworthy that the basal level of IL-6 was much higher in years. These data suggested that warfarin-associated cell senescence was more sensitive in younger patients.

3.8. Warfarin Enhanced the Susceptibility to Pi-Induced Calcification

To investigate whether warfarin promotes Pi-stimulation-induced vascular calcification ex vivo, rat aorta rings were cultured with 0.65 mM Pi and different concentrations of warfarin for 7 days. Under normal phosphate and low-grade warfarin conditions, no calcification was observed in cultured aortas, and vascular smooth muscle layers of the aorta showed calcium deposition in the presence of 4.0 μM or higher concentration of warfarin (Figure 4(a)). We used VSMC in early passage (passage 5 (P5)) to mimic young cells, and continuously passage cultured VSMC were used as the replicative senescent VSMC. In both young and senescent VSMC, incubation in calcification medium with a higher concentration of Pi (3.0 mM) resulted in calcium deposition within 7 days but not in a relatively low concentration (2.0 mM). However, calcification occurred when VSMC were cultured in the medium containing 2.0 mM high phosphate and 2.0 μM warfarin. Warfarin-induced calcification was not different between the senescent and the control young VSMC in day 7. Interestingly, young VSMC showed increased calcification in day 3, which suggested that young VSMC were more sensitive to warfarin-induced calcification in the early stage (Figures 4(b) and 4(c)).

4. Discussion

The present study demonstrated that the association of warfarin use with AC differs in different age groups of patients. Specifically, warfarin adversely affects younger (<65 years) patients more than older (≥65 years) patients, and this is possibly due to the fact that warfarin-associated senescence was more sensitive in younger patients. In vitro study also revealed that young VSMC are more sensitive to warfarin-induced low-grade calcification.

In this paper, we found that SASP was positively related to AC from the MESA cohort, and initially suggested that SASP level was inversely correlated with vitamin K intake. Vitamin K has been widely considered an inhibitory factor in the progression of AC. There is accumulating data suggesting the relationship between vitamin K status and AC, including observational studies and the MESA cohort, which supports that our conclusion is solid. For example, Geleijnse et al. showed that dietary intake of vitamin K was inversely associated with AC and all-cause mortality in a cohort of 4807 [19]. In the MESA, we discovered that the accumulated intake of vitamin K was negatively related to AC, which demonstrated that warfarin may aggravate AC by downregulating the level of vitamin K in plasma. In addition, IL-6 and TNF-α levels were positively correlated with AC. When we divided individuals depending on vitamin K intake and analyzed the relationship among vitamin K, AC, and SASP, results showed that Agatston score and SASP level decreased with the increase of vitamin K intake. It suggested that warfarin may participate in SASP regulation during AC process. In further subgroup analysis, we demonstrated that age but not warfarin use was independently related to AC in patients with AF. Furthermore, we divided the AF patients into two different age groups (over or less than 65) and noted that the dose-time-response of warfarin was positively correlated with AC score and IL-6 level in younger AF patients. Interestingly, Weijs et al. reported that the intake of vitamin K antagonists would increase the levels of coronary artery calcification (CAC) [5]. Similar findings were demonstrated in CKD and diabetic patients [2023]. However, our results demonstrated that warfarin treatment did not aggravate AC in AF patients. Compared to these studies, there are several distinctions in our study. One possibility is that, in our study, all the AF patients were at low risk without other underlying diseases. However, the participants from other studies were diagnosed with other AC-associated comorbidities such as CKD and DM [22, 24, 25]. With the disorder of the homeostasis, such as the vitamin D deficiency and hyperphosphatemia, AF patients might aggravate AC without the use of warfarin [26]. Besides, previous studies mainly concentrated on CAC in AF patients [27, 28] instead of AC in this article. It has been highlighted that AC preceded the occurrence of CAC in patients on hemodialysis and systemic chronic inflammation had more deteriorative effect on the AC than CAC [29, 30]. Furthermore, the present observation was a single-center study, and all of the AF patients in our study were Asians, while in most of the previous studies, the participants were Europeans [31, 32, 33]. According to the guidelines of the management of AF in 2010 European Society of Cardiology, the level of INR should be controlled from 2.0 to 3.0. However, in China, most of the AF patients have their INR controlled in the range of 1.6 to 2.5, which is much lower than the guideline suggested [25]. In our study, the average level of the INR was 1.39. Besides, anticoagulant responses to warfarin vary among patients, based on their genotype. Gaikwad et al. demonstrated that the genotype of CYP2C9 and VKORC1 in East Asian countries was different from that seen among Europeans [30]. These two genes play an important role in warfarin dose variation and its related adverse effects. Since the reactivity, the intake, and the drug metabolism of warfarin in East Asians are much different from those in Europeans, the effects of warfarin on AC might be different in these two areas.

Our data from Chinese AF patients and the MESA showed that aging is an independent risk factor of AC. Early vascular aging, which particularly focuses on vascular inflammation, stem cells, and calcification, has been of interest [34, 35]. Cellular senescence is a program of cell cycle which is initiated in response to a variety of stimuli, such as oncogene activation (oncogene-induced senescence (OIS)) or replicative exhaustion (replicative senescence (RS)) [36]. SASP was first revealed by Vijg and Campisi and their group in 2008 [37]. They demonstrated that senescent cells secreted a myriad of factors associated with inflammation and oncogenesis. Senescent cells are characterized by a proinflammatory, oxidant SASP that occurs within a few days of senescence. Classical SASP inflammatory mediators include IL-6, IL-8, TNF-α, and monocyte chemotactic protein 1 (MCP-1). According to MESA data, IL-6 and TNF-α level was positively related to AC. The proinflammatory cytokine IL-6 is upregulated in the SASP of cells and plays a causal role in disease pathology. Kuilman et al. demonstrated that IL-6 has a crucial role in the generation of OIS since its depletion abolished oncogene-induced cellular senescence. Our statistical result (Figures 1(c) and 1(d)) revealed that IL-6 showed more relevance to vitamin K intake than TNF-α. Thus, we investigated the relationship between IL-6 level and warfarin intake. Interestingly, we found that warfarin-treatment influenced plasma IL-6 concentration in AF patients, especially in the younger (years). It is suggested that within the vascular aging, senescent cells have been associated with warfarin treatment. Therefore, warfarin may induce vascular cell senescence and likely contributes to the spread of vascular inflammation and oxidative stress via SASP. Thus, warfarin is likely to induce chronic low-grade inflammation in nonsenescent VSMC and then enhance vascular aging and calcification.

Treatment with warfarin enhanced calcification; compared with senescent VSMC, it worked sensitively in young VSMC at the early stage. Previous studies have reported that treatment with warfarin deteriorated VSMC calcification only in the young VSMC, not in the senescent VSMC [38], which is consistent with our finding. This different effect of warfarin is probably due to the loss of MGP in the senescent VSMC.

There are several limitations that should be highlighted in this study. This study is a small case-control study so a cause-and-effect relationship could not be concluded; more large-scale randomized clinical trials and basic researches are needed to corroborate the current observations. Besides, the genotype of the patients might be involved in the interaction between warfarin and AC. Further long-term large-scale randomized study is needed to corroborate the current observations.

5. Conclusions

Warfarin impacts AC more in younger patients (<65 years) than in older patients (≥65 years), and this is possibly due to the increased senescence by warfarin in this population. These findings suggest that warfarin use should be considered more cautiously in younger AF patients.

Data Availability

All data generated or analyzed during this study are included in the article. The original clinical and experimental data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare no conflict of interests.

Authors’ Contributions

Ningle Wei (NLW), Liuyi Lu (LYL), Ming Gao (MG), and Hui Huang (HH) designed the research. NLW, LYL, Huanji Zhang (HJZ), and MG performed the research. NLW, LYL, MG, Jie Chen (JC), Junhua Qi (JHQ), Zhaoyu Liu (ZYL), Jingfeng Wang (JFW), and HH interpreted results of the experiments. NLW, LYL, Sounak Ghosh, and HH drafted and revised the manuscript. HH reviewed the manuscript and made final changes. Ningle Wei, Liuyi Lu, Huanji Zhang, and Ming Gao contributed equally to this work.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (NSFC) (81870506, 81670676, and 81422011); project of Traditional Chinese Medicine Bureau of Guangdong Province (20201062); the Natural Science Foundation of Guangdong Province (2014A030313035); Basic Research Project of Shenzhen Science and Technology Innovation Committee (JCYJ20180306174648342 and JCYJ20190808102005602); the Public Health Scientific Research Project in Futian District, Shenzhen City (FTWS2019003), to Hui Huang; and NSFC (81500563) to Jie Chen.