Sander Pajusalu
Mari-Anne Vals
Laura Mihkla
Ustina Šamarina
Tiina Kahre1,2
Katrin Õunap- 1Department of Clinical Genetics, Institute of Clinical Medicine, University of Tartu, Tartu, Estonia
- 2Department of Clinical Genetics, United Laboratories, Tartu University Hospital, Tartu, Estonia
- 3Children’s Clinic, Tartu University Hospital, Tartu, Estonia
Congenital disorders of glycosylation (CDG) are a widely acknowledged group of metabolic diseases. PMM2-CDG is the most frequently diagnosed CDG with a prevalence as high as one in 20,000. In contrast, the prevalence of other CDG types remains unknown. This study aimed to analyze the estimated prevalence of different N-linked protein glycosylation disorders. We extracted allele frequencies for diverse populations from The Genome Aggregation Database (gnomAD), encompassing variant frequency information from 141,456 individuals. To identify pathogenic variants, we used the ClinVar database as a primary source. High confidence loss-of-function variants as defined by the LOFTEE algorithm were also classified as pathogenic. After summing up population frequencies for pathogenic alleles, estimated disease birth prevalence values with confidence intervals were calculated using the Bayesian method. We first validated our approach using two more common recessive disorders (cystic fibrosis and phenylketonuria) by showing that the estimated prevalences calculated from population allele frequencies were in accordance with previously published epidemiological studies. Among assessed 27 autosomal recessive N-glycosylation disorders, the only disease with estimated birth prevalence higher than one in 100,000 was PMM2-CDG (in both, all gnomAD individuals and those with European ancestry). The combined prevalence of 27 different N-glycosylation disorders was around one in 22,000 Europeans but varied considerably across populations. We will show estimated prevalence data from diverse populations and explain the possible pitfalls of this analysis. Still, we are confident that these data will guide CDG research and clinical care to identify CDG across populations.
Introduction
Congenital disorders of glycosylation (CDG) are a growing group of metabolic diseases with at least 137 defects (Ondruskova et al., 2021). According to the International Classification of Inherited Metabolic Disorders (ICIMD) database, 32 are classified as N-linked protein glycosylation defects (Ferreira et al., 2021). Based on population allele frequencies, the expected birth prevalence of the most common PMM2-CDG could be as high as 1:20,000 (Schollen et al., 2000), and in later reports, 1:77,000 to 1:286,000 (Vals et al., 2018; Yildiz et al., 2020). Mainly at the expense of PMM2-CDG, type 1 defects are more frequently diagnosed than type 2 defects (Peanne et al., 2017; Medrano et al., 2019) but compared to PMM2-CDG, other defects occur much seldom. The estimated combined prevalence of CDG among the Saudi population is 14 per million (Alsubhi et al., 2017), whereas, in Poland, the prevalence is approximately one case per million (Lipinski et al., 2021). The worldwide individual prevalence of different defects remains unknown.
The emergence of large-scale population-based exome and genome sequencing studies and the aggregation into even more extensive cross-population databases like gnomAD makes allele frequency data available from more than 100,000 individuals across different populations (Karczewski et al., 2020). Disease prevalence can be estimated when these data are combined with clinical variant classification databases like ClinVar (Landrum et al., 2018) and variant loss-of-function predictions. This approach has been previously used on assessing disease prevalence of different limb-girdle muscular dystrophy subtypes across populations (Liu et al., 2019). Although direct calculations using the Hardy-Weinberg equation can be used, the methods developed by Liu et al. (2019) make use of more advanced Bayesian statistics to add confidence intervals on prevalence estimators.
This study aimed to analyze the estimated prevalence of N-linked protein glycosylation defects across different populations based on allele frequencies in general population databases.
Materials and Methods
Only N-linked protein glycosylation defects with autosomal recessive inheritance patterns (27/32) were included in the analysis. The list was based on the ICIMD database (Ferreira et al., 2021). We extracted allele frequencies for different populations from The Genome Aggregation Database (gnomAD v2.1.1) (Karczewski et al., 2020), encompassing variant frequency information from 141,456 individuals. We defined variants as pathogenic if they were classified pathogenic or likely pathogenic in the ClinVar database (version 20210119) (Landrum et al., 2018) or were marked as high confidence loss-of-function (Moller et al., 2016) variants by the LOFTEE algorithm (Karczewski et al., 2020). In the case of conflicting annotations, high confidence LoF variants marked as (likely) benign in ClinVar were considered as non-pathogenic. Also, variants not passing all quality filters in the gnomAD variant call set were eliminated from counting toward the sum of pathogenic alleles.
The calculations were carried out across all seven gnomAD populations (African/African American, Latino/Admixed American, non-Finnish European, Finnish, Ashkenazi Jewish, East Asian, and South Asian). The number of individuals per population ranged from 5,185 (Ashkenazi Jewish) to 64,603 (non-Finnish European). In addition, we calculated prevalence estimates for the Estonian population (2,418 individuals), which is a subpopulation of non-Finnish Europeans. After obtaining population frequencies for pathogenic alleles, estimated disease birth prevalence values with confidence intervals were calculated using Bayesian statistics adapted from the previously published study by Liu et al. (2019). No manual variant curation was performed, and no gene-based exceptions were made to assure reproducibility and robustness (i.e., for PMM2-CDG, the p.Arg141His variant homozygotes were not excluded although known to be embryonically lethal). For validation, we used two common recessive disorders (cystic fibrosis and phenylketonuria) and compared the estimated disease prevalence with the data from epidemiological studies conducted in different populations.
Results
The estimated prevalence from population allele frequencies of cystic fibrosis and phenylketonuria was in accordance with previously published data (Supplementary Table 1). Results from published epidemiological studies fell into 95% confidence intervals in non-Finnish European, Estonian, and Finnish populations, where reliable population prevalence is known.
The selected list of N-linked protein glycosylation defects, their estimated prevalence, and the number of reported cases is presented in Table 1. The full list of all defects in all assessed populations is shown in Supplementary Table 2.
Table 1. N-linked protein glycosylation defects with highest estimated prevalence.
PMM2-CDG showed the highest estimated prevalence counted from 71 different pathogenic variants, and it is more prevalent in European (1 in 27,000), Ashkenazi Jewish (1 in 20,000), and admixed American populations (1 in 64,000). All other N-linked protein glycosylation defects had a much lower prevalence.
FUK-CDG and MAN2B2-CDG, which are both only recently reported and classified as N-linked protein glycosylation defects, showed high estimated prevalence, especially in East-Asians where FUK-CDG prevalence was estimated at 1:12,000 and MAN2B2-CDG at 1:11,000.
The estimated prevalence of combined N-linked protein glycosylation defects in Europeans was one in 22,000 or one in 24,000 after excluding FUK-CDG and MAN2B2-CDG. In East-Asians, however, the total CDG prevalence was estimated at 1:5,613, dropping to 1:121,935 after excluding FUK-CDG and MAN2B2-CDG. CDG seems to be more common in Finnish (1:16,000) and Ashkenazi Jewish populations (1:14,000). In Estonians, the combined CDG birth prevalence is estimated at around 1:50,000.
Discussion
This study presents the estimated prevalence of different N-linked protein glycosylation defects calculated from population allele frequencies. As the CDG group involves at least 137 defects, we emphasize that we only included 27 autosomal recessive protein N-glycosylation affecting defects and excluded defects that affect multiple glycosylation pathways.
The used methods make several assumptions and do not take variant- and gene-specific properties into account. For example, the calculations assume that pathogenic variants are independent and not on the same allele. Also, Hardy-Weinberg equilibrium is assumed but may not be present in some populations with more consanguinity. All diseases are considered to result from the biallelic loss-of-function mechanism; however, we cannot exclude the possibility that for some CDG, only specific variants (e.g., certain missense variants) cause the disease.
Also, this method did not consider the possibility that some variant recombinants might not be compatible with life. Since we cannot exclude this phenomenon in these genes, the chance of p.Arg141His homozygosity was also included when calculating the estimated prevalence of PMM2-CDG. However, as seen from the subanalysis of the Estonian population, it does not change the prevalence that much. In Estonia, the estimated prevalence of PMM2-CDG with the exclusion of p.Arg141His homozygosity was 1:77,000, and with inclusion, 1:62,000 (Vals et al., 2018). This example shows that the estimated prevalence of different defects might be somewhat overrated.
On the other hand, some aspects can cause an underestimation of disease prevalence. For example, currently, the data relies on gnomAD variants, and thus ultra-rare pathogenic variants not seen in gnomAD are not counted toward prevalence estimation. Also, only small variants (single nucleotide substitutions and small indels) were included, and thus pathogenic copy-number variants (deletions) may raise the disease prevalence estimators. Regarding missense variants, only those classified as pathogenic in ClinVar were included, while some other missense, as well as synonymous, intronic, and regulatory variants, may truly be pathogenic but missed due to lack of such information in databases. Regarding missense variants, one possibility to improve the estimates would be to include variants with high pathogenicity predictions. However, as shown by Liu et al. (2019) and consistent with our experience (data not shown), this leads to an overestimation of prevalence. Thus, we did not include any missense variant prediction tools in our final calculations. To be on the conservative side, we eliminated all low-quality variants in gnomAD; however, some of them, especially indels, may be true and pathogenic. Of note, presented results rely on data extracted from two publicly available databases, ClinVar and gnomAD, at a single time point. If other databases would have been used, we expect to see slightly different results; however, confidence intervals are shown to address this variability. Moreover, as the databases update both allele frequencies and pathogenicity classifications, a follow-up study in the future is planned to investigate whether the prevalence estimates will change after some years.
Despite many pitfalls, the validation with two well-known recessive diseases—cystic fibrosis and phenylketonuria—proved the method’s accuracy and usability for estimating the disease prevalence. From our previous epidemiological studies in Estonia, the prevalence of phenylketonuria is known to be 1:6,700 (Lillevali et al., 2018) correlating with the estimate from gnomAD data is 1:6,604 (95% CI 1:4,269-1:12,073). For cystic fibrosis, the birth prevalence in Estonia is 1:7,743 (Kahre, 2004) comparing well with the calculated estimate 1:8,482 (95% CI 1:5,351-1:16,269). For other well-studied populations like Europeans and Finnish, the data is in good correlation as well (Supplementary Table 1).
As shown in Table 1, compared to PMM2-CDG, all other N-linked protein glycosylation defects are less common. Still, our reported prevalence of different defects shows compliance with published clinical observations. In 2016, an informal inquiry was conducted in various European laboratories to evaluate the number of screening positive and molecularly confirmed CDG-I and CDG-II patients (Peanne et al., 2017). Based on their results, the most common type 1 N-glycosylation defect in Europe was PMM2-CDG, followed by ALG6-CDG, ALG1-CDG, and MPI-CDG. If we compare this distribution with non-Finnish Europeans, a similar pattern is seen with type 1 N-glycosylation defects. The most common type 2 N-glycosylation defects, according to the inquiry, were MAN1B1-CDG, MGAT2-CDG, and B4GALT1-CDG. In the non-Finnish European population, MAN1B1-CDG is followed by MOGS-CDG, which does not show an abnormal serum transferrin profile with routine CDG screening (De Praeter et al., 2000), and therefore was not reported by laboratories. Compared to these two defects, MGAT2-CDG and B4GALT1-CDG show a much lower prevalence.
Unexpectedly, two recently described defects, FUK-CDG and MAN2B2-CDG, showed the prevalence of higher than one per million. To date, FUK-CDG is described only in two individuals, whereas MAN2B2-CDG only in one (Ng et al., 2018; Verheijen et al., 2020). For FUK, we found 119, and for MAN2B2, 84 different possibly disease-causing alleles matching our definition of pathogenicity. Disagreement between relatively high estimated prevalence and the low number of reported cases could be explained by the phenomenon where homozygosity or compound heterozygosity of some variants may not be viable, as is the case with p.Arg141His homozygosity in PMM2 (Matthijs et al., 1998). Also, FUK and MAN2B2 are not included in smaller panels and are only identified by whole-exome sequencing. Theoretically, it may be the cause for the underdiagnosis of these defects to some extent. Notably, one of the pathogenic missense variants reported by Ng et al. (2018), NM_145059.3:c.2980A > C (p.Lys994Gln), reaches allele frequency of 0.002 in the East Asian population, and there is one homozygous individual present in gnomAD. Although mild phenotypes cannot be excluded in gnomAD individuals, also non-pathogenicity of this allele, at least in the homozygous state, should be considered. Of note, FUK and MAN2B2 both lack homozygous LoF variant carriers in gnomAD, which hints that biallelic LoF may be an actual pathogenic mechanism for those genes.
Unlike FUK-CDG and MAN2B2-CDG, other recently described defects (GFUS-CDG, OSTC-CDG, and FUT8-CDG) show lower prevalence.
Although the individual prevalence of each N-linked protein glycosylation defect is very low, the combined prevalence for this group of CDG is notable. If all 27 defects are included, the combined prevalence in non-Finnish Europeans is one in 22,000. If we exclude FUK-CDG and MAN2B2-CDG, the prevalence in Europeans is slightly lower, one in 24,000. Compared to available data (Vals et al., 2018; Lipinski et al., 2021), the expected and observed prevalence differ. The lower observed prevalence might be somewhat underestimated because of undiagnosed patients due to different causes such as unrecognized phenotypes, negative screening results or limited access to different metabolic and molecular studies. In the Estonian population, the combined prevalence of 25 defects is 1:50,000. Therefore, with a yearly birth rate of 13,000–14,000, one individual with protein N-glycosylation defect should be born every 4 years correlating with our clinical data.
In summary, we broaden the knowledge about N-linked protein glycosylation disorders, which hopefully helps raise awareness about the prevalence of this CDG subgroup.
Data Availability Statement
Publicly available datasets were analyzed in this study. These data can be found at: https://gnomad.broadinstitute.org/ and https://www.ncbi.nlm.nih.gov/clinvar/.
Author Contributions
SP planned and conducted the study, carried out the statistical analysis, and compiled the manuscript. M-AV planned, conducted the study, analyzed the data, and compiled the manuscript. LM conducted literature review, analyzed the data, and reviewed the manuscript. UŠ and TK analyzed the data and reviewed the manuscript. KÕ planned and conducted the study and reviewed the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was supported and funded by the Estonian Research Council grants PRG471 and MOBTP175.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2021.719437/full#supplementary-material
References
Abiramalatha, T., Arunachal, G., Muthusamy, K., and Thomas, N. A. (2019). Family with floppy neonates with severe respiratory insufficiency: A lethal phenotype of RFT1-CDG due to a novel mutation. Eur. J. Med. Genet. 62, 248–253. doi: 10.1016/j.ejmg.2018.07.023
Alsubhi, S., Alhashem, A., Faqeih, E., Alfadhel, M., Alfaifi, A., Altuwaijri, W., et al. (2017). Congenital disorders of glycosylation: The Saudi experience. Am. J. Med. Genet. Part A 173, 2614–2621.
Altassan, R., Peanne, R., Jaeken, J., Barone, R., Bidet, M., Borgel, D., et al. (2019). International clinical guidelines for the management of phosphomannomutase 2-congenital disorders of glycosylation: Diagnosis, treatment and follow up. J. Inherit. Metab. Dis. 42, 5–28.
Anzai, R., Tsuji, M., Yamashita, S., Wada, Y., Okamoto, N., Saitsu, H., et al. (2021). Congenital disorders of glycosylation type IIb with MOGS mutations cause early infantile epileptic encephalopathy, dysmorphic features, and hepatic dysfunction. Brain Dev. 43, 402–410. doi: 10.1016/j.braindev.2020.10.013
Balasubramanian, M., Johnson, D. S., and Study, D. D. D. (2019). MAN1B-CDG: Novel variants with a distinct phenotype and review of literature. Eur. J. Med. Genet. 62, 109–114. doi: 10.1016/j.ejmg.2018.06.011
Bastaki, F., Bizzari, S., Hamici, S., Nair, P., Mohamed, M., Saif, F., et al. (2018). Single-center experience of N-linked congenital disorders of glycosylation with a summary of molecularly characterized cases in arabs. Ann. Hum. Genet. 82, 35–47. doi: 10.1111/ahg.12220
Cechova, A., Altassan, R., Borgel, D., Bruneel, A., Correia, J., Girard, M., et al. (2020). Consensus guideline for the diagnosis and management of mannose phosphate isomerase-congenital disorder of glycosylation. J. Inherit. Metab. Dis. 43, 671–693. doi: 10.1002/jimd.12241
De la Morena Barrio, M. E., Sabater, M., de la Morena Barrio, B., Ruhaak, R. L., Minano, A., Padilla, J., et al. (2020). ALG12-CDG: An unusual patient without intellectual disability and facial dysmorphism, and with a novel variant. Mol. Genet. Genomic Med. 8:e1304.
De Praeter, C. M., Gerwig, G. J., Bause, E., Nuytinck, L. K., Vliegenthart, J. F., Breuer, W., et al. (2000). A novel disorder caused by defective biosynthesis of N-linked oligosaccharides due to glucosidase I deficiency. Am. J. Hum. Genet. 66, 1744–1756. doi: 10.1086/302948
Ferreira, C. R., Rahman, S., Keller, M., Zschocke, J., and Group, I. A. (2021). An international classification of inherited metabolic disorders (ICIMD). J. Inherit. Metab. Dis. 44, 164–177.
Hock, M., Wegleiter, K., Ralser, E., Kiechl-Kohlendorfer, U., Scholl-Burgi, S., Fauth, C., et al. (2015). ALG8-CDG: novel patients and review of the literature. Orphanet J. Rare Dis. 10:73.
Karczewski, K. J., Francioli, L. C., Tiao, G., Cummings, B. B., Alfoldi, J., Wang, Q., et al. (2020). The mutational constraint spectrum quantified from variation in 141,456 humans. Nature 581, 434–443.
Landrum, M. J., Lee, J. M., Benson, M., Brown, G. R., Chao, C., Chitipiralla, S., et al. (2018). ClinVar: improving access to variant interpretations and supporting evidence. Nucleic Acids Res. 46, D1062–D1067.
Lekka, D. E., Brucknerova, J., Salingova, A., Sebova, C., Ostrozlikova, M., Ziburova, J., et al. (2021). Congenital disorders of glycosylation - an umbrella term for rapidly expanding group of rare genetic metabolic disorders - importance of physical investigation. Bratisl. Lek. Listy. 122, 190–195. doi: 10.4149/bll_2021_030
Lillevali, H., Reinson, K., Muru, K., Simenson, K., Murumets, U., Mols, T., et al. (2018). Hyperphenylalaninaemias in estonia: genotype-phenotype correlation and comparative overview of the patient cohort before and after nation-wide neonatal screening. JIMD Rep. 40, 39–45. doi: 10.1007/8904_2017_61
Lipinski, P., Bogdanska, A., and Tylki-Szymanska, A. (2021). Congenital disorders of glycosylation: Prevalence, incidence and mutational spectrum in the Polish population. Mol. Genet. Metab. Rep. 27:100726. doi: 10.1016/j.ymgmr.2021.100726
Liu, W., Pajusalu, S., Lake, N. J., Zhou, G., Ioannidis, N., Mittal, P., et al. (2019). Estimating prevalence for limb-girdle muscular dystrophy based on public sequencing databases. Genet. Med. 21, 2512–2520. doi: 10.1038/s41436-019-0544-8
Lo Barco, T., Osanni, E., Bordugo, A., Rodella, G., Iascone, M., Tenconi, R., et al. (2021). Epilepsy and movement disorders in CDG: Report on the oldest-known MOGS-CDG patient. Am. J. Med. Genet. Part A 185, 219–222. doi: 10.1002/ajmg.a.61916
Matthijs, G., Schollen, E., Van Schaftingen, E., Cassiman, J. J., and Jaeken, J. (1998). Lack of homozygotes for the most frequent disease allele in carbohydrate-deficient glycoprotein syndrome type 1A. Am. J. Hum. Genet. 62, 542–550. doi: 10.1086/301763
Medrano, C., Vega, A., Navarrete, R., Ecay, M. J., Calvo, R., Pascual, S. I., et al. (2019). Clinical and molecular diagnosis of non-phosphomannomutase 2 N-linked congenital disorders of glycosylation in Spain. Clin. Genet. 95, 615–626. doi: 10.1111/cge.13508
Moller, R. S., Larsen, L. H., Johannesen, K. M., Talvik, I., Talvik, T., Vaher, U., et al. (2016). Gene panel testing in epileptic encephalopathies and familial epilepsies. Mol. Syndromol. 7, 210–219. doi: 10.1159/000448369
Ng, B. G., Rosenfeld, J. A., Emrick, L., Jain, M., Burrage, L. C., Lee, B., et al. (2018). Pathogenic variants in fucokinase cause a congenital disorder of glycosylation. Am. J. Hum. Genet. 103, 1030–1037. doi: 10.1016/j.ajhg.2018.10.021
Ng, B. G., Shiryaev, S. A., Rymen, D., Eklund, E. A., Raymond, K., Kircher, M., et al. (2016). ALG1-CDG: clinical and molecular characterization of 39 unreported patients. Hum. Mutat. 37, 653–660.
Ondruskova, N., Cechova, A., Hansikova, H., Honzik, T., and Jaeken, J. (2021). Congenital disorders of glycosylation: Still “hot” in 2020. Biochim. Biophys. Acta Gen. Subj. 1865:129751. doi: 10.1016/j.bbagen.2020.129751
Peanne, R., de Lonlay, P., Foulquier, F., Kornak, U., Lefeber, D. J., Morava, E., et al. (2017). Congenital disorders of glycosylation (CDG): Quo vadis? Eur. J. Med. Genet. 61, 643–663. doi: 10.1016/j.ejmg.2017.10.012
Quelhas, D., Jaeken, J., Fortuna, A., Azevedo, L., Bandeira, A., Matthijs, G., et al. (2019). RFT1-CDG: absence of epilepsy and deafness in two patients with novel pathogenic variants. JIMD Rep. 43, 111–116. doi: 10.1007/8904_2018_112
Schollen, E., Kjaergaard, S., Legius, E., Schwartz, M., and Matthijs, G. (2000). Lack of Hardy-Weinberg equilibrium for the most prevalent PMM2 mutation in CDG-Ia (congenital disorders of glycosylation type Ia). Eur. J. Hum. Genet. 8, 367–371. doi: 10.1038/sj.ejhg.5200470
Sturiale, L., Bianca, S., Garozzo, D., Terracciano, A., Agolini, E., Messina, A., et al. (2019). ALG12-CDG: novel glycophenotype insights endorse the molecular defect. Glycoconj. J. 36, 461–472. doi: 10.1007/s10719-019-09890-2
Tahata, S., Gunderson, L., Lanpher, B., and Morava, E. (2019). Complex phenotypes in ALG12-congenital disorder of glycosylation (ALG12-CDG): Case series and review of the literature. Mol. Genet. Metab. 128, 409–414. doi: 10.1016/j.ymgme.2019.08.007
Vals, M. A., Pajusalu, S., Kals, M., Magi, R., and Ounap, K. (2018). The Prevalence of PMM2-CDG in estonia based on population carrier frequencies and diagnosed patients. JIMD Rep. 39, 13–17. doi: 10.1007/8904_2017_41
Verheijen, J., Wong, S. Y., Rowe, J. H., Raymond, K., Stoddard, J., Delmonte, O. M., et al. (2020). Defining a new immune deficiency syndrome: MAN2B2-CDG. J. Allergy Clin. Immunol. 145, 1008–1011. doi: 10.1016/j.jaci.2019.11.016
Vuillaumier-Barrot, S., Schiff, M., Mattioli, F., Schaefer, E., Dupont, A., Dancourt, J., et al. (2019). Wide clinical spectrum in ALG8-CDG: clues from molecular findings suggest an explanation for a milder phenotype in the first-described patient. Pediatr. Res. 85, 384–389. doi: 10.1038/s41390-018-0231-5
Keywords: CDG, congenital disorders of glycosylation, N-glycosylation, prevalence, population allele frequencies
Citation: Pajusalu S, Vals M-A, Mihkla L, Šamarina U, Kahre T and Õunap K (2021) The Estimated Prevalence of N-Linked Congenital Disorders of Glycosylation Across Various Populations Based on Allele Frequencies in General Population Databases. Front. Genet. 12:719437. doi: 10.3389/fgene.2021.719437
Received: 02 June 2021; Accepted: 21 July 2021;
Published: 10 August 2021.
Edited by:
Karolina Stepien, Salford Royal NHS Foundation Trust, United KingdomReviewed by:
Gregory M. Pastores, Mater Misericordiae University Hospital, IrelandJames O’Byrne, Mater Misericordiae University Hospital, Ireland
Copyright © 2021 Pajusalu, Vals, Mihkla, Šamarina, Kahre and Õunap. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Katrin Õunap, a2F0cmluLm91bmFwQGtsaWluaWt1bS5lZQ==
†These authors have contributed equally to this work and share first authorship