Next Article in Journal
Friend or Foe? Exploring the Role of Cytomegalovirus (HCMV) Infection in Head and Neck Tumors
Previous Article in Journal
Mechanical and Material Analysis of 3D-Printed Temporary Materials for Implant Reconstructions—A Pilot Study
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Biomedicines
Volume 12
Issue 4
10.3390/biomedicines12040871
biomedicines-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Decoding Immuno-Competence: A Novel Analysis of Complete Blood Cell Count Data in COVID-19 Outcomes

by
Prakasha Kempaiah
1,†,
Claudia R. Libertin
2,†,
Rohit A. Chitale
3,
Islam Naeyma
4,
Vasili Pleqi
1,
Johnathan M. Sheele
5,
Michelle J. Iandiorio
6,
Almira L. Hoogesteijn
7,
Thomas R. Caulfield
4,8 and
Ariel L. Rivas
9,*
1
Department of Medicine, Division of Infectious Diseases, Mayo Clinic, Jacksonville, FL 32224, USA
2
Department of Critical Care, Mayo Clinic, Jacksonville, FL 32224, USA
3
Mayo Clinic Alix School of Medicine, Mayo Clinic, Jacksonville, FL 32224, USA
4
Department of Neuroscience, Division of QHS Computational Biology, Mayo Clinic, Jacksonville, FL 32224, USA
5
Department of Emergency Medicine, Mayo Clinic, Jacksonville, FL 32224, USA
6
Department of Internal Medicine, School of Medicine, University of New Mexico, Albuquerque, NM 87131, USA
7
Human Ecology, Centro de Investigaciones Avanzadas, Merida 97310, Mexico
8
Department of Biochemistry and Molecular Biology, Mayo Clinic, Rochester, MN 55905, USA
9
Center for Global Health, Department of Internal Medicine, School of Medicine, University of New Mexico, Albuquerque, NM 87131, USA
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Biomedicines 2024, 12(4), 871; https://doi.org/10.3390/biomedicines12040871
Submission received: 10 February 2024 / Revised: 14 March 2024 / Accepted: 23 March 2024 / Published: 15 April 2024
(This article belongs to the Section Immunology and Immunotherapy)
Browse Figures
Review Reports Versions Notes

Abstract

:
Background: While ‘immuno-competence’ is a well-known term, it lacks an operational definition. To address this omission, this study explored whether the temporal and structured data of the complete blood cell count (CBC) can rapidly estimate immuno-competence. To this end, one or more ratios that included data on all monocytes, lymphocytes and neutrophils were investigated. Materials and methods: Longitudinal CBC data collected from 101 COVID-19 patients (291 observations) were analyzed. Dynamics were estimated with several approaches, which included non-structured (the classic CBC format) and structured data. Structured data were assessed as complex ratios that capture multicellular interactions among leukocytes. In comparing survivors with non-survivors, the hypothesis that immuno-competence may exhibit feedback-like (oscillatory or cyclic) responses was tested. Results: While non-structured data did not distinguish survivors from non-survivors, structured data revealed immunological and statistical differences between outcomes: while survivors exhibited oscillatory data patterns, non-survivors did not. In survivors, many variables (including IL-6, hemoglobin and several complex indicators) showed values above or below the levels observed on day 1 of the hospitalization period, displaying L-shaped data distributions (positive kurtosis). In contrast, non-survivors did not exhibit kurtosis. Three immunologically defined data subsets included only survivors. Because information was based on visual patterns generated in real time, this method can, potentially, provide information rapidly. Discussion: The hypothesis that immuno-competence expresses feedback-like loops when immunological data are structured was not rejected. This function seemed to be impaired in immuno-suppressed individuals. While this method rapidly informs, it is only a guide that, to be confirmed, requires additional tests. Despite this limitation, the fact that three protective (survival-associated) immunological data subsets were observed since day 1 supports many clinical decisions, including the early and personalized prognosis and identification of targets that immunomodulatory therapies could pursue. Because it extracts more information from the same data, structured data may replace the century-old format of the CBC.

1. Introduction

To prevent the spread of rapidly disseminating infections, such as COVID-19, personalized assessments of immunity have been proposed [1]. Yet, immuno-competence lacks a standardized method that could routinely be used in clinical medicine [2]. Tests that measure immuno-competence should be both informative and biologically valid. Ideally, they should also predict disease outcomes and help select personalized therapies [3,4]. To develop such methods, several biomedical and methodological concepts need to be considered and integrated into the chosen operationalizations, including (i) dynamics, (ii) non-binary (polychotomous) conditions, (iii) the risk of confounding, (iv) personalized factors, (v) biological complexity, (vi) data structuring and (vii) data analysis.

1.1. Dynamic, Polychotomous and Personalized Methods That Prevent Confounding

The dynamics of the immune responses, such as COVID-19, are poorly known [5]. One reason for such a cognitive gap is that, typically, classic methods have been static. Static methods assume only two conditions: “disease-positive” and “disease-negative” results. Binary and static methods are not valid because biological functions are dynamic processes that can reveal three or more biomedical conditions [6,7]. For example, inflammatory (immunological) responses exhibit three or more presentations: no inflammation, early inflammation, late inflammation (followed by recovery) and/or late inflammation (followed by chronic disease). Therefore, to assess immuno-competence, new methods should handle non-dichotomous and dynamic conditions [8,9,10]. By subtracting earlier from later data values, the potential bias of extremely low (or high) values of earlier tests can be prevented, and dynamics can be assessed [11].
In addition, assessments of immuno-competence should be personalized. Because patients differ in co-morbidities and medical histories, to prevent confounding, methods that attempt to capture dynamic and non-binary responses should not assume that populations are homogenous [12].
Because most published research on major diseases, such as COVID-19, has been cross-sectional and observational, the risk of bias is not trivial [13,14,15,16]. One source of confounding refers to the fact that patients may be hospitalized at different disease stages, and therefore, early and late disease stages could erroneously be grouped together. Another variation of the same error can occur when data with different disease stages are aggregated. Aggregate data do not provide personalized information but population averages. Thus, methods that handle non-binary, dynamic and personalized data are needed [17,18].
A third common source of confounding emerges when data distributions of different biological conditions—such as “cases” and “non-cases”—overlap [19,20]. Data overlapping (different outcomes associated with similar data intervals) should be avoided because it prevents medical discrimination even when statistical significance is achieved [21,22]. The previous considerations result in one consequence that novel methods should also address: given the dynamic, non-binary and personalized aspects of disease processes, there is no standard patient, and therefore, there is no control either [23].

1.2. Complex and Dynamic Interactions That Structure the Data and Show Patterns

To measure personalized immune dynamics and to extract hidden patterns, the assessment of immunological complexity has been recommended [24,25,26]. To explore complexity, interactions involving two or more biological elements—such as lymphocytes and monocytes—may describe biological functions that influence outcomes [27]. Such complex functions change over time and may exhibit feedforward and feedback loops [9,28].
Structured data—not isolated variables, such as the counts or percentages of a single cell type—can investigate relationships. Data structuring may provide interpretable and medically useful information and, in addition, avoid the confounding associated with unstructured data [29].
Structured data may promote pattern recognition. The shapes of data distributions may inform [30,31]. One specific data distribution (positive kurtosis or L-shaped data) has been viewed, for many years, as statistically intractable [32]. Yet, such data presentations are not rare and may be informative in several medical fields [33,34,35].
Interpretations also depend on the format used to present and analyze the data. Because complex and dynamic processes involve interactions, ratios—which assess relationships among two or more elements—are preferred over counts. Novel methods that utilize the data generated by the century-old CBC could inform more or better if they are not limited to analyses of counts [29,36].

1.3. COVID-19-Related Dynamic and Personalized Immuno-Competence

The previous considerations may influence the way immuno-competence is explored. In several clinical entities (e.g., COVID-19, sepsis, hantavirus), the early disease stage is not always associated with death [27,37,38,39,40].
Immuno-competence cannot be evaluated with a single (static) assessment. This concept may matter particularly in sepsis where dynamic changes can occur and both excessive inflammation and immuno-suppression may induce fatalities [41,42]. Therefore, to estimate immuno-competence, temporal assessments are needed.
In adopting a personalized approach, here, longitudinal data on COVID-19 are retrospectively explored with the purpose of assessing immunological dynamics and their associated outcomes. Assuming that survival is an unambiguous outcome of immuno-competent responses, this study was set to answer whether, soon after admission, (i) structured and temporal blood-related data can extract more information from the same data than unstructured data and (ii) whether non-overlapping and/or distinct data distributions may, rapidly and easily, predict disease outcomes in COVID-19 patients.

2. Materials and Methods

2.1. Clinical Data

Two-hundred and ninety-one complete blood count (CBC/differential) samples collected from 101 anonymized individuals hospitalized at the Mayo Clinic of Jacksonville, Florida before August of 2020 with a COVID-19 diagnosis are reported in Table S1. Inclusion criteria considered subjects of 18 or more years of age, treated at or before August 2020, with SARS-CoV-2 positive test results conducted within 72 h of admission, radiographic changes consistent with COVID-19 and deemed to be at risk of severe illness. Discharge was based on negative testing and the alleviation of life-threatening conditions. Individuals treated for immunosuppression, malignancy, pregnancy and/or hospitalized for at least three weeks in the previous six months were excluded [27]. The outcome (survival or no survival) was determined at 30 days after hospitalization.

2.2. Data Structure and Analysis

Time was measured with seven approaches that analyzed uni-, bi- and three-dimensional data. The first approach analyzed the data collected on the first post-hospitalization day (Day 1). The second approach considered the data collected at days 1 and 2. The third approach evaluated the difference between day 2 and day 1 observations. These analyses were uni-dimensional and focused on a single (unstructured or structured) variable.
The fourth and fifth approaches were bi-dimensional, considering two (unstructured or structured) variables. The sixth and seventh approaches were three-dimensional, focusing on personalized and therapy-oriented information, respectively.
In all approaches, the presence or absence of L-shaped data distributions (positive kurtosis) was determined. Kurtosis (non-normally distributed [not bell-shaped] data) was documented when observations differed 5 or more standard deviations from the mean, and consequently, most observations were concentrated on one end of the distribution (hence, the shape of the data resembled an ‘L’). The statistical analysis was limited to indicating whether outcomes revealed non-kurtotic data intervals, and, when kurtosis was observed, whether L-shaped data distributions were revealed by one or many blood-related variables.
In addition, the informative effect of non-structured and structured data was also investigated. Observations were analyzed as either separate counts or relative percentages of single-cell types (unstructured approach) or complex multicellular interactions that included two or more ratios (structured approach). The structured method has been described before when several diseases that affect human and non-human species were investigated [21,24,27,28,40]. The structured method is based on dimensionless indicators (DIs)—temporary guides used to detect distinct data subsets. DIs include numerous relationships, and they are complex: they have two or more interactions that, individually, could include two or more elements, e.g., the triple ratio resulting from (i) generating a ratio between neutrophils and monocytes (the N% divided the M% or N/M ratio), (ii) creating a separate ratio between mononuclear cells (L%+M% or MC cells) and neutrophils (or MC/N ratio), and (iii) the overall relationship resulting from the first divided the second ratio. Because DLs are used to display three-dimensional patterns, the number of possible combinations observed is very high (if not infinite) at least because each axis may influence and/or be influenced by any other axis. Given such internal and external possible variability, DLs are identified with letters lacking any biological meaning (such as AA, AAA, BBB…). Biologically interpretable contents are explicitly investigated only after data patterns are detected in 3D space.
The rationale of this method can be found in Niels Jerne’s work [43]. He emphasized the language-like properties of immunology, i.e., a set of few or very few elements (comparable to the 26 letters of the English alphabet) which, once combined and recombined, results not only in a very large number of ‘words’, but also in an even larger number of ‘sentences’ and ‘paragraphs’. When we refer to Biology, this combinatorial ‘language’ can generate millions of antibodies and/or proteins that can recognize antigens even when the immune system has not previously encountered such antigens.
This means that analyzing only the elemental immunological units is not informative: in considering the language metaphor, no knowledge can be derived from memorizing all the letters of the alphabet. Similarly, useful biomedical knowledge may be missed, unless and until basic leukocyte data are (a) reformatted as information, (b) information is converted into knowledge, and (c) finally, knowledge is applied as concrete (context- and resource-dependent) decisions. This process has been described as the ‘DIKW (data–information–knowledge–‘wisdom’) pyramid’ [29].
While a proprietary software can, in principle, conduct all these steps in real time, any method that focuses on pattern recognition and explores three-dimensional data combinations (structures in which the information emerges only after the three variables converge in space) is likely to detect usable information, as previously described [24].
This method was conducted with a proprietary software package [29,39,44]. To generate visual patterns and to assess statistical properties, a commercial software package was used (Minitab 21, Minitab LLC, State College, PA, USA).

3. Results

Due to overlapping data intervals of different outcomes, no medical discrimination was possible, at any time when either counts or percentages of cell types were assessed in isolation (rectangles, Figure 1A–D). Overlapping data distributions occurred even when statistically significant differences were observed. For example, when all the temporal data were considered, the median lymphocyte percentage of non-survivors was 12.3% (n = 37), while that of survivors was 17.2% (n = 254), a difference that reached statistical significance (p < 0.01, Mann–Whitney test, Table S2). This statistical difference occurred even when most survivor and non-survivor observations overlapped (Figure 1D).
The lack of discrimination between outcomes was not limited to leukocyte-related data. Other hematological variables—such as IL-6 and hemoglobin—also showed, over time, overlapping distributions between outcomes (Figure 2A,B). In contrast, structured data that simultaneously captured many interactions among all leukocytes displayed not only statistically significant differences between outcomes, but also a range of observations with non-overlapping outcomes (Figure 2C).
When, regardless of outcome, leukocytes were assessed as separate (non-structured) cell types at hospitalization day 1, no L-shaped (non-kurtotic) intervals were observed (Table S2, Figure 3A,B). However, when other hematological variables—such as IL-6 and Hb—were assessed, survivors revealed kurtosis, while non-survivors did not (Table S2, Figure 3C–F). Similarly, structured leukocyte data showed kurtosis, although only in survivors (AAT, Figure 3G,H).
Regardless of outcome, kurtosis was not observed when all longitudinal observations were assessed as non-structured data, i.e., when the percentages of lymphocytes, neutrophils or monocytes were investigated (Table S2, Figure 4A–F). In contrast, both non-leukocyte and structured leukocyte data displayed kurtosis when survivors were tested and all temporal data points were evaluated (Figure 5A–F). Hence, it was concluded that kurtosis was associated with both (i) the structure of the data (complex ratios) and (ii) the outcome (survival).
To prevent biased assessments of dynamics, an additional comparison was conducted, which considered the net difference between two consecutive observations. Dynamics were then evaluated in two ways: (a) considering the net difference of values shown between day 2 and day 1 by either survivors or non-survivors and (b) considering the dispersion of such values when survivors and non-survivors were measured simultaneously. When day 1 values were subtracted from day 2 values, survivors displayed kurtosis, while non-survivors did not (Table S2, Figure 6A–F). Some variables—such as IL-6 and complex indicators that assessed leukocyte multicellular interactions—oscillated over time in competent immunological responses. In contrast, the result of the subtraction, for most survivors, approached zero (Figure 6B,F).
The last comparison simultaneously assessed both outcomes. Complementing the assessments of kurtosis reported in Figure 6, dynamic assessments supported the hypothesis that survivors displayed numerous and large temporal oscillations around zero, while non-survivors did not (Figure 7A–F).
To further explore the hypothesis that survivors and non-survivors differed immunologically, additional time-related, bi-dimensional comparisons between pairs of cell types were conducted with non-structured data. While survivors exhibited a broad range of observations that could facilitate relationships between neutrophils and monocytes at days 1 and 2, non-survivors lacked such data ranges (Figure 8A–D). Likewise, non-survivors did not show high values of lymphocyte percentages that could interact with monocytes at such times, while survivors did (Figure 8E–H).
More information was extracted from the same data when bi-dimensional, temporal and structured data were investigated. Two data ranges occupied only by survivors were revealed when complex leukocyte interactions were investigated at days 1 and 2 (Figure 9A–D).
The two protective data ranges of immunological interactions were not influenced by time (Figure 10A–C). The two immunological relationships associated with survival (here named ‘A’ and ‘B’) were characterized by non-overlapping intervals of non-structured data (the percentages of lymphocytes [L], neutrophils [N] or monocytes [M]) and structured data (the small leukocytes [L and N]/M, the N/M and the N/L ratios, Figure 10D,E).
When all temporal observations were considered (n = 291), three-dimensional assessments of structured data retrieved even more information from the data: three data subsets—perpendicular to one another—included only survivors when the first three hospitalization days were assessed (Figure 11A,B). Therefore, immuno-competence (responses that promoted survival) was an early process—it was measurable in the first three post-admission days.
In classifying the data into four groups (three groups composed of only survivors [here named I, II and III] and one group that included survivors and non-survivors, Figure 11C), survivors revealed three non-overlapping data intervals. The data groups that included only survivors differed from one another: group I showed the lowest values of lymphocyte percentages and highest neutrophil/lymphocyte (N/L) ratios; group II showed the lowest monocyte values intermediate and N/L ratio values between those of groups I and III; and group III displayed the lowest neutrophil percentages (Figure 11D).
When day 1 and day 2 structured data were investigated (n = 195) as recorded (i.e., when the same patient contributed two data points), the findings supported at least four types of inferences. For example, between one and three data ranges concerning only a single outcome were detected at days 1 and 2 (Figure 12A,B). Because survival-related findings were based on two different data structures that considered two temporal points (a double redundant approach), inferences revealed internal validity (Figure 12A,B).
The three data subsets containing only survivors were named ‘left’, ‘top’ and ‘right’ (Figure 12C). They differed immunologically from one another: the same non-structured and structured variables that distinguished similar subsets when all temporal data were analyzed (reported in Figure 11) also differentiated the subsets identified with the data collected on the first two hospitalization days (Figure 12D).
In addition, day 1 and day 2 structured data displayed trajectory: the directionality of the temporal data (where each data point came from/went to) was objectively displayed using an easily detectable structure (a single, one data point-wide, line of observations that simultaneously captures time and immunological interactions, Figure 12D). In this example, the arrows generated by this minimal (only 24 h long) temporal assessment identified two prognostic inferences: two movements away from zero predicted survival at a personalized level, i.e., even a single individual could be prognosticated.
Therefore, the findings reveal that (i) assessments limited to day 1 data are adequate for predicting outcomes; (ii) the data collected at day 1 and day 2 revealed trajectory; (iii) the assessment of dynamics (day 2 minus day 1 data) demonstrated that survivors’ data oscillated, while the data of non-survivors did not; and (iv) numerous hematological, non-cellular variables also revealed kurtosis-related differences between non-survivors and survivors.

4. Discussion

4.1. Caveats

The fundamental limitation of this study refers to the fact that there is no consensus on the meaning of immuno-competence. To address this challenge, this study analyzed both humoral and cell-mediated responses. Because both approaches generated similar inferences (survivors exhibited L-shaped [kurtotic] data intervals while non-survivors did not), the stated limitation was ameliorated.
A second potential limitation is that, given the combinatorial nature of leukocyte-related multicellular interactions, the number of possible ratios to be assessed is very high, and consequently, any investigated ratio is not necessarily representative. The analysis of two or more leukocyte-derived complex indicators (redundant analysis) diminished this risk.

4.2. Seven Major Findings

The data supported seven inferences: (1) over time, survivors displayed non-randomly distributed (circular-like) immunological responses that revealed feedback/feedforward loops; (2) L-shaped data distributions (kurtosis) characterized survivors; (3) non-survivors were associated with static immunological responses; (4) immuno-competence (responses that, over time, differed in magnitude and/or internal composition) was expressed early, after hospitalization; (5) data ranges of immunological data associated with survival may provide targets for immunomodulatory strategies; (6) unstructured, simple (reductionist) indicators—such as counts and percentages—missed information conveyed by complex (non-reductionist) indicators; and (7) the anticipatory testing of immuno-competence (evaluations conducted in the absence of a disease) may be desirable. These inferences indicate that the time-related detection of non-randomly distributed data intervals may offer medically useful information.
These findings corroborate the non-random (oscillatory) data features of blood leukocytes. The circadian cycle of blood leukocytes has been described before—a process associated with immuno-competence [45]. For example, the circadian cycles of CD4 and CD8+ cells are opposed to one another [46]. The circadian cycle that regulates the hypothalamic–pituitary–thyroid axis is altered in immuno-suppression: in lung cancer, the ratio of melatonin/cortisol is diminished (cortisol is increased), and the count and relative proportions of CD8+ cells are decreased [47]. Such feedback-like cycles are not randomly distributed. For instance, lymphocyte data exhibit non-normally distributed intervals [48,49]. Because effective (immuno-competent) responses are non-random [50,51], a lack of such dynamics may reflect immuno-suppression or deficiency.
Here, we followed the new recommended research emphasis that promotes explorations of L-shaped data distributions, i.e., kurtosis [52,53]. In L-shaped distributions, data oscillations around zero (a quasi-stationary mode) induce a left-side peak [54]. Chronic infections have been characterized by a quasi-stationary mode in which lymphocytes are inactive [55]. Because non-normal (L-shaped or kurtotic) data distributions were mainly exhibited by survivors, the hypothesis that kurtosis describes cyclic immunological oscillations associated with survival was not rejected [9,24]. This continuous ‘fine-tuning’ resembles the way cars are driven: to drive straight, an open-ended series of minor turns to the left and right are needed. Other biological (non-immunological) functions also show similar cyclic adjustments [56].
Numerous studies have suggested that immuno-suppression is associated with COVID-19 [57,58,59,60,61,62,63]. This study also corroborated earlier reports that indicated immune-competence, in COVID-19, predicts survival [3]. The fact that, in patients with cancer, COVID-19 is associated with increased mortality supports the view that immuno-suppression is a negative predictor of disease outcome [64,65].
The findings also support the hypothesis that immuno-competence, in COVID-19, involves the cellular immune response [66,67,68,69]. Patients presenting with severe COVID-19 have frequently been reported with lymphopenia [70,71].
This study also showed differences between survivors and non-survivors when non-cellular factors (including hemoglobin, IL-6 and hematocrit) were analyzed [72]. Because numerous humoral indicators differed between the survivors and non-survivors and such differentiation was robust over time, the involvement of both arms of the immune system on disease outcomes was documented, and the approach used to explore immuno-competence exhibited robust internal and external validity [73].
Non-survivors showed differences in their immune responses as early as hospitalization day 1. Because the kurtosis associated with survivors was not detected after day 3, the findings support the hypothesis in that immuno-competence should be evaluated early [74]. Because immuno-suppression was already detected at hospitalization day 1, the hypothesis that, in non-survivors, immunosuppression could predate the infection cannot be ruled out [3]. This hypothesis has been sustained before: patients tested two or more weeks before SARS-CoV-2 infections took place have shown reduced numbers of eosinophils—a finding associated with poorer outcomes [75]. Very early immunological differences between survivors and non-survivors have also been reported in other studies [76,77,78].
Data ranges of immunological data associated with survival may provide targets for immunomodulatory therapies. Because, to be effective, immunomodulation should be implemented early, the findings provide unambiguous targets for such therapies [79].
In agreement with earlier reports, unstructured, simple indicators—such as counts and percentages—revealed weaknesses associated with reductionist approaches [21,80]. In contrast, structured (complex, non-reductionist) indicators provided biologically interpretable and medically useful information [24,26]. When complex data structures reveal distinct patterns—such as single, one data point-wide, lines of observations—the directionality of temporal data promotes earlier evaluations [81].
Because the risk of COVID-19-related death may be high even in vaccinated immuno-suppressed patients, the findings supported the implementation of the anticipatory testing of immuno-competence, i.e., before diseases are diagnosed [82]. This possibility may be synonymous with updating the format of the CBC/differential: after a century in which no biological functions (no interactions among cells and molecules that participate in immune responses) have been assessed, it is now possible to extract more biologically interpretable and medically useful information from the same data.

5. Conclusions

This study suggests that the temporal assessment of non-randomly distributed, oscillatory cycles of immune responses may facilitate early and personalized prognosis as well as identification of targets to be potentially used by immunomodulatory-oriented therapies. To further explore this proposition, additional studies on dynamics of host–microbial interactions are recommended.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/biomedicines12040871/s1: Table S1: Longitudinal assessments of 101 COVID-19 patients (n = 291); Table S2: Temporal assessments of L-shaped distributions.

Author Contributions

Conceptualization, T.R.C.; Methodology, M.J.I.; Validation, P.K.; Resources, C.R.L., R.A.C. and V.P.; Data curation, I.N.; Writing—original draft, A.L.R.; Writing—review & editing, J.M.S.; Visualization, A.L.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

This study was approved by the Review Board of Mayo Clinic, Jacksonville, FL, USA (protocol D:21-002778, 21 April 2021).

Informed Consent Statement

Patient consent was waived due to the retrospective and anonymous nature of the study, i.e., no patient was hospitalized at the time of data collection and all identifiers were removed.

Data Availability Statement

Data are contained within the article and Supplementary Materials.

Acknowledgments

The technical support provided by José Luis Febles (Dept. of Human Ecology, Centro de Investigaciones Avanzadas [CINVESTAV], Merida, Mexico) is appreciated.

Conflicts of Interest

A.L.R. and A.L.H. are co-inventors of the software utilized. All other authors declare no conflict of interest.

References

  1. Bland, J.S. COVID-19 Risk: Clinical Tools for Assessing and Personalizing Immunity. Integr. Med. 2021, 20, 18–23. [Google Scholar]
  2. Chang, S.; Kohrt, H.; Maecker, H.T. Monitoring the immune competence of cancer patients to predict outcome. Cancer Immunol. Immunother. 2014, 63, 713–719. [Google Scholar] [CrossRef] [PubMed]
  3. Lee, G.C.; Restrepo, M.I.; Harper, N.; Manoharan, M.S.; Smith, A.M.; Meunier, J.A.; Sanchez-Reilly, S.; Ehsan, A.; Branum, A.P.; Winter, C.; et al. Immunologic resilience and COVID-19 survival advantage. J. Allergy Clin. Immunol. 2021, 148, 1176–1191. [Google Scholar] [CrossRef] [PubMed]
  4. Grimm, C.; Dickel, S.; Grundmann, J.; Payen, D.; Schanz, J.; Zautner, A.E.; Tampe, B.; Moerer, O.; Winkler, M.S. Case Report: Interferon-g Rescues Monocytic Human Leukocyte Antigen Receptor (mHLA-DR) Function in a COVID-19 Patient With ARDS and Superinfection With Multiple MDR 4MRGN Bacterial Strains. Front. Immunol. 2021, 12, 753849. [Google Scholar] [CrossRef]
  5. Li, Q.; Cao, Y.; Chen, L.; Wu, D.; Yu, J.; Wang, H.; He, W.; Chen, L.; Dongv, F.; Chen, W.; et al. Hematological features of persons with COVID-19. Leukemia 2020, 34, 2163–2172. [Google Scholar] [CrossRef] [PubMed]
  6. Grossman, Z.; Berke, G. Tumor escape from immune elimination. J. Theor. Biol. 1980, 83, 267–296. [Google Scholar] [CrossRef] [PubMed]
  7. Sontag, E.D. A Dynamic Model of Immune Responses to Antigen Presentation Predicts Different Regions of Tumor or Pathogen Elimination. Cell Syst. 2017, 22, 231–241.e11. [Google Scholar] [CrossRef] [PubMed]
  8. Assmus, H.E.; Herwig, R.; Cho, K.-H.; Wolkenhauer, O. Dynamics of biological systems: Role of systems biology in medical research. Expert Rev. Mol. Diagn. 2006, 6, 891–902. [Google Scholar] [CrossRef] [PubMed]
  9. Lesterhuis, W.J.; Bosco, A.; Millward, M.J.; Small, M.; Nowak, A.K.; Lake, R.A. Dynamic versus static biomarkers in cancer immune checkpoint blockade: Unravelling complexity. Nat. Rev. Drug Discov. 2017, 16, 264–272. [Google Scholar] [CrossRef] [PubMed]
  10. Giannoula, A.; Gutierrez-Sacristán, A.; Bravo, Á.; Sanz, F.; Furlong, L.I. Identifying temporal patterns in patient disease trajectories using dynamic time warping: A population-based study. Sci. Rep. 2018, 8, 4216. [Google Scholar] [CrossRef] [PubMed]
  11. Taber, D.J.; Palanisamy, A.P.; Srinivas, T.R.; Titte, R.; Gebregziabher, M.; Odeghe, J.; Chavin, K.D.; Egede, L.E.; Baliga, P.K. Inclusion of dynamic clinical data improves the predictive performance of a 30-day readmission risk model in kidney transplantation. Transplantation 2015, 99, 324–330. [Google Scholar] [CrossRef] [PubMed]
  12. Christaki, E.; Giamarellos-Bourboulis, E.J. The beginning of personalized medicine in sepsis: Small steps to a bright future. Clin. Genet. 2014, 86, 56–61. [Google Scholar] [CrossRef] [PubMed]
  13. Tleyjeh, I.M.; Kashour, T.; Mandrekar, J.; Petitti, D.B. Overlooked Shortcomings of Observational Studies of Interventions in Coronavirus Disease 2019: An Illustrated Review for the Clinician. Open Forum Infect. Dis. 2021, 8, ofab317. [Google Scholar] [CrossRef] [PubMed]
  14. Millard, L.A.C.; Fernández-Sanlés, A.; Carter, A.R.; Hughes, R.A.; Tilling, K.; Morris, T.P.; Major-Smith, D.; Griffith, G.J.; Clayton, G.L.; Kawabata, E.; et al. Exploring the impact of selection bias in observational studies of COVID-19: A simulation study. Int. J. Epidemiol. 2023, 52, 44–57. [Google Scholar] [CrossRef] [PubMed]
  15. Fung, K.; Jones, M.; Doshi, P. Sources of bias in observational studies of COVID-19 vaccine effectiveness. J. Eval. Clin. Pract. 2024, 30, 30–36. [Google Scholar] [CrossRef] [PubMed]
  16. Accorsi, E.K.; Qiu, X.; Rumpler, E.; Kennedy-Shaffer, L.; Kahn, R.; Joshi, K.; Goldstein, E.; Stensrud, M.J.; Niehus, R.; Cevik, M.; et al. How to detect and reduce potential sources of biases in studies of SARS-CoV-2 and COVID-19. Eur. J. Epidemiol. 2021, 36, 179–196. [Google Scholar] [CrossRef] [PubMed]
  17. Huang, Y.; Tang, J.; Tam, W.W.; Mao, C.; Yuan, J.; Di, M.; Yang, Z. Comparing the Overall Result and Interaction in Aggregate Data Meta-Analysis and Individual Patient Data Meta-Analysis. Medicine 2016, 95, e3312. [Google Scholar] [CrossRef] [PubMed]
  18. Lyman, G.H.; Kuderer, N.M. The strengths and limitations of meta-analyses based on aggregate data. BMC Med. Res. Methodol. 2005, 5, 14. [Google Scholar] [CrossRef] [PubMed]
  19. Mueller, Y.M.; Schrama, T.J.; Ruijten, R.; Schreurs, M.W.; Grashof, D.G.; van de Werken, H.J.; Lasinio, G.J.; Álvarez-Sierra, D.; Kiernan, C.H.; Castro Eiro, M.D.; et al. Stratification of hospitalized COVID-19 patients into clinical severity progression groups by immuno-phenotyping and machine learning. Nat. Commun. 2022, 13, 915. [Google Scholar] [CrossRef] [PubMed]
  20. Altman, D.G.; Royston, P. What do we mean by validating a prognostic model? Statist. Med. 2000, 19, 453–473. [Google Scholar] [CrossRef]
  21. Leitner, G.; Blum, S.; Rivas, A.L. Visualizing the indefinable: Three-dimensional complexity of ‘infectious diseases. PLoS ONE 2015, 10, e0123674. [Google Scholar] [CrossRef] [PubMed]
  22. Lo, A.; Chernoff, H.; Zheng, T.; Lo, S.H. Why significant variables aren’t automatically good predictors. Proc. Natl. Acad. Sci. USA 2015, 112, 13892–13897. [Google Scholar] [CrossRef] [PubMed]
  23. Anjum, R.L.; Copeland, S.; Rocca, E. Rethinking Causality and Evidence for the Unique Patient; Springer: Berlin/Heidelberg, Germany, 2020; Available online: http://library.oapen.org/handle/20.500.12657/39574 (accessed on 3 March 2024).
  24. Rivas, A.L.; Leitner, G.; Jankowski, M.D.; Hoogesteijn, A.L.; Iandiorio, M.J.; Chatzipanagiotou, S.; Ioannidis, A.; Blum, S.E.; Piccinini, R.; Antoniades, A.; et al. Nature and consequences of biological reductionism for the immunological study of infectious diseases. Front. Immunol. 2017, 8, 612. [Google Scholar] [CrossRef]
  25. Agur, Z.; Elishmereni, M.; Foryś, U.; Kogan, Y. Accelerating the Development of Personalized Cancer Immunotherapy by Integrating Molecular Patients’ Profiles with Dynamic Mathematical Models. Clin. Pharmacol. Ther. 2020, 108, 515–527. [Google Scholar] [CrossRef] [PubMed]
  26. Relecom, A.; Merhi, M.; Inchakalody, V.; Uddin, S.; Rinchai, D.; Bedognetti, D.; Dermime, S. Emerging dynamics pathways of response and resistance to PD-1 and CTLA-4 blockade: Tackling uncertainty by confronting complexity. J. Exp. Clin. Cancer Res. 2021, 40, 74. [Google Scholar] [CrossRef] [PubMed]
  27. Verma, J.S.; Libertin, C.R.; Gupta, Y.; Khanna, G.; Kumar, R.; Arora, B.S.; Krishna, L.; Fasina, F.O.; Hittner, J.B.; Antoniades, A.; et al. Multi-Cellular Immunological Interactions Associated With COVID-19 Infections. Front. Immunol. 2022, 13, 794006. [Google Scholar] [CrossRef]
  28. Rivas, A.L.; Jankowski, M.D.; Piccinini, R.; Leitner, G.; Schwarz, D.; Anderson, K.L.; Fair, J.M.; Hoogesteijn, A.L.; Wolter, W.; Chaffer, M.; et al. Feedback-based, system-level properties of vertebrate-microbial interactions. PLoS ONE 2013, 8, e53984. [Google Scholar] [CrossRef] [PubMed]
  29. Libertin, C.R.; Kempaiah, P.; Gupta, Y.; Fair, J.M.; van Regenmortel, M.H.V.; Antoniades, A.; Rivas, A.L.; Hoogesteijn, A.L. Data structuring may prevent ambiguity and improve personalized medical prognosis. Mol. Aspects Med. 2022, 15, 101142. [Google Scholar] [CrossRef] [PubMed]
  30. Greiff, V.; Bhat, P.; Cook, S.C.; Menzel, U.; Kang, W.; Reddy, S.T. A bioinformatic framework for immune repertoire diversity profiling enables detection of immunological status. Genome Med. 2015, 7, 49. [Google Scholar] [CrossRef]
  31. Miho, E.; Yermanos, A.; Weber, C.R.; Berger, C.T.; Reddy, S.T.; Greiff, V. Computational Strategies for Dissecting the High-Dimensional Complexity of Adaptive Immune Repertoires. Front. Immunol. 2018, 9, 224. [Google Scholar] [CrossRef] [PubMed]
  32. Bradley, J.V. The insidious L-shaped distribution. Bull. Psychon. Soc. 1982, 20, 85–88. [Google Scholar] [CrossRef]
  33. Raab, P.; Hattingen, E.; Franz, K.; Zanella, F.E.; Lanfermann, H. Cerebral gliomas: Diffusional kurtosis imaging analysis of microstructural differences. Radiology 2010, 254, 876–881. [Google Scholar] [CrossRef]
  34. Hittner, J.B.; Fasina, F.O.; Hoogesteijn, A.L.; Piccinini, R.; Maciorowski, D.; Kempaiah, P.; Smith, S.D.; Rivas, A.L. Testing-related and geo-demographic indicators strongly predict COVID-19 deaths in the United States during March of 2020. Biomed. Environ. Sci. 2021, 34, 734–738. [Google Scholar] [CrossRef]
  35. Van Cauter, S.; Veraart, J.; Sijbers, J.; Peeters, R.R.; Himmelreich, U.; De Keyzer, F.; Van Gool, S.W.; Van Calenbergh, F.; De Vleeschouwer, S.; Van Hecke, W.; et al. Gliomas: Diffusion kurtosis MR imaging in grading. Radiology 2012, 263, 492–501. [Google Scholar] [CrossRef] [PubMed]
  36. Houwen, B. The Differential Cell Count. Lab. Hematol. 2001, 7, 89–100. [Google Scholar]
  37. Hamers, L.; Kox, M.; Pickkers, P. Sepsis-induced immunoparalysis: Mechanisms, markers, and treatment options. Minerva Anestesiol. 2015, 81, 426–439. [Google Scholar] [PubMed]
  38. Peters van Ton, A.M.; Kox, M.; Abdo, W.F.; Pickkers, P. Precision Immunotherapy for Sepsis. Front. Immunol. 2018, 9, 1926. [Google Scholar] [CrossRef] [PubMed]
  39. Rivas, A.L.; Hoogesteijn, A.L.; Antoniades, A.; Tomazou, M.; Buranda, T.; Perkins, D.J.; Fair, J.M.; Durvasula, R.; Fasina, F.O.; Tegos, G.P.; et al. Assessing the dynamics and complexity of disease pathogenicity using 4-dimensional immunological data. Front. Immunol. 2019, 10, 1258. [Google Scholar] [CrossRef] [PubMed]
  40. Chatzipanagiotou, S.; Ioannidis, A.; Trikka-Graphakos, E.; Charalampaki, N.; Sereti, C.; Piccinini, R.; Higgins, A.M.; Buranda, T.; Durvasula, R.; Hoogesteijn, A.L.; et al. Detecting the hidden properties of immunological data and predicting the mortality risks of infectious syndromes. Front. Immunol. 2016, 7, 217. [Google Scholar] [CrossRef] [PubMed]
  41. Venet, F.; Monneret, G. Advances in the understanding and treatment of sepsis-induced immunosuppression. Nat. Rev. Nephrol. 2018, 14, 121–137. [Google Scholar] [CrossRef] [PubMed]
  42. Pradhan, K.; Yi, Z.; Geng, S.; Li, L. Development of exhausted memory monocytes and underlying mechanisms. Front. Immunol. 2021, 12, 778830. [Google Scholar] [CrossRef] [PubMed]
  43. Jerne, N.K. The generative grammar of the immune system. EMBO J. 1985, 4, 847–852. [Google Scholar] [CrossRef] [PubMed]
  44. Satashia, P.H.; Franco, P.M.; Rivas, A.L.; Isha, S.; Hanson, A.; Narra, S.A.; Singh, K.; Jenkins, A.; Bhattacharyya, A.; Guru, P.; et al. From numbers to medical knowledge: Harnessing combinatorial data patterns to predict COVID-19 resource needs and distinguish patient subsets. Front. Med. 2023, 10, 240426. [Google Scholar] [CrossRef] [PubMed]
  45. Mazzoccoli, G.; Sothern, R.B.; Greco, A.; Pazienza, V.; Vinciguerra, M.; Liu, S.; Cai, Y. Time-related dynamics of variation in core clock gene expression levels in tissues relevant to the immune system. Int. J. Immunopathol. Pharmacol. 2011, 24, 869–879. [Google Scholar] [CrossRef] [PubMed]
  46. Mazzoccoli, G.; De Cata, A.; Greco, A.; Damato, M.; De Pinto, D.G.; Montrano, M.; Marzulli, N.; Dagostino, M.P.; Carughi, S. Opposing circadian rhythms of CD3+, CD4+ and CD3+, CD8+ lymphocyte subpopulations in healthy humans. Biol. Rhythm. Res. 2011, 42, 111–118. [Google Scholar] [CrossRef]
  47. Mazzoccoli, G.; Vendemiale, G.; De Cata, A.; Carughi, S.; Tarquini, R. Altered time structure of neuro-endocrine-immune system function in lung cancer patients. BMC Cancer 2010, 10, 314. [Google Scholar] [CrossRef]
  48. Fricke, G.M.; Letendre, K.A.; Moses, M.E.; Cannon, J.L. Persistence and Adaptation in Immunity: T Cells Balance the Extent and Thoroughness of Search. PLoS Comput. Biol. 2016, 12, e1004818. [Google Scholar] [CrossRef] [PubMed]
  49. Macfarlane, F.R.; Lorenzi, T.; Chaplain, M.A.J. Modelling the Immune Response to Cancer: An Individual-Based Approach Accounting for the Difference in Movement Between Inactive and Activated T Cells. Bull. Math. Biol. 2018, 80, 1539–1562. [Google Scholar] [CrossRef] [PubMed]
  50. Potthoff, R.F.; Tudor, G.E.; Pieper, K.S.; Hasselblad, V. Can one assess whether missing data are missing at random in medical studies? Stat. Methods Med. Res. 2006, 15, 213–234. [Google Scholar] [CrossRef]
  51. Lakshmikanth, T.; Muhammad, S.A.; Olin, A.; Chen, Y.; Mikes, J.; Fagerberg, L.; Gummesson, A.; Bergström, G.; Uhlen, M.; Brodin, P. Human Immune System Variation during 1 Year. Cell Rep. 2020, 32, 107923. [Google Scholar] [CrossRef] [PubMed]
  52. Cain, M.K.; Zhang, Z.; Yuan, K.H. Univariate and multivariate skewness and kurtosis for measuring nonnormality: Prevalence, influence and estimation. Behav. Res. Methods 2017, 49, 1716–1735. [Google Scholar] [CrossRef] [PubMed]
  53. Grossman, Z. Immunological Paradigms, Mechanisms, and Models: Conceptual Understanding Is a Prerequisite to Effective Modeling. Front. Immunol. 2019, 10, 2522. [Google Scholar] [CrossRef] [PubMed]
  54. Panda, A.; Chen, S.; Shaw, A.C.; Allore, H.G. Statistical approaches for analyzing immunologic data of repeated observations: A practical guide. J. Immunol. Methods 2013, 398–399, 19–26. [Google Scholar] [CrossRef] [PubMed]
  55. Grossman, Z.; Greenblatt, C.L.; Cohen, I.R. Parasite immunology and lymphocyte population dynamics. J. Theor Biol. 1986, 121, 129–139. [Google Scholar] [CrossRef] [PubMed]
  56. Jurema, M.W.; Bracero, N.J.; Garcia, J.E. Fine tuning cycle day 3 hormonal assessment of ovarian reserve improves in vitro fertilization outcome in gonadotropin-releasing hormone antagonist cycles. Fertil. Steril. 2003, 80, 1156–1561. [Google Scholar] [CrossRef] [PubMed]
  57. Corey, L.; Beyrer, C.; Cohen, M.S.; Michael, N.L.; Bedford, T.; Rolland, M. SARS-CoV-2 Variants in Patients with Immunosuppression. N. Engl. J. Med. 2021, 385, 562–566. [Google Scholar] [CrossRef] [PubMed]
  58. Thng, Z.X.; De Smet, M.D.; Lee, C.S.; Gupta, V.; Smith, J.R.; McCluskey, P.J.; Thorne, J.E.; Kempen, J.H.; Zierhut, M.; Nguyen, Q.D.; et al. COVID-19 and immunosuppression: A review of current clinical experiences and implications for ophthalmology patients taking immunosuppressive drugs. Br. J. Ophthalmol. 2021, 105, 306–310. [Google Scholar] [CrossRef] [PubMed]
  59. Goldman, J.D.; Robinson, P.C.; Uldrick, T.S.; Ljungman, P. COVID-19 in immunocompromised populations: Implications for prognosis and repurposing of immunotherapies. J. Immunother. Cancer 2021, 9, e002630. [Google Scholar] [CrossRef] [PubMed]
  60. Tassone, D.; Thompson, A.; Connell, W.; Lee, T.; Ungaro, R.; An, P.; Ding, Y.; Ding, N.S. Immunosuppression as a risk factor for COVID-19: A meta-analysis. Intern. Med. J. 2021, 51, 199–205. [Google Scholar] [CrossRef] [PubMed]
  61. Baek, M.S.; Lee, M.-T.; Kim, W.-Y.; Choi, J.C.; Jung, S.-Y. COVID-19-related outcomes in immunocompromised patients: A nationwide study in Korea. PLoS ONE 2021, 16, e0257641. [Google Scholar] [CrossRef] [PubMed]
  62. Abbasi, J. Researchers Tie Severe Immunosuppression to Chronic COVID-19 and Virus Variants. JAMA 2021, 325, 2033–2035. [Google Scholar] [CrossRef] [PubMed]
  63. Kim, A.H.J.; Sparks, J.A. Immunosuppression and SARS-CoV-2 breakthrough infections. Lancet Rheumatol. 2022, 4, e379–e380. [Google Scholar] [CrossRef] [PubMed]
  64. Di Felice, G.; Visci, G.; Teglia, F.; Angelini, M.; Boffetta, P. Effect of cancer on outcome of COVID-19 patients: A systematic review and meta-analysis of studies of unvaccinated patients. eLife 2022, 11, e74634. [Google Scholar] [CrossRef]
  65. Bilich, T.; Roerden, M.; Maringer, Y.; Nelde, A.; Heitmann, J.S.; Dubbelaar, M.L.; Peter, A.; Hörber, S.; Bauer, J.; Rieth, J.; et al. Preexisting and Post-COVID-19 Immune Responses to SARS-CoV-2 in Patients with Cancer. Cancer Discov. 2021, 11, 1982–1995. [Google Scholar] [CrossRef] [PubMed]
  66. Primorac, D.; Vrdoljak, K.; Brlek, P.; Pavelič, E.; Molnar, V.; Matišić, V.; Erceg Ivkošić, I.; Parčina, M. Adaptive Immune Responses and Immunity to SARS-CoV-2. Front. Immunol. 2022, 13, 848582. [Google Scholar] [CrossRef] [PubMed]
  67. Bobcakova, A.; Petriskova, J.; Vysehradsky, R.; Kocan, I.; Kapustova, L.; Barnova, M.; Diamant, Z.; Jesenak, M. Immune Profile in Patients With COVID-19: Lymphocytes Exhaustion Markers in Relationship to Clinical Outcome. Front. Cell. Infect. Microbiol. 2021, 11, 646688. [Google Scholar] [CrossRef] [PubMed]
  68. Bange, E.M.; Han, N.A.; Wileyto, P.; Kim, J.; Gouma, S.; Robinson, J.; Greenplate, A.R.; Hwee, M.A.; Porterfield, F.; Owoyemi, O.; et al. CD8+ T cells contribute to survival in patients with COVID-19 and hematologic cancer. Nat. Med. 2021, 27, 1280–1289. [Google Scholar] [CrossRef] [PubMed]
  69. Başcı, S.; Ata, N.; Altuntaş, F.; Yiğenoğlu, T.N.; Dal, M.S.; Korkmaz, S.; Namdaroğlu, S.; Baştürk, A.; Hacıbekiroğlu, T.; Doğu, M.J.; et al. Patients with hematologic cancers are more vulnerable to COVID-19 compared to patients with solid cancers. Intern. Emerg. Med. 2022, 17, 135–139. [Google Scholar] [CrossRef] [PubMed]
  70. Martínez-Urbistondo, M.; Gutiérrez-Rojas, A.; Andrés, A.; Gutiérrez, I.; Escudero, G.; García, S.; Gutiérrez, A.; Sánchez, E.; Herráiz, J.; De La Fuente, S.; et al. Severe Lymphopenia as a Predictor of COVID-19 Mortality in Immunosuppressed Patients. J. Clin. Med. 2021, 10, 3595. [Google Scholar] [CrossRef] [PubMed]
  71. Garbo, R.; Valent, F.; Gigli, G.L.; Valente, M. Pre-Existing Lymphopenia Increases the Risk of Hospitalization and Death after SARS-CoV-2 Infection. Infect. Dis. Rep. 2022, 14, 20–25. [Google Scholar] [CrossRef]
  72. Giacaman, A.; Henriquez, W.; Tolosa, G.; Prado, A.; Jerez, R.; Reveco, Y.; Martínez, C.; Baumert, C.; Rodríguez, B.; Sanhueza, B.; et al. Hematological abnormalities in immunosuppressed patients with COVID-19: Evidence from a single center. A cross sectional study. Int. Immunopharmacol. 2022, 109, 108862. [Google Scholar] [CrossRef] [PubMed]
  73. Mara, C.A.; Peugh, J.L. Validity of Data Collected from Randomized Behavioral Clinical Trials During the COVID-19 Pandemic. J. Pediatr. Psychol. 2020, 45, 971–976. [Google Scholar] [CrossRef] [PubMed]
  74. Rubina, K.; Shmakova, A.; Shabanov, A.; Andreev, Y.; Borovkova, N.; Kulabukhov, V.; Evseev, A.; Popugaev, K.; Petrikov, S.; Semina, E. Novel prognostic determinants of COVID-19-related mortality: A pilot study on severely-ill patients in Russia. PLoS ONE 2022, 17, e0264072. [Google Scholar] [CrossRef] [PubMed]
  75. Zein, J.G.; Strauss, R.; Attaway, A.H.; Hu, B.; Milinovich, A.; Jawhari, N.; Chamat, S.S.; Ortega, V.E. Eosinophilia Is Associated with Improved COVID-19 Outcomes in Inhaled Corticosteroid-Treated Patients. J. Allergy Clin. Immunol. Pract. 2022, 10, 742–750.e14. [Google Scholar] [CrossRef] [PubMed]
  76. Zhao, Y.; Nie, H.X.; Hu, K.; Wu, X.-J.; Zhang, Y.-T.; Wang, M.-M.; Wang, T.; Zheng, Z.-S.; Li, X.-C.; Zeng, S.-L. Abnormal immunity of non-survivors with COVID-19: Predictors for mortality. Infect. Dis. Poverty 2020, 9, 108. [Google Scholar] [CrossRef] [PubMed]
  77. Qin, R.; He, L.; Yang, Z.; Jia, N.; Chen, R.; Xie, J.; Fu, W.; Chen, H.; Lin, X.; Huang, R.; et al. Identification of Parameters Representative of Immune Dysfunction in Patients with Severe and Fatal COVID-19 Infection: A Systematic Review and Meta-analysis. Clinic. Rev. Allerg. Immunol. 2022, 64, 33–65. Available online: https://doi-org.libproxy.unm.edu/10.1007/s12016-021-08908-8 (accessed on 3 March 2024). [CrossRef] [PubMed]
  78. Limon-de la Rosa, N.; Cervantes-Alvarez, E.; Méndez-Guerrero, O.; Gutierrez Gallardo, M.A.; Kershenobich, D.; Navarro-Alvarez, N. Time-Dependent Changes of Laboratory Parameters as Independent Predictors of All-Cause Mortality in COVID-19 Patients. Biology 2022, 11, 580. [Google Scholar] [CrossRef]
  79. Ingraham, N.E.; Lotfi-Emran, S.; Thielen, B.K.; Techar, K.; Morris, R.S.; Holtan, S.G.; Dudley, R.A.; Tignanelli, C.J. Immunomodulation in COVID-19. Lancet Respir. Med. 2020, 8, 544–546. [Google Scholar] [CrossRef] [PubMed]
  80. Van Regenmortel, M.H. Reductionism and complexity in molecular biology. Scientists now have the tools to unravel biological and overcome the limitations of reductionism. EMBO Rep. 2004, 5, 1016–1020. [Google Scholar] [CrossRef] [PubMed]
  81. Fair, J.M.; Rivas, A.L. Systems Biology and ratio-based, real-time disease surveillance. Transb. Emerg. Dis. 2015, 62, 437–445. [Google Scholar] [CrossRef] [PubMed]
  82. Heudel, P.; Favier, B.; Solodky, M.L.; Assaad, S.; Chaumard, N.; Tredan, O.; Bachelot, T.; Ray-Coquard, I.; Russias, B.; Fournier, M.L.; et al. Survival and risk of COVID-19 after SARS-CoV-2 vaccination in a series of 2391 cancer patients. Eur. J. Cancer 2022, 165, 174–183. [Google Scholar] [CrossRef]
Biomedicines 12 00871 g001
Figure 1. Non-structured data analysis. Total leukocyte counts (WBC) and relative percentages of lymphocytes, neutrophils and monocytes were assessed considering only the data available at hospitalization day 1 (A,B) and all temporal observations (C,D). In all comparisons, non-survivors and survivors displayed overlapping data distributions, which prevented their differentiation (rectangles, (A–D)).
Figure 1. Non-structured data analysis. Total leukocyte counts (WBC) and relative percentages of lymphocytes, neutrophils and monocytes were assessed considering only the data available at hospitalization day 1 (A,B) and all temporal observations (C,D). In all comparisons, non-survivors and survivors displayed overlapping data distributions, which prevented their differentiation (rectangles, (A–D)).
Biomedicines 12 00871 g001
Biomedicines 12 00871 g002
Figure 2. Day-specific, non-structured data analysis. The blood concentrations of IL-6 and hemoglobin (Hb) (A,B) as well as the values of a complex indicator that captured numerous relationships among leukocytes (C) were determined at each post-hospitalization day for non-survivors and survivors. The overlapping data distributions of IL-6 and Hb prevented the differentiation of outcomes (A,B). However, for the non-cellular, non-structured (IL-6, Hb) and leukocyte-related (AAT) structured variables in the first 3 days, survivors displayed a range of values for the complex indicator not including non-survivors (rectangle, (C)).
Figure 2. Day-specific, non-structured data analysis. The blood concentrations of IL-6 and hemoglobin (Hb) (A,B) as well as the values of a complex indicator that captured numerous relationships among leukocytes (C) were determined at each post-hospitalization day for non-survivors and survivors. The overlapping data distributions of IL-6 and Hb prevented the differentiation of outcomes (A,B). However, for the non-cellular, non-structured (IL-6, Hb) and leukocyte-related (AAT) structured variables in the first 3 days, survivors displayed a range of values for the complex indicator not including non-survivors (rectangle, (C)).
Biomedicines 12 00871 g002
Biomedicines 12 00871 g003
Figure 3. Day 1-specific assessments of kurtosis. Histograms of variables collected on the first hospitalization day did not reveal kurtosis (values equal to or less than 0.13) when, regardless of outcome, non-structured, leukocyte-related data were investigated (A,B). While non-survivors did not show kurtosis when other variables were investigated (values equal to or less than 0.29, (C,E,G)), the same variables displayed L-shaped patterns when survivors were evaluated (values equal to or higher than 16.89, (D,F,H)).
Figure 3. Day 1-specific assessments of kurtosis. Histograms of variables collected on the first hospitalization day did not reveal kurtosis (values equal to or less than 0.13) when, regardless of outcome, non-structured, leukocyte-related data were investigated (A,B). While non-survivors did not show kurtosis when other variables were investigated (values equal to or less than 0.29, (C,E,G)), the same variables displayed L-shaped patterns when survivors were evaluated (values equal to or higher than 16.89, (D,F,H)).
Biomedicines 12 00871 g003
Biomedicines 12 00871 g004
Figure 4. Longitudinal assessments of kurtosis—I. When all longitudinal observations were assessed together (n = 291), no kurtosis was observed (all estimates were equal to or less than 1.22) when, regardless of outcome, the percentages of lymphocytes, monocytes or neutrophils were separately investigated (A–F).
Figure 4. Longitudinal assessments of kurtosis—I. When all longitudinal observations were assessed together (n = 291), no kurtosis was observed (all estimates were equal to or less than 1.22) when, regardless of outcome, the percentages of lymphocytes, monocytes or neutrophils were separately investigated (A–F).
Biomedicines 12 00871 g004
Biomedicines 12 00871 g005
Figure 5. Longitudinal assessments of kurtosis—II. When non-cellular indicators as well as leukocyte-related complex indicators were investigated and all longitudinal observations were considered (n = 291), survivors showed kurtosis estimates ranging from 13.41 to 253.78. In contrast, non-survivors displayed smaller kurtosis estimates (A–F).
Figure 5. Longitudinal assessments of kurtosis—II. When non-cellular indicators as well as leukocyte-related complex indicators were investigated and all longitudinal observations were considered (n = 291), survivors showed kurtosis estimates ranging from 13.41 to 253.78. In contrast, non-survivors displayed smaller kurtosis estimates (A–F).
Biomedicines 12 00871 g005
Biomedicines 12 00871 g006
Figure 6. Longitudinal assessments of kurtosis—III (dynamics of paired observations, separate analyses). When the net difference between day 2 and day 1 values of 92 pairs of observations was considered, non-survivors did not show evidence of kurtosis (values ranging between −0.44 and 8.37), while survivors exhibited kurtosis (values between 18.37 and 80.99, (A–F)).
Figure 6. Longitudinal assessments of kurtosis—III (dynamics of paired observations, separate analyses). When the net difference between day 2 and day 1 values of 92 pairs of observations was considered, non-survivors did not show evidence of kurtosis (values ranging between −0.44 and 8.37), while survivors exhibited kurtosis (values between 18.37 and 80.99, (A–F)).
Biomedicines 12 00871 g006
Biomedicines 12 00871 g007
Figure 7. Longitudinal assessments of kurtosis—IV (dynamics of paired observations, integrated analysis). Survivors displayed a larger dispersion of values than non-survivors when the difference between the day 2 and day 1 data values of survivors and non-survivors was simultaneously investigated in a single plot (92 paired observations). Both non-cellular indicators (such as IL-6, D-dimer and PCT, (A–C)) and complex (structured) indicators derived from leukocyte data (AAT, BBF, BBT, (D–F)) revealed that non-survivors did not vary much in values over time (the difference of day 2 minus day 1 approached zero). In contrast, all other variables (except PCT) revealed broad oscillations over time when survivors were considered, which expressed both above- and below-zero values. This assessment may be interpreted as a single and perpendicular perspective of the data reported separately for survivors and non-survivors in Figure 6.
Figure 7. Longitudinal assessments of kurtosis—IV (dynamics of paired observations, integrated analysis). Survivors displayed a larger dispersion of values than non-survivors when the difference between the day 2 and day 1 data values of survivors and non-survivors was simultaneously investigated in a single plot (92 paired observations). Both non-cellular indicators (such as IL-6, D-dimer and PCT, (A–C)) and complex (structured) indicators derived from leukocyte data (AAT, BBF, BBT, (D–F)) revealed that non-survivors did not vary much in values over time (the difference of day 2 minus day 1 approached zero). In contrast, all other variables (except PCT) revealed broad oscillations over time when survivors were considered, which expressed both above- and below-zero values. This assessment may be interpreted as a single and perpendicular perspective of the data reported separately for survivors and non-survivors in Figure 6.
Biomedicines 12 00871 g007
Biomedicines 12 00871 g008
Figure 8. Detection of simple, bi-dimensional relationships associated with survival. While on days 1 and 2, survivors exhibited a broad data interval that facilitated relationships between non-structured variables (neutrophils and monocytes), non-survivors lacked such a data range (rectangles, (A–D)). A similar absence was observed when lymphocytes and monocytes were assessed (rectangles, (E–H)).
Figure 8. Detection of simple, bi-dimensional relationships associated with survival. While on days 1 and 2, survivors exhibited a broad data interval that facilitated relationships between non-structured variables (neutrophils and monocytes), non-survivors lacked such a data range (rectangles, (A–D)). A similar absence was observed when lymphocytes and monocytes were assessed (rectangles, (E–H)).
Biomedicines 12 00871 g008
Biomedicines 12 00871 g009
Figure 9. Detection of complex, bi-dimensional relationships associated with survival. When structured, leukocyte-related variables were investigated at days 1 and 2 in survivors and non-survivors, survivors exhibited two data ranges that were lacking when non-survivors were evaluated (A–D). Because the analysis of non-structured data did not reveal two separate data ranges associated with survival (shown in Figure 8), the analysis of structured (complex) indicators extracted more information from the same data.
Figure 9. Detection of complex, bi-dimensional relationships associated with survival. When structured, leukocyte-related variables were investigated at days 1 and 2 in survivors and non-survivors, survivors exhibited two data ranges that were lacking when non-survivors were evaluated (A–D). Because the analysis of non-structured data did not reveal two separate data ranges associated with survival (shown in Figure 8), the analysis of structured (complex) indicators extracted more information from the same data.
Biomedicines 12 00871 g009
Biomedicines 12 00871 g010
Figure 10. Influence of time on the complexity of bi-dimensional assessments associated with survival. A bi-dimensional, outcome-related, immunological and temporal assessment of structured data is not influenced by time: each of the two data ranges occupied only by survivors included both day 1 and day 2 observations (A–C). The two subsets of survivor-only observations (here named A and B) were biologically valid and distinguishable: subsets A and B differed from one another both in terms of non-structured data (non-overlapping percentages of lymphocytes, neutrophils and monocytes) and in several ratios, including the small leukocyte (lymphocyte and neutrophil)/monocyte (SL/M), neutrophil/monocyte (N/M) and neutrophil/lymphocyte (N/L) ratios (D,E).
Figure 10. Influence of time on the complexity of bi-dimensional assessments associated with survival. A bi-dimensional, outcome-related, immunological and temporal assessment of structured data is not influenced by time: each of the two data ranges occupied only by survivors included both day 1 and day 2 observations (A–C). The two subsets of survivor-only observations (here named A and B) were biologically valid and distinguishable: subsets A and B differed from one another both in terms of non-structured data (non-overlapping percentages of lymphocytes, neutrophils and monocytes) and in several ratios, including the small leukocyte (lymphocyte and neutrophil)/monocyte (SL/M), neutrophil/monocyte (N/M) and neutrophil/lymphocyte (N/L) ratios (D,E).
Biomedicines 12 00871 g010
Biomedicines 12 00871 g011
Figure 11. Three-dimensional, long-term, biologically validated detection of structured data ranges associated with survival. More information was retrieved from the same data when three-dimensional relationships were investigated with data collected on all testing days (n = 291). Three non-randomly distributed data subsets—perpendicular to one another—were populated by survivors (A). Data oscillations (values far from zero) indicated protective responses if they were observed in the first three hospitalization days (B). The combination of spatial and outcome (quantitative and qualitative) data identified three data subsets (named I–III) composed of only survivors (C). All three (I–III) data subsets were validated: they differed from one another by non-overlapping intervals of neutrophil or monocyte percentages or monocyte/lymphocyte ratios (D). Therefore, immuno-competence was an early response that involved, at least, three complex immunological functions that could involve all cell types.
Figure 11. Three-dimensional, long-term, biologically validated detection of structured data ranges associated with survival. More information was retrieved from the same data when three-dimensional relationships were investigated with data collected on all testing days (n = 291). Three non-randomly distributed data subsets—perpendicular to one another—were populated by survivors (A). Data oscillations (values far from zero) indicated protective responses if they were observed in the first three hospitalization days (B). The combination of spatial and outcome (quantitative and qualitative) data identified three data subsets (named I–III) composed of only survivors (C). All three (I–III) data subsets were validated: they differed from one another by non-overlapping intervals of neutrophil or monocyte percentages or monocyte/lymphocyte ratios (D). Therefore, immuno-competence was an early response that involved, at least, three complex immunological functions that could involve all cell types.
Biomedicines 12 00871 g011
Biomedicines 12 00871 g012
Figure 12. Three-dimensional, two-day long, biologically validated detection of structured data ranges associated with survival. Structured data collected on days 1 and 2 (n = 195) displayed several 3D patterns that facilitated numerous inferences. For example, every data point detected within a data subset perpendicular to the remaining data prognosticated survival (rectangle, (A)). Personalized prognoses were supported through redundant analyses: the predictions generated by plot (A) were corroborated by the three data subsets distinguished in plot (B). The three data subsets identified in (B) were biologically valid: each subset associated with survival (‘left’, ‘top’ and ‘right’) differed from one another by two or more variables (C). Because this differentiation and validation corroborated the findings reported in Figure 11, it is concluded that the assessment of either day 1 only or days 1 and 2 conveyed similar information. Other structured indicators showed additional potential targets of immunomodulatory therapies and facilitated earlier evaluations of treatments. For example, the temporal data directionality described by a single (one data point-wide) line of observations showed wo desirable trajectories (temporal data movements that followed the arrows, (D)). Because this assessment (which was based on data collected over 48 h) corroborated the findings reported both in the analysis of all temporal data and the one that considered only day 1 data, it is concluded than the analysis of complex (structured) leukocyte-related, alone or together with the analysis of other non-cellular hematological variables, may assess immuno-competence, facilitate personalized prognosis and (when temporal data directionality is considered) evaluate therapies earlier.
Figure 12. Three-dimensional, two-day long, biologically validated detection of structured data ranges associated with survival. Structured data collected on days 1 and 2 (n = 195) displayed several 3D patterns that facilitated numerous inferences. For example, every data point detected within a data subset perpendicular to the remaining data prognosticated survival (rectangle, (A)). Personalized prognoses were supported through redundant analyses: the predictions generated by plot (A) were corroborated by the three data subsets distinguished in plot (B). The three data subsets identified in (B) were biologically valid: each subset associated with survival (‘left’, ‘top’ and ‘right’) differed from one another by two or more variables (C). Because this differentiation and validation corroborated the findings reported in Figure 11, it is concluded that the assessment of either day 1 only or days 1 and 2 conveyed similar information. Other structured indicators showed additional potential targets of immunomodulatory therapies and facilitated earlier evaluations of treatments. For example, the temporal data directionality described by a single (one data point-wide) line of observations showed wo desirable trajectories (temporal data movements that followed the arrows, (D)). Because this assessment (which was based on data collected over 48 h) corroborated the findings reported both in the analysis of all temporal data and the one that considered only day 1 data, it is concluded than the analysis of complex (structured) leukocyte-related, alone or together with the analysis of other non-cellular hematological variables, may assess immuno-competence, facilitate personalized prognosis and (when temporal data directionality is considered) evaluate therapies earlier.
Biomedicines 12 00871 g012
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

Kempaiah, P.; Libertin, C.R.; Chitale, R.A.; Naeyma, I.; Pleqi, V.; Sheele, J.M.; Iandiorio, M.J.; Hoogesteijn, A.L.; Caulfield, T.R.; Rivas, A.L. Decoding Immuno-Competence: A Novel Analysis of Complete Blood Cell Count Data in COVID-19 Outcomes. Biomedicines 2024, 12, 871. https://doi.org/10.3390/biomedicines12040871

AMA Style

Kempaiah P, Libertin CR, Chitale RA, Naeyma I, Pleqi V, Sheele JM, Iandiorio MJ, Hoogesteijn AL, Caulfield TR, Rivas AL. Decoding Immuno-Competence: A Novel Analysis of Complete Blood Cell Count Data in COVID-19 Outcomes. Biomedicines. 2024; 12(4):871. https://doi.org/10.3390/biomedicines12040871

Chicago/Turabian Style

Kempaiah, Prakasha, Claudia R. Libertin, Rohit A. Chitale, Islam Naeyma, Vasili Pleqi, Johnathan M. Sheele, Michelle J. Iandiorio, Almira L. Hoogesteijn, Thomas R. Caulfield, and Ariel L. Rivas. 2024. "Decoding Immuno-Competence: A Novel Analysis of Complete Blood Cell Count Data in COVID-19 Outcomes" Biomedicines 12, no. 4: 871. https://doi.org/10.3390/biomedicines12040871

APA Style

Kempaiah, P., Libertin, C. R., Chitale, R. A., Naeyma, I., Pleqi, V., Sheele, J. M., Iandiorio, M. J., Hoogesteijn, A. L., Caulfield, T. R., & Rivas, A. L. (2024). Decoding Immuno-Competence: A Novel Analysis of Complete Blood Cell Count Data in COVID-19 Outcomes. Biomedicines, 12(4), 871. https://doi.org/10.3390/biomedicines12040871

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