1. Introduction
When the World Health Organization (WHO) declared on May 5, 2023, the end of the emergency for Covid-19, there was still an increase in confirmed COVID-19 deaths at several places in African, American, European, Southeastern Asia, and Western Pacific regions [
1]. Currently, over 7 million deaths and over 775 million confirmed cases have been reported globally [
2]. It should be noted that these numbers are underestimations since many countries have stopped or changed the frequency of reporting [
2].
Some people cannot fully recover after the COVID-19 disease, presenting long-term symptoms which are usually named “Long-COVID-19”, “Long COVID syndrome”, “post-COVID-19 condition”, or “post-acute COVID-19 syndrome” [
3,
4]. In short, the term “Long COVID” is commonly used to describe signs and symptoms that persist after acute COVID-19 and this term is going to be used here.
COVID-19 disease affects a multitude of organs and systems, and Long COVID, as a post-infection condition, relates to various systems such as the respiratory, cardiovascular, and nervous systems, underpinning the need for elucidating the exact pathophysiological mechanism [
5]. COVID-19 not only reduces significantly axial blood microvessel velocity but also has the devastating effect of causing extensive microthrombosis, as was shown by recent quantitative work on post-COVID patients [
6].
Soon after, a normative velocity model in the exchange microvessels was proposed and tested successfully as a disease discriminator on post-COVID conjunctival data [
7]. However, microvascular loss (ML), due to microthrombosis, was not taken into account in the normative velocity model because only functional microvessels (with blood flow) were included in the velocity measurements. In reality, the tissue blood supply was even lower due to micro-occlusion. In addition, data in those papers [
6,
7] were taken only from the conjunctival tissue and a limited geographical region.
In this work, Long COVID symptoms are described first, according to the available data. Then, since most of these symptoms can be attributed to microvascular deprivation of proper blood supply, the concept of microvascular blood supply (S) per flat tissue area is defined. With S, the contribution of both the microvascular hemodynamic decrease (HD) and the ML, is taken into account in post-COVID patients. Tissue blood supply reduction (SR) for a long time is the cause of persistent lower mass diffusion rates and therefore multiple tissue hypoxia and undernutrition which can explain most of the Long COVID symptoms. Evidence in support of this mechanism is presented from previously published case-control microcirculatory data from multiple tissues and geographical areas around the world.
2. Background on Long COVID and Symptoms
Paul Garner was the first to describe in a blog of the British Medical Journal (May 5, 2020) his severe post-COVID symptoms, ongoing 7 weeks after infection. The term “Long COVID” was first used two weeks later by Elisa Perego (May 20, 2020), to describe the cyclical, and multiphasic Long-term Covid illness [
8]. Carfi et al. [
9] were the first to report symptoms from a cohort of 143 patients that persisted 36 days (on average) after hospital discharge and 87% of the patients suffered persistence of at least one symptom, especially fatigue (53%), dyspnea (43%), joint pain (27%), and chest pain (22%). In a multistate survey of 270 COVID-19 adults with milder outpatient illness [
10], only 65% had returned to their usual health state when interviewed 2-3 weeks after positive testing. In comparison, more than 90% of influenza outpatients recover within 2 weeks after positive testing [
10]. Puntmann et al. [
11] were the first to report a cardiovascular magnetic resonance imaging cohort study of 100 patients who recently recovered from COVID-19, and in their work, cardiac involvement and ongoing myocardial inflammation were demonstrated in 78% and 60% of the patients, respectively. In addition, they found a statistically significant discrimination between controls and post-COVIDs using C-reactive protein (CRP), a blood biomarker that positively correlates with the incidence of thrombosis.
According to a 6-month retrospective cohort study of 273,618 survivors of COVID-19 [
12], 57% had one or more Long COVID features recorded during the whole 6-month period. The most commonly reported symptoms by order of percentage were anxiety/depression (23%), abnormal breathing (19%), abdominal symptoms (16%), fatigue/malaise (13%), chest/throat pain (13%), other pain (12%), headache (9%), cognitive symptoms (8%), and myalgia (3%). Huang et al. [
13] reported that 76% of 1733 COVID-19 patients discharged from the hospital had at least one symptom 6 months after the acute infection. There were similar reports for milder cases of COVID-19 without hospitalization [
3], where 61% of patients older than 46 years had persistent symptoms at a 6-month follow-up.
Rezel-Potts et al. [
14] conducted a large case-control study with 428,650 COVID-19 patients with a maximum follow-up time of 52 weeks after diagnosis and found that cardiovascular disease was increased in the early COVID phase (first 4 weeks) with principal manifestations of pulmonary embolism, atrial arrhythmias, and venous thromboses (5-fold increase). The 5-fold increase in venous thromboses aligns closely with the 6-fold increase in conjunctival microthrombosis reported in COVID-19 survivors soon after hospital discharge [
6]. Katsoularis et al. [
15] conducted a nationwide matched cohort study and found that COVID-19 is an independent risk factor for deep vein thrombosis, pulmonary embolism, and bleeding and that the risk of these outcomes is increased for three, six, and two months after COVID-19, respectively.
Other researchers have reported fatigue as a major manifestation of Long COVID with a high prevalence (60-70%) [16-19]. Poor sleep quality is another common manifestation of Long COVID [
16,
19,
20], and there were demonstrations of chronic pain either as a newly developed condition or a worsening of preexisting chronic pain [
21]. The important neurological and psychiatric symptoms observed during Long COVID greatly puzzle scientists worldwide [
22]. In addition, a recent retrospective study [
16] of 287 patients with Long COVID symptoms found that fatigue, sleepiness, and sleep quality were SARS-CoV-2 variant independent. This supports the view of a common Long COVID pathophysiological mechanism for all virus variants, which is presented in the next section.
3. Materials and Methods
3.1. The Normative Range (NR) Diagram
The recently introduced velocity-diffusion (V-J) equation [
23] is based on the well-known flow-diffusion (Q-J) equation [
24], which is one of the fundamental equations of vascular physiology. The V-J equation describes the relation of the axial blood microvessel velocity (V) with the mass diffusion rate (J) of each blood solute through the walls of a microvessel (
Figure 1).
After measuring axial velocities at many microvessels with the same diameter, from the same tissue of many healthy subjects (control group), an average axial velocity (AV) can be estimated statistically corresponding to an average mass diffusion rate (AJ). Since the AV is determined statistically, a normative range (NR) of the AV can be defined for microvessels of this diameter, corresponding to normal resting conditions, in a normative range (NR) diagram (
Figure 1). The UR (Underemic Range) is a range of average axial velocities corresponding to underemic conditions on the left of the NR, and the HR (Hyperemic Range) is the range of average axial velocities corresponding to hyperemic conditions on the right of the NR. Underemic conditions are indicative of hypoxia, malnutrition, and disease. The NR diagram is the backbone of the normative microvascular velocity model and can be extended to many diameters [
7].
3.2. Hemodynamic Decrease (HD)
In case-control hemodynamic studies, when a statistically significant reduction in the average blood velocity is detected in the case (post-COVID) group, it is direct proof that the average axial velocity (AV) of the case group is shifted to the left, in the underemic range (
Figure 1). To quantify this shift, “Hemodynamic Decrease” (HD) is defined here as:
Where AV
CONTROL and AV
COVID are
the average axial blood velocities in the control group and the post-COVID
group, respectively. HD ranges between 0% and 100% (for AV
CONTROL >
AV
COVID). HD can be quantified with various techniques such as
Conjunctival Video Capillaroscopy (CVC) [
25,
26],
Laser Doppler Flowmetry (LDF) [
27], and
special Magnetic Resonance Imaging (MRI) [
28].
3.3. Microvascular Loss (ML)
“Microvascular loss” (ML) is defined generally as a deprivation of microvessels in a given tissue. ML can be detected in various tissues with various techniques such as Optical Coherence Tomography Angiography (OCTA), CVC, sublingual video capillaroscopy (SVC), and Nailfold Video Capillaroscopy (NVC).
In the following subsections, the well-known indices of vessel density (VD), foveal avascular zone (FAZ), capillary density (CD), and percentage of perfused vessels (PPV) are utilized to make quantitative estimates of ML from case-control studies.
3.3.1. Vessel Density Reduction (VDR)
Vessel density (VD) is defined as the total functional microvessel area per unit area in the region of measurement and is measured in percent (%). Functional microvessels are considered those with blood flow since most techniques require blood flow to detect a microvessel [
29]. VD reduction (VDR) is defined as:
Where VDCONTROL and VDCOVID
are the averages in the control and the post-COVID group, respectively. For VDCONTROL
> VDCOVID, VDR ranges between 0% and 100%.
3.3.2. Foveal Avascular zone Enlargement (FAZE)
The foveal avascular zone (FAZ) is an approximately circular area without any blood vessels with a diameter of about 0.5 mm, located at the center of the fovea in the eye fundus. FAZ however, gets proper supply due to the special tissue architecture of that retinal section (foveal pit or foveola about 0.35 mm in diameter). FAZ enlargement (FAZE) is defined as:
Where FAZCONTROL and FAZCOVID are the averages in the control and the post-COVID group, respectively. FAZE can be considered a kind of microvascular loss in the most critical area of vision.
3.3.3. Capillary Density Reduction (CDR)
Capillary density (CD) is defined as the total capillary number per linear millimeter and is measured by nailfold video capillaroscopy (NVC). Capillary density reduction (CDR) is defined as:
Where CDCONTROL and CDCOVID are the averages in the control and the post-COVID group, respectively. CDR ranges between 0% and 100 (for CDCONTROL > CDCOVID).
3.3.4. Percentage of Perfused Vessels Reduction (PPVR)
The percentage of perfused vessels (PPV) is defined as the percentage of the number of perfused (functional) vessels over the total number of vessels, at a given microvascular area (field of view). The PPV is an index inversely related to the microvascular damage (thrombotic microoclusion) sustained by a tissue. The percentage of perfused vessels reduction (PPVR) is defined as:
Where PPVCONTROL and PPVCOVID are the averages in the control and the post-COVID group, respectively. PPVR ranges between 0% and 100% (for PPVCONTROL > PPVCOVID).
3.4. Blood Supply Reduction (SR)
In case-control microvascular studies, when a statistically significant ML is detected in the tissue of the case group, it is evident that there is an oxygen and nutrient undersupply in that tissue area. The current difficulty with the normative velocity model of
Figure 1 is that ML is not taken into account. It is as if
Figure 1 has many “zero points” (V = 0, J = 0) which are not measured. However, the higher the average zero points, the higher the tissue undersupply.
To overcome this difficulty, the blood supply (S) per flat tissue area is defined. For a given flat tissue surface area (field of view) of a subject, the blood supply S is the sum of all the blood volume flows (Q1 + Q2 + … + Qn) of the n functional exchange microvessels. This can be written as S = n Q, where Q is the average blood volume flow of the n microvessels. This can be written as S = n b V, where b is a parameter for the average cross-sectional area (m2) of the n microvessels and V is the average axial blood velocity (m/s) of the n microvessels. In this way, changes in the average number of functional microvessels can be incorporated into the S calculations.
In a case-control study with post-COVID patients, the average blood supply per flat tissue area for the control group can be expressed as SCONTROL = nCONTROL•bCONTROL•VCONTROL, where SCONTROL, nCONTROL, bCONTROL, and VCONTROL are the control group averages for the S, n, b, and V, respectively. In the same manner, for the post-COVID group: SCOVID = nCOVID•bCOVID• VCOVID.
In a situation where an average reduction of alpha (α) percent was observed in the number of microvessels in the post-COVID group, then n
COVID could be expressed as [(100-α)/100]n
CONTROL. If on surplus an average hemodynamic reduction of HD percent was observed in the remaining functional microvessels, then V
COVID could be expressed as [(100-HD)/100]V
CONTROL. Therefore, the average post-COVID group blood supply is written as:
With the condition that the average cross-sectional area of the remaining functional microvessels in the post-COVID group is the same as that in the control group (b
CONTROL = b
COVID), equation 6 is written:
Finally, the blood supply reduction (SR) in percent is defined as:
After inserting equation 7 into equation 8:
Using equation 9, the total supply reduction SR from both the microvascular hemodynamic decrease (HD) and the microvascular loss (α) is taken into account. Of the four microvascular loss quantitative measures (VDR, FAZE, CDR, and PPVR), only two were defined on a flat tissue area: VDR and PPVR. Here it is considered that α can be approximated by VDR or PPVR depending on the available data.
3.5. The Proposed Pathophysiological Microcirculatory Mechanism for Long COVID
The steps and possible outcomes of the COVID-19 disease progress are shown in numbered rectangular boxes (
Figure 2, boxes 1-9). The common framework for the proposed Long COVID mechanism is represented schematically by boxes 6, 7, 8, and 9, which are included in a dashed black line (
Figure 2). Box 6 represents the subgroup of COVID-19 survivors with microvascular loss (ML) and hemodynamical decrease (HD) for several months after testing negative. Box 7 shows that both ML and HD, lead to tissue blood supply reduction (SR). The significant SR has consequently much lower diffusion rates (J) hence, multiple tissue hypoxia and undernutrition (box 8), which can be considered the cause of the Long COVID symptoms (box 9).
3.6. Statistical Analysis
Data were gathered from case-control microcirculatory studies with patients in the recovery phase of the COVID-19 disease (0 to 6 months after hospitalization or after having tested negative). These studies have reported statistically significant HD and ML metrics (VDR, FAZE, CDR, and PPVR). The Microsoft Office EXCEL 2016 software (professional edition) was used for the estimation of median, mean, and standard error of the mean (SEM). The statistical estimations were performed on different tissues since HD and ML metrics were all defined in relative units (percentages) to the control group reference.
5. Discussion
Life is sustained by the diffusion of oxygen and nutrients towards the intracellular space. In animals with a circulatory system, this diffusive process takes mainly place in the exchange microvessels where the adequate supply of blood is of uttermost importance. When the normal tissue blood supply is disrupted at the microvascular level, there are serious side effects.
Early in the COVID-19 pandemic, Polak et al. [
47] reported the importance of microvascular histological patterns and their persistence throughout the clinical course in their systematic review including 131 lung samples from either antemortem or postmortem COVID-19 patients (
Figure 2, boxes 3-5). They identified 76 cases (59%) with microvascular damage (microthrombi) and proposed directed anti-inflammatory, anticoagulant, and/or anticomplement agents. Pretorius et al. [
48] detected amyloid microclots in twenty COVID-19-positive blood samples before patient treatment, and Jung et al. [
49] reported abdominal microcirculatory disorders (
Figure 2, box 3). At the beginning of 2022, coagulopathy in patients hospitalized with COVID-19 was well-documented [
50] and there were reports of enhanced hypercoagulability [
51] and impaired fibrinolysis [
52]. Endothelial dysfunction [
53] is another characteristic of COVID-19 related to microvascular disease and vascular aging [
54].
In this work, a common framework (
Figure 2, boxes 6-9) was presented for the pathophysiological mechanism of the Long COVID symptoms reported by COVID-19 patients after hospital exit or after testing negative. After the acute phase (
Figure 2, box 3), a multisystemic microvasculopathy persists for several months (
Figure 2, box 6), as was presented in the results, for about 800 post-COVID patients (
Table 1,
Table 2,
Table 3,
Table 4 and
Table 5). This can explain the Long COVID symptoms through the quantitative pathophysiological mechanism of the tissue blood supply reduction (SR) (
Figure 2, boxes 7-9).
Regarding post-COVID patients, Patterson et al. [
55] in the middle of 2020, reported that SARS-CoV-2 infection is associated with a wide spectrum of neurological syndromes, and D-dimers were markedly elevated in all patient subgroups. In a cohort study of 100 patients who recovered from COVID-19 [
11], statistically significant discrimination between controls and post-COVIDs was reported using C-reactive protein (CRP), a blood biomarker that positively correlates with the incidence of thrombosis. Pretorius et al. [
56] found microclots that were resistant to fibrinolysis, and a substantial increase in α2-antiplasmin, in blood samples from 11 post-COVID patients at least 2 months after recovery. The presence of amyloid microclots was also reported in blood samples from 80 post-COVID patients [
57] (but without controls). Scheim et al. [
58] proposed a mechanism of microclot formation with a central role of sialylated glycan attachments between SARS-CoV-2 spike protein and red blood cells. Kell and Pretorius [
59] proposed ischemic injury from fibrin amyloid microclots as the primary factor for the Long COVID condition, and in line with this, Astin et al. [
60] proposed lower tissue oxygen availability (chronic hypoxia) as a pathophysiological mechanism of Long COVID. The microclot blocking of capillaries in post-COVID patients was confirmed
in vivo by Koutsiaris et al. [
6] who also proposed microthrombosis as a possible explanation for Long COVID syndrome and speculated the existence of a possibly unknown coagulation factor.
The propositions mentioned above are well-fitted into the pathophysiological context of the present work, with the difference that here, case-control
in vivo quantitative data were gathered, from multiple human tissues and multiple geographical places, with statistically significant differences between post-COVID cases and controls (
Tables 1 to 5,
Figure 2, box 6).
The major finding of a case-control study [
61] with 120 Long COVID individuals 3-4 months after the acute infectious phase was the association of Long COVID with decreased antioxidant defenses as indicated by the lowered total antioxidant capacity of plasma. In addition, a high percentage of the variance in the severity of the Long COVID neuropsychiatric symptoms was explained by the increased C-reactive protein and the oxidative stress toxicity to antioxidant ratio. These findings are in support of the pathophysiological model resulting in tissue hypoxia in this work (
Figure 2, box 8).
In a study with 87 COVID-19 survivors after hospitalization [
62], exercise intolerance was reported, and a hypothesis of lung microvascular injury was made, as a pathophysiological mechanism leading to increased dead space as a fraction of tidal volume during exercise in post-COVID-19 patients. This hypothesis fits very well in the context of this work.
Gareau and Barrett [
63] proposed a role of the impaired microbiota-gut-brain axis signaling in the development of Long COVID but this does not explain the extensive microthrombosis. However, a potential link to microthrombosis may exist in the gut-lung axis through the degradation of the intestinal epithelial cell junctional proteins and hence of the biochemical barrier to the microvessels.
Reiss et al. [
64] explored possible mechanistic pathways between Long COVID and nervous system inflammation. In their review, it was noted that impaired blood flow in the brain due to viral invasion of the microvascular endothelium may be a neuropathological mechanism of Long COVID which agrees with this work.
In three recent Long COVID reviews [
5,
65,
66] the need for understanding the pathophysiological mechanism of Long COVID was underscored, and various possible cellular and molecular mechanisms were described. Among them, the general concept of endothelial dysfunction-damage, microclot formation, persistent microvascular injury, and impairment of oxygen transfer, is in good agreement with this work, which supports it as the most prevalent mechanism.
For an indication of when the microvascular status returns to normal, a recent work [
67] with retinal evaluations at an average of 15.2 ± 6.9 months post SARS-CoV-2 infection, reported no significant differences from controls.
Assuming conservatively that only 10% of the COVID-19-infected people developed Long COVID [
68], at least 77.5 million subjects worldwide have shown Long COVID symptoms [
2]. With insufficient current Long COVID diagnostic and treatment options [
68] more basic and clinical research is needed for understanding Long COVID pathophysiological mechanisms.