Skip to main content

Establishing an appropriate Z score regression equation for Chinese pediatric coronary artery echocardiography: a multicenter prospective cohort study

Abstract

Background

Z score utility is emphasized in classifying coronary artery lesions in Kawasaki disease patients. The present study is the largest such multicenter Chinese pediatric study about coronary artery diameter reference values and Z score regression equation to date. It is useful in Chinese pediatric echocardiography.

Methods

A multicenter cohort was assembled, which consisted of 852 healthy children between 1 month and 17 years of age, ten children were excluded because their ultrasound images were not clear, or lost in following up. Diameters of the right coronary artery, left coronary artery, and left anterior descending coronary artery were assessed using echocardiography. Data were body surface area (BSA)-corrected using BSA calculated via either the Stevenson BSA formula or the Haycock BSA formula. Coronary artery diameter reference values and Z score regression equations were established for use in the Chinese pediatric population.

Results

No difference was observed between coronary artery diameter data corrected using BSAste or BSAhay. Of the five assessed regression models, the exponential model exhibited the best fit and was therefore selected as the basis for derivation of the SZ method. When comparing Z scores, those produced by the SZ method conformed to the standard normal distribution, while those produced by the D method did not. In addition, there was a statistically significant difference between Z scores produced by the SZ and D methods (P < 0.05).

Conclusions

Coronary artery diameter reference values for echocardiography were successfully established for use in the Chinese pediatric population, and a Z score regression equation more suitable for clinical use in this population was successfully developed.

Peer Review reports

Background

Echocardiography is commonly used to evaluate pediatric patients for coronary artery disease (CAD). Kawasaki disease is the major cause of childhood-acquired CAD. Accurate pediatric reference values for echocardiographic monitoring of changes in coronary artery diameter contribute to establishing diagnosis and prognosis. The most recent American Heart Association guidelines recommend body surface area (BSA)-corrected coronary artery diameter as the gold standard, and emphasize Z score utility in evaluating coronary artery injury risk and in classifying coronary artery lesions in Kawasaki disease patients. The American Society of Echocardiography also recommends the use of Z scores in pediatric cardiology [1]. These scores indicate the distribution of measurements about the mean within a healthy population, and facilitate comparison of datasets exhibiting differential means and distributions (such as datasets deriving from children of different ages and sizes).

Coronary artery Z score regression equations have been derived for multiple healthy pediatric populations, including those of the United States [2], Canada [3, 4], Japan [5], Singapore [6], Turkey [7], Korea [8], and China [9]. With the advent of reference value predictive models, many manufacturers have incorporated such models into echocardiography equipment-integrated software in order to allow rapid computation of coronary artery measurement Z scores. However, whether models developed in North America and/or using different BSA formulae are suitable for use in the Chinese pediatric population remains unclear.

The present study therefore aimed to: (1) establish coronary artery diameter reference values in a healthy Chinese pediatric population, and (2) compare suitability for use in this population of a Z score regression equation derived from these data with that of one derived from data of mostly-Caucasian pediatric population [4].

Methods

Materials

A multicenter healthy Chinese pediatric cohort was assembled. Pediatric cardiologists and ultrasound experts from six provinces in Southern and Northern China participated in this study. All researchers reached a consensus and unified the inclusion and exclusion criteria. The cohort consisted of 862 participants between July 1, 2016, and December 30, 2019.The children within this cohort were followed up for 6 months. Ten children were excluded because their ultrasound images were not clear, or lost in following up. There were 546 males and 296 females. Height (H), weight (W) and systolic / diastolic blood pressure (SBP / DBP) were recorded. At each center, echocardiographic data of each child were evaluated.

Inclusion criteria

  1. (1)

    Age from 1 month to 18 years old.

  2. (2)

    No cardiovascular history.

  3. (3)

    Blood pressure was normal [10].

  4. (4)

    Physical examination showed no abnormality of cardiovascular system.

  5. (5)

    Heart structure were normal by clinical examination and echocardiography.

  6. (6)

    Echocardiography examination for the following reasons: asymptomatic heart murmur, chest pain, suspicious abnormal chest X-ray or ECG performance, family history of congenital heart disease.

Exclusion criteria

  1. (1)

    Congenital heart disease (The presence of a patent foramen ovale in infants was considered normal).

  2. (2)

    History of KD or surgical catheterization for heart disease.

  3. (3)

    Myocardial disease.

  4. (4)

    Chromosomal abnormalities (i.e. trisomy 21, Turner syndrome, noonan syndrome, Marfan syndrome, Williams syndrome, 22q11.2 deletion, etc.).

  5. (5)

    Major potential diseases (such as malignant tumor, leukemia, metabolic or endocrine diseases, neuromuscular diseases, chronic renal failure, etc.).

  6. (6)

    Professional athletes.

This study was approved by the ethics committee of Shenzhen children’s hospital. Signed the informed consent by children’s parents were obtained.

Equipment and methods

The age, weight, and height of each patient were recorded at the time of the echocardiography evaluation. Data were body surface area (BSA)-corrected using BSA calculated via either BSA Stevenson formula (BSAste) [11] or the BSA Haycock formula (BSAhay) [12].

$$ \mathrm{BSASte}\left[{\mathrm{m}}^2\right]=0.0061\times \mathrm{H}\left[\mathrm{cm}\right]+0.0128\times \mathrm{W}\left[\mathrm{kg}\right]-0.1529 $$
$$ \mathrm{BSAHay}\left[{\mathrm{m}}^2\right]=0.024265\times \mathrm{H}{\left[\mathrm{cm}\right]}^{0.3964}\times \mathrm{W}{\left[\mathrm{kg}\right]}^{0.5378} $$

The studies were performed using GE Vivid (GE Healthcare, Milwaukee, WI) or Philips iE33 echocardiograph (Philips Medical Systems, Bothell, WA) echocardiography system. All Echocardiographic data were digitally recorded, allowing offline measurements. The frequency of the probe used should be as high as possible, the frame frequency > 60fps. Diameters of the right coronary artery (RCA), left coronary artery (LCA), and left anterior descending (LAD) coronary artery were assessed using echocardiography [13], which were defined as raw values (Fig. 1).

Fig. 1
figure 1

Short-axis echocardiographic view used for coronary artery diameter measurements. A RCA, right coronary artery; B LCA, left coronary artery; LAD, left anterior descending coronary artery

Statistics analysis

General data: SPSS 11.0 software (SPSS, Inc., Chicago, IL, UnitedStates) was used for statistical analysis. Continuous variables (age, height, weight, blood pressure, heart rate, BSA) were expressed as mean (SD) and / or medians.

The Anderson darling test is used to determine whether the data conform to the normal distribution. The student’s t test was used to compare the difference between two groups of data which conformed to a standard normal distribution; Wilcoxon signed-rank test was used to compare the difference between two groups of data which did not. P values< 0.05 were considered statistically significant.

Regression analysis

Five regression models (linear, logarithmic, power, exponential, and squareroot) were assessed for goodness of fit using R language,based on adjustedR2, residual standard error (RSE), and Akaike information criterion (AIC) values. The best-performing model (maximum adjusted R2and minimumRSEandAIC) was carried forward to derive a novelreference regression equation, termed the ShenZhen (SZ) method.

Comparison of two regression equations

Both the SZmethod and the North AmericanDallaire(D)method [4] were applied to the same set of raw data, and resultantZ scores were compared.

Results

Cohort demographic data

Overall, a pediatric cohort of 842 healthy subjects was assembled, comprising 67% males and 33% females, and spanning an age range of 1 month to 17 years (Table 1).

Table 1 Age demographic information of healthy Participants

Echocardiography-determined coronary artery diameters corrected using BSAs calculated via two different formulae.

Both BSAste- and BSAhay-corrected coronary artery diameters were normally distributed (P > 0.05) and did not differ significantly from each other (P > 0.05) (Table 2).

Table 2 Comparison of coronary artery diameters corresponding to BSAste and BSAhay

Regression model performance comparison

The R project for statistical computing was used to compare five regression models, including linear, logarithmic, power, exponential, and square root. Of these, the exponential model exhibited the best fit and was therefore selected as the basis for derivation of the SZ method. Exponential model performance based on BSAste-corrected values was as follows: lnyLCA = 1.0042 + 0.3494lnx (R2 = 0.678, RSE = 0.114, AIC = − 932.57), lnyLAD = 0.8013 + 0.3448lnx (R2 = 0.641, RSE = 0.120, AIC = − 849.8), and lnyRCA = 0.8107 + 0.3098lnx (R2 = 0.4806, RSE = 0.150, AIC = − 604) (Table 3). Exponential model performance based on BSAhay-corrected values was as follows: lnyLCA = 1.0084 + 0.3740lnx (R2 = 0.682,RSE = 0.115, AIC = − 938.7), lnyLAD = 0.8063 + 0.3697lnx (R2 = 0.650, RSE = 0.119, AIC = − 858.8), and lnyRCA = 0.8040 + 0.3291lnx (R2 = 0.4815, RSE = 0.150, AIC = − 598.2) (Table 3). Because the exponential model based on BSAhay-corrected values yielded the maximum adjusted R2 and minimum RSE and AIC, it was selected to construct a novel regression equation, herein referred to as the SZ method.

Table 3 Reference regression equation of coronary artery Echocardiography based on the BSAste and BSAhay formulas

Relationship of coronary artery diameter to BSAhay, and Z score distribution

Coronary artery diameters increased non-linearly with BSAhay (Fig. 2A). Values were normally distributed (Fig. 2B), and the normal distribution was standard (P > 0.05) (Fig. 2C).

Fig. 2
figure 2

A1-A3: Trend charts and predicted mean curves for LCA, LAD, and RCA. The black line is the predicted value, and the red line is the prediction interval with 95% confidence. B1 ~ B3: studentized residual diagrams for LCA, LAD, and RCA. The black line is the normal curve, and the red line is the estimated curve. C1 ~ C3: Z score scatter plots for LCA, LAD, and RCA, which represent the fitted values of the regression models, along with upper and lower bounds, i.e. 2 and − 2 z-scores, respectively

Comparison of SZ versus D method Z scores

When comparing Z scores resulting from the SZ versus D methods, SZ-derived Z scores conformed to a standard normal distribution, while D-derived Z scores (within certain BSAhay ranges) did not, and the difference between Z scores resulting from the two methods was significant (P < 0.05) (Table 4). Based on raw LCA diameter data, SZ method Z scores were larger than those of the D method when BSA ≤ 1.0 m2, but smaller than those of the D method when BSA > 1.0 m2 (P < 0.05). Similarly, based on raw LAD artery diameter data, SZ method Z scores were larger than those of the D method when BSA ≤ 1.3 m2, but smaller than those of the D method when BSA > 1.3 m2 (P < 0.05). Finally, based on raw RCA artery diameter data, SZ method Z scores were larger than those of the D method when BSA ≤ 0.7 m2, but smaller than those of the D method when BSA > 0.7 m2 (P < 0.05).

Table 4 Based on the same set of raw data, two groups result of Z scores calculated by SZ method and D method respectively

D method Z scores

Within certain BSAhay ranges, D method Z scores were non-normally distributed (Fig. 3).

Fig. 3
figure 3

A1 ~ A3: D method Z score scatter plots for LCA, LAD, and RCA with 2Z score boundaries. B1 ~ B3: studentized residual diagrams for LCA, LAD, and RCA. The black line is the predicted value, and the red line is the prediction interval with 95% confidence

Performance of the SZ versus D methods

Within certain BSAhay ranges, the SZ and D methods predicted different values (P < 0.05) (Fig. 4). Specifically, D method-predicted LCA diameter means were larger than those predicted by the SZ method when BSAhay> 0.9 m2, D method-predicted LAD artery diameter means were smaller than those predicted by the SZ method when BSAhay ≤1.2 m2, and D method-predicted means of both diameters were larger than those predicted by the SZ method when BSAhay> 0.6 m2.

Fig. 4
figure 4

A, B, C Means±2Z scores predicted by the SZ and D methods for LCA, LAD, and RCA. The solid black line represents values predicted by the SZ method, the dotted blue line represents values predicted by the D method, and the solid red line represents the BSA in which the two predicted mean curves intersect. The red vertical line indicates the point where the SZ and D models yield the same z-score prediction for that specific BSA

Discussion

Because China is large and exhibits regional diversity in individual lifestyle, height, and weight, a multicenter study was conducted to more accurately ascertain pediatric coronary artery diameter reference values. The present study is the largest such multicenter Chinese pediatric study to date, and includes coronary artery measurement data from six regions in Northern and Southern China.

Study findings demonstrate that coronary artery diameters of a healthy Chinese pediatric population increase non-linearly with an increase in BSA, and differ significantly from those of the Dallaire cohort. Specifically, Chinese pediatric coronary artery diameters are larger at lower BSA ranges, but smaller at higher BSA ranges, than those of North American children. In addition, relative to the D method, the SZ method mean reference value-predictive regression model exhibits lower Z scores at lower BSA ranges and higher Z scores at higher BSA ranges. This suggests that the SZ method may be more accurate than the D method in the Chinese pediatric population.

Previous studies suggested the existence of a linear correlation between equal variance of coronary artery diameter and BSA [14,15,16,17]. However, multiple larger studies have found a non-linear correlation between coronary artery diameter and BSA in healthy pediatric populations [2,3,4,5,6,7,8,9, 18]. Coronary artery diameter standard deviation from the mean differs within distinct BSA ranges, resulting in population-inappropriate linear regression methods introducing bias during Z score calculation. Therefore, researchers have sought to identify regression models able to more accurately represent the non-linear relationship between coronary artery diameter mean and standard error. Specifically concerning heteroscedastic non-linear relationships, prior studies have evaluated the goodness of fit of various regression models, including quadratic, cubic polynomial, logarithmic, exponential, and square root [2,3,4,5,6,7,8,9].

In the present study, an exponential regression model exhibited the best fit (based on maximum R2 and minimum RSE and AIC) and the standard test of normality was therefore applied to model residuals and Z scores. Results indicate that the model performs reliably. The present study therefore contributes to establishment of Chinese pediatric coronary artery diameter reference values, and provides data which may help overcome the lack of universality associated with single-center studies.

Up to 25 BSA calculation methods are available, including the Stevenson, Haycock, and Du Bois EF formula [11, 12, 19, 20]. The American Society of Echocardiography guidelines for quantitative pediatric echocardiography recommends use of the Haycock formula when calculating Z scores for cardiovascular structure measurements [21, 22]. The Haycock formula may also be a better estimator of BSA for smaller children [4, 22], although different BSA formulae do not result in different model Z scores in the present study. However, the Stevenson formula is commonly used within the Chinese healthcare system. Therefore, in order to determine the impact of formula choice on model Z scores, the present study derived regression models from a single original dataset (using BSA calculated via either the Stevenson or Haycock formulae), and compared resultant Z scores. Because Z scores did not differ significantly, either BSA formula is appropriate for use in quantitative evaluation of echocardiography data in a Chinese pediatric population. The present study ultimately incorporated the Haycock formula into the SZ method regression model. This is consistent with the approach of the D method, which also uses the Haycock formula, and facilitates comparison of these two methods.

In comparing Z scores resulting from the SZ and D methods, it was found that the SZ method produced Z scores closer to zero. This indicates that the SZ method may be more suitable than the D method for use in the Chinese pediatric population. Furthermore, predicted coronary artery diameter mean values and standard deviations differed between North American and Chinese pediatric populations, and the predicted mean value curve provided by the D method had a higher gradient than that provided by the SZ method. When BSA is low, predicted mean values for North American children were lower than those for Chinese children, and when BSA is high, predicted mean values for North American children were higher than those for Chinese children. This indicates that using predictive regression models based on North American pediatric reference values may over- or underestimate coronary artery measurement Z scores in the Chinese pediatric population, which would negatively impact CAD diagnosis and treatment.

We suggest that the D method regression model used for the North American pediatric population is unsuitable for use in the Chinese pediatric population. Factors such as region and race should be taken into account when incorporating automatic calculation functions into echocardiographic equipment and when selecting predictive regression models for clinical applications. The SZ method regression model established during the present study may be more suitable for use in the Chinese pediatric population.

However, we acknowledge certain study limitations. For example, determination coefficients (R2 values) are marginally lower than those determined by previous studies, especially for the RCA [4, 9]. Determination coefficients for the LCA, LAD artery, and RCA obtained using the D method were higher than those obtained using the SZ method. Two possible reasons may account for this observation: differences between multiple participating study sites, or an unbalanced cohort gender ratio. Furthermore, coronary artery diameters were too small to avoid errors in measurement when enlarging images, and it is challenging to ensure inter-site consistency during a multicenter study. Mitigation strategies for such limitations should be considered when designing future studies.

Furthermore, Study cohort was skewed toward male and there are few subjects older than 12 years of age which limit generalization of this data to children > 12 years of age.

Conclusions

This study established new coronary artery diameter reference values for echocardiography and a Z score regression equation that hold clinical value for the evaluation of Chinese pediatric populations. A Z score regression equation derived from data of a Caucasian pediatric population may be not suitable for use in chinese population.

Availability of data and materials

The data that support the findings of this study are available from the authors upon reasonable request.

Abbreviations

BSA:

Body surface area

BSAste:

BSA Stevenson

BSAhay:

BSA Haycock

CAD:

Coronary artery disease

H:

Height

W:

Weight

SBP / DBP:

Systolic / diastolic blood pressure

RCA:

Right coronary artery

LCA:

Left coronary artery

LAD:

Left anterior descending

RSE:

Residual standard error

AIC:

Akaike information criterion

SZ:

ShenZhen

D:

Dallaire

References

  1. McCrindle BW, Rowley AH, Newburger JW, et al. Diagnosis, treatment, and long- term management of Kawasaki disease: a scientific statement for health professionals from the American Heart Association. Circulation. 2017;135(17):e927–99.

    Article  Google Scholar 

  2. Olivieri L, Arling B, Friberg M, et al. Coronary artery Z score regression equations and calculators derived from a large heterogeneous population of children undergoing echocardiography. J Am Soc Echocardiogr. 2009;22(2):159–64.

    Article  Google Scholar 

  3. McCrindle BW, Li J, Minich L, Colan S, Atz A, Takahashi M, et al. Coronary artery involvement in children with Kawasaki disease: risk factors from analysis of serial normalized measurements. Circulation. 2007;116:174–9.

    Article  Google Scholar 

  4. Dallaire F, Dahdah N. New equations and a critical appraisal of coronary artery Z scores in healthy children. J Am Soc Echocardiogr. 2011;24(1):60–74.

    Article  Google Scholar 

  5. Kobayashi T, Fuse S, Sakamoto N, et al. A new Z-score curve of the coronary arterial internal diameter using the lambda-mu-sigma method in a pediatric population. J Am Soc Echocardiogr. 2016;29(8):794–801.e29.

    Article  Google Scholar 

  6. Tan TH, Wong KY, Cheng TK, et al. Coronary normograms and the coronary-aorta index: objective determinants of coronary artery dilation. Pediatr Cardiol. 2003;24:328–35.

    Article  CAS  Google Scholar 

  7. Karagol BS, Orun UA, Zenciroglu A, et al. The diameter of coronary arteries in healthy newborns at birth,1 and 6 months of ages. J Matern Fetal Neonatal Med. 2012;25(12):2729–34.

    Article  Google Scholar 

  8. Ryu K, Yu JJ, Jun HO, et al. A comparative study of established Z score models for coronary artery diameters in 181 healthy Korean children. Kor J Pediatr. 2017;60(11):373–8.

    Article  Google Scholar 

  9. Zhang YQ, Chen SB, Huang GY, et al. Coronary artery indexed diameter and z score regression equations in healthy Chinese Han children. J Clin Ultrasound. 2015;43(1):39–46.

    PubMed  Google Scholar 

  10. Mi J, Wang TY, Meng LH, et al. Development of blood pressure reference standards for Chinese children and adolescents (in Chinese). Chin J Evid Based Pediatr. 2010;5(1):4–14.

    Google Scholar 

  11. Stevenson PH. Height weight surface formula for the estimation or surface area in the Chinese subjects. Chin J Physiol.1937;3:327–30.

  12. Haycock GB, Schwartz GJ, Wisotsky DH. Geometric method for measuring body surface area: a height-weight formula validated in infants, children, and adults. J Pediatr. 1978;93:62–66.

  13. Fuse S, Kobayashi T, Arakaki Y, et al. Standard method for ultrasound imaging of coronary artery in children. Pediatr Int. 2010;52:876–82.

  14. Oberhoffer R, Lang D, Feilen K. The diameter of coronaryarteries in infants and children without heart disease. Eur. JPediatr. 1989;148:389–92.

  15. DeZorzi A, Colan SD, Gauvreau K, et al. Coronary dimensions may be misclassified asnormal in Kawasaki disease. J. Pediatr. 1998;133:254–8.

  16. Arjunan K, Daniels SR, Meyer RA, et al. Coronary artery caliber in normal children andpatients with Kawasaki disease but without aneurysm: Anechocardiographic and angiographic study. JAmCollCardiol. 1986;8(5):1119–24.

  17. Kurotobi S, Nagai T, Kawakami N, Sano T. Coronary diameter in normal infants, children and patients with Kawasaki disease. PediatrInt. 2002;44:1–4.

  18. Lin MT, Chang CH, Hsieh WC, et al. Coronary Diameters in Taiwanese Children Younger than 6 Years Old: Z-Score Regression Equations Derived from Body Surface Area. Acta Cardiol Sin. 2014;30(4):266–73.

  19. Redlarski G, Palkowski A. Krawczuk. Body surface area formulae: an alarming ambiguity. M Sci Rep. 2016;6:27966.

    Article  CAS  Google Scholar 

  20. Du Bois D, Du Bois EF. A formula to estimate the approximate surface area if height and weight be known. Arch Intern Med. 1916;17(6):863–71.

    Article  Google Scholar 

  21. Lopez L, Colan S, Stylianou M, et al. Relationship of echocardiographic Z scores adjusted for body surface area to age, sex, race, and ethnicity: the pediatric heart network Normal echocardiogram database. Circ Cardiovasc Imaging. 2017;10(11):e006979.

    Article  Google Scholar 

  22. Lopez L, Colan SD, Frommelt PC, et al. Recommendations for quantification methods during the performance of a pediatric echocardiogram: a reportfrom the pediatric measurements writing Group of the American Society of echocardiography pediatric and congenital heart disease council. JAmSocEchocardiogr. 2010;23:465–95.

    Google Scholar 

Download references

Acknowledgments

We would like to thank the children and parents included in this study. Thanks to the continuous and excellent support of all authors and colleagues. Thanks to Editage (www.editage.cn) for English language editing.

Funding

National Natural Science Foundation of China (62071309).

Author information

Authors and Affiliations

Authors

Contributions

Bei Xia had primary responsibility for the study, protocol development, patient enrollment and outcome assessment and for writing the manuscript. Shu-min Fan contributed to data analyses, outcome assessment and manuscript writing. Wei Yu, Shu-Bao Chen, Qing-Hua Liu, Wen-juan Chen, Shan-liang Zhu, Mei Jin, Jie-Qing Min, Yang Xu were responsible for patient enrollment and data collection. Wei-xiang Liu, Zhi-xia Wu performed final data analyses. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Bei Xia.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Ethical committee for medical research at Shenzhen Children’s Hospital. Written consent from each parent was retrieved. Reference Number: 2020012 158–31.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Fan, Sm., Xia, B., Liu, Wx. et al. Establishing an appropriate Z score regression equation for Chinese pediatric coronary artery echocardiography: a multicenter prospective cohort study. BMC Pediatr 21, 429 (2021). https://doi.org/10.1186/s12887-021-02877-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12887-021-02877-9

Keywords