To explore the association of obesity comprehensively defined by Body Mass Index (BMI) and Body Fat percentage (BF%) with osteopenia.
MethodsIn this cross-sectional study, data of adult men and postmenopausal women aged ≥ 50 years old were obtained from the National Health and Nutrition Examination Surveys (NHANES) database. Weighted logistic regression analysis was conducted to investigate the association of BF% obesity with osteopenia in participants who had different gender and BMI obesity conditions. The association of obesity comprehensively evaluated by BMI and BF% with osteopenia was also explored in the total population and in gender subgroups.
ResultsAmong 1720 eligible subjects, 1054 had osteopenia. Multivariate analysis suggested that in males, BMI obesity combined with BF% obesity was associated with higher osteopenia odds compared to BMI obesity only (Odds Ratio [OR = 4.01], 95% Confidence Interval [95% CI 1.43‒11.27]). Compared to participants with both BMI and BF% obesity, those with BMI obesity have lower osteopenia odds (OR = 0.46, 95% CI 0.28‒0.76), whereas those with BF% obesity have higher odds of osteopenia (OR = 2.03, 95% CI 1.35‒3.05, p = 0.002). In females, compared to BMI obesity combined with BF% obesity, BF% obesity (OR = 3.37, 95% CI 1.47‒7.73) or non-obesity (OR = 2.11, 95% CI 1.18‒3.75) was respectively associated with higher osteopenia odds. In males, BMI obesity was linked to lower osteopenia odds compared to both BMI and BF% obesity (OR = 0.25, 95% CI 0.10‒0.62).
ConclusionsThe comprehensively assessed obesity by BMI and BF% may be more meaningful in the evaluation of potential osteopenia risk, as well as further prevention and intervention of osteoporosis.
Osteoporosis is a most common chronic disease in the middle-aged and elderly population, characterized by low bone mass, bone strength impairment, and increased fragility fracture risk, leading to heavy family and social burdens.1 Osteopenia is widely recognized to be the preclinical stage of osteoporosis, with a global prevalence of approximately 40.4%.2,3 Hence, early screening and intervening in osteopenia are significant for reducing the burdens of osteoporosis and related diseases.
The association between obesity and bone health is complicated, and current evidence is inconsistent. Although most previous studies suggested that obesity might be a protective factor related to osteoporosis and fractures,4,5 there is also a demonstration that obesity could not prevent osteoporosis and may even have harmful effects.6 Body Mass Index (BMI) is always a classical index used to measure obesity. In fact, one possible explanation for these inconsistent findings is the limitation that the BMI does not distinguish between fat mass and muscle mass of the human body.7 Besides, the Body Fat percentage (BF%) is a commonly used indicator that directly reflects the body fat mass. In recent years, researchers have proposed that BF% may be superior to BMI on predicting the risk of multiple diseases.8,9 In a prospective follow-up study, Lin et al.10 considered that the comprehensive evaluation of obesity through both BMI and BF% is better than using a single index in evaluating mortality risk in Chronic Kidney Disease (CKD) patients. Nevertheless, no study has explored the association of obesity that was measured comprehensively via the BMI and BF% with osteopenia.
Herein, the present study based on data extracted from the National Health and Nutrition Examination Survey (NHANES) database aims to explore the association of obesity defined by BMI and BF% comprehensively with osteopenia, and to provide information for early screening and intervention in the osteopenia population.
Materials and methodsStudy subjectsData in this cross-sectional study were extracted from the NHANES database in 2005‒2006, 2013‒2014, and 2017‒2018. The NHANES survey is conducted by both the Centers for Disease Control and Prevention (CDC) and the National Center for Health Statistics (NCHS), and assesses the United States noninstitutionalized population’s nutritional and health status. This database integrated a comprehensive, multistage stratified probability sample from selected counties, blocks, households, and individuals. The NCHS well-trained professionals conducted interviews at individuals’ homes, and performed extensive physical examinations at the Mobile Exam Centers (MECs). For detailed information on: https://www.cdc.gov/nchs/nhanes/index.htm.
A total of 6414 individuals were initially included. The inclusion criteria were: Adults aged ≥ 50-years; Males or postmenopausal females (postmenopausal status was defined as self-reported cessation of menstruation for ≥ 12-months); Availability of complete data on BMI, Body Fat Percentage (BF%), Bone Mineral Density (BMD), and covariates (sleep duration, marital status, waist circumference, total 25-hydroxyvitamin D, total energy intake, cotinine, and blood pressure).
The exclusion criteria were: BMI < 18.5 kg/m2 (to exclude underweight individuals, as low body weight may independently affect bone health through mechanisms unrelated to obesity); Missing data on osteopenia diagnosis or BF%; Missing values in any of the study covariates listed above.
Finally, 1720 of them were eligible for further analysis. This public database has obtained ethical approval from the Institutional Review Board (IRB) of the NCHS. Since information is de-identified, and participants have provided informed consent, the ethical approval has been waived by the IRB of Affiliated Hospital of Nanjing University of Chinese Medicine. This study follows the STROBE Statement.
Assessment of obesity through BMI and BF%Calculation of BMI was according to the formula: BMI = weight (kg) ÷ height2 (m2), and obesity was recognized as BMI of ≥ 30 kg/m2 referring to the World Health Organization (WHO) standard.11
In NHANES, whole-body dual-energy X-Ray Absorptiometry (DXA) scans were performed at the MEC through the advanced Hologic QDR 4500 A fan beam X-Ray bone densitometer manufactured by Hologic Inc. Before the DXA scans, persons underwent tests of nuclear medicine or radiographic contrast material during last 72 hours or 3 days, or reported a weight exceeding 300 pounds or a height over 6′·5′ were excluded.12 The DXA scan data were collected and analyzed by the Hologic Discovery software, version 12.1 after a rigorous quality control process. It can provide precise measurements of body composition in both total and regional, such as BF%. In this study, obesity was defined as BF% ≥ sex-specific median, where the cutoff values for males and females were respectively 28.2% and 41.3%. In addition, subjects were categorized into four groups based on the comprehensive assessment of obesity using BMI and BF%, including BMI obesity + BF% obesity, BMI obesity + BF% non-obesity, BMI non-obesity + BF% obesity, and BMI non-obesity + BF% non-obesity.
Diagnosis of osteopeniaBone Mineral Density (BMD) was measured using DXA scans at four skeletal sites: femoral neck, total femur, trochanter, and intertrochanter. In accordance with the World Health Organization (WHO) diagnostic criteria for osteopenia,14 the T-score for each site was calculated as follows:
The reference group comprised healthy adults aged 20–29 years from the NHANES database in 2005–2006, 2013–2014, and 2017–2018 cycles.13 Osteopenia was defined as a T-score between −2.5 and −1.0 at any of the measured skeletal sites (femoral neck, total femur, trochanter, or intertrochanter). This approach aligns with the WHO recommendation that osteopenia can be diagnosed based on T-scores from either the femoral neck or lumbar spine.14 Although lumbar spine BMD data were available in NHANES, this study prioritized hip-related measurements (femoral neck, total femur, etc.) to minimize potential confounding effects of degenerative changes in the spine common in older adults.
Variables selectionVariables could be potential confounding factors were also extracted, including age, race, marital status, gender, educational level, smoking, drinking, Poverty Income Ratio (PIR), physical activity, CKD, Cardiovascular Disease (CVD), anemia, anti-osteoporosis therapy, sleep duration, fracture history, waist circumference, height, weight, Systolic Blood Pressure (SBP), Diastolic Blood Pressure (DBP), fasting glucose, total 25-hydroxyvitamin D, cotinine, total energy intake, Calcium (Ca) intake and caffeine intake.
The information on physical activity was collected by the NHANES questionnaire and converted using the following formula: weekly energy expenditure (MET·min/week) = recommended Metabolic Equivalent (MET) × weekly exercise time of corresponding activity (min), and was categorized with the cut-off value 750 MET·min/week. CKD is recognized as the Urinary Albumin to Creatinine Ratio (UACR) of ≥ 30 mg/g and/or estimated Glomerular Filtration Rate (eGFR) < 60 mL/min/1.73 m2. The eGFR is calculated as eGFR = 175 × standardized Serum creatinine (Scr) −1.154 × age −0.203 × 1.212 (if black) × 0.742 (if female), in which the units for GFR and Scr are respectively mL/min/1.73 m2 of body surface area and mg/dL. Diagnosis of CVDs (including angina, heart attack, heart failure, coronary heart disease, congestive heart failure, stroke, or cardiovascular drug use) was self-reported. Anemia was diagnosed as serum hemoglobin of < 12 g/dL (males) or < 11 g/dL (females) or having taken treatment for anemia in the past 3-months. The NHANES collected dietary intake information through two 24-hour dietary recall surveys. The first one was conducted in person at the MEC, and the second one was through telephone or mail later in 3‒10 days.15 The total energy intake was calculated via dietary intake plus dietary supplement.
Statistical analysisData of continuous variables were presented using mean ± standard error (mean ± SE). Student t-test was utilized for the comparison of characteristics of participants between the osteopenia group and non-osteopenia group. Data of categorical variables were presented by frequency and constituent ratio [n (%)]. The Chi-Square test (χ2) was performed for comparison. In accordance with the NHANES recommendation, the “Full Sample 2 Year MEC Exam Weight (WTMEC2YR)” should be used for combination analysis of three 2-cycle data.
The covariates associated with osteopenia were screened through weighted univariate logistic regression analysis. Weighted multivariate logistic regression was utilized to investigate the association between BF% obesity and osteopenia in different gender and BMI obesity groups, which was evaluated by Odds Ratios (ORs) with 95% Confidence Intervals (95% CIs). Multivariate models adjusted for age, race, anemia, anti-osteoporosis therapy, fracture history and waist circumference. Restricted Cubic Spline (RCS) curves were drawn to fit the nonlinear association between BF% and osteopenia. The associations of obesity comprehensively evaluated by BMI and BF% with osteopenia were also explored in the total study population and in gender subgroups. Additionally, the Receiver Operating Character (ROC) curve with Area Under the Curve (AUC) was drawn to explore the predictive value of the comprehensive index on osteopenia. Statistical analyses were conducted with SAS 9.4 software (SAS Institute, Cary, NC, USA).
ResultsCharacteristics of study subjectsThe study process is shown in Fig. 1. To be specific, 6414 adult men and postmenopausal women aged ≥ 50 years were initially included. Then participants with BMI of < 18.5 kg/m2 (n = 394), missing information on osteopenia (n = 1036), BF% (n = 3110), marital status (n = 2), sleep duration (n = 2), waist circumference (n = 10), total 25-Hydroxyvitamin-D (n = 43), total energy intake (n = 64), cotinine (n = 5) and blood pressure (n = 28) were excluded 1720 of them were finally eligible.
Among the participants, 1054 had osteopenia. Comparison of characteristics between the non-osteopenia group and osteopenia group was shown in Table 1. More than half of the participants were aged 50‒59 years old (n = 1242). Male adults accounted for 61.10% (n = 1040). The numbers of individuals with anemia (34 vs. 16), anti-osteoporosis therapy (49 vs. 11), and fracture history (33.87 vs. 23.92) in the osteopenia group were significantly higher than those in the non-osteopenia group. Also, compared to those with non-osteopenia, the median height and weight was significantly lower in osteopenia patients, while the mean waist circumference and cotinine was significantly higher (all p < 0.05). 51.56% of participants had BMI obesity in the non-osteopenia group, whereas only 27.85% of participants had BMI obesity in the osteopenia group. Differently, the proportion of BF% obesity in the non-osteopenia group was 56.94%, and that in the osteopenia group was 49.23%.
Characteristics of participants between non-osteopenia group and osteopenia group.
t, t-test, χ2, Chi-Square test.
PIR, Poverty Income Ratio; CKD, Chronic Kidney Disease; CVD, Cardiovascular Disease; SE, Standard Error; SBP, Systolic Blood Pressure; DBP, Diastolic Blood Pressure; Ca, Calcium; BMI, Body Mass Index; BF%, Body Fat percentage.
Before exploring the association of BF% obesity with osteopenia in different gender and BMI obesity groups, covariates associated with osteopenia were screened (Supplementary Table 1). The authors observed that age, race, anemia, anti-osteoporosis therapy, fracture history, and waist circumference are significantly linked to the odds of osteopenia (all p < 0.05). After adjusting for selected covariates, the authors only observed that in males, having both BMI obesity and BF% obesity was linked to higher odds of osteopenia (OR = 4.01, 95% CI 1.43‒11.27, p = 0.011), compared to having BMI obesity only (Table 2). According to the RCS curve of the relationship between BF% and osteopenia odds in men with BMI obesity showed that when BF% > 28.2, the OR values of osteopenia > 1 (Fig. 2).
Association of BF% obesity with osteopenia in individuals with different gender and BMI obesity conditions.
| Groups | BF% obesity (outcome/total) | Unadjusted model | Adjusted modela | ||
|---|---|---|---|---|---|
| OR (95% CI) | p | OR (95% CI) | p | ||
| Male & BMI non-obesity | No (313/446) | Ref | Ref | ||
| Yes (158/235) | 0.77 (0.45‒1.32) | 0.332 | 1.60 (0.94‒2.74) | 0.082 | |
| Male & BMI obesity | No (23/74) | Ref | Ref | ||
| Yes (146/285) | 2.62 (1.08‒6.38) | 0.034 | 4.01 (1.43‒11.27) | 0.011 | |
| Female & BMI non-obesity | No (206/289) | Ref | Ref | ||
| Yes (102/127) | 1.42 (0.73‒2.74) | 0.290 | 1.64 (0.71‒3.77) | 0.229 | |
| Female & BMI obesity | No (18/50) | Ref | Ref | ||
| Yes (88/214) | 1.51 (0.48‒4.74) | 0.467 | 1.00 (0.24‒4.11) | 0.998 | |
BF%, Body Fat percentage; BMI, Body Mass Index; OR, Odds Ratio; CI, Confidence Interval; Ref, Reference.
The authors further assessed the association of obesity defined by BMI and BF% with osteopenia (Table 3). After adjusting for covariates, comparing to participants with BMI obesity combined with BF% obesity, those with BMI obesity only have lower odds of osteopenia (OR = 0.46, 95% CI 0.28‒0.76, p = 0.004), whereas those with BF% obesity only have higher odds (OR = 2.03, 95% CI 1.35‒3.05, p = 0.002). Similarly, Fig. 3 clearly shows the association of comprehensively evaluated obesity with osteopenia.
Association of BMI and BF% compositive defined obesity with osteopenia.
| Variables | Unadjusted model | Adjusted modela | ||
|---|---|---|---|---|
| OR (95% CI) | p | OR (95% CI) | p | |
| BMI obesity & BF% obesity | Ref | Ref | ||
| BMI obesity & BF% non-obesity | 0.47 (0.25‒0.86) | 0.017 | 0.46 (0.28‒0.76) | 0.004 |
| BMI non-obesity & BF% obesity | 2.30 (1.54‒3.45) | <0.001 | 2.03 (1.35‒3.05) | 0.002 |
| BMI non-obesity & BF% non-obesity | 2.47 (1.88‒3.26) | <0.001 | 1.40 (1.00‒1.97) | 0.053 |
BMI, Body Mass Index; BF%, Body Fat percentage; OR, Odds Ratio; CI, Confidence Interval; Ref, Reference.
Besides, the association of BMI and BF% comprehensively evaluated obesity with osteopenia was explored in gender subgroups (Table 4). In females, compared to BMI obesity combined with BF% obesity, BF% obesity only (OR = 3.37, 95% CI 1.47‒7.73, p = 0.007) or non-obesity (OR = 2.11, 95% CI 1.18‒3.75, p = 0.014) was respectively linked to higher odds of osteopenia. Among males, BMI obesity only was linked to lower odds of osteopenia compared to both BMI and BF% obesity (OR = 0.25, 95% CI 0.10‒0.62, p = 0.005).
Association of BMI and BF% compositive defined obesity with osteopenia in gender subgroups.
BMI, Body Mass Index; BF%, Body Fat percentage; OR, Odds Ratio; CI, Confidence Interval; Ref, Reference.
Adjusted for age, race, anemia, anti-osteoporosis therapy, fracture history and waist circumference.
Moreover, Fig. 4 shows the ROC curve of the predictive value of the comprehensive index on osteopenia. It clearly suggested that the AUC was 0.71 (0.68‒0.74). Also, the true positive rate of obesity evaluated comprehensively by BMI and BF% was 0.8284, indicating that it had a relatively high sensitivity in predicting osteopenia (Supplementary Table 2).
DiscussionIn the present study, the association of obesity comprehensively assessed by BMI and BF% with osteopenia in adults aged ≥50 years old was investigated based on the NHANES database. The study results suggested that men with both BMI and BF% obesity seemed to have higher odds of osteopenia compared to those with BMI obesity only. Also, comparing with individuals with both BMI and BF% obesity, those who with only BMI obesity have lower odds of osteopenia while those with only BF% obesity have higher odds, after adjusting covariates associated with osteopenia. Moreover, in females, BF% obesity only or non-obesity was respectively associated with higher odds of osteopenia compared to both BMI and BF% obesity.
There is a complex relationship between obesity and bone. Although the broadly accepted notion that obesity has a positive effect on bone health, evidence has challenged this view.6,16 In recent years, the coexistence of osteoporosis and sarcopenia has been considered as a syndrome.17 However, due to the BMI is a traditional measurement index for obesity that cannot reflect muscle mass, it may have lower priority than BF% in the assessment of the association between obesity and multiple diseases.8,9 Besides, sarcopenic obesity is a newly identified pathological entity that is characterized by an increased body fat mass with an associated sarcopenia, which was also associated with bone health.18 Based on these theoretical foundations, the present study explored the association of obesity measured comprehensively through BMI and BF% with osteopenia, with the aim of providing some new idea for prevention and prevention in the pre-disease stage of osteoporosis. Lin et al.10 performed a prospective cohort study that included CKD patients from a medical center in Taiwan and found that BMI-defined obesity was related to significantly lower mortality risk, while the relationship was reversed when obesity was calculated using BF%. The authors similarly observed that persons with BMI obesity have lower odds of osteopenia while those with BF% obesity have higher odds, compared with those who had obesity assessed by both BMI and BF%. Lin et al.10 considered this diagnostic discordance may explain the obesity paradox to some extent. Although no study has reported the same topic before, subjects in this study were from the NHANES database, which includes large United States representative populations, so these findings were relatively reliable and may supplement some information on the association between comprehensive measurement of obesity and osteoporosis. However, causal associations of BMI, BF%, and comprehensively measured obesity with bone health need further clarified.
Adipose tissue and bone tissue interact with each other. Adipose tissue plays a role in the musculoskeletal system through a series of substances, including definable hormones such as leptin19,20 and adiponectin.21 Besides, obesity as a chronic dysfunction condition with characteristics of a low-grade, systemic inflammatory status. The adiponectin concentration is also inversely proportional to numerous inflammatory cytokines, such as Interleukin (IL)-6, C-Reactive Protein (CRP), and Tumor Necrosis Factor (TNF)-α. Adiponectin can stimulate osteoblastic proliferation, with an increased activity of alkaline phosphatase, as well as the formation of type I collagen and osteocalcin, which are markers of osteoblasts’ differentiation and maturation.7 The concentration of adiponectin is usually presented at a low level in obesity.22 Herein, it can be presumable that in obesity, chronic inflammatory status can express a high concentration of these inflammation markers, which potently inhibit the expression of adiponectin that may further suppress differentiation and maturation of osteoblasts. The effect of Oxidative Stress (OS) in the development of osteoporosis is also considerable.23 OS in obesity, especially sarcopenic obesity, is a crucial factor in muscle function alteration and metabolic dysfunction development.24 In addition, a recent study has revealed the interaction between ferroptosis and TNF-α and its impact on osteogenesis and angiogenesis, which might be related to the pathogenesis and regenerative therapy of obesity-related osteoporosis.25 Herein, obesity comprehensively measured by both BMI and BF% may partly take full account of the effects of changes in both fat mass and muscle mass. However, due to the observational nature of this study, the mechanism that obesity comprehensively evaluated by both BMI and BF% may more comprehensively reflect the correlation between obesity and osteopenia is still unclear.
The authors further investigated these associations in different gender subgroups. Notably, compared to both BMI and BF% obesity, non-obesity was associated with higher odds of osteopenia in postmenopausal females, which may be related to the phenomenon “obesity paradox”.26 However, the potential protective effect of increased BMI and obesity on bone health based on age and gender has not reached a consensus.7,27 In fact, osteoporosis and fracture risk increases along with age, independently from body weight, and, in menopause, the abrupt decrease of estrogen is the primary cause of osteoporosis. Therefore, the authors only included postmenopausal females to avoid this biasing factor. Besides, although obesity measured by BMI and BF% had a positive association with osteopenia risk, it seemed that BF% evaluated obesity played a more negative role. On the contrary, regarding to males, BMI obesity rather than BF% obesity was a potential protective factor for osteopenia, which was consistent with the conclusions of previous studies.28,29
The present findings are further supported by epidemiological evidence highlighting the interplay between obesity and metabolic dysfunction in postmenopausal women. A cross-sectional study of 5027 Brazilian postmenopausal women demonstrated that obesity (BMI > 30 kg/m²) was prevalent in 30% of the cohort and significantly associated with adverse metabolic profiles, including hypertension, elevated triglycerides, fasting glucose, and reduced HDL-C levels.30 These metabolic disturbances, such as insulin resistance and dyslipidemia, may exacerbate bone loss by promoting chronic inflammation and altering adipokine secretion (e.g., reduced adiponectin),31 potentially explaining the stronger association between BF% obesity and osteopenia observed in the female subgroup. This aligns with the hypothesis that obesity, particularly when characterized by high body fat, reflects a systemic pro-inflammatory state that accelerates bone resorption.
Moreover, the authors have compared multiple influencing factors related to osteopenia between osteopenia patients and non-osteopenia persons, including total energy intake, Ca intake, caffeine intake and total 25-hydroxyvitamin-D concentration.32,33 Unfortunately, no significant difference has been observed. The predictive value of obesity evaluated by both BMI and BF% on osteopenia was assessed, and the AUC of the ROC curve showed a high sensitivity of this comprehensive index. Our results indicated that in addition to the usual recommendation for osteoporosis screening, focusing on comprehensively evaluated obesity may have the potential to be a significant link in osteoporosis prevention. In clinical, it could be an easier and low-cost tool for early screening of populations with high risk of osteopenia. Also, following medical advice to adjust dietary nutrition and take appropriate physical activity to keep both healthy BMI and BF% may be beneficial to reduce the risk of osteopenia and subsequence osteoporosis.
As far as the authors know, this is the first to explore the association of obesity comprehensively evaluated by BMI and BF% with osteopenia in old adults. In this study, BMI and BF% were combined to more accurately distinguish different types of persons with obesity and more comprehensively reflect the correlation between obesity and osteopenia. These results suggested that BF% may be a more valuable indicator of obesity for osteopenia, and provide relevant evidence support for subsequent preventive interventions in osteoporosis. However, there are still some limitations in the present research. Due to the cross-sectional study design, the authors are unable to clarify the causal association of obesity with osteopenia. In addition, the classification cutoff values of BF% were based on the median values of males and females respectively in the study participants, which limited the extrapolation of these results in other populations. Therefore, prospective cohort studies are needed to further validate and supplement the conclusions in the future.
ConclusionObesity evaluated by both BMI and BF% may more comprehensively reflect the association between obesity and osteopenia. In clinical practice, not only BMI but also BF% worth attention in the intervention and prevention of osteoporosis in old adults.
DeclarationsEthics approval and consent to participate: Not applicable, because NHANES belongs to public databases, the patients involved in the database have obtained ethical approval, users can download relevant data for free for research and publish relevant articles, and the present study is based on open-source data, and the Affiliated Hospital of Nanjing University of Chinese Medicine do not require research using publicly available data to be submitted for review to their ethics committee, so there are no ethical issues and other conflicts of interest.
Consent for publicationNot applicable, because this paper did not reveal any personal information of patients.
Availability of data and materialsThe datasets used and/or analyzed during the current study were publicly available from the NHANES database.
Authors’ contributionsXin Liu, Yan Lou and Guangquan Sun conceived and designed the study.
Xin Liu, Yan Lou, Zhiyong Chang, Changyuan Gu and Bin Du collected the data.
Xin Liu, Yan Lou, Zhiyong Chang, Changyuan Gu and Bin Du analyzed and interpreted the data.
Xin Liu and Yan Lou writing the manuscript.
Guangquan Sun provided critical revisions that are important for the intellectual content.
All authors approve the final version of the manuscript.
FundingThis study was supported by the National Natural Science Foundation of China (n° 82,074,471) and the General Foundation of Jiangsu Province Traditional Chinese Medicine Bureau (n° MS2022020).
The authors declare no conflicts of interest.
This study was supported by the National Natural Science Foundation of China (n° 82074471) and the General Foundation of Jiangsu Province Traditional Chinese Medicine Bureau (n° MS2022020).










