Skip to main content

ORIGINAL RESEARCH article

Front. Mol. Biosci., 31 October 2023
Sec. Lipids, Membranes and Membranous Organelles

Lipidomic landscape of lipokines in adipose tissue derived extracellular vesicles

Yan Zhang,,
&#x;Yan Zhang1,2,3*Tingyan Dong,&#x;Tingyan Dong2,3Muyao Wang,Muyao Wang1,2
  • 1Department of Oral and Maxillofacial Surgery, Tianjin Stomatological Hospital, School of Medicine, Nankai University, Tianjin, China
  • 2Key Laboratory of Bioactive Materials (Ministry of Education), State Key Laboratory of Medicinal Chemical Biology, College of Life Sciences, Nankai University, Tianjin, China
  • 3Tianjin Key Laboratory of Oral and Maxillofacial Function Reconstruction, Tianjin, China
  • 4Department of Periodontology, Hexi Subsidiary Clinical-service of Tianjin Stomatological Hospital, The Affiliated Stomatological Hospital of Nankai University, Tianjin, China

Introduction: Adipose tissue-derived extracellular vesicles (EVs-AT) are recognized as critical mediators of metabolic alterations in obesity-related diseases. However, few studies have focused on the role of lipids within EVs-AT in the development of obesity-related diseases.

Methods: In this study, we performed a targeted lipidomic analysis to compare the lipidome of EVs secreted by inguinal white adipose tissue (EVs-iWAT), epididymal white adipose tissue (EVs-eWAT), and interscapular brown adipose tissue (EVs-BAT) in lean and obese mice.

Results: We uncovered a comprehensive lipidomic map, revealing the diversity and specific lipid sorting in EVs-iWAT, EVs-eWAT, and EVs-BAT in obesity. Biological function analyses suggested that lipids encapsulated within EVs-AT of obese individuals might correlate with metabolism, pro-inflammatory response, and insulin resistance. These effects were particularly pronounced in EVs-eWAT and EVs-BAT.

Conclusion: Our findings indicated that EVs-AT serves as novel carriers for lipokines, thereby mediating the biological functions of EVs-AT. This study holds promise for the identification of new biomarkers for obesity-related diseases and the development of new strategies to combat metabolic diseases.

GRAPHICAL ABSTRACT

1 Introduction

The prevalence of obesity has increased substantially over the past four decades, imposing an enormous burden on public health (Kusminski et al., 2016). Obesity leads to chronic and systemic inflammation, contributing to various complications, including insulin resistance (IR), diabetes, non-alcoholic fatty liver disease (NAFLD), cardiovascular disease, and nephropathy (Reilly and Saltiel, 2017; Piché et al., 2020). The expansion of adipose tissue during obesity is a determinant of metabolic outcomes, prompting significant interest in adipose tissue biology (Sake et al., 2022).

Adipose tissue is commonly categorized into white, beige, and brown adipose tissues, each possessing distinct functions and phenotypes. White adipose tissue (WAT) primarily serves as an energy reservoir in the form of triglycerides. WAT is distributed across various depots and is primarily classified as subcutaneous WAT (sWAT) and visceral WAT (vWAT). The distribution and morphology of adipose tissue play a key role in determining the extent of adverse metabolic effects. Previous studies have demonstrated that sWAT exerts a protective effect on energy homeostasis (Marcelin et al., 2022). However, the dysregulation of vWAT is associated with increased risks of metabolic diseases (Cypess, 2022). Brown adipose tissue (BAT) is responsible for thermogenesis, and increasing BAT mass has been regarded as a promising strategy to combat obesity and related diseases. However, a recent study has shown that BAT from obese and hyperglycemic mice exhibits higher levels of inflammation (macrophages and T cell infiltration), endoplasmic reticulum (ER) stress, oxidative damage, and enhanced mitochondrial respiration activity (Villarroya et al., 2018; Carpentier et al., 2023). These studies highlight the involvement of various adipose tissue types in the development of systemic inflammation.

Adipose tissue secretes a multitude of protein factors, known as adipokines, to regulate the distal target organs (Scheja and Heeren, 2019). These adipokines have pro-inflammatory or anti-inflammatory properties, contributing to the development of systemic inflammation (Ouchi et al., 2011). Beyond polypeptide adipokines, the adipose tissue secretome comprises several chemically distinct classes of bioactive molecules, including bioactive lipids (Cao et al., 2008). These lipids are actively secreted into circulation and are termed as “lipokines” (Li et al., 2020). Recent studies have proposed lipokines to function as endocrine mediators facilitating cross-talk between adipose tissue and distant organs, including the pancreas, liver, and muscle, through diverse cell-surface and intracellular mechanisms of action (Li et al., 2020). For instance, obese adipose tissue secretes lipokines that induce inflammation in preadipocytes through a combined paracrine/autocrine manner, involving TLR4 and increased ROS, thus creating a paracrine loop that facilitates the differentiation of adipocytes with a proinflammatory profile (Renovato-Martins et al., 2020).

Extracellular vesicles (EVs) have emerged as a novel means of adipokine release and intercellular communication. However, few studies have focused on the endocrine roles of lipids within EVs-AT. Previous studies have reported that adipose tissue from lean mice releases ∼1% of its lipid content per day via EVs, a rate that more than doubles in obese animals (Flaher et al., 2019). Adipocytes release lipid-laden exosomes, delivering TAG locally to macrophages, thereby influencing their differentiation (Antonyak et al., 2019; Flaher et al., 2019). This indicates that lipokines are not only transported by EVs but also play vital roles in regulating EV-mediated functions. EVs provide a novel way for lipids to transfer and carry various lipid types (not limited to the secreted lipids). Therefore, characterization and quantification of the lipids in the EVs-AT will not only expand the list of lipokines but also provide insights into the health status of adipose tissue and the development of obesity and related diseases. To this purpose, we performed a targeted lipidomic analysis to compare the lipidome of EVs secreted by inguinal white adipose tissue (EVs-iWAT), epididymal white adipose tissue (EVs-eWAT), and interscapular brown adipose tissue (EVs-BAT) in lean and obese mice. Our work presents a comprehensive lipidomic landscape, revealing the lipid diversity and their specific sorting in EVs-iWAT, EVs-eWAT, and EVs-BAT. We also identified distinct EVs-AT lipid classes or species closely related to the obese state. Particularly, enrichment in some EV lipid species may represent novel biomarker candidates or mediators in metabolic dysfunctions associated with obesity.

2 Material and methods

2.1 Animals

The study was approved by the Research Ethics Committee of Tianjin Stomatological Hospital. Male 10-week-old wild type (WT) C57BL/6 mice and male 10-week-old obese C57BL/6 (ob/ob) mice were purchased from Nanjing Junke Bioengineering Co., Ltd. Animals were housed in a plexiglass cage (5 per cage) at a temperature (22ºC ± 3°C) and humidity (55% ± 15%). Animals were provided with food and sterile water ad libitum and kept on a 12-h light-dark cycle acclimated for 1 week before the study.

2.2 Isolation of adipose tissue-derived EVs

For this study, 5 g of adipose tissue (iWAT, eWAT, BAT) were separately collected from 8-week-old male wild-type or obese C57BL/6 mice. The tissue samples were washed extensively with sterile phosphate-buffered saline (PBS) to remove the debris and red blood cells. The tissue was finely sectioned into small pieces (1–2 mm (Piché et al., 2020)) under aseptic conditions. These explants were evenly distributed with a 5 mm spacing in three T75 flasks. After the explants adhered to the flasks (2 h), 10 mL a-MEM (without FBS) was gently added to each flask without disturbing the explants. The culture medium was collected after 24 h. Subsequently, the medium was centrifuged at 300 g for 10min, 2,000 g for 20 min, and filtered (0.22 μm filter) to remove the debris and large vesicles of cells. The resulting supernatant was further centrifuged at 100,000 g for 120 min 1 mL of the supernatant at the bottom of the ultracentrifuge tube was transferred to a sterile, and EVs isolation was performed using an exoEasy Maxi kit (QIAGEN) according to the manufacturer’s protocol. Briefly, 1 mL of the supernatant was mixed with an equal volume of “buffer XBP” and centrifuged at 500 g for 1 min at room temperature (RT) in the exoEasy spin column. Subsequently, the filter was washed with “buffer XWP” and the EVs were eluted in 400 µL elution “buffer XE”. The mix was then centrifuged at 500 g for 5 min at RT, reapplied to the filter, and centrifuged at 3,000 g for 5 min at RT. All samples were adjusted to a final volume of 500 µL with elution buffer and stored at −80°C (Supplementary Figure S1).

2.3 Extraction of EVs-AT lipids

Lipids were extracted according to the methyl-tert-butyl ether (MTBE) method. Briefly, the samples were spiked with an appropriate amount of internal lipid standards and then homogenized with 200 µL water and 240 µL methanol. Subsequently, 800 µL of MTBE was added, and the mixture was subjected to 20 min of ultrasound at 4°C followed by 30 min of settling at RT. The solution was then centrifuged at 14,000 g for 15 min at 10°C, yielding an upper organic solvent layer, which was collected and dried under nitrogen.

2.4 LC-MS/MS analysis

Reverse-phase chromatography was employed for LC separation using CSH C18 column (1.7 µm, 2.1 mm × 100 mm, Waters). The lipid extracts were re-dissolved in 200 µL 90% isopropanol/acetonitrile and centrifuged at 14,000 g for 15 min. Solvent A consisted of acetonitrile-water (6:4, v/v) with 0.1% formic acid and 0.1 Mm ammonium formate, while solvent B comprised acetonitrile isopropanol (1:9, v/v) with 0.1% formic acid and 0.1 Mm ammonium formate. The initial mobile phase was 30% solvent B at a flow rate of 300 μL/min, maintained for 2 min. Subsequently, solvent B was linearly increased to 100% in 23 min, followed by equilibration at 5% solvent B for 10 min. Mass spectra were acquired using Q-Exactive Plus in both positive and negative modes. ESI parameters were optimized and preset for all measurements as follows: Source temperature at 300°C, capillary temperature at 350°C, ion spray voltage at 3000V, S-Lens RF Level at 50%, and a scan range of m/z 200.

2.5 Data processing

Lipid species were identified using LipidSearch software version 4.2 (Thermo Scientific™) to process the raw data, peak alignment, retention time correction, and peak area extraction. LipidSearch contains data on more than 30 lipid classes, with information on more than 1,700,000 ion fragments. Adducts of +H and +NH4 were selected for positive mode searches, while -H and +CH3COO were selected for negative mode searches since ammonium acetate was used in the mobile phases. In the data extracted from LipidSearch, ion peaks with more than 50% missing values were removed from the group. After normalization and integration using the Perato scaling method, the processed data were imported into SIMPCA-P 16.1 (Umetrics, Umea, Sweden) for multivariate statistical analysis, including principal component analysis (PCA), partial least squares discriminant analysis (PLS-DA), and orthogonal partial least squares discriminant analysis (OPLS-DA). Lipids with significant differences were identified based on a combination of statistically significant thresholds of variable influence on projection (VIP>1) values obtained from the OPLS-DA model (multi-dimensional statistical analysis) and two-tailed Student’s t-test (p-value <0.05) applied to the raw data (unidimensional statistical analysis) using PASS 16 (https://www.ncss.com/software/pass/) before experiments.

2.6 Bioinformatics analysis

LINT-web (http://www.lintwebomics.info/single) is a web-based tool designed for comprehensive lipidomic dataset analysis (Li et al., 2021). The entire database was uploaded in LINT, and results of the lipid classification information, acyl-chain length analysis, and principal component analysis were exported as various plots, including stacked histograms of cumulative lipid classes, heatmaps, or volcano plots of lipid differences in EVs-iWAT, EVs-eWAT, EVs-BAT, and PCA score plots. Enrichment analysis for biological pathways was performed using LipidSig (http://chenglab.cmu.edu.tw/lipidsig) (Lin et al., 2021). The statistical significance of enrichment was assessed using the corrected p-value using the Benjamini–Hochberg method.

2.7 Statistical analysis

Each experiment was repeated at least three times. Data were expressed as mean values ± SD. Statistical analysis was performed with Student’s paired t-test. p < 0.05 was considered as statistically significant.

3 Results

3.1 Characteristics of adipose tissue-derived EVs (EVs-AT)

Firstly, we conducted a comprehensive characterization of EVs isolated from both lean and obese iWAT (EVs-iWAT), eWAT (EVs-eWAT), and BAT (EVs-BAT) using transmission electron microscopy (TEM), nanoparticle tracking analysis (NTA), and Western blot (WB) techniques. We found that obesity decreased the secretion of sEVs in all three adipose tissue types (Figure 1A). These vesicles exhibited a cup-shaped morphology (Figure 1B) with a diameter ranging from 50 to 200 nm (Figure 1C). The WB analysis demonstrated that the isolated EVs carried multiple positive EV markers, including CD63 and TSG101 (Figure 1D). The combined results from NTA, WB, and TEM analyses indicated no obvious differences in vesicle size distribution or appearance between lean or obese EVs-AT.

FIGURE 1
www.frontiersin.org

FIGURE 1. Characterization of EVs isolated from lean (wild type, WT) or obese iWAT (EVs-iWAT), eWAT (EVs-eWAT) and BAT (EVs-BAT). (A) Quantification analysis of EVs-AT in the lean or obese mice, Data represent mean ± SEM, **p < 0.01, Student’s t-test (Zhang et al., 2020); (B) Electron micrograph of EVs-iWAT, EVs-eWAT and EVs-BAT. Scale bar: 200 nm; (C) The size distribution of EVs-iWAT, EVs-eWAT and EVs-BAT as determined by NTA (D) Protein levels of CD63, TSG101 in EVs-iWAT, EVs-iWAT and EVs-BAT were detected by Western blot.

3.2 Lipidome analysis of normal/lean EVs-AT

We conducted an in-depth absolute-quantitative mass spectrometric analysis to investigate the differences in lipid composition in EVs-iWAT, EVs-eWAT, and EVs-BAT. We identified 1,697 distinct lipid molecules from adipose tissue-derived EVs (AT-EVs), further classified into 39 lipid classes (Figure 2A; Supplementary Tables S1–S3). Notably, several lipid classes, including triglycerides (TG), diacylglycerols (DG), phosphatidylcholines (PC), phosphatidylethanolamines (PE), ceramides (CER), Sphingomyelin (SM) (Table 1) were the most abundant lipid types in AT-EV (Figure 2B). Furthermore, we observed specific enrichment of zymosterol (ZyE), wax ester (WE), and sterol lipids (ST) in AT-EVs when compared to the ExoCarta database (http://www.exocarta.org/). Importantly, AT-EVs were not abundant in cholesterol ester (CE) lipids, indicating limited lipoprotein contamination and validating the robustness of our EV isolation methodology from tissue. However, variations in lipid class concentrations were evident among EVs-iWAT, EVs-eWAT, and EVs-BAT (Figure 2C).

FIGURE 2
www.frontiersin.org

FIGURE 2. Lipidome analysis of normal/lean EVs-AT. (A) Total lipid molecules and lipid subclasses identified in EVs-AT; (B) Statistical map of lipid subclasses and lipid molecules; (C) A stacked histogram of cumulative lipid classes in EVs-iWAT, EVs-eWAT and EVs-BAT; (D) Heatmaps of lipid classes differences in EVs-iWAT, EVs-eWAT and EVs-BAT; (E) Hierarchical clustered heatmaps of lipid species in EVs-iWAT, EVs-eWAT and EVs-BAT; (F) KEGG analysis of the enriched lipid classes in EVs-eWAT compared with EVs-iWAT; (G) KEGG analysis of the enriched lipid classes in EVs-BAT compared with EVs-iWAT.

TABLE 1
www.frontiersin.org

TABLE 1. Lipids abbreviation.

We observed a significantly higher total lipid concentration in EVs-iWAT and EVs-eWAT than in EVs-BAT. Furthermore, EVs-iWAT, EVs-eWAT, and EVs-BAT exhibited a similar lipid pattern. In Figure 1D, a bar chart was presented to highlight visual differences in lipid variance between EVs-iWAT, EVs-eWAT, and EVs-BAT. WAT-derived EVs were enriched in Cer, FA, Hex1Cer, ZyE, GD1a, and GD2 compared to EVs-BAT (Supplementary Figure S2A). However, several glycerophospholipids classes, such as LPC, LPE, LPG, PI, and CL, were higher in EVs-eWAT than in EVs-iWAT (Supplementary Figure S2B). Conversely, CerG3GNAc2 were specifically enriched in EVs-iWAT (Supplementary Figure S2C), while PG and Co were enriched in Evs-BAT relative to EVs-eWAT (Supplementary Figure S2D).

Our lipidomic analysis was next extended to the individual molecular species of the identified lipid subclasses. Specifically, we identified ten PG species, Cer(m38:2), sm(t42:5) as enriched lipids and six PE species and seven SM species as less abundant lipids in EVs-BAT compared with EVs-iWAT and EVs-eWAT (Figure 2D). In contrast, six PE species and five SM species were specifically enriched in EVs-eWAT. Subsequent enrichment analysis for the differential lipids revealed that the metabolism, PI3K-AKT, MAPK signaling, and autophagy pathways were the main KEGG pathways in EVs-eWAT (Figure 2E), while EVs-BAT demonstrated enrichment in metabolism and cancer-related pathways (Figure 2F). Collectively, these findings highlighted the variance in lipidomic profiles and biological functions in EVs derived from distinct adipose tissue types.

3.3 Obesity alters EVs-AT lipidome

To elucidate the effects of obesity on the lipid composition of EVs-AT, we performed a comparative lipidomic analysis between the obese and lean state on the leptin-deficient (ob/ob) mouse. We observed a significant increase in the total lipid concentration in obese EVs-iWAT. Principal component analysis (PCA) of lipidomic datasets revealed a distinct separation between WT EVs-iWAT and ob/ob EVs-iWAT (Figure 3A, Supplementary Figure S3A), highlighting substantial lipidomic changes in the ob/ob model. Prominent increases in several lipid species, including SPHP, StE, TG, FA, PE, ST, PIP, PIP2, PIP3, OAHFA, PG, LPC, LPE, and LPG were decreased in obese EVs-iWAT compared with the lean controls (Figure 3B, Supplementary Figure S3B). Then, we analyzed the subclasses based on the values of predictive variable importance (VIP) obtained from the OPLS-DA model for VIP> 1, p < 0.05 as differential lipid molecules visualized using bubble plots. The results also confirmed the significant upregulation of several TG, PE, and ST species and the downregulation of PC and CL species in obese EVs-iWAT (Figure 3C; Supplementary Table S4). Furthermore, we observed a significant decrease of carbon atoms in the acyl chain composition of several lipid classes, including CL, CerG3NAc2, and GD2, and an increase of carbon atoms in FA with 20 carbon atoms and Hex2Cer with 36–42 carbon atoms (Supplementary Figure S3C).

FIGURE 3
www.frontiersin.org

FIGURE 3. Obesity alters EVs-AT lipidome. (A) Principal component analysis (PCA) of EVs-iWAT lipidomes in lean and obese mice; (B) A stacked histogram of cumulative lipid classes in lean or obese EVs-iWAT; (C) The relative percentage difference in concentration of all quantified VIP lipid species between lean (left panel) or obese EVs-iWAT (right panel), each dot represents a lipid species, and the dot size indicates significance; (D) PCA of EVs-eWAT lipidomes in lean and obese mice; (E) A stacked histogram of cumulative lipid classes in lean or obese EVs-eWAT; (F) The relative percentage difference in concentration of all quantified VIP lipid species between lean (left panel) or obese EVs-eWAT (right panel), each dot represents a lipid species, and the dot size indicates significance; (G) PCA of EVs-BAT lipidomes in lean and obese mice; (H) A stacked histogram of cumulative lipid classes in lean or obese EVs-BAT; (I) The relative percentage difference in concentration of all quantified VIP lipid species between lean (left panel) or obese EVs-BAT (right panel), each dot represents a lipid species, and the dot size indicates significance.

We also observed notable lipidomic changes in EVs-eWAT in the ob/ob model (Figure 3D, Supplementary Figure S4A). Several lipid species, including FA, PC, PE, DG, WE and SM, were significantly increased, while PIP, PIP2, and PIP3 were decreased in obese EVs-eWAT. Interestingly, LPC and LPE, which were decreased in obese EVs-iWAT, exhibited an increase in obese EVs-eWAT (Figure 3E, Supplementary Figure S4B). The PC/PE ratio, inversely related to insulin sensitivity in adipocytes, decreased in obese EVs-eWAT but increased in obese EVs-iWAT. Further VIP lipid analysis reconfirmed the significant upregulation of several PC, PE, DG, Hex2Cer and downregulation of CerP, Cer, and OAHFA species in obese EVs-eWAT (Figure 3F, Supplementary Table S5). Furthermore, we observed a significant decrease of carbon atoms in the acyl chain composition of CO, MG, and DG, while carbon atoms of several species of LPC, LPG, and OAHFA decreased in obese EVs-eWAT (Supplementary Figure S4C).

In the EVs-BAT, PCA analysis also revealed substantial lipidomics changes in the ob/ob model (Figure 3G, Supplementary Figure S5A). Interestingly, despite the significant increase in total lipid concentration in obese EVs-eWAT, the lipidome profile in obese EVs-BAT tended to shift towards a more EVs-WAT-like composition (Figure 3H, Supplementary Figure S5B). Several lipid species, including TG, FA, PC, PE, LPE, LPG, and SM were significantly upregulated, while PIP2 and PG were downregulated in obese EVs-BAT. Furthermore, VIP lipid analysis confirmed the significant upregulation of several PC, PE, and SM species and the downregulation of PG and OAHFA species in obese EVs-BAT (Figure 3I; Supplementary Table S6). Furthermore, unlike EVs-iWAT and EVs-eWAT, carbon numbers of several lipids, including Hex1Cer, CL, Cer, FA, and PS, were significantly increased in obese EVs-BAT (Supplementary Figure S5C). These results underscore the profound impact of obesity on the composition, structure, and lipidomic profiles of EVs-AT.

3.4 Function analysis identifies altered lipid expression involved in metabolism and inflammation in obese EVs-AT

We performed an enrichment analysis to elucidate the functional implications of the differentially expressed lipids in obese EVs-AT. In EVs-iWAT, the differential lipids were associated with metabolic pathways, glycerophospholipid metabolism, sphingolipid signaling pathway, and choline metabolism. Additionally, pathways related to type II diabetes mellitus, adipocytokine signaling pathway, and adrenergic signaling in cardiomyocytes were also involved, indicating an association of the upregulated lipids with endocrine regulation (Figure 4A). Further analysis revealed an enrichment of the B cell receptor signaling pathway, T cell receptor signaling pathway, chemokine signaling pathway and PI3K-Akt signaling pathway in EVs-eWAT, indicating an association of the upregulated lipids with inflammation regulation, which was in contrast with EVs-iWAT (Figure 4B). In EVs-BAT, we found that altered lipids were involved in inflammatory mediator regulation of TRP channels, PI3K-Akt signaling, and TNF signaling pathways (Figure 4C).

FIGURE 4
www.frontiersin.org

FIGURE 4. KEGG pathway analysis of the changed lipids in obese EVs-iWAT (A), EVs-eWAT (B) and EVs-BAT (C).

We further analyzed the functions of differentially expressed lipid classes in obese EVs-AT, summarizing their involvement in the regulation of cell metabolism, inflammation, and signaling transduction. In obese EVs-iWAT, we observed an increase in lipids associated with metabolism, such as FA, TG, and Co, while most lipids involved in inflammation (LPC, LPE, PC, PE, SM, CerP) and signaling transduction (PIP, PIP2, GD1a, GD2, GM3) were decreased (Figure 5A). Moreover, the PC/PE ratio, a contributor to the development of insulin resistance, was increased (Figure 5B). We also used multivariate statistical analysis to identify the top 75 differential lipid substances, revealing 20 PC species, 7 PE species, 2 LPE species, and 1 LPC species to be downregulated in obese EVs-iWAT (Figure 5C).

FIGURE 5
www.frontiersin.org

FIGURE 5. Functional analysis of the changed lipids in obese EVs-AT. (A) A heatmap of changed lipid classes in obese EVs-AT, the fold change of lipid classes compared to the lean EVs-AT was represented by color ranges from blue to red, as indicated by the scale bar; (B) The concentration of PC and PE in the lean and obese EVs-iWAT, EVs-eWAT and EVs-BAT; (C) Hierarchical clustered heatmaps of lipid species in the lean and obese EVs-iWAT; (D) Hierarchical clustered heatmaps of lipid species in the lean and obese EVs-eWAT; (E) Hierarchical clustered heatmaps of lipid species in the lean and obese EVs-BAT.

In contrast to obese EVs-iWAT, obese EVs-eWAT exhibited a significant increase in lipids associated with inflammation, accompanied by a decrease in lipids related to signaling transduction (Figure 5A). The PC/PE ratio was significantly increased in obese EVs-eWAT, indicating a potential development of insulin resistance (Figure 5B). Furthermore, analysis of lipid substances revealed significant upregulation of 14 PC species, 11 PE species, and 10 Cer species. Furthermore, the other inflammation-related lipids, including 3 LPC species and 1 LPE species, were also upregulated (Figure 5D), collectively suggesting that obese EVs-eWAT may contribute to the development of low inflammation states in obesity.

Similar to obese EVs-eWAT, inflammation-associated lipids were significantly increased, while those linked to signal transduction decreased in obese EVs-BAT (Figure 5A). Analysis of lipid substances showed that 10 PG species were significantly decreased (Figure 5F). Since PG is predominantly located in mitochondrial and microsomal membranes, this decrease may indicate a decrease of mitochondria in EVs or an impaired metabolic homeostasis regulation. Furthermore, we also noted an increase in several inflammation-related lipids, including 9 SM species. These findings indicate that lipid species enriched in obese EVs-AT could serve as potential biomarkers or mediators of obesity-associated metabolic dysfunctions.

By generating chord plots, we can clearly observe the differences among the obese EVs-iWAT (Supplementary Figure S6A) and EVs-eWAT (Supplementary Figure S6B) and EVs-BAT (Supplementary Figure S6C). Our results suggested that three groups presented similar lipid correlation patterns. However, the three groups also had substantial patterned differences. For instance, the positive correlations between PC and PE are stronger in EVs-eWAT than EVs-iWAT and EVs-BAT. On the contrary, the positive correlations between PG and DG are stronger in EVs-iWAT and EVs-BAT than EVs-eWAT. Therefore, lipidomic correlation patterns provided more comprehensive molecular insights into assessing the influence of obesity on the regulation of metabolism and inflammation of EVs-AT.

4 Discussion

Previous studies highlighted the significance of EVs-AT as critical mediators of metabolic alterations associated with obesity-related diseases (Zhang et al., 2016; Gesmundo et al., 2021; Zhang et al., 2021). However, limited studies have focused on the role of EVs-AT lipids in the development of obesity-related diseases. Our study provides comprehensive lipid maps of EVs derived from subcutaneous WAT, visceral WAT, and BAT, revealing the alterations of specific lipid classes or species within obese EVs-AT, which may contribute to the pathogenesis of obesity-related diseases.

We found that SPH, PE, DG, TG, ZyE, Cer, and SM were the commonly enriched lipids in EVs-iWAT, EVs-eWAT, and EVs-BAT. Notably, enrichment of ZyE, Cer, and SM have been reported in other EV types, which might provide a higher EV membrane order degree and contribute to increasing the rigidity of EVs and increase EV resistance to non-ionic detergents. (Skotland et al., 2020; Su et al., 2021; Blandin et al., 2023). Moreover, we found that SPH, PE, DG, and TG were specifically enriched in EVs-AT compared to EVs derived from other tissues. The accumulation of TG and DG (associated with TG turnover) in EVs-AT indicates that EVs-AT could provide energy to recipient cells and regulate energy metabolism. Furthermore, both sphingosines (SPH) and PE play an important role in various complex biological processes such as signal transduction, endocytosis, cell metabolism, and apoptosis (van der Veen et al., 2017; Ogretmen, 2018), indicating that EVs-AT may facilitate vital cellular communication processes.

Interestingly, we observed distinct differences in lipid class composition between EVs-iWAT, EVs-eWAT, and EVs-BAT. For instance, SM, LPE, LPC, and LPG, which were closely related to inflammation, were enriched in EVs-eWAT (Zhang et al., 2018). This indicates that EVs-eWAT may mediate the regulation of inflammation through lipid transfer. However, PG, primarily located in mitochondrial and microsomal membranes, was enriched in EVs-BAT compared to the white adipose tissue-derived EVs (EVs-iWAT and EVs-eWAT) (Morita and Terada, 2015). Previous studies have reported that PG can be incorporated into cardiolipin to improve mitochondrial activity and inhibit inflammation (Chen et al., 2018). Therefore, EVs-BAT may transfer mitochondria and regulate mitochondrial activity. These results suggest that EVs-eWAT might be more active in the regulation of inflammation, while EVs-BAT might be involved in thermogenic and metabolic activities from the basal lipidomic standpoint.

Our investigation revealed that obesity-induced significant changes in the lipid profiles of EVs-AT. Firstly, we noticed an overall increase of longer lipids species in obese EVs-AT, especially in EVs-BAT. Recent studies have established that a change in acyl chain length affects the pathogenesis of various human diseases, including cancer, cardiomyopathy, and metabolic diseases (Grösch et al., 2012; Park and Park, 2015). For instance, the long-chain ceramides were increased in irritable bowel syndrome and cystic fibrosis, where it induces apoptosis (Grösch et al., 2012). Consistently, we also found an increase in the long-chain ceramides in obese EVs-eWAT and EVs-BAT. Although the chain length of the other lipids was also increased, limited studies have focused on their pathological roles in the adipose tissue or the EVs-AT. Understanding the significance of such changes in obesity could potentially provide novel markers for the diagnosis of obesity-related disease.

Moreover, we observed a marked decrease in lipid classes associated with signal transduction in obese EVs-AT, indicating impaired biological functions. Specifically, PIP2, which plays a crucial role in regulating many biological processes such as signaling transduction, intracellular trafficking, membrane dynamics, and cell-matrix adhesion, was decreased in obese EVs-iWAT and EVs-BAT compared to obese EVs-eWAT (Mandal, 2020). PIP2 is central to these processes. PIP2 binds to protein kinase AKT/PKB and is implicated in the PIP3-Akt signaling pathway (Shen et al., 2022). The PI3K/AKT signaling pathway promotes lipid biosynthesis and inhibits lipolysis in adipose tissue and liver (Li et al., 2011; Huang et al., 2018). During obesity, inactivation of the PI3K/AKT signaling pathway can inhibit lipolysis, leading to impaired glucose utilization (On et al., 2008; Yue and Lam, 2012). This process triggers the ectopic accumulation of lipids and subsequently causes insulin resistance. Moreover, decreased activation of the PI3K/AKT pathway can also reduce insulin secretion and impair β cell function (Huang et al., 2018), further aggravating insulin resistance in multiple tissues. These findings indicated that, in addition to miRNAs and adipokines, lipids within EVs-AT also mediate the development of insulin resistance (Kahn et al., 2019).

Furthermore, our study suggests that lipids may mediate EVs-AT-induced inflammation in obesity. Several glycerophospholipids were significantly increased in obese EVs-eWAT and EVs-BAT. Previous studies have shown the correlation between glycerophospholipids and inflammation, particularly in obesity (Bansal et al., 2016; Darabi and Kontush, 2016). For instance, LPC serves as a pro-inflammatory lipid and enhances the release of interleukin-1β (IL-1β), IL-6, and tumor necrosis factor-α (TNF-α) from adipocytes, enhances secretion of interferon-γ from peripheral blood mononuclear leucocytes, and increases activation of B cells (Darabi and Kontush, 2016; Liu et al., 2020). Further studies have also demonstrated an increased PC expression in obese intestinal epithelial cells. The EVs-derived PC binds to and activates AhR (Aryl Hydrocarbon Receptor), subsequently inhibiting the expression of genes essential for activation of the PI3K/AKT signaling pathway (Kumar et al., 2021). Consequently, this can negatively affect insulin signaling when the EVs are taken up by macrophages and hepatocytes, leading to inhibition of the insulin signaling pathway. We also observed a decrease in the PC/PE ratio in EVs-iWAT and EVs-BAT. Previous studies have indicated that this is critical for insulin signaling, with a reduced PC/PE ratio improving insulin signaling in hepatocytes (van der Veen et al., 2017; van der Veen et al., 2019). This indicates that the ratio of PC and PE in EVs-AT could serve as a potential marker for insulin resistance, and obese EVs-AT may alter the PC/PE ratio of targeted cells, including the adipocytes and liver cells, thereby mediating the development of insulin resistance in obesity. In recent years, there has been increased awareness of the importance of adipose tissue macrophages in obesity. As a response to specific microenvironmental stimuli, levels of proinflammatory M1-like ATM subtypes are elevated, whereas antiinflammatory M2-like macrophages are reduced. The M1-like ATMs have been shown to impair insulin signaling in both mice and humans. Lipids regulate signal transduction and gene regulation during macrophage activation. Dysregulated lipid metabolism is implicated in aberrant macrophage functions in obesity. Macrophage lipid synthesis regulates bioenergetics, phagocytosis, and inflammatory cytokine production. Considering that macrophages were important EV providers and enriched in obese AT. Therefore, we speculate that macrophage-derived EVs could be involved in obesity-associated EVs-AT lipid fingerprints. In summary, our study highlights the role of lipids in obese AT-EVs as inflammatory mediators, a major contributor to decreased insulin sensitivity.

5 Conclusion

In summary, the current study presented a comprehensive analysis of lipidome in EVs derived three different types of adipose tissue on scale of the entire lipidome analysis, providing new insight in comprehensive understanding of the heterogeneity of EVs-AT. In this study, we reported that the lipids of obese EVs-AT might also correlated with lipogenesis, inflammatory response and insulin resistance. This study implied that EVs-AT could be served as a novel carrier for lipokines and mediated the biological functions of EVs-AT. However, further research is needed to explore the lipid profiles in the lean/obese human EVs dervied from iWAT, eWAT and BAT and verified their physiological or the pathological roles. Taken together, our work will facilitate identifying new therapeutic targets and developing new strategies to combat metabolic diseases.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

Ethics statement

The animal studies were approved by the Research Ethics Committee of Tianjin Stomatological Hospital. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent was obtained from the owners for the participation of their animals in this study.

Author contributions

YZ: Conceptualization, Supervision, Writing–original draft, Writing–review and editing. TD: Methodology, Validation, Writing–original draft, Writing–review and editing. MW: Investigation, Methodology, Project administration, Validation, Writing–review and editing.

Funding

The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This research was supported by the China Postdoctoral Science Foundation (2022M721714), Tianjin Stomatological Hospital Key Project (2021KLMS07), and Tianjin Key Medical Discipline (Specialty) Construction Project (WKZKT02).

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/fmolb.2023.1281244/full#supplementary-material

References

Antonyak, M. A., Lukey, M. J., and Cerione, R. A. (2019). Lipid-filled vesicles modulate macrophages. Science 363 (6430), 931–932. doi:10.1126/science.aaw6765

PubMed Abstract | CrossRef Full Text | Google Scholar

Bansal, P., Gaur, S. N., and Arora, N. (2016). Lysophosphatidylcholine plays critical role in allergic airway disease manifestation. Sci. Rep. 6, 27430. doi:10.1038/srep27430

PubMed Abstract | CrossRef Full Text | Google Scholar

Blandin, A., Dugail, I., Hilairet, G., Ponnaiah, M., Ghesquière, V., Froger, J., et al. (2023). Lipidomic analysis of adipose-derived extracellular vesicles reveals specific EV lipid sorting informative of the obesity metabolic state. Cell Rep. 42 (3), 112169. doi:10.1016/j.celrep.2023.112169

PubMed Abstract | CrossRef Full Text | Google Scholar

Cao, H., Gerhold, K., Mayers, J. R., Wiest, M. M., Watkins, S. M., and Hotamisligil, G. S. (2008). Identification of a lipokine, a lipid hormone linking adipose tissue to systemic metabolism. Cell 134 (6), 933–944. doi:10.1016/j.cell.2008.07.048

PubMed Abstract | CrossRef Full Text | Google Scholar

Carpentier, A. C., Blondin, D. P., Haman, F., and Richard, D. (2023). Brown adipose tissue-A translational perspective. Endocr. Rev. 44 (2), 143–192. doi:10.1210/endrev/bnac015

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, W. W., Chao, Y. J., Chang, W. H., Chan, J. F., and Hsu, Y. H. (2018). Phosphatidylglycerol incorporates into cardiolipin to improve mitochondrial activity and inhibits inflammation. Sci. Rep. 8 (1), 4919. doi:10.1038/s41598-018-23190-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Cypess, A. M. (2022). Reassessing human adipose tissue. N. Engl. J. Med. 386 (8), 768–779. doi:10.1056/NEJMra2032804

PubMed Abstract | CrossRef Full Text | Google Scholar

Darabi, M., and Kontush, A. (2016). Phosphatidylserine in atherosclerosis. Curr. Opin. Lipidol. 27 (4), 414–420. doi:10.1097/MOL.0000000000000298

PubMed Abstract | CrossRef Full Text | Google Scholar

Flaherty, S., Grijalva, A., Xu, X., Ables, E., Nomani, A., and Ferrante, A. W. (2019). A lipase-independent pathway of lipid release and immune modulation by adipocytes. Sci. (New York, N.Y.) 363 (6430), 989–993. doi:10.1126/science.aaw2586

PubMed Abstract | CrossRef Full Text | Google Scholar

Gesmundo, I., Pardini, B., Gargantini, E., Gamba, G., Birolo, G., Fanciulli, A., et al. (2021). Adipocyte-derived extracellular vesicles regulate survival and function of pancreatic β cells. JCI Insight 6 (5), e141962. doi:10.1172/jci.insight.141962

PubMed Abstract | CrossRef Full Text | Google Scholar

Grösch, S., Schiffmann, S., and Geisslinger, G. (2012). Chain length-specific properties of ceramides. Prog. Lipid Res. 51 (1), 50–62. doi:10.1016/j.plipres.2011.11.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, X., Liu, G., Guo, J., and Su, Z. (2018). The PI3K/AKT pathway in obesity and type 2 diabetes. Int. J. Biol. Sci. 14 (11), 1483–1496. doi:10.7150/ijbs.27173

PubMed Abstract | CrossRef Full Text | Google Scholar

Kahn, C. R., Wang, G., and Lee, K. Y. (2019). Altered adipose tissue and adipocyte function in the pathogenesis of metabolic syndrome. J. Clin. Invest. 129 (10), 3990–4000. doi:10.1172/JCI129187

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumar, A., Sundaram, K., Mu, J., Dryden, G. W., Sriwastva, M. K., Lei, C., et al. (2021). High-fat diet-induced upregulation of exosomal phosphatidylcholine con tributes to insulin resistance. Nat. Commun. 12 (1), 213. doi:10.1038/s41467-020-20500-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Kusminski, C. M., Bickel, P. E., and Scherer, P. E. (2016). Targeting adipose tissue in the treatment of obesity-associated diabetes. Nat. Rev. Drug Discov. 15 (9), 639–660. doi:10.1038/nrd.2016.75

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, F., Song, J., Zhang, Y., Wang, S., Wang, J., et al. (2021). LINT-web: a web-based lipidomic data mining tool using intra-omic integrative correlation strategy. Small Methods 5 (9), e2100206. doi:10.1002/smtd.202100206

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, S., Ogawa, W., Emi, A., Hayashi, K., Senga, Y., Nomura, K., et al. (2011). Role of S6K1 in regulation of SREBP1c expression in the liver. Biochem. Biophys. Res. Commun. 412 (2), 197–202. doi:10.1016/j.bbrc.2011.07.038

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, V. L., Kim, J. T., and Long, J. Z. (2020). Adipose tissue lipokines: recent progress and future directions. Diabetes 69 (12), 2541–2548. doi:10.2337/dbi20-0012

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, W. J., Shen, P. C., Liu, H. C., Cho, Y. C., Hsu, M. K., Lin, I. C., et al. (2021). LipidSig: a web-based tool for lipidomic data analysis. Nucleic Acids Res. 49 (1), W336–w345. doi:10.1093/nar/gkab419

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, P., Zhu, W., Chen, C., Yan, B., Zhu, L., Chen, X., et al. (2020). The mechanisms of lysophosphatidylcholine in the development of diseases. Life Sci. 247, 117443. doi:10.1016/j.lfs.2020.117443

PubMed Abstract | CrossRef Full Text | Google Scholar

Mandal, K. (2020). Review of PIP2 in cellular signaling, functions and diseases. Int. J. Mol. Sci. 21 (21), 8342. doi:10.3390/ijms21218342

PubMed Abstract | CrossRef Full Text | Google Scholar

Marcelin, G., Gautier, E. L., and Clément, K. (2022). Adipose tissue fibrosis in obesity: etiology and challenges. Annu. Rev. Physiol. 84, 135–155. doi:10.1146/annurev-physiol-060721-092930

PubMed Abstract | CrossRef Full Text | Google Scholar

Morita, S. Y., and Terada, T. (2015). Enzymatic measurement of phosphatidylglycerol and cardiolipin in cultured cells and mitochondria. Sci. Rep. 5, 11737. doi:10.1038/srep11737

PubMed Abstract | CrossRef Full Text | Google Scholar

Ogretmen, B. (2018). Sphingolipid metabolism in cancer signalling and therapy. Nat. Rev. Cancer 18 (1), 33–50. doi:10.1038/nrc.2017.96

PubMed Abstract | CrossRef Full Text | Google Scholar

Ono, H., Pocai, A., Wang, Y., Sakoda, H., Asano, T., Backer, J. M., et al. (2008). Activation of hypothalamic S6 kinase mediates diet-induced hepatic ins ulin resistance in rats. J. Clin. Invest. 118 (8), 2959–2968. doi:10.1172/JCI34277

PubMed Abstract | CrossRef Full Text | Google Scholar

Ouchi, N., Parker, J. L., Lugus, J. J., and Walsh, K. (2011). Adipokines in inflammation and metabolic disease. Nat. Rev. Immunol. 11 (2), 85–97. doi:10.1038/nri2921

PubMed Abstract | CrossRef Full Text | Google Scholar

Park, W. J., and Park, J. W. (2015). The effect of altered sphingolipid acyl chain length on various disease models. Biol. Chem. 396 (6-7), 693–705. doi:10.1515/hsz-2014-0310

PubMed Abstract | CrossRef Full Text | Google Scholar

Piché, M. E., Tchernof, A., and Després, J. P. (2020). Obesity phenotypes, diabetes, and cardiovascular diseases. Circ. Res. 126 (11), 1477–1500. doi:10.1161/CIRCRESAHA.120.316101

PubMed Abstract | CrossRef Full Text | Google Scholar

Reilly, S. M., and Saltiel, A. R. (2017). Adapting to obesity with adipose tissue inflammation. Nat. Rev. Endocrinol. 13 (11), 633–643. doi:10.1038/nrendo.2017.90

PubMed Abstract | CrossRef Full Text | Google Scholar

Renovato-Martins, M., Moreira-Nunes, C., Atella, G. C., Barja-Fidalgo, C., and Moraes, J. A. (2020). Obese adipose tissue secretion induces inflammation in preadipocytes: role of toll-like receptor-4. Nutrients 12 (9), 2828. doi:10.3390/nu12092828

PubMed Abstract | CrossRef Full Text | Google Scholar

Sakers, A., De Siqueira, M. K., Seale, P., and Villanueva, C. J. (2022). Adipose-tissue plasticity in health and disease. Cell 185 (3), 419–446. doi:10.1016/j.cell.2021.12.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Scheja, L., and Heeren, J. (2019). The endocrine function of adipose tissues in health and cardiometabolic disease. Nat. Rev. Endocrinol. 15 (9), 507–524. doi:10.1038/s41574-019-0230-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Shen, Q., Han, Y., Wu, K., He, Y., Jiang, X., Liu, P., et al. (2022). MrgprF acts as a tumor suppressor in cutaneous melanoma by restraining PI3K/Akt signaling. Signal Transduct. Target Ther. 7 (1), 147. doi:10.1038/s41392-022-00945-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Skotland, T., Sagini, K., Sandvig, K., and Llorente, A. (2020). An emerging focus on lipids in extracellular vesicles. Adv. Drug Deliv. Rev. 159, 308–321. doi:10.1016/j.addr.2020.03.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Su, H., Rustam, Y. H., Masters, C. L., Makalic, E., McLean, C. A., Hill, A. F., et al. (2021). Characterization of brain-derived extracellular vesicle lipids in Alzheimer's disease. J. Extracell. Vesicles 10 (7), e12089. doi:10.1002/jev2.12089

PubMed Abstract | CrossRef Full Text | Google Scholar

van der Veen, J. N., Kennelly, J. P., Wan, S., Vance, J. E., Vance, D. E., and Jacobs, R. L. (2017). The critical role of phosphatidylcholine and phosphatidylethanolamine metabolism in health and disease. Biochim. Biophys. Acta Biomembr. 1859 (9), 1558–1572. doi:10.1016/j.bbamem.2017.04.006

PubMed Abstract | CrossRef Full Text | Google Scholar

van der Veen, J. N., Lingrell, S., McCloskey, N., LeBlond, N. D., Galleguillos, D., Zhao, Y. Y., et al. (2019). A role for phosphatidylcholine and phosphatidylethanolamine in hepatic insulin signaling. FASEB J. official Publ. Fed. Am. Soc ieties Exp. Biol. 33 (4), 5045–5057. doi:10.1096/fj.201802117R

CrossRef Full Text | Google Scholar

Villarroya, F., Cereijo, R., Gavaldà-Navarro, A., Villarroya, J., and Giralt, M. (2018). Inflammation of brown/beige adipose tissues in obesity and metabolic disease. J. Intern Med. 284 (5), 492–504. doi:10.1111/joim.12803

PubMed Abstract | CrossRef Full Text | Google Scholar

Yue, J. T. Y., and Lam, T. K. T. (2012). Lipid sensing and insulin resistance in the brain. Cell metab. 15 (5), 646–655. doi:10.1016/j.cmet.2012.01.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, C., Wang, K., Yang, L., Liu, R., Chu, Y., Qin, X., et al. (2018). Lipid metabolism in inflammation-related diseases. Analyst 143 (19), 4526–4536. doi:10.1039/c8an01046c

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., Yu, M., Dong, J., Wu, Y., and Tian, W. (2020). Identification of novel adipokines through proteomic profiling of small extracellular vesicles derived from adipose tissue. Journal of Proteome Research 19 (8), 3130–3142.

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., Yu, M., and Tian, W. (2016). Physiological and pathological impact of exosomes of adipose tissue. Cell Prolif. 49 (1), 3–13. doi:10.1111/cpr.12233

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., Zhang, C., Wang, J., Liu, H., and Wang, M. (2021). Bone-adipose tissue crosstalk: role of adipose tissue derived extracellular vesicles in bone diseases. J. Cell Physiol. 236 (11), 7874–7886. doi:10.1002/jcp.30414

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: extracellular vesicles, lipokines, adipose tissue, obesity, lipidome

Citation: Zhang Y, Dong T and Wang M (2023) Lipidomic landscape of lipokines in adipose tissue derived extracellular vesicles. Front. Mol. Biosci. 10:1281244. doi: 10.3389/fmolb.2023.1281244

Received: 22 August 2023; Accepted: 17 October 2023;
Published: 31 October 2023.

Edited by:

Anna Sansone, Consiglio Nazionale delle Ricerche, Italy

Reviewed by:

Sheetal Gandotra, Council of Scientific and Industrial Research (CSIR), India
Cristina Mas Bargues, University of Valencia, Spain

Copyright © 2023 Zhang, Dong and Wang. 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: Yan Zhang, dentistzy@126.com

These authors have contributed equally to this work

Disclaimer: 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.