Introduction
Symptoms and signs of dysfunction of multiple organ systems may be present during both the acute and post-viral stages of coronavirus disease 2019 (COVID-19) (Al-Aly, Xie, & Bowe, Reference Al-Aly, Xie and Bowe2021; Bell et al., Reference Bell, Catalfamo, Farland, Ernst, Jacobs, Klimentidis and Pogreba-Brown2021; Blomberg et al., Reference Blomberg, Mohn, Brokstad, Zhou, Linchausen, Hansen and Langeland2021; Davis et al., Reference Davis, Assaf, McCorkell, Wei, Low, Re'em and Akrami2021; Huang et al., Reference Huang, Huang, Wang, Li, Ren, Gu and Cao2021a, Reference Huang, Yao, Gu, Wang, Ren, Wang and Cao2021b; Lopez-Leon et al., Reference Lopez-Leon, Wegman-Ostrosky, Perelman, Sepulveda, Rebolledo, Cuapio and Villapol2021; Menges et al., Reference Menges, Ballouz, Anagnostopoulos, Aschmann, Domenghino, Fehr and Puhan2021; Nasserie, Hittle, & Goodman, Reference Nasserie, Hittle and Goodman2021; Sudre et al., Reference Sudre, Murray, Varsavsky, Graham, Penfold, Bowyer and Steves2021; Writing Committee for the COMEBAC Study Group et al., Reference Morin, Savale, Pham, Colle, Figueiredo and Monnet2021). Cognitive deficits, anxiety, depression, sleep disturbances, and post-traumatic stress disorder (PTSD) have been frequently shown to persist for weeks to several months after COVID-19 (Abel et al., Reference Abel, Carr, Ashcroft, Chalder, Chew-Graham, Hope and Pierce2021; Becker et al., Reference Becker, Lin, Doernberg, Stone, Navis, Festa and Wisnivesky2021; Bourmistrova, Solomon, Braude, Strawbridge, & Carter, Reference Bourmistrova, Solomon, Braude, Strawbridge and Carter2021; Damiano et al., Reference Damiano, Guimarães, Cincoto, Rocca, Serafim, Bacchi and Forlenza2022; Hampshire et al., Reference Hampshire, Trender, Chamberlain, Jolly, Grant, Patrick and Mehta2021; Ismael et al., Reference Ismael, Bizario, Battagin, Zaramella, Leal, Torales and Castaldelli-Maia2021; Mazza et al., Reference Mazza, De Lorenzo, Conte, Poletti, Vai, Bollettini and Benedetti2020, Reference Mazza, Palladini, De Lorenzo, Magnaghi, Poletti, Furlan and Benedetti2021; Søraas et al., Reference Søraas, Bø, Kalleberg, Støer, Ellingjord-Dale and Landrø2021; Taquet, Geddes, Husain, Luciano, & Harrison, Reference Taquet, Geddes, Husain, Luciano and Harrison2021), highlighting the relevance of central nervous system (CNS) manifestations in the long COVID syndrome (or post-acute sequelae of SARS-CoV-2 infection, PASC).
As in many disorders that present multiple clinical manifestations (Ahmad et al., Reference Ahmad, Pencina, Schulte, O'Brien, Whellan, Piña and Felker2014; Aktas, Walsh, & Rybicki, Reference Aktas, Walsh and Rybicki2010; Manca, De Marco, Ince, & Venneri, Reference Manca, De Marco, Ince and Venneri2021), the detection of specific patterns of symptom clustering may help to elucidate the pathophysiology of PASC, which remains largely unknown. The identification of patterns of co-occurrence of psychiatric and cognitive symptoms with manifestations pertaining to other organ systems could guide hypotheses as to whether CNS symptoms in PASC may emerge under the influence of unifying pathological mechanisms, such as persistent systemic inflammatory and prothrombotic states (Bornstein et al., Reference Bornstein, Voit-Bak, Donate, Rodionov, Gainetdinov, Tselmin and Straube2021; Mackay, Reference Mackay2021; Perrin et al., Reference Perrin, Riste, Hann, Walther, Mukherjee and Heald2020; Phillips & Williams, Reference Phillips and Williams2021), or chronic dysfunction of the lungs (Sasannejad, Ely, & Lahiri, Reference Sasannejad, Ely and Lahiri2019) or kidneys (Desmond et al., Reference Desmond, Lindner, Troost, Held, Callaway, Oh and Massengill2021).
Based on online surveys applied to 3762 subjects with confirmed or suspected COVID-19, one investigation of temporal profiles of PASC symptoms highlighted a cluster combining CNS manifestations with multiple other organ system symptoms, which shared patterns of persistence up to 7 months after acute disease (Davis et al., Reference Davis, Assaf, McCorkell, Wei, Low, Re'em and Akrami2021). Another study applied statistical methods to directly assess patterns of clustering of different PASC manifestations in subjects with moderate or severe COVID-19, based on data from a large sample of adults (n = 1077) evaluated during in-person visits between 2 and 7 months after in-hospital discharge (Evans et al., Reference Evans, McAuley, Harrison, Shikotra, Singapuri and Sereno2021). The authors of the latter study described four clusters reflecting different levels of PASC severity, and in all clusters there was a close co-occurrence of three CNS symptoms (i.e. depression, anxiety, and PTSD) with manifestations pertaining to other organ systems (i.e. dyspnea, fatigue, and poor physical performance). Cluster severity was directly associated with blood levels of C-reactive protein measured at the same timepoint, possibly reflecting persistent systemic inflammation (Mackay, Reference Mackay2021; Perrin et al., Reference Perrin, Riste, Hann, Walther, Mukherjee and Heald2020). However, this investigation was limited by the inclusion of only three non-CNS manifestations in the cluster analysis design. Moreover, no analyses were conducted to directly account for the association found between C-reactive protein levels and increased body weight (Evans et al., Reference Evans, McAuley, Harrison, Shikotra, Singapuri and Sereno2021), or the potentially confounding effects of other pro-inflammatory comorbidities.
In this study, we investigated the patterns of co-occurrence of 30 multi-organ system symptoms (including six psychiatric manifestations and cognitive complaints) in a relatively large sample of adults (n = 749) assessed by multidisciplinary teams at approximately 6–11 months after hospitalization due to COVID-19. We initially conducted an exploratory factor analysis (EFA) on the ratings of all symptoms, with the prediction that one factor would emerge combining cognitive and psychiatric complaints with multiple symptoms pertaining to other organ systems (Davis et al., Reference Davis, Assaf, McCorkell, Wei, Low, Re'em and Akrami2021; Evans et al., Reference Evans, McAuley, Harrison, Shikotra, Singapuri and Sereno2021). Subsequently, we generated a latent variable combining symptoms of PASC using Item Response Theory (IRT), a statistical approach suitable for scaling multiple health outcomes along one single continuum of severity (latent trait modeling) (Hays, Morales, & Reise, Reference Hays, Morales and Reise2000; Krueger et al., Reference Krueger, Nichol, Hicks, Markon, Patrick, Iacono and McGue2004). This latent trait of PASC symptoms was used to investigate associations with objective signs of persistent physical disability, pulmonary dysfunction, and kidney function impairment. Finally, we investigated whether latent trait severity would be directly related to levels of blood markers of persistent inflammatory and prothrombotic states measured at the same follow-up assessment (C-reactive protein and D-dimer), accounting for the influence of obesity and other pro-inflammatory comorbidities.
Methods
Study design, participants, and procedures
We consecutively invited for a follow-up visit all adult (⩾18 years) patients that had been admitted for at least 24 h as inpatients to Hospital das Clínicas, Faculdade de Medicina, Universidade de São Paulo (HCFMUSP), Brazil, due to laboratory-confirmed COVID-19, over a period of 5 months (between March and August 2020). Follow-up data collection occurred from October 2020 to April 2021.
Comorbid conditions prior to COVID-19 were identified using a database of information for all cases with suspected COVID-19 during their admission as inpatients at HCFMUSP (Busatto et al., Reference Busatto, de Araújo, Duarte, Levin, Guedes and Kallas2021), and we excluded patients with a previous diagnosis of dementia or end-stage cancer. Additional exclusion criteria were nosocomial COVID-19 infection, subjects living in nursing homes or long-term care facilities, and insufficient physical mobility to leave home. Finally, as re-infection by SARS-CoV-2 is possible, we enquired subjects about the emergence of symptoms and signs of infection (e.g. fever) during telephone calls on the day before the follow-up visit and again upon their arrival for the assessments. Any subjects with suspected re-infection had their research visit postponed, and they were referred to the infectious disease outpatient clinic at HCFMUSP dedicated to the diagnosis and management of acute COVID-19.
Online Supplementary material 1 (Methods S1) provides further details on the study period, procedures, and laboratory confirmation of COVID-19. A full description of general procedures for invitations and assessments was provided elsewhere (Busatto et al., Reference Busatto, de Araújo, Duarte, Levin, Guedes and Kallas2021).
Data from interviews, scales, and complementary examinations were captured and stored at real-time using web-based case report forms developed on a Research Electronic Data Capture (REDCap) system hosted at HCFMUSP (Harris et al., Reference Harris, Taylor, Thielke, Payne, Gonzalez and Conde2009).
Protocols were approved by the local ethics committee (numbers 4.270.242, 4.502.334, 4.524.031, 4.302.745, and 4.391.560). Informed consent was obtained from all participants or their proxy prior to study procedures.
Evaluation of psychiatric and cognitive symptoms at the follow-up
To assess symptoms of PTSD (Weathers, Litz, Herman, Huska, & Keane, Reference Weathers, Litz, Herman, Huska and Keane1993), anxiety and depression (Zigmond & Snalth, Reference Zigmond and Snalth1983), insomnia (Bastien, Vallières, & Morin, Reference Bastien, Vallières and Morin2001), and subjective memory impairment (Vale, Balieiro, & Silva-Filho, Reference Vale, Balieiro and Silva-Filho2012), multi-item scales were applied by specialized teams (Busatto et al., Reference Busatto, de Araújo, Duarte, Levin, Guedes and Kallas2021). Validated cutoffs were used to generate categorical ‘yes-no’ variables (see online Supplementary Methods S2 and Table S1 for details). Complaints of impaired concentration were documented using the Clinical Interview Schedule-Revised (CIS-R) (Lewis, Pelosi, Araya, & Dunn, Reference Lewis, Pelosi, Araya and Dunn1992). To identify potentially incident cases, we estimated the date of onset of symptoms using structured questions from the CIS-R; ‘no’ ratings were applied when symptoms were reported by patients as already present prior to SARS-CoV-2 infection.
Other self-reported symptoms at the follow-up
Additional symptoms previously highlighted as relevant to PASC (Lopez-Leon et al., Reference Lopez-Leon, Wegman-Ostrosky, Perelman, Sepulveda, Rebolledo, Cuapio and Villapol2021; Nasserie et al., Reference Nasserie, Hittle and Goodman2021) (i.e. fatigue, dyspnea, muscle/joint pain, and taste and/or olfaction changes) were also assessed using standardized scales (Bestall et al., Reference Bestall, Paul, Garrod, Garnham, Jones and Wedzicha1999; Boonstra, Preuper, Balk, & Stewart, Reference Boonstra, Preuper, Balk and Stewart2014; Brandão Neto et al., Reference Brandão Neto, Fornazieri, Dib, Di Francesco, Doty, Voegels and Pinna2021; Webster, Cella, & Yost, et al., Reference Webster, Cella and Yost2003) (see online Supplementary Table S1 for details). An adaptation of the WHO screening tool for neuroepidemiology investigations in low-and-middle-income countries (WHO, 1982) was applied by a team of neurologists (Busatto et al., Reference Busatto, de Araújo, Duarte, Levin, Guedes and Kallas2021) to assess neurological manifestations including headache, weakness, gait problems, episodes of loss of consciousness, and paresthesia. Finally, subjects were inquired about the presence of 14 persistent post-COVID symptoms covering additional organ-system domains (Busatto et al., Reference Busatto, de Araújo, Duarte, Levin, Guedes and Kallas2021) (see full list in online Supplementary Table S1).
The symptoms above, added to those cited in item 2.2, comprised a total of 30 items (online Supplementary Table S1).
Objective assessments of functioning and blood laboratory indices at the follow-up
Objective indices of organ system dysfunction were generated from the following assessments conducted during the same in-person visit (Busatto et al., Reference Busatto, de Araújo, Duarte, Levin, Guedes and Kallas2021): anthropometric measurements (including body mass index; BMI); 1 min sit-to-stand and handgrip strength tests (to measure physical disability) (Bohannon, Reference Bohannon2015; Litmanovich, Chung, Kirkbride, Kicska, & Kanne, Reference Litmanovich, Chung, Kirkbride, Kicska and Kanne2020; Strassmann et al., Reference Strassmann, Steurer-Stey, Lana, Zoller, Turk, Suter and Puhan2013; Vianna, Oliveira, & Araújo, Reference Vianna, Oliveira and Araújo2007); measurements of resting oxygen saturation, decreased oxygen saturation during the sit-to-stand test, a chest X-ray (Litmanovich et al., Reference Litmanovich, Chung, Kirkbride, Kicska and Kanne2020), and a spirometry test (Pereira, Sato, & Rodrigues, Reference Pereira, Sato and Rodrigues2007) (all addressing pulmonary dysfunction); and blood creatinine levels (to estimate glomerular filtration rate, as a measure of kidney function impairment) (Levey et al., Reference Levey, Stevens, Schmid, Zhang, Castro and Feldman2009).
Details regarding the criteria for rating objective indices of dysfunction as present are provided in online Methods S2 and Table S2 (Supplementary material 1). BMI values were used to detect the presence of clinically relevant weight change after COVID-19 (online Supplementary Table S2), and also to classify study subjects as obese (⩾30 kg/m2) or non-obese (⩽30 kg/m2) (Table 1).
Table 1. Baseline and hospitalization characteristics of the sample

s.d., standard deviation; IQR, interquartile range; ICU, intensive care unit.
a Six categories assessed in accordance to current criteria of the Associação Brasileira de Empresas de Pesquisa (ABEP, 2020).
b Based on self-reported mixed black and white races.
c Rated as present if the body mass index of individuals was equal or higher than 30 kg/m2, based on self-reported estimates of body weight prior to COVID-19 given by study subjects at the time of the follow-up assessments.
d WHO scale categories: 3–4, no continuous supplemental oxygen needed; 5, continuous supplemental oxygen only; 6, continuous positive airway pressure ventilation, bi-level positive airway pressure or high flow nasal oxygen; 7–9, invasive mechanical ventilation and/or extra-corporeal membrane oxygenation (ECMO). WHO Working Group on the Clinical Characterisation and Management of COVID-19 infection (2020).
Missing values for objective signs of dysfunction are provided in online Supplementary Table S3; most of the missing data were due to the fact that some of the assessment protocols were not ready for use at the starting date of data collection.
Methods for measuring C-reactive protein and D-dimer serum levels and information regarding missing data are provided in online Supplementary Methods S3.
Statistical analyses
Identifying symptom dimensions through exploratory factor analysis
Using the binary ratings for the 30 self-reported PASC symptoms outlined in items 2.2 and 2.3 above, we conducted an initial EFA including all symptoms. A total of 203 patients had to be excluded from this analysis (27.1% of the total sample) due to missing values for at least one symptom (see online Supplementary Table S4). Further methodological details and information on missing values are provided in online Methods S4 (Supplementary material 1).
Generation of a unidimensional latent variable through item response theory
Based on previous literature findings (Davis et al., Reference Davis, Assaf, McCorkell, Wei, Low, Re'em and Akrami2021; Evans et al., Reference Evans, McAuley, Harrison, Shikotra, Singapuri and Sereno2021), we expected that the initial EFA could indicate symptom data to fit an unidimensional latent trait space. Based on this prediction, we chose to generate a latent variable of PASC symptoms using two-parameter logistic IRT modeling (Hays et al., Reference Hays, Morales and Reise2000; Krueger et al., Reference Krueger, Nichol, Hicks, Markon, Patrick, Iacono and McGue2004), scaling the 30 symptoms of PASC along one single continuum. The IRT approach is suitable for scaling multiple items related to health outcomes along a unidimensional continuum of severity, and it has the additional advantage of being highly flexible for handling missing values (Krueger et al., Reference Krueger, Nichol, Hicks, Markon, Patrick, Iacono and McGue2004). In order to classify the 30 multi-organ symptoms included in the latent trait of PASC, we used the property of discrimination as our parameter of interest, referring to the increment in the probability a given symptom to be scored as present as the latent dimension score increased. Further information about the calculation of discrimination indices and IRT details is provided in online Supplementary Methods S5.
We run a sensitivity analysis assessing the properties of all symptoms included in the IRT model after leaving out subjects with comorbidities known to be themselves associated with CNS manifestations (see online Supplementary Methods S6). We also run additional sensitivity analyses including only one psychiatric or cognitive symptom at a time, to test the unique role of each self-reported psychiatric or cognitive symptom within the latent variable (see online Supplementary Methods S6).
Relationships between the latent variable of PASC symptoms, objective signs of organ system dysfunction and laboratory test results
Using two-parameter logistic regression models, we investigated significant associations between the latent variable described above and eight objective signs of organ system dysfunction (see online Supplementary Methods S7).
Additionally, we used a Differential Item Functioning (DIF) analysis approach to investigate relationships between elevated levels of blood tests (C-reactive protein and D-dimer) and the latent variable of PASC symptoms. The DIF approach allowed us to test whether significant relationships between abnormal test results and the latent dimension of PASC were due to different probabilities of endorsing any of the symptoms included in the latent variable (Saha et al., Reference Saha, Compton, Pulay, Stinson, Ruan, Smith and Grant2010; Penfield & Camilli, Reference Penfield, Camilli, Rao and Sinharay2007). Initially, we investigated the relationship between blood test results and the latent dimension of symptoms under investigation (i.e. syndrome level), with statistical significance assessed using logistic coefficients and 95% confidence intervals. If any significant relationship with the severity of a latent PASC trait was detected, we then tested differential effects at the level of the individual PASC symptoms (i.e. manifestation level) through Mantel–Haenszel tests (see details in online Supplementary Methods S7).
We also conducted sensitivity analyses (logistic regressions) investigating relationships between the latent variable of PASC and blood test results after leaving out subgroups of subjects stratified by sex, age, race, and socio-economic status, in order to verify if the above associations were substantially affected by the influence of any demographic variables on levels of C-reactive protein or D-dimer (see details in online Methods S8, Supplementary material 1). Finally, to take account of the influence of confounding comorbidities on blood laboratory markers, we conducted additional sensitivity analyses after stratifying the sample by the presence of the following diagnoses: obesity; diabetes; chronic lung or heart disease; and additional comorbidities known to affect the blood indices investigated (rheumatic disease, chronic liver or kidney disease, hematological disease; active cancer; and organ transplantation).
Results
Demographics and baseline characteristics
Figure 1 provides the flowchart of potential participants and patient selection (including data on refusals and exclusions). A total of 749 eligible individuals attended the in-person follow-up assessments. Table 1 provides their details regarding demographics, comorbidities, hospitalization events, and time between symptom onset and hospital admission. The mean duration of symptoms characterizing acute COVID-19 prior to hospitalization (as referred by patients and/or relatives upon admission) equaled 9.0 ± 6.5 days (Table 1).

Fig. 1. Potentially eligible subjects and final sample included in the follow-up investigation.
Online Supplementary Table S5 provides results of comparisons between the patients who attended the in-person assessments v. the remaining surviving individuals who did not participate (see details also in online Results S1, Supplementary material 1).
Exploratory factor analysis and latent variable modeling
The frequencies of each of the 30 symptoms in the overall sample are provided in online Supplementary Table S4.
One single factor was extracted in the initial EFA (eigenvalue = 4.65), explaining 64.8% of variance. All symptoms presented positive loadings for this factor; loadings >0.35 were found for 16 symptoms, with fatigue, insomnia, psychiatric, and cognitive complaints presenting the highest loadings (>0.5) (see online Supplementary Table S6). One additional factor approached the eigenvalue of one (0.99) and explained an additional amount of 13.8% of data variance, including the symptoms of weakness, gait impairment, and paresthesia, all with loadings above 0.45. Further details related to the EFA results are provided in online Results S2 (Supplementary material 1).
The unidimensional nature of the PASC symptom data, demonstrated by the above EFA, confirmed the validity of using the IRT approach to generate one single latent dimension of PASC. This analysis was conducted with all 749 subjects, taking advantage of the flexibility of the IRT approach for handling missing values. Discrimination properties for each of the 30 manifestations are displayed in Table 2. All 30 symptoms displayed discrimination coefficients from 0.2 to 2.3 (p ⩽ 0.005), with highest values for fatigue (2.3) and the six psychiatric and cognitive symptoms (1.5–2.1) (Table 2).
Table 2. Latent variable modeling of post-acute sequelae of SARS-CoV-2 infection (PASC) using Item Response Theory

CI, confidence intervals.
The sensitivity analysis assessing the properties of all symptoms included in the IRT model after leaving out subjects with comorbidities known to be themselves associated with CNS manifestations showed that all psychiatric and cognitive manifestations retained a high degree of discrimination (see online Supplementary Table S7).
In the six IRT analyses including only one psychiatric or cognitive symptom at a time, each of those variables persisted as highly discriminative within the resulting latent variable (see online Supplementary Table S8).
Significant associations between the latent variable of PASC symptoms, objective signs of organ system dysfunction, and laboratory test results
The three following objective signs of organ system dysfunction were significantly associated with the IRT-based latent variable of PASC symptoms: reduced BMI, abnormal sit-to-stand performance, and decreased handgrip strength (Table 3).
Table 3. Associations between the latent variable of PASC symptoms and objective signs of organ system dysfunction

CI, confidence intervals; BMI, body mass index; FVC, forced vital capacity; eGFR, estimated glomerular filtration rate.
a Difference between current BMI measurement and calculated BMI based on self-reported estimates of body weight prior to COVID-19.
b Oxygen saturation reduction of ⩾4 points during sit-to-stand test (not undertaken in subjects presenting resting pulse oximetry ratings lower than 90%).
c Estimated glomerular filtration rate calculated according to CKD EPI30.
Blood levels of C-reactive protein and D-dimer at the follow-up evaluation showed significant direct relationships with the latent PASC dimension, with highest significance for the former (Table 4). There were no differential effects on any of the individual manifestations included in the PASC trait either for C-reactive protein or D-dimer levels (Table 4).
Table 4. Significant findings from the differential item functioning analysis investigating blood laboratorial predictors of post-acute sequelae of SARS-CoV-2 infection (PASC)

CI, confidence interval; aOR, odds ratio adjusted for the relationship between the latent trait score and the other PASC symptoms; npc, not possible to calculate odds ratio and 95%CI.
Results of the sensitivity analyses investigating relationships between the latent variable of PASC and blood test results after exclusion of each subgroup stratified by demographic variables are detailed in online Supplementary Results S3; statistical significance was retained in all analyses involving levels of C-reactive protein levels, but not D-dimer. Associations between the latent trait of PASC symptoms and blood biomarkers in each subgroup stratified by demographic variables are detailed in online Supplementary Table S9.
Finally, the sensitivity analyses conducted after excluding subjects with obesity or other comorbidities known to affect values of C-reactive protein and D-dimer showed that associations of the latent variable of PASC remained statistically significant with those two blood biomarkers (online Supplementary Table S10).
Discussion
In this study, we applied statistical methods to document patterns of co-occurrence of several PASC symptoms in a relatively large sample of patients with moderate or severe COVID-19 evaluated 6–11 months after hospitalization, investigating how such symptom co-occurrence was associated with signs of physical disability and persistent inflammation.
We demonstrated that fatigue, insomnia, psychiatric, and cognitive complaints were the symptoms with highest loadings (>0.5) in the single factor extracted by EFA, which explained a large proportion of data variance (64.8%). Several other symptoms presented meaningful loadings in the same EFA factor, and each of them also showed significant discrimination (p ⩽ 0.005) along the latent variable of PASC generated using IRT (which allowed the recovery of cases with missing values). These findings confirm the notion that there is a significant degree of co-occurrence of multiple clinical manifestations in COVID-19 patients evaluated several months after hospitalization, suggesting that common underlying pathological mechanisms might influence on the persistence of different organ system symptoms in PASC.
The discriminative role of psychiatric and cognitive symptoms on the IRT-based latent variable of PASC remained high after the exclusion of cases with comorbid conditions themselves known to be associated with psychiatric and cognitive manifestations independently of COVID-19. Moreover, each of those CNS symptoms retained the same level of high discrimination when evaluated in isolation within the IRT model, together with fatigue. The fact that all other symptoms also displayed significant discrimination values in the same latent variable of PASC symptoms indicates that the likelihood of those additional symptoms to be referred by patients was higher in proportion to the presence of fatigue and psychiatric/cognitive complaints. These findings are consistent with the results of previous investigations that traced the trajectory of multiple individual symptoms of PASC over time, as these have indicated an incremental number of symptoms affecting multiple organ systems from approximately 10 weeks after acute COVID-19 onwards, with a prominence of CNS manifestations (Davis et al., Reference Davis, Assaf, McCorkell, Wei, Low, Re'em and Akrami2021).
The latent variable of PASC symptoms was directly related to three objective signs of persistent physical disability, namely low performance on the sit-to-stand test, impaired handgrip strength, and reduced BMI. Conversely, all other indices generated from objective tests of organ system functioning were unrelated to the latent trait of PASC symptoms, including the four signs of persistent respiratory dysfunction. The latter negative findings are consistent with previous observations that fatigue and ill health may be reported by patients several months after acute COVID-19 independently of whether objective signs of respiratory dysfunction are detected at the same timepoint (Townsend et al., Reference Townsend, Dowds, O'Brien, Sheill, Dyer, O'Kelly and Bannan2021a Reference Townsend, Fogarty, Dyer, Martin-Loeches, Bannan, Nadarajan and O'Donnell2021b). The pattern of high discrimination of fatigue, psychiatric, and cognitive symptoms in our latent trait of PASC symptoms, together with its specific association with physical disability, is consistent with the proposed similarity of PASC with myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS) (Bornstein et al., Reference Bornstein, Voit-Bak, Donate, Rodionov, Gainetdinov, Tselmin and Straube2021; Mackay, Reference Mackay2021; Perrin et al., Reference Perrin, Riste, Hann, Walther, Mukherjee and Heald2020; Phillips & Williams, Reference Phillips and Williams2021), a disorder thought to be often triggered by viral disease (Morris, Anderson, Galecki, Berk, & Maes, Reference Morris, Anderson, Galecki, Berk and Maes2013). Cognitive complaints are listed in the diagnostic criteria for ME/CFS (Committee on the Diagnostic Criteria for Myalgic Encephalomyelitis/Chronic Fatigue Syndrome, Board on the Health of Select Populations, & Institute of Medicine, 2015). The diagnosis of ME/CFS is also commonly associated with anxiety, depression, and weight loss (Afari & Buchwald, Reference Afari and Buchwald2003; Committee on the Diagnostic Criteria for Myalgic Encephalomyelitis/Chronic Fatigue Syndrome, Board on the Health of Select Populations, & Institute of Medicine, 2015).
The latent variable of PASC symptoms described herein was directly associated with persistently elevated blood levels of C-reactive protein at the follow-up assessment, supporting the proposed role for dysregulated inflammatory/immune mechanisms in the pathophysiology of PASC (Mackay, Reference Mackay2021; Perrin et al., Reference Perrin, Riste, Hann, Walther, Mukherjee and Heald2020). The latent trait of PASC symptoms was also significantly associated with elevated D-dimer levels, suggesting the persistence of a pro-thrombotic state possibly related to inflammatory mechanisms (Nalbandian et al., Reference Nalbandian, Sehgal, Gupta, Madhavan, McGroder, Stevens and Wan2021; Pasini et al., Reference Pasini, Corsetti, Romano, Scarabelli, Chen-Scarabelli, Saravolatz and Dioguardi2021; Townsend et al., Reference Townsend, Dowds, O'Brien, Sheill, Dyer, O'Kelly and Bannan2021a, Reference Townsend, Fogarty, Dyer, Martin-Loeches, Bannan, Nadarajan and O'Donnell2021b). These results extend the findings of Evans et al. (Reference Evans, McAuley, Harrison, Shikotra, Singapuri and Sereno2021), who reported elevated blood levels of C-reactive protein in proportion to the severity of a symptom cluster combining persistent psychiatric manifestations, dyspnea, fatigue, and poor physical performance in COVID-19 patients evaluated between 2 and 7 months after hospital discharge (Evans et al., Reference Evans, McAuley, Harrison, Shikotra, Singapuri and Sereno2021). Differently from that study, we conducted a set of sensitivity analyses which demonstrated that the association between PASC symptom severity and blood levels of either C-reactive protein or D-dimer remained significant in subgroups excluding patients with conditions known to affect levels of these biomarkers, including obesity, diabetes, and other pro-inflammatory comorbidities. Moreover, findings of the sensitivity analyses using subgroups stratified by demographic variables showed that the significant association between C-reactive protein levels and the latent trait of PASC symptoms was not determined by the influence of these factors (see details in online Discussion S1, Supplementary material 1).
The DIF analyses in our study showed no differential effects of either C-reactive protein or D-dimer levels on the individual symptoms included in the latent variable of PASC symptoms, suggesting that these biomarkers exert a uniform influence on psychiatric, cognitive, and all other symptoms along the overall PASC trait. Persistently elevated levels of these two indices and other inflammatory biomarkers were described in other studies that assessed subjects at variable time points after the onset of acute COVID-19 (Huang et al., Reference Huang, Yao, Gu, Wang, Ren, Wang and Cao2021b; Mandal et al., Reference Mandal, Barnett, Brill, Brown, Denneny, Hare and Hurst2021; Pasini et al., Reference Pasini, Corsetti, Romano, Scarabelli, Chen-Scarabelli, Saravolatz and Dioguardi2021; Townsend et al., Reference Townsend, Dowds, O'Brien, Sheill, Dyer, O'Kelly and Bannan2021a, Reference Townsend, Fogarty, Dyer, Martin-Loeches, Bannan, Nadarajan and O'Donnell2021b). However, associations with the severity of persistent clinical manifestations were inconsistently reported in those investigations, possibly due to modest sample sizes and the assessment of isolated PASC symptoms, rather than using a unifying latent variable as reported herein.
Peripheral biomarker findings suggestive of systemic inflammation were previously reported in association with symptoms of PTSD (Sankowski, Mader, & Valdés-Ferrer, Reference Sankowski, Mader and Valdés-Ferrer2015), mood disorders (Pariante, Reference Pariante2017), anxiety disorders (Michopoulos, Powers, Gillespie, Ressler, & Jovanovic, Reference Michopoulos, Powers, Gillespie, Ressler and Jovanovic2017), and ME/CFS (Montoya et al., Reference Montoya, Holmes, Anderson, Maecker, Rosenberg-Hasson, Valencia and Davis2017). A persistent state of systemic inflammation in PASC may lead to prominent psychiatric, cognitive, and physical dysfunctional manifestations via pro-inflammatory agents entering the CNS by circumventricular organs that have incomplete blood–brain barrier, or via abnormally permeable portions of blood–brain barrier damaged by cytokines (Mackay, Reference Mackay2021; Perrin et al., Reference Perrin, Riste, Hann, Walther, Mukherjee and Heald2020; Sankowski et al., Reference Sankowski, Mader and Valdés-Ferrer2015). Through those pathways, pro-inflammatory agents may affect limbic regions and the hypothalamus, reduce monoamine neurotransmission, and trigger microglia activation and neurotoxicity (Boldrini, Canoll, & Klein, Reference Boldrini, Canoll and Klein2021; Mackay, Reference Mackay2021; Perrin et al., Reference Perrin, Riste, Hann, Walther, Mukherjee and Heald2020).
The findings of this study highlight a multi-symptom dimension of PASC that is dominated by psychiatric manifestations and physical disability, and which may be related to persistent systemic inflammation. However, these results do not rule out the possible existence of additional clinical dimensions of relevance in PASC, combining other symptoms and objective signs of dysfunction. For instance, objectively assessed cognitive deficits (which were not addressed in the present study) may persist after COVID-19 (Becker et al., Reference Becker, Lin, Doernberg, Stone, Navis, Festa and Wisnivesky2021; Evans et al., Reference Evans, McAuley, Harrison, Shikotra, Singapuri and Sereno2021; Hampshire et al., Reference Hampshire, Trender, Chamberlain, Jolly, Grant, Patrick and Mehta2021) and have been suggested to be independent of psychiatric symptoms and physical disability (Evans et al., Reference Evans, McAuley, Harrison, Shikotra, Singapuri and Sereno2021); such objective cognitive deficits might be specifically associated with signs of pulmonary dysfunction (Sasannejad et al., Reference Sasannejad, Ely and Lahiri2019) or chemosensory symptoms (Pirker-Kees, Platho-Elwischger, Hafner, Redlich, & Baumgartner, Reference Pirker-Kees, Platho-Elwischger, Hafner, Redlich and Baumgartner2021), both of which were prevalent in our sample (online Supplementary Tables S3 and S4). It is also worth noting that our EFA analysis detected an additional trend pattern of symptom co-occurrence involving three neurological manifestations (weakness, gait impairment, and paresthesia) (WHO, 1982). This suggests that there may be different patterns of symptom co-occurrence in PASC involving neurological, psychiatric, and cognitive manifestations, which could emerge under the influence of separate risk factors and underlying brain mechanisms (Boldrini et al., Reference Boldrini, Canoll and Klein2021; Taquet et al., Reference Taquet, Geddes, Husain, Luciano and Harrison2021). These issues warrant further investigations in large samples of PASC sufferers.
One strength of the present study regards to the use of multi-item scales with validated cutoffs for the assessment of psychiatric manifestations ad other key symptoms (Zigmond & Snalth, Reference Zigmond and Snalth1983), rather than simpler ‘yes-no’ questions (Blomberg et al., Reference Blomberg, Mohn, Brokstad, Zhou, Linchausen, Hansen and Langeland2021; Huang et al., Reference Huang, Huang, Wang, Li, Ren, Gu and Cao2021a, Reference Huang, Yao, Gu, Wang, Ren, Wang and Cao2021b; Søraas et al., Reference Søraas, Bø, Kalleberg, Støer, Ellingjord-Dale and Landrø2021; Sudre et al., Reference Sudre, Murray, Varsavsky, Graham, Penfold, Bowyer and Steves2021; Writing Committee for the COMEBAC Study Group et al., Reference Morin, Savale, Pham, Colle, Figueiredo and Monnet2021). A study limitation is that subjects were asked at the time of the follow-up assessments about the presence of symptoms before COVID-19, which is prone to recall bias. Therefore, we could not confirm that the onset of symptoms occurred after COVID-19, since subjects did not undergo similar assessments prior to the SARS-CoV-2 infection. Other limitations include: the lack of a control group (Abel et al., Reference Abel, Carr, Ashcroft, Chalder, Chew-Graham, Hope and Pierce2021); the lack of chest computed tomography imaging and diffusion capacity testing to detect pulmonary involvement in COVID-19 with greater accuracy (Huang et al., Reference Huang, Wang, Li, Ren, Zhao, Hu and Cao2020, Reference Huang, Huang, Wang, Li, Ren, Gu and Cao2021a); the successful recruitment of only 43% of all patients potentially eligible for in-person visits; and the inclusion of a single hospital site, which may limit the generalizability of our findings. However, it should be noted that the size of our single-site sample was larger than that of most unicentric studies of COVID-19 patients that undertook in-person follow-up assessments to date (Mandal et al., Reference Mandal, Barnett, Brill, Brown, Denneny, Hare and Hurst2021; Mazza et al., Reference Mazza, De Lorenzo, Conte, Poletti, Vai, Bollettini and Benedetti2020; Menges et al., Reference Menges, Ballouz, Anagnostopoulos, Aschmann, Domenghino, Fehr and Puhan2021; Sonnweber et al., Reference Sonnweber, Boehm, Sahanic, Pizzini, Aichner, Sonnweber and Weiss2020; Townsend et al., Reference Townsend, Dyer, Jones, Dunne, Mooney, Gaffney and Conlon2020; Writing Committee for the COMEBAC Study Group et al., Reference Morin, Savale, Pham, Colle, Figueiredo and Monnet2021).
Conclusions and implications
The unidimensional profile of symptom clustering described in this study indicates that multiple persistent manifestations following moderate or severe COVID-19 may be united within one latent trait of PASC. This latent trait is dominated by fatigue, psychiatric, and cognitive symptoms, and is significantly associated with objective signs of physical disability and persistent systemic inflammation. Further longitudinal studies are needed to evaluate whether these PASC manifestations persist over longer follow-up periods. If confirmed in such additional investigations, our findings would reinforce the need for the development of healthcare services providing focused diagnostic assessments and evaluating the usefulness of tailored treatment programs for PASC sufferers (David, Reference David2021; Heightman et al., Reference Heightman, Prashar, Hillman, Marks, Livingston, Ridsdale and Banerjee2021). These services could identify patients that do not present long-term, organ-specific sequelae of COVID-19, but who may benefit from multidisciplinary interventions emphasizing physical rehabilitation and mental health-promoting strategies. The present findings also warrant investigations using panels of specific pro-inflammatory and anti-inflammatory markers in large-sized samples of PASC sufferers, in order to ascertain which inflammatory pathways may be most distinctly related to the symptom dimension highlighted herein; the identification of key biomarkers in such investigations may guide randomized clinical trials testing the efficacy of pharmacological interventions to reduce the burden of co-occurring physical disability and neuropsychiatric features of PASC. Finally, the findings of the present study highlight the relevance of full vaccination schemes against SARS-Cov-2, as these have been found to reduce not only the severity of acute COVID-19 but also the risk of PASC manifestations (Antonelli et al., Reference Antonelli, Penfold, Merino, Sudre, Molteni, Berry and Steves2022).
Supplementary material
The supplementary material for this article can be found at https://doi.org/10.1017/S0033291722001374.
Acknowledgements
We acknowledge the financial contribution to the study setup provided by donations from the general public under the HC-COMVIDA crowdfunding scheme (https://viralcure.org/c/hc) with funds managed by the Fundação Faculdade de Medicina. We are also grateful for the support from Patricia Manga Favaretto, Maria Cristina Coelho de Nadai, Vivian RB Saboya, and other members of the Diretoria Executiva dos Laboratórios de Investigação Médica at HCFMUSP in organizing the logistics for the follow-up assessments of COVID-19 subjects; the assistance of Rosa Maria Affonso Moyses in supervising clinical assessments, and Cristiana Rocca and Antonio de Pádua Serafim in supervising neuropsychological assessments; the help of Pedro Bacchi for the online programming of Clinical Interview Schedule-Revised scores; and the infrastructure support from the HCFMUSP COVID-19 task force (Antonio José Pereira, Rosemeire K Hangai, Danielle P Moraes, Renato Madrid Baldassare, Elizabeth de Faria, Gisele Pereira, Lucila Pedroso, Marcelo CA Ramos, Taciano Varro, and Vilson Cobello Junior) both during the baseline stage of in-hospital data collection and during the setting-up of the follow-up assessments. Finally, we thank the team led by Bruno Gualano in organizing the set-up for follow-up data collection regarding spirometry measurements.
Non-author contributors
HCFMUSP PASC Initiative. The members of the group (Claudia da Costa Leite, Cristiano Gomes, Guilherme Fonseca, Jorge Hallak, José Eduardo Krieger, Luis Yu, Luiz Henrique Martins Castro, Maria Cassia J Mendes Corrêa, Maria Elizabeth Rossi, Marilia Seelaender, Nelson Gouveia, Paulo A. Lotufo, Roger Chammas, Rossana Pulcinelli Francisco, Thais Mauad, Thiago Avelino-Silva, Wilson Jacob Filho) contributed in outlining the components of the multidisciplinary assessment protocols, conducting assessments of patients, and organizing the storage of demographic, clinical, and biomarker data.
Author contributions
G. F. B. had full access to all of the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis. Concept and design: G. F. B., A. L. d. A., J. M. C.-M., H. P. S. Acquisition, analysis, or interpretation of data: all authors. Drafting of the manuscript: G. F. B., A. L. d. A., J. M. C.-M., H. P. S., O. V. F., L. R. B. Critical revision of the manuscript for important intellectual content: all authors. Statistical analysis: J. M. C.-M. Obtained funding: E. G. K., C. R. R. d. C., G. F. B., V. G. R., G. G. C. Administrative, technical, or material support: A. L. d. A., A. J. d. S. D., N. M. S., K. R. d. S., M. I., M. L. G., M. V. Y. S., B. F. G., R. F. D.
Financial support
HCCOMVIDA crowdfunding scheme.
Conflict of interest
Dr Ferreira received speaker fees from Medtronic, outside of the submitted work. Dr Burdmann received speaker fees from AstraZeneca and Fresenius, outside of the submitted work. Dr Guedes holds stock in Fleury Ltd., a clinical analysis lab, which is not the provider of tests for this study. The other authors have no conflicts of interest to disclose.
Role of the funder/sponsor
The funders had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.
Group information
Non-author members of the HCFMUSP Long-COVID PASC Initiative (see full list in online Supplementary material 2).