SUMMARY
The application of metabolomics for studying modifications in host metabolism due to viral infections has proven to be a game-changing approach. Prior to our study, only one other ‘omics’ study has been carried out that investigates the interplay between the host and CCHFV and its subsequent pathogenesis. We employed NMR spectroscopy, given its advantages in terms of reproducibility, minimal sample preparation, and capability to analyze complex biofluids. Our methodology builds upon the proven success of metabolomics in the research of other viral hemorrhagic fevers such as Ebola, Marburg, and Dengue. Our research underlines the critical role of SAH, a metabolite involved in numerous biochemical reactions. We provide new insights into the metabolic alterations occurring in CCHF patients. These alterations not only shed light on the disease’s pathogenesis but also pave the way for potential biomarkers and therapeutic targets. Among all the metabolites detected, S-Adenosyl-L-homocysteine and Carnosine stood out as the most prevalent, warranting further exploration of their roles in CCHFV pathogenesis and their potential as therapeutic targets.
1. Introduction
Crimean-Congo Hemorrhagic Fever (CCHF) is a severe viral infection with substantial implications for public health in a wide range of countries in Asia, Africa, Southern Europe and the Middle East. CCHF is caused by CCHF virus (CCHFV) belonging to the genus of Orthonairovirus, and the family of Nairoviridae being one of deadliest viruses of its kind with reported mortality rate of 3-30% (1). CCHFV is transmitted to humans through tick bites of infected ticks, contact with blood or tissues of infected livestock, or contact with infected patients (2). The exact course of pathogenesis of CCHF is not clearly known however it is divided into four phases: incubation, pre-haemorrhagic, haemorrhagic and convalescence (2, 3). CCHF is classified as a severe hemorrhagic fever with a short incubation period of 1-3 days although longer incubation periods have been documented (3, 4). The onset of infection is often sudden and includes symptoms, fever, diarrhea, vomiting, nausea, myalgia, back and abdominal pain followed by an hemorrhagic phase where severe bruises, uncontrollable bleeding at the body orifices are observed and in severe cases, deterioration of kidneys, liver and lungs (5). Deaths associated with the infection mostly occur between 5-14 days from the start of the viremic phase (6, 7). In terms of treatment, early hospitalizations and early administration of therapeutics are shown to reduce both severity and mortality of CCHF (8). Hence, lack of early detection of CCHF is one of the leading factors causing the particular high mortality rate of CCHF.
There are several challenges of diagnosing CCHF infection particularly before the hemorrhagic phase of the infection and patients who are not suspected of being bitten by infected ticks or contacted with infected livestock. CCHF is relatively rare in certain areas, so healthcare providers may not consider CCHF as a possible diagnosis at first. Likewise, there are documented cases of difficulties diagnosing CCHF, both because of the latter but also the absence of a universally applicable diagnostic kit for surveillance and diagnosis of all CCHFV strains (9). Standard blood tests such as hemogram, biochemical analysis and physical examination at the beginning of hospitalization is applied on all patients who are suspected to be infected with CCHF although the results are often relevant for short term prognostic factors as biochemical values change often quickly. More comprehensive tests such as viral antigen and nucleic acid amplification tests are used as standard diagnostic methods (10). The initial symptoms of CCHF, such as fever, headache, and muscle aches, are rather non-specific hence diagnosing CCHF early by differential diagnosis can be difficult. (3,6) Biomarkers, essentially biomolecules, provide a measure of specific diseases or their stages due to their varying concentrations. They are instrumental in diagnosing and monitoring the progression of viral infections. The associated changes in their levels, often indicative of the disease, are typically attributed to the host’s immune reaction and the disturbance of key biochemical routes in reaction to the infectious process. In this sense, omics studies and biomarkers could be used for both analysis and diagnosis of CCHF and for characterizing better treatment strategies of hospitalized patients swiftly is crucial to pinpoint optimal treatment strategies for treating CCHF (11).
So far there is only one omics study which investigates host-viral response and pathogenesis of CCHF utilizing transcriptomics and proteomics methods (12). Led by this gap in the literature, we conducted a nationwide analysis of metabolomes of patients hospitalized due to CCHF. Turkey has over 10,000 cases of CCHF with an average fatality rate of 5%, making it a critical public health concern affecting people living in rural areas as ticks are widespread in these regions (13). In our present investigation, we employed Nuclear Magnetic Resonance (NMR) spectroscopy exploratory metabolomics to investigate the overall temporal variations in plasma metabolites during a seasonal outbreak of CCHF infection. We employed PLS-DA statistical analysis of the blood serum metabolome of CCHF patients and categorized certain metabolomes linked to metabolic dysregulation caused by CCHF. Preliminary results from our study suggest that specific metabolic markers can be identified in the serum of CCHF patients pointing to metabolic dysregulation, which may allow for earlier diagnosis and more targeted treatment strategies. Additionally, being the first study to categorize alterations of patient metabolome during CCHF viremic phase may be valuable for efforts to develop therapeutics or targeted treatment strategies to reduce the severity and high mortality rate of CCHF.
2. Materials and methods
Monitoring of CCHF diagnosed patients and collection of the samples from the patients
Patients who were diagnosed with CCHF were selected from Çorum Hitit University Hospital. The first sample collection was conducted in April 2022, and blood samples were collected daily using Ethylenediaminetetraacetic acid (EDTA) tubes. In current metabolomics studies, opinions regarding the appropriate determination of sample size can considerably vary. However, numerous statistical analyses have underscored that a substantial sample size, for achieving meaningful results, around 30 samples (14). In this study, serum specimens were procured from 29 patients (n=29) diagnosed with CCHFV infection, as well as from 10 healthy control group. CCHF patients also admitted to the hospital with an infection diagnosis were subclassified into two categories based on their blood test values and symptoms: moderate (n=24) and severe (n=5). Four blood serum samples from patients each consecutive day (n=116 samples) and a single sample from the control group (n=10 samples) was taken. These blood samples were then subjected to centrifugation for 5 minutes at 3000g to separate the sample into plasma, white blood cell, and red blood cell phases. The plasma phase was extracted and subjected to metabolite extraction using multiple approaches: single methanol extraction, triple alcohol extraction, Methanol-chloroform, Acetone, Acetonitrile, and Ultrafiltration. Cold methanol-chloroform was chosen as the most effective extraction and used for extraction of all samples. Only the polar metabolites in the plasma were investigated, while proteins and apolar compounds were removed from the plasma samples. Overall in this study, we utilized a rigorous approach to sample collection, preparation, and analysis to investigate polar metabolites in the plasma of CCHF patients (15, 16, 17).
Ethical Statement
The study was approved by Koç University Ethics Committee and all the procedures performed in this study involving human participants were in accordance with the ethical standards of the institutional research committee ethical standards. A written informed consent was obtained from all patients.
NMR sample preparation
4 mL ice-cold methanol-chloroform (1:1) mixture was added to 2 mL serum for methanol-chloroform extraction. The mixture was vortexed for 30 seconds and incubated for 10 minutes on ice. After incubation the mixture was centrifuged at 4500g at 4°C for 30 minutes. The methanol phase was collected and dried using a vacuum concentrator. The dried samples were dissolved in 550 uL D2O based NMR sample solution (50 mM PBS (pH 7.4), 20mM NaCl, 1 mM DSS) for standardized sample preparation.
NMR data collection, processing and statistical analysis
500 MHz Bruker Ascend magnet with BBO paired resonance probe and Avance NEO console was used for NMR data collection. 1D NOESY-presat (noesygppr1d) pulse sequence was used for data collection. Each NMR data spectrum is composed of 4K screening and 32K complex data points. Spectrum widths were set to 9615.4 Hz. Bruker Topspin 4.2.0 software was used for NMR data processing. Data was divided into 0.02 ppm data packages along with their normalization coefficients. The dataset, which comprises data packets with a resolution of 0.02 ppm, was analyzed using the MetaboAnalyst 5.0 online metabolomics statistical analysis software. Henceforth the data will be referred to as (Bin.x.xx [ppm]) data packets and the day of when sample was taken. All data points were normalized using the average centering normalization method. Following this normalization, the dataset underwent statistical analysis using Partial Least Squares Discriminant Analysis (PLS-DA). PLS-DA is a classification and discrimination technique based on the Partial Least Squares (PLS) regression method. This method is widely utilized to determine the differences between classes, particularly in high-dimensional and multivariate datasets. PLS-DA is a commonly employed method in analyzing complex biological systems, such as metabolomics studies. The VIP Projection variable importance score plot is a graph that is utilized to assess the outcomes of PLS-DA and determine the most significant variables in the analysis. VIP scores quantify the importance of each variable (e.g., metabolites) in classification and aid in identifying the most critical features. VIP scores are computed based on the contribution of each variable to the components (latent variables) in the PLS-DA model. The values begin at 1, and higher VIP scores indicate that the variable is more important for classification. Variables with VIP scores greater than 1 are generally deemed significant, although this threshold may vary in practice. The VIP score plot displays the VIP scores of the variables on the vertical axis, while the variables themselves or their indices are shown on the horizontal axis. This graph facilitates the identification of important variables visually and helps focus on the variables that require prioritization in the analysis. In the VIP score graph, the peaks at the relevant ppm values that make up the data packages (Bin.x.xx) have been examined in more detail and the metabolites to which they belong have been identified. For this operation, NMR spectra have been reopened, and the metabolites to which the peak in the relevant ppm region belongs have been determined using the Chenomx software. It is thought that some peaks might belong to metabolites not found in the database, and the molecules these peaks belong to have not been identified. Further investigation and characterization may be required to fully understand these unidentified peaks and their role in the overall metabolic profile, ensuring that the final analysis provides an accurate reflection of the biological system under investigation.
Metabolomic pathway visualization
All metabolic pathways visualized using Metastate software Version BETA (https://metastate.bio) (Metastate Bio Inc.). Metastate algorithm employs the Kyoto Encyclopedia of Genes and Genomes (KEGG) database as its foundational input source. Software systematically retrieves details pertaining to biological pathways, chemical compounds, and molecular reactions of interest. Software curates and assembles a dynamic graphical representation of the data.
3. Results
The demographic characteristics of the participants in the study, includes 29 patients diagnosed with CCHF and 10 healthy individuals serving as controls. The mean age of the control group was 50.1 years, with an age range between 40 to 64 years, while the CCHF patient group had a slightly higher mean age of 50.5 years, though with a wider age range of 22 to 77 years. In total 126 blood samples were collected, their serum extracted by methanol extraction, and NMR samples were prepared for all samples. Following metabolite extraction, NMR analyses were conducted on the samples. Finally, a statistical analysis approach was applied. In the statistical analyses, all samples collected over a period of 4 days. All groups were subjected to statistical evaluation as moderate, severe, and control groups as described in materials and methods.
In PLS-DA analysis, a model is developed to include the primary components (latent variables) of the data. The first five components are considered important in classifying the data and explain the majority of the dataset. The matching score plots exhibit these first five components in two-dimensional graphs and compare each component with another (Figure 1a). Subsequently, the next step involved obtaining the Partial Least Squares Discriminant Analysis (PLS-DA) score plot, as shown in Figure 1b demonstrates a high discrimination between healthy and diseased individuals. PLS-DA effectively deals with highly dimensional data, highlighting the relevant variables responsible for the distinction, thereby providing a clear differentiation between health states based on the variation and correlation within the biomarker data. In order to examine the significant data packets responsible for this discrimination, the PLS-DA VIP Projection variable importance score plot is illustrated in Figure 1c. This technique, graphically demonstrated in Figure 1c, reveals the variables that have substantial influence on the discrimination. VIP scores signify the importance of each variable in the model, enabling us to identify and focus on the key biomarkers that are significantly altering between healthy and diseased states. This plot aids in the in-depth interpretation of the complex multivariate data and fosters a more nuanced understanding of the disease markers. Figure 1c showcases the top 15 variables that have the greatest impact on discrimination among these variables. As shown in Figure 1c, it is noticeable that several data packets exhibit low concentration values (blue) in the control group and gradually increase towards higher concentrations with the passage of time. Ultimately as depicted in Figure 1d, the listed compounds exhibited a notable increase in their concentration, as substantiated by increased signals from their corresponding data buckets. Using a statistical analysis approach consisting of three different groups. Samples taken on four consecutive days from the two patient groups were successfully differentiated from the control group. However, no discernible distinction was observed between the moderate and severe groups. Overall these compounds were detected to be elevated compared to the control: SAH, GTP, Carnosine, maleate, 2-Deoxyuridine, IMP, AMP and NADP+.
As a secondary approach to data analysis, data collected on the day 1 and day 2 were used as severe patient and moderate patient groups to further investigate the metabolite profile of the patients during the initial stages of hospitalization. Subsequently, only samples obtained on day 1 and day 2 were categorized as diseased samples, and the blood samples obtained from healthy individuals were designated as the control group for statistical analysis (Fig 2). In addition to compounds detected in our initial analysis novel compounds with significant increase during day 1 and 2 detected that may have a more active role in the patogenesis.
Finally, in a comprehensive metabonomic analysis of CCHF patients’ blood serum, distinct patterns of metabolites were observed between the severe and moderate infection levels. Specifically, through the PLS-DA score plots and Variable Importance Score (VIP) projections detailed in Figure 3, notable variations in the metabolite profiles were discerned. In the case of severe CCHF infections, the compounds AMP, IMP, and NAAD were identified as present within the serum, whereas they were not as significantly detected in the samples from patients with moderate infection levels. Conversely, GTP was found to be significantly increased in the samples of patients with moderate CCHF infection levels, but not in those categorized as severe. Overall the results outlined in our final analysis reveal a marked metabolic variety between the moderate and severe CCHF infection levels which may indicate the critical pathways in CCHF patogenesis.
4. Discussion
In this study we demonstrated our results of the first non-targeted CCHF metabolomics study to improve the current understanding of the pathogenesis of CCHF for diagnostic purposes and aid studies on development of potential therapies against CCHF. Our main findings revealed considerable increases in the levels of SAH, GTP, Carnosine, Maleate, 2-Deoxyuridine, IMP, AMP and NADP+ in the blood serum of CCHF patients, suggesting these metabolites could play crucial roles in the pathogenesis of CCHF and therefore, potentially serve as important biomarkers for early detection and monitoring of the disease progression.
Metabolomics has emerged as a pivotal approach for inspecting the modifications in host metabolism elicited by viral infections. Up until the present, only a single omics study has been conducted by MS that explores the interaction between the host and CCHFV and the ensuing pathogenesis (12). NMR and MS are the two commonly utilized methods used for metabolomics studies with each having certain strengths and weaknesses depending on their purpose of use. NMR spectroscopy offers better reproducibility, minimal sample preparation, and the ability to analyze complex biofluids, whereas MS has relatively higher sensitivity and better for targeted metabolomics studies. Although reproducibility and superior comparability of NMR spectra across different instruments and laboratories compared to MS makes it preferable for non targeted metabolome analysis for diagnostic and prognostic analysis of several diseases (18). NMR spectroscopy in metabolomics is used in analysis of infection on host metabolism and its implications for the pathogenesis of viruses such as HIV (20), Dengue virus (DENV) (21), Chikungunya virus (22) Metabolomics analysis has been employed in the study of other viral hemorrhagic fevers, such as Ebola, Marburg, and Dengue, with promising results for early diagnosis and prognosis prediction (23). Our study demonstrates a new path by harnessing the power of NMR metabolomics analysis to further decipher the intricate interplay between the host and CCHFV and to elucidate the mechanisms of CCHFV pathogenesis.
S-adenosyl homocysteine (SAH) is an intermediate metabolite, a precursor to homocysteine and adenosine(24). Moreover, SAH is the substrate of SAH hydrolase enzyme which is a crucial enzyme involved in the S-adenosylmethionine (SAM/AdoMet) regeneration cycle (Figure 4a) (25). This cycle is noteworthy in the context of CCHF patients, as evidenced by increased SAH levels from the beginning of hospitalization. Furthermore, a genomic study investigating the effect of Methylenetetrahydrofolate reductase (MTHFR) polymorphism showed that MTHFR polymorphism creates a predisposition to milder CCHF (26). MTHFR’s primary role as an enzyme in folate metabolism extends to functions related to methylation as well. Thus, we hypothesized that alterations in methylation pathways during the viremic phase of CCHF may serve as valid indicators of prognosis in CCHF patients. Furthermore, SAM is understood to function as a primary methyl donor in various cellular methylation reactions, including those necessary for the 5’ RNA capping process (27). Similar to the methylation process in Flaviviruses (28, 29), the Ebola virus’ L protein has been found to exhibit methyltransferase activity, affecting co-transcriptional modifications at the cap structure and internal adenosine-2’-O-methylation (30). This additional complexity in RNA methylation within some viral families could provide a framework to consider the implications of SAM-related domains in replication complexes as potential drug targets.
GTP, AMP and inosine monophosphate (IMP) are nucleotides involved in purine metabolism, which play a vital role in energy metabolism (Figure 4b), and nucleic acid synthesis (31). Increase in the level of GTP and IMP are known indicators for infection of several viruses belonging to Orthornavirae kingdom (32). In fact, increase in purine compounds can be attributed to infection of CCHFV considering known effective broad spectrum antivirals such as Ribavirin as a nucleotide analog targets the substrate-binding site of the IMPDH enzyme and prevents binding of IMP causing reduction of RNA synthesis in infected cells (33, 34) through down regulating GTP synthesis. Monitoring the levels of these metabolites may provide vital insights into the efficacy of antiviral treatments, especially nucleotide analogs like Ribavirin. Identifying shifts in GTP and IMP concentrations could be a good addition to therapeutic strategies for clinical trials, optimizing the use of antiviral agents that target purine metabolism, and thereby enhancing the response to CCHF infection. (35, 36) Carnosine (beta-alanyl-L-histidine) (Figure 4c) is a dipeptide with antioxidant, free radical scavenging; anti-glycation, and antiinflammatory properties which is predominantly found in skeletal muscle and in the brain (29). Increase in carnosine levels may reflect changes in the antioxidant defense system, reduction of muscle mass due to hospitalization and infection related cell death and as a possible anti-inflammatory response to the CCHFV infection. Carnosine is considered as a potential therapy for Zika and Dengue viruses, (37) and SARS-CoV-2 (38) as it showed significant reduction in viral replication and ease the symptoms of these diseases, however the effectiveness of carnosine as a therapeutic against CCHF requires further investigation. 2’-Deoxyuridine (2’-dU) is an intermediate in the synthesis of thymidylate (Figure 4d), which acts as a precursor for DNA synthesis and Edoxudine as an antiviral therapeutic. (39). Increase in the levels of 2-deoxyuridine can have significant implications for viral pathogenesis as 2’-dU variants such as BVDU is effective against Herpes simplex virus type 1 (HSV-1) and varicella-zoster virus (VZV) (40). Nevertheless 2’-dU is an antiviral of DNA viruses and used as an indicator to determine presence of viruses which have DNA as their genetic material (41). Presumably 2’-dU levels increased in the blood serum may indicate a dysregulation in the pyrimidine and energy metabolism in accordance to our other findings which signal metabolic disruption of these pathways. Nicotinamide Adenine Dinucleotide Phosphate (NADP+), another key metabolite, serves as a crucial cofactor in various enzymatic reactions, primarily within the pentose phosphate pathway, where it facilitates the synthesis of NADPH, essential for fatty acid synthesis and the regeneration of reduced glutathione. Additionally, NADP+ plays a vital role in supporting the body’s antioxidant defense system and is involved in the reduction-oxidation (redox) reactions. It is difficult to definitively pinpoint why NADP+ increases in CCHF patients, as the underlying cause could be multifaceted, encompassing a range of biological, metabolic, and external factors. Furthermore in relation to our other findings it is possible that elevation of NADP+ also overlaps with increase of GTP, IMP, and AMP which take part in energy metabolism and nucleotide synthesis, as a possible indication of a coordinated increase in use of energy and nucleic acid synthesis.. This is further supported by our final analysis between severe and moderate infection CCHF in Fig 3, considering AMP, IMP and NAAD levels increased more significantly in severe cases signaling increased activity in metabolic pathways including these compounds. Maleate (cis-butenedioic acid), a dicarboxylic acid, and trans-isomer of fumaric acid is a metabolite taking role in nicotinate and nicotinamide metabolism. (42) Previous studies on rats have demonstrated that the injection of sodium maleate can induce a generalized renal transport defect resembling Fanconi syndrome. Furthermore, maleate has been shown to reduce cellular levels of coenzyme A, inhibit the tricarboxylic acid cycle, lower ATP concentrations, and affect various enzymatic activities. Its impact even extends to the ultrastructure of kidney cells. The elevated presence of maleate in the blood serum of CCHF patients could signal underlying metabolic and cellular disruptions, warranting further investigation into its potential role and implications in the disease. (43,44)
Conclusion
In this study we conducted a nationwide analysis of blood serum metabolites of patients admitted to hospital due to CCHF. We employed NMR spectroscopy as a novel approach to characterize host viral response and pathogenesis of CCHFV. Among detected compounds, most prevalent of them were detected to be SAH and carnosine and notable increases in compounds relating to TCA cycle, nucleic acid synthesis and redox related coenzymes in the cell.
Data Availability
All data produced in the present study are available upon reasonable request to the authors
Author Contributions
Çağdaş Dağ: Conceptualization, Methodology, Supervision, Funding acquisition, Writing-Reviewing and Editing, Writing-Original draft preparation Kerem Kahraman: Investigation, Formal analysis, Visualization, Writing-Reviewing and Editing Oktay Göcenler: Investigation, Formal analysis, Writing-Original Draft, Visualization Derya Yapar: Investigation Yaren Kahraman: Investigation, Formal analysis Cengizhan Büyükdağ: Formal analysis, Writing-Original Draft Gülen Esken: Methodology, Investigation Serena Ozabrahamyan: Methodology, Investigation Tayfun Barlas: Investigation Yüksel Karadağ: Investigation Aysel Kocagül Çelikbaş: Investigation Füsun Can: Supervision, Investigation Nurcan Baykam: Supervision, Investigation Mert Kuşkucu: Supervision, Writing-Reviewing and Editing Önder Ergönül: Supervision, Funding acquisition, Conceptualization, Methodology
ACKNOWLEDGMENT
CD acknowledges support from TÜBİTAK (Project No: 221S353). The authors acknowledge the use of the services and facilities of n2STAR-Koç University Nanofabrication and Nanocharacterization Center for Scientific and Technological Advanced Research. The authors gratefully acknowledge use of the services and facilities of the Koç University Is Bank Infectious Disease Center (KUIS-CID). We are also immensely grateful to Fırat Kahya, Boran Saruhan and Oğuzcan Ünver for constructive feedback and for their comments.