Next Article in Journal
Vitamin D Alleviates Anxiety and Depression in Elderly People with Prediabetes: A Randomized Controlled Study
Next Article in Special Issue
Identification of Exhaled Metabolites in Children with Cystic Fibrosis
Previous Article in Journal
Impaired Neurovascular Function Underlies Poor Neurocognitive Outcomes and Is Associated with Nitric Oxide Bioavailability in Congenital Heart Disease
Previous Article in Special Issue
The Human Skin Volatolome: A Systematic Review of Untargeted Mass Spectrometry Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Non-Invasive Differential Diagnosis of Cervical Neoplastic Lesions by the Lipid Profile Analysis of Cervical Scrapings

1
National Medical Research Center for Obstetrics Gynecology and Perinatology Named after Academician V.I., Kulakov of the Ministry of Healthcare of Russian Federation, 117997 Moscow, Russia
2
Moscow Institute of Physics and Technology, 141700 Moscow, Russia
3
V.L. Talrose Institute for Energy Problems of Chemical Physics, Russia Academy of Sciences, 119991 Moscow, Russia
4
Translational Medicine Department, Siberian State Medical University, 634050 Tomsk, Russia
*
Author to whom correspondence should be addressed.
Metabolites 2022, 12(9), 883; https://doi.org/10.3390/metabo12090883
Submission received: 22 August 2022 / Revised: 12 September 2022 / Accepted: 13 September 2022 / Published: 19 September 2022
(This article belongs to the Special Issue Personalized Metabolomics)

Abstract

:
Cervical cancer is one of the most common cancers in women with pronounced stages of precancerous lesions. Accurate differential diagnosis of such lesions is one of the primary challenges of medical specialists, which is vital to improving patient survival. The aim of this study was to develop and test an algorithm for the differential diagnosis of cervical lesions based on lipid levels in scrapings from the cervical epithelium and cervicovaginal canal. The lipid composition of the samples was analyzed by high-performance chromato-mass spectrometry. Lipid markers were selected using the Mann–Whitney test with a cutoff value of 0.05 and by projections to latent structures discriminant analysis, where a projection threshold of one was chosen. The final selection of variables for binomial logistic regressions was carried out using the Akaike information criterion. As a result, a final neoplasia classification method, based on 20 logistic regression sub-models, has an accuracy of 79% for discrimination NILM/cervicitis/LSIL/HSIL/cancer. The model has a sensitivity of 83% and a specificity of 88% for discrimination of several lesions (HSIL and cancer). This allows us to discuss the prospective viability of further validation of the developed non-invasive method of differential diagnosis.

1. Introduction

Cervical cancer (CC) ranks fourth among women’s cancers and second among women aged 15–44. According to a report by the International Agency for Research on Cancer (IARC), more than 604,000 new cases of cervical cancer were diagnosed worldwide in 2020, more than half of which were fatal [1].
The following tests are used to detect cervical neoplasia in clinical practice: cytological examination (liquid-based cytology (LBC), Pap test), HPV test, colposcopy, and biopsy, followed by histological examination. In some cases, analysis of p16 and Ki-67 markers by immunocytochemistry and immunohistochemistry is used to assess individual risk. The main purpose of a cytological examination is to identify precancerous diseases of the cervix for their timely treatment [2].
The correlation of the results of cytological and histological morphological studies is one of the most important criteria for the quality of CC screening. A study by Alanbay et al. showed that cytology results were inconsistent with subsequent histological diagnosis in 31% of cases [3]. R. Gupta et al. also analyzed the agreement between cytology and histology and found major discrepancies (incorrect differentiation absence of intraepithelial lesion and HSIL or cancer) in 6.4% of cases and minor discrepancies (one-step difference between cervical smear and cervical biopsy) in 20.4% of cases [4]. In the study of B.A. Crothers et al. discrepancies were noted in 20.4% of cases [5]. The multicenter KGOG 1040 study revealed 7.9% discrepancies between cytological classification into groups of negative for intraepithelial lesion or malignancy (NILM), atypical squamous cells of undetermined significance (ASCUS), low-grade squamous intraepithelial lesion (LSIL) and histological classification into groups of the high-grade squamous intraepithelial lesion (cervical intraepithelial neoplasia (CIN II–III) and cervical cancer) [6]. These results oblige researchers to search for and develop new, more accurate non-invasive clinical and diagnostic approaches to the diagnosis of neoplastic processes in the cervix. Today, there is much discussion about the involvement of new and rapidly developing omics studies such as metabolomics, proteomics, and transcriptomics [7]. Thus, significant efforts are being focused on developing a panel of genetic markers, which methylation may be associated with cervical neoplastic lesions [8,9]. Duraipandian et al. showed a difference in the molecular composition of tissues affected by precancerous lesions and normal tissues [10]. Paraskevaidi demonstrated LA-REIMS potential to differentiate CIN2+ grade lesions [11]. Furthermore, Porcaci et al. discovered metabolites that can accurately detect HSIL+ tissues [12].
Lipid metabolism disorder is one of the most prominent metabolic changes observed in cancer. Increased synthesis or uptake of lipids contributes to the rapid growth of cancer cells and the formation of tumors [13,14,15,16,17]. Thus, it is extremely important to study and evaluate the diagnostic potential of cervical lipidome changes in tissues with HPV-associated lesions, including CC, using mass spectrometry. Our group has previously created an accurate, rapid diagnostic method based on lipidome (ESI-MS analysis of a cervical epithelial tissue direct injection) [16]. This method can be further developed as an adjunct to histological analysis. However, cervical cancer screening requires a non-invasive approach, in particular, using samples of scrapings for cytological examination. The creation of a panel of lipid markers may allow accurate non-invasive differential diagnosis of the severity of cervical lesions and predict the outcomes of neoplastic transformation. The aim of this investigation was to study the possibility of non-invasive diagnosis of the cervix HPV-associated lesions based on the analysis of the lipid composition of the cervix epithelium.

2. Materials and Methods

The study included 111 patients aged 21 to 45 years (mean age 34 ± 9 years).
Inclusion Criteria:
  • Age from 21 to 45 years;
  • High-risk HPV;
  • Low and high-grade squamous intraepithelial lesions, cervical cancer confirmed by histology examination;
  • Regular menstrual cycle;
  • Ability to comply with protocol requirements;
  • Providing written informed consent to participate in the study.
Exclusion Criteria:
  • Pregnancy;
  • Lactation period;
  • Hormone therapy;
  • Acute inflammation;
  • Dysfunction of the kidneys, liver, lungs in the stage of decompensation;
  • Psychoneurological conditions.
We performed an examination of all patients, including cytological examination, cervical biopsy, histological examination of biopsy material, and lipidomic analysis of scrapings of the cervical epithelium. In cases of CC, cytological evaluation of smears was carried out according to the Bethesda system (2014).
Based on the results of biopsy material histological examination, 5 groups were formed:
  • Group 1 (control)—NILM and HPV-positive (n = 8; 7.2%);
  • Group 2—Chronical cervicitis and HPV-positive (n = 29; 26.1%);
  • Group 3—LSIL (n = 32; 28.8%);
  • Group 4—HSIL (n = 19; 17.1%);
  • Group 5—Cervical cancer (n = 23; 20.7%).
Histological verification of the diagnosis was based on a two-level histopathological classification of CC precancerous processes. According to it, the term “mild epithelial dysplasia” (CIN I) corresponds to LSIL, and the term “moderate and severe cervical intraepithelial neoplasia” (CIN II–III) corresponds to high-grade squamous intraepithelial lesion (HSIL).
Scrapings from the cervical epithelium and cervical canal were taken with a cervical brush. The extraction of lipids from the collected biomaterial was carried out as follows: the brush was placed in an Eppendorf with 500 μL of water/methanol (1/1, v/v) and kept for 5 min in a Vortex, and then for 5 min in an ultrasonic bath. Then it was removed from the Eppendorf. After adding 1 mL of chloroform to the Eppendorf, it was kept in a Vortex for 10 min and then centrifuged for 5 min at 15,000 rpm. The lower organic layer with a volume of 950 µL was taken into a separate vial. After drying in a stream of nitrogen, the sample was redissolved in 200 μL of isopropanol/acetonitrile (1/1, v/v).
The lipid composition of the tissue was analyzed by LC/MS according to a previously developed protocol [18]. Samples were separated on a Dionex UltiMate 3000 chromatograph (Thermo Scientific, Waltham, MA, USA) coupled with Maxis Impact qTOF mass spectrometer (Bruker Daltonics, Bremen, Germany) with electrospray ionization. The analysis was performed in positive and negative ion modes [18]. Data preprocessing (Table S1) and compound identification were carried out according to the protocol developed by J.P. Koelmel et al. by characteristic ion fragments (Table S2, Figure S4) [19]. The lipid nomenclature corresponded to Lipid Maps [20]. The resulting LC/MS data were normalized using autoscaling (Table S3) [21].
Statistical analysis was performed using scripts in the R language (4.1.1) [22] in the RStudio environment (1.383 GNU) [23] with packages “ropls” and “pROC”(Supplementary File S1 contains scripts).
The lipid features for logistic models for classification into any two of the clinical groups were selected in two ways. In the positive ion mode, at the first stage, lipids with p-value < 0.05 according to the Mann–Whitney statistical test were selected. In the negative ion mode, lipids with variable importance in projection (VIP) greater than 1 according to the orthogonal projection on latent structures discriminant analysis (OPLS-DA) (Supplementary File S2 contains statistical summary and S-plot, loading plot, and VIP-plot of OPLS-DA models) were selected. The second stage was the same for both positive and negative ion modes. The lipids selected in the first stage were used in logistic regression models and filtered by the Akaike information criterion [24]. The response variable was the group diagnosis. The milder case was assigned a value of “0” and the more severe case was assigned a value of “1”.
The final classification model was built based on 20 logistic regressions. There were 10 logistic regressions with lipids detected in positive ion mode for all pairwise combinations of the five studied diagnoses and 10 logistic regressions with lipids detected in negative ion mode. The final diagnosis was determined based on which of them received the most “votes”(Figure 1).
Each intermediate model was tested by cross-validation to determine its sensitivity, specificity, and threshold. The accuracy of each diagnosis was assessed as the ratio of the number of true cases of a diagnosis to the total number of cases with such a diagnosis given by the model.

3. Results

3.1. Lipidomic Profile

Histology is the gold standard for diagnosing cervical lesions. All other diagnostic methods are validated by histology. Table 1 demonstrates matches and differences between the cervix state classification results by cytology and histology for the samples studied. There were 29 patients with chronic cervicitis (further cervicitis) according to histological examination, and only 8 of them (27.6%) were revealed during a cytological examination. Similarly, the cytological examination revealed 18 patients with LSIL out of 32 patients with this diagnosis according to histology (56.3%), 10 patients with HSIL out of 19 (52.6%), and 18 patients with CC out of 23 (78.3%). Overall, cytological accuracy was 56%, sensitivity 76%, and specificity 88%.
The diagnostic potential of cervix lipid composition was studied. Lipid extracts from scrapings of the cervical epithelium were analyzed by liquid chromatography with mass spectrometry (HPLC/MS) in positive and negative ion modes. Characteristic chromatograms obtained for the studied groups are shown in Figure S2, and characteristic integrated mass spectra are shown in Figures S3 and S4.
In order to select lipids for diagnostic logistic regression models for discriminant analysis between two clinical groups, a pairwise comparison was performed using Mann–Whitney statistical test. In positive ion mode (Figures S2a and S3), 26 lipids were selected. Most lipids belonged to cholesterol esters, ceramides, lyso- and phosphatidylcholines, phosphatidylethanolamines, and triacylglycerols with an ether bond (Figure 2). Differences in the lipid composition of cervix epithelium when comparing the cervicitis control group with the LSIL and HSIL groups were expressed in a statistically significant decrease in the level of ceramides Cer-NDS d16:0/16:0 and Cer-NDS d16:0/18:0. In the CC group, compared to the control group, increased levels of phosphatidylcholine PC 16:0_18:1, PC 14:0_16:0, PC 16:0_20:4, PC 16:1_18:0, cholesterol esters CE 18:1, CE 24:1 and triacylglycerol TG 10:0_8:0_8:0 were observed. The analysis of cholesterol esters in cervicitis revealed a decrease in their levels in LSIL and HSIL for CE 18:1 and CE 24:1 and an increase in CE 18:1 in CC. The levels of phosphatidylcholines PC 16:0_22:6, PC 16:1_18:2 and sphingomyelins SM d18:0/16:0, SM d18:1/22:0 were also found to decrease in LSIL compared to cervicitis, while phosphatidylcholine PC 14:0_16:0, PC 16:0_16:0 and triacylglycerol TG 16:0_16:0_18:0, TG 14:0_16:0_16:1 levels increase in CC compared to cervicitis (Figure 2, Table S2).
In negative ion mode (Figures S2b and S4), 36 lipids were selected for inclusion in diagnostic logistic regression models. The selected lipids belonged to lyso- and phosphatidylcholines, sphingomyelins, and oxidized phospholipids (Figure 3). There was a statistically significant increase in the levels of phosphatidylcholines PC O-16:0/16:0, PC 16:0_18:2, and sphingomyelin SM d16:0/18:2 in CC compared to those in the NILM group. The cervicitis and LSIL groups had higher levels of lysophosphatidylcholine LPC 18:2 compared to the control group (NILM). There is a significant increase in the levels of phosphatidylcholines and sphingomyelins SM d22:0/20:3, SM d24:1/18:1, SM d20:0/14:0 in CC compared to cervicitis (Figure 3, Table S4).
With malignant transformation of the epithelium of the cervix, the level of particular phosphatidylcholines and sphingomyelins increased (PC 16:0_22:6, PC 14:0_16:0, SM d18:1/20:0 in HSIL vs LSIL groups, and PC 16:0_18:1, LPC 22:3, PC 16:0_20:5, SM d20:0/22:0 in CC vs LSIL). Lipid markers that differentiate precancerous HSIL and cervical cancer (CC) included lipid esters LPC O-16:0, PE P-16:0/20:4, and sphingomyelins SM d16:0/18:1, SM d22:0/20:3, SM d18:0/16:0, levels of which increased in cancer (Figure 2 and Figure 3, Table S4).
Based on the above lipids, a set of discriminative models was developed. Detailed parameters of the models are given in Table 2 and Table S2, including the area under the curve (AUC) with its confidence interval (CI), the threshold of the model, and its sensitivity and specificity. The discriminative models were used to build a classification algorithm, the accuracy of which was 79% (Table 3). The method demonstrated high accuracy in differentiating cervicitis and CC from NILM. The diagnosis of severe lesions (HSIL+) had a sensitivity of 83% and specificity of 88%. The lipidome-based model can be used to refine the diagnosis, and the agreement was higher for histology results than cytological diagnosis (Table 4).

3.2. Clinical Cases

A set of clinical cases were analyzed in more detail to demonstrate the application of the developed models (Figure 1). We re-classified cervical lesions into two groups: NILM, cervicitis, and LSIL—low lesion; HSIL and cancer—high lesion. The probability of belonging to each new class was calculated as division votes for the class (for example, votes for the high lesion is the sum of votes for HSIL and cancer) to sum all votes. The results are presented in Table 5.

3.2.1. Clinical Case No. 1

Patient II (Table 4) had an HSIL cytology result, and HPV type 51 was identified. Colposcopy results revealed mild abnormalities in the epithelium of the cervix (Figure S5a). The epithelial scraping lipidome analysis indicated the absence of neoplasia (Table 5). A biopsy of the cervix was performed. The histological diagnosis resembled chronic cervicitis. This clinical case demonstrates a significant discrepancy between cytology and histology results, with good agreement between histology and lipid model. The latter gives an 85% probability of low-grade lesion (cervicitis or LSIL) and a 15% probability of high-grade lesion (HSIL or CC) (Table 5).

3.2.2. Clinical case No. 2

The patient designated VIII in Table 4 had an LSIL cytology result. HPV type 16 was detected. Colposcopy results revealed mild abnormalities in the epithelium of the cervix (Figure S5b). The epithelial scraping lipidome analysis revealed HSIL (Table 5). A biopsy of the cervix was performed. The histological diagnosis resembled HSIL (CIN II-III). The treatment was carried out. According to Table 5, the probability of high lesions (HSIL and CC) in this patient was also higher than low ones—60% and 40%, respectively, which is consistent with the histological diagnosis.

4. Discussion

It is known that lipid metabolism disorders are closely associated with carcinogenesis. Cholesterol esters and cholesterol accumulate in macrophages [25], which engulf dead cells. Previously, based on studies of kidney cancer, it was shown that cholesterol metabolism is disturbed in cancer cells [26]. Neutral sphingomyelinase decomposes sphingomyelins into ceramides and choline phosphate. It is involved in the activation of signals associated with inflammatory processes and acts as a tumor suppressor [27,28]. Phosphatidylcholines are associated with the carcinogenesis process [29,30]. Similarly, A.M. Porcari et al. demonstrated differences in levels of ceramides and sphingosine metabolites in normal cervical tissues and tissues with severe neoplasia, as well as the potential utility of these markers for cervical tissue differentiation [12].
Our final model has a higher specificity for CC and high predictive value for predicting cervicitis and cervical cancer, as well as an average value for predicting LSIL relative to the NILM, compared to models based on p16/Ki-67 proteins immunocytochemistry, microRNA, and HPV DNA [31]. It should be noted that the above-mentioned methods are based on compounds associated with HPV and neoplastic transformations, which makes it difficult to differentiate lesions with low and no viral load (LSIL/cervicitis/NILM). The model, based on information obtained by Raman spectroscopy, demonstrates high accuracy in the differentiation of no lesions/CIN1/CIN 2-3 lesions/CC; however, this method is more invasive since it requires tissue biopsy [32,33]. Our previous study demonstrated 88% sensitivity and 71% specificity of the OPLS-DA model used to classify mild neoplastic lesions of the cervix (LSIL) from severe precancerous and malignant diseases of the cervix (HSIL and SCC) [16]. This model was based on the lipid profile of cervical tissues (biopsy samples). The non-invasive approach proposed in this study has an even higher diagnostic potential (83% sensitivity and 88% specificity).
Models based on lipid markers in serum also have high diagnostic potential for differentiating NILM and SIL [34]. Nam M et al. conducted a comparative study of plasma lipidome in HSIL, CC, LSIL, and control groups by mass spectrometry and found a difference in plasma lipid profile for each of the four groups [35]. Khan I, Nam M, and Kwon M found a difference in the polar metabolite levels in blood plasma in cases with no neoplastic cervical lesions and CIN1, CIN 2/3, CC [36].

5. Conclusions

The developed model, based on the lipid profile obtained by HPLC-MS analysis of cervical epithelium scrapings, makes it possible to differentiate HSIL and CC from LSIL, cervicitis, and NILM. The lipidome diagnostic panel in patients with HSIL and cervical cancer is characterized by the predominance of lipid classes of phosphatidylcholines and sphingomyelin, as well as cholesterol esters. The non-invasive method used in the work demonstrated a high diagnostic value for severe lesions (HSIL+) with 83% sensitivity and 88% specificity. The diagnostic potential for the lipid profile models was higher than cytological examination, shown in two clinical cases. This suggests the potential applicability of the developed method in cervical cancer screening.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/metabo12090883/s1. Figure S1. Tandem mass spectra characteristic for the identified lipid markers. Figure S2. Characteristic chromatograms of samples from each clinical group obtained in (a) positive ion mode; (b) negative ion mode. Figure S3. Characteristic mass spectra of samples from each clinical group, obtained in positive ion mode by averaging over the analysis time. Figure S4. Characteristic mass spectra of samples from each clinical group, obtained in negative ion mode by averaging over the analysis time. Figure S5. The result of a colposcopy of a sample with 3% solution of acetic acid: (a) Clinical case # 1, (b) Clinical case # 2. Table S1. The resulting abundance-m/z table for compounds detected by LC/MS. Table S2. Identified lipid species with the corresponding ion mass and mass of the characteristic fragment ion. Table S3. Levels of each lipid in each sample after auto scaling normalization. Table S4. Box plot of levels of lipid markers in samples analyzed for each clinical case. Table S5. Description of each discriminative logistic regression model: lipid markers, their coefficients, standard deviation and confidence interval of coefficient value, Wald criteria, P-value. Supplementary File S1. R-scripts, using for processing, statistical analysis and model validation. Supplementary File S2. Summary and loading plot, S-plot and VIP-plot of OPLS-DA models, based on lipids in negative ion mode analysis.

Author Contributions

Data curation, A.T., D.A., K.G. and E.K.; Formal analysis, A.T. and N.S.; Funding acquisition, N.N., V.F. and G.S.; Investigation, V.C., E.K. and N.S.; Methodology, V.C., A.T., N.N. and V.F.; Project administration, N.S., N.N., V.F. and G.S.; Resources, D.A., K.G. and E.K.; Software, A.T. and V.C.; Supervision, V.C., N.N., V.F. and G.S.; Validation, D.A., K.G. and E.K.; Visualization, A.T. and D.A.; Writing—original draft, A.T., D.A., K.G. and N.S.; Writing—review and editing, V.C., N.N., V.F. and G.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the state assignment of the Ministry of Healthcare of the Russian Federation (Registration No. 121040600125-0).

Institutional Review Board Statement

The study was approved by the Ethical Committee of the National Medical Research Center for Obstetrics, Gynecology and Perinatology named after Academician V.I. Kulakov (protocol No. 4 of October 2019).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Data are contained within the Supplementary Material.

Acknowledgments

This work was supported by the state assignment of the Ministry of Healthcare of the Russian Federation. The authors are grateful to the laboratory for the collection and storage of biological material (Biobank) for providing tissue samples.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bruni, L.; Albero, G.; Serrano, B.; Mena, M.; Collado, J.J.; Gómez, D.; Muñoz, J.; Bosch, F.X.; Sanjosé, S. Human Papillomavirus and Related Diseases in the World. 2021. Available online: https://hpvcentre.net/statistics/reports/XWX.pdf (accessed on 22 October 2021).
  2. Lee, M.H.; Finlayson, S.J.; Gukova, K.; Hanley, G.; Miller, D.; Sadownik, L.A. Outcomes of Conservative Management of High Grade Squamous Intraepithelial Lesions in Young Women. J. Low. Genit. Tract Dis. 2018, 22, 212–218. [Google Scholar] [CrossRef] [PubMed]
  3. Alanbay, I.; Öztürk, M.; Fıratlıgil, F.B.; Karaşahin, K.E.; Yenen, M.C.; Bodur, S. Cytohistological discrepancies of cervico-vaginal smears and HPV status. Ginekol. Polska 2017, 88, 235–238. [Google Scholar] [CrossRef] [PubMed]
  4. Gupta, R.; Hariprasad, R.; Dhanasekaran, K.; Sodhani, P.; Mehrotra, R.; Kumar, N.; Gupta, S. Reappraisal of cytology-histology correlation in cervical cytology based on the recent American Society of Cytopathology guidelines (2017) at a cancer research centre. Cytopathology 2019, 31, 53–58. [Google Scholar] [CrossRef] [PubMed]
  5. Crothers, B.A.; Ghofrani, M.; Zhao, C.; Dodd, L.G.; Goodrich, K.; Husain, M.; Kurtycz, D.F.; Russell, D.K.; Shen, R.Z.; Souers, R.J.; et al. Low-Grade Squamous Intraepithelial Lesion or High-Grade Squamous Intraepithelial Lesion? Concordance Between the Interpretation of Low-Grade Squamous Intraepithelial Lesion and High-Grade Squamous Intraepithelial Lesion in Papanicolaou Tests: Results from the College of American Pathologists PAP Education Program. Arch. Pathol. Lab. Med. 2018, 143, 81–85. [Google Scholar] [CrossRef]
  6. Ouh, Y.-T.; Park, J.J.; Kang, M.; Kim, M.; Song, J.Y.; Shin, S.J.; Shim, S.-H.; Yoo, H.J.; Lee, M.; Lee, S.-J.; et al. Discrepancy between Cytology and Histology in Cervical Cancer Screening: A Multicenter Retrospective Study (KGOG 1040). J. Korean Med Sci. 2021, 36, e164. [Google Scholar] [CrossRef] [PubMed]
  7. Güzel, C.; Hoff, J.v.S.-V.; de Kok, I.M.; Govorukhina, N.I.; Boychenko, A.; Luider, T.M.; Bischoff, R. Molecular markers for cervical cancer screening. Expert Rev. Proteom. 2021, 18, 675–691. [Google Scholar] [CrossRef]
  8. Boers, A.; Wang, R.; van Leeuwen, R.W.; Klip, H.G.; de Bock, G.H.; Hollema, H.; van Criekinge, W.; de Meyer, T.; Denil, S.; van der Zee, A.G.J.; et al. Discovery of new methylation markers to improve screening for cervical intraepithelial neoplasia grade 2/3. Clin. Epigenetics 2016, 8, 2262. [Google Scholar] [CrossRef]
  9. del Pino, M.; Sierra, A.; Marimon, L.; Delgado, C.M.; Rodriguez-Trujillo, A.; Barnadas, E.; Saco, A.; Torné, A.; Ordi, J. CADM1, MAL, and miR124 Promoter Methylation as Biomarkers of Transforming Cervical Intrapithelial Lesions. Int. J. Mol. Sci. 2019, 20, 2262. [Google Scholar] [CrossRef]
  10. Duraipandian, S.; Mo, J.; Zheng, W.; Huang, Z. Near-infrared Raman spectroscopy for assessing biochemical changes of cervical tissue associated with precarcinogenic transformation. Analyst 2014, 139, 5379–5386. [Google Scholar] [CrossRef]
  11. Paraskevaidi, M.; Cameron, S.J.; Whelan, E.; Bowden, S.; Tzafetas, M.; Mitra, A.; Semertzidou, A.; Athanasiou, A.; Bennett, P.R.; MacIntyre, D.A.; et al. Laser-assisted rapid evaporative ionisation mass spectrometry (LA-REIMS) as a metabolomics platform in cervical cancer screening. eBioMedicine 2020, 60, 103017. [Google Scholar] [CrossRef]
  12. Porcari, A.M.; Negrão, F.; Tripodi, G.L.; Pitta, D.R.; Campos, E.A.; Montis, D.M.; Martins, A.M.A.; Eberlin, M.N.; Derchain, S.F.M. Molecular Signatures of High-Grade Cervical Lesions. Front. Oncol. 2018, 8, 99. [Google Scholar] [CrossRef] [PubMed]
  13. Brügger, B. Lipidomics: Analysis of the Lipid Composition of Cells and Subcellular Organelles by Electrospray Ionization Mass Spectrometry. Annu. Rev. Biochem. 2014, 83, 79–98. [Google Scholar] [CrossRef] [PubMed]
  14. Zhang, J.; Feider, C.L.; Nagi, C.; Yu, W.; Carter, S.A.; Suliburk, J.; Cao, H.S.T.; Eberlin, L.S. Detection of Metastatic Breast and Thyroid Cancer in Lymph Nodes by Desorption Electrospray Ionization Mass Spectrometry Imaging. J. Am. Soc. Mass Spectrom. 2017, 28, 1166–1174. [Google Scholar] [CrossRef] [PubMed]
  15. Snaebjornsson, M.T.; Janaki-Raman, S.; Schulze, A. Greasing the Wheels of the Cancer Machine: The Role of Lipid Metabolism in Cancer. Cell Metab. 2019, 31, 62–76. [Google Scholar] [CrossRef]
  16. Starodubtseva, N.L.; Chagovets, V.V.; Nekrasova, M.E.; Nazarova, N.M.; Tokareva, A.O.; Bourmenskaya, O.V.; Attoeva, D.I.; Kukaev, E.N.; Trofimov, D.Y.; Frankevich, V.E.; et al. Shotgun Lipidomics for Differential Diagnosis of HPV-Associated Cervix Transformation. Metabolites 2022, 12, 503. [Google Scholar] [CrossRef] [PubMed]
  17. Iurova, M.V.; Chagovets, V.V.; Pavlovich, S.V.; Starodubtseva, N.L.; Khabas, G.N.; Chingin, K.S.; Tokareva, A.O.; Sukhikh, G.T.; Frankevich, V.E. Lipid Alterations in Early-Stage High-Grade Serous Ovarian Cancer. Front. Mol. Biosci. 2022, 9. [Google Scholar] [CrossRef] [PubMed]
  18. Tonoyan, N.M.; Chagovets, V.V.; Starodubtseva, N.L.; Tokareva, A.O.; Chingin, K.; Kozachenko, I.F.; Adamyan, L.V.; Frankevich, V.E. Alterations in lipid profile upon uterine fibroids and its recurrence. Sci. Rep. 2021, 11, 1–11. [Google Scholar] [CrossRef]
  19. Koelmel, J.P.; Kroeger, N.M.; Ulmer, C.Z.; Bowden, J.A.; Patterson, R.E.; Cochran, J.A.; Beecher, C.W.W.; Garrett, T.J.; Yost, R.A. LipidMatch: An automated workflow for rule-based lipid identification using untargeted high-resolution tandem mass spectrometry data. BMC Bioinform. 2017, 18, 1–11. [Google Scholar] [CrossRef]
  20. Sud, M.; Fahy, E.; Cotter, D.; Brown, A.; Dennis, E.A.; Glass, C.K.; Merrill, A.; Murphy, R.C.; Raetz, C.R.H.; Russell, D.; et al. LMSD: LIPID MAPS structure database. Nucleic Acids Res. 2006, 35, D527–D532. [Google Scholar] [CrossRef]
  21. Tokareva, A.O.; Chagovets, V.V.; Kononikhin, A.S.; Starodubtseva, N.L.; Nikolaev, E.N.; Frankevich, V.E. Normalization methods for reducing interbatch effect without quality control samples in liquid chromatography-mass spectrometry-based studies. Anal. Bioanal. Chem. 2021, 413, 3479–3486. [Google Scholar] [CrossRef]
  22. R Core Team. R: A Language and Environment for Statistical Computing. Version 4.1.1; R Foundation for Statistical Computing: Vienna, Austria, 2021. [Google Scholar]
  23. RStudio Team. RStudio: Integrated Development for R. Version 1.1.463; RStudio, Inc.: Boston, MA, USA, 2016. [Google Scholar]
  24. Tokareva, A.O.; Chagovets, V.V.; Kononikhin, A.S.; Starodubtseva, N.L.; Nikolaev, E.N.; Frankevich, V.E. Comparison of the effectiveness of variable selection method for creating a diagnostic panel of biomarkers for mass spectrometric lipidome analysis. J. Biol. Mass Spectrom. 2021, 56, e4702. [Google Scholar] [CrossRef] [PubMed]
  25. Chistiakov, D.A.; Bobryshev, Y.V.; Orekhov, A.N. Macrophage-mediated cholesterol handling in atherosclerosis. J. Cell. Mol. Med. 2016, 20, 17–28. [Google Scholar] [CrossRef]
  26. Gebhard, R.L.; Clayman, R.V.; Prigge, W.F.; Figenshau, R.; Staley, N.A.; Reesey, C.; Bear, A. Abnormal cholesterol metabolism in renal clear cell carcinoma. J. Lipid Res. 1987, 28, 1177–1184. [Google Scholar] [CrossRef]
  27. Shamseddine, A.A.; Airola, M.V.; Hannun, Y.A. Roles and regulation of neutral sphingomyelinase-2 in cellular and pathological processes. Adv. Biol. Regul. 2014, 57, 24–41. [Google Scholar] [CrossRef] [PubMed]
  28. Hannun, Y.A.; Obeid, L.M. Sphingolipids and their metabolism in physiology and disease. Nat. Rev. Mol. Cell Biol. 2017, 19, 175–191. [Google Scholar] [CrossRef] [PubMed]
  29. Zeisel, S.H.; Canty, D.J. Choline phospholipids: Molecular mechanisms for human diseases: A meeting report. J. Nutr. Biochem. 1993, 4, 258–263. [Google Scholar] [CrossRef]
  30. Podo, F. Tumour phospholipid metabolism. NMR Biomed. 1999, 12, 413–439. [Google Scholar] [CrossRef]
  31. Frega, A.; Pavone, M.; Sesti, F.; Leone, C.; Bianchi, P.; Cozza, G.; Colombrino, C.; Lukic, A.; Marziani, R.; Sanctis, L.D.E.; et al. Sensitivity and specificity values of high-risk HPV DNA, p16/ki-67 and HPV mRNA in young women with atypical squamous cells of undetermined significance (ASCUS) or low-grade squamous intraepithelial lesion (LSIL). Eur. Rev. Med. Pharmacol. Sci. 2019, 23, 10672–10677. [Google Scholar] [CrossRef]
  32. Traynor, D.; Duraipandian, S.; Bhatia, R.; Cuschieri, K.; Tewari, P.; Kearney, P.; D’Arcy, T.; O’Leary, J.J.; Martin, C.M.; Lyng, F.M. Development and Validation of a Raman Spectroscopic Classification Model for Cervical Intraepithelial Neoplasia (CIN). Cancers 2022, 14, 1836. [Google Scholar] [CrossRef]
  33. Wang, J.; Zheng, C.-X.; Ma, C.-L.; Zheng, X.-X.; Lv, X.-Y.; Lv, G.-D.; Tang, J.; Wu, G.-H. Raman spectroscopic study of cervical precancerous lesions and cervical cancer. Lasers Med Sci. 2021, 36, 1855–1864. [Google Scholar] [CrossRef]
  34. Cheng, F.; Wen, Z.; Feng, X.; Wang, X.; Chen, Y. A serum lipidomic strategy revealed potential lipid biomarkers for early-stage cervical cancer. Life Sci. 2020, 260, 118489. [Google Scholar] [CrossRef] [PubMed]
  35. Nam, M.; Seo, S.-S.; Jung, S.; Jang, S.Y.; Lee, J.; Kwon, M.; Khan, I.; Ryu, D.H.; Kim, M.K.; Hwang, G.-S. Comparable Plasma Lipid Changes in Patients with High-Grade Cervical Intraepithelial Neoplasia and Patients with Cervical Cancer. J. Proteome Res. 2020, 20, 740–750. [Google Scholar] [CrossRef] [PubMed]
  36. Khan, I.; Nam, M.; Kwon, M.; Seo, S.-S.; Jung, S.; Han, J.S.; Hwang, G.-S.; Kim, M.K. LC/MS-Based Polar Metabolite Profiling Identified Unique Biomarker Signatures for Cervical Cancer and Cervical Intraepithelial Neoplasia Using Global and Targeted Metabolomics. Cancers 2019, 11, 511. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. The structure of the final diagnostic model. Panelx D1/D2 indicate lipids included in the logistic regression to discriminate diagnosis D1 from diagnosis D2 in the x ion mode; Modelx D1/D2 indicate model based on logistic regression to discriminate diagnoses D1 and D2 in the x ion mode, Thrx D1/D2 indicate a threshold of a logistic regression result to set a D1 or D2 diagnosis.
Figure 1. The structure of the final diagnostic model. Panelx D1/D2 indicate lipids included in the logistic regression to discriminate diagnosis D1 from diagnosis D2 in the x ion mode; Modelx D1/D2 indicate model based on logistic regression to discriminate diagnoses D1 and D2 in the x ion mode, Thrx D1/D2 indicate a threshold of a logistic regression result to set a D1 or D2 diagnosis.
Metabolites 12 00883 g001
Figure 2. Box plot of lipid levels detected in positive ion mode and included in the discriminative models. CE—cholesterol ester; Cer-NDS—no-hydroxy fatty acid dihydrosphingosineceramide; CerP—ceramide-1-phosphate; LPC—lysophosphatidylcholine; LPC O—plasmanyllysophosphatidylcholine; MG—monoacylglycerol; PC—phosphatidylcholine; PE—phosphatidylethanolamine, PE P—plasmenylphosphatidylethanolamine; PEtOH—phosphatidylethanol; SM—sphingomyelin; TG—triacylglycerol; TG O—plasmanyltriacylglycerol.
Figure 2. Box plot of lipid levels detected in positive ion mode and included in the discriminative models. CE—cholesterol ester; Cer-NDS—no-hydroxy fatty acid dihydrosphingosineceramide; CerP—ceramide-1-phosphate; LPC—lysophosphatidylcholine; LPC O—plasmanyllysophosphatidylcholine; MG—monoacylglycerol; PC—phosphatidylcholine; PE—phosphatidylethanolamine, PE P—plasmenylphosphatidylethanolamine; PEtOH—phosphatidylethanol; SM—sphingomyelin; TG—triacylglycerol; TG O—plasmanyltriacylglycerol.
Metabolites 12 00883 g002
Figure 3. Box plot of lipid levels included in the discriminative models in the negative ion mode. Cer-NDS—non-hydroxy fatty acid dihydrosphingosineceramide; Cer-NP—non-hydroxy fatty acid phytosphingosineceramide; CL—cardiolipin; HexCer-AP—alpha-hydoxy fatty acid phytosphingosineglucosylceramide; LPC—lysophosphatidylcholine; OxCL—oxidized cardiolipin; OxPC—oxidized phosphatidylcholine; PC—phosphatidylcholine; PC O—plasmanylphosphatidylcholine; PE P—plasmenylphosphatidylethanolamine; SM—sphingomyelin.
Figure 3. Box plot of lipid levels included in the discriminative models in the negative ion mode. Cer-NDS—non-hydroxy fatty acid dihydrosphingosineceramide; Cer-NP—non-hydroxy fatty acid phytosphingosineceramide; CL—cardiolipin; HexCer-AP—alpha-hydoxy fatty acid phytosphingosineglucosylceramide; LPC—lysophosphatidylcholine; OxCL—oxidized cardiolipin; OxPC—oxidized phosphatidylcholine; PC—phosphatidylcholine; PC O—plasmanylphosphatidylcholine; PE P—plasmenylphosphatidylethanolamine; SM—sphingomyelin.
Metabolites 12 00883 g003
Table 1. Contingency table for diagnoses of cervix lesions established by cytology and histology. Number of patients and respective accuracy (%). The accuracy of each diagnosis is indicated in parentheses, calculated as the ratio of the number of true cases of any cytological diagnosis to the total number of cases with such a cytological diagnosis.
Table 1. Contingency table for diagnoses of cervix lesions established by cytology and histology. Number of patients and respective accuracy (%). The accuracy of each diagnosis is indicated in parentheses, calculated as the ratio of the number of true cases of any cytological diagnosis to the total number of cases with such a cytological diagnosis.
Histological Diagnosis Cytological Diagnosis
NILMCervicitisLSILHSILCancer
NILM8 (29)0000
Cervicitis128 (53)720
LSIL3518 (64)60
HSIL42310 (45)0
Cancer100418 (100)
Table 2. Parameters of discriminative models based on logistic regression obtained during cross-validation.
Table 2. Parameters of discriminative models based on logistic regression obtained during cross-validation.
Ion ModeModelAUC (CI AUC)Thr. 1Sens. 2, %Spec. 3, %
Positive modeNILM/cervicitis0.89 (0.79–1.00) 0.9710067
NILM/LSIL0.93 (0.85–1.00) 0.9610057
NILM/HSIL0.68 (0.46–0.89)0.8110044
NILM/cancer0.94 (0.87–1.00) 0.9910080
Cervicitis/LSIL0.81 (0.70–0.93)0.137495
Cervicitis/HSIL0.63 (0.46–0.81)0.466173
Cervicitis/Cancer0.90 (0.76–1.00)0.990085
LSIL/HSIL0.68 (0.52–0.84) 0.376775
LSIL/Cancer0.89 (0.71–0.98) 0.019188
HSIL/Cancer0.95 (0.85–1.00) 0.9910086
Negative modeNILM/cervicitis0.81 (0.56–1.00) 0.49 93 86
NILM/LSIL0.78 (0.56–1.00) 0.50 91 71
NILM/HSIL0.81 (0.66–1.00) 0.48 89 75
NILM/cancer0.75 (0.40–1.00) 0.99 87 63
Cervicitis/LSIL0.65 (0.51–0.79) 0.51 67 63
Cervicitis/HSIL0.82 (0.67–0.95) 0.15 64 87
Cervicitis/Cancer0.94 (0.86–1.00) 0.01 92 96
LSIL/HSIL0.64 (0.48–0.81) 0.39 56 73
LSIL/Cancer0.89 (0.73–0.99) 0.99 95 89
HSIL/Cancer0.77 (0.63–0.92) 0.27 76 76
1 Thr.—threshold; 2 Sens.—sensitivity; 3 Spec.—specificity.
Table 3. Lipidome diagnostic contingency table for differential diagnosis of the severity of neoplastic lesions. Values in parentheses indicate the accuracy of each diagnosis.
Table 3. Lipidome diagnostic contingency table for differential diagnosis of the severity of neoplastic lesions. Values in parentheses indicate the accuracy of each diagnosis.
NILMCervicitisLSILHSILCancer
NILM 8 (0.89) 0 0 0 0
Cervicitis0 23 (0.82) 2 4 0
LSIL 0 5 23 (0.74) 4 0
HSIL 1 0 6 12 (0.57) 0
Cancer0 0 0 1 22 (1.00)
Table 4. Comparison of lipidome, cytological and histological analysis of several patients. Bold type indicates patients for whom histology and lipidomic diagnosis coincided.
Table 4. Comparison of lipidome, cytological and histological analysis of several patients. Bold type indicates patients for whom histology and lipidomic diagnosis coincided.
PatientsHistologyCytologyLipidome
IHSILCervicitisLSIL
IICervicitisHSILCervicitis
IIIHSILNILMLSIL
IVHSILLSILHSIL
VHSILNILMCervicitis
VIHSILNILMLSIL
VIIHSILLSILCervicitis
VIIIHSILLSILHSIL
Table 5. Distribution of votes for each diagnosis in clinical cases for patients II and VIII.
Table 5. Distribution of votes for each diagnosis in clinical cases for patients II and VIII.
PatientNILMCervicitisLSILHSILCancerProbability of
Low Lesions, %
Probability of
High Lesions, %
II386218515
VIII431754060
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tokareva, A.; Chagovets, V.; Attoeva, D.; Starodubtseva, N.; Nazarova, N.; Gusakov, K.; Kukaev, E.; Frankevich, V.; Sukhikh, G. Non-Invasive Differential Diagnosis of Cervical Neoplastic Lesions by the Lipid Profile Analysis of Cervical Scrapings. Metabolites 2022, 12, 883. https://doi.org/10.3390/metabo12090883

AMA Style

Tokareva A, Chagovets V, Attoeva D, Starodubtseva N, Nazarova N, Gusakov K, Kukaev E, Frankevich V, Sukhikh G. Non-Invasive Differential Diagnosis of Cervical Neoplastic Lesions by the Lipid Profile Analysis of Cervical Scrapings. Metabolites. 2022; 12(9):883. https://doi.org/10.3390/metabo12090883

Chicago/Turabian Style

Tokareva, Alisa, Vitaliy Chagovets, Djamilja Attoeva, Natalia Starodubtseva, Niso Nazarova, Kirill Gusakov, Eugenii Kukaev, Vladimir Frankevich, and Gennady Sukhikh. 2022. "Non-Invasive Differential Diagnosis of Cervical Neoplastic Lesions by the Lipid Profile Analysis of Cervical Scrapings" Metabolites 12, no. 9: 883. https://doi.org/10.3390/metabo12090883

APA Style

Tokareva, A., Chagovets, V., Attoeva, D., Starodubtseva, N., Nazarova, N., Gusakov, K., Kukaev, E., Frankevich, V., & Sukhikh, G. (2022). Non-Invasive Differential Diagnosis of Cervical Neoplastic Lesions by the Lipid Profile Analysis of Cervical Scrapings. Metabolites, 12(9), 883. https://doi.org/10.3390/metabo12090883

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop