- 1Department of Otolaryngology-Head and Neck Surgery, First Affiliated Hospital of Guangxi Medical University, Nanning, China
- 2State Key Laboratory of Non-Food Biomass and Enzyme Technology, Guangxi Key Laboratory of Bio-refinery, National Engineering Research Center for Non-Food Biorefinery, Guangxi Biomass Engineering Technology Research Center, Guangxi Academy of Sciences, Nanning, China
- 3Key Laboratory of Early Prevention and Treatment for Regional High-Frequency Tumor, Guangxi Key Laboratory of High-Incidence-Tumor Prevention and Treatment, Ministry of Education, Guangxi Medical University, Nanning, China
- 4Guangxi Key Laboratory of Marine Natural Products and Combinatorial Biosynthesis Chemistry, Beibu Gulf Marine Research Center, Guangxi Academy of Sciences, Nanning, China
- 5Guangxi Health Commission Key Laboratory of Molecular Epidemiology of Nasopharyngeal Carcinoma, Wuzhou Red Cross Hospital, Wuzhou, China
Plasma lipidomics has been commonly used for biomarker discovery. Studies in cancer have suggested a significant alteration of circulating metabolite profiles which is correlated with cancer characteristics and treatment outcome. However, the lipidomics characteristics of nasopharyngeal carcinoma (NPC) have rarely been studied. We previously described the phenomenon of lipid droplet accumulation in NPC cells and showed that such accumulation could be regulated by latent infection of Epstein–Barr virus (EBV). Here, we compared the plasma lipidome of NPC patients to that of healthy controls by liquid chromatography-tandem mass spectrometry (LC–MS/MS). We found 19 lipids (e.g., phosphatidylinositols 18:0/20:4 and 18:0/18:2 and free fatty acid 22:6) to be remarkably decreased, whereas 2 lipids (i.e., diacylglycerols 16:0/16:1 and 16:0/20:3) to be increased, in the plasma of NPC patients, compared with controls. Different lipid profiles were also observed between patients with different titers of EBV antibodies (e.g., EA-IgA and VCA-IgA) as well as between patients with and without lymph node or distant organ metastasis. In conclusion, plasma lipidomics might help to differentiate NPC cases from controls, whereas EBV infection might influence the risk and prognosis of NPC through modulating lipid metabolism in both tumor cells and peripheral blood.
Introduction
As reported in GLOBOCAN 2020, there were 133,354 newly diagnosed cases of nasopharyngeal carcinoma (NPC) in 2020, accounting for 0.7% of the total human cancer globally (Sung et al., 2021). NPC is prevalent in Southeast Asia, especially China, and relatively rare elsewhere (Yu and Yuan, 2002; Zhang et al., 2015; Chen et al., 2019). Multiple factors contribute to the initiation and development of NPC, among which, the crosstalk of Epstein–Barr virus (EBV), genetic factors, and environmental carcinogens likely constitutes the main contributor.
Epstein–Barr virus infection is the most established risk factor for NPC, especially for non-keratinizing NPC: the most common histological type of NPC in endemic regions (Pathmanathan et al., 1995; Yu and Hussain, 2009). EBV infection may exhibit a latent or lytic phase (Hoch et al., 1988; Moss et al., 2001). In the latent phase, the EBV genomic DNA exists like host chromosomal DNA. Expression of several EBV proteins facilitates this phase, including latent membrane proteins (LMP1, LMP2A, and LMP2B) and EBV nuclear antigens (EBNA1, EBNA2, EBNA3A, EBNA3B, and EBNA3C). The expression of the immediate-early proteins, Zta and Rta (encoded by BZLF1 and BRLF1, respectively), initiates on the other hand the lytic phase of EBV infection (Robertson and Kieff, 1995; Binné et al., 2002; Hui and Chiang, 2010). EBV reactivation may be identified by detection of IgA antibody to EBV early antigen (EA) and late antigen (VCA; Hoch et al., 1988). As early as 37 months before the clinical onset of NPC, an increased serological level of EBV VCA-IgA might already be detected (Ji et al., 2007). To date, serological detection of EBV-related antibodies including VCA-IgA, EA-IgA, and EBV nuclear antigen 1 (EBNA1)-IgA antibodies has been widely applied in NPC screening programs (Zeng et al., 1982; Liu et al., 2012; Coghill et al., 2014; Gao et al., 2017).
Accumulating evidence suggests the involvement of lipid metabolism reprogramming in cancer (Luo et al., 2017). Cancer cells increase de novo synthesis to better utilize lipids in energy metabolism to facilitate malignant behaviors like cell death resistance and replicative immortality enablement (Carracedo et al., 2013). In a previous study, we showed that there was an aberrant accumulation of lipid droplets in NPC cells, likely attributable to the epigenetic downregulation of the ISG15-conjugating enzyme UbcH8, leading to altered stabilization of adipose triglyceride lipase (Zhou et al., 2015). High-density lipoprotein cholesterol has also been shown to enhance proliferation, migration, colony formation, and chemoresistance ability of NPC cells, whereas a higher level of high-density lipoprotein cholesterol was found to predict a poor survival rate of NPC patients (Liu et al., 2016).
It is known that there is an interplay between host cell lipid metabolism and infection of oncogenic viruses including EBV (Daker et al., 2013). EBV infection may therefore promote the development of NPC through modulating lipid metabolism. EBV-encoded latent membrane protein 1 may promote the proliferation and progression of NPC by activating the sterol regulatory element-binding protein 1 and its downstream molecular fatty acid synthase, leading to an increased de novo lipid synthesis and droplet formation (Lo et al., 2018). Our previous study also showed that EBV-encoded latent membrane protein 2A (LMP2A) could drive the downregulation of adipose triglyceride lipase, leading to lipid accumulation and migration of NPC cells (Zheng et al., 2020). EBV-encoded RNAs are found highly expressed in EBV-positive NPC cells (Raab-Traub, 2015). Although they do not seem to affect the migration or invasion of cancer cells, EBV-encoded RNAs might confer a low-density lipoprotein-dependent proliferation in NPC cells (Daker et al., 2013). Furthermore, EBV could target adipocytes to build a cancer-promoting environment by rewiring lipid metabolism, leading to increased levels of free fatty acids, glycerol, and proinflammatory cytokines (Liu et al., 2021). Indeed, increased serum levels of free fatty acids have been shown to indicate a poor prognosis in NPC patients (Liu et al., 2021).
Lipidomics identifies changes in lipid metabolism, transport and homeostasis and has played an important role in understanding biochemical mechanisms of lipid-related diseases (Han, 2016). Recent advances in mass spectrometry and other spectroscopic methods have further promoted the development and use of lipidomics in research (Hu et al., 2009). For instance, liquid chromatography-tandem mass spectrometry (LC–MS/MS) has been used widely in studying complex metabolomes (Leung and Fong, 2014). In terms of cancer, lipidomics has been used in studies of lung cancer (Zhang et al., 2020), bladder cancer (Shahid et al., 2020), esophageal cancer (Liu et al., 2013), colorectal cancer (Tan et al., 2013), thyroid cancer (Jiang et al., 2021), and breast cancer (Fichtali et al., 2020). However, to our best knowledge, no study has examined lipidomics in NPC, to date.
To this end, we used LC–MS/MS to determine the lipid profiles in the plasma of patients with NPC, through comparing different lipid metabolites between NPC patients and healthy controls in an endemic area of Southeast China. We also studied the correlations between lipid profiles with EBV infection, measured through EBV antibodies, and clinical characteristics of NPC, namely TMN and clinical stages.
Materials and Methods
Study Design
We recruited 100 NPC patients and 100 age- and sex-matched healthy controls in Wuzhou Red Cross Hospital (Wuzhou, China) between 2018 and 2020. The diagnosis of NPC was cross validated by two experienced pathologists according to the World Health Organization classification (Stelow and Wenig, 2017). The patients were required to have not received radiation- or chemotherapy and have no hyperlipidemia, diabetes, or other conditions affecting lipid metabolism. Patients with other malignancy or recurrent NPC were also excluded. The healthy controls were required to have no metabolic disorders and no EBV VCA-IgA and EA-IgA antibodies. We collected information on clinical characteristics of the NPC patients, including age, sex, EBV antibodies (VCA-IgA, EA-IgA, EBNA1-IgA, and Zta-IgA), as well as Tumor Node Metastasis (TNM) stage and clinical stage defined according to the eighth edition of the American Joint Committee on Cancer (AJCC) staging system.
Ethics approval of this study was obtained from the Ethical Evaluation Committee of Wuzhou Red Cross Hospital (No. LL2017-19, Wuzhou, China). All participants signed informed consent before recruitment to the study.
EBV Antibody Detection in Serum
A morning fasting peripheral blood (5 ml) sample was collected using a non-anticoagulant-treated tube. The blood samples were allowed to clot at room temperature for half an hour. The samples were then centrifuged at 2,000 rmp at 4°C for 10 min, with the supernatant serum transferred, aliquoted, and stored at –70°C for further analysis of the serological levels of VCA-IgA and EA-IgA antibodies using an immunoenzymatic assay as previously described (Cai et al., 2014). All controls had negative results on these antibodies, whereas NPC patients were classified as with high or low titers of these antibodies based on the median titers of the antibodies among all patients with NPC. Information on EBV Zta-IgA and EBNA-IgA antibodies (classified as negative or positive) was obtained from medical records review and only available for NPC patients. EBNA1-IgA and Zta-IgA antibodies were measured in Wuzhou Red Cross Hospital using commercial kits (Zhongshan Biotech, Zhongshan, China) based on enzyme-linked immunosorbent assay (ELISA).
Lipid Extraction and LC–MS/MS
We collected another sample of remaining morning fasting blood in the anticoagulant-treated tubes at 2,000 rmp at 4°C for 10 min. The supernatant plasma was transferred, aliquoted, and stored at −80°C before the samples were sent to the Shanghai Institute of Clinical Bioinformatics where 20 μl plasma was added into a glass tube to extract the total lipids using the isopropanol precipitation method (Sarafian et al., 2014). Briefly, 9 μl internal standard cocktails (AB SCIEX, 5040146-5,040,155) and 350 μl pre-cold isopropanol were added to each tube which was then shaken vigorously and incubated at −20°C overnight. On the next day, samples were centrifuged at 14,000 rpm at 4°C for 20 min with the supernatant transferred and centrifuged again under the same condition. Finally, the supernatant was ready for lipid extraction.
All samples were measured in the normal-phase liquid chromatography-coupled triple-quadrupole mass spectrometers (Q-trap® 5,500 LC–MS/MS, Sciex, Framingham, MA, United States of America). The determination conditions used in LC–MS/MS were as follows: Acquity UPLC BEH HILIC (Waters Co., Milford, MA, US) column (100 mm × 2.1 mm, 1.7 μm) was used as stationary phase, acetonitrile solution (containing 10 mmol/l ammonium acetate) was used as mobile phase A, and 50% acetonitrile solution (containing 10 mmol/l ammonium acetate) was used as mobile phase B. The gradient elution of mobile phase B: 0–10 min, raised from 0.1 to 20%; 10–11 min, raised from 20 to 98%; 11–13 min, maintained at 98%; 13.1 min, reduced to 0.1%; and 13.2–16 min, maintained at 0.1%. Other conditions mainly included a flow rate of 0.5 ml/min, a column temperature of 35°C, an injection volume of positive/negative ion mode at 4 μl, a gas curtain pressure at 35 psi, a nebulizer pressure at 50 psi, an auxiliary gas pressure at 60 psi, and a heating temperature at 500°C. The positive and negative ion mode ion spray voltages were − 5,500 V and 5,500 V, the positive and negative ion mode optimized de-clustering voltages were 80 V and − 80 V, whereas the positive and negative injection voltages were 10 V and − 10 V, respectively.
An equal amount of liquid harvested from each sample was mixed and used as a quality control (QC) mixture. The QC sample was treated the same as the study samples. One QC sample was added to each ten study samples. The stability of LC–MS/MS was evaluated by principal component analysis of QC samples and other experimental samples. On the principal component diagram, the projection points of QC samples were relatively clustered, demonstrating that LC–MS/MS had good stability and repeatability. The results also showed that the differences among the test samples were mainly due to the variations in the metabolites of the test samples instead of the variations of the test methods. The production was scanned by Q-Trap in multiple reaction monitoring modes.
Statistical Analysis
Among all detected lipids, we retained the ones with > 30% of non-zero values among either the cases or controls and disregarded the rest. After this data filtering, we retained 1,125 of the total 1,191 lipid molecules detected in LC–MS/MS. Using the filtered data, we performed a case–control comparison as well as subgroup analyses of the NPC patients by EBV antibodies, tumor stage, and sex. To address whether the sex-specific lipid profiles were specific to NPC, we also compared the lipids between male and female controls. In the subgroup analyses, if a lipid molecule had < 5% non-zero values in either of the comparison groups, we discarded such comparison. We used unpaired Welch t test in all comparisons, considering the homogeneity of variance of the samples. Both fold-change value and p value were calculated in these analyses. A fold-change value was computed for each lipid molecule through comparing the averages of different groups. A p value <0.05 was considered statistically significant. Volcano plots were used to present the lipid molecules that differed statistically significantly between groups (with a fold-change > 2 and a p value < 0.05).
We also performed multivariable analyses, namely partial least squares discriminant analysis (PLS-DA) and orthogonal partial least squares discriminant analysis (OPLS-DA) to examine whether the lipid profiles can be used to discriminate NPC patients from healthy controls. PLS-DA and OPLS-DA are supervised methods to establish a model of the relationship between the predictor (i.e., a lipid profile) and outcome (i.e., case–control status) and to subsequently estimate the usefulness of the predictor (Worley and Powers, 2013). Parameters R2 and Q2 were used to evaluate the interpretability and predictability of the PLS-DA and OPLS-DA models, respectively. A VIP value was used to indicate the contribution of the predictor to the OPLS-DA model to discriminate cases from healthy controls. The greater the contribution, the higher the VIP value. We used a 7-fold cross-validation and permutation test with 200 repeats when constructing the PLS-DA model. R version 4.0.5 was used in all data analyses and for graphic visualization.
Results
Table 1 shows the characteristics of the NPC patients and their age- and sex-matched controls. There were 27 female and 73 male among the NPC patients, with a mean age of 49.4 years (range: 23–87 years), and 27 female and 73 male among the controls, with a mean age of 48.6 years (range: 17–91 years).
Table 1. Characteristics of nasopharyngeal carcinoma (NPC) patients and their age- and sex-matched healthy controls.
We identified 1,191 lipids in both the NPC patients and controls, including 3 cholesterols, 21 cholesterol esters, 2 ceramides, 12 sphingomyelins, 1 sphingosine, 1 sphingosine-1-phosphate, 17 monoacylglycerols, 50 diacylglycerols, 445 triglycerides, 77 phosphatic acids, 141 phosphatidylethanolamines, 77 phosphatidylglycerols, 76 phosphatidylinositols, 77 phosphatidylserines, 78 phosphatidylcholine, 13 lyso-phosphatic acids, 16 lyso-phosphatidylethanolamines, 16 lyso-phosphatidylglycerols, 16 lyso-phosphatidylinositols, 16 lyso-phosphatidylserines, 16 lyso-phosphatidylcholine and 20 free fatty acids. After data filtering, we retained 1,125 of these lipids in the analysis.
Lipid Profiles Between NPC Patients and Controls
Statistically significant differences were noted for 349 lipids between NPC patients and controls (p value < 0.05). Figure 1 shows the lipid molecules with a fold-change > 2 (alternatively < 0.5) of these 349 lipids, when comparing NPC patients to controls. Figure 2A shows the results of PLS-DA model and Figure 2B shows the results of OPLS-DA model. The R2 and Q2 values of the score plots suggest that the models are stable with good predictive ability. A statistically significant difference was noted for 21 lipid molecules between NPC patients and controls, including three classes, namely free fatty acids, glycerolipids, and glycerophospholipids, as well as seven sub-classes (Table 2). Among these, diacylglycerols 16:0/20:3 and 16:0/16:1 showed a higher level (Figure 3), whereas the other 19 lipids, including diacylglycerol 20:0/20:0, free fatty acids 16:1, 18:0, 18:1, 18:3, 20:0, 20:1, 20:4, and 22:6, lyso-phosphatic acid 18:0, phosphatic acid 20:0/20:5, phosphatidylglycerol 18:0/20:5, phosphatidylinositols 16:0/18:2, 18:0/18:2, 18:0/20:3, 18:0/20:4, and 18:1/18:1, as well as phosphatidylserines 14:0/18:3 and 16:0/14:0 showed a lower level (Figure 4), among NPC patients than controls.
Figure 1. Volcano plot for the 21 lipid molecules showing statistically significantly different levels between nasopharyngeal carcinoma (NPC) patients and controls. The x-axis shows Log2 (FC), the y-axis shows –Log10 (p value) and each dot represents one lipid. The big red and green dots above the straight line [−Log10(0.05)] show lipids with statistically significantly different levels between NPC patients and controls as well as a fold-change (FC) greater than 2 (alternatively smaller than 0.5). PI, phosphatidylinositol; FFA, free fatty acid; DAG, diacylglycerol.
Figure 2. Distinguishing NPC patients from controls using supervised statistical analyses. (A) NPC patients and controls are partially distinguished based on 2-dimensional scores made by PLS-DA, where Q2 > 0.5 and R2 > 0.5 indicate the predictive ability of the model. A blue square represents a healthy control whereas a red square represents an NPC patient. (B) Score plot of OPLS-DA showing the separation of NPC patients and controls, where Q2 > 0.5 and R2 > 0.5 indicate the predictive ability of the model. A blue square represents a healthy control whereas a red square represents an NPC patient.
Figure 3. Violin plots showing the two lipids with a higher level among NPC patients compared with controls. * means p < 0.05. NPC means nasopharyngeal carcinoma. Control means healthy controls. DAG, diacylglycerol.
Figure 4. Violin plots showing the 19 lipids with a lower level among NPC patients compared with controls. * means p < 0.05, ** means p < 0.01, *** means p < 0.001, **** means p < 0.0001, and “ns” means p > 0.05. DAG, diacylglycerol; FFA, free fatty acid; PI, phosphatidylinositol; PG, phosphatidylglycerol; PS, phosphatidylserine; PA, phosphatidic acid; LPA, lyso- phosphatidic acid.
Lipid Profiles in NPC Patients by EBV Antibodies
There was a lower level of one lipid (phosphatidic acid 14:0/20:1) whereas a higher level of 19 lipids, including 18 glycerolipids, among NPC patients with a lower titer of EA-IgA antibodies, compared with NPC patients with a higher titer of EA-IgA antibodies (Table 3). There was a higher level of 52 lipids, all glycerolipids, among NPC patients with a lower titer of VCA-IgA antibodies, compared with patients with a higher titer of VCA-IgA (Table 4). NPC patients with the presence of Zta-IgA had a higher level of 2 lipids (diacylglycerol16:1/20:2 and triacylglycerol54:0-FA16:0), whereas a lower level of 1 lipid (phosphatidic acid 16:0/14:0), compared with NPC patients without Zta-IgA. Finally, NPC patients with the presence of EBNA1-IgA had a lower level of 2 lipids whereas a higher level of 143 lipids, compared with patients without EBNA1-IgA. Table 5 shows the top 20 lipids of these.
Lipid Profiles in NPC Patients by Tumor Stage
No different levels of lipids were found when comparing NPC patients with a T1 or T2 tumor to those with a T3 or T4 tumor (all p values > 0.05). However, NPC patients with a N0 tumor had a higher level of 1 lipid whereas a lower level of another 112 lipids, compared with patients with a N1, N2, or N3 tumor (Table 6). NPC patients with a M0 tumor had a higher level of four lipids, whereas a lower level of another four lipids, compared with patients of a M1 tumor (Table 7). Finally, there was a lower level of 93 lipids, whereas a higher level of 1 lipid, among NPC patients with an early stage (stage II) cancer, compared with patients with an advanced stage (stage III, IVA or IVB) cancer (Table 8).
Table 6. Top 20 lipids with different levels between NPC patients with a N0 tumor and NPC patients with N1, N2, or N3 tumor.
Table 7. Lipids with different levels between NPC patients with a M0 stage tumor and NPC patients with a M1 stage tumor.
Table 8. Top 30 lipids with different levels between NPC patients with a stage II tumor and NPC patients with a stage III or IV tumor.
Lipid Profiles in NPC Patients and Controls by Sex
Male NPC patients had a lower level of 8 lipids, whereas a higher level of 1 lipid, when compared to female NPC patients (Table 9). Statistically significant differences were also noted for 49 lipids between male and female controls (Table 10). Male controls had a lower level of 42 glycerophospholipids, whereas a higher level of seven free fatty acids, compared with female controls.
Discussion
In a study of 100 patients with NPC and 100 age- and sex-matched healthy controls, we used non-targeted lipid metabolomics techniques based on LC–MS/MS to compare the plasma lipid profiles of NPC patients to that of healthy controls and to analyze lipid profiles of the NPC patients by EBV antibodies, tumor stage, and sex. We found that 21 lipids, including 10 glycerophospholipids, 3 glycerolipids, and 8 free fatty acids, demonstrated differential levels between patients and controls. Among these lipids, majority showed a lower level, whereas 2 glycerolipids showed a higher level, among NPC patients than controls. We also found the lipid profiles to differ by the titer and presence of different EBV antibodies, as well as by tumor stage and sex.
Although plasma lipidomics has rarely been studied in NPC, reprogramming of energy metabolism has been recognized as a hallmark (Hanahan and Weinberg, 2011) and alteration of lipid metabolism is known as a prominent metabolic feature (Beloribi-Djefaflia et al., 2016; Cheng et al., 2018; Snaebjornsson et al., 2020) in cancer. Further, increased uptake of exogenous lipids is known to promote cancer cell proliferation, growth, invasion, and migration (Snaebjornsson et al., 2020), whereas serum lipids are energy sources for tumor progression (Ward and Thompson, 2012; Beloribi-Djefaflia et al., 2016). Similarly, although the exact roles of different lipids on NPC are unknown, existing literature has suggested potential mechanisms linking together lipids to carcinogenesis in general. For instance, the mTORC2 signaling pathway promotes carcinogenesis whereas mTORC2 stimulates de novo lipid synthesis, especially synthesis of sphingolipids and glycerophospholipids (Guri et al., 2017). Glycerophospholipids, including phosphatic acids, phosphatidylglycerols, phosphatidylinositols, phosphatidylserines and lyso-phosphatic acids, are the most abundant lipids observed in our study. Phosphatidylinositols are the precursor of the second messenger inositol 1,4,5-trisphosphate and diacylglycerol; the former controls cytosolic calcium levels and the latter activates the protein kinase C kinase family (Goncalves et al., 2018). An increase in phosphatidylinositols has been shown in renal clear cell carcinoma (Neumann et al., 2020) and colon cancer (Hofmanova et al., 2020). On the other hand, inhibition of lyso-phosphatic acids was shown to reduce tumor formation in mouse ovaries (Carrasco and Merida, 2007). Phospholipids are the most prominent membrane lipids in mammalian cells (van Meer and Sprong, 2004) and have been shown to be increased in tumor cells (Hofmanova et al., 2020), suggesting that elevated phospholipids may be essential for tumor cells to assemble membrane structure to sustain proliferation. In addition to phospholipids, intracellular lipids, especially fatty acids, and cholesterols, are also involved in tumor cell proliferation, apoptosis, stemness, invasiveness, and epithelial-mesenchymal transition (Koundouros and Poulogiannis, 2020). In a previous study, we found that fatty acids were upregulated whereas the ISG15-conjugating enzyme UbcH8 was epigenetically inactivated in NPC cells (Zhou et al., 2015). In another study, we further showed that the re-expression of SLC27A6 might increase fatty acid uptake in NPC cells (Zhong et al., 2021).
EBV infection plays an important role in the development of NPC, especially in the endemic areas such as Southeast Asia. The titers of serum EBV antibodies are known to increase dramatically in NPC patients (Chen et al., 2015). Interestingly, EBV has also been shown to be able to infect primary human adipocytes in vitro, leading to increased release of free fatty acids, glycerol, and expression of proinflammatory adipokines (Liu et al., 2021). Further, EBV-encoded latent membrane protein 1 and EBV-encoded RNAs were shown to promote cell proliferation by upregulating fatty acid synthase in NPC cells (Lo et al., 2018) whereas EBV-encoded latent membrane protein 2A (LMP2A) were shown to promote lipid accumulation by blocking adipose triglyceride lipase in NPC cells, resulting in an enhanced cell migration (Zheng et al., 2020). Taken together, our results of different lipid profiles between NPC patients with different status of EBV antibodies suggest that EBV infection might modulate the risk and prognosis of NPC through the reprogramming of lipid metabolism.
In the host, EBV can establish two alternative modes of life cycle, known as latent infection and lytic infection. The switch from latent infection to lytic infection is known as EBV reactivation (Li et al., 2016). EBV found in NPC cells is mostly in latent cycle, with the viral genome existing as an episome in the nucleus and only a few viral genes, including EBNA1, expressed (Kempkes and Robertson, 2015). However, it is speculated that it is a small number of cells with lytic infection of EBV that promote NPC carcinogenesis and progression, because antibodies against lytic infection, such as EA-IgA and VCA-IgA, were found to be elevated several years before clinical onset of NPC, whereas levels of these antibodies were correlated with disease progression (Chien et al., 2001; Cao et al., 2011). It is therefore highly plausible that reactivation of EBV is critical in the initiation and progression of NPC. EBV antibodies are effective markers for different phases of EBV infection. In the present study, we found that NPC patients with positive EBNA1-IgA antibodies had a higher level of 143 lipids, but a lower level of only two lipids, compared to NPC patients with negative EBNA1-IgA antibodies. However, NPC patients with a higher titer of either EA-IgA or VCA-IgA antibodies had a lower level of lipids, compared to NPC patients with a lower titer of these antibodies. The different characteristics of lipid metabolism during the latent (indicated by positive EBNA1-IgA) and lytic (indicated by high titers of EA-IgA and VCA-IgA) phases of EBV infection suggest that altered lipid metabolism in relation to or as a result of EBV reactivation might be an important mechanism linking together EBV infection, especially EBV reactivation, to the development of NPC. When studying lipid profiles in relation to tumor stage of the NPC patients, we found that there was a lower level of many lipids but a higher level of only few lipids in early stage (e.g., N0 and stage II), compared with advanced stage (e.g., N1 and stage III, IVA or IVB) NPC, suggesting a dynamic change of lipid metabolism during tumor progression. It is interesting to note that these lipids are not the same as the lipids showing a different level between NPC patients and healthy controls. These results might therefore suggest that different characteristics of lipid metabolism are present in the initiation and progression of NPC. Further, NPC patients with advanced stage tumor demonstrated only a lower level of 1 lipid, namely phosphatidylinositol 18:0/20:1, compared with NPC patients with early-stage tumor. In the comparison between NPC patients and healthy controls, we also found that NPC patients demonstrated a lower level of many phosphatidylinositols. It is therefore possible that these phosphatidylinositols are originated from NPC cells and are regulated during NPC progression. If this hypothesis is true, phosphatidylinositols might be of potential use as biomarkers for the diagnosis and prognostic indication of NPC. However, as the pattern was less clear in terms of T stage and metastasis status, these findings need to be validated in future studies.
The sex-specific results among NPC patients and healthy controls suggest distinct lipid profiles between male and female, which are independent of NPC status. Although further research is clearly needed in unraveling the underlying mechanisms for such pattern, we speculate that sex-specific factors such as sex hormones might contribute.
The main strengths of the study are the use of lipidomics analysis based on LC–MS/MS and the rich data on clinical characteristics for NPC patients such as EBV antibodies and TNM and clinical stages. Given the cross-sectional design of the present study, a temporal relationship cannot be examined between EBV infection, alterations in lipid profile, and onset of NPC, precluding the possibility to discuss causation. The observed lipid profile alterations among patients with NPC, compared with healthy controls, might therefore be attributable to different reasons. First, the alterations might be due to epiphenomena, i.e., the lipids change as a result of NPC. Second, the alterations are secondary to EBV infection, especially EBV reactivation, which leads subsequently to the development of NPC. Third, the alterations are causally related to NPC development, i.e., changed lipid profiles, independently of or jointly with EBV infection, lead to the development of NPC. Further, as our primary aim was to understand the potential role of altered lipid metabolism in the causation of NPC, we only analyzed lipid profiles before any treatment was given, precluding the possibility to examine the role of cancer treatment on lipids. Prospective studies with longitudinal plasma samples collected from long before onset of NPC (preferably at the time of EBV reactivation) until after cancer diagnosis (preferably after cancer treatment) will help to better understand the sequence of these events as well as the impact of cancer treatment on lipid profiles.
Another limitation of the study is the lack of separate assay to validate the observed lipid profiles using LC–MS/MS. This is however not unique to the present study as such validation is currently non-existent in lipidomics studies. A further limitation of the study is the lack of independent samples to validate the observed alterations in lipid profiles among patients with NPC, compared to healthy controls. As a result, we deemed the present study as a pilot effort and call for future studies to confirm (or refute) our findings. We did not focus on the prediction ability of individual lipids in the study, but rather the prediction ability of a composite lipid profile. The prediction ability of the model is likely influenced by sample size as well as noise due to confounding factors and the accuracy of lipid measurements. Future studies with larger sample size and better control of confounding are therefore needed to further assess the prediction ability of our model. Finally, as our study was based on an endemic area of NPC, whether our findings are generalizable to NPC in non-endemic areas remains to be examined.
Conclusion
In conclusion, plasma lipidomics might help to differentiate NPC cases from controls whereas EBV infection might contribute to the risk and prognosis of NPC through modulating lipid metabolism in tumor cells as well as in the periphery.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics Statement
The studies involving human participants were reviewed and approved by ethics approval of this study was obtained from the Ethical Evaluation Committee of Wuzhou Red Cross Hospital (No. LL2017-19, Wuzhou, China). Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin.
Author Contributions
ZZ and YC conceived the idea and designed the study. YH and JL performed the experiments. WH, XP, YL, XZ, and YY contributed to the data analysis. XX and WZ drafted the manuscript. XZ contributed to the result interpretation, critically reviewed the manuscript, and supervised experiments. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the National Natural Science Foundation of China (81960490, 81760489, 82060511, and 81860601), the Guangxi Natural Science Foundation of China (2020GXNSFAA297105), the Youth Program of Guangxi Natural Science Foundation of China (2018GXNSFBA281158, 2018GXNSFBA281028), the High-level Talent Introduction Plan of the First Affiliated Hospital of Guangxi Medical University (the fifth level), and the fund of Guangxi Key Laboratory of Early Prevention in Regional High Incidence Cancer (GKE-ZZ 202119).
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.
Abbreviations
NPC, Nasopharyngeal Carcinoma; EBV, Epstein–Barr virus; EA, Early Antigen; VCA, Viral Capsid Antigen; EBNA1, EBV Nuclear Antigen 1; Zta-IgA, BZLF1 Transcription Activator Protein; PLS-DA, Partial Least Squares Discriminant Analysis; OPLS-DA, Orthogonal Partial Least Squares Discriminant Analysis; LC–MS/MS, Liquid Chromatography-Tandem Mass Spectrometry.
References
Beloribi-Djefaflia, S., Vasseur, S., and Guillaumond, F. (2016). Lipid metabolic reprogramming in cancer cells. Oncogenesis 5:e189. doi: 10.1038/oncsis.2015.49
Binné, U. K., Amon, W., and Farrell, P. J. (2002). Promoter sequences required for reactivation of Epstein-Barr virus from latency. J. Virol. 76, 10282–10289. doi: 10.1128/JVI.76.20.10282-10289.2002
Cai, Y. L., Li, J., Lu, A. Y., Zheng, Y. M., Zhong, W. M., Wang, W., et al. (2014). Diagnostic significance of combined detection of Epstein-Barr virus antibodies, VCA/IgA, EA/IgA, Rta/IgG and EBNA1/IgA for nasopharyngeal carcinoma. Asian Pac. J. Cancer Prev. 15, 2001–2006. doi: 10.7314/apjcp.2014.15.5.2001
Cao, S. M., Liu, Z., Jia, W. H., Huang, Q. H., Liu, Q., Guo, X., et al. (2011). Fluctuations of epstein-barr virus serological antibodies and risk for nasopharyngeal carcinoma: a prospective screening study with a 20-year follow-up. PLoS One 6:e19100. doi: 10.1371/journal.pone.0019100
Carracedo, A., Cantley, L. C., and Pandolfi, P. P. (2013). Cancer metabolism: fatty acid oxidation in the limelight. Nat. Rev. Cancer 13, 227–232. doi: 10.1038/nrc3483
Carrasco, S., and Merida, I. (2007). Diacylglycerol, when simplicity becomes complex. Trends Biochem. Sci. 32, 27–36. doi: 10.1016/j.tibs.2006.11.004
Chen, Y. P., Chan, A. T. C., Le, Q. T., Blanchard, P., Sun, Y., and Ma, J. (2019). Nasopharyngeal carcinoma. Lancet 394, 64–80. doi: 10.1016/S0140-6736(19)30956-0
Chen, Y., Zhao, W., Lin, L., Xiao, X., Zhou, X., Ming, H., et al. (2015). Nasopharyngeal Epstein-Barr virus load: an efficient supplementary method for population-based nasopharyngeal carcinoma screening. PLoS One 10:e0132669. doi: 10.1371/journal.pone.0132669
Cheng, C., Geng, F., Cheng, X., and Guo, D. (2018). Lipid metabolism reprogramming and its potential targets in cancer. Cancer Commun. (Lond) 38:27. doi: 10.1186/s40880-018-0301-4
Chien, Y. C., Chen, J. Y., Liu, M. Y., Yang, H. I., Hsu, M. M., Chen, C. J., et al. (2001). Serologic markers of Epstein-Barr virus infection and nasopharyngeal carcinoma in Taiwanese men. N. Engl. J. Med. 345, 1877–1882. doi: 10.1056/NEJMoa011610
Coghill, A. E., Hsu, W. L., Pfeiffer, R. M., Juwana, H., Yu, K. J., Lou, P. J., et al. (2014). Epstein-Barr virus serology as a potential screening marker for nasopharyngeal carcinoma among high-risk individuals from multiplex families in Taiwan. Cancer Epidemiol. Biomark. Prev. 23, 1213–1219. doi: 10.1158/1055-9965.EPI-13-1262
Daker, M., Bhuvanendran, S., Ahmad, M., Takada, K., and Khoo, A. S. (2013). Deregulation of lipid metabolism pathway genes in nasopharyngeal carcinoma cells. Mol. Med. Rep. 7, 731–741. doi: 10.3892/mmr.2012.1253
Fichtali, K., Bititi, A., Elghanmi, A., and Ghazi, B. (2020). Serum Lipidomic profiling in breast Cancer to identify screening, diagnostic, and prognostic biomarkers. Biores. Open Access. 9, 1–6. doi: 10.1089/biores.2018.0022
Gao, R., Wang, L., Liu, Q., Zhang, L. F., Ye, Y. F., Xie, S. H., et al. (2017). Evaluation of seven recombinant VCA-IgA ELISA kits for the diagnosis of nasopharyngeal carcinoma in China: a case-control trial. BMJ Open 7:e013211. doi: 10.1136/bmjopen-2016-013211
Goncalves, M. D., Hopkins, B. D., and Cantley, L. C. (2018). Phosphatidylinositol 3-kinase, growth disorders, and Cancer. N. Engl. J. Med. 379, 2052–2062. doi: 10.1056/NEJMra1704560
Guri, Y., Colombi, M., Dazert, E., Hindupur, S. K., Roszik, J., Moes, S., et al. (2017). mTORC2 promotes tumorigenesis via lipid synthesis. Cancer Cell 32:e12. doi: 10.1016/j.ccell.2017.11.011
Han, X. (2016). Lipidomics for studying metabolism. Nat. Rev. Endocrinol. 12, 668–679. doi: 10.1038/nrendo.2016.98
Hanahan, D., and Weinberg, R. A. (2011). Hallmarks of cancer: the next generation. Cell 144, 646–674. doi: 10.1016/j.cell.2011.02.013
Hoch, J., Dryjski, M., Jarrell, B. E., Carabasi, R. A., and Williams, S. K. (1988). In vitro endothelialization of an aldehyde-stabilized native vessel. J. Surg. Res. 44, 545–554. doi: 10.1016/0022-4804(88)90160-6
Hofmanova, J., Slavik, J., Ovesna, P., Tylichova, Z., Dusek, L., Strakova, N., et al. (2020). Phospholipid profiling enables to discriminate tumor- and non-tumor-derived human colon epithelial cells: Phospholipidome similarities and differences in colon cancer cell lines and in patient-derived cell samples. PLoS One 15:e0228010. doi: 10.1371/journal.pone.0228010
Hu, C., van der Heijden, R., Wang, M., van der Greef, J., Hankemeier, T., and Xu, G. (2009). Analytical strategies in lipidomics and applications in disease biomarker discovery. J. Chromatogr. B Analyt. Technol. Biomed. Life Sci. 877, 2836–2846. doi: 10.1016/j.jchromb.2009.01.038
Hui, K. F., and Chiang, A. K. S. (2010). Suberoylanilide hydroxamic acid induces viral lytic cycle in Epstein-Barr virus-positive epithelial malignancies and mediates enhanced cell death. Int. J. Cancer 126, 2479–2489. doi: 10.1002/ijc.24945
Ji, M. F., Wang, D. K., Yu, Y. L., Guo, Y. Q., Liang, J. S., Cheng, W. M., et al. (2007). Sustained elevation of Epstein-Barr virus antibody levels preceding clinical onset of nasopharyngeal carcinoma. Br. J. Cancer 96, 623–630. doi: 10.1038/sj.bjc.6603609
Jiang, N., Zhang, Z., Chen, X., Zhang, G., Wang, Y., Pan, L., et al. (2021). Plasma Lipidomics profiling reveals biomarkers for papillary thyroid cancer diagnosis. Front. Cell Dev. Biol. 9:682269. doi: 10.3389/fcell.2021.682269
Kempkes, B., and Robertson, E. S. (2015). Epstein-Barr virus latency: current and future perspectives. Curr. Opin. Virol. 14, 138–144. doi: 10.1016/j.coviro.2015.09.007
Koundouros, N., and Poulogiannis, G. (2020). Reprogramming of fatty acid metabolism in cancer. Br. J. Cancer 122, 4–22. doi: 10.1038/s41416-019-0650-z
Leung, K. S., and Fong, B. M. (2014). LC-MS/MS in the routine clinical laboratory: has its time come? Anal. Bioanal. Chem. 406, 2289–2301. doi: 10.1007/s00216-013-7542-5
Li, H., Liu, S., Hu, J., Luo, X., Li, N., Bode, A. M., et al. (2016). Epstein–Barr virus lytic reactivation regulation and its pathogenic role in carcinogenesis. Int. J. Biol. Sci. 12, 1309–1318. doi: 10.7150/ijbs.16564
Liu, Y., Huang, Q., Liu, W., Liu, Q., Jia, W., Chang, E., et al. (2012). Establishment of VCA and EBNA1 IgA-based combination by enzyme-linked immunosorbent assay as preferred screening method for nasopharyngeal carcinoma: a two-stage design with a preliminary performance study and a mass screening in southern China. Int. J. Cancer 131, 406–416. doi: 10.1002/ijc.26380
Liu, Y. Y., Lin, S. J., Chen, Y. Y., Liu, L. N., Bao, L. B., Tang, L. Q., et al. (2016). High-density lipoprotein cholesterol as a predictor of poor survival in patients with nasopharyngeal carcinoma. Oncotarget 7, 42978–42987. doi: 10.18632/oncotarget.7160
Liu, R., Peng, Y., Li, X., Wang, Y., Pan, E., Guo, W., et al. (2013). Identification of plasma metabolomic profiling for diagnosis of esophageal squamous-cell carcinoma using an UPLC/TOF/MS platform. Int. J. Mol. Sci. 14, 8899–8911. doi: 10.3390/ijms14058899
Liu, S. C., Tsang, N. M., Lee, P. J., Sui, Y. H., Huang, C. H., and Liu, T. T. (2021). Epstein-Barr virus induces adipocyte dedifferentiation to modulate the tumor microenvironment. Cancer Res. 81, 3283–3294. doi: 10.1158/0008-5472.CAN-20-3121
Lo, A. K., Lung, R. W., Dawson, C. W., Young, L. S., Ko, C. W., Yeung, W. W., et al. (2018). Activation of sterol regulatory element-binding protein 1 (SREBP1)-mediated lipogenesis by the Epstein-Barr virus-encoded latent membrane protein 1 (LMP1) promotes cell proliferation and progression of nasopharyngeal carcinoma. J. Pathol. 246, 180–190. doi: 10.1002/path.5130
Luo, X., Cheng, C., Tan, Z., Li, N., Tang, M., Yang, L., et al. (2017). Emerging roles of lipid metabolism in cancer metastasis. Mol. Cancer 16:76. doi: 10.1186/s12943-017-0646-3
Moss, D. J., Burrows, S. R., Silins, S. L., Misko, I., and Khanna, R. (2001). The immunology of Epstein-Barr virus infection. Philosophical transactions of the Royal Society of London. Ser. B Biol. Sci. 356, 475–488. doi: 10.1098/rstb.2000.0784
Neumann, C. K. A., Silver, D. J., Venkateshwari, V., Zhang, R., Traughber, C. A., Przybycin, C., et al. (2020). MBOAT7-driven phosphatidylinositol remodeling promotes the progression of clear cell renal carcinoma. Mol. Metab. 34, 136–145. doi: 10.1016/j.molmet.2020.01.011
Pathmanathan, R., Prasad, U., Chandrika, G., Sadler, R., Flynn, K., and Raab-Traub, N. (1995). Undifferentiated, nonkeratinizing, and squamous cell carcinoma of the nasopharynx. Variants of Epstein-Barr virus-infected neoplasia. Am. J. Pathol. 146, 1355–1367.
Raab-Traub, N. (2015). Nasopharyngeal carcinoma: an evolving role for the Epstein-Barr virus. Curr. Top. Microbiol. Immunol. 390, 339–363. doi: 10.1007/978-3-319-22822-8_14
Robertson, E., and Kieff, E. (1995). Reducing the complexity of the transforming Epstein-Barr virus genome to 64 kilobase pairs. J. Virol. 69, 983–993. doi: 10.1128/jvi.69.2.983-993.1995
Sarafian, M. H., Gaudin, M., Lewis, M. R., Martin, F. P., Holmes, E., Nicholson, J. K., et al. (2014). Objective set of criteria for optimization of sample preparation procedures for ultra-high throughput untargeted blood plasma lipid profiling by ultra performance liquid chromatography-mass spectrometry. Anal. Chem. 86, 5766–5774. doi: 10.1021/ac500317c
Shahid, M., Yeon, A., and Kim, J. (2020). Metabolomic and lipidomic approaches to identify biomarkers for bladder cancer and interstitial cystitis (review). Mol. Med. Rep. 22, 5003–5011. doi: 10.3892/mmr.2020.11627
Snaebjornsson, M. T., Janaki-Raman, S., and Schulze, A. (2020). Greasing the wheels of the Cancer machine: the role of lipid metabolism in Cancer. Cell Metab. 31, 62–76. doi: 10.1016/j.cmet.2019.11.010
Stelow, E. B., and Wenig, B. M. (2017). Update From The 4th edition of the World Health Organization classification of head and neck tumours: Nasopharynx. Head Neck Pathol. 11, 16–22. doi: 10.1007/s12105-017-0787-0
Sung, H., Ferlay, J., Siegel, R. L., Laversanne, M., Soerjomataram, I., Jemal, A., et al. (2021). Global Cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin. 71, 209–249. doi: 10.3322/caac.21660
Tan, B., Qiu, Y., Zou, X., Chen, T., Xie, G., Cheng, Y., et al. (2013). Metabonomics identifies serum metabolite markers of colorectal cancer. J. Proteome Res. 12, 3000–3009. doi: 10.1021/pr400337b
van Meer, G., and Sprong, H. (2004). Membrane lipids and vesicular traffic. Curr. Opin. Cell Biol. 16, 373–378. doi: 10.1016/j.ceb.2004.06.004
Ward, P. S., and Thompson, C. B. (2012). Metabolic reprogramming: a cancer hallmark even Warburg did not anticipate. Cancer Cell 21, 297–308. doi: 10.1016/j.ccr.2012.02.014
Worley, B., and Powers, R. (2013). Multivariate analysis in metabolomics. Curr. Metabolomics 1, 92–107. doi: 10.2174/2213235X11301010092
Yu, W. M., and Hussain, S. S. (2009). Incidence of nasopharyngeal carcinoma in Chinese immigrants, compared with Chinese in China and South East Asia: review. J. Laryngol. Otol. 123, 1067–1074. doi: 10.1017/S0022215109005623
Yu, M. C., and Yuan, J. M. (2002). Epidemiology of nasopharyngeal carcinoma. Semin. Cancer Biol. 12, 421–429. doi: 10.1016/s1044579x02000858
Zeng, Y., Zhang, L. G., Li, H. Y., Jan, M. G., Zhang, Q., Wu, Y. C., et al. (1982). Serological mass survey for early detection of nasopharyngeal carcinoma in Wuzhou City, China. Int. J. Cancer 29, 139–141. doi: 10.1002/ijc.2910290204
Zhang, L. F., Li, Y. H., Xie, S. H., Ling, W., Chen, S. H., Liu, Q., et al. (2015). Incidence trend of nasopharyngeal carcinoma from 1987 to 2011 in Sihui County, Guangdong Province, South China: an age-period-cohort analysis. Chin. J. Cancer 34, 350–357. doi: 10.1186/s40880-015-0018-6
Zhang, L., Zhu, B., Zeng, Y., Shen, H., Zhang, J., and Wang, X. (2020). Clinical lipidomics in understanding of lung cancer: opportunity and challenge. Cancer Lett. 470, 75–83. doi: 10.1016/j.canlet.2019.08.014
Zheng, S., Matskova, L., Zhou, X., Xiao, X., Huang, G., Zhang, Z., et al. (2020). Downregulation of adipose triglyceride lipase by EB viral-encoded LMP2A links lipid accumulation to increased migration in nasopharyngeal carcinoma. Mol. Oncol. 14, 3234–3252. doi: 10.1002/1878-0261.12824
Zhong, X., Yang, Y., Li, B., Liang, P., Huang, Y., Zheng, Q., et al. (2021). Downregulation of SLC27A6 by DNA Hypermethylation promotes proliferation but suppresses metastasis of nasopharyngeal carcinoma Through modulating lipid metabolism. Front. Oncol. 11:780410. doi: 10.3389/fonc.2021.780410
Keywords: nasopharyngeal carcinoma, Epstein–Barr virus, lipidomics, LC–MS/MS, cancer stage
Citation: Huang Y, Liang J, Hu W, Liang Y, Xiao X, Zhao W, Zhong X, Yang Y, Pan X, Zhou X, Zhang Z and Cai Y (2022) Integration Profiling Between Plasma Lipidomics, Epstein–Barr Virus and Clinical Phenomes in Nasopharyngeal Carcinoma Patients. Front. Microbiol. 13:919496. doi: 10.3389/fmicb.2022.919496
Edited by:
Jianhong Lu, Central South University, ChinaCopyright © 2022 Huang, Liang, Hu, Liang, Xiao, Zhao, Zhong, Yang, Pan, Zhou, Zhang and Cai. 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: Zhe Zhang, emhhbmd6aGUxOTc1QGdtYWlsLmNvbQ==; Yonglin Cai, Y3lsemVuQDE2My5jb20=
†These authors have contributed equally to this work and share first authorship