- 1National Medical Research Center for Obstetrics Gynecology and Perinatology Named after Academician V.I. Kulakov of the Ministry of Healthcare of Russian Federation, Moscow, Russia
- 2Department of Chemical Physics, The Moscow Institute of Physics and Technology, Moscow, Russia
- 3Laboratory of Translational Cellular And Molecular Biomedicine, National Research Tomsk State University, Tomsk, Russia
- 4Cancer Research Institute, Tomsk National Research Medical Center, Russian Academy of Sciences, Tomsk, Russia
- 5Laboratory of Genetic Technologies, Siberian State Medical University, Tomsk, Russia
- 6Institute of Transfusion Medicine and Immunology, Mannheim Faculty of Medicine, University of Heidelberg, Heidelberg, Germany
- 7German Red Cross Blood Service Baden-Württemberg–Hessen, Mannheim, Germany
- 8Laboratory of Translational Medicine, Siberian State Medical University, Tomsk, Russia
Introduction: Immunometabolism is essential factor of tumor progression, and tumor-associated macrophages are characterized by substantial changes in their metabolic status. In this study for the first time, we applied targeted amino acid LC-MS/MS analysis to compare amino acid metabolism of circulating monocytes isolated from patients with breast, ovarian, lung, and colorectal cancer.
Methods: Monocyte metabolomics was analyzed by liquid chromatography-mass spectrometry/mass spectrometry (LC-MS/ MS) analysis of amino acid extracts. The targeted analysis of 26 amino acids was conducted by LCMS/MS on an Agilent 6460 triple quadrupole mass spectrometer equipped with an electrospray ionization source and an Agilent 1260 II liquid chromatograph.
Results: Comparison of monocytes of cancer patients with monocytes of healthy control individuals demonstrated that in breast cancer most pronounced changes were identified for tryptophan (AUC = 0.76); for ovarian cancer, aminobutyric acid was significantly elevated (AUC= 1.00); for lung cancer significant changes we indented for citrulline (AUC = 0.70). In order to identify key amino acids that are characteristic for monocytes in specific cancer types, we compared each individual cancer with other 3 types of cancer. We found, that aspartic acid and citrulline are specific for monocytes of patients with colorectal cancer (p<0.001, FC = 1.40 and p=0.003, FC = 1.42 respectively). Citrulline, sarcosine and glutamic acid are ovarian cancer-specific amino acids (p = 0.003, FC = 0.78, p = 0.003, FC = 0.62, p = 0.02, FC = 0.78 respectively). Glutamine, methionine and phenylalanine (p = 0.048, FC = 1.39. p = 0.03, FC = 1.27 and p = 0.02, FC = 1.41) are lung cancer-specific amino acids. Ornithine in monocytes demonstrated strong positive correlation (r = 0.63) with lymph node metastasis incidence in breast cancer patients. Methyl histidine and cysteine in monocytes had strong negative correlation with lymph node metastasis in ovarian cancer patients (r = -0.95 and r = -0.95 respectively). Arginine, citrulline and ornithine have strong negative correlation with tumor size (r = -0.78, citrulline) and lymph node metastasis (r = -0.63 for arginine and r = -0.66 for ornithine).
Discussion: These alterations in monocyte amino acid metabolism can reflect the reaction of systemic innate immunity on the growing tumor. Our data indicate that this metabolic programming is cancer specific and can be inhibiting cancer progression. Cancer-specific differences in citrulline, as molecular link between metabolic pathways and epigenetic programing, provide new option for the development and validation of anti-cancer therapies using inhibitors of enzymes catalyzing citrullination.
1 Introduction
Monocytes in the adults represent a dynamic population of precursors of key innate immune cells in tumor microenvironment - tumor associate macrophages (TAMs). Circulating monocytes are differentiated constantly from bone marrow progenitors. Monocytes can spend 1.5 to 4 days in the bloodstream, and can be recruited to the tissues to give rise to macrophages (1–5). Monocyte-derived macrophages in tumor microenvironments can support proliferation of cancer cells, tumor angiogenesis and cancer cells migration and extravasation into blood and lymphatic vessels (6). Different types of monocyte-derived TAMs can have diverse effects on tumor development, growth, angiogenesis, and metastatic progression, as well as suppression of adaptive immune responses (7). In majority of cancers, including breast and lung cancer, both total amount of TAMs as well as their M2 polarization positively correlate with primary tumor growth, metastasis and poor prognosis and decreased survival of patients (8).
Tumor progression is associated with the initiation of a cascade of immune reactions. These pathological processes are linked to changes in both the overall level of peripheral monocytes and their transcriptome and metabolome (9–14). However, majority of the studies were focus on transcriptome.
Monocytes can encounter cancer-related factors in the bloodstream, potentially acquiring programs detrimental to the patient, which are then transferred to tumor tissue, affecting the differentiation of tumor-associated macrophages (TAMs), dendritic cells (DC), and myeloid-derived suppressor cells (MDSC) (12, 15–17). Recent studies utilizing high-throughput methods have shown that transcriptional changes in peripheral blood monocytes can serve as diagnostic, predictive, and prognostic biomarkers in various cancers, including renal, colorectal, breast, cervical, skin, thyroid, hepatocellular, and lung cancers (3, 12, 16–20). The study of the unique features of amino acid metabolism in monocytes across different cancer types deserves special attention. Alterations in amino acid levels can profoundly impact the development of an effective immune response (21–24) and may disrupt the migration, division, and maturation of immune cells. Amino acids play a crucial role in regulating several pathways within immune cells, including mTOR signaling and nitric oxide (NO) production (23, 25). Competition for metabolites and signaling interactions among host immune cells and pathogens can influence disease development (22).
Advancements in analytical techniques, particularly liquid chromatography-mass spectrometry (LC-MS), have led to the emergence of metabolomics (26–28). This field is actively utilized to identify metabolic peculiarities in oncological processes, uncover pathogenesis mechanisms, and identify new therapeutic targets (26, 29). Unlike TAMs, DCs, and MDSCs, which require sophisticated isolation techniques from heterogeneous tumor tissues without compromising highly adherent subpopulations, monocytes are more amenable to examination. Understanding the heterogeneity of monocytes and their roles at various stages of cancer progression is paramount for investigating how cancer cells evade immune surveillance and elude immune cell aggression. Intervening in the pathways and mechanisms involved in monocyte and TAMs reprogramming holds promising potential as a novel therapeutic strategy for personalized cancer immunotherapy (14). We and others have previously identified that monocytes undergo cancer-specific programming that can be detected on the level of phenotype (surface markers) or on the level transcriptome (3, 30, 31). In breast cancer, CD14lowCD16+HLA-DR+ monocytes were indicative for the good response to neoadjuvant chemotherapy (3). In colon, but not in rectal cancer, overexpression of activator of glycolysis PFKFB3 in monocytes, was correlated with tumor elimination after therapy (30). Metabolism of immune cells is critical for the ability to stimulate or suppress tumor growth. However, despite the monocytes are massively recruited into growing tumor, their metabolic status is almost neglected.
In light of these considerations, our study was designed to analyze the monocyte metabolome in patients with tumors of diverse origins and to examine correlation of cancer-specific amino acid profiles in circulating monocytes with clinical parameters of tumor progression.
2 Materials and methods
2.1 Study design
The prospective multicenter study included 43 patients with cancer of different localizations (invasive breast carcinoma of no special type, colon and rectal adenocarcinoma, high-grade serous ovarian carcinoma, and lung adenocarcinoma). Patients with breast cancer (n=11) and colorectal cancer (n=16) were treated in the Cancer Research Institute, Tomsk National Research Medical Centre (Tomsk, Russia). Patients with ovarian cancer (n=6) and lung cancer (n=10) were treated in both Cancer Research Institute and Tomsk Regional Oncology Center (Tomsk, Russia). All patients lived in the same geographic area. Patients had no acute pathologies, no infectious disorders, and did not have a history of any other types of cancer in addition to the cancer. Patients did not receive any anticancer treatment prior to sample collection. Inclusion criteria were: primary cancer incident; no history other types of cancers; no anticancer treatment. Exclusion criteria were: acute or chronic infections disorders; pregnancy; alcohol and drug dependence; early-onset of cancer; terminal stage of cancer; cachexia; asthma or autoimmune disorders.
Healthy volunteers (5 women and 5 men) were enrolled in this study as a control group. The inclusion criteria for the healthy cohort were as follows: (a) age from 40 to 77 years, (b) no active medical conditions, (c) not taking immunomodulatory medication (over the counter or prescription) within 30 days of study, and (d) no current or past history of an oncology disease.
Detailed demographic and clinical data of 53 participants included in this study are included in the Supplementary Table 1 in Supplementary Material.docx.
The study received approval from the Local Committee for Medical Ethics and was conducted in accordance with the guidelines of the Declaration of Helsinki and the International Conference on Harmonisation’s Good Clinical Practice Guidelines (ICH GCP). Written informed consent was obtained from all subjects, including the patients/participants, who willingly provided their consent to participate in the study.
2.2 Monocyte isolation
Peripheral whole-blood samples were collected from both healthy volunteers (n=10) and patients with cancer of different localizations (n=43). Two millions of monocytes were isolated from 10 ml of blood of a patient. The isolation of monocytes from peripheral blood was performed using density gradients followed by positive magnetic selection with CD14+ MACS beads (#130-050-201, Miltenyi Biotech, Germany) as previously described (30). The isolation quality of monocytes was determined by flow cytometry. Monocyte samples were obtained by positive magnetic separation via СD14 microbeads. Peripheral blood mononuclear cells were the first step in the process. After enrichment on a density gradient, we observed a minimal number of granulocytes in the samples compared to the baseline condition (Supplementary Figures 1A and Supplementary Figure 1B in Supplementary material.docx). Magnetic separation led to enrichment of monocyte content (Supplementary Figure 1C in Supplementary Material.docx). Detection by flow cytometry showed a purity of CD14 monocyte isolation of more than 98% (Supplementary Figures 1D-G in Supplementary Material.docx). After monocyte isolation, the samples were washed twice with DPBS without calcium and magnesium at 300 g for 7 minutes. Subsequently, the cell precipitate obtained was dried under a nitrogen atmosphere and stored at -80°C until further experiments (32).
2.3 Amino acid extraction
Monocyte metabolomics was analyzed by liquid chromatography-mass spectrometry/mass spectrometry (LC-MS/MS) analysis of amino acid extracts (32). A 480 µL chloroform-methanol mixture (2:1, v/v) was added at 4°C to the monocyte sample, containing 2*106 cells; sonicated for 10 minutes, following addition of 150 µL of water; centrifuged at 13000 g for 5 minutes at room temperature. The upper aqueous-methanol layer (150 µL) was collected and dried under a stream of nitrogen for 30 minutes at 60°C following addition of 200 µL of 0.1 M hydrochloric acid in butanol for derivatization of amino acids. The solution was mixed for 3 minutes, centrifuged at 13000 g for 15 seconds at room temperature and kept at 60°C for 15 minutes to carry out the derivatization reaction. Afterward, centrifugation was performed at 13000 g for 15 seconds at room temperature. The sample was dried under a stream of nitrogen for 30 minutes at 60°C and reconstituted in 300 µL of acetonitrile/water solution (1:1, v/v). The mixture was stirred for 5 minutes. Centrifugation was carried out at 13000 g for 15 seconds at room temperature for 10 minutes. Finally, 200 µL of the resulting sample was transferred to a vial with an insert for further LC-MS/MS analysis.
2.4 Targeted amino acid LC-MS/MS analysis
The targeted analysis of 26 amino acids was conducted by LC-MS/MS on an Agilent 6460 triple quadrupole mass spectrometer (Agilent) equipped with an electrospray ionization source and an Agilent 1260 II liquid chromatograph (Agilent). The system featured a binary high-pressure pump, column thermostat, and an autosampler with 122 vials on an Amino acids-HPLC Column (Jasem) at a temperature of 30°C. The elution mode is provided in Supplementary Table 2 in Supplementary Material.docx, with a flow rate of 0.15 mL/min and a sample volume of 3 µL. Mass spectrometric analysis was performed in positive ion mode at a drying gas temperature and flow rate of 150°C and 10 L/min, respectively. The nebulizer gas pressure was 2.76 bar, the curtain gas temperature and flow rate were 400°C and 10 L/min, respectively, and the voltage was set at 2 kV. MRM transitions for detection are listed in Supplementary Table 3 in Supplementary Material.docx.
No internal standards were used in this study. The reproducibility of analyte extraction was controlled by preparing each sample in triplicate. The difference between preparations did not exceed 15%. Equal amounts of all samples were pooled as a quality control sample to monitor the stability of the LC-MS system and the reproducibility of the analysis. The difference between QC injections did not exceed 10%. The peak areas of amino acids are provided in the Supplementary_Data.xlxs, the TIC and EIC of analite were demonstrated in Supplementary Figure 2 in the Supplementary Material.docx.
2.5 Statistical analysis
The statistical analysis of experimental data was performed using scripts written in the R language (33) within RStudio (34).
Comparison of amino acid levels between control and cancer groups and comparison between one cancer group and other three was performed by the non-parametric Wilcoxon-Mann-Whitney test. For describing quantitative data, medians (Me) and quartiles Q1 and Q3 were used. The significance level (uncorrected p-value) was set to 0.05. In addition, powers of tests were calculated based on effect size (35) and included in Supplementary Table 4 in the Supplementary Material.docx. Correlations between parameters describing cancer progression (tumor size, number of LNM and hematogenous metastasis) and amino acid levels were assessed using the Spearman rank test, with a significance threshold of p = 0.05.
Logistic regression models for classifying “cancer” vs. “control” were built through stepwise selection of amino acids that exhibited statistically significant differences until there was a reduction in the Akaike Information Criterion (AIC) value. Model quality was evaluated using leave-one-out cross-validation to choose optimal thresholds, sensitivity, and specificity based on maximizing the sum of sensitivity and specificity.
The analysis of the involvement of the cancer specific amino acids in metabolic pathways was conducted using the MetaboAnalyst 5.0 resource (https://www.metaboanalyst.ca/home.xhtml). The pathway impact is calculated as the sum of significance scores of the corresponding metabolites, normalized by the sum of significance scores of all metabolites in each pathway. Enrichment analysis of metabolic pathways was performed using Over Representation Analysis (ORA) employing the hypergeometric test.
3 Results
3.1 Comparison cancer patients to healthy individuals.
The levels of asparagine, aspartic acid, and tryptophan statistically significantly decreased in patients with breast cancer compared to control group (women) (see Supplementary Figure 3 and Supplementary Table 5 in Supplementary Material.docx). The levels of citrulline, sarcosine and aminobutyric acid statistically significant decrease in patients with ovarian cancer compared to control group (women) (see Supplementary Figure 4 and Supplementary Table 6 in Supplementary Material.docx). The level of citrulline statistically significant decrease with the development of lung cancer compared to control (Supplementary Figure 5 and Supplementary Table 7 in Supplementary Material.docx). In case of colorectal cancer, no statistically significant differences in amino acid profile were identified compared to the healthy donor controls.
Logistic regression models were created for each cancer type classification versus appropriate control group. The formulas are presented below:
y is the dependent variable representing the probability of cancer, and Ix is the intensity of the amino acid peak X.
The diagnostic models exhibited an AUC of at least 0.80 and a sensitivity exceeding 85% (see Figure 1 and Table 1).
In breast and ovarian cancers the most pronounced disturbances in amino acid pathways included arginine biosynthesis and alanine, aspartate, and glutamate metabolism (Supplementary Figure 6, Supplementary Tables 8, 9 in Supplementary Material in Supplementary Material.docx).
3.2 Key differences in amino acid profiles in monocytes of patients identified by direct comparison of monocytes profiles between cancer groups.
Next, we analyzed specific changes in amino acid profiles by comparing mass spectrometry data for monocytes isolated from patients with 4 different cancer types. The comparison was performed by analysis of each individual cancer type with the group of three other cancer types. Statistically significant increase of aspartic acid and citrulline levels were detected in case of colorectal cancer (Figure 2A; Supplementary Table 10 in Supplementary Material.docx), statistically significant decrease of citrulline, glutamic acid and sarcosine were detected in case of ovarian cancer (Figure 2B; Supplementary Table 11 in Supplementary Material.docx), statistically significant increase of glutamine, methionine and phenylalanine levels detected in case of lung cancer (Figure 2C; Supplementary Table 12 in Supplementary Material.docx). In case of breast cancer statistically significant differences were not detected.
Figure 2 Boxplot of cancer-specific amino acids levels for: (A) colorectal cancer; (B) ovarian cancer; (C) lung cancer. *p (without multiply correction)<0.01, **p (without multiply correction)<0.001. Points represent the outliers.
3.3 Correlation of amino acid levels in peripheral blood monocytes correlation with clinical parameters of cancer
Next, we have analyzed whether amino acids profiles in monocytes of cancer patients correlate with the clinical parameters of disease progression: tumor size, number of lymph node metastasis and number of hematogenous metastasis. In breast cancer, a strong positive correlation was observed between the levels of ornithine and the number of LNM (r=0.63) (see Figure 3A). In ovarian cancer, a very strong negative correlation was found between the number of LNM and the levels of methyl-histidine (r=-0.95) and cystine (r=-0.95) (see Figure 3B). Citrulline, arginine and ornithine levels in monocytes of patients with lung cancer negatively correlated with tumor size and LNM (see Figure 3C).
Figure 3 Amino acids exhibiting statistically significant correlations with cancer progression in (A) breast cancer, (B) ovarian cancer, and (C) lung cancer. T, tumor size; N, lymph node metastases.
4 Discussion
In this study, we undertook an effort to unveil the tumor-induced reprogramming of amino acid metabolism in peripheral blood monocytes. We have previously demonstrated that samples suitable for analysis should contain a minimum of 2 × 106 cells (32).
The observed identified differences in monocyte amino acid profiles among patients with different tumors (breast, ovarian, lung, and colorectal cancer) could reflect the long-term impact of tumor cells developing tumor on the metabolism of body in general and, in particular, the immune system (3, 12, 16–21). In cancer patients, the systemic control of monocytes is likely can be achieved through various factors, such as blood cytokines, chemokines, inflammation mediators, exosomes, and lipid and carbohydrate metabolites produced by the tumor (3, 12). Consequently, tumors have the potential to influence both the amount and phenotype of circulating monocytes. Subsequently modified monocyte subpopulations may limit the effectiveness of tumor therapy (36–38). Nevertheless, the role of monocytes as regulators of the tumor’s response to chemotherapeutic agents remains extremely limited (39).
A notable strength of this study lies in the meticulous patient selection and the collection of samples prior to any treatment, effectively mitigating the influence of therapy. It is well-established that various treatment modalities, such as neoadjuvant chemotherapy, surgical procedures, chemotherapy, and radiation therapy, can all exert an impact on blood monocyte levels and their cellular composition (5). Elevated monocyte levels can stem from either increased mobilization from the bone marrow or heightened monopoiesis, both of which have been observed in cancer patients (12). Elevated serum levels of cytokines induced by cancer tumors, combined with low-grade systemic inflammation, may contribute to the distant reprogramming of myelopoiesis (40). Importantly, in several cancer types, patients with higher blood monocyte counts have exhibited a correlation with a less favorable disease prognosis (41–45).
Immune cells face a deprivation of essential amino acids such as tryptophan, cysteine, and arginine (46). In this study, the amino acid profile of peripheral monocytes was dramatically substantially altered in patients with various cancers during the oncogenic process.
In the case of lung and breast cancer, arginine, ornithine and citrulline exhibit strong correlations (r>0.6, p<0.05) with the number of LNM and tumor size. In cancer, arginine metabolism violation is related to immune evasion (promotion of M2-type macrophages), impaired NO synthesis, and increased polyamine production (stimulation of cancer cell proliferation) (47–49). The hypoxic TME increases HIF-1α expression in tumor cells, which in turn stimulate inducible nitric oxide synthase (iNOS) and arginase (ARG1 and ARG2) synthesis (50–52). Enhanced arginase expression is a hallmark of many tumors, such as head and neck, kidney, breast, hepatocellular, and prostate cancers (53–57). Furthermore, anti-inflammatory cytokines released by tumors, such as IL-4, IL-10, and TGF-β1, trigger the activation of ARGs while suppressing iNOS in macrophages (58–60). This process fosters the development of pro-tumoral M2-like macrophages, characterized by heightened ARG activity, reduced iNOS expression, increased proline and polyamine production, and decreased citrulline and NO levels (47, 61, 62). Consistent with these localized pathological mechanisms, a significant reduction in citrulline is observed in peripheral monocytes obtained from patients with lung, and ovarian cancer. Strategies involving the use of ARG1 inhibitors to reinvigorate M1 immunity may, in turn, reactivate tumor-specific Th1 immunity and restore cytotoxic activity (62).
Cancer-induced changes in tryptophan and asparagine in monocytes are of particular interest. Disruptions in tryptophan/asparagine metabolic pathways are associated with tumor processes (21, 22, 26, 63–67). Tryptophan is an essential amino acid, and its availability is a critical factor in the strength and quality of the immune response. The proliferation and activation of human T cells were significantly suppressed in environments lacking tryptophan compared to standard nutrient-rich environments (68, 69). Cancer cells, myeloid-derived suppressor cells (MDSCs), and tumor-associated fibroblasts can reduce tryptophan levels through the enzymatic activity of indoleamine 2,3-dioxygenase (IDO) (21). MDSCs and tumor cells upregulate IDO, inducing an immunosuppressive microenvironment through at least two mechanisms: tryptophan depletion and the accumulation of tryptophan catabolites, such as kynurenine, 3-hydroxyanthranilate, and quinolinate (70, 71). Tryptophan depletion inhibits the proliferation of activated T cells, while catabolites derived from tryptophan act as ligands for aromatic hydrocarbon receptors (72). Kynurenine acts as a suppressor of T-cell immunity, directing the differentiation of naive T cells towards regulatory T cells and suppressing the differentiation of Th17 cells (73).
Intratumoral metabolic status significantly depends on the availability of the nutrients. Tumor reprograms metabolic pathways, promoting metabolic autonomy (74–77). The discovery of oncogenes and tumor suppressors that regulate nutrient uptake and utilization has revealed that nutrients themselves play a key role in cell growth and proliferation (77–80). Significant efforts are being made to develop agents that deplete other amino acids and target central metabolic pathways disrupted in cancer cells, including glycolysis, the tricarboxylic acid cycle, and lipogenesis. Many of these drugs are still in preclinical stages, but some are currently validated in clinical trials (26). The cancer-specific difference in citrulline levels in circulating monocytes are of particular interest since citrulline can link metabolic and epigenetic programs. Citrulline levels can reflect the levels of citrullinated histones. Citrullination of histones is an essential epigenetic mechanism acting both via relaxation of chromatin and as well by interaction with other histone modifications. Enzymes that control citrullination are critical factors of cancer progression, and their expression as well as level of citrullinated proteins in cancer cells and in tumor-associated neutrophils are bad prognostic factors. Inhibitors of these enzymes are under development as cancer therapy drugs (81, 82) However, the role of these enzymes and citrullination of histones in monocytes and TAMs is almost unknown in cancer. Here for the first time we found that levels of citrulline are good prognostic factors for lung cancer. Our data suggest that at least in lung cancer patients citrullination can have a protective role. Thus, the application of inhibitors of the enzymes responsible for the citrullination has to consider their effect also on monocyte-derived TAMs in cancer-specific context.
In conclusion, we found statistically significant cancer-specific changes in the amino acid profile in peripheral blood monocytes in patients who have not received any anti-cancer therapy. Our data suggest that cancer localization has significant impact on the specific changes in the aminoacid profiles in monocytes of patients. The most pronounced changes in amino acid metabolism in monocytes of cancer patients compared to healthy donors were found for breast (decrease in tryptophan and aspartic acid) and ovarian (decrease in citrulline) cancers, and these changes can be indicative for the immunosuppressive programming of monocytes. In contrast, in colorectal cancer monocytes compared to monocytes of patients with other types of cancer, significant increase in aspartic acid and citrulline was identified. This finding provides novel explanation for the anti-tumor activity of TAMs in colorectal cancer, that, in contrast to other types of cancer, retain the ability to inhibit tumor growth. Thus, changes in the aminoacid metabolism in monocytes are involved in the regulation of the balance between pro- and anti-inflammatory activity of macrophages, and can contribute to the pro-tumoral TAM programming during their differentiation from monocytes. Our current study aimed to identify major changes in amino acid metabolism in CD14+ monocytes population that cover over 90% of circulating monocytes, and is composed of both CD14+CD16- and CD14+CD16+ monocytes (83). Due to the limited amount of blood that can be obtained from cancer patients due to ethical restriction, we would be unable to collect sufficient amounts of CD14-CD16+ monocytes to analyse them separately. The other limitations of the study have to be also considered, that include relatively small number of patients for each individual type of cancer, focus of the study on the primary tumors, where majority of patients had no metastasis, absence of analysis of potential gender or age factor contribution, that have to be separately assessed in the large cohort groups in future. However, the statistical reliability of the data suggests the significance of the findings, which can be further validated in larger patient’s cohorts. Correlation of metabolic, transcriptional and epigenetic changes in monocytes will enable identification of the mechanistic links between these major levels of TAMs programming (84). Such new mechanistic links are needed to develop a new generation of TAM-targeted anti-cancer therapy.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics statement
The studies involving humans were approved by Ethic Committee of Tomsk National Research Medical Centre, Tomsk, Russia; Ethic Committee of the Tomsk Regional Oncology Center, Tomsk, Russia. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
VC: Data curation, Investigation, Software, Writing – review & editing. NSt: Conceptualization, Methodology, Writing – original draft, Writing – review & editing. AT: Data curation, Formal analysis, Software, Writing – review & editing. AN: Data curation, Formal analysis, Validation, Writing – review & editing. MP: Resources, Validation, Writing – review & editing. IL: Conceptualization, Validation, Writing – review & editing. EP: Investigation, Resources, Writing – review & editing. MR: Data curation, Visualization, Writing – review & editing. AK: Methodology, Resources, Writing – review & editing. ET: Investigation, Resources, Writing – review & editing. NSh: Conceptualization, Project administration, Writing – review & editing. JK: Conceptualization, Funding acquisition, Writing – review & editing. VF: Conceptualization, Funding acquisition, Writing – review & editing, Writing – original draft. GS: Funding acquisition, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. The study was financially supported by the Russian Federation represented by the Ministry of Science and Higher Education of the Russian Federation (agreement dated 29 September 2021 № 075-15-2021-1073 on the topic “Genetic and epigenetic editing of tumor cells and the microenvironment in order to block metastasis “).
Acknowledgments
Authors wish to thank the Core Facility “Clinical Mass-spectrometry” of the National Medical Research Center for Obstetrics Gynecology and Perinatology Named after Academician V.I. Kulakov of the Ministry of Healthcare of Russian Federation for providing HPLC-MS/MS equipment for sample analysis. Authors are grateful to Bragina O.D., Tarasova A.S., Dobrodeev A.Yu., Kostromitskiy D.N., Afanasiev S.G., Villert A.B., and Kolomiets L.A. from the clinical departments of Cancer Research Institute, Tomsk National Research Medical Center, Russian Academy of Sciences, Tomsk, Russia for providing blood samples from patients with breast, ovarian and colorectal cancers.
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/fimmu.2023.1332043/full#supplementary-material
References
1. Kzhyshkowska J, Gudima A, Moganti K, Gratchev A, Orekhov A. Perspectives for monocyte/macrophage-based diagnostics of chronic inflammation. Transfus Med Hemother (2016) 43:66–77. doi: 10.1159/000444943
2. Nikitina E, Larionova I, Choinzonov E, Kzhyshkowska J. Monocytes and macrophages as viral targets and reservoirs. Int J Mol Sci (2018) 19:2821. doi: 10.3390/ijms19092821
3. Patysheva M, Larionova I, Stakheyeva M, Grigoryeva E, Iamshchikov P, Tarabanovskaya N, et al. Effect of early-stage human breast carcinoma on monocyte programming. Front Oncol (2022) 11:800235. doi: 10.3389/fonc.2021.800235
4. Matuschik L, Riabov V, Schmuttermaier C, Sevastyanova T, Weiss C, Klüter H, et al. Hyperglycemia induces inflammatory response of human macrophages to CD163-mediated scavenging of hemoglobin-haptoglobin complexes. Int J Mol Sci (2022) 23:1385. doi: 10.3390/ijms23031385
5. Patysheva M, Frolova A, Larionova I, Afanas’ev S, Tarasova A, Cherdyntseva N, et al. Monocyte programming by cancer therapy. Front Immunol (2022) 13:994319. doi: 10.3389/fimmu.2022.994319
6. Larionova I, Kazakova E, Gerashchenko T, Kzhyshkowska J. New angiogenic regulators produced by tams: Perspective for targeting tumor angiogenesis. Cancers (Basel) (2021) 13:1–40. doi: 10.3390/cancers13133253
7. Hinshaw DC, Shevde LA. The tumor microenvironment innately modulates cancer progression. Cancer Res (2019) 79:4557–67. doi: 10.1158/0008-5472.CAN-18-3962
8. Larionova I, Tuguzbaeva G, Ponomaryova A, Stakheyeva M, Cherdyntseva N, Pavlov V, et al. Tumor-associated macrophages in human breast, colorectal, lung, ovarian and prostate cancers. Front Oncol (2020) 10:566511. doi: 10.3389/fonc.2020.566511
9. Qian BZ, Li J, Zhang H, Kitamura T, Zhang J, Campion LR, et al. CCL2 recruits inflammatory monocytes to facilitate breast-tumour metastasis. Nature (2011) 475:222–5. doi: 10.1038/nature10138
10. Hanna RN, Cekic C, Sag D, Tacke R, Thomas GD, Nowyhed H, et al. Patrolling monocytes control tumor metastasis to the lung. Science (80-) (2015) 350:985–90. doi: 10.1126/science.aac9407
11. Shi C, Pamer EG. Monocyte recruitment during infection and inflammation. Nat Rev Immunol (2011) 11:762–74. doi: 10.1038/nri3070
12. Kiss M, Caro AA, Raes G, Laoui D. Systemic reprogramming of monocytes in cancer. Front Oncol (2020) 10:1399. doi: 10.3389/fonc.2020.01399
13. Gabrilovich DI, Ostrand-Rosenberg S, Bronte V. Coordinated regulation of myeloid cells by tumours. Nat Rev Immunol (2014) 12:253–68. doi: 10.1038/nri3175.Coordinated
14. Olingy CE, Dinh HQ, Hedrick CC. Monocyte heterogeneity and functions in cancer. J Leukoc Biol (2019) 106:309–22. doi: 10.1002/JLB.4RI0818-311R
15. Engblom C, Pfirschke C, Pittet MJ. The role of myeloid cells in cancer therapies. Nat Rev Cancer (2016) 16:447–62. doi: 10.1038/nrc.2016.54
16. Cassetta L, Fragkogianni S, Sims AH, Swierczak A, Forrester LM, Zhang H, et al. Human tumor-associated macrophage and monocyte transcriptional landscapes reveal cancer-specific reprogramming, biomarkers, and therapeutic targets. Cancer Cell (2019) 35:588–602.e10. doi: 10.1016/j.ccell.2019.02.009
17. Hamm A, Prenen H, Van Delm W, Di Matteo M, Wenes M, Delamarre E, et al. Tumour-educated circulating monocytes are powerful candidate biomarkers for diagnosis and disease follow-up of colorectal cancer. Gut (2016) 65:990–1000. doi: 10.1136/gutjnl-2014-308988
18. Chittezhath M, Dhillon MK, Lim JY, Laoui D, Shalova IN, Teo YL, et al. Molecular profiling reveals a tumor-promoting phenotype of monocytes and macrophages in human cancer progression. Immunity (2014) 41:815–29. doi: 10.1016/j.immuni.2014.09.014
19. Wang L, Simons DL, Lu X, Tu TY, Avalos C, Chang AY, et al. Breast cancer induces systemic immune changes on cytokine signaling in peripheral blood monocytes and lymphocytes. EBioMedicine (2020) 52:102631. doi: 10.1016/j.ebiom.2020.102631
20. Zhang W, Ruan J, Zhou D, Han X, Zhang Y, Wang W, et al. Predicting worse survival for newly diagnosed T cell lymphoma based on the decreased baseline CD16–/CD16 + monocyte ratio. Sci Rep (2020) 10:7757. doi: 10.1038/s41598-020-64579-z
21. Leone RD, Powell JD. Metabolism of immune cells in cancer. Nat Rev Cancer (2020) 20:516–31. doi: 10.1038/s41568-020-0273-y
22. Viola A, Munari F, Sánchez-Rodríguez R, Scolaro T, Castegna A. The metabolic signature of macrophage responses. Front Immunol (2019) 10:1462. doi: 10.3389/fimmu.2019.01462
23. Schairer DO, Chouake JS, Nosanchuk JD, Friedman AJ. The potential of nitric oxide releasing therapies as antimicrobial agents. Virulence (2012) 3:271–9. doi: 10.4161/viru.20328
24. Van den Bossche J, Baardman J, Otto NA, van der Velden S, Neele AE, van den Berg SM, et al. Mitochondrial dysfunction prevents repolarization of inflammatory macrophages. Cell Rep (2016) 17:684–96. doi: 10.1016/j.celrep.2016.09.008
25. Qualls JE, Subramanian C, Rafi W, Smith AM, Defreitas AA, Shirey KA, et al. Sustained generation of nitric oxide and control of mycobacterial infection requires argininosuccinate synthase 1. Cell Host Microbe (2012) 12:313–23. doi: 10.1016/j.chom.2012.07.012.Sustained
26. Schmidt DR, Patel R, Kirsch DG, Lewis CA, Vander Heiden MG, Locasale JW. Metabolomics in cancer research and emerging applications in clinical oncology. CA Cancer J Clin (2021) 71:333–58. doi: 10.3322/caac.21670
27. Odom JD, Sutton VR. Metabolomics in clinical practice: improving diagnosis and informing management. Clin Chem (2021) 67:1606–17. doi: 10.1093/clinchem/hvab184
28. Clish CB. Metabolomics: an emerging but powerful tool for precision medicine. Mol Case Stud (2015) 1:a000588. doi: 10.1101/mcs.a000588
29. Wishart DS. Emerging applications of metabolomics in drug discovery and precision medicine. Nat Rev Drug Discovery (2016) 15:473–84. doi: 10.1038/nrd.2016.32
30. Larionova I, Patysheva M, Iamshchikov P, Kazakova E, Kazakova A, Rakina M, et al. PFKFB3 overexpression in monocytes of patients with colon but not rectal cancer programs pro-tumor macrophages and is indicative for higher risk of tumor relapse. Front Immunol (2023) 13:1080501. doi: 10.3389/fimmu.2022.1080501
31. Chen DP, Ning WR, Jiang ZZ, Peng ZP, Zhu LY, Zhuang SM, et al. Glycolytic activation of peritumoral monocytes fosters immune privilege via the PFKFB3-PD-L1 axis in human hepatocellular carcinoma. J Hepatol (2019) 71:333–43. doi: 10.1016/j.jhep.2019.04.007
32. Frankevich VE, Novoselova AV, Starodubtseva NL, Patysheva MR, Larionova IV, Rakina MA, et al. Methodology of determining the metabolomic profile of tumor-associated macrophages and monocytes in oncological diseases. Bull Russ State Med Univ (2022) 5:58–64. doi: 10.24075/brsmu.2022.049
33. CoreTeam R. R: A language and environment for statistical computing (2021). Available at: https://www.r-project.org.
34. R team. R studio: integrated development for R. Boston, MA: RStudio, Inc (2016). Available at: http://www.rstudio.com/.
35. Cohen J. Statistical power analysis for the behavioral science. 2nd ed. New York: Lawrence Erlbaum Associates All (1988).
36. Bron S, Henry L, Faes-van’t Hull E, Turrini R, Vanhecke D, Guex N, et al. TIE-2-expressing monocytes are lymphangiogenic and associate specifically with lymphatics of human breast cancer. Oncoimmunology (2016) 5:e1073882. doi: 10.1080/2162402X.2015.1073882
37. Zhang B, Cao M, He Y, Liu Y, Zhang G, Yang C, et al. Increased circulating M2-like monocytes in patients with breast cancer. Tumor Biol (2017) 39:1010428317711571. doi: 10.1177/1010428317711571
38. Butler KL, Clancy-Thompson E, Mullins DW. CXCR3 + monocytes/macrophages are required for establishment of pulmonary metastases. Sci Rep (2017) 7:45593. doi: 10.1038/srep45593
39. Larionova I, Cherdyntseva N, Liu T, Patysheva M, Rakina M, Kzhyshkowska J. Interaction of tumor-associated macrophages and cancer chemotherapy. Oncoimmunology (2019) 8:e1596004. doi: 10.1080/2162402X.2019.1596004
40. Manz MG, Boettcher S. Emergency granulopoiesis. Nat Rev Immunol (2014) 14:302–14. doi: 10.1038/nri3660
41. Feng F, Zheng G, Wang Q, Liu S, Liu Z, Xu G, et al. Low lymphocyte count and high monocyte count predicts poor prognosis of gastric cancer. BMC Gastroenterol (2018) 18:148. doi: 10.1186/s12876-018-0877-9
42. Hayashi T, Fujita K, Nojima S, Hayashi Y, Nakano K, Ishizuya Y, et al. Peripheral blood monocyte count reflecting tumor-infiltrating macrophages is a predictive factor of adverse pathology in radical prostatectomy specimens. Prostate (2017) 77:1383–8. doi: 10.1002/pros.23398
43. Sasaki A, Iwashita Y, Shibata K, Matsumoto T, Ohta M, Kitano S. Prognostic value of preoperative peripheral blood monocyte count in patients with hepatocellular carcinoma. Surgery (2006) 139:755–64. doi: 10.1016/j.surg.2005.10.009
44. Lee YY, Choi CH, Sung CO, Do IG, Huh S, Song T, et al. Prognostic value of pre-treatment circulating monocyte count in patients with cervical cancer: Comparison with SCC-Ag level. Gynecol Oncol (2012) 124:92–7. doi: 10.1016/j.ygyno.2011.09.034
45. Shigeta K, Kosaka T, Kitano S, Yasumizu Y, Miyazaki Y, Mizuno R, et al. High absolute monocyte count predicts poor clinical outcome in patients with castration-resistant prostate cancer treated with docetaxel chemotherapy. Ann Surg Oncol (2016) 23:4115–22. doi: 10.1245/s10434-016-5354-5
46. Chen G, Wu K, Li H, Xia D, He T. Role of hypoxia in the tumor microenvironment and targeted therapy. Front Oncol (2022) 12:961637. doi: 10.3389/fonc.2022.961637
47. Capuano G, Rigamonti N, Grioni M, Freschi M, Bellone M. Modulators of arginine metabolism support cancer immunosurveillance. BMC Immunol (2009) 10:1. doi: 10.1186/1471-2172-10-1
48. Somasundaram V, Basudhar D, Bharadwaj G, No JH, Ridnour LA, Cheng RYS, et al. Molecular mechanisms of nitric oxide in cancer progression, signal transduction, and metabolism. Antioxid Redox Signal (2019) 30:1124–43. doi: 10.1089/ars.2018.7527
49. Matos A, Carvalho M, Bicho M. Arginine and arginases modulate metabolism, tumor microenvironment and prostate cancer progression. Nutrients (2021) 13:4503. doi: 10.3390/nu13124503
50. Hendrickson M, Poyton R. Crosstalk between nitric oxide and hypoxia-inducible factor signaling pathways: an update. Res Rep Biochem (2015) 5:147–61. doi: 10.2147/rrbc.s58280
51. Morris SM. Recent advances in arginine metabolism: Roles and regulation of the arginases. Br J Pharmacol (2009) 157:922–30. doi: 10.1111/j.1476-5381.2009.00278.x
52. Mori M, Gotoh T. Arginine metabolic enzymes, nitric oxide and infection. J Nutr (2004) 134:2820–5. doi: 10.1093/jn/134.10.2820s
53. Bron L, Jandus C, Andrejevic-Blant S, Speiser DE, Monnier P, Romero P, et al. Prognostic value of arginase-II expression and regulatory T-cell infiltration in head and neck squamous cell carcinoma. Int J Cancer (2013) 132:E85–93. doi: 10.1002/ijc.27728
54. Mumenthaler SM, Yu H, Tze S, Cederbaum SD, Pegg AE, Seligson DB, et al. Expression of arginase II in prostate cancer. Int J Oncol (2008) 32:357–65. doi: 10.3892/ijo.32.2.357
55. Tate DJ, Vonderhaar DJ, Caldas YA, Metoyer T, Patterson JR IV, Aviles DH, et al. Effect of arginase II on L-arginine depletion and cell growth in murine cell lines of renal cell carcinoma. J Hematol Oncol (2008) 1:14. doi: 10.1186/1756-8722-1-14
56. Mao H, Gao W, Lu G, Fang F, Teng L. Clinicopathological and prognostic implications of arginase expression in hepatocellular carcinoma. Clin Lab (2013) 59:37–43. doi: 10.7754/Clin.Lab.2012.120210
57. de Boniface J, Mao Y, Schmidt-Mende J, Kiessling R, Poschke I. Expression patterns of the immunomodulatory enzyme arginase 1 in blood, lymph nodes and tumor tissue of early-stage breast cancer patients. Oncoimmunology (2012) 1:1305–12. doi: 10.4161/onci.21678
58. Satriano J. Arginine pathways and the inflammatory response: Interregulation of nitric oxide and polyamines: Review article. Amino Acids (2004) 26:321–9. doi: 10.1007/s00726-004-0078-4
59. Grillo MA, Colombatto S. Arginine revisited: Minireview article. Amino Acids (2004) 26:345–51. doi: 10.1007/s00726-004-0081-9
60. Munder M, Eichmann K, Morán JM, Centeno F, Soler G, Modolell M. Th1/Th2-regulated expression of arginase isoforms in murine macrophages and dendritic cells. J Immunol (1999) 163:3771–7. doi: 10.4049/jimmunol.163.7.3771
61. Massi D, Marconi C, Franchi A, Bianchini F, Paglierani M, Ketabchi S, et al. Arginine metabolism in tumor-associated macrophages in cutaneous Malignant melanoma: evidence from human and experimental tumors. Hum Pathol (2007) 38:1516–25. doi: 10.1016/j.humpath.2007.02.018
62. Sica A, Allavena P, Mantovani A. Cancer related inflammation: The macrophage connection. Cancer Lett (2008) 267:204–15. doi: 10.1016/j.canlet.2008.03.028
63. Andrejeva G, Rathmell JC. Similarities and distinctions of cancer and immune metabolism in inflammation and tumors. Cell Metab (2017) 26:49–70. doi: 10.1016/j.cmet.2017.06.004
64. Pan S, Fan M, Liu Z, Li X, Wang H. Serine, glycine and one-carbon metabolism in cancer (Review). Int J Oncol (2021) 58:158–70. doi: 10.3892/ijo.2020.5158
65. Zhu Y, Li X, Wang L, Hong X, Yang J. Metabolic reprogramming and crosstalk of cancer-related fibroblasts and immune cells in the tumor microenvironment. Front Endocrinol (Lausanne) (2022) 13:988295. doi: 10.3389/fendo.2022.988295
66. Neuman RE, McCoy TA. Dual requirement of walker carcinosarcoma 256 in vitro for asparagine and glutamine. Science (80-) (1956) 124:124–5. doi: 10.1126/science.124.3212.124
67. Clavell LA, Gelber RD, Cohen HJ, Suzanne Hitchcock-Bryan RN, Cassady JR, Tarbel NJ, et al. Four-agent induction and intensive asparaginase therapy for treatment of childhood acute lymphoblastic leukemia. N Engl J Med (1986) 315:657–63. doi: 10.1056/NEJM198609113151101
68. Munn DH, Sharma MD, Baban B, Harding HP, Zhang Y, Ron D, et al. GCN2 kinase in T cells mediates proliferative arrest and anergy induction in response to indoleamine 2,3-dioxygenase. Immunity (2005) 22:633–42. doi: 10.1016/j.immuni.2005.03.013
69. Munn BDH, Shafizadeh E, Attwood JT, Bondarev I, Pashine A, Mellor AL. Inhibition of T cell proliferation by macrophage. J Exp Med (1999) 189:1363–72. doi: 10.1084/jem.189.9.1363
70. Platten M, Wick W, Van Den Eynde BJ. Tryptophan catabolism in cancer: Beyond IDO and tryptophan depletion. Cancer Res (2012) 72:5435–40. doi: 10.1158/0008-5472.CAN-12-0569
71. Moffett JR, Namboodiri MA. Tryptophan and the immune response. Immunol Cell Biol (2003) 81:247–65. doi: 10.1046/j.1440-1711.2003.t01-1-01177.x
72. Opitz CA, Litzenburger UM, Sahm F, Ott M, Tritschler I, Trump S, et al. An endogenous tumour-promoting ligand of the human aryl hydrocarbon receptor. Nature (2011) 478:197–203. doi: 10.1038/nature10491
73. Stephens GL, Wang Q, Swerdlow B, Bhat G, Kolbeck R, Fung M. Kynurenine 3-monooxygenase mediates inhibition of Th17 differentiation via catabolism of endogenous aryl hydrocarbon receptor ligands. Eur J Immunol (2013) 43:1727–34. doi: 10.1002/eji.201242779
74. DeBerardinid RJ, Sayed N, Ditsworth D, Thompson CB. Brick by brick: metabolisn and tumor growth. Curr Opin Genet Dev (2008) 18:54–61. doi: 10.1016/j.gde.2008.02.003
75. Galluzzi L, Kepp O, Heiden MGV, Kroemer G. Metabolic targets for cancer therapy. Nat Rev Drug Discovery (2013) 12:829–46. doi: 10.1038/nrd4145
76. Vander Heiden MG. Exploiting tumor metabolism: Challenges for clinical translation. J Clin Invest (2013) 123:3648–51. doi: 10.1172/JCI72391
77. Ward PS, Thompson CB. Metabolic reprogramming: A cancer hallmark even warburg did not anticipate. Cancer Cell (2012) 21:297–308. doi: 10.1016/j.ccr.2012.02.014
78. Csibi A, Lee G, Yoon SO, Tong H, Ilter D, Elia I, et al. The mTORC1/S6K1 pathway regulates glutamine metabolism through the eif4b-dependent control of c-Myc translation. Curr Biol (2014) 24:2274–80. doi: 10.1016/j.cub.2014.08.007
79. Wise DR, Deberardinis RJ, Mancuso A, Sayed N, Zhang XY, Pfeiffer HK, et al. Myc regulates a transcriptional program that stimulates mitochondrial glutaminolysis and leads to glutamine addiction. Proc Natl Acad Sci U.S.A. (2008) 105:18782–7. doi: 10.1073/pnas.0810199105
80. Yuneva M, Zamboni N, Oefner P, Sachidanandam R, Lazebnik Y. Deficiency in glutamine but not glucose induces MYC-dependent apoptosis in human cells. J Cell Biol (2007) 178:93–105. doi: 10.1083/jcb.200703099
81. Teng Y, Chen Y, Tang X, Wang S, Yin K. PAD2: A potential target for tumor therapy. Biochim Biophys Acta - Rev Cancer (2023) 1878:188931. doi: 10.1016/j.bbcan.2023.188931
82. Zhu D, Zhang Y, Wang S. Histone citrullination: a new target for tumors. Mol Cancer (2021) 20:1–17. doi: 10.1186/s12943-021-01373-z
83. Ziegler-Heitbrock L, Hofer TPJ. Toward a refined definition of monocyte subsets. Front Immunol (2013) 4:23. doi: 10.3389/fimmu.2013.00023
Keywords: mass spectrometry, metabolomics, oncology, monocytes, tumor-associated macrophages
Citation: Chagovets V, Starodubtseva N, Tokareva A, Novoselova A, Patysheva M, Larionova I, Prostakishina E, Rakina M, Kazakova A, Topolnitskiy E, Shefer N, Kzhyshkowska J, Frankevich V and Sukhikh G (2024) Specific changes in amino acid profiles in monocytes of patients with breast, lung, colorectal and ovarian cancers. Front. Immunol. 14:1332043. doi: 10.3389/fimmu.2023.1332043
Received: 02 November 2023; Accepted: 19 December 2023;
Published: 08 January 2024.
Edited by:
Ramalingam Bethunaickan, National Institute of Research in Tuberculosis (ICMR), IndiaReviewed by:
Nagaraj Manickam, Madras Diabetes Research Foundation, IndiaParamasivam Natarajan, AbbVie, United States
Balaji Pathakumari, Mayo Clinic, United States
Copyright © 2024 Chagovets, Starodubtseva, Tokareva, Novoselova, Patysheva, Larionova, Prostakishina, Rakina, Kazakova, Topolnitskiy, Shefer, Kzhyshkowska, Frankevich and Sukhikh. 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: Julia Kzhyshkowska, SnVsaWEua3poeXNoa293c2thQG1lZG1hLnVuaS1oZWlkZWxiZXJnLmRl; Vladimir Frankevich, dmZyYW5rZXZpY2hAZ21haWwuY29t
†These authors have contributed equally to this work
‡These authors share senior authorship