Literature DB >> 33841339

Poor Metabolic Health Increases COVID-19-Related Mortality in the UK Biobank Sample.

Filip Morys1, Alain Dagher1.   

Abstract

Previous studies link obesity and components of metabolic health, such as hypertension or inflammation, to increased hospitalizations and mortality of patients with COVID-19. Here, in two overlapping samples of over 1,000 individuals from the UK Biobank we investigate whether metabolic health as measured by waist circumference, dyslipidemia, hypertension, type 2 diabetes, and systemic inflammation is related to increased COVID-19 infection and mortality rate. Using logistic regression and controlling for confounding variables such as socioeconomic status, age, sex or ethnicity, we find that individuals with worse metabolic health (measured on average eleven years prior to 2020) have an increased risk for COVID-19-related death (adjusted odds ratio: 1.75). We also find that specific factors contributing to increased mortality are increased serum glucose levels, systolic blood pressure and waist circumference.
Copyright © 2021 Morys and Dagher.

Entities:  

Keywords:  COVID-19; diabetes; dyslipidemia; hypertension; obesity

Mesh:

Substances:

Year:  2021        PMID: 33841339      PMCID: PMC8027509          DOI: 10.3389/fendo.2021.652765

Source DB:  PubMed          Journal:  Front Endocrinol (Lausanne)        ISSN: 1664-2392            Impact factor:   5.555


Introduction

Since the beginning of the COVID-19 pandemic, mounting evidence supports an association between obesity and poor outcomes (1–12). The association holds for both obesity and obesity-associated metabolic health – hypertension, diabetes, dyslipidemia, and systemic inflammation (13). Similarly, obesity and excess adipose tissue have also been associated with higher risk of SARS-CoV-2 infection (10, 14, 15). Previous studies, however, have tended to use small sample sizes, focus predominantly on the effects of body mass index (BMI) as a measure of obesity, or not account for confounding factors, such as ethnicity or socioeconomic status (14, 16–18). Since ethnicity and socioeconomic status are themselves associated with obesity and metabolic health (19, 20), they could confound interpretation of analyses in patients with COVID-19. Here, we aim to present a comprehensive evaluation of obesity-associated metabolic risk factors that might be related to poor health outcomes in SARS-CoV-2 infected patients while controlling for confounding variables and limiting potential collider bias, which has previously resulted in incorrect epidemiological conclusions (21–23). We investigate whether metabolic health is related to higher chance for SARS-CoV-2 infection, but also COVID-19-related death.

Materials and Methods

Participants

In this study, we used the UK Biobank dataset – a large scale study with extensive phenotyping carried out in the United Kingdom (24). This study was performed under UK Biobank application ID 35605. SARS-CoV-2 real-time PCR test results in the UK Biobank dataset are derived from the Public Health England microbiology database Second Generation Surveillance System that is dynamically linked to the UK Biobank database (25). Here we only included individuals who were recorded as tested for SARS-CoV-2. We distinguished between two samples for two aims of our project: Sample 1 – a larger sample (n=12,659) of all individuals who were tested for SARS-CoV-2 between 16th March 2020 and 24th August 2020, to investigate the risk of COVID-19 infection and how it is related to metabolic health; and Sample 2, a subset of Sample 1 consisting of individuals who tested positive for SARS-CoV-2 (n=1,152). We also obtained data on mortality from COVID-19 for all the individuals included in our study population. Sample characteristics can be found in . The discrepancies in number of patients who tested positive in both samples is due to outlier exclusions (see section: ‘Measures used in the study’). Participants were classified as having metabolic syndrome based on the criteria in (26). All participants signed written informed consent prior to participating in the UK Biobank study, which was approved by the North-West Multi-centre Research Ethics Committee (11/NW/0382). All UK Biobank actions are overseen by the UK Biobank Ethics Advisory Committee.
Table 1

Participants characteristics.

MeasureSample 1 (n=12,659)Sample 2 (n=1,148)
SARS-CoV-2 positive patientsmean (SD)SARS-CoV-2 negative patientsmean (SD)p-valueSARS-CoV-2 positive survivorsmean (SD)SARS-CoV-2 positive deceasedmean (SD)p-value
n1,13811,521970178
COVID-19 mortality rate1.42%15.51%
Age in years68 (8)70 (8)<0.00167 (9)74 (6)<0.001
Sex distribution47.36% women52.32% women0.00249.07% women36.52%0.003
Waist circumference (cm)*93.38 (13.89)92.07 (13.69)0.00392.56 (13.63)98.49 (14.31)<0.001
Body mass index (kg/m2)*28.40 (5.15)27.87 (4.90)<0.00128.22 (4.96)29.57 (6.06)0.006
Obesity (BMI>30kg/m2)*38.14%28.01%0.13629.38%35.96%0.096
Metabolic syndrome*1.58%1.56%1.0001.24%4.49%0.006
Serum HDL (mmol/l)*1.36 (0.34)1.43 (0.38)<0.0011.37 (0.45)1.34 (0.31)0.366
Serum TG (mmol/l)*1.78 (1.04)1.77 (1.02)0.5451.78 (1.05)1.88 (1.01)0.207
HbA1c (mmol/mol)*35.83 (4.51)35.76 (4.33)0.63735.73 (4.63)36.94 (5.19)0.004
Serum glucose (mmol/l)*4.98 (0.64)5.00 (0.63)0.3184.95 (0.63)5.20 (0.74)<0.001
Serum C-reactive protein (mg/l)*3.00 (5.23)2.84 (4.45)0.3082.82 (4.73)4.01 (7.25)0.036
Systolic blood pressure (mmHg)*137.20 (18.59)138.55 (18.55)0.019135.68 (18.15)145.74 (18.43)<0.001
Diastolic blood pressure (mmHg)*82.46 (10.25)82.24 (9.94)0.48482.19 (10.30)83.99 (10.06)0.030
Diabetes*9.23%7.61%0.0558.66%15.17%0.010
Hypertension*41.65%42.32%0.68538.14%62.36%<0.001
Liver disease2.46%2.63%0.8072.68%1.69%0.605
Pulmonary disease25.3%25.22%0.97924.85%29.21%0.256
Cardiovascular disease14.94%16.22%0.27913.40%25.28%<0.001
Kidney disease7.64%6.86%0.3497.11%10.67%0.137
Townsend deprivation index*0.14 (1.04)-0.07 (0.96)<0.001-0.06 (0.98)0.04 (1.03)0.215
Cigarette smokers*12.57%11.89%0.19712.58%12.36%0.108

SD, standard deviation; BMI, body mass index; HDL, high density lipoprotein cholesterol; TG, triglycerides; HbA1C, hemoglobin A1c. Uncorrected p-values reflect significance of within samples differences calculated with t-tests in case of numerical variables and chi-squared tests in case of categorical variables.

Note that measures with an asterisk were obtained on average 11 years prior to 2020, while the other measures are from the time of testing.

Bold values indicate significant between-group differences.

Participants characteristics. SD, standard deviation; BMI, body mass index; HDL, high density lipoprotein cholesterol; TG, triglycerides; HbA1C, hemoglobin A1c. Uncorrected p-values reflect significance of within samples differences calculated with t-tests in case of numerical variables and chi-squared tests in case of categorical variables. Note that measures with an asterisk were obtained on average 11 years prior to 2020, while the other measures are from the time of testing. Bold values indicate significant between-group differences.

Measures Used in the Study

To investigate how metabolic health is related to SARS-CoV-2 infection and mortality rate, we used the following measures: waist circumference, serum triglyceride (TG), serum high density lipoprotein cholesterol (HDL), glycated hemoglobin (HbA1c), serum glucose (corrected for fasting times prior to blood drawing), serum C-reactive protein (27), previous type 2 diabetes diagnosis, resting systolic and diastolic blood pressure (mean of two measurements each), and hypertension diagnosis. In our analyses, we also controlled for age, sex, socioeconomic status [Townsend deprivation index (28)], highest achieved educational qualifications, mean family income, smoking status, and ethnic background, the laboratory where COVID-19 testing was performed, and the origin of the sample used for COVID-19 test (e.g. nose, throat etc.). Based on the previous literature, we also controlled for past liver disease (ICD10 identifiers K7), cardiovascular disease (ICD10 identifiers I2 and I3), pulmonary disease (ICD10 identifiers J), and kidney disease (ICD10 identifiers N1) in all our analyses, as they are all related to the severity of COVID-19 (29–32). For Sample 2, we used mortality data provided by the UK Biobank – COVID-19-related death was described using the ICD10 identifier U07.1. All variables related to metabolic health and all confounding variables were collected on average 11 years prior to COVID-19 tests. Prior to the analyses, all numeric variables were z-scored, all serum level values were log-transformed, and we excluded outliers from the samples (2.2 interquartile range below 1st or above 3rd quartile). No other exclusions were applied. If participants were tested for SARS-CoV-2 more than once, they were considered positive if at least one test result was positive.

Statistical Analyses

The same analyses were performed for each of the two population samples. Data were analyzed using R (v. 3.6.0). First, using confirmatory factor analysis in lavaan (v. 0.6-7), we estimated the fit of a latent variable ‘metabolic health’, which consisted of waist circumference, serum C-reactive protein levels, a latent variable ‘dyslipidemia’ [serum HDL and TG levels; with loadings constrained to be equal between the two variables (33)], and two other latent variables, diabetes and hypertension. The latent variable ‘diabetes’ consisted of serum glucose and HbA1c levels, and type 2 diabetes diagnosis, while the latent variable ‘hypertension’ consisted of blood pressure measurements (systolic and diastolic blood pressure), and hypertension diagnosis (). The model was estimated using robust maximum likelihood estimation and model’s fit was evaluated using common indices: comparative fit index (CFI), root mean square error of approximation (RMSEA), and standardized root mean square residual (SRMR). Acceptable fit was defined as CFI>0.9, RMSEA<0.1, and SRMR<0.08.
Figure 1

Confirmatory factor analysis model used to derive the latent variable ‘metabolic health’. Squares indicate measured variables; circles indicate latent variables. HbA1c, glycated hemoglobin A1c; HDL, high density lipoprotein cholesterol; TG, triglycerides; CRP, C-reactive protein.

Confirmatory factor analysis model used to derive the latent variable ‘metabolic health’. Squares indicate measured variables; circles indicate latent variables. HbA1c, glycated hemoglobin A1c; HDL, high density lipoprotein cholesterol; TG, triglycerides; CRP, C-reactive protein. Next, we extracted components of the latent variable ‘metabolic health’ for each participant and entered them in a logistic regression. The outcome variable in logistic regression for Sample 1 was COVID-19 test result, while for Sample 2 the outcome variable was COVID-19-related death. In the analyses we used a set of confounding variables to calculate adjusted odds ratio: age, sex socioeconomic status, smoking status, ethnic background, test laboratory, and sample origin (6, 34–36). Finally, for Sample 2 we explored how individual factors contributed to the COVID-19-related mortality by using a logistic regression with individual components of metabolic health, instead of the latent variable ‘metabolic health’. Overall, analyses using a latent variables approach allowed us to maximize information shared between different measured variables within the same metabolic domains, e.g. serum glucose levels, HbA1C levels and type 2 diabetes diagnosis for the latent variable ‘diabetes’. Furthermore, adding continuous measures in this study, such as glucose levels or blood pressure, instead of only using type 2 diabetes or hypertension diagnosis, enabled us to use more information that is available in the dataset and investigate in depth how continuous metabolic health measures are related to SARS-CoV-2 infection and COVID-19 mortality. A script for the analysis of the data as well as the output of statistical software can be found at https://github.com/FilipMorys/COVID_MetS.

Results

Metabolic Health and the Risk of COVID-19 Infections

In Sample 1, the confirmatory factor analysis provided an acceptable model fit (CFI=0.932, RMSEA=0.065, SRMR=0.040). In the logistic regression the relation between metabolic health and the chance of having a positive SARS-CoV-2 test did not reach our pre-set statistical significance threshold (p=0.059; odds ratio 1.10; 95% confidence intervals (CI): 1.00-1.22).

Metabolic Health Is Related to an Increased COVID-19-Related Mortality

Among the individuals who tested positive, 178 people (16%) died from COVID-19, allowing us to investigate how metabolic health influences COVID-19-related mortality. Here, the confirmatory factor analysis provided an acceptable model fit (CFI=0.919, RMSEA=0.068, SRMR=0.049). Logistic regression showed that obesity-associated metabolic impairment was related to an increased mortality rate among COVID-19 positive individuals (p<0.001) – adjusted odds ratio: 1.75 (95% CI: 1.30-2.37), pointing to a 75% increase for each unit increase on the metabolic health latent variable (). Each unit increase on the metabolic health variable means increased serum glucose levels by 0.21 mmol/l, HbA1c levels by 3.21 mmol/mol, C-reactive protein levels by 2.18 mg/l, triglyceride levels by 0.67 mmol/l, systolic blood pressure by 7.17 mmHg, diastolic blood pressure by 4.27 mmHg, waist circumference by 18.72 cm, and decreased HDL cholesterol levels by 0.24 mmol/l.
Figure 2

Odds ratio for COVID-19-related death. Circles represent adjusted odds ratios, error bars represent 95% confidence intervals.

Odds ratio for COVID-19-related death. Circles represent adjusted odds ratios, error bars represent 95% confidence intervals. Exploratory analysis with individual components of metabolic health revealed that higher systolic blood pressure, higher serum glucose, and higher waist circumference increased the chance of COVID-19-related mortality – the adjusted odds ratio were 1.32 (95% CI: 1.02-1.70), 1.27 (95% CI: 1.04-1.55), and 1.34 (95% CI: 1.05-1.72), respectively (systolic blood pressure: p=0.035, glucose levels: p=0.020, waist circumference: p=0.019; ).

Discussion

We investigated whether obesity-associated poor metabolic health, here defined by increased waist circumference, increased TG levels, decreased HDL levels, systemic inflammation, increased glucose and HbA1C levels, increased systolic and diastolic blood pressure, and the presence of type 2 diabetes and hypertension, is a risk factor for COVID-19 infection and mortality. We were able to show that a one unit increase of the latent variable ‘metabolic health’ results in a 75% higher risk of death because of COVID-19. In contrast, the relation between metabolic health and the likelihood of test positivity was weaker, with an odds ratio of only 1.1 and a confidence interval that included the null effect. In sum, poor metabolic health contributes little to no risk of test positivity but substantially increases the odds of an adverse outcome. A strength of our analysis lies in the fact that the risk factor variables were measured prior to COVID-19 diagnosis and therefore were not influenced by the disease itself or possible treatment. Another strength is that we controlled for known shared risk factors between obesity and COVID-19 outcomes, such as socioeconomic status, ethnicity, and sex. We also used history of pulmonary, cardiovascular, kidney, and liver disease as covariates in our analyses, as these have been linked to COVID-19 severity and mortality (29–32). The results of our study are in line with previous reports linking obesity and individual components of metabolic health with poor COVID-19 outcomes and death (1, 5, 7, 10, 12, 18, 37). In particular, the fact that waist circumference was related to COVID-19 mortality is consistent with a study showing that visceral adipose tissue is significantly related to the risk of admission to an intensive care unit in COVID-19 patients (38). Here, the authors argue that this might be in part due to increased local and systemic inflammation. In their review, Bansal et al. discuss some of the physiological mechanisms that might mediate the relationship between obesity, related comorbidities and worse outcomes of COVID-19 (39). These include an enhanced expression of the angiotensin converting enzyme 2 (ACE2), diabetes-related microvascular dysfunction, increased expression of pro-inflammatory cytokines, or activation of the renin-angiotensin-aldosterone system related to hypertension (39–41). Recent reports especially highlight the role of interleukin-6 (IL-6) in the pathogenesis of COVID-19, but also the role of IL-6 inhibitors in COVID-19 therapy (42, 43). This is important given that adiposity is associated with chronic systemic inflammation, which generally delays immune response to pathogens and can also lead to worse outcomes in COVID-19 patients (37, 39, 40, 44). In general, our results suggest that a number of these mechanisms might also contribute to increased mortality from COVID-19. Previous studies investigating COVID-19 in the UK Biobank used positive test results obtained between March 16th and April 26th as a proxy of severe COVID-19 (18, 45). The rationale for this is that, during this time, only patients admitted to hospitals and with COVID-19-like symptoms were tested for SARS-CoV-2. Here, we decided to not use this approach for several reasons. First, it is possible that individuals with COVID-19-like symptoms for which they were admitted to a hospital but who did not have COVID-19 were only infected after being admitted. In those cases, positive tests would not reflect COVID-19 severity. Second, it is not possible to determine the exact reason for which inpatients were tested for SARS-CoV-2; positive test results might therefore not only reflect severity of COVID-19 disease, but also testing in anticipation of isolating patients admitted for other reasons. For example, patients having to undergo unrelated medical procedures might have received precautionary testing. We therefore recommend that studies that used this approach as a proxy for COVID-19 severity be interpreted with caution. For the interpretation of our and similar results from the UK Biobank, it is important to note that the UK Biobank is not a sample representative of the entire UK population and therefore the findings might not be generalizable (46). Infection and mortality rates calculated here should not be used as an indicator of true prevalence and mortality rate in the general population. Furthermore, observational studies such as ours are prone to collider bias, which has already been identified in UK Biobank COVID-19 investigations (23). Collider bias occurs when the sample population is conditioned on a variable that correlates with the variables of interest. For example, during the time period of the current study, it is thought that health workers were more likely to be tested for SARS-CoV-2, which may have contributed to incorrect conclusion that cigarette smoking is protective, as health care workers have a lower incidence of smoking (23). Our study population includes individuals who were seen in a health care setting, who may therefore have a higher incidence of obesity and poorer metabolic health than the general population. However, it is difficult to see how this would account for the effect of obesity on death from COVID-19 among people who tested positive for the virus. In addition, COVID-19 testing in the UK during the time period of this study was generally restricted to individuals with symptoms such as fever, cough, or loss of smell or taste – asymptomatic COVID-19 individuals were less likely to be tested. This might further increase the extent of collider bias in this and similar studies investigating the predictors of COVID-19 severity or mortality. Current strategies to account for such bias or measure the extent thereof rely on models and strong assumptions that might be incorrect. We therefore suggest that our findings be interpreted with caution. Nonetheless, our model did account for potential collider variables or confounds such as socioeconomic status, ethnicity and age. Finally, in this study we only used limited data pertaining to COVID-19 diagnosis or mortality. Information on the seriousness of illness, treatment administered, or actual cause of death were not available. In sum, we used the UK Biobank dataset to confirm that, in individuals who tested positive for COVID-19 in the early stages of the pandemic, metabolic health, and especially visceral adiposity, hypertension, and serum glucose levels, were associated with an increased risk of death.

Data Availability Statement

The data analyzed in this study is subject to the following licenses/restrictions: The dataset can be obtained from the UK Biobank. Requests to access these datasets should be directed to https://www.ukbiobank.ac.uk.

Ethics Statement

The studies involving human participants were reviewed and approved by North-West Multi-centre Research Ethics Committee. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

Both authors contributed equally to all stages of the research. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by a Foundation Scheme award to AD from the Canadian Institutes of Health Research. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
  40 in total

1.  The UK Biobank sample handling and storage protocol for the collection, processing and archiving of human blood and urine.

Authors:  Paul Elliott; Tim C Peakman
Journal:  Int J Epidemiol       Date:  2008-04       Impact factor: 7.196

2.  UK biobank: an open access resource for identifying the causes of a wide range of complex diseases of middle and old age.

Authors:  Cathie Sudlow; John Gallacher; Naomi Allen; Valerie Beral; Paul Burton; John Danesh; Paul Downey; Paul Elliott; Jane Green; Martin Landray; Bette Liu; Paul Matthews; Giok Ong; Jill Pell; Alan Silman; Alan Young; Tim Sprosen; Tim Peakman; Rory Collins
Journal:  PLoS Med       Date:  2015-03-31       Impact factor: 11.069

3.  Comparison of Sociodemographic and Health-Related Characteristics of UK Biobank Participants With Those of the General Population.

Authors:  Anna Fry; Thomas J Littlejohns; Cathie Sudlow; Nicola Doherty; Ligia Adamska; Tim Sprosen; Rory Collins; Naomi E Allen
Journal:  Am J Epidemiol       Date:  2017-11-01       Impact factor: 4.897

4.  Chronic kidney disease is associated with severe coronavirus disease 2019 (COVID-19) infection.

Authors:  Brandon Michael Henry; Giuseppe Lippi
Journal:  Int Urol Nephrol       Date:  2020-03-28       Impact factor: 2.370

5.  Obesity and risk of COVID-19: analysis of UK biobank.

Authors:  Thomas Yates; Cameron Razieh; Francesco Zaccardi; Melanie J Davies; Kamlesh Khunti
Journal:  Prim Care Diabetes       Date:  2020-05-28       Impact factor: 2.459

6.  Early administration of interleukin-6 inhibitors for patients with severe COVID-19 disease is associated with decreased intubation, reduced mortality, and increased discharge.

Authors:  Pranay Sinha; Anahita Mostaghim; Catherine G Bielick; Angela McLaughlin; Davidson H Hamer; Lee M Wetzler; Nahid Bhadelia; Maura A Fagan; Benjamin P Linas; Sabrina A Assoumou; Michael H Ieong; Nina H Lin; Ellen R Cooper; Karrine D Brade; Laura F White; Tamar F Barlam; Manish Sagar
Journal:  Int J Infect Dis       Date:  2020-07-25       Impact factor: 3.623

7.  Association Between Midlife Obesity and Its Metabolic Consequences, Cerebrovascular Disease, and Cognitive Decline.

Authors:  Filip Morys; Mahsa Dadar; Alain Dagher
Journal:  J Clin Endocrinol Metab       Date:  2021-09-27       Impact factor: 5.958

8.  Metabolic syndrome and COVID-19: An update on the associated comorbidities and proposed therapies.

Authors:  Fernanda Farias Costa; Wilian Reis Rosário; Ana Cláudia Ribeiro Farias; Ramon Guimarães de Souza; Roberta Sabrine Duarte Gondim; Wermerson Assunção Barroso
Journal:  Diabetes Metab Syndr       Date:  2020-06-11

9.  Metabolic Syndrome and COVID-19 Mortality Among Adult Black Patients in New Orleans.

Authors:  John Xie; Yuanhao Zu; Ala Alkhatib; Thaidan T Pham; Frances Gill; Albert Jang; Stella Radosta; Gerard Chaaya; Leann Myers; Jerry S Zifodya; Christine M Bojanowski; Nassir F Marrouche; Franck Mauvais-Jarvis; Joshua L Denson
Journal:  Diabetes Care       Date:  2020-08-25       Impact factor: 19.112

Review 10.  COVID-19 Susceptibility in chronic obstructive pulmonary disease.

Authors:  Jordi Olloquequi
Journal:  Eur J Clin Invest       Date:  2020-09-02       Impact factor: 5.722

View more
  5 in total

Review 1.  The Endocrine Adipose Organ: A System Playing a Central Role in COVID-19.

Authors:  Francesca Cinti; Saverio Cinti
Journal:  Cells       Date:  2022-07-04       Impact factor: 7.666

2.  Relationship between hyperglycemia, waist circumference, and the course of COVID-19: Mortality risk assessment.

Authors:  Mykola Khalangot; Nadiia Sheichenko; Vitaly Gurianov; Viola Vlasenko; Yulia Kurinna; Oksana Samson; Mykola Tronko
Journal:  Exp Biol Med (Maywood)       Date:  2021-10-20

Review 3.  Lessons From COVID-19: Physical Exercise Can Improve and Optimize Health Status.

Authors:  Dario Cerasola; Christiano Argano; Salvatore Corrao
Journal:  Front Med (Lausanne)       Date:  2022-05-13

4.  Visceral fat inflammation and fat embolism are associated with lung's lipidic hyaline membranes in subjects with COVID-19.

Authors:  Georgia Colleluori; Laura Graciotti; Mauro Pesaresi; Angelica Di Vincenzo; Jessica Perugini; Eleonora Di Mercurio; Sara Caucci; Patrizia Bagnarelli; Cristina M Zingaretti; Enzo Nisoli; Stefano Menzo; Adriano Tagliabracci; Annie Ladoux; Christian Dani; Antonio Giordano; Saverio Cinti
Journal:  Int J Obes (Lond)       Date:  2022-01-26       Impact factor: 5.551

5.  Comparing body mass index and obesity-related comorbidities as predictors in hospitalized COVID-19 patients.

Authors:  Michael W Tsoulis; Victor L Garcia; Wei Hou; Chrisa Arcan; Joshua D Miller
Journal:  Clin Obes       Date:  2022-02-22
  5 in total

北京卡尤迪生物科技股份有限公司 © 2022-2023.