ORIGINAL RESEARCH article

Front. Cell. Infect. Microbiol., 24 July 2024

Sec. Molecular Bacterial Pathogenesis

Volume 14 - 2024 | https://doi.org/10.3389/fcimb.2024.1413103

The predictive value of peripheral blood cell mitochondrial gene expression in identifying the prognosis in pediatric sepsis at preschool age

  • 1. Key Laboratory of Birth Defects and Related Diseases of Women and Children of MOE, Department of Pediatrics, West China Second University Hospital, Sichuan University, Chengdu, Sichuan, China

  • 2. Department of Emergency, Chengdu Women’s and Children’s Central Hospital, School of Medicine, University of Electronic Science and Technology of China, Chengdu, Sichuan, China

Abstract

Background:

Sepsis represents a severe manifestation of infection often accompanied by metabolic disorders and mitochondrial dysfunction. Notably, mitochondrial DNA copy number (mtDNA-CN) and the expression of specific mitochondrial genes have emerged as sensitive indicators of mitochondrial function. To investigate the utility of mitochondrial gene expression in peripheral blood cells for distinguishing severe infections and predicting associated outcomes, we conducted a prospective cohort study.

Methods:

We established a prospective cohort comprising 74 patients with non-sepsis pneumonia and 67 cases of sepsis induced by respiratory infections, aging from 2 to 6 years old. We documented corresponding clinical data and laboratory information and collected blood samples upon initial hospital admission. Peripheral blood cells were promptly isolated, and both total DNA and RNA were extracted. We utilized absolute quantification PCR to assess mtDNA-CN, as well as the expression levels of mt-CO1, mt-ND1, and mt-ATP6. Subsequently, we extended these comparisons to include survivors and non-survivors among patients with sepsis using univariate and multivariate analyses. Receiver operating characteristic (ROC) curves were constructed to assess the diagnostic potential.

Results:

The mtDNA-CN in peripheral blood cells was significantly lower in the sepsis group. Univariate analysis revealed a significant reduction in the expression of mt-CO1, mt-ND1, and mt-ATP6 in patients with sepsis. However, multivariate analysis did not support the use of mitochondrial function in peripheral blood cells for sepsis diagnosis. In the comparison between pediatric sepsis survivors and non-survivors, univariate analysis indicated a substantial reduction in the expression of mt-CO1, mt-ND1, and mt-ATP6 among non-survivors. Notably, total bilirubin (TB), mt-CO1, mt-ND1, and mt-ATP6 levels were identified as independent risk factors for sepsis-induced mortality. ROC curves were then established for these independent risk factors, revealing areas under the curve (AUCs) of 0.753 for TB (95% CI 0.596–0.910), 0.870 for mt-CO1 (95% CI 0.775–0.965), 0.987 for mt-ND1 (95% CI 0.964–1.000), and 0.877 for mt-ATP6 (95% CI 0.793–0.962).

Conclusion:

MtDNA-CN and mitochondrial gene expression are closely linked to the severity and clinical outcomes of infectious diseases. Severe infections lead to impaired mitochondrial function in peripheral blood cells. Notably, when compared to other laboratory parameters, the expression levels of mt-CO1, mt-ND1, and mt-ATP6 demonstrate promising potential for assessing the prognosis of pediatric sepsis.

Introduction

Currently, infectious diseases remain a predominant threat to children’s health. Alarming statistics reveal that over 13 million children worldwide succumb to infectious diseases annually, constituting 63% of child fatalities (). Among these, sepsis, a rapidly advancing condition triggered by severe infection, not only stands as the leading cause of hospitalization and mortality in pediatric intensive care units (PICUs) but also represents a significant source of childhood disability. Consequently, a compelling clinical imperative exists to actively identify biomarkers characterized by high specificity and sensitivity for predicting the prognosis of patients with pediatric sepsis. Research has unearthed several candidates, such as serum-soluble urokinase-type plasminogen activator receptor (SuPAR), neutrophil CD64, pro-adrenomedullin, and cell-free plasma DNA (cfDNA), which exhibit the potential to forecast outcomes sepsis patient. However, out of more than 5,000 studies scrutinizing sepsis biomarkers from 2009 to 2019, merely 189 biomarkers linked to sepsis prognosis have been identified. Among these, only pro-adrenomedullin demonstrated a notably high predictive value, with an area under the curve (AUC) exceeding 0.8 (). Consequently, the quest for effective biomarkers to evaluate the prognosis of patients with sepsis remains ongoing. Therefore, it is critical to demonstrate a more efficient sepsis biomarker and assess its prognosis.

Mitochondria are not only the primary cellular energy source but also pivotal contributors to the immune response. They actively engage in the antiviral immune response through Toll-like and RIG-1-like receptor signaling pathways. They participate in the immune response against bacterial infections by regulating the production of reactive oxygen species (ROS) and orchestrating metabolic adaptations in phagocytic cells. Moreover, they play a critical role in facilitating the metabolic changes essential for differentiating T and B lymphocytes, significantly contributing to the adaptive immune response (). It is imperative to note that mitochondrial dysfunction can impair immune cell functionality, and this immune dysfunction is intricately linked to the onset and progression of sepsis. The mitochondrial respiratory chain is the primary executor of mitochondrial function, comprising complexes I to IV. Its principal role involves the conveyance of hydrogen and electrons, with the energy generated in this process harnessed by complex V for the synthesis of ATP, which serves as the body’s energy currency. Within this intricate system, the mitochondrial genes mt-ND1, mt-CO1, and mt-ATP6 function as subunits of complex I, complex IV, and complex V, respectively. Both mitochondrial DNA copy number (mtDNA-CN) and mitochondrial gene expression serve as valuable indicators of mitochondrial function. Moreover, in the previous studies of our group, mtDNA-CN and specific mitochondrial genes expressions were associated with mitochondrial function, and several studies applied such parameters into mitochondrial functional assessment (; ). According to previous studies, mtDNA-CN was associated with several types of diseases’ prognosis via the critical regulation mechanisms of mitophagy in mitochondrial quality control (). However, the specific genes’ expressions of mtDNA directly reflected the function of respiratory chain. Thus, assessing the exact expression of particular genes, such as mt-ND1, mt-CO1, and mt-ATP6, might serve as a more precious role in determining related clinical outcomes of infection. Moreover, in clinical practice, it was difficult to apply rapid ATP production measurement, as the mitochondria might be damaged during ex vivo transportation and storage. Thus, the examination of mtDNA-CN and specific genes expression would be much easier in clinical application, while numerous studies have explored the correlation between mtDNA-CN, mitochondrial gene expression, and infection.

Herein, we carried out this research to illustrate the changes of mtDNA-CN and mitochondrial gene expression (mt-ND1, mt-CO1, and mt-ATP6) in peripheral blood cells. Then, we attempted to identify the potential role of mtDNA-CN and mitochondrial gene expression in measuring the severity of infectious diseases and in predicting the prognosis of some patients with pediatric sepsis. The research revealed the advantages of mitochondrial DNA transcriptional levels in septic management compared with general laboratory parameters.

Materials and methods

Ethics statement

This study was conducted following the ethical guidelines and was approved by the Ethics Committee of West China Second University Hospital of Sichuan University (approval number: 2021–069). Informed consent was obtained from all parents or legal guardians of the patients, who provided their consent for including their child’s clinical and imaging details in the manuscript for publication purposes. This is a prospectively designed research. All participants in this study were hospitalized in either the respiratory department or the general intensive care unit of West China Second University Hospital of Sichuan University between January 2022 and December 2022. There were two groups involved in the observational cohort: the non-sepsis respiratory infection (non-S) group and the respiratory infection-induced sepsis (S) group.

Inclusion and exclusion criteria

The inclusion criteria included the following (1): the diagnosis of sepsis meeting the criteria from the 2016 International Guidelines for Sepsis and Septic Shock, which should be induced by or have originated from initial respiratory infection (sepsis 3.0 diagnosis criteria); (2) the non-sepsis respiratory infection was recognized as bronchial pneumonia without other complications based on the Guidelines for the Management of Community-Acquired Pneumonia in Children in China; (3) the age of enrolled patients ranged from 2 to 6 years; and (4) the peripheral blood samples should be collected at the time of initial hospital admission before antibiotic agent supplementation.

The exclusion criteria were as follows: (1) patients with viral infection-associated sepsis; (2) patients with birth defects; (3) bronchopulmonary dysplasia history; (4) chest trauma history; (5) patients with immunodeficiency; (6) diagnosed with inherited metabolic disorders; (7) genetic test confirmed mitochondrial diseases; (8) patients with neoplastic diseases; (9) recent surgery history; (10) beyond respiratory or combined other system infections; (11) incomplete medical archive or laboratory tests; (12) lost to follow-up; and (13) severe chronic diseases.

The observational end point of the cohort had been designed as short-term follow-up of 1 month after hospital discharge or patient death. The clinical survival was considered as primary outcome.

Blood sample collection and nucleotide extraction

Peripheral blood samples (4 mL) were collected from the patients on admission, of which 2 mL was stored in EDTA tubes and 2 mL was stored in RNA tubes; then, such samples were placed in a −40°C refrigerator. The experimental samples were kept and registered by special personnel, and the extraction time of RNA and DNA shall not exceed 6 months. Blood DNA Kit (D3392–02) and Blood RNA Kit (R6814–02, Omega, USA) were used to extract DNA and RNA from peripheral blood samples, and the concentration and purity of DNA and RNA were determined and stored at −20°C.

Plasmid construction and mtDNA-CN measurement

To obtain the DNA sequence amplification of mt-ND1, design primers (Forward—ATACCCATGGCCAACCTCCTA; Reverse—TAGGTTTGAGGGGGAATGCTG) were used to amplify the whole mitochondrial mt-ND1 gene to quantify the mitochondrial genome. Fragments were amplified using NEB Q5 High-Fidelity DNA polymerase. Then, the mt-ND1 sequence was connected to the T vector according to the Mighty TA-cloning Reagent Set for Prime STAR® (Code N0.6019 TaKaRa Company) kit instructions for research and the plasmid was linearized to make a gradient concentration standard. The number of DNA copies was calculated according to the following formula:

DNA (copies/μL) = (DNA concentration (ng/μL) × 6.022 × 1023)/(Length (bp) × 1 × 109 × 660) [6.022 × 1023 = Avogadro’s constant, 1 × 109 = mole amount converted to ng, 660 = average mass of a pair of double-stranded DNA bases (g), Length = vector + length of insert]. Dilute the linearized plasmid 10-fold with ddH2O to a concentration of 103–108 copies/μL to make a standard.

RT-qPCR quantification of MT-CO1, MT-ND1, and MT-ATP6

RNA was reverse-transcribed into complementary deoxyribonucleic acid [complementary DNA (cDNA)] (according to HiScriptIII AII-in-one RT SuperMix Perfect for qPCR), and then RT-PCR reaction was performed (Qiagen SYBR Green Mixture reagent protocol). The relative gene expression was calculated using the 2−ΔΔCt method. The sequence of each gene was amplified, as shown in Table 1.

Table 1

GenePrimer sequence (5′-3′)
MT-CO1 ForwardTCTCAGGCTACACCCTAGACCA
MT-CO1 ReverseATCGGGGTAGTCCGAGTAACGT
MT-ND1 ForwardGGCTATATACAACTACGCAAAGGC
MT-ND1 ReverseGGTAGATGTGGCGGGTTTTAGG
ATP6 ForwardGAAGCGCCACCCTAGCAATA
ATP6 ReverseGCTTGGATTAAGGCGACAGC
18S ForwardTTGACGGAAGGGCACCACCAG
18S ReverseGCACCACCACCCACGGAATCG

Gene amplification primer sequence.

Statistical analysis

Statistical analysis was performed by GraphPad Prism software (Version 8, San Diego, CA). Normality test would be completed initially, and the data were presented as mean ± SD once they passed normality test. Otherwise, the continuous data would be presented in median with range, while the categorical variables were expressed as n (%). For continuous variable data passed normality test, an independent-sample t-test was used between two groups; for categorical variables, a chi-square test was used for comparison. When p < 0.05, it was judged as a significant statistical difference. Univariate and multivariate analyses had been applied to analyze the clinical and laboratory indicators and the expression of mt-CO1, mt-ND1, and mt-ATP6 between survivors and non-survivors in sepsis, and identify particular independent risk factors. Then, SPSS 22.0 software had been used to illustrate the ROC curve of independent risk factors to predict adverse clinical outcome of sepsis, and AUC had been used to clarify the test efficiency.

Results

Characters of included patients

Ninety-three patients were initially enrolled in the non-S group. Seven patients were subsequently excluded due to lost follow-ups, while three patients were identified as having congenital disorders, seven had concomitant injuries to other systems, and two patients did not complete coagulation tests. Thus, the final analysis included 74 patients in the non-S group. In the S group, 97 patients were initially included. Among them, 8 patients had congenital birth defects, 5 patients had positive genetic analysis results, 3 were diagnosed with immunodeficiency diseases, 2 were previously diagnosed with bronchopulmonary dysplasia, and 12 patients had concurrent COVID-19 either before or after the sepsis diagnosis. Consequently, 67 cases with complete medical records were included for further analysis. Baseline information was compared between the two groups, revealing no significant differences in gender distribution, average age, or body weight. There were 51 patients (76%) in the S group identified with positive results by blood culture or metagenomics analysis, including 17 cases with S. aureus, 15 cases with S. pneumoniae, 7 cases with P. aeruginosa, 5 cases with H. influenzae, 3 cases with K. pneumoniae, 3 cases with E. coli, and 1 case with E. faecium. However, the S group had a substantially longer average hospital stay duration (23.68 ± 7.74 days) than the non-S group (4.74 ± 1.37 days). Additionally, the S group experienced 11 in-hospital non-survivors, whereas no deaths were recorded in the non-S group (Table 2).

Table 2

VariableNon-S (n = 74)S (n = 67)p
Female/male (n)27/4736/31
Age (m)45.5 ± 14.2340.25 ± 17.60
Weight (kg)16.24 ± 3.3115.39 ± 4.48
Hospital duration (days)4.74 ± 1.3723.68 ± 7.74<0.001
Death (n)0 (0%)5 (7.46%)0.022
Blood culture or metagenomics analysisNot involved51 (76%)
S. aureus (17/25.4%)
S. pneumoniae (15/22.4%)
P. aeruginosa (7/10.4%)
H. influenzae (5/7.5%)
K. pneumoniae (3/4.5%)
E. coli (3/5.4%)
E. faecium (1/1.5%)
WBC (109/L)15.77 ± 22.4611.54 ± 7.810.145
N (%)62.33 ± 17.7958.06 ± 24.310.232
L (%)27.34 ± 17.0629.61 ± 19.470.461
HGB (g/L)119.05 ± 9.04101.13 ± 27.60<0.001
CRP (mg/L)35.21 ± 30.9972.45 ± 69.06<0.001
PCT (μg/L)1.18 ± 1.0412.59 ± 27.31<0.001
ALT (U/L)21.97 ± 9.38192.72 ± 436.93<0.001
AST (U/L)42.76 ± 10.47301.22 ± 754.10<0.001
TB (mmol/L)5.24 ± 2.5116.95 ± 25.92<0.001
ALB (g/L)42.83 ± 2.8139.20 ± 22.340.167
GLB (g/L)22.48 ± 3.4419.65 ± 7.340.003
γ-GT (U/L)9.98 ± 3.6262.07 ± 96.96<0.001
LDH (U/L)342.37 ± 90.971943.11 ± 9117.900.133
UN (mmol/L)7.04 ± 26.396.77 ± 9.040.938
Cr (μmol/L)21.94 ± 4.8364.07 ± 112.430.001
PT (s)13.28 ± 0.6516.12 ± 16.950.150
APTT (s)30.28 ± 2.9435.21 ± 11.79<0.001
Fg (g/L)4.05 ± 0.54304.72 ± 171.36<0.001
INR1.44 ± 1.111.32 ± 0.350.385
DDI (μg/mL)1.04 ± 0.257.39 ± 10.44<0.001
FDP (μg/mL)3.05 ± 0.5931.78 ± 80.250.002
2-ΔΔCT-CO160.89 ± 128.211.41 ± 0.96<0.001
2-ΔΔCT-ND1101.90 ± 210.261.58 ± 0.98<0.001
2-ΔΔCT-ATP664.90 ± 132.771.33 ± 0.90<0.001

Univariate analysis of laboratory tests and mitochondrial function between non-S and S groups.

Non-S, non-sepsis respiratory infection; S, respiratory infection-induced sepsis. Numerical variables were expressed as mean ± SD; categorical variables were presented as percentage.

Severe infection induced reduction of mtDNA-CN quantification

Both the plasmid standard product with a concentration of 103–108 copies/μL and the DNA samples from the two groups were simultaneously subjected to fluorescent qPCR. Following the completion of multiple batches, a standard curve was constructed. The quantitative templates with various concentrations of standards exhibited a linear correlation with the Ct value. The regression equations for the mt-ND1 standard curve were derived as follows: y = −4.146x + 39.575, where y represents Ct, x signifies the base 10 logarithm of the plasmid template’s copy number, and the correlation coefficients reached 0.9939. Subsequently, the Ct value was substituted into the regression equation of the standard curve to calculate the absolute quantity of mtDNA-CN. Based on the Ct values obtained from qPCR assessment, the average mtDNA-CN was 105.95 ± 0.32 in the S group and 106.11 ± 0.15 in the non-S group. These findings highlight a significant decrease in mtDNA-CN among patients with sepsis. Consequently, it becomes crucial to assess the expression of specific mitochondrial genes transcription from mtDNA and analyze their predictive value in early identification of adverse prognosis associated with sepsis.

Mitochondrial DNA expression level failed to serve as an independent factor for severe infection

The clinical characteristics of patients in the non-S group and S group are presented in Table 2. Univariate analysis had been used to analyze the clinical and laboratory parameters between the two groups. It was found that the values of HBG, CRP, PCT, ALT, AST, TB, GLB, γ-GT, Cr, APTT, Fg, DDI, FDP, 2-ΔΔCT-CO1, 2-ΔΔCT-ND1, and 2-ΔΔCT-ATP6 were significantly different between two groups (p < 0.05), as shown in Table 2. Then, logistic regression had been used for multivariate analysis to evaluate whether such parameters would contribute as independent factors in distinguishing severe conditions among patients with respiratory infection. Among the included variables that were identified based on univariate analysis, HGB (OR = 0.78, p = 0.031), TB (OR = 1.447, p = 0.003), γ-GT (OR = 2.531, p = 0.009), Cr (OR = 1.392, p = 0.034), and FDP (OR = 2.647, p = 0.002) demonstrated their potential in serving as independent risk factors for distinguishing severe respiratory infection, although the expression of mt-CO1, mt-ND1, and mt-ATP6 had been recorded to be extremely reduced in the S group. However, the logistic regression analysis of the values among 2-ΔΔCT-CO1, 2-ΔΔCT-ND1, and 2-ΔΔCT-ATP6 presented no significant meaning in separating sepsis and non-severe infection (Table 3). Therefore, although the expression of mitochondrial DNA coding gene is significantly downregulated in patients with sepsis, the expression of mitochondrial DNA coding genes failed to serve as an independent factor to distinguish sepsis.

Table 3

VariableBS.E.WaldORp
HGB−0.2480.1154.6340.7800.031
CRP0.0320.0183.0241.0330.082
PCT1.5050.8003.5374.5040.060
ALT−0.0170.0550.0900.9830.764
AST0.0740.0790.8761.0770.349
TB0.3700.1268.5631.4470.003
GLB−0.3170.1942.6740.7280.102
GT0.9290.3556.8412.5310.009
Cr0.3310.1564.4901.3920.034
APTT0.0320.0410.5991.0320.439
Fg0.44726.7870.0001.5640.987
DDI−0.0960.2020.2230.9090.637
FDP0.9740.3199.3402.6470.002
2-ΔΔCT-CO10.4620.7960.3371.5870.561
2-ΔΔCT-ND1−0.6510.6970.8730.5210.350
2-ΔΔCT-ATP60.2060.3370.3731.2290.541

Multivariate analysis to identify independent risk factors between Non-S and S groups.

Downregulated mitochondrial gene transcription was associated with septic lethality

In the next step analysis, we also performed univariate analysis to demonstrate potential factors in predicting adverse clinical outcomes among patients with sepsis. According to the initial results, the values of TB, UN, PT, APTT, FDP, 2-ΔΔCT-CO1, 2-ΔΔCT-ND1, and 2-ΔΔCT-ATP6 presented significant differences between the septic survivors and non-survivors (p < 0.05, Table 4). Then, logistic regression has been used to evaluate the potential role of the initially identified parameters that serve as an independent factor to predict septic lethality (Table 5). According to the results, the expressions of TB (OR = 1.042, p = 0.037), mt-ND1 (OR = 0.00, p = 0.021), mt-CO1 (OR = 0,004, p = 0.004), and mt-ATP6 (OR = 0.007, p = 0.008) revealed persistent strength in assessing clinical outcomes of sepsis, which could serve as independent factors in early identifying non-survivors.

Table 4

VariableSurvivors (n = 56)Non-survivors (n = 11)p
WBC (109/L)11.74 ± 13.7810.53 ± 6.200.642
N (%)56.62 ± 44.6965.84 ± 24.580.243
L (%)30.27 ± 32.3926.65 ± 24.640.575
N/L (%)4.04 ± 12.256.03 ± 6.050.352
HGB (g/L)99.98 ± 59.18107 ± 25.890.445
CRP (mg/L)80.33 ± 131.6535.18 ± 59.080.055
PCT (μg/L)12.05 ± 65.0714.17 ± 33.890.827
ALT (U/L)172.12 ± 948.46295.72 ± 327.500.395
AST (U/L)286.81 ± 1638.05373.27 ± 422.360.731
TB (mmol/L)12.31 ± 22.5722.4 ± 15.310.025
ALB (g/L)35.09 ± 11.4332.85 ± 6.850.268
GLB (g/L)19.95 ± 17.0318.13 ± 6.270.456
γ-GT (U/L)63.32 ± 259.5055.72 ± 25.010.814
LDH (U/L)2040.5 ± 27463.451447.36 ± 1931.940.845
UN (mmol/L)5.34 ± 18.2511.31 ± 7.750.013
Cr (μmol/L)51.35 ± 229.23101.54 ± 104.980.106
PT (s)13.77 ± 6.3529.08 ± 40.480.005
APTT (s)33.97 ± 17.6342.02 ± 18.570.038
Fg (g/L)310.67 ± 367.22272.3 ± 165.720.501
INR1.30 ± 0.671.46 ± 0.420.191
DDI (μg/mL)6.81 ± 15.9110.59 ± 14.690.276
FDP (μg/mL)21.23 ± 47.3390.08 ± 182.020.008
2-ΔΔCT-CO11.56 ± 1.300.69 ± 0.240.005
2-ΔΔCT-ND11.79 ± 1.600.52 ± 0.130.000
2-ΔΔCT-ATP61.48 ± 1.550.54 ± 0.180.001

Univariate analysis between survivors and non-survivors in patients with sepsis.

Numerical variables were expressed as mean ± SD; categorical variables were presented as percentage.

Table 5

VariableBS.E.WaldORp
TB0.0410.0242.9341.0420.037
UN0.0730.0472.3931.0750.122
PT0.0700.1270.3001.0720.584
APTT0.0160.0350.2021.0160.653
FDP0.0040.0130.1021.0040.750
MT-CO1−5.6281.9488.3500.0040.004
MT-ND1−29.31212.6595.3620.0000.021
MT-ATP6−4.8961.8417.0760.0070.008

Multivariate analysis reveals the independent risk factors in assessing non-survivors in sepsis.

In order to further clarify the predictive value of mitochondrial DNA expression assessment in the clinical outcome of sepsis, ROC curves had been illustrated to determine the predictive value of different indicators. The AUC of the ROC among the four parameters had been calculated as 0.753 for TB (95% CI 0.596–0.910), 0.870 for mt-CO1 (95% CI 0.775–0.965), 0.987 for mt-ND1 (95% CI 0.964–1.000), and 0.877 for mt-ATP6 (95% CI 0.793–0.962) (Figure 1; Table 6).

Figure 1

Table 6

VariableAUCS.E.95% CI
TB0.7530.0800.596–0.910
MT-CO10.8700.0480.775–0.965
MT-ND10.9870.0120.964–1.000
MT-ATP60.8770.0430.793–0.962

AUC of particular ROC curves in determining survivor in patients with sepsis.

AUC, area under the curve; CI, confidence interval; S.E., standard error; TB, total bilirubin.

Discussion

The pathogenesis of sepsis is inherently intricate. Traditionally, it has been attributed to the body’s response to infection becoming imbalanced. Upon entering the body, pathogens engage with pathogen-related molecular patterns, triggering immune cell activation and inflammatory responses. These responses, in turn, set off complex cascade reactions within and between cells, intensifying the inflammatory processes. This ultimately results in heightened capillary permeability and vasodilation, leading to decreased blood pressure, tissue underperfusion, and abnormal blood coagulation. Notably, research has demonstrated a reciprocal promotion between inflammatory responses and coagulation abnormalities in sepsis. Further investigations have revealed that the initial stages of sepsis are characterized by inflammation and immune effector cell activation aimed at pathogen clearance (; ). As sepsis progresses, an anti-inflammatory response gains prominence, accompanied by diminished counts and compromised function of immune cells, including monocytes, lymphocytes, and dendritic cells. This phase also witnesses reduced expression and production of pro-inflammatory cytokines, coupled with increased anti-inflammatory factors, all working together to mitigate inflammation and facilitate tissue repair. Evidently, immune dysfunction plays a pervasive role throughout the onset and progression of sepsis, and the intertwined dynamics of inflammatory response and coagulation dysfunction closely contribute to its development.

Mitochondria possess their own genome, known as mt-DNA, which stands as the sole DNA molecule residing in the human cytoplasm independently from the nuclear chromosomes. The quantification of these genomes, termed mtDNA-CN, serves as a reflection of mitochondrial size, quantity, and functional status (). Mitochondria, dynamic organelles crucial for cellular homeostasis, serve as the primary site for ATP production. They also engage in fatty acid β-oxidation and the tricarboxylic acid cycle, regulate cellular calcium metabolism and signal transduction, and play a role in programmed cell death. Moreover, mitochondria actively contribute to the innate immune response during pathogen infections through various pathways. Concurrently, mitochondria are essential for the metabolic transformations required for the differentiation of T and B lymphocytes and play a pivotal role in the adaptive immune response (). Research indicates that during sepsis, immune cells experience mitochondrial dysfunction, characterized by decreased mitochondrial membrane potential in peripheral blood mononuclear cells (PBMCs) and energy depletion (; ). Additionally, the activity of mitochondrial complexes is variably diminished (; ; ; ). Dysfunctional mitochondria profoundly impact immune cell activity, ultimately leading to dysregulated immune responses.

Altered mitochondrial function has been linked to a spectrum of diseases, encompassing neoplastic diseases (), Parkinson’s disease (), cardiovascular diseases (; ), hepatitis C (), acquired immunodeficiency syndrome (), and novel coronavirus pneumonia (). In the context of sepsis, Pyle et al. discovered that the blood cells of surviving patients with sepsis exhibited significantly higher mt-DNA levels compared to those in the non-surviving group. The increased mt-DNA levels were closely associated with the prognosis of patients with sepsis (). Another researcher reported the predictive value of plasma mt-DNA for the early diagnosis and prognosis in patients with sepsis (; ). Busani et al. conducted a study on the expression of mt-DNA genes in peripheral blood and its association with clinical outcomes in patients with sepsis in the PICU. Their findings indicated that mt-DNA can serve as a biomarker for disease severity and even mortality in both adults and children receiving critical care (). MT-CO1 is a functional subunit of cytochrome C oxidase (COX). Mutations in the mt-CO1 coding region lead to structural and functional changes in COX, resulting in increased production of ROS, apoptosis, and genetic susceptibility to sepsis. Moreover, the mt-DNA T6459C mutation is responsible for mitochondrial damage in the early stages of sepsis, leading to significant pathological alterations associated with sepsis (). Rahmel et al. observed reduced abundance of mitochondrial TFAM in PBMCs of patients with sepsis. This reduction was accompanied by diminished mt-DNA copy number, decreased expression of MT-ND1, and a decline in cellular ATP content ().

In our research, we detected the mtDNA-CN of peripheral blood cells between non-S and S groups. The results showed that the mt-DNA copy number of children in the severe infection group was significantly lower. It suggested that the mitochondrial function of blood cells had been impaired in severe infection. However, the expression of mt-CO1, mt-ND1, and mt-ATP6 failed to demonstrate any statistical significance and could not work as an independent factor for distinguishing severe infection. Furthermore, we identified that the expression of mt-ND1, mt-CO1, and mt-ATP6 was associated with sepsis-related death, which had been confirmed by AUC under ROC calculations. According to the results, the expression of mitochondrial genes presented a much higher predictive value than previous applied clinical parameters. Therefore, it was believed that mitochondrial function revealed an efficient predictive value in determining the adverse clinical outcome of children with sepsis.

To date, numerous biomarkers for predicting sepsis prognosis have been extensively investigated. Studies have revealed that elevated levels of CRP and PCT are indicative of poor outcomes in patients with sepsis. Specifically, CRP has demonstrated a favorable discriminative capacity for predicting mortality in children with sepsis beyond 6 months in the PICU. However, the calculated AUC of CRP was not satisfied enough, and other newly defined parameters were warranted (; ). Early studies also highlighted the prognostic potential of cfDNA in predicting sepsis mortality. Dwivedi et al. demonstrate that the AUC for predicting sepsis mortality of cfDNA is notably high at 0.97 (95% CI 0.93–1.00), presenting a perfect predictive value for septic lethality (). Subsequently, Forsblom et al. identified cf-DNA as an independent prognostic marker for death in patients with sepsis (). However, more recent literature reports with a higher sample size have suggested that cfDNA, with an average AUC around 0.76, exhibited less favorable specificity and sensitivity in predicting death among patients with general sepsis (). Moreover, other studies demonstrate that neutrophil CD64 expression shows a relatively weak association with sepsis mortality (AUC < 0.8) (; ), while the AUCs of CYSTM1, MMP8, and CD177, which are strongly associated with a variety of immune cells, were 0.988, 0.973, and 0.986, respectively, in predicting pediatric sepsis (). Notably, several studies have emphasized the sensitivity and specificity of plasma-soluble urokinase-type plasminogen activator receptor (suPAR) as an independent prognostic biomarker in patients with sepsis (; ; ). While it surpasses PCT and CRP in predicting patient mortality, its AUC remains relatively low, with correspondingly low specificity and sensitivity (; ; ). In a recent meta-analysis focused on sepsis mortality, proadrenomedullin was the sole biomarker to achieve an AUC exceeding 0.8, accompanied by high specificity (92% specificity; 75% sensitivity) (). Moreover, in recent years, more and more studies have focused on the predictive value of different biomarker combinations for pediatric sepsis; the AUC of these combinations can even reach 1 (; ; ). Remarkably, in our study, mitochondrial function-related indicators, including mt-CO1, mt-ND1, and mt-ATP6, exhibited significantly higher discriminative capabilities for predicting the clinical outcomes of sepsis-associated lethality compared to traditional laboratory indicators. Therefore, we contend that mitochondrial function indices hold substantial predictive value for the clinical outcomes of sepsis in children.

This study has several limitations that warrant consideration: The sample size in this study is relatively small, and it would benefit from validation through the expansion of the sample cohort. Meanwhile, the patients we included ranged in age from 2 to 6 years, which should be expanded to 0–18 years in future studies. However, the results from this prospective cohort study are deemed credible and hold significant scientific value due to the detailed statistical analyses conducted. The study assessed mitochondrial function by measuring mtDNA-CN and mitochondrial gene expression. It lacks direct methods to evaluate mitochondrial function, which should be addressed in subsequent investigations.

Conclusion

The morbidity and mortality rates associated with pediatric sepsis continue to pose a significant challenge. Early prognosis prediction serves as a valuable tool for clinicians, aiding in the assessment of potential disease progression and informing therapeutic decisions. Both mtDNA-CN and the expression of mitochondrial genes have shown associations with the severity and clinical outcomes of infectious diseases. Moreover, severe infections have been observed to impair the mitochondrial function of peripheral blood cells. In comparison to other laboratory parameters, the expression levels of mt-CO1, mt-ND1, and mt-ATP6 demonstrate notable potential in predicting the prognosis of pediatric sepsis.

Statements

Data availability statement

The original contributions presented in the study are included in the article/supplementary material. Further inquiries can be directed to the corresponding authors.

Ethics statement

This study was conducted following the ethical guidelines and was approved by the Ethics Committee of West China Second University Hospital of Sichuan University (approval number: 2021–069). Informed consent was obtained from all parents or legal guardians of the patients, who provided their consent for including their child’s clinical and imaging details in the manuscript for publication purposes. This is a prospectively designed research. All participants in this study were hospitalized in either the respiratory department or the general intensive care unit of West China Second University Hospital of Sichuan University between January 2022 and December 2022. There were two groups involved in the observational cohort: the non-sepsis respiratory infection (non-S) group and the respiratory infection-induced sepsis (S) group.

Author contributions

SJ: Conceptualization, Data curation, Investigation, Methodology, Resources, Validation, Writing – original draft. YZ: Data curation, Formal analysis, Investigation, Methodology, Writing – original draft. WZ: Investigation, Methodology, Resources, Software, Validation, Writing – original draft. YL: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Software, Supervision, Validation, Writing – original draft, Writing – review & editing. YW: Conceptualization, Data curation, Investigation, Methodology, Project administration, Resources, Validation, Writing – review & editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by grants from the Technology Project of Sichuan Province of China (2021YFQ0061) and the National Natural Science Foundation of China (82270249). The funding sources had no role in the design of the study and collection, analysis, and interpretation of data or in the writing of the manuscript.

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.

References

  • 1

    AsharF. N.ZhangY.LongchampsR. J.LaneJ.MoesA.GroveM. L.et al. (2017). Association of mitochondrial DNA copy number with cardiovascular disease. JAMA Cardiol.2, 1247–1255. doi: 10.1001/jamacardio.2017.3683

  • 2

    BrealeyD.BrandM.HargreavesI.HealesS.LandJ.SmolenskiR.et al. (2002). Association between mitochondrial dysfunction and severity and outcome of septic shock. Lancet.360, 219–223. doi: 10.1016/S0140-6736(02)09459-X

  • 3

    BusaniS.De BiasiS.NasiM.PaoliniA.VenturelliS.TosiM.et al. (2020). Increased plasma levels of mitochondrial DNA and normal inflammasome gene expression in monocytes characterize patients with septic shock due to multidrug resistant bacteria. Front. Immunol.11. doi: 10.3389/fimmu.2020.00768

  • 4

    ChengS. C.SciclunaB. P.ArtsR. J.GresnigtM. S.LachmandasE.Giamarellos-BourboulisE. J.et al. (2016). Broad defects in the energy metabolism of leukocytes underlie immunoparalysis in sepsis. Nat. Immunol.17, 406–413. doi: 10.1038/ni.3398

  • 5

    CuiN.ZhangH.ChenZ.YuZ. (2019). Prognostic significance of PCT and CRP evaluation for adult ICU patients with sepsis and septic shock: retrospective analysis of 59 cases. J. Int. Med. Res.47, 1573–1579. doi: 10.1177/0300060518822404

  • 6

    DimoulaA.PradierO.KassengeraZ.DalcomuneD.TurkanH.VincentJ. L. (2014). Serial determinations of neutrophil CD64 expression for the diagnosis and monitoring of sepsis in critically ill patients. Clin. Infect. Dis.58, 820–829. doi: 10.1093/cid/cit936

  • 7

    DonadelloK.ScollettaS.CovajesC.VincentJ. L. (2012). SuPAR as a prognostic biomarker in sepsis. BMC Med.10, 2. doi: 10.1186/1741-7015-10-2

  • 8

    DonadelloK.ScollettaS.TacconeF. S.CovajesC.SantonocitoC.CortesD. O.et al. (2014). Soluble urokinase-type plasminogen activator receptor as a prognostic biomarker in critically ill patients. J. Crit. Care29, 144–149. doi: 10.1016/j.jcrc.2013.08.005

  • 9

    DwivediD. J.ToltlL. J.SwystunL. L.PogueJ.LiawK. L.WeitzJ. I.et al. (2012). Canadian Critical Care Translational Biology G. Prognostic utility and characterization of cell-free DNA in patients with severe sepsis. Crit. Care16, R151. doi: 10.1186/cc11466

  • 10

    ForsblomE.AittoniemiJ.RuotsalainenE.HelmijokiV.HuttunenR.JylhavaJ.et al. (2014). High cell-free DNA predicts fatal outcome among Staphylococcus aureus bacteraemia patients with intensive care unit treatment. PLoS One9, e87741. doi: 10.1371/journal.pone.0087741

  • 11

    GarrabouG.MorenC.LopezS.TobiasE.CardellachF.MiroO.et al. (2012). The effects of sepsis on mitochondria. J. Infect. Dis.205, 392–400. doi: 10.1093/infdis/jir764

  • 12

    HuangQ.XiongH.YanP.ShuaiT.LiuJ.ZhuL.et al. (2020). The diagnostic and prognostic value of suPAR in patients with sepsis: A systematic review and meta-analysis. Shock.53, 416–425. doi: 10.1097/SHK.0000000000001434

  • 13

    HuttunenR.SyrjanenJ.VuentoR.HurmeM.HuhtalaH.LaineJ.et al. (2011). Plasma level of soluble urokinase-type plasminogen activator receptor as a predictor of disease severity and case fatality in patients with bacteraemia: a prospective cohort study. J. Intern. Med.270, 32–40. doi: 10.1111/j.1365-2796.2011.02363.x

  • 14

    JapiassuA. M.SantiagoA. P.d'AvilaJ. C.Garcia-SouzaL. F.GalinaA.Castro Faria-NetoH. C.et al. (2011). Bioenergetic failure of human peripheral blood monocytes in patients with septic shock is mediated by reduced F1Fo adenosine-5'-triphosphate synthase activity. Crit. Care Med.39, 1056–1063. doi: 10.1097/CCM.0b013e31820eda5c

  • 15

    KochA.VoigtS.KruschinskiC.SansonE.DuckersH.HornA.et al. (2011). Circulating soluble urokinase plasminogen activator receptor is stably elevated during the first week of treatment in the intensive care unit and predicts mortality in critically ill patients. Crit. Care15, R63. doi: 10.1186/cc10037

  • 16

    KochR. E.JosefsonC. C.HillG. E. (2017). Mitochondrial function, ornamentation, and immunocompetence. Biol. Rev. Camb Philos. Soc92, 1459–1474. doi: 10.1111/brv.12291

  • 17

    LarsenF. F.PetersenJ. A. (2017). Novel biomarkers for sepsis: A narrative review. Eur. J. Intern. Med.45, 46–50. doi: 10.1016/j.ejim.2017.09.030

  • 18

    LeonardS.GuertinH.OdoardiN.MillerM. R.PatelM. A.DaleyM.et al. (2024). Pediatric sepsis inflammatory blood biomarkers that correlate with clinical variables and severity of illness scores. J. Inflamm. (London England).21, 7. doi: 10.1186/s12950-024-00379-w

  • 19

    LiX.WeiY.XuZ.LiT.DongG.LiuX.et al. (2023). Lymphocyte-to-C-reactive protein ratio as an early sepsis biomarker for neonates with suspected sepsis. Mediators Inflamm.2023, 9077787. doi: 10.1155/2023/9077787

  • 20

    LiuH.ZhenC.XieJ.LuoZ.ZengL.ZhaoG.et al. (2024). TFAM is an autophagy receptor that limits inflammation by binding to cytoplasmic mitochondrial DNA. Nat. Cell Biol.26, 878–891. doi: 10.1038/s41556-024-01419-6

  • 21

    LivaditiO.KotanidouA.PsarraA.DimopoulouI.SotiropoulouC.AugustatouK.et al. (2006). Neutrophil CD64 expression and serum IL-8: sensitive early markers of severity and outcome in sepsis. Cytokine.36, 283–290. doi: 10.1016/j.cyto.2007.02.007

  • 22

    MalikA. N.CzajkaA. (2013). Is mitochondrial DNA content a potential biomarker of mitochondrial dysfunction? Mitochondrion.13, 481–492. doi: 10.1016/j.mito.2012.10.011

  • 23

    PierrakosC.VelissarisD.BisdorffM.MarshallJ. C.VincentJ. L. (2020). Biomarkers of sepsis: time for a reappraisal. Crit. Care24, 287. doi: 10.1186/s13054-020-02993-5

  • 24

    PlunkettA.TongJ. (2015). Sepsis in children. BMJ.350, h3017. doi: 10.1136/bmj.h3017

  • 25

    PyleA.AnugrhaH.Kurzawa-AkanbiM.YarnallA.BurnD.HudsonG. (2016). Reduced mitochondrial DNA copy number is a biomarker of Parkinson's disease. Neurobiol. Aging.38, 216.e7–216.e10. doi: 10.1016/j.neurobiolaging.2015.10.033

  • 26

    PyleA.BurnD. J.GordonC.SwanC.ChinneryP. F.BaudouinS. V. (2010). Fall in circulating mononuclear cell mitochondrial DNA content in human sepsis. Intensive Care Med.36, 956–962. doi: 10.1007/s00134-010-1823-7

  • 27

    QuickA. P.WangQ.PhilippenL. E.Barreto-TorresG.ChiangD. Y.BeaversD.et al. (2017). SPEG (Striated muscle preferentially expressed protein kinase) is essential for cardiac function by regulating junctional membrane complex activity. Circ. Res.120, 110–119. doi: 10.1161/CIRCRESAHA.116.309977

  • 28

    RahmelT.MarkoB.NowakH.BergmannL.ThonP.RumpK.et al. (2020). Mitochondrial dysfunction in sepsis is associated with diminished intramitochondrial TFAM despite its increased cellular expression. Sci. Rep.10, 21029. doi: 10.1038/s41598-020-78195-4

  • 29

    RannikkoJ.SeiskariT.HuttunenR.TarkiainenI.JylhavaJ.HurmeM.et al. (2018). Plasma cell-free DNA and qSOFA score predict 7-day mortality in 481 emergency department bacteraemia patients. J. Intern. Med.284, 418–426. doi: 10.1111/joim.12766

  • 30

    ReznikE.MillerM. L.SenbabaogluY.RiazN.SarungbamJ.TickooS. K.et al. (2016). Mitochondrial DNA copy number variation across human cancers. Elife.5, e10769. doi: 10.7554/eLife.10769

  • 31

    RoscaM.MinklerP.HoppelC. L. (2011). Cardiac mitochondria in heart failure: normal cardiolipin profile and increased threonine phosphorylation of complex IV. Biochim. Biophys. Acta1807, 1373–1382. doi: 10.1016/j.bbabio.2011.02.003

  • 32

    ShenX.HanG.LiS.SongY.ShenH.ZhaiY.et al. (2018). Association between the T6459C point mutation of the mitochondrial MT-CO1 gene and susceptibility to sepsis among Chinese Han people. J. Cell Mol. Med.22, 5257–5264. doi: 10.1111/jcmm.13746

  • 33

    ShenoyS. (2020). Coronavirus (Covid-19) sepsis: revisiting mitochondrial dysfunction in pathogenesis, aging, inflammation, and mortality. Inflamm. Res.69, 1077–1085. doi: 10.1007/s00011-020-01389-z

  • 34

    SinghP.ParajuliN.MayeuxP. R.MacMillan-CrowL. A. (2016). Renal mitochondrial lipid peroxidation during sepsis. J. Kidney.2, 116. doi: 10.4172/jok.1000116

  • 35

    TonialC. T.CostaC. A. D.AndradesG. R. H.CrestaniF.BrunoF.PivaJ. P.et al. (2021). Performance of prognostic markers in pediatric sepsis. J. Pediatr. (Rio J).97, 287–294. doi: 10.1016/j.jped.2020.07.008

  • 36

    WangX.LiR.QianS.YuD. (2023). Multilevel omics for the discovery of biomarkers in pediatric sepsis. Pediatr. Invest.7, 277–289. doi: 10.1002/ped4.12405

  • 37

    WeissS. L.SelakM. A.TulucF.Perales VillarroelJ.NadkarniV. M.DeutschmanC. S.et al. (2015). Mitochondrial dysfunction in peripheral blood mononuclear cells in pediatric septic shock. Pediatr. Crit. Care Med.16, e4–e12. doi: 10.1097/PCC.0000000000000277

  • 38

    WeissS. L.ZhangD.BushJ.GrahamK.StarrJ.MurrayJ.et al. (2020). Mitochondrial dysfunction is associated with an immune paralysis phenotype in pediatric sepsis. Shock.54, 285–293. doi: 10.1097/SHK.0000000000001486

  • 39

    YanH. P.LiM.LuX. L.ZhuY. M.Ou-YangW. X.XiaoZ. H.et al. (2018). Use of plasma mitochondrial DNA levels for determining disease severity and prognosis in pediatric sepsis: a case control study. BMC Pediatr.18, 267. doi: 10.1186/s12887-018-1239-z

  • 40

    YangY.YangJ.YuB.LiL.LuoL.WuF.et al. (2019). Association between circulating mononuclear cell mitochondrial DNA copy number and in-hospital mortality in septic patients: A prospective observational study based on the Sepsis-3 definition. PLoS One14, e0212808. doi: 10.1371/journal.pone.0212808

  • 41

    YueP.JingS.LiuL.MaF.ZhangY.WangC.et al. (2018). Association between mitochondrial DNA copy number and cardiovascular disease: Current evidence based on a systematic review and meta-analysis. PLoS One13, e0206003. doi: 10.1371/journal.pone.0206003

  • 42

    YueP.ZhangY.LiuL.ZhouK.XiaS.PengM.et al. (2022). Yap1 modulates cardiomyocyte hypertrophy via impaired mitochondrial biogenesis in response to chronic mechanical stress overload. Theranostics.12, 7009–7031. doi: 10.7150/thno.74563

  • 43

    ZhangD.LiY.Heims-WaldronD.BezzeridesV.GuatimosimS.GuoY.et al. (2018). Mitochondrial cardiomyopathy caused by elevated reactive oxygen species and impaired cardiomyocyte proliferation. Circ. Res.122, 74–87. doi: 10.1161/CIRCRESAHA.117.311349

  • 44

    ZhangW. Y.ChenZ. H.AnX. X.LiH.ZhangH. L.WuS. J.et al. (2023). Analysis and validation of diagnostic biomarkers and immune cell infiltration characteristics in pediatric sepsis by integrating bioinformatics and machine learning. World J. pediatrics: WJP.19, 1094–1103. doi: 10.1007/s12519-023-00717-7

Summary

Keywords

sepsis, mitochondrial DNA copy number, MT-CO1, Mt-ND1, MT-ATP6, prognosis

Citation

Jing S, Zhang Y, Zhao W, Li Y and Wen Y (2024) The predictive value of peripheral blood cell mitochondrial gene expression in identifying the prognosis in pediatric sepsis at preschool age. Front. Cell. Infect. Microbiol. 14:1413103. doi: 10.3389/fcimb.2024.1413103

Received

06 April 2024

Accepted

08 July 2024

Published

24 July 2024

Volume

14 - 2024

Edited by

Sharmistha Banerjee, University of Hyderabad, India

Reviewed by

Dan Yu, Capital Medical University, China

Xiaoling Guo, Wenzhou Medical University, China

Siyi He, General Hospital of Western Theater Command, China

Updates

Copyright

*Correspondence: Yan Wen, ; Yifei Li,

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics