Article Text

Download PDFPDF

Association between ATP2B1 and CACNB2 polymorphisms and high blood pressure in a population of Lithuanian children and adolescents: a cross-sectional study
  1. Sandrita Simonyte,
  2. Renata Kuciene,
  3. Virginija Dulskiene,
  4. Vaiva Lesauskaite
  1. Institute of Cardiology of the Medical Academy, Lithuanian University of Health Sciences, Kaunas, Lithuania
  1. Correspondence to Dr Sandrita Simonyte; sandrita.simonyte{at}lsmuni.lt

Abstract

Objectives Recently, genome-wide associated studies have identified several genetic loci that are associated with elevated blood pressure and could play a critical role in intracellular calcium homeostasis. The aim of this study was to assess the associations of ATP2B1 rs2681472 and CACNB2 rs12258967 gene polymorphisms with high blood pressure (HBP) among Lithuanian children and adolescents aged 12–15 years.

Study design and participants This was a cross-sectional study of a randomly selected sample of 646 12–15-year-old adolescents who participated in the survey ‘The Prevalence and Risk Factors of HBP in 12–15 Year-Old Lithuanian Children and Adolescents (from November 2010 to April 2012)’. Anthropometric parameters and BP were measured. The participants with HBP were screened on two separate occasions. Subjects were genotyped ATP2B1 rs2681472 and CACNB2 rs12258967 gene polymorphisms using real-time PCR method.

Results The prevalence of HBP was 36.7%, significantly higher for boys than for girls. In the multivariate analysis, after adjustment for body mass index and waist circumference, boys with CACNB2 CG genotype, CACNB2 GG genotype and CACNB2 CG +GG genotype had higher odds of having HBP in codominant (adjusted OR (aOR)=1.92; 95% CI 1.16 to 3.18, p=0.011; and aOR=2.64; 95% CI 1.19 to 5.90, p=0.018) and in dominant (aOR=2.05; 95% CI 1.27 to 3.30, p=0.003) inheritance models. Girls carrying CACNB2 CG genotype and CACNB2 CG +GG genotype had increased odds of HBP in codominant (aOR=1.82; 95% CI 1.02 to 3.24, p=0.044) and in dominant (aOR=1.89; 95% CI 1.09 to 3.28, p=0.023) inheritance models. Furthermore, significant associations were found in additive models separately for boys (aOR=1.72; 95% CI 1.20 to 2.46, p=0.003) and girls (aOR=1.52; 95% CI 1.05 to 2.20, p=0.027). No significant association was found between ATP2B1 gene polymorphism and the odds of HBP.

Conclusions Our results indicate that CACNB2 gene polymorphism was significantly associated with higher odds of HBP in Lithuanian adolescents aged 12–15 years.

  • cardiac epidemiology
  • hypertension
  • molecular biology
  • paediatric cardiology

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 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

  • To our knowledge, this is the first study to determine the association of the ATP2B1 rs2681472 polymorphism with high blood pressure (HBP) in adolescents aged 12–15 years in the Baltic states.

  • Multivariate logistic regression analyses were conducted separately for boys and girls to evaluate the associations between ATP2B1 rs2681472 and CACNB2 rs12258697 gene polymorphisms and HBP under different inheritance models.

  • As our study was a cross-sectional study in its design, we cannot determine a cause–effect relationship.

Introduction 

Hypertension is one of the major risk factors of cardiovascular disease, and does not only affect middle-aged and elderly people, but is also increasing in prevalence among children and adolescents worldwide.1 2 A well-known fact is that overweight and obesity, like family history of hypertension and low socioeconomic status, are related to elevated blood pressure (BP) in childhood, which is an important predictor of hypertension in adult life.3 4 According to various epidemiological studies, the prevalence of elevated BP in children and adolescents ranges from 0.8% up to 5%.5–7 Studies conducted in Lithuania have shown a higher prevalence of an increased BP among 3–7-year-old children (21.4%)8 and a higher prevalence of prehypertension and hypertension among 12–15-year-old adolescents (12.6% and 22.5%, respectively).9 The mortality rate from cardiovascular diseases in Lithuania is one of the highest in Europe.10 Hypertension is a highly heritable trait, and around 40%–60% of individual differences in BP have a genetic basis.11 Therefore, early identification of not only potential environmental but also genetic risk factors for the development of hypertension is essential, as they are important diagnostic and prognostic molecular markers.

Genome-wide associated studies have identified several genetic loci that are associated with elevated systolic (SBP) or diastolic blood pressure (DBP),12–14 but the mechanism by which a certain gene variant influences the risk of hypertension remains unclear. In our study, we selected genes (ATP2B1 and CACNB2) that play a critical role in intracellular calcium homeostasis15 and could be also involved in the pathogenic processes of hypertension. The ATP2B1 gene encodes the plasma membrane calcium ATPase isoform 1, which removes bivalent calcium ions from eucaryotic cells against very large concentration gradients16 and is responsible for controlling the contraction and dilation of vascular smooth muscles.17 The CACNB2 gene encodes the beta-2 subunit of a voltage-gated calcium channel and is a member of high voltage-gated calcium channel genes. The beta-2 subunit could interact with a pore-forming subunit of the calcium channel and could modulate calcium channel activity and blood pressure.18

The aim of our study was to investigate the associations of ATP2B1 rs2681472 and CACNB2 rs12258967 gene polymorphisms with HBP among Lithuanian children and adolescents aged 12–15 years.

Methods

Study design and population

Blood pressure and anthropometric measurements performed in this study have been described in our previous publication.19

This was a cross-sectional study of children and adolescents who had participated in the baseline survey ‘The Prevalence and Risk Factors of HBP in 12–15 Year-Old Lithuanian Children and Adolescents (Study 1, 2010–2012)’ in Kaunas City and Kaunas District, which are the second-largest city and district in Lithuania.19 The baseline survey that was based on a two-stage cluster sampling design enrolled 7638 children and adolescents aged 12 to 15 years who at the time of the examination (from November 2010 to April 2012) attended gymnasiums or secondary schools of all classes (grades 6, 7, 8 and 9). Data on clinically verified health disorders were collected from the subjects’ medical records (Form No. 027–1/a). The present study examined data of a randomly selected sample of 646 participants (313 boys and 333 girls) aged 12–15 years who underwent BP and anthropometric examinations as well as a genetic examination of the saliva.

Measurements

Blood pressure measurement

BP was measured in the morning hours by a physician who was not wearing a white coat in a quiet environment. The adolescents were advised to avoid coffee, tea, energy drinks and physical exercises until the measurements were taken. Before the BP measurement, the participants were asked to sit still for 10 min. BP measurements were performed with an automatic BP monitor (Omron M6; Omron Halthcare Co., Kyoto, Japan) using a cuff of the appropriate size. BP was measured three times with a 5 min rest interval between the measurements, with the participant being in a sitting position. The average of three BP measurements was calculated and used in the analysis. According to ‘The Fourth Report on the Diagnosis, Evaluation, and Treatment of High Blood Pressure in Children and Adolescents’, normal BP (NBP) was defined as SBP and DBP less than the 90th percentile for sex, age and height, and HBP was defined as SBP and DBP ≥90th percentile for sex, age and height.20 Prehypertension was defined as SBP or DBP levels between the 90th and the 95th percentile for sex, age and height, or if BP levels were ≥120/80 mm Hg even if below the 90th percentile.20 All subjects with high BP (BP was in the ≥90th percentile) during the first screening underwent a second evaluation of BP measurements within 2–3 weeks.

Anthropometric measurement

All participants underwent anthropometric measurements while being barefoot and wearing only light clothes. Body weight and height of the participants were measured to the nearest 0.1 kg and 0.1 cm, respectively, by using a balance beam scale (Seca) and a portable stadiometer. The participants were grouped into three categories: normal weight, overweight, and obese, using the cut-off points of body mass index (BMI) by age and sex for children and adolescents proposed by the International Obesity Task Force (IOTF).21

Waist circumference (WC) was measured to the nearest 0.5 cm by using a flexible measuring tape (Seca) at a level midway between the lower rib margin and the iliac crest. Using the cut-off values of the percentiles of the WC according to the criteria of the National Health and Nutrition Examination Survey III,22 the participants were grouped into the categories on the basis of their WC: below the 75th percentile (normal waist value), 75th<90th percentile and ≥90th percentile. High waist value was defined as WC ≥75th percentile.

Both BP and anthropometric measurements were performed at the subjects’ schools by the same team of trained research personnel (physicians and research assistants).

Genetic analysis

Saliva samples were collected into tubes from each individual for DNA extraction. Genomic DNA was extracted using a commercial DNA isolation kit—the ‘Genomic DNA Purification Kit’ (Thermo Fisher Scientific, Lithuania) according to the manufacturer’s instructions.

ATP2B1 rs2681472 and CACNB2 rs12258967 gene polymorphisms in the subjects were genotyped using the real-time PCR technique with TaqMan allelic discrimination Assay-By-Design genotyping assays C_30872739_10 and C_31302374_10, according to the manufacturer’s instructions (Applied Biosystems, Foster City, California, USA). Allele-specific fluorescence was analysed on the ABI 7900HT Sequence Detection System with SDS V.2.1 (Applied Biosystems).

Statistical analysis

Statistical analysis was performed using the statistical software package SPSS V.20 for Windows. Categorical variables were expressed as numbers and percentage. The z-test was used to compare differences between the groups. Continuous variables were presented as mean values±SD. The normality of the distribution of the continuous variables was assessed by applying the Kolmogorov-Smirnov test. Comparisons between the groups were performed by applying the χ2 test, Student’s t-test and analysis of variance with Bonferroni post hoc test. The z-test was used to compare the proportions. The χ2 test was used for the assessment of the Hardy-Weinberg equilibrium (HWE) for the distribution of genotypes. Logistic regression analyses were used to test for the associations of ATP2B1 (rs2681472) and CACNB2 (rs12258967) polymorphisms with HBP under five inheritance models based on the Akaike information criterion (AIC) to calculate the OR (95% CI) for CACNB2 and ATP2B1 polymorphism. ORs were adjusted for BMI and WC. The data were analysed separately for boys and girls. P values <0.05 were considered statistically significant.

Patient and public involvement

The schoolchildren and their parents or guardians were not involved in the design and conduct of the study.

Results

The characteristics of the study population are given in table 1. Among 646 study subjects aged 12–15 years (mean±SD age: 13.27±1.14 years), 48.5% (n=313) were boys, and 51.5% (n=333) were girls. Boys were significantly heavier and taller, and had a significantly greater mean WC. They had significantly higher mean SBP and significantly lower mean DBP, compared with girls. There was no significant difference in mean age and BMI between these groups.

Table 1

Demographic, anthropometric and BP characteristics of the study subjects by sex

The characteristics of the study participants according to BP levels are presented in table 2. The overall prevalence of HBP was 36.7%. HBP was found more frequently in boys than in girls (45.7% vs 28.2%; p<0.001). Subjects aged 14 to 15 years were more likely to have HBP than subjects aged 12 to 13 years, but a significant difference was found only between the boys (56.3% vs 38.5%, p<0.05) (data not shown). The subjects with HBP had significantly higher mean values of weight, height, BMI and WC, compared with subjects with NBP. The prevalence of overweight, obesity and large WC (≥75th percentile) were much more common in participants with HBP than in the NBP group.

Table 2

Characteristics of the study population according to the BP level

The distribution of the analysed genotypes was in the HWE (p>0.05) (table 3). No significant differences in the frequencies of ATP2B1 and CACNB2 genotypes or alleles between boys and girls were found.

Table 3

Distribution of the ATP2B1 and CACNB2 genotypes and frequency of alleles in the study population

The characteristics of the study subjects according to the ATP2B1 genotypes are presented in table 4. Mean values of SBP and DBP were similar in groups with different genotypes separately for both sexes. Statistically significant differences between ATP2B1 genotypes with respect to the mean values of age, weight, height and WC were observed only among boys. The subjects with ATP2B1 AA genotype had higher mean values of anthropometric variables and SBP than ATP2B1 AG and GG genotype carriers.

Table 4

Characteristics of the population (means and SD) according to ATP2B1 genotypes

The characteristics of the study subjects according to CACNB2 genotypes are presented in table 5. Mean values of weight, height, BMI, WC and DBP were similar in groups with different genotypes for both sexes separately, except for age and SBP among the boys. Furthermore, boys with the CACNB2 GG genotype had the highest level of SBP.

Table 5

Characteristics of the population (means and SD) according to CACNB2 genotypes

The mean SBP and DBP levels of the subjects according to genotypes are given in table 6. No significant differences were observed comparing the mean values of SBP and DBP between the ATP2B1 genotypes (AG vs AA, GG vs AA or GG vs AG), the CACNB2 genotypes (CG vs CC, GG vs CC or GG vs CG) separately for both sexes within each BP group.

Table 6

Distribution of ATP2B1 and CACNB2 genotypes according to systolic and diastolic blood pressure in the study population

The boys who carried one or two copies of the CACNB2 G allele were by 1.64 times more likely to have HBP than boys who had one or two copies of the CACNB2 C allele (table 7).

Table 7

Associations between ATP2B1 and CACNB2 polymorphisms and high blood pressure by sex (allelic models)

The results of multivariate logistic regression analyses are presented in table 8. To examine the associations with HBP, we calculated adjusted ORs (aORs) adjusted by the BMI and WC for boys and girls separately. No significant association between ATP2B1 genotypes and HBP was found either for boys or girls.

Table 8

Associations between ATP2B1 and CACNB2 polymorphisms and high blood pressure by sex (multivariate analyses)

The results showed that boys who carried the CACNB2 CG genotype, the CACNB2 GG genotype and the CACNB2 CG +GG genotypes had higher odds of having HBP than carriers of the CACNB2 CC genotype did (in the codominant model: aOR=1.92; p=0.011 and aOR=2.64; p=0.018; in the dominant model: aOR=2.05; p=0.003, respectively). In girls, carriers of the CACNB2 CG genotype and the CACNB2 CG +GG genotypes had significantly higher odds of HBP than girls with the CACNB2 CC genotype did (in the codominant model: aOR=1.82; p=0.044, and in dominant model: aOR=1.89; p=0.023). In addition, significant associations were found in additive models separately for both sexes.

Discussion

Though hypertension in adolescents is less common than in adults, it has become a major concern in some countries over the past decade because of the evidence suggesting that HBP tracks from childhood to adulthood.23 Children become hypertensive in adult years, particularly if they are obese as children or become obese as young adults.3 4 In adults, primary hypertension is more common than secondary hypertension in the paediatric population.24 In order to promote cardiovascular health in adults, it is strongly recommended to measure BP in all children ≥3 years of age.25 It was examined that the model of prediction of adult hypertension consists of known physical and environmental childhood risk factors, family history of hypertension and novel genetic variants.4 However, new findings of genome-wide association studies on adult hypertension are still little reflected in paediatric research. In the present study, we examined two newly identified genes together with measured hypertension risk factors among Lithuanian children and adolescents aged 12–15 years.

The prevalence of hypertension in the paediatric population varies because BP relates to sex, age and height during childhood.20 26 27 Nevertheless, many epidemiological studies reported the increasing prevalence of hypertension (from 11% to 22%) in children and adolescents between the ages of 3 and 18 years.1 28–31 In our study, the prevalence of HBP was 36.7%, and was significantly higher in boys than in girls. Our results are consistent with those of other studies showing sex differences in BP.32 33 Many studies have shown that obesity was associated with elevated BP.23 34–36 In this study, we found that the prevalence of overweight, obesity and high WC (≥75th percentile) was much more common in participants with HBP than in the NBP group, which is in line with the findings of other researchers.

There is no doubt that the genetic element in the evaluation of hypertension is an important factor. However, the pathophysiological mechanisms of the association between gene polymorphisms and variations of BP are poorly understood. We examined ATP2B1 rs2681472 and CACNB2 rs12258967 variants located in an intron region of these genes. According to numerous Genome-wide association studies (GWAS) reports,12–14 ATP2B1 and CACNB2 are involved in calcium transportation. In addition to well-known biological functions, calcium triggers muscle contraction and is the second messenger of hormones and growth factors.37 A defect in the regulation of calcium and calcium signalling plays an important role in hypertension-associated vascular dysfunction.38 The greatest amount of cellular calcium resides in the sarcoplasmic reticulum and mitochondria, resulting in its lower intracellular than extracellular concentration.39 After depolarisation, calcium enters cardiomyocytes via L-type calcium channels. It activates ryanodine receptors on the sarcoplasmic reticulum to release calcium, and triggers muscle contraction. As the contraction ends, intracellular calcium returns to the sarcoplasmic reticulum via calcium ATPase or plasma membrane Ca2+ ATPase isoform 1 (PMCA1), which is encoded by ATP2B1.40 The plasma membrane calcium/calmodulin-dependent ATPase is involved in calcium pumping from the cytosol to the extracellular compartment41 and regulates the homeostasis of cellular calcium levels, which is important in controlling the contraction and dilatation of vascular smooth muscles.42

CACNB2 encodes the intracellular beta-2 subunit of a calcium channel (voltage-dependent L-type calcium channel), which is a member of the high voltage-gated calcium channel genes and serves as a target of calcium channel blockers.43 These blockers, also known as calcium antagonists, act both on the cardiac tissue and on vascular smooth muscles. In cardiac tissues, a decrease in calcium levels directly results in a decrease of the force of the contraction of the heart, and thus a decrease in blood pressure. It was suggested that the beta-2 subunit could interact with the pore-forming alpha-1 subunit to modulate activity through conformational changes in the calcium channel. Mutations in genes could result in channel activation at more hyperpolarised membrane potentials, implicating an increased calcium influx in disease pathogenesis.44 This may explain why mutations in CACNB2 can influence intracellular calcium homeostasis and alter blood pressure.18 45 46

Studies have shown that CACNB2 was associated with DBP,12 13 23 47 with systolic pressure,13 48–52 mean arterial pressure53 and hypertension.13 49–52 The Kaunas Cardiovascular Risk cohort study found that CACNB2 CC carriers had the highest values of DBP in childhood. However, the logistic regression analysis failed to show significant associations between CACNB2 (rs12258967) and the risk of hypertension.23 CACNB2 polymorphism (rs4373814) was found to be significantly associated with a decreased risk (OR=0.70, 95% CI 0.51 to 0.95) of hypertension in the Han Chinese population.51 The association of the CACNB2 polymorphism (rs11014166) with hypertension (OR=0.79, 95% CI 0.65 to 0.97) was also observed in a study of China.52 A CACNB2 intronic single-nucleotide polymorphism (SNP), rs1571787, had the most significant association with pulse pressure (minor allele frequency (MAF)=0.27, p=0.0003) in a large cohort study.54 Our data showed that boys with the CACNB2 GG genotype had the highest level of SBP. The dominant model of CACNB2 adjusted for the BMI and WC showed that carriers of CG +GG genotypes had a significantly higher odds of HBP than carriers of the CACNB2 CC genotype did (in boys and girls: aOR=2.05; p=0.003 and aOR=1.89; p=0.023, respectively).

It remains unclear how ATP2B1 could result in an increased risk of hypertension, but a significant association of a variant of the ATP2B1 gene with blood pressure and hypertension was identified among Asian and European adults in GWAS.12 13 55 56 The Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium GWAS analysis showed that for rs2681472 in ATP2B1, the strongest signal for diastolic pressure – p=3.7×10-8 and with an OR for hypertension of 1.17 per risk allele.13 A meta-analysis involving 15 909 cases and 18 529 controls confirmed that rs2681472 was significantly associated with the risk of hypertension in East Asians (OR=1.18, 95% CI 1.10 to 1.27, p<0.0001).17 Another study among Chinese children confirmed a significant association of ATP2B1 rs17249754 with the risk of hypertension (allelic OR=1.25, 95% CI 1.08 to 1.44, p=0.003).57 However, five studies conducted in Korean,58 Chinese59 and Japanese populations16 60 61 yielded inconsistent results. Our study found no significant association of ATP2B1 with BP or with the risk of hypertension.

BP could be increased during puberty, which is related to rapid physical growth, the endocrine system and hormonal changes.62 Epidemiological studies have demonstrated that adolescents’ age was significantly associated with HBP.63–65 A cross-sectional study in China including children and adolescents aged 5–18 years reported that subjects aged 12–14 and 15–18 years had a significantly increased risk of prehypertension, compared with those aged 5–8 years.63 Data from the Healthy Kids Project showed that 12–17-year-old both boys and girls separately had a greater risk for prehypertension, hypertension and prehypertension/hypertension compared with 5–11-year-old participants.64 Results from the ‘Beijing children and adolescents (BP) study’ showed that adolescents with hypertension during puberty had a higher risk of hypertension in adulthood than children with hypertension in prepuberty did.65

Elevated blood pressure at a young age predicts HBP in adulthood, which is the leading risk factor for cardiovascular diseases. Studies have shown that HBP in childhood correlates with carotid intima-media thickness, atherosclerosis, left ventricular hypertrophy and kidney failure in adulthood.66–68 Childhood BP, overweight and obesity showed the largest OR for adult hypertension, followed by parental hypertension status in different longitudinal studies.4 23 69 Screening children and adolescents for elevated blood pressure could identify hypertension at an early stage, decreasing the rate of progression of hypertension from childhood to adulthood and reducing the risk of cardiovascular disease in adulthood.70

Strengths and limitations

To our knowledge, this is the first study investigating the relationships between ATP2B1 rs2681472 polymorphism and HBP among children and adolescents aged 12–15 years in the Baltic countries. Furthermore, we used multivariate logistic regression analyses to evaluate the associations under different inheritance models, separately for boys and girls. However, the present research has some limitations. The subject selection was limited to 12–15-year-old children and adolescents. Therefore, these results need to be confirmed in a large-scale study involving a larger population of younger children and older adolescents. In this study, BP readings were obtained by an automatic oscillometric BP monitor, although, according to the Fourth Report of the National High Blood Pressure Education Program (NHBPEP), HBP readings should be repeated by using auscultation.20 The present study did not evaluate pubertal status and biochemical parameters. In addition, there was no adjustment for family history of hypertension because information obtained from self-reports was lacking. The design of our study did not allow us to determine the cause–effect relationship.

Conclusion

In conclusion, CACNB2 rs12258967 gene polymorphism was significantly associated with higher odds of HBP in Lithuanian adolescents aged 12–15 years. ATP2B1 rs2681472 gene polymorphism was not associated with the odds of HBP. Further studies, particularly those evaluating the effect of genetic–anthropometric–environmental interactions on blood pressure, could help us to understand new pathophysiological mechanisms of the regulation of BP and to find potential targets for treatments.

Acknowledgments

The authors would like to thank Jurate Medzioniene for carrying out the statistical analysis. The authors also would like to thank the participants of the present study.

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.
  25. 25.
  26. 26.
  27. 27.
  28. 28.
  29. 29.
  30. 30.
  31. 31.
  32. 32.
  33. 33.
  34. 34.
  35. 35.
  36. 36.
  37. 37.
  38. 38.
  39. 39.
  40. 40.
  41. 41.
  42. 42.
  43. 43.
  44. 44.
  45. 45.
  46. 46.
  47. 47.
  48. 48.
  49. 49.
  50. 50.
  51. 51.
  52. 52.
  53. 53.
  54. 54.
  55. 55.
  56. 56.
  57. 57.
  58. 58.
  59. 59.
  60. 60.
  61. 61.
  62. 62.
  63. 63.
  64. 64.
  65. 65.
  66. 66.
  67. 67.
  68. 68.
  69. 69.
  70. 70.

Footnotes

  • Contributors SS, RK and VD contributed to the conception or design of the work. SS, RK and VL contributed to the acquisition, analysis or interpretation of data for the work. SS and RK drafted the manuscript. VL critically revised the manuscript. All authors gave final approval and agree to be accountable for all aspects of work ensuring integrity and accuracy.

  • Funding This research was funded by a grant (no. LIG-02/2011) from the Research Council of Lithuania.

  • Competing interests None declared.

  • Patient consent Parental/guardian consent obtained.

  • Ethics approval The study was approved by Kaunas Regional Biomedical Research Ethics Committee at the Lithuanian University of Health Sciences (protocol No. BE-2-69).

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

  • Data sharing statement No additional data are available.