- Research article
- Open Access
- Open Peer Review
Epidemiology and molecular characterization of the re-emerging measles virus among children and adults in the Haut-Ogooue, Gabon
BMC Infectious Diseasesvolume 19, Article number: 90 (2019)
Measles is one of the most infectious diseases with a high mortality rate worldwide. It is caused by the measles virus (MeV) which is a single stranded RNA virus with genetic diversity based on the nucleoprotein gene, including 24 genotypes. In Gabon, several outbreaks occurred in the past few years, especially in 2016 in Libreville and Oyem. A surveillance network of infectious diseases highlighted a measles outbreak which occurred in the south of Gabon from April to June 2017.
Clinical specimens of urine, blood, throat and nasal swabs were collected in the two main cities of the Haut-Ogooue province, Franceville and Moanda. Virological investigations based on real-time polymerase chain reaction for molecular diagnosis and conventional PCR for genotype identification were done.
Specimens were collected from 139 suspected measles patients. A total of 46 (33.1%) children and adults were laboratory-confirmed cases among which 16 (34.8%) were unvaccinated, 16 (34.8%) had received one dose, and 11 (23.9%) had received two doses of the measles vaccine. Phylogenetic analysis revealed that all the sequences of the nucleoprotein gene belonged to genotype B3.
Measles infection was more commonly confirmed among those with one recorded dose compared to suspect cases with none, unknown or two recorded doses. The molecular characterization of the strains showed the circulation of the B3 genotype which is endemic on the African continent, thirty years after the B2 genotype was described in an outbreak in Libreville, the capital of Gabon. These findings highlight that surveillance and molecular investigation of measles should be continued in Gabon.
Measles is a highly contagious infectious disease caused by the measles virus (MeV) which is a negative-sense single stranded RNA virus, belonging to the Morbillivirus genus, Measles morbillivirus species, and the Paramyxoviridae family. After an incubation period of 10–14 days, measles is characterized by respiratory infection symptoms such as fever (38 °C), cough, coryza, and conjunctivitis, followed by a maculopapular rash and high fever (39–40 °C). Measles complications which include pneumonia, diarrhea, vomiting and encephalitis are more common among children under five years old and adults . The measles virus’ genetic phylogeny is divided into 8 clades (A-H) and subdivided into 24 genotypes based on the diversity of the C-terminal hypervariable domain (a 456 nucleotide sequence) of the nucleoprotein gene .
Measles is one of the most infectious diseases with a high mortality rate worldwide. Therefore, the World Health Organization (WHO) recommended a molecular surveillance of measles strains and established a vaccination program for measles elimination [3, 4]. In 2010, the World Health Assembly planned to increase routine coverage with the first dose of measles vaccine for children aged one year old to ≥90% nationally and ≥ 80% in every district, reduce and maintain annual measles incidence at < 5 cases per million and reduce mortality by 95% from a 2000 estimate in order to reach eradication by 2015 . From 2000 to 2015, the number of measles cases reported annually worldwide decreased by 70%, from 853,479 to 254,928 cases, and measles incidence decreased by 75%, from 146 to 36 cases per million population . From 2000 to 2016, measles vaccination prevented an estimated 20.4 million deaths. Global measles deaths have decreased by 84% from an estimated 550,100 in 2000 to 89,780 in 2016 . Accelerated immunization activities have had a major impact on reducing measles deaths, and yet from 2016 to 2017, more than 66,000 laboratory-confirmed measles cases have been reported from the different regions of the world . Hence, several reports concluded that measles continued to spread across many countries because the vaccination coverage was sub-optimal .
In Gabon, several outbreaks of measles were reported, especially in 1984, 1992 and 2001, of which 7234, 1245 and 5129 reported cases occurred, respectively. More recently, 1509 cases were notified in 2016 . In 2012, WHO set up supplementary measles immunization activities nationwide in Gabon in children aged 9 months − 5 years old in order to protect children who did not respond to the first measles vaccination. This program reached 168,769 (67%) children in the targeted age group . A study performed on clinical specimens collected during the 1984 measles outbreak in Libreville, the capital of Gabon, showed that the strains investigated belonged to genotype B2 . No data has been available in the past years on the molecular characterization of measles strains in Gabon. In this study we describe the epidemiology and molecular investigation of an outbreak of measles and report the genotype strain which circulated in the Haut-Ogooue province, Gabon, in 2017.
Patients and samples
A measles outbreak occurred from April to June 2017 in the south-east of Gabon. A surveillance network of infectious diseases was set up according to the guidelines of the Gabonese Ministry of Health. Suspected measles cases were patients who visited health centers in the two main cites (Franceville and Moanda) of the Haut-Ogooue’s province in Gabon, with clinical symptoms such as fever (≥ 38 °C), maculopapular rash, cough, coryza and conjunctivitis. Clinical specimens of urine, blood, throat and nasal swabs were collected and sent to the Centre International de Recherches Médicales de Franceville (CIRMF). Demographical data such as the patient’s name, age, sex and neighborhood of residence were recorded. Furthermore, clinical data (fever, cough, coryza, rash, conjunctivitis, diarrhea, vomiting), the date of the onset of the rash and the vaccination status were also reported. All subjects gave their informed consent for inclusion before they participated in the study. These suspected cases were confirmed in the laboratory of the CIRMF with virological and molecular investigations.
Nucleic acid extraction and amplification
RNA was extracted, from serum, urine, nasal and throat swabs, on a BioRobot EZ1 workstation (Qiagen) using the EZ1 Virus Mini Kit version 2.0 (Qiagen) according to the manufacturer’s instructions. Fast Real-time PCR was performed using the primers and the probe for the 188 bp virus nucleoprotein (NP) target . This probe was labeled at the 5′ ends with the FAM quencher and at the 3’ends with the Black Hole Quencher 1. After a first step of reverse transcription, each 15 μl of reaction mixture contained 5 μl of eluted DNA, 10 μl of Master Mix Fast (Applied Biosystems), 0.25 μM of primers and probe. A 7500 Fast Real Time PCR system was run for 2 min at 50 °C, for 20 s at 95 °C, followed by 45 cycles at 95 °C for 15 s and 60 °C for 30 s.
Measles virus genotyping and phylogenetic analysis
The measles virus positive samples were genotyped by sequencing the 450 nucleotides coding for the carboxyl terminus of the nucleocapsid protein. The NP gene fragment was amplified using a conventional ReverseTranscriptase-Polymerase Chain Reaction (RT-PCR) method (SuperScript III Reverse Transcriptase kit; Invitrogen) with the forward primer (5’-CGGGCAAGAGATGGTAAGGAGGTCAG -3′) and the reverse primer (5′- AGGGTAGGCGGATGTTGTTCTGG -3′) . Each 25-μl reaction mixture contained 1X buffer, 0.1 mM of magnesium sulfate, 40 ng/μl of bovine serum albumin, 0.4 μM of primers, 1 μl of Platinium Taq (Invitrogen) and 5 μl of eluted RNA. Cycling conditions were reverse transcription for 30 min at 50 °C, denaturation for 3 min at 95 °C, followed by 45 cycles of amplification at 95 °C for 30 s, 62 °C for 45 s, 72 °C for 1 min, and a final extension for 10 min at 72 °C. The generated amplicons were sequenced by the Sanger method with a 3500 Genetic Analyzer (Applied Biosystems).
Phylogenetic analyses were performed using a multiple sequence alignment of 10 MeV sequences obtained from patients of our study with a selection of reference strains divided into 6 genotypes among the 24 available in the GenBank database. The sequences were compared with measles strain sequences on the BLAST (Basic Local Alignment Search Tool) and MeaNS (Measles Nucleotide Surveillance) databases for obtaining the genotype of the strains. Sequences were aligned using ClustalX (version 1.81). Phylogenetic relationships were reconstructed using the best-fitting ML model based on the Akaike information criterion, GTR (general time reversible) + ModelTest. The phylogenetic tree was built using the maximum-likelihood method with the PhyML algorithm [13,14,15,16,17,18] and drawn using FigTree v.1.4.0.
Statistical analysis was performed using Statview V5.0 software. Data were analyzed with the Pearson’s chi-squared and Fisher’s exact tests. Statistical tests were used to compare measles symptoms and the vaccination status among the laboratory-confirmed measles cases. A two-tailed critical alpha value of 0.05 was used. A p-value below 0.05 was considered to indicate statistical significance.
Sequence accession number
The sequences of the nucleoprotein N gene of MeV are available in the DDBJ/EMBL/GenBank database under the accession number MK393355-MK393364.
Demographic, clinical and vaccination parameters
Among the 139 patients with suspected measles, 77 were males (55.4%), 62 were females (44.6%) giving a sex ratio (M:F) of 1.24 (Table 1). The age ranged from 2 months to 51 years. The median age was 2.0 years and the mean age was 5.1 ± 8.0 years. The age group distribution was as follows: 17 patients (12.2%) were under nine months old, 39 patients (28.1%) were 9–17 months old, 51 patients (36.7%) were aged 18 months-5 years old, 21 (15.1%) patients were 6–14 years old and 11 (7.9%) were 15–51 years old (Table 1).
Among the suspected cases, 34 (24.5%) were unvaccinated patients whereas 67 (48.2%) and 34 (24.4%) had received one and two doses of vaccine, respectively. The vaccination status of 4 (2.9%) patients was unknown (Table 1).
Epidemiology of the measles virus outbreak
The measles virus was detected in 46/139 (33.1%) patients. The weekly distribution of the measles virus outbreak revealed that the number of suspected cases and the incidence of the measles virus peaked at week S16 and gradually decreased until week S25 (Fig. 1).
The distribution of laboratory-confirmed measles cases among age groups was as follows: 6 (13%) in the < 9 months-old age group, 4 (8.7%) in the 9–17-months age group, 18 (39.1%) in the 18 months-5 years age group, 11 (23.9%) in the 6–14-years age group and 7 (15.2%) in the ≥15 years old age group (Table 1). Sixty per cent (28/46) of laboratory-confirmed measles cases were children ≤5 years old, among which 14 (50%) were vaccinated (Table 2).
Among the 46 patients infected with the measles virus, 27 (58.7%) were vaccinated, 16 (34.8%) were unvaccinated and no data was available on the vaccination status of 3 (6.5%) patients (Table 2).
The rate of patients with the measles virus who received one and two doses of vaccine was 34.8% (16/46) and 23.9% (11/46), respectively. Among the < 9 months and 9–17 months age groups, only two children had received one dose of vaccine and another had received two doses. Children in the 9–17 months age group had the lowest positive rate of measles (10%) (Table 1). Confirmation of measles infection was highest amongst suspected cases over 15 years of age at 63.6% (Table 1). In this group, among the 7 positive cases, two patients were unvaccinated, 4 had received one dose of vaccine and 1 had an unknown vaccination status (Table 2).
The most common clinical symptoms were cough (97.8%), coryza (78.3%) and rash (69.6) (Table 2). There was a significant difference in fever between patients vaccinated with one dose, patients vaccinated with two doses of vaccine and unvaccinated patients (X2 = 11.5, p = 0.0095) (Table 2). Statistical analyses showed that 81% of unvaccinated laboratory-confirmed cases had a high fever ≥39 °C. Half of the patients who had received one dose of vaccine had high fever. However, 72% of patients who had received two doses of vaccine had fever < 39 °C. The two patients with conjunctivitis were positive measles cases, but their vaccination status was unknown. Likewise, the percentage of patients with diarrhea was highest among unvaccinated laboratory-confirmed cases and lowest among patients who had received two doses (X2 = 7.7, p = 0.0518) (Table 2).
The phylogenetic analysis based on 456 nucleotides sequence of the nucleoprotein N gene was performed using 10 sequences of this measles outbreak and sequences of the GenBank database. According to the WHO nomenclature used to describe the genotypes of measles virus strains , our sequences were named as follows: MVs (extracted from clinical material), city and country where the case occurred, number of the week and years the clinical symptoms appeared, followed by the isolate number (Table 3, Fig. 2).
We obtained 36 sequences from patients infected with the measles virus among which 22 were similar to the sequence from isolate 3, 2 were similar to the sequence from isolate 23, 1 was similar to the sequence from isolate 17 and another was similar to the sequence from isolate 71 (Fig. 2). All the 8 sequences obtained from the specimens from Moanda’s health center, were similar to the sequence from isolate 3 from Franceville.
The 10 obtained sequences were compared with WHO reference sequences of measles on database of MeaNS which are recommended for genotype identification. Genotyping based on the C-terminal 450 nucleotides of the N gene of measles virus revealed that the 10 strains displayed 99 and 98% identity at the nucleotide level with the strains MVi/Ibadan.NGA/0.97/ and MVi/New_York.USA/0.94/, respectively, belonging to the genotype B3. Morever, the obtained strains displayed 99% identity with a strain isolated in France in 2016 (Accession no. KX838946) and clustered with genotype B3.
This outbreak of the measles virus which occurred April to June 2017 in the Haut-Ogooue province, in Gabon demonstrates the challenges with controlling measles worldwide. Measles has been re-emerging in several countries in Europe and Asia in the past two years, particularly in France since 2017 [19,20,21]. The report of the French health authorities in April 2018 reported 1527 cases of measles since the beginning of the year. In Central Africa, Cameroon reported the endemic circulation of the B3 MeV genotype between 2010 and 2011 and during a recent measles outbreak in 2015 [22, 23]. In Gabon, this study is the first to describe the circulation of genotype B3. Data on measles outbreaks showed circulation of genotype B2 in 1984 in the capital Libreville. Furthermore, several suspected patients were reported in Libreville and Oyem (Northern Gabon) in 2016 . However there was no data available on the molecular characterization of these strains. In this study, more than half of the positive cases were children under five years old, suggesting that vaccination coverage is not sufficient. The vaccination program in Gabon provides one dose of vaccine between 9 and 12 months old, and the second dose should be administered at 18 months. The 6 children under 9 months old accounted for 13% of positive measles cases. This finding suggests that these young children did not have maternal antibodies and they were susceptible to measles infection until they received the first dose of vaccine. Among these 6 children, 4 were aged 5–6 months. Considering that maternal immunity should last 4–6 months, it would seem that the mothers of these 4 children might also be unprotected and susceptible to measles. A study mentioned that the high rate among young infants suggested low prevalence of maternal measles antibodies due to a large number of measles-susceptible mothers . More than half of the children aged 6 to 14 years old and the adults infected with the measles virus had received only one dose of vaccine. Our results are compatible with studies describing outbreaks showing that measles is likely endemic until a vaccination coverage of at least 95% is reached in children and adults [25, 26]. Measles was most frequently confirmed in suspected cases aged 15 years and over. Our results are compatible with those of Namibian and Chinese studies showing that the highest age incidence occurred among young adults and infants [27, 28].
There were no significant differences in measles symptoms between vaccinated and unvaccinated laboratory-confirmed cases for cough, coryza, rash and vomiting. However, there was an inversely proportional relationship between high fever and vaccination status. High fever (39–40 °C) was most frequent in unvaccinated patients. Likewise, the percentage of patients with diarrhea was highest among unvaccinated cases and lowest among patients who had received two doses. A less severe symptom profile in these latter cases could probably be due to a secondary failure of the vaccination. The severity of the disease seemed to depend on the vaccination status and the immunity of the patient. One review mentioned that mild form of measles could be observed in subjects with higher immunity to the virus .
The fact that the number of laboratory–confirmed cases continue to spread despite vaccination could probably be due to the fact that vaccination coverage was sub-optimal, as is the case in many countries of the European Union . Failure in vaccine protection may not be due to the circulating genotype but to the number of vaccine doses received. Indeed, children who received only one dose of vaccine develop measles due to vaccine failures (as expected in 10% of cases). A second dose of vaccine to all children would reduce the expected vaccine failure rate to around 1%. The measles virus only has one serotype. As a result, the anti-measles immune response protects against the disease regardless of the genotype .
Among the 36 sequences from patients infected with the measles virus, 10 strains displayed one or two nucleotides difference. These 10 strains that have been compared with the WHO measles surveillance network (MeaNS) formed a cluster belonging to genotype B3. According to the World Health Organization, the global distribution of the measles genotype from 2010 to 2015 revealed that the B3 genotype was spread worldwide. However, this genotype is endemic and the most frequent reported on the African continent . Molecular characterization of MeV strain from Nigeria and Ghana suggested that clade B was prevalent in Sub-Saharan Africa . In Gabon, data showed that the B2 genotype was the only once circulating in 1984 . The B2 genotype seemed to disappear in the 2000s and was declared inactive by WHO in 2003 . The insufficient molecular surveillance of measles in Central Africa could explain these results . However, several studies reported evidence of the continued circulation of B2 on the African continent, in South Africa (2002), Angola (2003), and the Democratic Republic of the Congo (2004) [12, 32]. Furthermore, a co-circulation of the B2 and B3.1 genotypes of the measles virus in the Central African Republic was reported in the 2000 measles epidemic . In the last 10 years, most studies reported the circulation of genotype B3, particularly B3.1, in Central Africa, especially in Cameroon between 2010 and 2011, and then in 2014. Our finding was consistent with the results which reported the circulation of the B3 genotype in Central Africa (Fig. 1).
The syndromic surveillance of infectious diseases allowed us to describe a measles outbreak which occurred from April to June 2017 in the Haut- Ogooue province, in Gabon. Measles was commonly confirmed in patients who had received one dose of vaccine, due to the expected failure rate of around 10%. Molecular characterization of the strains showed the circulation of the B3 genotype, the same strain circulating in several neighboring countries.
These data enhanced our knowledge on the measles virus strains circulating in Gabon and suggest that surveillance and molecular investigation of measles strains should be pursued to better understand the burden of measles outbreaks. Following this epidemic, the Gabonese Ministry of Health immediately set up a vaccination program.
World Health Organization
Freymuth F, Dina J, Parent du Chatelet I, Mourez B, Waku-Kouomou D, Vabret A. La rougeole et son virus. Virologie. 2011;15(1):6–22.
WHO. Measles virus nomenclature update: 2012. Wkly Epidemiol Rec. 2012;87(9):73–80.
Perry RT, Gacic-Dobo M, Dabbagh A, Mulders MN, Strebel PM, Okwo-Bele JM, et al. Global control and regional elimination of measles, 2000-2012. MMWR Morb Mortal Wkly Rep. 2014;63(5):103–7.
Kalaycioglu AT, Yolbakan S, Guldemir D, Korukluoglu G, Coskun A, Cosgun Y, et al. Towards measles elimination: phylogenetic analysis of measles viruses in Turkey (2012-2013) and identification of genotype D8. J Med Virol. 2016;88(11):1867–73.
WHO. Progress towards regional measles elimination – worldwide, 2000–2015. Wkly Epidemiol Rec. 2016;91(45):525–36.
WHO. Measles. 2018.
WHO. WHO measles and rubella surveillance data, 24 jul 2017. 2017.
Grobusch MP, Rodriguez-Morales AJ, Wilson ME. Measles on the move. Travel Med Infect Dis. 2017;18:1–2.
WHO. WHO vaccine-preventable diseases: monitoring system. 2017 global summary. In: Incidence time series for Gabon (GAB); 2018.
Waku-Kouomou D, Alla A, Blanquier B, Jeantet D, Caidi H, Rguig A, et al. Genotyping measles virus by real-time amplification refractory mutation system PCR represents a rapid approach for measles outbreak investigations. J Clin Microbiol. 2006;44(2):487–94.
Ammour Y, Faizuloev E, Borisova T, Nikonova A, Dmitriev G, Lobodanov S, et al. Quantification of measles, mumps and rubella viruses using real-time quantitative TaqMan-based RT-PCR assay. J Virol Methods. 2013;187(1):57–64.
Kremer JR, Nkwembe E, Bola Oyefolu AO, Smit SB, Pukuta E, Omilabu SA, et al. Measles virus strain diversity, Nigeria and Democratic Republic of the Congo. Emerg Infect Dis. 2010;16(11):1724–30.
Anisimova M, Gascuel O. Approximate likelihood-ratio test for branches: a fast, accurate, and powerful alternative. Syst Biol. 2006;55(4):539–52.
Castresana J. Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol Biol Evol. 2000;17(4):540–52.
Chevenet F, Brun C, Banuls AL, Jacq B, Christen R. TreeDyn: towards dynamic graphics and annotations for analyses of trees. BMC Bioinformatics. 2006;7:439.
Dereeper A, Guignon V, Blanc G, Audic S, Buffet S, Chevenet F, et al. Phylogeny.fr: robust phylogenetic analysis for the non-specialist. Nucleic Acids Res. 2008;36(Web Server):W465–9.
Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32(5):1792–7.
Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol. 2010;59(3):307–21.
Caseris M, Burdet C, Lepeule R, Houhou N, Yeni P, Yazdanpanah Y, et al. An update on measles. Rev Med Interne. 2015;36(5):339–45.
Currie J, Davies L, McCarthy J, Perry M, Moore C, Cottrell S, et al. Measles outbreak linked to European B3 outbreaks, Wales, United Kingdom 2017. Euro Surveill. 2017;22(42).
Kaida A, Iritani N, Kanbayashi D, Yamamoto SP, Hirai Y, Hakui N, et al. Ten-year surveillance of Measles virus from 2007-2016 in Osaka City, Japan. Jpn J Infect Dis. 2018;71(2):152–4.
Demanou M, Ratsitoharana R, Yonga M, Dosseh A, Anya B, Kobela M, et al. Molecular characterization of measles viruses that circulated in Cameroon between 2010 and 2011. Virol J. 2013;10:71.
Njim T, Agyingi K, Aminde LN, Atunji EF. Trend in mortality from a recent measles outbreak in Cameroon: a retrospective analysis of 223 measles cases in the Benakuma Health District. Pan Afr Med J. 2016;23:135.
Hagan JE, Takashima Y, Sarankhuu A, Dashpagma O, Jantsansengee B, Pastore R, et al. Risk factors for Measles virus infection among adults during a large outbreak in Postelimination era in Mongolia, 2015. J Infect Dis. 2017;216(10):1187–95.
Freymuth F, Vabret A. Measles, a re-emerging disease in France? Clin Microbiol Infect. 2011;17(6):793.
Gay NJ. The theory of measles elimination: implications for the design of elimination strategies. J Infect Dis. 2004;189(Suppl 1):S27–35.
Ma C, Yan S, Su Q, Hao L, Tang S, An Z, et al. Measles transmission among adults with spread to children during an outbreak: implications for measles elimination in China, 2014. Vaccine. 2016;34(51):6539–44.
Ogbuanu IU, Muroua C, Allies M, Chitala K, Gerber S, Shilunga P, et al. Measles outbreak reveals measles susceptibility among adults in Namibia, 2009 - 2011. S Afr Med J. 2016;106(7):715–20.
Stokes J Jr, Reilly CM, Buynak EB, Hilleman MR. Immunologic studies of measles. Am J Hyg. 1961;74:293–303.
WHO. Policy statement on data sharing by WHO in the context of public health emergencies (as of 13 April 2016). Wkly Epidemiol Rec. 2016;91(18):237–40.
Hanses F, Truong AT, Ammerlaan W, Ikusika O, Adu F, Oyefolu AO, et al. Molecular epidemiology of Nigerian and Ghanaian measles virus isolates reveals a genotype circulating widely in western and Central Africa. J Gen Virol. 1999;80(Pt 4):871–7.
Smit SB, Hardie D, Tiemessen CT. Measles virus genotype B2 is not inactive: evidence of continued circulation in Africa. J Med Virol. 2005;77(4):550–7.
Gouandjika-Vasilache I, Waku-Kouomou D, Menard D, Beyrand C, Guye F, Ngoay-Kossy JC, et al. Cocirculation of measles virus genotype B2 and B3.1 in Central African Republic during the 2000 measles epidemic. J Med Virol. 2006;78(7):964–70.
We acknowledge Heïdi Lançon for the English revision of the manuscript. We thank all the hospital staff of the health centers. We thank Alain-Prince Okouga for his technical assistance, Lekolo Gwladys Lekolo Sangoye a medicine intern and Kassa Kassa Roland Fabrice a doctor of the CIRMF.
Availability of data and materials
The sequences of the nucleoprotein N gene of MeV are available in the DDBJ/EMBL/GenBank database (https://www.ncbi.nlm.nih.gov/pubmed) under accession numbers MK393355-MK393364.
Ethics approval and consent to participate
This study was approved by the Gabonese Ministry of Health which granted its ethics approval (number 000292/VP/MSPSSN/SG/DGS). It was set up in response to the measles outbreak which occurred in Gabon. We obtained informed verbal consent from patients and their parents or guardians (for minors) for inclusion before they participated in the study, which was approved by the Ethics Committee.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.