Abstract Metabolic reprogramming is a distinctive characteristic of SARS-CoV-2 infection, which refers to metabolic changes in hosts triggered by viruses for their survival and spread. It is current urgent to understand the metabolic health status of COVID-19 survivors and its association with long-term health consequences of infection, especially for the predominant non-severe patients. Herein, we show systemic metabolic signatures of survivors of non-severe COVID-19 from Wuhan, China at six months after discharge using metabolomics approaches. The serum amino acids, organic acids, purine, fatty acids and lipid metabolism were still abnormal in the survivors, but the kynurenine pathway and the level of itaconic acid have returned to normal. These metabolic abnormalities are associated with liver injury, mental health, energy production, and inflammatory responses. Our findings identify and highlight the metabolic abnormalities in survivors of non-severe COVID-19, which provide information on biomarkers and therapeutic targets of infection and cues for post-hospital care and intervention strategies centered on metabolism reprogramming. Keywords: COVID-19 survivors, Metabolic abnormalities, Metabolomics, Health consequences __________________________________________________________________ COVID-19 survivors; Metabolic abnormalities; Metabolomics; Health consequences. 1. Introduction Coronavirus disease 2019 (COVID-19), a global pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has resulted in 523,786,368 confirmed cases with 6,279,667 deaths, as of 24 May 2022 ([45]https://covid19.who.int/). About 80% of confirmed COVID-19 patients presented with non-severe symptoms, and the proportion is increasing as vaccination and treatment protocols mature ([46]Huang et al., 2020; [47]Costanzo et al., 2022a). Patients were generally discharged within one month, but they may not have “truly recovered” due to the long-term sequelae and debilitation. The epidemiological cohort studies have reported that COVID-19 survivors mainly suffered from fatigue or muscle weakness, sleep difficulties, anxiety or depression and dyspnoea at six months after infection, and these symptoms retained even one year ([48]Huang et al., 2021a, [49]2021b). Besides, confirmed persistent symptoms in survivors of SARS and Middle East respiratory syndrome (MERS) also call attention on “long COVID-19”. Metabolic reprogramming is initiated by virus for their replication and survival and also modulated by the immune response after SARS-CoV-2 infection ([50]Shen and Wang, 2021; [51]Xiao et al., 2021). Several metabolic pathways, mainly including amino acids, energy, lipids and immune-related metabolism, were reprogrammed in hospitalized COVID-19 patients ([52]Su et al., 2020; [53]Shen et al., 2020; [54]Costanzo et al., 2022b). These metabolic changes had not completely returned to normal even the discharge criteria of the official care protocol for COVID-19 were met ([55]Wu et al., 2020). Metabolites are more sensitive than clinical indicators to reflect the condition of the body. Considering the fact that metabolic changes usually progress to a chronic state, we propose that the effects of SARS-CoV-2 infection on metabolism may be long-lasting, especially in the absence of intervention and treatment after discharge. Metabolic health plays a crucial role in regulating the long-term health consequences and convalescent duration for COVID-19 patients, given the close relationship between metabolism and pathophysiological states. However, metabolic health status of COVID-19 survivors after discharge has not been completely elucidated. In the present study, we characterized the serum metabolic signatures of survivors of non-severe COVID-19 at six months after discharge from the first pandemic wave in Wuhan, China, and further explored the potential adverse outcomes associated with these metabolic abnormalities. The findings are expected to provide more information on understanding of coronavirus pathology and biomarkers and the health cares and intervention targets of COVID-19 survivors after discharge. 2. Results 2.1. Survivors exhibit differences in demographics, clinical characteristics and metabolic profile from controls Twenty-one past-infected volunteers (survivors) with non-severe COVID-19 symptoms and 22 control volunteers without SARS-CoV-2 infection history were recruited at Renmin Hospital of Wuhan University, Wuhan, China. Their clinical data are presented in Table S1. Survivors of non-severe COVID-19 had a median age of 50.2 years (young: 31.0 ± 4.4; elder: 56.2 ± 6.7) and 28.5% were male, and control volunteers had a median age of 54.5 years (young: 29.6 ± 2.9; elder: 61.9 ± 8.3) in which 36.4% were male. Compared to controls, red blood cell counts, hemoglobin, albumin, total protein and aspartate aminotransferase (AST) were significantly decreased in COVID-19 survivors, but high-density lipoprotein cholesterol (HDL-c) decreased only in elder survivors. Changes in the first two indices suggested the hematopoiesis in the bone marrow and the oxygen transport in blood was inhibited to a certain extent, and that in others might be associated with hepatic dysfunction. We integrated non-targeted and targeted metabolomic approaches to characterize the difference of serum metabolic profile of survivors of non-severe COVID-19 and controls. Total of 116 polar metabolites and 455 lipids were detected in serum samples (Table S2). Orthogonal partial least squares discriminant analysis (OPLS-DA) found a dramatic difference between the serum metabolic profiles of survivors and controls ([56]Figure 1a). Our results suggest that the metabolic-health effects of SARS-CoV-2 infection still existed six months after discharge. According to the screening criteria for significantly differential metabolites (in the Methods section), 83 of 165 metabolites with |log2(fold change)| ≥ 0.25, and p < 0.05 in volcano plot were significantly varied by meeting variable importance in the projection (VIP) > 1.0 and area under curve (AUC) > 0.75, mainly including amino acids, organic acids, fatty acids and lipids ([57]Figure 1b and Table S3). As shown in [58]Figure 1c and 1d, the numbers of differential metabolites in elder group were higher than that of young group (159 vs 30, Table S3), suggesting the metabolic abnormality in elder survivors might be more profound although the sample size of young survivors is small. Change trends of the differential metabolites from each class in both age groups were basically consistent. The top differential metabolites (VIP ≥2.0) in all, young and elder samples were integrated and analyzed, and the change trends showed that most of them were reduced in serum from survivors of non-severe COVID-19 ([59]Figure 1e). The significantly differential metabolites in all survivors were further performed pathway analysis, and the result showed that the primary abnormal metabolic pathways of survivors of non-severe COVID-19 were amino acid metabolism, TCA cycle, lipid metabolism and fatty acid metabolism ([60]Figure 1f and Table S4). Moreover, similar metabolic abnormalities were also found in blood of in-hospital COVID-19 patients ([61]Mahmud and Garrett, 2020; [62]Aggarwal et al., 2021), demonstrating that metabolic reprogramming induced by SARS-CoV-2 infection had not returned to normal at six months after discharge. The persistent metabolomic abnormalities in COVID-19 survivors may require a longer recovery period ([63]Li et al., 2022). Figure 1. [64]Figure 1 [65]Open in a new tab Abnormal serum metabolites associated with the history of SARS-CoV-2 infection. (a) OPLS-DA score plot (R^2Y = 0.96, p = 0.003; Q^2 = 0.60, p < 0.001 for permutation test) obtained from modelling metabolic profile in the survivors of non-severe COVID-19 (red) and controls (blue). Volcano plots for differential metabolites in all (b), young (c), and elder (d) samples. (e) Heat maps of top differential metabolites in all, young, and elder samples. The color of each column represents the logarithm of fold change of metabolites with base 2, and p values were calculated using a student's t-test, ∗p < 0.05, ∗∗p < 0.01, and ∗∗∗p < 0.001. (f) The main KEGG pathways disturbed in survivors of non-severe COVID-19. 2.2. Abnormal amino acid and purine metabolism indicate hepatic injury and mental disorders Amino acids are mainly synthesized in the liver, and play important roles in the immune system and maintaining redox status. Hepatic injury is frequently present in hospitalized COVID-19 patients ([66]Zhang et al., 2020; [67]Caterino et al., 2021). The main manifestations of hepatic injury, elevated alanine aminotransferase and aspartate aminotransferase levels, were not observed in discharged COVID-19 patients in our cohort, but the decreased albumin, total protein and AST might be associated with hepatic dysfunction. As shown in [68]Figure 2, serum ornithine, aspartate, asparagine, alanine, glycine, serine, proline, histidine and threonine were significantly down-regulated (p < 0.05) in COVID-19 survivors compared with controls, while arginine, glutamate, valine, taurine, phenylalanine, lysine, glutamine and methionine were significantly up-regulated. Moreover, the change trends were similar in young and elder survivors (Figure S1). Further pathway enrichment analysis results show these significantly altered amino acids in COVID-19 survivors are involved in the urea cycle, ammonia recycling, glycine and serine metabolism, arginine and proline metabolism, aspartate metabolism and glutamate metabolism (Figure S2). Urea cycle converts toxic ammonia into urea in the liver, and its disorder also is the common feature of hepatic injury, such as acute and chronic liver dysfunction and liver carcinogenesis ([69]Bigot et al., 2017). In COVID-19 patients, most of urea cycle-related metabolites were decreased ([70]Shen et al., 2020; [71]Song et al., 2020). We found that the urea cycle was still disordered with the decrease of ornithine and increase of arginine in survivors of non-severe COVID-19. Therefore, hepatic injury caused by SARS-CoV-2 infection has not been recovered six months after discharge, and might mainly result from the urea cycle blocked in the arginine node. Figure 2. [72]Figure 2 [73]Open in a new tab Disorders in targeted metabolites in survivors of non-severe COVID-19. Heatmap representation of Z-score transformed concentrations of the amino acids and amino acid derivatives, organic acids and neurotransmitters. Each column represents a single sample, and the dark red in color scale is for the highest value, white is for the midpoint value and dark blue is for the lowest value. Anxiety, depression, disturbed sleep and post-traumatic stress disorder (PTSD) are common adverse mental health symptoms to the COVID-19 pandemic ([74]Rajkumar, 2020; [75]Taquet et al., 2021). Glutamate is the main excitatory neurotransmitter in the central nervous system, and its disorder affects mental health ([76]Nasir et al., 2020). The serum level of glutamate was significantly increased in survivors of non-severe COVID-19, and can contribute to the abovementioned mental symptoms through neuron damage ([77]Zeredo et al., 2019). PTSD was the most common condition among COVID-19 survivors, and the high occurrence rate (42.9%) of denial or avoidance of infection history was found in our cohort (Table S1), and the higher glutamate levels was found in elder survivors with PTSD (Figure S3). Elevated glutamate in survivors also suggested that the energy metabolism of their central nervous system was weakened in survivors, which could explain their remaining symptoms of fatigue and sleep disturbance ([78]Tarazona et al., 2021; [79]Drews et al., 2020). Furthermore, purine metabolism was inhibited in survivors of non-severe COVID-19 with the reduction of inosine, deoxyinosine and hypoxanthine. As end-product of purine metabolism, low serum levels of uric acid in elder survivors (Table S3) might increase the risks of multiple sclerosis, neurodegenerative diseases and optic neuritis ([80]Kutzing and Firestein, 2008). These results suggested that the raising awareness of long-term mental health care in COVID-19 survivors is urgent. 2.3. Altered organic acid and carbohydrate metabolites indicate inhibition of cellular energy production SARS-CoV-2 infection can decrease serum levels of organic acids, and further severely inhibit energy metabolism ([81]Song et al., 2020). We found that organic acids involved in the TCA cycle, succinic acid and 2-ketoglutaric acid remained significantly down-regulated in serum samples from survivors of non-severe COVID-19 ([82]Figure 2 and Figure S1), suggesting that the damaged energy metabolism has not returned to normal. Notably, the levels of citric acid, malic acid, fumaric acid and aconitic acid showed no significant difference in survivors compared to controls, demonstrating that the TCA cycle inhibited by SARS-CoV-2 infection was recovering. Glycolysis is adopted to facilitate virus replication and cytokine production after SARS-CoV-2 infection, and the elevated glucose was observed in COVID-19 patients ([83]Wu et al., 2020). Glycolysis remained inhibited in young survivors of non-severe COVID-19, with significant increase in glucose-6-phosphate (G6P), and which might be a major risk factor for new-onset diabetes and kidney injury in COVID-19 survivors ([84]Rubino et al., 2020; [85]Sundaram et al., 2021). The TCA cycle and glycolysis are metabolic pathways for the generation of ATP, but the former is more highly efficient pathway. The significant down-regulations of lactic acid and AST were also evidence for the disrupted energy metabolism (Figure S1 and [86]Figure 2e). Therefore, ATP production of COVID-19 survivors was mainly attenuated by weakened TCA cycle, further resulted in fatigue or muscle weakness of the COVID-19 survivors (41.7%, in Table S1). 2.4. Disturbed glycerophospholipid and glycerolipid metabolism are more significant in elder survivors Lipids are envelope components of SARS-CoV-2, and play important roles in host-pathogen interactions and pathogenesis ([87]Abu-Farha et al., 2020). Dysregulated lipid metabolism in COVID-19 patients was ubiquitous and closely related to the levels of disease severity ([88]Su et al., 2020). Even when patients has met clinical discharge criteria, their lipid metabolism still remains abnormal ([89]Wu et al., 2020). In our study, lipid metabolic dysregulations sustained in survivors of non-severe COVID-19 six months after discharge (Figure S4), and it may be more obvious in elder cases (Table S3). As shown in [90]Figure 3, the total levels of lysophosphatidylcholines (LysoPC), phosphatidic acids (PA), phosphatidylcholines (PC), phosphatidylethanolamines (PE), phosphatidylserines (PS), and ceramides (Cer) in elder survivors were significantly down-regulated, except for triglycerides (TG). In young cases, no lipid classes presented significant change (p > 0.05). These initial observations suggest that glycerophospholipid and sphingolipid metabolism could be inhibited in elder survivors, but glycerolipid metabolism could be activated. The serum lipid metabolic disruptions in survivors of non-severe COVID-19 differed from that of SARS survivors twelve years after infection, with increase of phosphatidylinositols (PI), lysophosphatidylinositols (LysoPI) and LysoPC ([91]Wu et al., 2017), and a possible explanation for this is the high-dose steroid treatment for SARS patients. Figure 3. [92]Figure 3 [93]Open in a new tab Metabolic network of significantly changed lipid classes in all, young and elder survivors of non-severe COVID-19. Bars in all graphs represent mean abundance of controls (light blue column) and survivors (pink column) in different subgroups, error bars represent standard deviation. The normalized total peak abundance of individual lipid class is shown on the Y-axis. Red, green, black and gray colors are significantly up-regulated, down-regulated, unchanged and undetected lipids or metabolites, and the arrows indicate the direction of metabolism. The dotted-line box includes all free fatty acids (FFA), which convert to Acyl-CoA for participate the synthesis of LPA and PA. p values were calculated using a student's t-test, ∗p < 0.05, ∗∗p < 0.01, and ∗∗∗p < 0.001. MUFA: monounsaturated fatty acids, PUFA: polyunsaturated fatty acids, MHFA: mono-hydroxyl fatty acids; DHFA: di-hydroxyl fatty acids, G6P: glucose-6-phosphate, Glyeraldehyde-3-P: glyeraldehyde-3-phosphate, Glycerol-3-P: glycerol-3-phosphate, PC: phosphatidylcholines, LysoPC: lysophosphatidylcholines, PA: phosphatidic acids, LysoPA: lyso phosphatidic acids, PE: phosphatidylethanolamines, LysoPE: lisophosphatidylethanolamines, PS: phosphatidylserines, MG: monoglycerides, DG: diglycerides, and TG: triglyceride. As important components of bio-membranes, glycerophospholipids are essential for intracellular metabolism and transmembrane signaling. Total abundance of PC, PE, and PS were significantly down-regulated in elder survivors, indicating that their transformation was markedly inhibited in the whole glycerophospholipid metabolism, and cell membrane structure could be destructed. As major precursors of PS and PC, the significant decrease of serine and choline in the serum of elder survivors also proved this inhibition (Figure S5). Cytosolic phospholipase A2α (PLA2α), which hydrolyzes glycerophospholipids to free fatty acids (FFAs) and lysophospholipids, is activated after SARS-CoV-2 infection ([94]Shen et al., 2020). The down-regulation of total lysoPC and FFA in elder survivors demonstrated the abnormal activation of PLA2α had returned to normal. Most of PC and PE species with at least one unsaturated fatty acyl chain were significantly reduced in elder COVID-19 survivors six months after discharge, and several PA, PI and PS species were also reduced (Table S3). As important energy storage units, total TG levels were significantly up-regulated in elder survivors, demonstrating the biosynthesis of TG might be accelerated along with the reduction of MG, to maintain the balance of the energy metabolism in the body. It was closely associated with fatigue or muscle weakness syndrome of the COVID-19 survivors in our cohort. The accumulation of TG was also found in the plasma of COVID-19 survivors after discharge ([95]Xu et al., 2021; [96]Li et al., 2022). Besides, elevated TG and reduced HDL-c might increase the risk of coronary heart diseases in elder survivors ([97]Krauss and Siri, 2004). 2.5. Alterations in inflammation mediators are associated with age As bioactive signaling mediators, FFAs and their oxygenated products are involved in numerous functions, including inflammation, coagulation and apoptosis ([98]Dennis and Norris, 2015). Fatty acid metabolism was significantly inhibited in survivors of non-severe COVID-19, and varied with the ages of survivors ([99]Figure 3). The total levels of mono-hydroxyl fatty acids (MHFAs) were significantly down-regulated in both young and elder survivors, while FFAs reduced only in elder cases, and di-hydroxyl fatty acids (DHFAs) reduced only in young cases. MHFAs are active anti-inflammatory and pro-resolution mediators ([100]Serhan et al., 2008), two and five of MHFA species were respectively reduced in young and elder discharged patients, whose reduction demonstrates a potential deficit in the resolve of inflammatory in survivors of non-severe COVID-19. [101]Li et al. (2022) also found that two MHFA species were downregulated in plasma of severe and non-severe COVID-19 survivors 6 months after discharge from the hospital. In elder survivors, 38.8% of FFA species, including FA 16:0, FA 18:0, FA 18:1, FA 18:2, FA 20:1, FA 20:4 and FA 22:6, were significantly down-regulated (Table S3), demonstrating that the elder survivors might have a poor resistance to infection and/or disease. It was different from the elevation of serum FFAs in recovered SARS patients twelve years after infection ([102]Wu et al., 2017). Unsaturated fatty acids can inactivate enveloped viruses and inhibit their proliferation, but their changes were inconsistent in human blood after SARS-CoV-2 infection ([103]Su et al., 2020; [104]Shen et al., 2020). Further, DHFAs with eighteen carbon atoms were significantly reduced in young survivors (Table S3). Among them, DHFA from FA 18:2 can inhibit inflammatory and immune responses, regulate the vascular development and repair process, negatively regulate the pain, enhance increase exercise and skeletal muscle regulation ([105]Hildreth et al., 2020). Thus, a high fatty acid diet is suggested for survivors of non-severe COVID-19 for the elevation of the levels of FFAs and their oxygenated products both in blood and tissues to boost their immune system. 2.6. Recovered immune-related metabolites indicate the macrophages in normal state Immune dysregulation is common in COVID-19 patients. The kynurenine pathway participates in the regulation of the body's immune system, and is significantly activated in SARS-CoV-2 patients ([106]Thomas et al., 2020). Activation of this pathway can help to enhance the phagocytic ability of macrophages and turn on the switch for macrophage effector responses by autonomic synthesis of NAD^+ ([107]Minhas et al., 2019). We found the serum levels of metabolites in the kynurenine pathway, including tryptophan, kynurenine, kynurenic acid, xanthurenic acid and 5-hydroxy-tryptamine were not significantly altered from those of the controls, suggesting that kynurenine pathway in survivors of non-severe COVID-19 has been returned to normal. As an anti-inflammatory metabolite secreted by macrophages, itaconic acid is significantly up-regulated in the early stages of SARS-CoV-2 infection ([108]Song et al., 2020). However, itaconic acid in survivors of non-severe COVID-19 was also not significantly different from that in the control cases. In short, recovered kynurenine pathway metabolites and itaconic acid indicated macrophages activated by SARS-CoV-2 infection might have returned to normal. 3. Discussion SARS-CoV-2 infection triggers immune response and inflammatory storm accompanied by metabolic reprogramming. Up to now, most studies on metabolic changes in COVID-19 patients focused on acute infection and recovery periods for exploring the pathogenesis, prognosis, and potential biomarkers of infection ([109]Su et al., 2020; [110]Shen et al., 2020; [111]Wu et al., 2020; [112]Hasan et al., 2021). Reprogramming of amino acids, glucose, organic acids, lipids and fatty acids is common in hospitalized COVID-19 patients. These reprogrammed metabolisms not only offer energy and essential substrates for SARS-CoV-2 replication, but also regulate immune and inflammation response. Notably, the long-lasting effects of infection on metabolic health in post-convalescence may be very considerable, not just during acute infection and recovery periods ([113]Ayres, 2020). Although the metabolic activities were increased in COVID-19 patients after clinical treatment relative to acute infection period, the partial metabolic changes continued until discharge, maybe even months after discharge ([114]Wu et al., 2020; [115]Zheng et al., 2021; [116]Shi et al., 2021; [117]Li et al., 2022). This study aim to illustrate the systemic metabolic abnormalities of survivors of non-severe COVID-19 at six months after discharge by serum metabolomics analyses. We found that amino acids, organic acids, purine, fatty acids and lipid metabolism of the survivors still are abnormal compared to controls ([118]Figure 4), but the kynurenine pathway metabolites and itaconic acid had returned to normal levels. It was indicated that the effects of SARS-CoV-2 infection on health are profound. As reported, plasma glycerophospholipid, purine, and phenylalanine metabolism, and TCA cycle in the moderate/severe/critical patients without previous underlying diseases at three months after discharge were different from those of controls ([119]Zhang et al., 2021). Moreover, metabolic profiles of survivors with organ dysfunctions, such as abnormal pulmonary function, decreased glomerular filtration rate or cerebral white matter alterations also have been changed ([120]Xu et al., 2021; [121]Zhou et al., 2021; [122]Yang et al., 2021). Figure 4. [123]Figure 4 [124]Open in a new tab Abnormalities of body metabolism in survivors of non-severe COVID-19. Metabolic reprogramming is not only a prominent feature of SARS-CoV-2 infection, but also an important contributor to the progression and prognosis. Sequelae are very common among COVID-19 survivors, fatigue or muscle weakness and PTSD were dominant in our study cohort. The metabolic features of hepatic injury, mental disorders, fatigue or muscle weakness, and inflammation in survivors have been found. Abnormal energy metabolism for ATP production is a typical characteristic of fatigue or muscle weakness, and also include amino acids, nitrogen, oxidative stress and hormone metabolic imbalance ([125]Armstrong et al., 2014). Coronaviruses are neurotropic, neurovirulent, and neuroinvasive, and the point prevalence of PTSD was 32.2% in the post-illness stage of SARS and MERS ([126]Rogers et al., 2020). Glutamate is involved in the pathophysiology of PTSD via its actions on the initiation and maintenance of the hypothalamic-pituitary-adrenal response and learning and memory formation ([127]Averill et al., 2017). Energy metabolism and inflammatory pathways associated with stress responses, such as urea cycle, TCA cycle, purine metabolism and glycolysis, also occurred in PTSD mice ([128]Miller et al., 2020). Furthermore, lipid and fatty acid metabolism are closely associated with PTSD-related psychopathology ([129]Karabatsiakis et al., 2015). In conclusion, the systemic metabolic changes induced by infection are still present and might further affect the long-term health outcomes of COVID-19 survivors. The long-term physiological and mental health care to promote metabolic health are needed for COVID-19 survivors. To this end, targeting these metabolic pathways may be an important approach to prevent and treat complications and further improve health-related quality of life. 4. Limitations of the study The chemical signatures of COVID-19 survivors were identified using metabolomics in the present study, but there were a number of limitations. First, we cannot rule out the possibility that the survivors had taken medication, especially for elders with previous underlying diseases. Second, we recruited as many recovered patients as possible, but the sample size was small. Third, serum samples reflect the systemic metabolic abnormalities of survivors, being the sum of the metabolic contributes of many organs and tissues. Fourth, these findings are from patients infected by the original strain in China, and different virus variant, vaccination, racial, ethnic, and geographical cohorts were not considered in our study. Finally, confirmatory studies is needed for extending the current understanding of metabolic health in COVID-19 survivors. Declarations Author contribution statement Fang Li and Lei Fu: Performed the experiments; Analyzed and interpreted the data; Wrote the paper. Xiaoxiong Liu, Yueguang Lv, Zhiyu Chen and Wenbo Li: Performed the experiments. Xin-an Liu: Conceived and designed the experiments; Performed the experiments. Yong Liang, Qian Luo and Fan Pan: Conceived and designed the experiments. Zhiyi Yang and Ang Guo: Analyzed and interpreted the data. Funding statement Fang Li was supported by National Natural Science Foundation of China [22076197]. Prof. Qian Luo was supported by National Natural Science Foundation of China [82127801], Scientific Instrument Developing Project of the Chinese Academy of Sciences [YJKYYQ20200034], Shenzhen Science and Technology Innovation Commission [JCYJ20200109115405930, JCYJ20210324115811031]. Data availability statement Data will be made available on request. Declaration of interests statement The authors declare no conflict of interest. Additional information No additional information is available for this paper. Acknowledgements