key: cord-1036659-hw3cm3el authors: Torres Rives, Bárbara; Zúñiga Rosales, Yaíma; Mataran Valdés, Minerva; Roblejo Balbuena, Hilda; Martínez Téllez, Goitybell; Rodríguez Pérez, Jacqueline; Caridad Marín Padrón, Lilia; Rodríguez Pelier, Cira; Sotomayor Lugo, Francisco; Valdés Zayas, Anet; Carmenate Portilla, Tania; Sánchez Ramírez, Belinda; Carlos Silva Aycaguer, Luis; Ángel Portal-Miranda, José title: Assessment of changes in immune status linked to COVID-19 convalescent and its clinical severity in patients and uninfected exposed relatives date: 2022-04-12 journal: Immunobiology DOI: 10.1016/j.imbio.2022.152216 sha: 15648e989642d8004d1c1bc9e8db9b212c505365 doc_id: 1036659 cord_uid: hw3cm3el INTRODUCTION: The immune response during and after SARS-CoV-2 infection can be complex and heterogeneous, and it can be affected by the severity of the disease. It can also contribute to an unfavorable evolution and bring about short and long term effects. The aim of this study was to characterize the lymphocyte composition according to the severity of COVID-19, as well as its degree of relationship to the specific humoral response to SARS-CoV-2 in convalescents up to 106 days after the infection and in their exposed relatives. METHODS: An applied research was carried out with a cross-section analytical design, from March 11 to June 11, 2020 in Cuba. The sample consisted of 251 convalescents from COVID-19 over 18 years of age and 88 exposed controls who did not become ill. The B and T cell subpopulations, including memory T cells, as well as the relationship with the humoral immune response against SARS-CoV-2, were identified by flow cytometry and enzyme immunoassay. RESULTS: Convalescent patients, who evolved with severe forms, showed a decrease in frequency and a greater proportion of individuals with values ​​lower than the minimum normal range of B cells, CD3 + CD4 + cells and the CD4 + / CD8 + ratio, as well as a higher frequency and a greater proportion of individuals with values ​​above the normal maximum range of CD3 + CD8 + and NK cells. Convalescent patients with severe forms of COVID-19 that exhibited IgG / RBD titers ≥ 1/200 had a lower frequency of TEMRA CD8 + cells (p = 0.0128) and TEMRA CD4 + (p = 0.0068). IgG / RBD titers were positively correlated with the relative frequency of CD4 + CM T memory cells (r = 0.4352, p = 0.0018). CONCLUSIONS: The identified alterations of B and T lymphocytes suggest that convalescent patients with the severe disease could be vulnerable to infectious, autoimmune or autotinflammatory processes; therefore, these individuals need medical follow-up after recovering from the acute disease. Furthermore, the role of T cells CD4 + CM in the production of antibodies against SARS-CoV-2 is confirmed, and it is noted that the defect of memory T cells CD8 + TEMRA could contribute to the development of severe forms of COVID-19. The infection by the novel coronavirus SARS-CoV-2, (severe acute respiratory syndrome coronavirus 2) first reported by , is a severe health problem worldwide that has produced the death of more than 5 million people until October 29, 2021 (Lu et al., 2020; MINSAP, R., 2021a) . Within the immuno-pathology of the disease, it has not only been stated that the damage may be mediated by the virus itself, but that other factors are involved. These include a hyperinflammatory immune response, Kuri-Cervantes et al., 2020; , the exhaustion or dysfunction of T cells Zheng et al., 2020) , and the development of a cytokine storm Mathew et al., 2020) . These elements are responsible for the large clinical spectrum of COVID-19 (Kuri-Cervantes et al., 2020; Mathew et al., 2020) . A post COVID-19 syndrome having a significant clinical repercussion has been described (Greenhalgh et al.,2020) . It is expressed in that the dysfunction of the immune system may subsist for up to two years after viral infections (Wiedemann et al., 2020) , the magnitude of the antibody and T cell response may be diverse, discordant, and it may be influenced by the severity of COVID-19 (Wiedemann et al., 2020) . Several authors have shown that the response of the memory cells to SARS-CoV-2 will play an important role in the infection, pathogenesis and protection against the disease Sekine et al., 2020) . The study focus on the identification of the degree of alterations of the immune system that may favor the presence of sequels at the short-or long-term, in convalescents from SARS-CoV-2 (Greenhalgh et al.,2020; Shuwa et al., 2020) . The aim of this study was to characterize the lymphocyte composition according to the severity of COVID-19, as well as its degree of relationship to the specific humoral response to SARS-CoV-2 in convalescents up to 106 days after the infection and in their exposed relatives. defined as the presence of a negative RT-PCR to SARS-CoV-2, 14 days after the first negative RT-PCR (SARS-CoV -2) that granted the clinical discharge (MINSAP, R., 2020b , World Health Organization, 2021 ) (see Supplementary material). The inclusion and exclusion criteria of the subjects to study were exposed in Table 1 . The convalescent study time was from epidemiological discharge to blood collection for the evaluation with an average of 68 days (interquartile range of 55-77 days, minimum 20 days, and maximum 106 days) ( Table 2 ). The COVID-19 convalescents individuals were grouped according to the clinical severity of the infection by SARS-CoV-2: Asymptomatic convalescent (A): they experienced asymptomatic COVID-19; Moderate convalescent (M): they experienced mild or moderate COVID-19; Severe convalescent (S): they experienced a severe or critical form of COVID-19. The COVID-19 disease severity was classified, according to the guidelines for clinical management of COVID-19 of the World Health Organization (World Health Organization, 2021) and the guidelines of the national action protocol for COVID-19 (MINSAP, R., 2020b) as: Asymptomatic disease: when patients infected with SARS-CoV-2 (positive to RT-PCR), present no signs or symptoms of the disease (see symptoms associated to COVID-19 and Table S1 in Supplementary material); Mild disease: when patients (positive to RT-PCR) present symptoms of COVID-19, without evidence of viral pneumonia, hypoxia or others complications; Moderate disease: when patients infected with SARS-CoV-2 (positive to RT-PCR) present clinical signs and imagines of pneumonia, but without signs of severity and with oxygen saturation as measured by pulse oximetry (SpO2) ≥90% on room air; Severe disease: when patients (positive to RT-PCR) present clinical signs of pneumonia (fever, cough, dyspnea, fast breathing) plus one of the following: respiratory rate > 30 breaths/min; severe respiratory distress; or SpO2 < 90% on room air; Critical disease: when patients (positive to RT-PCR) present the Acute Respiratory Distress Symptom (ARDS), sepsis, septic shock, multiple organ dysfunction or other severe complications (Table S1 in Supplementary material) A group of 88 uninfected exposed individuals (Exposed), who were living in close contact to a first-degree relative (mother, father, children) having COVID-19, was included (Table 1 and Supplementary material). Blood samples for the study of the uninfected exposed individuals were collected the same day it was collected from their first-degree relative convalescent (minimum 42 days, and maximum 107 days from the diagnoses the convalescent relatives). A face-to-face interview was carried out for the collection of clinical, epidemiological and social data of the patients, given directly by the patients or by their legal tutors, in the case of intellectual disability. the red blood cell lysis without washing. To each volume of the conjugate, we added 100μL of blood, mixed for 3 s, and incubated it in a dark chamber for 15 min. at room temperature. Then we added 1mL of the lysis buffer VersaLyse TM (Beckman Coulter, France) and incubated it for 10 min under the same conditions as the previous step. Finally, we immediately proceeded to the acquisition of the sample by the cytometer. The acquisition of the data was carried out through the Kaluza Acquisition v1.0 software by which we obtained a minimum of 50000 total events. For the analysis and the report of the results we used Kaluza Analysis v1.5a. The absolute cell counts of the lymphocyte populations were performed through a dual-platform. We designed a manual and sequential window selection strategy with bi-parametric graphs ( Figure S1 ). For the reference values of the analyses of cellular sub-populations we used previous studies in the Cuban population (Kokuina et al, 2019) . We quantified the CD 19+ lymphocytes (B cells), T CD3+ lymphocytes, T CD3+CD4+ lymphocytes, T CD3+CD8+ lymphocytes and CD56+CD3-cells (NK cells). A polychromatic flow cytometry tube was used for peripheral lymphocyte immunotyping developed at the immunology laboratory of the National Medical Genetics Center (Zúñiga et al., 2020) . The monoclonal antibodies conjugated with fluorochromes from MACS MiltenyiBiotec (Germany) included anti-CD45 APC-Vio770 (Clone 5B1), anti-CD19 PE-Vio700 (Clone LT19), anti-CD3 FITC (Clone BW264/56), anti-CD4 PerCP-Vio700 (Clone M-T466), anti-CD8 APC (Clone BW135/80), anti-CD56 PE (Clone REA196), ( Figure S1 ). We identified the memory cells as: central memory: CM, CD45RA-CD27+, effector memory: EM, CD45RA-CD27−, terminally differentiated T effector cells (TEMRA, CD45RA+CD27−) and naïve cells (CD45RA+CD27+). The monoclonal antibodies used were CD8-PE-Cy7 (invitrogen, eBioscience, clone SK1), CD45 RA APC-eF780 (invitrogen, eBioscience, clone HI100), CD3/FITC (Clone BW264/56), Anti-CD27 APC, eBioscience, clone 0323, San Diego, CA, anti-CD127-PE, BD Pharmingen, Clone HIL-7R-M21, BD Biosciences ( Figure S1 ). We carried out the determination in the serum of total antibodies against SARS-CoV-2 through a double antigen sandwich-type ultra-immune-enzymatic assay (UMELISA ANTI SARS-CoV-2) that was standardized and validated at the Immuno-assay Center of Cuba (CIE, according to its Spanish acronym). The SARS-CoV-2 antigens were fragments from the spike protein (S) and the nucleocapsid (N) of SARS-CoV-2. (Supplementary material) At the Cuban Center of Molecular Immunology (CIM, according to its Spanish acronym), we quantified the total IgG specific RBD antibodies in the serum of patients using an enzymelinked immunosorbent assay (ELISA). The plates were coated with RBD-mFc and incubated with serial dilutions of serum samples, starting at 1:100. The experimental titers of IgG were determined (Supplementary material). See the supplementary material for particulars of the methods. The normal distribution of the quantitative variables was verified using the Shapiro-Wilk test. To describe the quantitative variables, the estimates were made through the median and interquartile ranges (IQR) or the mean and standard deviation, as appropriate. The 95% confidence intervals were also calculated. To assess the statistical significance of the association between qualitative variables and the comparison of proportions between each convalescent group and between those and the exposed group the Fisher's exact test was used. The Mann-Whitney U test was used for comparisons between two groups for the analysis of cell subpopulations by flow cytometry and anti-SARS-CoV-2 antibody levels. The correlation between the flow cytometric variables (cell subpopulations CD3, CD4, CD8, CD19, NK and memory cells CM, EM and TEMRA) and the IgG / RBD titers and total antibodies against SARS-CoV-2, was performed using the Spearman's rank correlation Using the IBM SPSS Statistics software (version 22), we carried out multivariate logistic regression analyses to evaluate the influence of age, severity and duration of the disease (we adjusted age) on the variables: CD19+, CD3+, CD3+CD4+, CD3+CD8+, NK. We also used the GraphPad Prism 7 (GraphPad Software, California, USA). We consider that there is statistical significance when p <0.05. The research was carried out under the compliance of the regulations of the Helsinki Declaration of 2013 (World Medical Association, 2013) . All cases participants in the research signed the informed consent before accepting their participation. This study is part of a research project approved by the Ethics and Research Committee of the National Medical Genetics Center, and by the advisory committee of the Ministry of Public Health of Cuba. The sample was of 251 individuals who had been ill with COVID-19, and in the group of exposed persons we included 88 first-degree relatives who were exposed to the virus and did not become ill ( Table 2 ). The clinical forms of COVID-19, from mild to moderate, were more frequent in the convalescents (48.6%), followed by individuals with asymptomatic forms of the disease (67%) (p<0.0001; 95% CI: 13.5-29.9) ( Table 2 ). Females were predominant (n=142, 56.6%, p=0.0385, 95% CI: 0.73-25.1) within all patients having COVID-19 ( Table 2 ). The median of age was higher (p<0.0001) in patients with severe forms of the disease compared to asymptomatic individuals (p<0.0001) and moderate (p<0.0001) ( Table 2 ). In convalescents that had severe forms of COVID-19, the time lapse between the diagnosis by RT-PCR of SARS-CoV-2 infection and the first negative PCR of the disease was slightly longer compared to patients with moderate COVID-19 symptoms (p= 0.0313) ( Table 2 ). The Convalescent study time (adopted for the study was from epidemiological discharge until day of blood sampling collection) was of 68 days (IQR: 55.0-77.0 days, minimum 20 days, and maximum 106 days) ( Table 2) . The absolute lymphocyte count was lower in convalescent people that had asymptomatic forms (A) and moderate forms (M) of the disease compared to the Exposed group that was not infected (Figure 1 ). The severe convalescent patients (S) had a higher proportion of individuals with an absolute lymphocyte count (14.5%) higher than the normal reference value compared to the A (3.0%, p=0.0184) and M (5.7%, p=0.0464) groups ( Figure 2 ). In the S group compared to the non-severe convalescent groups (A and M) and to those Exposed, we observed a lower relative and absolute frequency of CD19+, a lower median of relative frequency of CD3+CD4+ and of the ratio CD4+/CD8+, as well as an increase of the median of the relative frequency of CD3+CD8+ and NK (Figure 1 ). The S group also showed (compared to A and M) an increase of the median of the absolute frequency of CD3+CD8+ and NK (Figura1). The A and M groups had lower median of the absolute frequency of total lymphocytes, CD3+, CD3+CD4+ and CD3+CD8+ compared to those exposed ( Figure 1 ). In another analysis we observed a higher proportion of S (21.0%) that showed values below the minimum normal of the median of the relative frequency of CD19+ compared to those of A (4.5%, p=0.0047), M (5.7%, p=0.0018) and the exposed group (Exposed: 5.7%, p=0.0047). We also observed that the S group had a higher proportion of individuals (24.0%) with a median of the absolute frequency of CD19+, lower than the median reference value compared to the M (11.5%, p=0.0257) and to the exposed (5.7%, p=0.0047) groups. We also identified a higher proportion of S with a median of the relative frequency (8.1%, p=0.0182) and absolute frequency (6.5%, p=0.0345) of CD3+CD4+, lower than the median minimum normal reference value compared to the A (0%) group. At the same time, the median of the absolute frequencies of CD3+ and CD3+CD4+ were lower than the median minimum reference value in a larger proportion of the S group (8.1%, p=0.0047) compared to the exposed persons (0%), with the same percentages in both sub-populations ( Figure 2 ). We also observed in group S, compared to those Exposed who were not infected, a higher percentage of individuals with values higher than that of the normal established value of the median absolute frequency (severe: 21.0% vs exposed: 6.82%, p=0.0106) and the relative frequency of NK+ (severe: 8.1% vs Exposed: 1.14%, p=0.0325) ( Figure 2 ). Interestingly, the S group had a higher proportion of individuals with values above the maximum range established as the normal value of the absolute frequency of CD3+ (12.90%), CD3+CD8+ (12.9%) and of NK cells (21.%), compared to the A (CD3+: 1.49%, p=0.0109; CD3+CD8+: 0%, p=0.0025; NK+: 3,0 p=0.0014. Similarly, the M group had a higher proportion of individuals (compared to the S group) with absolute frequencies higher than the normal maximum value of CD3+: 3.3% p=0.0127; CD3+CD8+: 3.3%, p=0.0127 and NK+: 5.7%, p=0.0018 ( Figure 2 ). We analyzed T CD4+ and CD8+ memory cells (TEMRA: CD45RA+CD27+, CM: CD45RA−CD27+, EM: CD45RA−CD27-) and naive (virgin) cells from 85 (convalescent) at an average of 82 days (42 days as the minimum and 107 days as the maximum) after viral clearance for SARS-CoV-2 identified by RT-PCR, who had clinical and epidemiological discharge. We also studied 29 exposed individuals. The median of the relative frequency of total CD8+ memory cells (CM, EM, TEMRA) was significantly greater compared to the median of the T CD4+ memory cells (CM, EM, TEMRA) in all the convalescents studied (p<0.0001) ( Figure 3A ), and this occurred in a similar manner in each one of the groups according to the clinical forms (asymptomatic: p<0.0001, moderate: p=0.0213 and severe: p<0.0001) of COVID-19 ( Figure 3A ). In the analysis of total memory cells in all convalescents (p=0.0397) and severe cases (p=0.0030) studied, we identified a higher median of the relative frequency of total T CD8+ memory cells in comparison to the Exposed cases ( Figure 3A ). Nonetheless, we observed that in the entire group of convalescents studied, the medians of the frequencies of CM, EM, TEMRA and naive cells were similar between convalescents and exposed for T CD4+ and CD8+ memory cells (figure 3B). Within the subtypes of T CD8+ memory cells, the highest median of the relative frequency corresponded to TEMRA CD8+ (34.0%), followed by memory cells CM T CD8+ (24.8%), although no statistical differences were found between the median of the frequencies (p=0.1889, 95% CI: -4.47 to 22.4) ( Figure 3 ). No statistical differences were observed between the median of the relative frequencies in the subpopulations of T CD8+ memory cells (TEMRA: CD45RA+CD27+, CM: CD45RA−CD27+, EM: CD45RA-CD27-) from the convalescent and Exposed groups ( Figure 3) . Nevertheless, the median of the frequency of virgin cells (CD45RA+CD27+) in the S group was lower compared to that of the Exposed group (9.5% vs 11.6%; p=0.0379) ( Figure 3A ). The median of the relative frequency of the total T CD4+ memory cells in relation to the Exposed group, was similar in the total convalescent group (p=0.3263) and in those presenting asymptomatic (p=0.7282), moderate (p=0.2666) and severe (p=0.4101) forms of the disease ( Figure 3B ). The median of the relative frequency of the total T CD4+ memory cells were also similar between the individuals with different degrees of COVID-19 severity ( Figure 3A ). The phenotype of T CD4+ memory cells that prevailed were the cells specialized in central memory (51.1 %), followed by the EM T CD4+ (21.3%) (p<0.0001, 95% CI: 16.6 to 43.6504) ( Figure 3A ). The medians of the relative frequencies of T CD4+ CM, T CD4+ EM memory cells and naive cells from convalescent individuals of all the clinical forms analyzed in this study, were similar to those exposed. However, the median of the relative frequency of the TEMRA T CD4+ cells was higher in the M group (6.9%) compared to the Exposed group (Exposed: 3.9%, p=0.0158) and to the S group (4.3%, p=0.0182) ( Figure 3A ). Using multivariate analysis, we detected that the risk of presenting an increase in the absolute frequency of NK+ in convalescents was 4.0 times greater in those presenting the severe forms of the disease (severe) (adjusted OR: 34.0; 95% IC: 1,47-10,8; p=0.007). Furthermore, it was 1.06 times higher for each day the acute disease was extended (duration of the disease) (adjusted OR: 1.06; 95% IC: 1.0-1.11; p=0.036). We also found that the convalescents that progressed with severe forms of the disease showed a trend towards an increase in the relative frequency of the NK+ cells (adjusted OR: 2.32; 95% IC: 0.58-9.28; p=0.2320) and a rise in the relative frequency (adjusted OR: 1.49; 95% IC: 0.37-6.08; p=0.5790), and absolute frequency of CD3+CD8+ (adjusted OR: 3.23; 95% IC: 0.95-11.0; p=0.0610) ( Figure 4 , Table 1S ). Convalescents, of 60 or more years of age, had a higher risk of decreasing the relative frequency (adjusted OR: 2.70; 95% IC: 1.07-6.78; p=0.0390) and absolute frequency (adjusted OR: 2.84; 95% IC: 1.34-6.03; p=0.007) of CD19+ ( Figure 4 , Table 1S ). Age was positively correlated to the total number of CD8+ memory cells in convalescent patients (r=0.5748. p<0.0001, Spearman correlation) and a similar behavior was found in the exposed persons (r= -0.5688, p=0.0013, Spearman correlation). Furthermore, age was negatively correlated with T CD8+ naive cells in the S and M groups (r= -0.5213, p<0.0001, Spearman correlation, data not shown) and in those exposed +(r= -0.5491, p=0.0020, Spearman correlation, data not shown). The duration of the disease was found to be positively correlated with the T CD8+ EM cells (r=0.3480, p=0.0192 Spearman correlation, data not shown). 3.5 Correlation of B, T, NK, memory and naive cells with the response of specific antibodies to SARS-CoV-2 in convalescent Cubans according to the clinical forms of COVID-19 IgG/RBD titers are positively correlated with the relative frequency of the CD8+ CM memory cells (r=0.3132, p=0.0320) in the S group ( Figure 5-A) . TEMRA CD8+ showed a tendency to correlate negatively, but there were no statistical differences (r= -0.2405, p=0.1035) ( Figure 5 In line with the previous results of IgG/RBD, we identified a lower frequency of TEMRA T CD8+ cells (53.3% vs 27,3 %, p=0.0137) in the peripheral blood of the positive S group for total antibodies against fragments of the N and S protein of SARS-CoV-2 (compared to those that were negative to this antitotal antibody) ( Figure 5 -E). Among the individuals with a presence or absence of antitotal antibodies to SARS-CoV-2, we identified similar relative frequencies of CD19+, CD3+, CD3+CD4+, CD3+CD8+, NK+ and the ratio CD4/CD8+ in all ranges of severity studied (results not shown). The alterations of the B, T and NK cells, as well as the more and more frequent presence of signs and symptoms, have been reported in COVID-19 convalescents, (Greenhalgh et al.,2020; Shuwa et al., 2020) . There are, however, discrepancies in the magnitude and the protective or pathogenic role of the following immune response to SARS-CoV-2 infection Chen et al., 2020; Mathew et al., 2020) . As observed by other authors, in this study we did not observe lymphopenia in COVID-19 convalescents Rodriguez et al., 2020) . A higher proportion of severe convalescents with values above the reference range of total lymphocytes could be influenced by the larger number of T CD8+ lymphocytes identified as a response to clear the persistence and the greater antigenic magnitude of SARS-CoV-2 and because of the use of immunomodulators according to the national protocols (MINSAP, R., 2020a; Hernández et al., 2021), since it has been reported that the administration of biological therapies in COVID-19 patients, produces an increase of circulating lymphocytes (Giamarellos et al., 2020; Hernández et al., 2021) . The normalization of CD19+ cells in COVID-19 convalescents is reported in the literature Sherina Sherina) . A decrease of CD19+ cells has also been observed in patients with severe forms of the disease, Deng et al., 2020) as those found in this study. The decrease of CD19+ lymphocytes could suggest that sub-populations of this compartment, such as B regulator cells with anti-inflammatory functions, are low, and as a consequence the convalescents show a delay in their complete recovery, and are vulnerable to auto-immune, auto-inflammatory and infectious processes (Mauri et al., 2017) . Several mechanisms may explain the lymphopenia of CD19+ and CD3+CD4+ observed in this study, such as the presence of a greater viral load and exposure time to SARS-CoV-2 in more severely ill patients. This leads to an increase in the direct action of SARS-CoV-2 on the cells, and the damage mediated by the immune system, the sequester of cells from the lung or peripheral lymphoid organs induced by the cytokine storm, apoptosis and the suppression of the bone marrow and the thymus (Wen et al., 2020) . The literature also reports, however, that convalescent individuals have similar values of CD3+CD4+ to those of individuals who did not become ill (Shuwa et al., 2021; Townsend et al., 2021; Wen et al., 2020) . The absence of asymptomatic convalescent patients with a decrease of the T CD4+ lymphocyte values and a higher frequency of these cells in relation to the convalescents with symptoms (moderate and severe), may correspond with the evidence that asymptomatic individuals have a higher secretion of INF -ᵞ e IL-12, as well as a proportional secretion of IL-10 and of pro-inflammatory cytokine (IL-6, TNF-α e IL-1β). This fact suggested that the asymptomatic patients have the ability of developing a less intense inflammatory process, but their antiviral response is protective, efficient, balanced and specific, so that it protects the host and does not produce any apparent pathology (Le Bert et al., 2021) . The discrete increase in T CD3+CD8+ lymphocytes in convalescents from the severe illness supports the role of these cells when facing a greater antigenic exposure and an exaggerated immune response to achieve effective viral clearance (Wen et al., 2020; Thieme et al., 2020) . These results agree with the expansion of T CD8+ lymphocytes in convalescents reported by other authors (Shuwa et al., 2021; Wen et al., 2020) . However, the recovery from lymphopenia of T CD8+ characteristic of the acute phase of the disease is also reported. This has led to the idea that the virus produces this alteration and that the effective anti-viral therapy leads to the recovery of T CD8+ cells (Zheng et al., 2020) . The increase of T CD8+ lymphocytes in convalescents may have implications in the development of later infections or the perpetuation of inflammatory processes, depending on the capacity of the cytokine secretion of these cells (Shuwa et al., 2021) . Similarly, other studies have reported an increase in NK cells during the convalescent stages (Rodríguez et al., 2020; Wen et al., 2020) . It has been reported that an effective therapy for SARS-CoV2 is accompanied by an increase of NK cells (Zheng et al., 2020) . The increase in NK+ cells is considered to be a valuable biomarker for monitoring the progression of the acute phase of the disease toward recovery stages in severe patients (Rodríguez et al., 2020) . In contrast, values of NK+ cells are also reported to be similar between convalescents and persons who are not infected by SARS-CoV-2 (Townsend et al., 2021; Liu et al., 2021) . This research also showed that an increase in the frequency of these cells is associated to the severity and duration of the disease. This supports the antiviral role and the participation in the immunopathology of these cells on the severity of the disease and on the stages of inflammation that may persist (Townsend et al., 2021; Fox et al., 2012; Market et al., 2020) . The immunologic memory is considered to be of great importance in preventing the recurrence of severe forms of COVID-19 in individuals who are seronegative to SARS-CoV-2, whether they are exposed or not (Sekine et al., 2020) . It is possible that a small part of the population infected with SARS-CoV-2, having a poor immunological memory, will be susceptible to reinfection shortly after recovering from the acute process . The similarity of the memory and naive cells among all convalescents and uninfected persons identified in this study has been reported in the literature . The prevalence of central memory cells in the compartment of T CD4+ memory in the present study, as well as the identification in severe convalescents of a positive correlation between the T CD4+ CM cells and the RBD titers, corresponds to the capacity of these cells of extravasation and migration to secondary tissues. These show a high proliferative capacity and a low dependence on co-stimulators, thus favoring the formation of specific antibodies against SARS-CoV-2, as observed (Wen et al., 2020; Peng et al., 2020; Weiskopf et al., 2020; Lugli et al., 2013; Neidleman et al., 2020; Mahnke et al., 2013) . Consistent with other reports from the literature, Wen et al., 2020; Dan et al., 2020; Peng et al., 2020; Weiskopf et al., 2020; Neidleman et al., 2020; Yang et al., 2007) , in this study there was a predominance of TEMRA cells in the sub-set of T CD8+ memory cells, which agrees with their function as potent producers of interferon -γ and perforins that mediate in the specific cytotoxicity of the antigen. This makes them highly important in viral infections (Sallusto et al.,2004) . The high frequency of TEMRA memory cells in convalescent individuals presenting moderate forms of the disease (compared to the severe forms), supports the protective role of these cells in the development of severe forms of COVID-19 and endorses the substantial role of T-cell immunity in SARS-CoV-2, (Sekine et al., 2020; Dan et al., 2020; Peng et al., 2020; Le Bert et al., 2020) , and other viral infections (Sridhar et al., 2013) . Hence, the formation of T memory cells has been associated to recovery from COVID-19, and it was reported that this response could predict severity and it could become a marker associated to the loss of the effectiveness of the anti-viral response (Odak et al., 2020) . Previous studies have demonstrated that the severity of the disease is inversely correlated with the immunity of T-cells (Ni et al., 2020) , and that the deficient T-cell response prevents the positive action of the immune system against SARS-CoV-2 (Odak et al., 2020; Wang et al., 2020) . Consistent with these reports, we observed that the convalescents who became ill with severe COVID-19 showed a lower frequency of TEMRA cells associated to higher titers of IgG/RBD and to the presence of total antibodies against fragments of the N and S proteins of SARS-CoV-2 (TEMRA CD8+). We also found a negative correlation between TEMRA CD4+ and the IgG/RBD titers. On the other hand, the similarity in the immune response (T CD3+CD4+, CD3+CD8+, NK and CD19+) in individuals who were seropositive or not to SARS-CoV-2, suggests that the state of protection evaluated through the detection of antibodies against SARS-CoV-2, may be underrated (Sekine et al., 2020) . Although we did not determine the specific SARS-CoV-2 T cells, we did observe that the frequency of the memory and naive cells in all convalescents was similar to that of the exposed individuals. It has been found a response of specific SARS-CoV-2 T cells in persons having close contact with COVID-19 patients, in which no positive RT-PCR was detected, nor the presence of antibodies anti-SARS-CoV-2; therefore, it suggested that there were no infections due to the limited exposure of the persons to viral particles, or the short exposure time (Sekine et al., 2020; Wang et al., 2021) . The decrease of T CD8+ naive cells in patients having severe forms of COVID-19 may be the result of the mobilization of effector cells, as the cytotoxic lymphocytes, in order to eradicate viral infection (Odak et al., 2020; Thieme et al., 2020) . Other authors observed a decrease in T CD8+ naive cells in individuals having a mild or moderate SARS-CoV-2 infection (Odak et al., 2020) . It should be mentioned that the group of severely ill patients were older, and that the decrease of naive T cells has been associated to immunosenescence, as a consequence of aging (Fulop et al., 2011; Verdecia et al., 2013; Saavedra et al., 2017) . It must be considered that older age could produce a hyper-inflammatory state, starting with the fact that during aging and immunosenescence there are changes occurring in the immune system (Sauce et al., 2011) . These include the effect on several cellular compartments, (Verdecia et al., 2013) deregulation of cytokine secretion and association to a state of "inflammaging" (low-degree chronic and sterile inflammation during aging), producing an increase in the frequency of infectious, neurodegenerative, and cardiovascular diseases and cancer (Fulop et al., 2018) . One limitation of this study was the fact that we did not analyze any specific cells for SARS-CoV-2. Also, nor we did not carry out a longitudinal study, with patients through the time (samples were collected between 14-and 106-days post-infection). The evaluation of specific responses to SARS-CoV-2 in a longitudinal study could have provided a more comprehensive view of the dynamics of the immune response during COVID-19 and in the convalescent study time, as reported by other authors (Rodriguez et al, 2020; Liu et al, 2021; Wen et al., 2021) . It would have been interesting to have included a group of patients experiencing acute COVID-19, as well as to correlate the results of all groups with the clinical condition present in the convalescence period. (Shuwa et al., 2021; Sekine et al., 2020; Peng et al, 2020) . Despite these limitations, it was observed that the identified alterations in the immune response in convalescents is influenced by the severity of the disease, similar to what other authors affirm (Peng et al, 2020; Shuwa et al., 2021; Sekine et al., 2020; Wen et al., 2021) . Besides, these patients are susceptible to subsequent complications mediated by the immune system (Shuwa et al., 2021) . The study of T cells specific to SARS-CoV-2 would have allowed the identification of the presence of cross reaction T cells against SARS-CoV-2 and other coronaviruses, which could have possibly been an element to explain the fact that uninfected exposed individual did not get sick Le Bert et al., 2021 . The immune status of COVID-19 convalescents is influenced by the severity of the disease. The alterations of the lymphocytes CD19+, CD8+, NK cells and of the specific antibodies against SARS-CoV-2 in severe convalescents suggest that these patients could be vulnerable to infectious, autoimmune or autotinflammatory processes. These findings could be associated with a more unfavorable recovery and the instauration of new sequels of the disease, thereby needing close medical supervision. The alterations in the effector memory cells may be related to the evolution towards severe forms of the disease. Zheng, M., Gao, Y., Wang, G., Song, G., Liu, S., Sun, D., Xu, Y., Tian, Z., 2020. Functional exhaustion of antiviral lymphocytes in COVID-19 patients. Cell Mol Immunol 17, 533-535. https://doi.org/10.1038/s41423-020-0402-2 Zúñiga Rosales, Y., Villegas Valverde, C., Torres Rives, B., Hernández Reyes, E., 2020. Diseño y optimización de un tubo policromático de citometría de flujo para inmunofenotipo linfocitario periférico. Revista Cubana de Hematología, Inmunología y Hemoterapia 36(3). URL http://www.revhematologia.sld.cu/index.php/hih/article/view/1151/1002 (accessed 7.20.21). Clinical discharge criteria Close contact Symptoms associated with COVID-19 Table S1 : COVID-19 disease severity Obtaining biological samples Flow Cytometry Qualitative determination of total antibodies anti-SARS-Cov-2 in the serum Detection of antibodies anti-RBD Figure. S1: Window strategy to identify immunophenotypes of lymphocytes in peripheral blood using flow cytometry. Table S2 : Influence of age ≥60 years old, the severe forms, and the duration of the disease on the alterations of the cellular sub-populations in COVID-19 convalescents. The patient's condition is stable and afebrile for more than 3 days without antipyretic medication, regular breathing and normal or markedly improved respiratory rate, clear consciousness, unaffected speech and normal diet, lung images show significant improvement without signs of organic dysfunction, two consecutive RT-PCR tests for SARS-CoV-2 with at least 24 hours between them. (MINSAP, R., 2020b; World Health organization, 2021). Defined as a person who lived and had been in contact or within a close distance (<1.0 m) to an individual(s) with COVID-19 in a confined space for > 24 h , from 2 days before and up to 14 days after the onset of symptoms. (MINSAP, R., 2020b; World Health Organization, 2020a). The symptoms and signs of COVID-19 can be variable. Fever, cough, fatigue, anorexia, dyspnea, myalgias. Other non-specific symptoms, such as sore throat, nasal congestion, headache, diarrhea, nausea and vomiting, loss of smell (anosmia) or loss of taste (ageusia). Elderly people and immunocompromised patients may present with non-classical symptoms such as fatigue, reduced alertness, reduced mobility, diarrhea, loss of appetite, neurological symptoms and absence of fever. (MINSAP, R., 2020b; World Health organization, 2021). Table S1 : COVID-19 disease severity criteria (MINSAP, R., 2020b; World Health organization, 2021) . Patients with laboratory-confirmed (RT-PCR to SARS-CoV-2) but who have no symptoms, which are consistent with COVID-19. See Symptoms associated with COVID-19 in Supplementary material. Patients with laboratory-confirmed (RT-PCR to SARS-CoV-2) and symptoms of COVID-19 (fever, cough, fatigue, anorexia, dyspnea, myalgias. Other nonspecific symptoms, such as sore throat, nasal congestion, headache, diarrhea, nausea and vomiting, loss of smell or loss of taste), without evidence of viral pneumonia, hipoxia or others complications. See Symptoms associated with COVID-19 in Supplementary material. Patients with a diagnose of infection by SARS-CoV-2 with clinical signs or Chest imaging (radiograph, computed tomography, or lung ultrasound) of pneumonia (fever, cough, dyspnea, fast breathing), but without signs of severe neumonia, with oxygen saturation as measured by pulse oximetry (SpO2) ≥90% on room air. Patients (RT-PCR to SARS-CoV-2) who experienced severe pneumonia (fever, cough, dyspnea, fast breathing) plus one of the following: respiratory rate > 30 breaths/min; severe respiratory distress; or (SpO2) < 90% on room air. Patients who experienced one of the following symptoms:  Acute respiratory distress syndrome (ARDS): pneumonia a week after beginning or worsening respiratory symptoms. Chest imaging: (radiograph, computed tomography: CT scan, or lung ultrasound): bilateral opacities, not fully explained by volume overload, lobar or lung collapse, or nodules. Pulmonary infiltrates: respiratory failure not fully explained by cardiac failure or fluid overload. • Mild ARDS: 200 mmHg < partial pressure arterial oxygen (PaO2)/ fraction of inspired oxygen (FiO2a) ≤ 300 mmHg (with positive end-expiratory pressure: PEEP or continuous positive airway pressure: CPAP ≥ 5 cmH2O). When PaO2 is not available, SpO2/FiO2 ≤ 315 suggests ARDS. • Moderate ARDS: 100 mmHg < PaO2/FiO2 ≤ 200 mmHg (with PEEP ≥ 5 cmH2O). When PaO2 is not available, SpO2/FiO2 ≤ 315 suggests ARDS • Severe ARDS: PaO2/FiO2 ≤ 100 mmHg (with PEEP ≥ 5 cmH2O). When PaO2 is not available, SpO2/FiO2 ≤ 315 suggests ARDS.  Sepsis: acute life-threatening organ dysfunction caused by a dysregulated host response have a suspicion or verified infection: fast breathing, low oxygen saturation, reduced urine output, altered mental status, fast heart rate, weak pulse, cold extremities or low blood pressure, laboratory evidence of coagulopathy, thrombocytopenia, acidosis, high lactate, or hyperbilirubinemia.  Septic shock: persistent hypotension despite volume resuscitation, requiring vasopressors to maintain mean arterial pressure (MAP) ≥ 65 mmHg and serum lactate level > 2 mmol/L. Complications in all levels can showed up and they should be considered as severe forms of the disease, and decompensated chronic diseases as well We took 7 mL of peripheral venous blood of the convalescents and controls through venipuncture at their forearm by specialized personnel. This was deposited (4mL) in a Vacutainer tube with ethylendiamine tetraacetic acid, in order to study the cellular immunotypes. The rest of the blood was used to obtain the serum for the quantification of specific antibodies against SARS-CoV-2. The serum was stored at −80 ˚C for its later analysis. Blood samples were collected from convalescents from epidemiological discharge until (68 days with an interquartile range of 55-77 days). Using the percentages of the cellular sub-populations obtained, we calculated the absolute values of each sub-population (cells/µL). The total number of lymphocytes /µL was obtained using the Mindray BC6800 hematologic counter. The following formula was applied: Cells/µL = percentage of the population of interest × (total lymphocytes/µL) 100 Quality control was performed with Flow-CheckTM Pro Fluorospheres, (Beckman Coulter, France) to verify the stability of the optic and fluid systems. Data acquisition was performed with the Kaluza Acquisition v1.0 software, and a minimum of 50000 total events was obtained. For the analysis and report of results, we used the Kaluza Analysis v1.5a software. The quantification of lymphocyte populations was carried out through a dual-platform. We designed a strategy of manual and sequential selection windows with bi-parametric graphs (Fig. S1 ). Graphs were created for the control of positive events of global fluorescence vs side scatter (SS), reflecting structural complexity. Using a discrimination graph for coinciding events, we created a gate for the elimination of detritus, so that we could finally obtain the window of the selection of lymphocytes of the graph of CD45+ against SS. From the CD45+ window, we constructed the panels (graphs) to quantify the populations and sub-populations of interest (Fig. S1 ). At the Immunoassays Center of Cuba (CIE, according to its Spanish acronym), we carried out the determination of total antibodies against SARS-CoV-2 in the serum, through a sandwich type ultra-immunoenzymatic double antigen assay (UMELISA ANTI SARS-CoV-2). For the solid phase, we used ultra-micro ELISA plates coated with the spike (S) and nucleocapsid (N) proteins of SARS-CoV-2. We added the controls and 10 µL of the diluted samples 1:2 with the buffer solution and the sheep serum. This was incubated for 30 minutes at 37 °C in a humid chamber. It was washed 4 times in a SUMA technology washer (SUMA, Cuba) adding 30 μL of the buffer solution in each well. We added 10 μL of the biotinilated antigen solution (CIE, Cuba) and incubated it for 30 minutes at 37 ºC in a humid chamber, and then we again made 4 washings with the buffer solution. From the flask of the Streptavidin Conjugate / alkaline phosphatase (CIE, Cuba) we added 10 μL of the mixture and incubated it for 30 minutes at 37 ºC in a humid chamber. It was washed 4 times and the 4-methylumbeliferir phosphate fluorigenic substrate was added and incubated for 30 minutes in a humid chamber at between 20 -25 ºC. Therefore, we guaranteed that the fluorescent signal of the positive control was between 60 and 180 units. The reading of the intensity of the fluorescence produced was carried out using a reader of the SUMA series (Havana, Cuba), with the UMELISA ANTI SARS-CoV-2 software. The samples are considered positive (reactive) when complying with the quality control parameters in the control specimens and obtaining a value of, or higher than, 2.9 as a result of the division of the fluorescence of the sample by the median of the negative control. The determination of Total RBD-specific IgG was made at the Molecular Immunology Center of Cuba (CIM). The RBD-specific IgG ELISA assay was performed in microtiter plates (Maxisorp, NUNC) that were coated with 0.5 µg/well of RBD-mFc for 1h at 37°C. Plates were blocked with 150 µL/well of 4% skim milk in PBS with Tween-20 0.05% (PBST) for 30 min, at 37°C. Next, the plates were incubated with serial dilutions of serum samples (100 µL/well, starting at 1:100) overnight at 4°C. After washing with PBST, plates were incubated with 100 µL/well of a peroxidase conjugated anti-human IgG monoclonal antibody (1:5 000) for 1h at room temperature (RT). Then, plates were washed and the reaction was visualized by the addition of o-phenylenediamine (0.5 mg/ml) in peroxidase substrate solution, and was stopped with H2SO4 (3 M). The absorbance at 490 nm was measured using a microwell system reader (Organon Teknica, Salzburg, Austria). Experimental antibody titers were considered as the inverse of the highest serum dilution giving optical density (OD) values four-fold the value of the negative control serum. ; the proportion of convalescents with severe clinical forms (S) (the red bars); those exposed (Exp) (gray bars). Where n is the total number of convalescents and those exposed with alterations in comparison to the total of each group (total number of convalescents with asymptomatic disease forms, 67; total number of convalescents with moderate disease forms, 122; total number of convalescents with severe disease forms, 62; total number of exposed individuals, 88). The cellular sub-populations are identified through flow cytometry. The 95% confidence intervals are shown in parenthesis, with which we identified the statistical significance. The between-groups significance was calculated through the proportions comparison; the statistical significance was established for p<0.05 and it was represented as: ****p<0.0001, ***p<0.001, **p<0.01, *p<0.05. Uninfected Exposed individuals -Age > 18 years -Of both sexes -Who accepted their participation in the study -That had been sick with COVID-19 confirmed by RT-PCR for SARS-CoV-2. -That were living in close contact** to a first-degree relative (mother, father, children) having COVID-19. -That had Epidemiological discharge.* -They did not test positive to two RT-PCR for SARS-CoV-2 done in the quarantine period (14 days) since their convalescent relative was diagnosed with COVID-19. Inclusion criteria -These exposed individuals were also negative to specific antibodies against SARS-CoV-2. -The individuals who were not Cuban residents -The individuals who were residing outside the health area at the time of inclusion in the study. -Deceased *Epidemiological discharge: Defined as the presence of a negative RT-PCR to (SARS-CoV -2), 14 days after the first negative RT-PCR to (SARS-CoV-2) that granted the clinical discharge (Ministry of Public Health, 2020, World Health Organization, 2021). **Close contact: Defined as a person who lived and had been in contact or within a close distance (<1.0 m) to an individual(s) with COVID-19 in a confined space for > 24 h , from 2 days before and up to 14 days after the onset of symptoms. (Ministry of Public Health, 2020a, World Health Organization, 2020). T cell responses in patients with COVID-19 Immunological memory to SARS-CoV-2 assessed for up to eight months after infection Dynamic changes in peripheral blood lymphocyte subsets in adult patients with COVID-19 Severe Pandemic H1N1 2009 Infection Is Associated with Transient NK and T Deficiency and Aberrant CD8 Responses The integration of inflammaging in age-related diseases Immunosenescence and gender: a study in healthy Cubans Complex Immune Dysregulation in COVID-19 Patients with Severe Respiratory Failure Management of post-acute covid-19 in primary care Targets of T Cell Responses to SARS-CoV-2 Coronavirus in Humans with COVID-19 Disease and Unexposed Individuals CIGB-258, a peptide derived from human heat-shock protein 60, decreases hyperinflammation in COVID-19 patients Clinical features of patients infected with 2019 novel coronavirus in Wuhan Clinical features of patients infected with 2019 novel coronavirus in Wuhan Normal Values of T, B and NK Lymphocyte Subpopulations in Peripheral Blood of Healthy Cuban Adults Comprehensive mapping of immune perturbations associated with severe COVID-19 SARS-CoV-2-specific T cell immunity in cases of COVID-19 and SARS, and uninfected controls Highly functional virus-specific cellular immune response in asymptomatic SARS-CoV-2 infection SARS-CoV-2-specific T cells in asymptomatic Analysis of the long-term impact on cellular immunity in COVID-19-recovered individuals reveals a profound nkt cell impairment Genomic characterisation and epidemiology of 2019 novel coronavirus: implications for virus origins and receptor binding The who's who of T-cell differentiation: Human memory T-cell subsets Flattening the COVID-19 Curve With Natural Killer Cell Based Immunotherapies Deep immune profiling of COVID-19 patients reveals distinct immunotypes with therapeutic implications Human regulatory B cells in health and disease: therapeutic potential Sitio oficial de gobierno del Ministerio de Salud Pública en Cuba Protocol of followed in Cuba for COVID-19. Version 1.4. Havana. Ministry of Public Health of Cuba SARS-CoV-2-Specific T Cells Exhibit Phenotypic Features of Helper Function, Lack of Terminal Differentiation, and High Proliferation Potential Detection of SARS-CoV-2-Specific Humoral and Cellular Immunity in COVID-19 Convalescent Individuals Reappearance of effector T cells is associated with recovery from COVID-19 Efficiency of Neyman Allocation Procedure over other Allocation Procedures in Stratified Random Sampling Broad and strong memory CD4 + and CD8 + T cells induced by SARS-CoV-2 in UK Systems-Level Immunomonitoring from Acute to Recovery Phase of Severe COVID-19 T Cell Subpopulations in Healthy Elderly and Lung Cancer Patients: Insights from Cuban Studies Central memory and effector memory T cell subsets: function, generation, and maintenance Altered thymic activity in early life: how does it affect the immune system in young adults? Current Opinion in Immunology, Host pathogens/Immune senescence 23 Severe Pandemic H1N1 2009 Infection Is Associated with Transient NK and T Deficiency and Aberrant CD8 Responses Cellular immune correlates of protection against symptomatic pandemic influenza Robust T Cell Response Toward Spike, Membrane, and Nucleocapsid SARS-CoV-2 Proteins Is Not Associated with Recovery in Critical COVID-19 Patients Longitudinal Analysis of COVID-19 Patients Shows Age-Associated T Cell Changes Independent of Ongoing Ill-Health. Front Recovery from severe H7N9 disease is associated with diverse response mechanisms dominated by CD8+ T cells Exposure to SARS-CoV-2 generates T-cell memory in the absence of a detectable viral infection COVID-19 Severity Correlates with Weaker T-Cell Immunity, Hypercytokinemia, and Lung Epithelium Injury Phenotype and kinetics of SARS-CoV-2-specific T cells in COVID-19 patients with acute respiratory distress syndrome Immune cell profiling of COVID-19 patients in the recovery stage by single-cell sequencing Long-lasting severe immune dysfunction in Ebola virus disease survivors Clinical management of COVID-19: interim guidance Considerations for quarantine of contacts of COVID-19 cases: interim guidance World Medical Association Declaration of Helsinki: ethical principles for medical research involving human subjects Persistent memory CD4+ and CD8+ T-cell responses in recovered severe acute respiratory syndrome (SARS) patients to SARS coronavirus M antigen sub-populations analyzed were the increase of NK+(%), increase of NK+ (AV), increase of CD8+(%), increase of CD8+(AV), decrease of CD19+(%) and decrease of CD19+(AV) Model 2: Influence of the adjusted severity of age on the cellular sub-populations. Model 3: Influence of the duration of the evolution of the acute disease adjusting for age and severity of the disease on the main alterations of the cellular sub-populations Legend: a : statistical significance using proportion comparison between both sexes in the total number of COVID-19 convalescent individuals. b1 : statistical significance between asymptomatic patients and those severely ill, b2 : statistical significance between moderately ill and severely ill convalescent patients, in both cases p <0.0001 identified through the Mann-Whitney test. c : statistical significance between the duration of the disease in persons with moderate and severe forms of COVID-19 through the Mann-Whitney test. d : statistical significance between the duration of convalescence in individuals with asymptomatic and severe forms of COVID-19 identified through the Mann-Whitney test. * Duration of the disease: was defined as the time lapse between the diagnosis made by PCR-RT of SARS-CoV-2 infection and the first negative PCR as part of the criteria for clinical discharge. ** Convalescent study time: was defined from epidemiological discharge until day of blood sampling. IQR: interquartile range. For all tests, statistical significance was considered as p<0.05 