Article Text

Download PDFPDF

Effect of different ranges of systolic blood pressure on left ventricular structure and diastolic function in a Chinese population: a cross-sectional population-based Shunyi study
  1. Liu Yongtai1,
  2. Lai Jinzhi1,
  3. Zhou Lixin2,
  4. Zhai Feifei2,
  5. Zhang Dingding3,
  6. Tian Zhuang1,
  7. Zhu Yanlin1,
  8. Chen Wei1,
  9. Bai Hua1,
  10. Wang Hui1,
  11. Zhu Yicheng2,
  12. Cui Liying2,
  13. Jin Zhengyu4,
  14. Shuyang Zhang1
  1. 1 Cardiology, Peking Union Medical College Hospital, Dongcheng-qu, Beijing, China
  2. 2 Neurology, Peking Union Medical College Hospital, Dongcheng-qu, Beijing, China
  3. 3 Central Research Laboratory, Peking Union Medical College Hospital, Dongcheng-qu, Beijing, China
  4. 4 Radiology, Peking Union Medical College Hospital, Dongcheng-qu, Beijing, China
  1. Correspondence to Professor Shuyang Zhang; shuyangzhang103{at}163.com

Abstract

Objectives To evaluate the effect of different ranges of systolic blood pressure (SBP) on left ventricular (LV) geometry and diastolic function in Chinese population.

Design Cross-sectional study.

Setting Peking Union Medical College Hospital in Beijing, China.

Participants All inhabitants aged 35 years or older, living in five villages of Shunyi were invited. Exclusion criteria included individuals who declined participation, presence of moderate to severe valvular heart disease, persistent atrial fibrillation and suboptimal echocardiograms.

Interventions The baseline data of 1051 participants were analysed. The relationship between SBP and LV geometric and diastolic function assessed by echocardiography was analysed after adjusting for conventional cardiac risk factors.

Results The adjusted value of SBP was independently associated with LV hypertrophy (LVH) and LV diastolic dysfunction (LVDDF) (all p<0.01). Setting individuals with SBP <120 mm Hg as the reference group (group 1), those with SBP between 120 mm Hg and 140 mm Hg (group 2) had higher risk odds of LVH and those with SBP ≥140 mm Hg (group 3) had higher risk odds of LVH and LVDDF (all p<0.01). With the increase of SBP, LV mass index (LVMI) and E/e′ stepwise increased and e′ stepwise decreased significantly from group 1 to 3 (all p<0.05). In the whole population, SBP was independently correlated with LVMI, LVEDD, Left Atrial Volume Index, e′, and E/e′ (all p<0.01).

Conclusions SBP was independently related to LVH and LVDDF, SBP between 120 and 140 mm Hg was independently related to worse LV remodelling and diastolic function, these findings indicated the potential benefit of intensive SBP control.

  • cardiac epidemiology
  • hypertension
  • echocardiography

This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.

Statistics from Altmetric.com

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

Strengths and limitations of this study

  • The strengths of the study include its large sample from the general Chinese population and its standardised acquisition of basic characteristics and echocardiography examinations.

  • A limitation is that this is a cross-sectional study and cannot prove a causal relationship between the geometric/functional changes and systolic blood pressure.

  • Another limitation is that this study was conducted in a Chinese population; studies that involve other ethnic populations are needed to replicate our results.

Introduction

It is well known that hypertension is a very important risk factor for cardiovascular disease.1–4 According to 2018 European Society of Cardiology guideline, hypertension was defined as office systolic blood pressure (SBP) values at least 140 mm Hg and/or diastolic BP (DBP) values at least 90 mm Hg.5 A complex series of changes in the geometry and function occur in the hearts of patients with hypertension.6 In these patients, the alterations in left ventricular (LV) geometry, delayed myocardial relaxation, increased passive stiffness of the sarcomere and extracellular matrix, and altered myocardial tone will cause LV remodelling and influence diastolic function.7 In the latest study, impairment of LV geometry and diastolic function was observed even in prehypertensive patients,8–10 who were defined as having a SBP ranging from 120 to 139 mm Hg and/or DBP from 80 to 89 mm Hg.

Elevations in SBP, particularly among older adults, are usually the result of stiffness in the large arteries and SBP appears to be a better predictor of events than DBP.11 12 Increasing emphasis has been placed on SBP alone as the most useful predictor of cardiovascular disease in hypertensive individuals and several guidelines have supported the use of SBP and strongly encourage the use of SBP goals in antihypertensive treatment.10

In the newest Chinese guideline,13 the ‘high normal’ SBP was defined as SBP between 120 and 140 mm Hg. Our hypothesis is that SBP between 120 and 140 mm Hg, regardless of the DBP level, increases the risk odds of LV remodelling and LV diastolic dysfunction (LVDDF) in Chinese adults.

The present study was a part of the ongoing community-based Shunyi study in China, designed to investigate the risk factors for cardiovascular and age-related diseases. In this study, we analysed the effect of different ranges of SBP on LV remodelling and diastolic function in the Chinese population.

Methods

Population

All inhabitants aged 35 years or older, living in five villages of Shunyi, a suburb district of Beijing, were invited, and the demographic data were recorded from June 2013 to April 2016. Exclusion criteria included individuals who declined participation, presence of moderate to severe valvular heart disease, persistent atrial fibrillation and suboptimal echocardiograms. Finally, 1586 individuals participated and standard baseline assessments were undertaken. All participants were invited for echocardiographic examination. Among those, 461 participants refused. Thus, the baseline echocardiography was performed in 1125 (70.9%) participants. The individuals with incomplete clinical data (59 individuals), moderate to severe valvular disease (6 individuals), persistent atrial fibrillation (4 individuals), and poor image quality (5 individuals) were excluded. The final analysis was performed based on 1051 subjects (figure 1).

Figure 1

The participants’ flow chart with respect to inclusion and exclusion. A total of 1586 inhabitants aged 35 years or older were invited. Exclusion criteria included individuals who declined participation, presence of moderate to severe valvular heart disease, persistent atrial fibrillation and suboptimal echocardiograms. Finally, 1051 individuals were included in this study.

Patient and public involvement

Neither patients nor public were involved in the development of the research question, study design, outcome measures, recruitment to and conduct of the study. The results of the study will be disseminated to participants through lectures presented by the investigators.

Assessment of the conventional risk factors affecting LV geometry and function

Conventional cardiovascular risk factors affecting LV geometry and function were selected and included in the analysis based on medical knowledge and previous reports included age, gender, heart rate (HR), body mass index (BMI), physical activity, smoking and alcohol drinking status, SBP, antihypertensive treatment, blood glucose level, serum low-density lipoprotein cholesterol (LDL-C) level, history of coronary heart disease (CHD) and renal function.

BMI was calculated by the weight in kilograms indexed to the height in metres squared. Smoking status was classified as current smoking (at least within the past 1 month) and non-current smoking. Alcohol drinking status was classified into current alcohol drinking (at least once per week within the past 1 month) and non-current alcohol drinking. We also collected the self-report degree of individual physical activity including five grades: very light, light, moderate, heavy and very heavy physical activity.

BP was measured by trained nurses and physicians on sitting position, supporting the limb used to measure BP, ensuring the BP cuff was at heart level, using a validated oscillometric device (Omron HEM-7211) three times after a 5 min rest and the mean value was used. Each BP measurement took place in at least 30 s interval, having the patient sit quietly for 5 min before a reading is taken. Hypertension was defined as self-reported hypertension, antihypertensive treatment, SBP ≥140 mm Hg or DBP ≥90 mm Hg. The device used in the biochemical test was Beckman AU 5800. Venous blood samples, routinely drawn after an overnight fast, were analysed for glycated haemoglobin A1c (HbA1c), plasma LDL-C and plasma glucose levels. Because of the large fluctuations in fasting blood glucose, HbA1c was chosen as an indicator of the blood glucose level. History of CHD was confirmed by two independent experienced cardiologists according to the patient’s medical history, symptoms, electrocardiography and other imaging data. If there was disagreement about the final diagnosis, a third cardiologist refereed. Renal function was expressed as estimated glomerular filtration rate (eGFR), which was calculated using the Chronic Kidney Disease Epidemiology Collaboration (CKD-EPI) equation, eGFR=141× min (Scr/κ, 1)α×max (Scr/κ, 1)–1.209×0.993Age×1.018 (if female), where Scr is serum creatinine, κ is 0.7 for females and 0.9 for males, α is –0.329 for females and –0.411 for males, min indicates the minimum of Scr/κ or 1, and max indicates the maximum of Scr/κ or 1.

Echocardiography

Echocardiography was performed using commercially available equipment (Vivid I; GE Vingmed Ultrasound, Horten, Norway). All studies were performed and reviewed by cardiologists with advanced training in echocardiography.

Left atrial anteroposterior dimension (LAD), LV end-diastolic dimension (LVEDD), interventricular septal wall thickness (IVST), LV posterior wall thickness (LVPWT), LV end-diastolic volume (LVEDV), LV end-systolic volume (LVESV), LV ejection fraction(), LV mass (LVM) and Left Atrial Volume Index (LAVI) were assessed according to the Guideline of the American Society of Echocardiography.14 LAD, LVEDD, IVST and LVPWT were measured in the parasternal long-axis view using M-mode echocardiography. Relative wall thickness (RWT) was calculated with the following formula: 2×LVPWT/LVEDD. LA volume was measured in a biplane calculation from both the apical four and two-chamber views, and LAVI was calculated by LA volume indexed to body surface area (BSA). LVM was calculated by linear dimension technique, and LVM index (LVMI) was the result of LVM indexed to BSA. LVEDV, LVESV and LVEF were assessed using the modified biplane Simpson’s equation in the apical four and two-chamber views.

Mitral inflow was assessed with pulsed-wave Doppler echocardiography from the apical four-chamber view.15 The Doppler beam was aligned parallel to the direction of flow, and a 1–2 mm sample volume was placed between the tips of mitral leaflets during diastole. From the mitral inflow profile, the E-wave (E) and A-wave (A) velocity and E/A ratio were measured.

Tissue Doppler imaging (TDI) was performed at the apical four-chamber view for the long-axis motion of the heart.15 16 Two-dimensional echocardiography with colour TDI was performed. The imaging angle was adjusted to ensure a parallel alignment of the sampling window with the myocardial segment of interest. Gain settings, filters, pulse repetitive frequency, sector size and depth were adjusted to optimise colour saturation. The frame rate was adjusted to greater than 100. Pulse Doppler sample volume was placed at the septal and lateral MV annulus to get the systolic velocity (s′-S and s′-L) and the early diastolic velocity (e′-S and e′-L), s′ was the mean value of s′-S and s′-L, e′ was the mean value of e′-S and e′-L and the ratio of E/e′ was calculated by the following formula: E/e′=(E/e′–S+E/e′–L)/2. At least three consecutive beats were stored, and the images were digitised and analysed offline by EchoPac software (EchoPac 6.3.6; Vingmed-General Electric) and the average value was taken.

Normal LV geometry was defined as RWT <0.42 and LVMI ≤95 g/m2 in women or LVMI ≤115 g/m2 in men.6 The presence of LVDDF was defined as LAVI ≥34 mL/m2 and e′-L ≤10 cm/s or e′-S ≤8 cm/s according to the recommendation from the American College of Cardiology.17

Statistical analysis

Continuous normally distributed data were displayed as mean±SD, and categorical data were shown as a total sample and proportion. The multivariate analysis selection criterion from the univariate analysis was p<0.05. Multivariate logistic and linear regression analysis was performed in two models with entry and retention at a significance level of 0.05. Model 1 was adjusted for age and gender, and model 2 was adjusted for covariates related to obesity, dyslipidaemia, diabetes, CHD, health-related behaviours and renal function (model 1+HR, BMI, physical activity, smoking and alcohol-drinking status, antihypertensive treatment, HbA1c level, serum LDL-C level, history of CHD and eGFR).

The participants were divided into two groups according to LV geometry and diastolic function. The effect of demographic data and conventional cardiovascular risk factors on LV geometry and LVDDF was analysed. The difference between the two groups was examined by independent t-test. Categorical variables were compared between the two groups by the X2 test. Multivariate logistic regression analysis was performed to investigate the independent relationship of various covariates and abnormal LV geometry and LVDDF.

To exam the effects of different ranges of SBP on LV geometry and diastolic function, all participants were divided on the basis of SBP values: group 1, SBP <120 mm Hg; group 2, SBP between 120 and 140 mm Hg; and group 3, SBP ≥140 mm Hg. Further multivariate logistic regression analysis was performed in two models. Model 1 was adjusted for age, gender, model 2 was adjusted for age, gender, HR, BMI, physical activity, smoking and alcohol drinking status, antihypertensive treatment, HbA1c level, serum LDL-C level, history of CHD and eGFR. In these two models, group 1 was set as the reference, and the OR and 95% CI of group 2 and group 3 were calculated.

The potential linear relationship between SBP and the echocardiographic parameters was analysed by univariate linear regression analysis in the whole population, and further multivariate linear regression analysis was performed in the above two models.

SPSS V.13.0 (SPSS) was used for the calculations. A p<0.05 was considered statistically significant.

Results

A total of 1051 subjects were included in the present analysis. As shown in table 1, the mean age of the participants was 56.6±9.3 years, and 37.9% were men. Among them, 561 (53.4%) individuals had hypertension, and the mean SBP and DBP were 133.4±19.2 and 78.5±10.7 mm Hg, respectively. The number of individuals with SBP <120 mm Hg, SBP between 12 and 140 mm Hg, and SBP ≥140 mm Hg was 278 (26.4%), 416 (39.6%) and 357 (34%), respectively. As for the echocardiographic data, the LVEDD was 47.2±4.3 mm, LVEDV was 69.6±15.9 mL and LVEF was 69.0%±6.4%. LVEF was less than 50% in 11 individuals, five of whom had old myocardial infarction. The LVMI of male and female participants was 77.9±17.9 g/m2 and 74.3±17.1 g/m2, respectively, and RWT was 35.1%±5.7%. The parameters of LV diastolic function were as follows: LAVI was 25.0±6.9 mL/m2, E was 69.2±17.6 cm/s and E/A ratio was 0.92±0.31. e′-S and e′-L were 6.9±1.9 cm/s and 9.8±2.7 cm/s, respectively, and the number of individuals with abnormal LV geometry and LVDDF was 191 (18.2%) and 91 (8.7%), respectively. Other clinical and echocardiographic characteristics are displayed in table 1.

Table 1

The clinical and echocardiographic data of the study population and the comparison between the groups divided by LV geometry and LVDDF

Association of SBP with LVDDF and abnormal LV geometry

As shown in table 1, participants with abnormal LV geometry and diastolic function were older (59.8±9.1 vs 55.9±9.2 years, and 60.1±10.4 vs 56.3±9.1 years, respectively, all p<0.01), had higher prevalence of hypertension and higher SBP (64.9% vs 50.8% and 138.4±19.3 vs 132.3±19.0 mm Hg, 69.2% vs 51.9% and 140.0±20.3 vs 132.8±18.9 mm Hg, respectively, all p<0.01), larger proportion of participants with higher SBP group and antihypertensive treatment (all p<0.05). The participants with abnormal LV geometry had higher BMI (27.3±3.8 vs 26.4±3.7 kg/m2, p<0.01), higher LDL-C (3.1±0.8 vs 2.9±0.8 mmol/L, p<0.01), lower eGFR (90.5±14.1 vs 93.5±14.4 mL/min/1.73 m2, p<0.01), less likely to be male (25.1% vs 40.7%, p<0.01) and less likely to be a current smoker and alcohol drinker (16.2% vs 24.5% and 15.7% vs 25.9%, respectively, all p<0.01). Other data were shown in table 1.

In the multivariate analysis (table 2), SBP was independently related to abnormal LV geometry and LVDDF (OR 1.01, 95% CI 1.00 to 1.02, and OR 1.02, 95% CI 1.01 to 1.03, respectively, all p<0.05). Other conventional risk factors independently related to abnormal LV geometry and LVDDF are shown in table 2.

Table 2

The multivariate analysis of the relationship between the risk factors and abnormal LV geometry and LVDDF

The effects of different ranges of SBP on LV geometry and diastolic function

The effects of different ranges of SBP on LV geometry and diastolic function were analysed in two models (figure 2). In the first model, the participants in group 3 had a higher risk odds of abnormal LVDDF and LV geometry than those in group 1 (OR 2.19, 95% CI 1.14 to 4.23 for LVDDF and OR 2.26, 95% CI 1.40 to 3.64 for abnormal LV geometry, all p<0.05). These differences remained statistically significant in the second model (OR 3.08, 95% CI 1.56 to 6.08 for LVDDF and OR 2.02, 95% CI 1.24 to 3.30 for abnormal LV geometry, all p<0.05).

Figure 2

The effects of different ranges of SBP on LV geometry and diastolic function, SBP <120 mm Hg was set as reference. Model 1was adjusted for age, gender, model 2 was adjusted for age, gender, HR, BMI, physical activity, smoking and alcohol drinking status, antihypertensive treatment, HbA1c level, serum LDL-C level, history of CHD and eGFR. BMI, body mass index; CHD, coronary heart disease; eGFR, estimated glomerular filtration rate; HbA1c, haemoglobin A1c; HR, heart rate; LDL-C, low-density lipoprotein cholesterol; LV, left ventricular; LVDDF, left ventricular diastolic dysfunction; SBP, systolic blood pressure.

Compared with the individuals of group 1, participants of group 2 also showed increased risk odds of abnormal LV geometry in both models (OR 1.92, 95% CI 1.19 to 3.08 in model 1 and OR 1.76, 95% CI 1.09 to 2.85 in model 2, all p<0.05) (figure 2). The participants of group 2 had greater LVMI than those of group 1 but lower than those from group 3 (table 3).

Table 3

The effects of different range of SBP on LV geometry and diastolic function

Although compared with the individuals of group 1, participants of group 2 did not show increased risk odds of LVDDF in both models (figure 2), the parameters of LV diastolic function (e′ and E/e′) were progressively worse in group 2 and group 3 in comparison with group 1 (table 3).

Linear association of SBP with cardiac remodeling and diastolic function assessed by echocardiography in the whole population

In the univariate analysis, SBP was positively correlated with LVMI, LVEDD, RWT, LAVI and E/e′ (all p<0.01) and was negatively correlated with e′ (p<0.01). In both models, after adjusting for the relevant parameters, SBP was still independently correlated with the echocardiographic parameters about cardiac remodelling (LVMI and LVEDD, p<0.01) and LV diastolic function (LAVI, e′, and E/e′, all p<0.01) (table 4).

Table 4

Associations of SBP with cardiac remodelling and diastolic function assessed by echocardiography in the whole population

We randomly selected 10 individuals for measurement of echocardiographic data by two experienced physicians, one of whom took the measurement again 1 day later and calculated the intraobserver and interobserver variability which were less than 4% and 5% in our echocardiographic laboratory.

Discussion

In this cross-sectional population-based study of general Chinese population, we assessed the relationship between the LV diastolic functional and structural changes and SBP, and we have two major findings.

First, SBP was independently related to the abnormal LV geometry and had an independent linear relationship with LVMI. Compared with individuals with SBP <120 mm Hg, participants with SBP between 120 and 140 mm Hg had higher LVMI, although this LV remodelling was less severe than what was seen in participants with SBP ≥140 mm Hg. LV remodelling is a well-known marker of cardiovascular risk,18 and BP is one of the major determinants of changes in LV geometry. As shown in this study, every 1 mm Hg increase in SBP related to 1% increase in the odds of abnormal LV geometry. In previous studies,9 19 prehypertensive patients were found to have higher LVM compared with normotensive subjects, which was consistent with our study. In our study, LV remodelling was based only on increasing LVMI, with no difference in RWT among different ranges of SBP and no difference in LVEDD between the individuals with SBP between 120 and 140 mm Hg and those with SBP ≥140 mm Hg. This finding was not consistent with previous research.9 The reason was unclear, but the participants in our study were older than those in the previous research (56.6±9.3 vs 40.3±8.1 years), and different age and sample size might be the possible explanation for the different results.

The second major finding was the effect of elevated SBP on subclinical LV diastolic function. Many previous studies8 9 20 21 suggested prehypertension was related to and accelerated the development of diastolic dysfunction of the heart. Our study showed that such association was continuous as the SBP increased. In our whole-study population, for every 1 mm Hg increase in SBP, the risk of diastolic dysfunction increased by 2% after adjusting for conventional cardiac risk factors. Further quantitative analysis of diastolic function (expressed as e′, E/e′) indicated that a dose–response relationship between SBP and LV diastolic parameters existed across the study population, from the lowest SBP to the highest SBP. Compared with individuals who had SBP <120 mm Hg, participants with SBP between 120 and 140 mm Hg had higher E/e′ and lower e′, although the impairment of LV diastolic function was less severe than what was seen in participants with SBP ≥140 mm Hg.

In a previous study,19 Santos et al found more frequent LVDDF in prehypertension than optimal BP. In the present study, there was no statistically significant difference in the prevalence of the diastolic dysfunction between the participants with SBP between 120 and 140 mm Hg and those with SBP <120 mm Hg. The reason was unclear, but the different definition of LVDDF might be the reason, and Santos et al used a more complex definition.

LV remodelling and LVDDF were the most critical factors in the pathogenesis of many cardiac diseases, such as heart failure with normal ejection fraction, and they started before the onset of symptoms. The newest guideline22 classified BP in adults into three categories: normal BP (SBP <120 mm Hg and DBP <80 mm Hg), elevated BP (SBP between 120 and 130 mm Hg and DBP <80 mm Hg), and hypertension (SBP ≥130 mm Hg or DBP ≥80 mm Hg). Our findings seemed to support this new classification. When the SBP exceeded the normal range (120 mm Hg), LV geometry ((LVMI) and LV diastolic function (expressed as e′, E/e′) were linearly, negatively correlated with SBP. From this perspective, more intensive BP control, especially controlling the SBP below 120 mm Hg, might reduce the impairment of LV structure and diastolic function. It also reminded us that in the population, elevated SBP alone was associated with impairment of LV structure and diastolic function, more attention should be paid to the control of SBP. Our findings were consistent with recent systolic blood pressure intervention trial (SPRINT) study.23

The strengths of the study include its large sample from the general Chinese population and its standardised acquisition of basic characteristics and echocardiography examinations. There are several limitations of this study: one limitation is that this is a cross-sectional study and cannot prove a causal relationship between the geometric/functional changes and SBP; another limitation is that this study was conducted in a Chinese population; studies that involve other ethnic populations are needed to replicate our results; third limitation is that we included self-reported hypertension into the definition of hypertension, a possibility of manipulating the BP readings exists, so end-digit preference for zero numbers and specific-value preference for readings just below the alert threshold exists.24

Conclusion

In the present population-based study, we found that in the Chinese population aged 35 years or older, SBP was independently related to abnormal LV geometry and LVDDF, SBP between 120 and 140 mm Hg was independently related to worse LV remodelling and diastolic function. Our findings indicate the potential benefit of intensive SBP control. But what is the threshold for SBP control, and whether the benefits of intensive SBP control on the heart will lead to improved long-term outcomes, all these issues warrant further investigation.

References

  1. 1.
  2. 2.
  3. 3.
  4. 4.
  5. 5.
  6. 6.
  7. 7.
  8. 8.
  9. 9.
  10. 10.
  11. 11.
  12. 12.
  13. 13.
  14. 14.
  15. 15.
  16. 16.
  17. 17.
  18. 18.
  19. 19.
  20. 20.
  21. 21.
  22. 22.
  23. 23.
  24. 24.

Footnotes

  • LY and LJ contributed equally.

  • Contributors LY, ZYi, CL, JZ and SZ designed the study. LY, LJ, ZL, ZF, TZ, ZYa, CW, BH and WH collected the data. LY and ZD analysed the data. All authors participated actively in the writing of the manuscript and approved the final draft.

  • Funding This work was supported by the Fundamental Research Funds for the Central Universities, CAMS Innovation Fund for Medical Sciences (CIFMS) (no.2017-12M-3-008, 2016-I2M-4-003 and no. 2016-12M-3-011), Peking Union Medical College (PUMC) Youth Fund, the National Natural Science Foundation of China (81500306) and the National Key Research and Development Program of China (grant no. 2016YFC0901500 and 2016YFC0901502).

  • Competing interests None declared.

  • Patient consent for publication Not required.

  • Ethics approval Ethics approval to conduct the study was obtained from by the Institutional Review Board of the Institute of Basic Medical Sciences, Chinese Academy of Medical Sciences (approval No. 009-2014, 031-2017).

  • Provenance and peer review Not commissioned; externally peer reviewed.

  • Data availability statement All data relevant to the study are included in the article or uploaded as online supplementary information.