Next Article in Journal
CAR-T-Cell-Based Cancer Immunotherapies: Potentials, Limitations, and Future Prospects
Next Article in Special Issue
Postoperative, but Not Preoperative, MELD-3.0 Prognosticates 3-Month Procedural Success in Patients Undergoing Orthotopic Heart Transplantation
Previous Article in Journal
Effects of Hypoxia and Inflammation on Hepcidin Concentration in Non-Anaemic COVID-19 Patients
Previous Article in Special Issue
BMI-Stratified Exploration of the ‘Obesity Paradox’: Heart Failure Perspectives from a Large German Insurance Database
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Exercise Immune Fitness Test to Unravel Disease Mechanisms—A Proof-of-Concept Heart Failure Study

by
Galyna Bondar
1,
Abhinandan Das Mahapatra
2,
Tra-Mi Bao
1,
Irina Silacheva
1,
Adrian Hairapetian
1,
Thomas Vu
1,
Stephanie Su
1,
Ananya Katappagari
1,
Liana Galan
1,
Joshua Chandran
1,
Ruben Adamov
1,
Lorenzo Mancusi
1,
Isabel Lai
1,
Anca Rahman
1,
Tristan Grogan
1,
Jeffrey J. Hsu
1,
Monica Cappelletti
1,
Peipei Ping
1,
David Elashoff
1,
Elaine F. Reed
1 and
Mario C. Deng
1,*
add Show full author list remove Hide full author list
1
David Geffen School of Medicine, University of California Los Angeles Medical Center, Los Angeles, CA 90095, USA
2
Strand NGS, Strand Life Sciences Pvt. Ltd., San Francisco, CA 94104, USA
*
Author to whom correspondence should be addressed.
J. Clin. Med. 2024, 13(11), 3200; https://doi.org/10.3390/jcm13113200
Submission received: 14 May 2024 / Revised: 22 May 2024 / Accepted: 23 May 2024 / Published: 29 May 2024

Abstract

:
Background: Cardiorespiratory fitness positively correlates with longevity and immune health. Regular exercise may provide health benefits by reducing systemic inflammation. In chronic disease conditions, such as chronic heart failure and chronic fatigue syndrome, mechanistic links have been postulated between inflammation, muscle weakness, frailty, catabolic/anabolic imbalance, and aberrant chronic activation of immunity with monocyte upregulation. We hypothesize that (1) temporal changes in transcriptome profiles of peripheral blood mononuclear cells during strenuous acute bouts of exercise using cardiopulmonary exercise testing are present in adult subjects, (2) these temporal dynamic changes are different between healthy persons and heart failure patients and correlate with clinical exercise-parameters and (3) they portend prognostic information. Methods: In total, 16 Heart Failure (HF) patients and 4 healthy volunteers (HV) were included in our proof-of-concept study. All participants underwent upright bicycle cardiopulmonary exercise testing. Blood samples were collected at three time points (TP) (TP1: 30 min before, TP2: peak exercise, TP3: 1 h after peak exercise). We divided 20 participants into 3 clinically relevant groups of cardiorespiratory fitness, defined by peak VO2: HV (n = 4, VO2 ≥ 22 mL/kg/min), mild HF (HF1) (n = 7, 14 < VO2 < 22 mL/kg/min), and severe HF (HF2) (n = 9, VO2 ≤ 14 mL/kg/min). Results: Based on the statistical analysis with 20–100% restriction, FDR correction (p-value 0.05) and 2.0-fold change across the three time points (TP1, TP2, TP3) criteria, we obtained 11 differentially expressed genes (DEG). Out of these 11 genes, the median Gene Expression Profile value decreased from TP1 to TP2 in 10 genes. The only gene that did not follow this pattern was CCDC181. By performing 1-way ANOVA, we identified 8/11 genes in each of the two groups (HV versus HF) while 5 of the genes (TTC34, TMEM119, C19orf33, ID1, TKTL2) overlapped between the two groups. We found 265 genes which are differentially expressed between those who survived and those who died. Conclusions: From our proof-of-concept heart failure study, we conclude that gene expression correlates with VO2 peak in both healthy individuals and HF patients, potentially by regulating various physiological processes involved in oxygen uptake and utilization during exercise. Multi-omics profiling may help identify novel biomarkers for assessing exercise capacity and prognosis in HF patients, as well as potential targets for therapeutic intervention to improve VO2 peak and quality of life. We anticipate that our results will provide a novel metric for classifying immune health.

1. Introduction

All living organisms have developed highly conserved and regulated stress response mechanisms. Strenuous activity leads to the acute activation of the hypothalamic–pituitary–adrenal axis including the adrenergic system [1]. Cardiorespiratory fitness positively correlates with longevity and inversely with all-cause mortality [2,3,4,5,6,7,8,9,10]. The World Health Organization and the American Heart Association are working towards the implementation of cardiorespiratory fitness measurements to improve risk classification and optimize disease prevention [11,12,13,14,15,16].
Cardiorespiratory fitness, defined as the ability of the circulatory and respiratory systems to supply oxygen to skeletal muscles during sustained physical activity correlates with immune health. Maintaining cardiorespiratory fitness exerts beneficial immunological effects [17,18]. It has the potential to delay, halt, or even reverse the age-associated decline in immune function [19,20]. Fitness is hypothesized to be linked to leukocyte gene expression in both healthy individuals and Heart Failure (HF) patients by regulating various physiological processes involved in oxygen uptake and utilization during exercise.
Various data suggest that regular exercise may provide health benefits by reducing systemic inflammation [21]. In chronic disease conditions, such as chronic HF and chronic fatigue syndrome, mechanistic links have been postulated between inflammation, muscle weakness, frailty, catabolic/anabolic imbalance, and aberrant chronic activation of immunity with monocyte upregulation [22,23,24,25,26,27,28]. In these conditions, in which the adrenergic system is permanently activated, levels of circulating cell-free mitochondrial DNA are chronically increased and correlate with secondary organ dysfunction and cell injury. This altered pattern is associated with a pro-inflammatory equilibrium [29,30,31,32].
Cardiorespiratory fitness is measured by peak oxygen uptake (VO2) in mL per kilogram body weight per minute during exercise by standardized cardiopulmonary exercise testing (CPX) [33]. CPX has become an important clinical tool to evaluate exercise capacity and predict outcomes in patients with HF and other cardiac conditions. It provides an assessment of the integrative exercise responses involving the pulmonary, cardiovascular, and skeletal muscle systems, which are not adequately reflected through the measurement of individual organ system function.
Our study aims to identify the gap in the detailed molecular signals induced by exercise that might benefit health and prevent disease. While the Molecular Transducers of Physical Activity Consortium is studying healthy individuals with no prior disease condition, there is a gap in focusing on subjects with specific diseases such as HF which this study addresses.
Based on our overall hypothesis that leukocyte gene expression in both healthy volunteers (HV) and Heart Failure (HF) is linked to cardiorespiratory fitness, this paper presents a proof-of-concept data analysis to characterize the relationship between exercise physiology as reflected in the capacity for functional exercise testing by CPX and the whole-transcriptome dynamics of mixed populations of peripheral blood mononuclear cells (PBMC) to characterize the phenotype of healthy persons and HF patients. We specifically hypothesize that (1) temporal changes in transcriptome profiles of PBMC during strenuous acute bouts of exercise using CPX are present in adult subjects, (2) these temporal dynamic changes are different between HV and HF patients and correlate with clinical CPX parameters and (3) they portend prognostic information. We believe that our results may be also applicable to various chronic disease conditions, such as other forms of organ failure, including immune system failure such as Post-Acute Sequelae of COVID (PASC), also called long-Covid, ref. [34] and chronic fatigue syndrome.

2. Methods and Design

2.1. Subject Population

Twenty persons were included in our proof-of-concept study (4 HV and 16 HF patients) who were evaluated for advanced cardiac care options and consented to UCLA IRB 12-000351. All participants underwent CPX. Blood samples were collected at three Time Points (TP) (Figure 1).
Subject Recruitment Inclusion criteria: On the basis of NIH-screening criteria (https://clinicalcenter.nih.gov/recruit/index.html accessed on 16 October 2019), we defined for the purpose of our proof-of-concept study HV as persons with no known significant health problems, specifically without any known cardiac, respiratory, or metabolic disease, not taking any chronic prescribed medication. An HF subject was defined as a subject with HF referred to our Advanced Heart Failure and Heart Transplant Program for evaluation of advanced cardiac therapies including heart transplantation.
Multivariate CPX panel All participants underwent CPX testing using a standardized bicycle ergometer Ramp protocol until each study participant achieved their individual maximum oxygen consumption or peak oxygen uptake (peak VO2). The multivariate CPX panel is a validated measure of cardiorespiratory fitness. Peak VO2 is reached when oxygen consumption remains at a steady state despite an increase in workload, measured as Respiratory Exchange Ratio (RER) > 1.0.
Clinical data were collected from questionnaires for HV and from medical records of HF on the day of CPX testing. We divided 20 participants into 3 clinically relevant groups of cardiorespiratory fitness, defined by peak VO2: HV (n = 4, VO2 ≥ 22 mL/kg/min), mild HF (HF1) (n = 7, 14 < VO2 < 22 mL/kg/min), and severe HF (HF2) (n = 9, VO2 ≤ 14 mL/kg/min).

2.2. Blood Processing

Samples were collected at 3 time points (TP): within 30 min before exercise (TP1), within 60 s of peak exercise (TP2), and within 1 h post-exercise (TP3) into 1 CPT (Becton Dickinson, Franklin Lakes, NJ, USA) tube for RNA-seq analyses.
RNA analyses Total RNA was isolated from PBMC for each blood sample. Purified RNA quality was verified on an Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA). RNA concentrations were determined using a NanoDrop ND-1000 spectrophotometer (NanoDrop Technologies, Wilmington, DE, USA). The mRNA library was prepared with Illumina TruSeq RNA kit (Illumina, San Diego, CA, USA). The cDNA libraries were quantified using Qubit and sequenced on Illumina HiSeq 2500 (Illumina, San Diego, CA, USA). Total mRNA was amplified and sequenced on the whole-genome Illumina HiSeq 3000.
PBMC sample processing and Gene Expression Profile protocol. Eight mLl of blood was drawn into a CPT tube. Peripheral Blood Mononuclear cells (PBMC) from each sample were purified within 2 h of phlebotomy. The collected blood was mixed and centrifuged at room temperature (22 °C) for 20 min at 2000× g RCF. Two mLl of plasma was separated without disturbing the cell layer into an Eppendorf tube (Sigma-Aldrich, St. Louis, MO, USA) and stored at −80 °C for future experiments. The cell layer was collected, transferred to 15 mL conical tubes, and re-suspended in cold Phosphate Buffer Saline (PBS) (Sigma-Aldrich, St. Louis, MO, USA) and centrifuged for 20 min at 300× g RCF at 4 °C. The supernatant was aspirated and discharged. The cell pellet was re-suspended in cold PBS, transferred into an Eppendorf tube, and centrifuged for 20 min at 300× g RCF at 4 °C. The supernatant was discharged. The pellet was re-suspended in 0.5 mL RNA Protect Cell Reagent (Qiagen, Valencia, CA, USA) and frozen at −80° C.
PBMC transcriptome RNA sequencing. All samples were processed using next-generation RNA sequencing transcriptome analysis at the UCLA Technology Center for Genomics and Bioinformatics. Briefly, the RNA was isolated from the PBMC using RNeasy Mini Kit (Qiagen, Valencia, CA, USA). The quality of the total RNA was assessed using NanoDrop® ND-1000 spectrophotometer (NanoDrop Technologies, Wilmington, DE, USA) and Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA) concentration above 50 ng/µL, purity—260/280~2.0, integrity—RIN > 9.0 and average > 9.5. Then, mRNA library was prepared with Illumina TruSeq RNA kit according to the manufacturer’s instructions (Illumina, San Diego, CA, USA). Library construction consists of random fragmentation of the polyA mRNA, followed by cDNA production using random polymers. The cDNA libraries were quantified using Qubit and size distribution was checked on Bioanalyzer 2100 (Agilent Technologies, Palo Alto, CA, USA). The library was sequenced on HiSeq 2500. Clusters were generated to yield approximately 725 K–825 K clusters/mm2. Cluster density and quality were determined during the run after the first base addition parameters were assessed. We performed single-end sequencing runs to align the cDNA sequences to the reference genome. Generated FASTQ files were transferred to the Advanced HF Research Data Center where Avadis NGS 1.5 (Agilent, Palo Alto, CA, USA and Strand Scientific, San Francisco, CA, USA) was used to align the raw RNA-Seq FASTQ reads to the reference genome. After RNA extraction, quantification and quality assessment, total mRNA was amplified and sequenced on the whole-genome Illumina HiSeq 2500. Data were then subjected to DeSeq normalization using NGS Strand/Avadis (v2.1 10 October 2014).

2.3. Statistical Analysis

The 60 single-end samples were aligned with Human Genome (hg38 build) and RNA Seq analysis was performed on the Strand NGS software V4.0. Reads were normalized using DESeq normalization method. We performed filter by expression (20th–100th percentile) on the samples to remove the genes that had very low normalized signal values in all the samples and retained 42,610 entities out of 57,923 entities. We performed fold change analysis with fold change cut off 2.0 and analyzed for all three TP based on all against the single condition. We identified 156 genes expressed in the conditions that satisfied the fold change cut-off. One-way ANOVA was performed on these entities with all TP and the 60 samples. Based on our extensive molecular prediction test development experience, we used the corrected p-value cut-off of 0.05, fold change cut-off of 2.0 and Benjamini–Hochberg as multiple testing corrections, retaining 11 genes expressed across all TPs. Subsequently, covariate regression analysis was performed using Pearson similarity metric. In addition, we performed ANOVA for all TP and oxygen uptake in HV, mild HF, and severe HF subjects separately. We were able to identify 5 out of these 11 genes expressed in both healthy and HF subjects.

3. Results

3.1. Differentially Expressed Genes (DEGs) before/at Max/1 h Post Max Exercise in the Combined Cohort

1.1: First, we were interested in the effect of CPX-based strenuous exercise on the temporal dynamics of PBMC-Gene Expression Profile in the combined cohort (n = 20). Based on the statistical analysis with 20–100% restriction, FDR-correction (p-value 0.05) and 2.0-fold change across the three time points (TP1, TP2, TP3) criteria, we obtained 11 DEGs using the analytical criteria detailed in the Section 2.3. The biological functions of these 11 genes are summarized in Table S1 (Supplementary Materials).
1.2: With respect to the directionality of gene expression, we observed that the 11 genes, when ranked by their median expression for each individual study participant and at each of the three TP, showed lower median signal values at TP2 and TP3 compared to TP1. Therefore, we were interested in the temporal dynamics of each of the 11 individual gene activity levels averaged over the cohort of all 20 subjects at each TP (Figure 2).
Out of these 11 genes, the median Gene Expression Profile value decreased from TP1 to TP2 in 10 genes. The only gene that did not follow this pattern was CCDC181 (circled in red). From TP2 to TP3, some GEP reversed their activity to TP1 level (Figure 3).
1.3: Next, we were interested in what extent the 11 strenuous-activity-bout-associated genes identified in the overall cohort were differentially expressed between HV (n = 4) and HF (n = 16) subjects, applying the same statistical criteria (20–100% restriction, FDR-correction (p-value 0.05) and 2.0-fold change across the three TPs). By performing 1-way ANOVA, we identified 8/11 genes in each of the two groups (HV versus HF) while 5 of the genes (TTC34, TMEM119, C19orf33, ID1, TKTL2) overlapped between the two groups as represented in the Venn diagram, suggesting some biological regulation differences in the PBMC pool (Figure 4). According to further clustering analyses with the 11 genes and their expression in all 20 subjects, we concluded that they did not have a close correlation to the clinical phenotypes (HV, HF1 and HF2 cohorts) at TP1 (resting state) and TP3 (Recovery state) but showed some correlation with HV/HF-state at peak exercise (Figure 5). At peak exercise (TP2), there were two main clusters differentiating severe HF from mild HF and HV and other cluster-mild HF and HV groups.

3.2. Differential Behavior of Genes in the HV versus Mild HF1 and Severe HF2 Cohorts (Peak VO2/Percent Predicted VO2)

2.1: Next, we were interested in understanding the correlation between the PBMC-transcriptome and the major clinical determinants of fitness, i.e., the condition of health, mild or severe HF, during acute bouts of strenuous exercise using CPX. The exploratory analysis suggested positive and inverse correlations. For exploratory reasons, we performed cluster entities on all 20 subjects for peak VO2 and percent predicted VO2 with the entity list that we derived from the previous fold change analysis (1.5/2) (Figure 6) which showed differences in TP2 samples between mild HF and severe HF groups. With the FC cut off 1.5, the samples for TP2 showed differences between mild HF and severe HF groups.

3.3. Correlations between the DEG/Clinical Profiles (Peak VO2/Percent Predicted VO2) and Survival/Heart Transplant Outcomes

3.1: Next, we were interested in understanding the correlation between genes associated with death in the severe HF group (n-9), analyzed by time-point (TP1, TP2, TP3) and as shared gene set represented by the Venn diagram (Figure 7). We found 265 genes which are differentially expressed between those who survived and those who died.
3.2: In order to provide additional support for our hypothesis that the currently identified genes and our previously described and UCLA-patented survival-related PBMC-genes are reflective of an underlying biological process, we compared the DEG-list for those three out of nine severe HF patients who died during follow-up (n = 265) with a DEG list derived from 29 AdHF-patient who had undergone MCS during 2012–2014 out of whom 11 died within a year (n = 105).
To reinforce our hypothesis regarding the association of identified genes, including our patented survival-related genes with an underlying biological process, we compared DEG lists. These lists were obtained from three out of nine severe HF patients who died during follow-up in the current study (n = 265) and 29 advanced HF patients, 11 of whom died within a year after undergoing mechanical circulatory support surgery in a previous study (n = 105) (Bondar, 2017 [35]). Out of these sets of DEG, six genes overlapped (SPOCD1, MUC20, OLFM1, CACNA1A, DDX3Y, BCORP1) (Table S2). The direction of regulation was the same in five out of six DEGs. The overlap of six DEGs suggests that there may be shared underlying immune-mediated processes related to the risk of dying from Advanced HF requiring further study.

4. Discussion

In our proof-of-concept study, we present data to support our hypotheses that (1) strenuous bouts of acute exercise during CPX induce longitudinal dynamic changes in transcriptome profiles of PBMC in adult subjects, (2) these longitudinal dynamic changes are different between HV and HF patients and correlate with clinical CPX-parameters and (3) they portend potential prognostic information.
First, our data show that longitudinal dynamic changes in PBMC-transcriptome profiles during acute bouts of strenuous exercise using CPX are present in adult subjects confirming prior studies conducted in younger people. Prior studies had suggested that brief exercise alters PBMC gene expression in early- and late-pubertal children. The pattern of change involves diverse genetic pathways, consistent with a global danger-type response, perhaps readying PBMCs for a range of physiological functions from inflammation to tissue repair that, to our knowledge, would be useful following a bout of physical activity [36].
Second, our results suggest that these longitudinal dynamic changes of each of the 11 DEGs may differ between HV and HF patients. The question regarding biological relevance needs to equally address (A) why some genes overlap and (B) why some genes do not overlap. We hypothesize, for future studies, that the key to discriminating adaptive versus maladaptive mechanisms of immune health may lie in understanding the patterns of dynamic cell free-mitochondrial DNA clearance and its subsequent effect on the transcriptome and proteome. Of our 11 DEGs, DPYD-AS1, TNFRSF12A, and CD300LD were found to play a role in inflammation and in the immune system [37,38,39]. TNFRSF12A was also found to be highly inducible and to play a key role in the development of cardiac hypertrophy followed by HF [40]. We postulate that the temporal dynamics of total cell free-DNA, specifically cell free-mitochondrial DNA, may either be adaptive or maladaptive, depending on the body’s ability to produce and eliminate these molecules. Therefore, in future studies underway in our lab, we propose to characterize the temporal dynamics of the immune system under exercise and resting conditions through the integrated analysis of cell free-DNA quantification (genomic), mRNA expression (transcriptomic), PBMC sub-population and cytokine production (proteomic), as well as clinical (phenomics) data from the multivariate CPX panel.
Third, we identified PBMC transcripts that might portend prognostic outcome prediction information. The fact that mortality/survival-related PBMC transcripts and CPX-related PBMC transcripts did only marginally overlap raises the interesting speculative hypothesis that these phenomena are governed by different biological pathways.
Fourth, differences in gene expression related to VO2 peak between healthy individuals and those with HF reflect the underlying pathophysiological changes associated with HF, including alterations in energy metabolism, mitochondrial function, oxygen transport, inflammatory responses, neurohormonal signaling, and structural remodeling. The difference in gene expression related to peak oxygen uptake (VO2 peak) between healthy individuals and those with HF is complex and multifactorial. Several studies have investigated gene expression profiles in cardiac and skeletal muscle tissues to understand the molecular mechanisms underlying exercise capacity and its alterations in HF.
We anticipate that our results will provide a novel metric for classifying immune health. We will use innovative multi-omics analyses to better understand this complex evolutionary physiology, which integrates data from different biological layers. Multi-omics profiling may help identify novel biomarkers for assessing exercise capacity and prognosis in HF patients, as well as potential targets for therapeutic intervention to improve VO2 peak and quality of life.

Limitations

First, our proof-of-concept study is very small with 4 healthy volunteers and 16 HF subjects. This limitation implies that the presented data should serve as building blocks for subsequent confirmatory studies. Second, the correlation analyses do not allow any inferences about causality relations between gene expression levels and clinical phenotypes. In addition to larger confirmatory clinical-translational studies, we therefore anticipate the necessity for mechanistic studies.

5. Conclusions

In summary, gene expression can influence VO2 peaks in both healthy individuals and HF patients by regulating various physiological processes involved in oxygen uptake and utilization during exercise. We found 11 differentially expressed genes that can illustrate changes in activity under exercise conditions and that also can characterize the clinical phenotypes at VO2 peak. We also found 265 genes associated with mortality, of which 6 overlap with genes that we found in a mortality cohort in a previous study conducted in our lab. As a whole, gene expression can influence VO2 peaks in both healthy individuals and HF patients by regulating various physiological processes involved in VO uptake and utilization during exercise. Dysregulation of gene expression in HF can contribute to reduced VO2 peak and exercise intolerance, highlighting the importance of understanding the molecular mechanisms underlying these alterations for developing targeted therapeutic strategies. Further research is needed to elucidate the specific gene expression signatures associated with exercise intolerance in HF to identify potential therapeutic targets for improving exercise capacity in these patients and to test the utility of an immune fitness test to predict outcomes in various disease conditions.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/jcm13113200/s1, Table S1: Gene ontology information on the 11 genes identified as described in the text; Table S2: Gene ontology information on the six overlapping genes. References [37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60] are cited in the Supplementary Materials.

Author Contributions

M.C.D. is the Principal Investigator of the Deng Advanced Heart Failure Research Lab and developed the hypothesis, designed the study, and participated in every stage of the manuscript development. G.B., A.D.M., T.-M.B. and I.S. co-designed the study and co-drafted the manuscript. D.E. led the plan for statistical analysis. All authors, including A.D.M., I.S., A.H., T.V., S.S., A.K., L.G., J.C., R.A., L.M., I.L., A.R., T.G., J.J.H., M.C., P.P. and E.F.R., contributed to drafting and reviewing the final manuscript. All authors approved the final version of the manuscript. All authors have agreed to be personally accountable for the author’s own contributions and to ensure that questions related to the accuracy or integrity of any part of the work, even ones in which the author was not personally involved, are appropriately investigated, resolved, and the resolution documented in the literature. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by UCLA National Institutes of Health [R21 1R21HL120040-01 to M.C.D. and G.B.], UCLA Department of Medicine Internal Funds and the Advanced Heart Failure Research Gift to Columbia University (Philip Geier, John Tocco and Robert Milo) and Advanced Heart Failure Research Gift to UCLA (Juan Mulder, Peter Schultz, Larry Layne, James and Candace Moose).

Institutional Review Board Statement

This study was approved by the Institutional Review Board of the University of California, Los Angeles (approval number: 12-000351) and was conducted in accordance with the Declaration of Helsinki.

Informed Consent Statement

Written informed consent was obtained from each patient.

Data Availability Statement

Data is contained within the article and Supplementary Materials.

Acknowledgments

We would like to acknowledge the contributions of Maral Bakir, David Beenhouwer, Talha Demirel, Diego Gonzales, Victoria Groysberg, Megan Kamath, Berke Korucu, Isabella Macias, Candace Moose, James Moose, Sammy Morales, Kaitlyn Ngai, Federica Raia, Joanna Schaenman, Kishan Shah, Karishma Shah, Natalie Simonian, Darko Vucicevic and Wei Wang.

Conflicts of Interest

The authors G.B. and M.C.D. are co-founding equity holders of LeukoLifeDx, Inc., the developer of MyLeukoMAPTM biomarker test. T.-M.B. was an employee of LeukoLifeDx. The authors have no other relevant affiliations or financial involvement with any organization or entity with a financial interest in or financial conflict with the subject matter or materials discussed in the manuscript apart from those disclosed.

References

  1. Ortega, E. The “bioregulatory effect of exercise” on the innate/inflammatory responses. J. Physiol. Biochem. 2016, 72, 361–369. [Google Scholar] [CrossRef] [PubMed]
  2. Clausen, J.S.R.; Marott, J.L.; Holtermann, A.; Gyntelberg, F.; Jensen, M.T. Midlife cardiorespiratory fitness and the long-term risk of mortality: 46 years of follow-up. J. Am. Coll. Cardiol. 2018, 72, 987–995. [Google Scholar] [CrossRef] [PubMed]
  3. Blair, S.N.; Kohl, H.W.; Paffenbarger, R.S.; Clark, D.G.; Cooper, K.H.; Gibbons, L.W. Physical fitness and all-cause mortality: A prospective study of healthy men and women. JAMA 1989, 262, 2395–2401. [Google Scholar] [CrossRef] [PubMed]
  4. Kodama, S.; Saito, K.; Tanaka, S.; Maki, M.; Yachi, Y.; Asumi, M.; Sugawara, A.; Totsuka, K.; Shimano, H.; Ohashi, Y.; et al. Cardiorespiratory fitness as a quantitative predictor of all-cause mortality and cardiovascular events in healthy men and women: A meta-analysis. JAMA 2009, 301, 2024–2035. [Google Scholar] [CrossRef] [PubMed]
  5. Ekelund, L.-G.; Haskell, W.L.; Johnson, J.L.; Whaley, F.S.; Criqui, M.H.; Sheps, D.S.; Lipid Research Clinics Mortality Follow-up, S. Physical fitness as a predictor of cardiovascular mortality in asymptomatic North American men. N. Engl. J. Med. 1988, 319, 1379–1384. [Google Scholar] [CrossRef] [PubMed]
  6. Blair, S.N.; Kohl, H.W.; Barlow, C.E.; Paffenbarger, R.S.; Gibbons, L.W.; Macera, C.A. Changes in physical fitness and all-cause mortality: A prospective study of healthy and unhealthy men. JAMA 1995, 273, 1093–1098. [Google Scholar] [CrossRef]
  7. Harber, M.P.; Kaminsky, L.A.; Arena, R.; Blair, S.N.; Franklin, B.A.; Myers, J.; Ross, R. Impact of cardiorespiratory fitness on all-cause and disease-specific mortality: Advances since 2009. Prog. Cardiovasc. Dis. 2017, 60, 11–20. [Google Scholar] [CrossRef] [PubMed]
  8. Kavanagh, T.; Mertens, D.J.; Hamm, L.F.; Beyene, J.; Kennedy, J.; Corey, P.; Shephard, R.J. Prediction of long-term prognosis in 12 169 men referred for cardiac rehabilitation. Circulation 2002, 106, 666–671. [Google Scholar] [CrossRef] [PubMed]
  9. Arem, H.; Moore, S.C.; Patel, A.; Hartge, P.; De Gonzalez, A.B.; Visvanathan, K.; Campbell, P.T.; Freedman, M.; Weiderpass, E.; Adami, H.O.; et al. Leisure time physical activity and mortality: A detailed pooled analysis of the dose-response relationship. JAMA Intern. Med. 2015, 175, 959–967. [Google Scholar] [CrossRef]
  10. Myers, J.; Prakash, M.; Froelicher, V.; Do, D.; Partington, S.; Atwood, J.E. Exercise capacity and mortality among men referred for exercise testing. N. Engl. J. Med. 2002, 346, 793–801. [Google Scholar] [CrossRef]
  11. Ross, R.; Blair, S.N.; Arena, R.; Church, T.S.; Després, J.-P.; Franklin, B.A.; Haskell, W.L.; Kaminsky, L.A.; Levine, B.D.; Lavie, C.J.; et al. Importance of assessing cardiorespiratory fitness in clinical practice: A case for fitness as a clinical vital sign: A scientific statement from the American Heart Association. Circulation 2016, 134, e653–e699. [Google Scholar] [CrossRef] [PubMed]
  12. Lee, I.M.; Shiroma, E.J.; Lobelo, F.; Puska, P.; Blair, S.N.; Katzmarzyk, P.T. Effect of physical inactivity on major non-communicable diseases worldwide: An analysis of burden of disease and life expectancy. Lancet 2012, 380, 219–229. [Google Scholar] [CrossRef] [PubMed]
  13. Kraus, W.E.; Bittner, V.; Appel, L.; Blair, S.N.; Church, T.; Després, J.-P.; Franklin, B.A.; Miller, T.D.; Pate, R.R.; Taylor-Piliae, R.E.; et al. The National Physical Activity Plan: A call to action from the American Heart Association: A science advisory from the American Heart Association. Circulation 2015, 131, 1932–1940. [Google Scholar] [CrossRef] [PubMed]
  14. Walsh, N.P.; Gleeson, M.; Pyne, D.B.; Nieman, D.C.; Dhabhar, F.S.; Shephard, R.J.; Oliver, S.J.; Bermon, S.; Kajeniene, A. Position statement part two: Maintaining immune health. Exerc. Immunol. Rev. 2011, 17, 64–103. [Google Scholar] [PubMed]
  15. Walsh, N.P.; Gleeson, M.; Shephard, R.J.; Gleeson, M.; Woods, J.A.; Bishop, N.; Fleshner, M.; Green, C.; Pedersen, B.K.; Hoffman-Goete, L.; et al. Position statement part one: Immune function and exercise. Exerc. Immunol. Rev. 2011, 17, 6–63. [Google Scholar] [PubMed]
  16. World Health Organization. Global Recommendations on Physical Activity for Health; World Health Organization: Geneva, Switzerland, 2010. [Google Scholar]
  17. Gjevestad, G.O.; Holven, K.B.; Ulven, S.M. Effects of exercise on gene expression of inflammatory markers in human peripheral blood cells: A systematic review. Curr. Cardiovasc. Risk Rep. 2015, 9, 34. [Google Scholar] [CrossRef] [PubMed]
  18. Huber, Y.; Gehrke, N.; Biedenbach, J.; Helmig, S.; Simon, P.; Straub, B.K.; Bergheim, I.; Huber, T.; Schuppan, D.; Galle, P.R.; et al. Voluntary distance running prevents TNF-mediated liver injury in mice through alterations of the intrahepatic immune milieu. Cell Death Dis. 2017, 8, e2893. [Google Scholar] [CrossRef] [PubMed]
  19. Spielmann, G.; McFarlin, B.K.; O’Connor, D.P.; Smith, P.J.W.; Pircher, H.; Simpson, R.J. Aerobic fitness is associated with lower proportions of senescent blood T-cells in man. Brain Behav. Immun. 2011, 25, 1521–1529. [Google Scholar] [CrossRef] [PubMed]
  20. Turner, J.E. Is immunosenescence influenced by our lifetime “dose” of exercise? Biogerontology 2016, 17, 581–602. [Google Scholar] [CrossRef]
  21. Shockett, P.E.; Khanal, J.; Sitaula, A.; Oglesby, C.; Meachum, W.A.; Castracane, V.D.; Kraemer, R.R. Plasma cell-free mitochondrial DNA declines in response to prolonged moderate aerobic exercise. Physiol. Rep. 2016, 4, e12672. [Google Scholar] [CrossRef]
  22. Mann, D.L. Inflammatory mediators and the failing heart: Past, present, and the foreseeable future. Circ. Res. 2002, 91, 988–998. [Google Scholar] [CrossRef] [PubMed]
  23. Nemeth, E.; Valore, E.V.; Territo, M.; Schiller, G.; Lichtenstein, A.; Ganz, T. Hepcidin, a putative mediator of anemia of inflammation, is a type II acute-phase protein. Blood J. Am. Soc. Hematol. 2003, 101, 2461–2463. [Google Scholar] [CrossRef] [PubMed]
  24. Shirazi, L.F.; Bissett, J.; Romeo, F.; Mehta, J.L. Role of inflammation in heart failure. Curr. Atheroscler. Rep. 2017, 19, 27. [Google Scholar] [CrossRef] [PubMed]
  25. Coats, A.J.S. The “muscle hypothesis” of chronic heart failure. J. Mol. Cell. Cardiol. 1996, 28, 2255–2262. [Google Scholar] [CrossRef] [PubMed]
  26. Doehner, W.; Frenneaux, M.; Anker, S.D. Metabolic impairment in heart failure: The myocardial and systemic perspective. J. Am. Coll. Cardiol. 2014, 64, 1388–1400. [Google Scholar] [CrossRef] [PubMed]
  27. Larun, L.; Brurberg, K.G.; Odgaard-Jensen, J.; Price, J.R. Exercise therapy for chronic fatigue syndrome. Cochrane Database Syst. Rev. 2019. [Google Scholar] [CrossRef] [PubMed]
  28. Loy, B.D.; O’Connor, P.J.; Dishman, R.K. Effect of acute exercise on fatigue in people with ME/CFS/SEID: A meta-analysis. Med. Sci. Sports Exerc. 2016, 48, 2003. [Google Scholar] [CrossRef] [PubMed]
  29. Grazioli, S.; Pugin, J. Mitochondrial damage-associated molecular patterns: From inflammatory signaling to human diseases. Front. Immunol. 2018, 9, 342896. [Google Scholar] [CrossRef] [PubMed]
  30. Levine, B.; Kalman, J.; Mayer, L.; Fillit, H.M.; Packer, M. Elevated circulating levels of tumor necrosis factor in severe chronic heart failure. N. Engl. J. Med. 1990, 323, 236–241. [Google Scholar] [CrossRef]
  31. Bierhaus, A.; Humpert, P.M.; Nawroth, P.P. Linking stress to inflammation. Anesthesiol. Clin. N. Am. 2006, 24, 325–340. [Google Scholar] [CrossRef]
  32. Dhabhar, F.S. Effects of stress on immune function: The good, the bad, and the beautiful. Immunol. Res. 2014, 58, 193–210. [Google Scholar] [CrossRef] [PubMed]
  33. Mancini, D.M.; Eisen, H.; Kussmaul, W.; Mull, R.; Edmunds, L.H., Jr.; Wilson, J.R. Value of peak exercise oxygen consumption for optimal timing of cardiac transplantation in ambulatory patients with heart failure. Circulation 1991, 83, 778–786. [Google Scholar] [CrossRef] [PubMed]
  34. Deng, M.C. An Exercise Immune Fitness Test to Unravel Mechanisms of Post-Acute Sequelae of COVID-19; Taylor & Francis: Abingdon, UK, 2023; Volume 19, pp. 693–697. [Google Scholar]
  35. Bondar, G.; Togashi, R.; Cadeiras, M.; Schaenman, J.; Cheng, R.K.; Masukawa, L.; Hai, J.; Bao, T.-M.; Chu, D.; Chang, E.; et al. Association between preoperative peripheral blood mononuclear cell gene expression profiles, early postoperative organ function recovery potential and long-term survival in advanced heart failure patients undergoing mechanical circulatory support. PLoS ONE 2017, 12, e0189420. [Google Scholar] [CrossRef] [PubMed]
  36. Radom-Aizik, S.; Zaldivar, F., Jr.; Leu, S.-Y.; Cooper, D.M. A brief bout of exercise alters gene expression and distinct gene pathways in peripheral blood mononuclear cells of early-and late-pubertal females. J. Appl. Physiol. 2009, 107, 168–175. [Google Scholar] [CrossRef] [PubMed]
  37. Singh, P.; Brock, C.O.; Volden, P.A.; Hernandez, K.; Skor, M.; Kocherginsky, M.; Park, J.E.; Brady, M.J.; Conzen, S.D. Glucocorticoid receptor C h IP-sequencing of subcutaneous fat reveals modulation of inflammatory pathways. Obesity 2015, 23, 2286–2293. [Google Scholar] [CrossRef]
  38. Solskov, L.; Magnusson, N.E.; Kristiansen, S.B.; Jessen, N.; Nielsen, T.T.; Schmitz, O.; Bøtker, H.E.; Lund, S. Microarray expression analysis in delayed cardioprotection: The effect of exercise, AICAR, or metformin and the possible role of AMP-activated protein kinase (AMPK). Mol. Cell. Biochem. 2012, 360, 353–362. [Google Scholar] [CrossRef]
  39. Hsu, W.B.; Hsu, W.H.; Hung, J.S.; Shen, W.J.; Hsu, R.W. Transcriptome analysis of osteoblasts in an ovariectomized mouse model in response to physical exercise. Bone Jt. Res. 2018, 7, 601–608. [Google Scholar] [CrossRef] [PubMed]
  40. Yerra, V.G.; Batchu, S.N.; Kabir, G.; Advani, S.L.; Liu, Y.; Siddiqi, F.S.; Connelly, K.A.; Advani, A. Empagliflozin disrupts a Tnfrsf12a-mediated feed forward loop that promotes left ventricular hypertrophy. Cardiovasc. Drugs Ther. 2022, 36, 619–632. [Google Scholar] [CrossRef] [PubMed]
  41. Ross, K.A. Coherent somatic mutation in autoimmune disease. PLoS ONE 2014, 9, e101093. [Google Scholar] [CrossRef]
  42. Saif, M.W.; Smith, M.; Maloney, A. The first case of severe takotsubo cardiomyopathy associated with 5-fluorouracil in a patient with abnormalities of both dihydropyrimidine dehydrogenase (DPYD) and thymidylate synthase (TYMS) genes. Cureus 2016, 8, e783. [Google Scholar] [CrossRef]
  43. Schwarz, T.; Prieler, B.; Schmid, J.A.; Grzmil, P.; Neesen, J. Ccdc181 is a microtubule-binding protein that interacts with Hook1 in haploid male germ cells and localizes to the sperm tail and motile cilia. Eur. J. Cell Biol. 2017, 96, 276–288. [Google Scholar] [CrossRef]
  44. Kanamoto, T.; Mizuhashi, K.; Terada, K.; Minami, T.; Yoshikawa, H.; Furukawa, T. Isolation and characterization of a novel plasma membrane protein, osteoblast induction factor (obif), associated with osteoblast differentiation. BMC Dev. Biol. 2009, 9, 70. [Google Scholar] [CrossRef]
  45. Zhang, M.; Zhu, J.; Zhang, P.; Li, L.; Min, M.; Li, T.; He, W. Development and validation of cancer-associated fibroblasts-related gene landscape in prognosis and immune microenvironment of bladder cancer. Front. Oncol. 2023, 13, 1174252. [Google Scholar] [CrossRef]
  46. Dinh, P.; Peng, J.-H.; Tran, C.; Tran, T.; Pan, S.-L. Construction of miRNAs and gene expression profiles associated with ischemic cardiomyopathy: Bioinformatics analysis. Res. Sq. 2022. [Google Scholar] [CrossRef]
  47. Almutairi, M.H.; Alrubie, T.M.; Alamri, A.M.; Almutairi, B.O.; Alrefaei, A.F.; Arafah, M.M.; Alanazi, M.; Semlali, A. Cancer-testis gene biomarkers discovered in colon cancer patients. Genes 2022, 13, 807. [Google Scholar] [CrossRef]
  48. Li, Y.; Wang, D.W.; Chen, Y.; Chen, C.; Guo, J.; Zhang, S.; Sun, Z.; Ding, H.; Yao, Y.; Zhou, L.; et al. Genome-wide association and functional studies identify SCML4 and THSD7A as novel susceptibility genes for coronary artery disease. Arterioscler. Thromb. Vasc. Biol. 2018, 38, 964–975. [Google Scholar] [CrossRef] [PubMed]
  49. Heidecker, B.; Kittleson, M.M.; Kasper, E.K.; Wittstein, I.S.; Champion, H.C.; Russell, S.D.; Hruban, R.H.; Rodriguez, E.R.; Baughman, K.L.; Hare, J.M.; et al. Transcriptomic biomarkers for the accurate diagnosis of myocarditis. Circulation 2011, 123, 1174–1184. [Google Scholar] [CrossRef] [PubMed]
  50. Gan, L.; Yang, C.; Zhao, L.; Wang, S.; Gao, Z.; Ye, Y. Prognostic Biomarker SPOCD1 and Its Correlation with Immune Infiltrates in Colorectal Cancer. Biomolecules 2023, 13, 209. [Google Scholar] [CrossRef] [PubMed]
  51. Liu, D.; Yang, Y.; Yan, A.; Yang, Y. SPOCD1 accelerates ovarian cancer progression and inhibits cell apoptosis via the PI3K/AKT pathway. OncoTargets Ther. 2020, 13, 351–359. [Google Scholar] [CrossRef]
  52. Lin, Z.; Chen, L.; Wu, T.; Zhang, Y.; Huang, X.; Chen, Y.; Chen, J.; Xu, Y. Prognostic Value of SPOCD1 in Esophageal Squamous Cell Carcinoma: A Comprehensive Study Based on Bioinformatics and Validation. Front. Genet. 2022, 13, 872026. [Google Scholar] [CrossRef]
  53. Li, M.; Feng, Z.; Han, R.; Hu, B.; Zhang, R.; Wang, H. Paclitaxel promotes mTOR signaling-mediated apoptosis in esophageal cancer cells by targeting MUC20. Thorac. Cancer 2023, 14, 3089–3096. [Google Scholar] [CrossRef] [PubMed]
  54. Xue, B.; Guo, W.-M.; Jia, J.-D.; Kadeerhan, G.; Liu, H.-P.; Bai, T.; Shao, Y.; Wang, D.-W. MUC20 as a novel prognostic biomarker in ccRCC correlating with tumor immune microenvironment modulation. Am. J. Cancer Res. 2022, 12, 695. [Google Scholar] [PubMed]
  55. Le-Niculescu, H.; Roseberry, K.; Gill, S.S.; Levey, D.F.; Phalen, P.L.; Mullen, J.; Williams, A.; Bhairo, S.; Voegtline, T.; Davis, H.; et al. Precision medicine for mood disorders: Objective assessment, risk prediction, pharmacogenomics, and repurposed drugs. Mol. Psychiatry 2021, 26, 2776–2804. [Google Scholar] [CrossRef] [PubMed]
  56. Li, R.; Diao, H.; Zhao, F.; Xiao, S.; El Zowalaty, A.E.; Dudley, E.A.; Mattson, M.P.; Ye, X. Olfactomedin 1 deficiency leads to defective olfaction and impaired female fertility. Endocrinology 2015, 156, 3344–3357. [Google Scholar] [CrossRef] [PubMed]
  57. Lipman, A.R.; Fan, X.; Shen, Y.; Chung, W.K. Clinical and genetic characterization of CACNA1A-related disease. Clin. Genet. 2022, 102, 288–295. [Google Scholar] [CrossRef] [PubMed]
  58. Le Roux, M.; Barth, M.; Gueden, S.; de Cepoy, P.D.; Aeby, A.; Vilain, C.; Hirsch, E.; de Saint Martin, A.; Des Portes, V.; Lesca, G.; et al. CACNA1A-associated epilepsy: Electroclinical findings and treatment response on seizures in 18 patients. Eur. J. Paediatr. Neurol. 2021, 33, 75–85. [Google Scholar] [CrossRef] [PubMed]
  59. Shen, H.; Yanas, A.; Owens, M.C.; Zhang, C.; Fritsch, C.; Fare, C.M.; Copley, K.E.; Shorter, J.; Goldman, Y.E.; Liu, K.F. Sexually dimorphic RNA helicases DDX3X and DDX3Y differentially regulate RNA metabolism through phase separation. Mol. Cell 2022, 82, 2588–2603. [Google Scholar] [CrossRef]
  60. Dicke, A.-K.; Pilatz, A.; Wyrwoll, M.J.; Punab, M.; Ruckert, C.; Nagirnaja, L.; Aston, K.I.; Conrad, D.F.; Di Persio, S.; Neuhaus, N.; et al. DDX3Y is likely the key spermatogenic factor in the AZFa region that contributes to human non-obstructive azoospermia. Commun. Biol. 2023, 6, 350. [Google Scholar] [CrossRef]
Figure 1. Study Design.
Figure 1. Study Design.
Jcm 13 03200 g001
Figure 2. Average expression of 11 genes across all three time points and 20 subjects. Color-coding identifies upper quartile (blue) and lower quartile (red).
Figure 2. Average expression of 11 genes across all three time points and 20 subjects. Color-coding identifies upper quartile (blue) and lower quartile (red).
Jcm 13 03200 g002
Figure 3. Gene Expression of 11 genes at Time Points 1, 2, and 3. The left columns represent the normalized signal values for Time Point 1, the middle columns for Time Point 2, and the right columns for Time Point 3.
Figure 3. Gene Expression of 11 genes at Time Points 1, 2, and 3. The left columns represent the normalized signal values for Time Point 1, the middle columns for Time Point 2, and the right columns for Time Point 3.
Jcm 13 03200 g003
Figure 4. The Venn diagram depicts the 11 genes described in Figure 3 and Table S1. While 5 of the 11 genes overlap between the group HV and HF, 6 of the 11 genes are specific to either HV (n = 4) or HF (n = 16).
Figure 4. The Venn diagram depicts the 11 genes described in Figure 3 and Table S1. While 5 of the 11 genes overlap between the group HV and HF, 6 of the 11 genes are specific to either HV (n = 4) or HF (n = 16).
Jcm 13 03200 g004
Figure 5. Cluster entities diagram for all 20 subjects and TP1, TP2 and TP3 with FC cut-off 2.0 and 11 genes. Top to bottom: TTC34, DPYD-AS1, CCDC181, TKTL2, THSD7A, TMEML19, TNRFSF12A, CD300LD, MGAT5B, C19orf33, ID1.
Figure 5. Cluster entities diagram for all 20 subjects and TP1, TP2 and TP3 with FC cut-off 2.0 and 11 genes. Top to bottom: TTC34, DPYD-AS1, CCDC181, TKTL2, THSD7A, TMEML19, TNRFSF12A, CD300LD, MGAT5B, C19orf33, ID1.
Jcm 13 03200 g005
Figure 6. Covariate regression analysis shows an expression of the 11 genes across all three TPs for HV and HF groups in correlation to VO2max and % predicted oxygen uptake. The color coding depicts a positive (red) and a negative (blue) correlation. The darker the color the stronger the respective correlation. Top to bottom: TTC34, DPYD-AS1, CCDC181, TKTL2, THSD7A, TMEM119, TNFRSF12A, CD300LD, MGAT5B, C19orf33, ID1.
Figure 6. Covariate regression analysis shows an expression of the 11 genes across all three TPs for HV and HF groups in correlation to VO2max and % predicted oxygen uptake. The color coding depicts a positive (red) and a negative (blue) correlation. The darker the color the stronger the respective correlation. Top to bottom: TTC34, DPYD-AS1, CCDC181, TKTL2, THSD7A, TMEM119, TNFRSF12A, CD300LD, MGAT5B, C19orf33, ID1.
Jcm 13 03200 g006
Figure 7. HF group (n-9), analyzed by time-point (TP1, TP2, TP3) and as a shared gene set. Comparison of Z test results for death-associated heart failure group and three time points separately.
Figure 7. HF group (n-9), analyzed by time-point (TP1, TP2, TP3) and as a shared gene set. Comparison of Z test results for death-associated heart failure group and three time points separately.
Jcm 13 03200 g007
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Bondar, G.; Mahapatra, A.D.; Bao, T.-M.; Silacheva, I.; Hairapetian, A.; Vu, T.; Su, S.; Katappagari, A.; Galan, L.; Chandran, J.; et al. An Exercise Immune Fitness Test to Unravel Disease Mechanisms—A Proof-of-Concept Heart Failure Study. J. Clin. Med. 2024, 13, 3200. https://doi.org/10.3390/jcm13113200

AMA Style

Bondar G, Mahapatra AD, Bao T-M, Silacheva I, Hairapetian A, Vu T, Su S, Katappagari A, Galan L, Chandran J, et al. An Exercise Immune Fitness Test to Unravel Disease Mechanisms—A Proof-of-Concept Heart Failure Study. Journal of Clinical Medicine. 2024; 13(11):3200. https://doi.org/10.3390/jcm13113200

Chicago/Turabian Style

Bondar, Galyna, Abhinandan Das Mahapatra, Tra-Mi Bao, Irina Silacheva, Adrian Hairapetian, Thomas Vu, Stephanie Su, Ananya Katappagari, Liana Galan, Joshua Chandran, and et al. 2024. "An Exercise Immune Fitness Test to Unravel Disease Mechanisms—A Proof-of-Concept Heart Failure Study" Journal of Clinical Medicine 13, no. 11: 3200. https://doi.org/10.3390/jcm13113200

APA Style

Bondar, G., Mahapatra, A. D., Bao, T. -M., Silacheva, I., Hairapetian, A., Vu, T., Su, S., Katappagari, A., Galan, L., Chandran, J., Adamov, R., Mancusi, L., Lai, I., Rahman, A., Grogan, T., Hsu, J. J., Cappelletti, M., Ping, P., Elashoff, D., ... Deng, M. C. (2024). An Exercise Immune Fitness Test to Unravel Disease Mechanisms—A Proof-of-Concept Heart Failure Study. Journal of Clinical Medicine, 13(11), 3200. https://doi.org/10.3390/jcm13113200

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

Article Metrics

Back to TopTop