BRIEF RESEARCH REPORT article
Genetic Variation of Kallikrein-Kinin System and Related Genes in Patients With Hereditary Angioedema
- 1Department of Biophysics, Federal University of São Paulo, São Paulo, Brazil
- 2Department of Dermatology and Allergy Centre, Odense University Hospital, Odense, Denmark
- 3Institute of Genetic Medicine, Johns Hopkins University, Baltimore, MD, United States
- 4Private Allergy and Immunology Clinic, Belo Horizonte, Brazil
- 5Department of Medicine, Allergy Clinic of Professor Polydoro Ernani de São Thiago University Hospital, Federal University of Santa Catarina, Florianopolis, Brazil
- 6Allergy Section, Internal Medicine Department, Hospital Universitari Vall d'Hebron, Barcelona, Spain
- 7Division of Clinical Immunology, Faculty of Medicine ABC, Santo André, Brazil
- 8Department of Internal Medicine, University of Campinas, Campinas, Brazil
Hereditary angioedema (HAE) is an autosomal dominant disease caused by C1-INH deficiency due to mutations in SERPING1 (C1-INH-HAE) in most of the cases, or by specific mutations in factor XII gene, F12 (F12-HAE). Identification of polymorphisms in the genes encoding proteins from key pathways driving HAE can help to understand how genetic diversity contributes to its phenotypic variability. Here, 15 genes related to the Kallikrein-Kinin System (KKS) were analyzed by next generation sequencing in 59 patients with C1-INH-HAE or F12-HAE from Brazil, Denmark and Spain, and 19 healthy relatives in a total of 31 families. We identified 211 variants, from which 23 occurred only in Danish subjects and 79 were found only in Brazilian individuals, resulting in 109/211 variations in common between European and Brazilian population in the HAE families analyzed. BDKRB2 and CPM presented a large number of variants in untranslated regions, 46/49 and 19/24, respectively; whereas ACE (n = 26), SERPING1 (n = 26), CPM (n = 24), and NOS3 (n = 16) genes presented the higher number of variants directly affecting amino acid sequence. Despite the large amount of variants identified, the lack of association between genotype and phenotype indicates that the modulation of HAE symptom requires a more complex regulation, probably involving pathways beyond the KKS, epigenetics and environmental factors. Considering the new HAE types recently described, molecules involved in the regulation of vasculature and in plasminogen activation become promising targets for future genetic studies.
Hereditary angioedema (HAE) is a rare autosomal dominant disease caused by deleterious mutations in SERPING1, leading to quantitative or functional C1-inhibitor deficiency (C1-INH-HAE) (1). Less frequently, specific mutations in F12 gene (F12-HAE) can be detected in patients with HAE with normal C1-inhibitor (2, 3). In addition, two new forms of HAE have been reported, HAE with mutations in angiopoietin-1 gene (ANGPT1) (4) and HAE with mutations in plasminogen gene (PLG) (5). A non-determined but substantial proportion of patients with HAE and normal C1-inhibitor cannot be explained by mutational findings at the moment and therefore is described as HAE of unknown cause (U-HAE) (1, 6). C1-INH-HAE and F12-HAE are related to the Kallikrein-Kinin System (KKS) activation and augmented production of bradykinin (BK), leading to vasodilation and angioedema episodes (1). In the KKS, factor XII (encoded by F12 gene) activates plasma kallikrein (KK; KLKB1 gene), which hydrolyzes high molecular weight kininogen (HMWK; KNG1 gene) releasing BK (7). BK binds the B2-receptor (BDKRB2 gene) leading to endothelial nitric oxide synthase activation (NOS3 gene). Carboxypeptidases N and M (genes CPN and CPM) metabolize BK into des-Arg9-BK, which preferentially activates B1-receptor (BDKRB1). Deleterious variations in genes responsible for the control of BK release, as SERPING1, or its metabolization, as CPN, CPM, angiotensin-converting enzyme (ACE), and neprilysin (MME), are able to increase the amount of BK released or its half-life, leading to or intensifying angioedema episodes. Whereas, deleterious variations in F12, KLKB1, and KNG1 could impair BK formation, working as a protective factor in BK-mediated angioedema (7). Although a small phenotype distinction can be noticed in HAE subtypes (2), there is a huge variation in clinical characteristics even among carriers of the same disease-causing mutation (1, 6).
A typical genome differs from the reference human genome at 4.1–5.0 million sites with >99.9% being single nucleotide polymorphisms and short indels (8). This large number of variants can possibly be disease-causing or may contribute to phenotypic variation and drug-response susceptibility. Indeed, a cumulative effect of multiple polymorphisms may lead to complex-trait genetic diseases (8). Here, we hypothesize that variants in the KKS could explain the phenotypic variation observed in HAE.
Materials and Methods
Seventy-eight subjects analyzed by targeted next-generation sequencing (NGS) belonged to 26 C1-INH-HAE families (45 mutation carriers, 15 healthy relatives) and 5 F12-HAE families (14 mutation carriers, 4 healthy relatives) (Table 1). Additional 56 HAE patients were included in the study and had single gene mutations sequenced by Sanger. Patients were previously diagnosed by biochemical C1-INH and C4 measurements and SERPING1 and/or F12 genotyping by Sanger sequencing.
This study was carried out in accordance with the recommendations of ethical standards of the 1964 Declaration of Helsinki with written informed consent from all subjects. The protocol was approved by the Ethics Committee in Research of Universidade Federal de São Paulo (n°56522).
Since only a few possible genetic variants have been suggested to be associated to the HAE phenotype (9–13), we explored 15 genes related to KKS in individuals from Brazil, Denmark, and Spain, by NGS as previously described (14). Briefly, the targeted genes (ACE, MME, KNG1, KLKB1, KLK1, F12, ENPEP, BDKRB1, BDKRB2, NOS3, TAC1, CPN1, CPM, SERPING1, PRCP) were amplified from patients' DNA extracted from blood samples by using Ion AmpliSeq™ Library Kit-2.0V (Life Technologies, Carlsbad, CA, USA). NGS was performed with the Ion 316™ Chip-v2 in Ion PGM™ Sequencer, which determined the amount of sequence data analyzed. To guarantee a high-quality coverage (minimum 20X) we limited the panel to the 15 selected genes, in a total of 366 amplicons. Torrent Suit v3.2.1 and Ion Reporter 4.0 (Life Technologies) were used to data analysis and interpretation of variants. Probably pathogenic variants, low coverage regions and SERPING1 and F12 coding regions were validated by Sanger sequencing. Variants annotation was performed using Ingenuity® Variant Analysis (Qiagen) and Metacore™ (Thomson Reuters).
In addition to the 78 subjects evaluated by NGS, in 56 HAE patients the exons 2 and 3 of BDKRB2 were also sequenced by Sanger method (11 C1-INH-HAE and 45 U-HAE). U-HAE patients were diagnosed according to Cicardi et al.'s guidelines (1). Distribution of variants between groups was calculated by Chi-square test.
Patients were grouped and compared according to nationality (Brazilian vs. European), HAE type (C1-INH-HAE vs. F12-HAE), edema localization (face, extremities, abdominal pain, upper respiratory tract, and genitals), mean frequency of attacks (per year) and duration of episodes (days). Patients were categorized in three groups according to the disease onset: (1) < 12-year-old, (2) from 12 to 25-year-old, and (3) >25-year-old. Mean frequency and duration of swelling episodes data were considered at baseline without prophylaxis or acute treatment. The frequency of attacks was divided in low (< 13 attacks per year), medium (from 13 to 24), and high (>24). Patients were also grouped according to the mean duration of swelling episodes in short duration (< 2 days), medium (2–3 days), and long duration (>3 days). Distribution of variants was compared within the selected groups.
In 26 C1-INH-HAE families, 23 mutations were considered responsible for HAE, and all F12-HAE patients carried p.Thr328Lys (Table 1). A total of 211 different variants were identified in the 15 genes analyzed. Nine alterations occurred in the 5′-UTR (4.3%) and 89 in the 3′-UTR (42.2%). Missense mutations corresponded to 22.7% of the total, followed by synonymous (21.8%), small insertions/deletions leading to frameshift (3.8%), non-sense (1.9%), splice site (1.4%), and intronic alterations (1.9%). A list containing all variants is found in Table S1 (Supplementary Material). BDKRB2 presented the highest number of variants (n = 49), followed by ACE (n = 26) and SERPING1 (n = 26), and CPM (n = 24). The distribution of variants per type or per gene was not significantly between C1-INH-HAE and F12-HAE groups (Chi-square test with Yates' correction). Considering only variants leading to change in the structure of the protein (missense, nonsense, indels, and predicted splicing sites), the genes in which more mutations were found were SERPING1 (n = 23) and ACE (n = 10). No probably damaging variants were identified in BDKRB1, CPN/M, KNG1, KLK1, ENPEP, PRCP, or TAC1 genes. Few studies suggest tissue kallikrein (KLK1 gene) may be involved in HAE manifestations (15). Tissue kallikrein releases Lys-BK from low molecular weight kininogen (also encoded by KNG1), which is converted into BK by aminopeptidases (16) or into Lys-des-Arg9-BK by carboxypeptidases.
The change p.Arg14Cys in BDKRB2 was found 15 individuals, and plGly354Glu in 16, without significant difference in distribution among HAE subtypes (p = 0.1203 and p = 0.1900, respectively, Chi-square test). The alterations p.Val376Met (rs141958164) and p.Val208Ile (rs139203012) were found in one U-HAE patient each.
Clinical Features and Phenotype-Genotype Association
We found 23 mutations exclusively in Danish individuals and 79 only in Brazilians, resulting in 109/211 mutations in common between European and Brazilian HAE families analyzed.
We compared the symptomatic patients regarding the presence of facial swellings (43/57), swellings of extremities (41/57), abdominal swellings (45/57), upper airways (28/57), and edema of the genitalia (21/57). The first attacks occurred before the age of 12 in 26 patients (24 C1-INH-INH, 2 F12-HAE) (mean = 6.7 years); between 12 and 25 years of age in 21 patients (20 C1-INH-HAE, 5 F12-HAE) (mean = 17), and after 25 years of age in eight (3 C1-INH-HAE, 5 F12-HAE) (mean = 33.8). Most of the patients reported to experience up to 12 swelling episodes/year without prophylaxis (31/47; mean = 4.6). Nine C1-INH-HAE patients reported 13–24 attacks/year (mean = 19.8), and seven reported more than 24 attacks/year (mean = 43.2).
Distribution of variants was compared within the selected groups, but no specific mutation was significantly associated to any of the clinical characteristic selected. However, five missense mutations were classified as probably damaging by all in silico predictive software (excluding SERPING1 variations) (Table 2).
HAE can be caused by C1-INH deficiency, in which more than 300 different mutations have already been considered responsible for the disease (17). A smaller proportion of HAE patients present specific mutations affecting a proline-rich region of factor XII that leads to loss of glycosylation and enhanced autoactivation mechanism of the zymogen activating the KKS (18). In spite of the large mutational spectrum of SERPING1 in HAE, the variation in the presentation of symptoms seems not to be related to the type of the disease-causing mutation or even with functional levels of C1-INH (19). Therefore, the search for new genetic biomarkers could not only explain the variation in HAE symptomatology, but also be used in diagnosis, prognosis, and treatment response. Proteins of the KKS are strong candidates to be regulating the severity and frequency of HAE attacks, since they are directly involved in the BK-release pathway.
In this study, a large number of variants was found in a regulatory portion of the 3′-untranslated BDKRB2 region (20). Also, 40/78 individuals presented the 9 bp deletion in the exon 1 of BDKRB2, described to increase B2-receptor expression (21). However, no correlation was found with these variants and HAE clinical characteristics suggesting that the regulation of B2-receptor expression through these regions is not strong enough to modulate HAE attacks itself, in agreement with previous data (22).
Variants classified as probably pathogenic in ACE (p.Gly354Arg, p.Tyr244Cys, p.Thr916Met) were found in Brazilian and Danish symptomatic patients, and p.Arg324Trp was found only in one non-affected relative (family 21). Since ACE degrades BK, deleterious alterations could increase BK half-life, favoring longer edema episodes. It has been demonstrated that ACE activity inversely correlates with severity scores in F12-HAE (23); in contrast, C1-INH-HAE severity is not depending on ACE, but highly depends on aminopeptidase P activity (24). However, the C1-INH-HAE patients carrying these mutations had distinguished clinical features from each other (Table 2).
The mutation p.Cys548Tyr in KLKB1, found in two patients from family 1 and one asymptomatic mutation carrier from family 27 (F12-HAE), has been associated with KK deficiency in a compound heterozygosity with the change p.Trp402Ter (25). KK deficiency is a rare autosomal recessive condition which causes an augmented activated partial thromboplastin time, but no coagulation anomalies. Therefore, coexistence of KK deficiency could attenuate or even abolish the symptoms of HAE, since it is the key enzyme responsible for BK release. The only missense mutation found together with p.Cys548Tyr was the high frequent polymorphism p.Ser143Asn (66/78). This polymorphism was shown to cause KK deficiency with the change p.Gly125Arg, both in homozygosity (26) and lying in the apple domain 2 of the heavy chain of KK, through which kallikrein interacts with HMWK domain 6 (27). However, p.Ser143Asn was found in heterozygosity in the three patients carrying the p.Cys548Tyr, discarding the hypothesis that this association could be responsible by itself for the lack of symptoms observed in the F12-HAE carrier. Furthermore, the Cys548 lies in the catalytic domain of KK, far from the apple domain 2.
Pathogenicity level of most of the new variants cannot be confirmed only by in silico analysis (28). For example, the variant p.Thr328Lys in F12 is responsible for most of F12-HAE cases (3, 18), but is classified as benign by many pathogenicity predictors. Therefore, the discovery of new and rare variants possibly involved in HAE pathogenesis requires not only experimental evidences to assure the pathogenicity degree of each variant but also the analysis of the overlap effect of different variations and genes. Thus, despite the large number of variants identified in KKS and related genes in this work, the small number of subjects included and the lack of biochemical and functional analysis on the variants are some limitations.
The studies of modulating factors in HAE must combine the discovery of new mutations, genotype-phenotype association including a large number of patients, and consider epigenetic and environmental factors. To our knowledge, this is the largest study exploring different genes in many HAE families. We expect that the disclosure of the variants here presented in HAE patients and the new recent findings in U-HAE will encourage further studies correlating the KKS and the regulation of vasculature, as well as larger international multicentric studies.
CV coordinated the analysis, and drafted the manuscript. AA and AB designed data collection. CV, PN, and AM performed all the experiments. CV, RM, and RF-S performed the bioinformatics analysis. AA, RG, JD, MG, AG, EM, and AB performed the data collection instruments, and coordinated data collection. AB conceptualized the study and supervised data collection. JP conceptualized and designed the study, coordinated and supervised all the experiments and analysis and revised the manuscript. All the authors drafted, reviewed and approved the manuscript.
This study received grants from Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP) (grant numbers 2013/02661-4, 2014/27198-8 and 2015/25494-1).
Conflict of Interest Statement
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.
We are thankful to the patients and their families for being part of this study, and to Dr. Claus Koch from the Department of Cancer & Inflammation Research, Institute of Molecular Medicine, University of Southern Denmark, Odense, Denmark.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmed.2019.00028/full#supplementary-material
Table S1. Genetic variants found by NGS in HAE families. All the 211 variants found in the 15 genes analyzed in the 78 subjects are listed in the table. Pathogenicity prediction from SIFT (D, damaging; T, tolerated), PolyPhen-2 (B, benign; P, possibly damaging; D, damaging), PROVEAN (D, damaging; N, non-damaging) and CADD (threshold > 20) indicate the in silico analysis results. †Found only in Danish individuals.
Table S2. Full dataset of genetic variants and clinical data. The sheet entitled Genetic Variants contains all the variants found by NGS before validation and patients filtering (raw data). Variants with the annotation ARTIFACT in the column Functional Class correspond to the variations not confirmed by Sanger sequencing validation and discarded in the final evaluation. The sheet Patients Raw contains the chip and barcode identification, and clinical data from all the patients initially sequenced, including U-HAE patients posteriorly excluded from NGS analysis.
ACE, angiotensin-converting enzyme; C1-INH-HAE, hereditary angioedema with C1 inhibitor deficiency; F12-HAE, hereditary angioedema with mutation in F12 gene; HMWK, high molecular weight kininogen; KK, Plasma Kallikrein; KKS, Kallikrein-Kinin System; U-HAE, hereditary angioedema of unknown cause.
1. Cicardi M, Aberer W, Banerji A, Bas M, Bernstein JA, Bork K, et al. Classification, diagnosis, and approach to treatment for angioedema: consensus report from the hereditary angioedema international working group. Allergy (2014) 69:602–16. doi: 10.1111/all.12380.
3. Veronez CL, Moreno AS, Constantino-Silva RN, Maia LSM, Ferriani MPL, Castro FFM, et al. Hereditary Angioedema with Normal C1 Inhibitor and F12 mutations in 42 Brazilian families. J Allergy Clin Immunol Pract. (2018) 6:1209–16.e8. doi: 10.1016/j.jaip.2017.09.025
4. Bafunno V, Firinu D, D'Apolito M, Cordisco G, Loffredo S, Leccese A, et al. Mutation of the angiopoietin-1 gene (ANGPT1) associates with a new type of hereditary angioedema. J Allergy Clin Immunol. (2018) 141:1009–17. doi: 10.1016/j.jaci.2017.05.020
6. Bouillet L, Boccon-Gibod I, Launay D, Gompel A, Kanny G, Fabien V, et al. Hereditary angioedema with normal C1 inhibitor in a French cohort: clinical characteristics and response to treatment with icatibant. Immun Inflamm Dis. (2017) 5:29–36. doi: 10.1002/iid3.137
7. Hofman ZL, Relan A, Zeerleder S, Drouet C, Zuraw B, Hack CE. Angioedema attacks in patients with hereditary angioedema: Local manifestations of a systemic activation process. J Allergy Clin Immunol. (2016) 138:359–66. doi: 10.1016/j.jaci.2016.02.041
9. Speletas M, Szilágyi Á, Csuka D, Koutsostathis N, Psarros F, Moldovan D, et al. F12-46C/T polymorphism as modifier of the clinical phenotype of hereditary angioedema. Allergy (2015) 70:1661–4. doi: 10.1111/all.12714
10. Piñero-Saavedra M, González-Quevedo T, Saenz de San Pedro B, Alcaraz C, Bobadilla-González P, Fernández-Vieira L, et al. Hereditary angioedema with F12 mutation: clinical features and enzyme polymorphisms in 9 Southwestern Spanish families. Ann Allergy Asthma Immunol. (2016) 117:520–6. doi: 10.1016/j.anai.2016.09.001
11. Gianni P, Loules G, Zamanakou M, Kompoti M, Csuka D, Psarros F, et al. Genetic Determinants of C1 inhibitor deficiency angioedema age of onset. Int Arch Allergy Immunol. (2017) 174:200–4. doi: 10.1159/000481987
12. Cumming SA, Halsall DJ, Ewan PW, Lomas DA. The effect of sequence variations within the coding region of the C1 inhibitor gene on disease expression and protein function in families with hereditary angio-oedema. J Med Genet. (2003) 40:e114. doi: 10.1136/jmg.40.10.e114
13. Duponchel C, Di Rocco C, Cicardi M, Tosi M. Rapid detection by fluorescent multiplex PCR of exon deletions and duplications in the C1 inhibitor gene of hereditary angioedema patients. Hum Mutat. (2001) 17:61–70. doi: 10.1002/1098-1004(2001)17:1 < 61::AID-HUMU7>3.0.CO;2-9
14. Veronez CL, da Silva ED, Lima Teixeira PV, Cagini N, Constantino-Silva RN, Grumach AS, et al. Genetic analysis of hereditary angioedema in a Brazilian family by targeted next generation sequencing. Biol Chem. (2016) 397:315–22. doi: 10.1515/hsz-2015-0212
15. Charest-Morin X, Hébert J, Rivard GÉ, Bonnefoy A, Wagner E, Marceau F. Comparing pathways of bradykinin formation in whole blood from healthy volunteers and patients with hereditary angioedema Due to C1 inhibitor deficiency. Front Immunol. (2018) 9:2183. doi: 10.3389/fimmu.2018.02183
16. Goto Y, Hattori A, Ishii Y, Mizutani S, Tsujimoto M. Enzymatic properties of human aminopeptidase a. Regulation of its enzymatic activity by calcium and angiotensin IV. J Biol Chem. (2006) 281:23503–13. doi: 10.1074/jbc.M603191200
18. Björkqvist J, de Maat S, Lewandrowski U, Di Gennaro A, Oschatz C, Schönig K, et al. Defective glycosylation of coagulation factor XII underlies hereditary angioedema type III. J Clin Invest. (2015) 125:3132–46. doi: 10.1172/JCI77139
20. Zamorano R, Suchindran S, Gainer JV. 3'-Untranslated region of the type 2 bradykinin receptor is a potent regulator of gene expression. Am J Physiol Renal Physiol. (2006) 290:F456–64. doi: 10.1152/ajprenal.00009.2005
21. Lung CC, Chan EK, Zuraw BL. Analysis of an exon 1 polymorphism of the B2 bradykinin receptor gene and its transcript in normal subjects and patients with C1 inhibitor deficiency. J Allergy Clin Immunol. (1997) 99:134–46.
22. Freiberger T, Vyskocilová M, Kolárová L, Kuklínek P, Krystufková O, Lahodná M, et al. Exon 1 polymorphism of the B2BKR gene does not influence the clinical status of patients with hereditary angioedema. Hum Immunol. (2002) 63:492–4. doi: 10.1016/S0198-8859(02)00397-X
23. Charignon D, Ghannam A, Defendi F, Ponard D, Monnier N, López Trascasa M, et al. Hereditary angioedema with F12 mutation: factors modifying the clinical phenotype. Allergy (2014) 69:1659–65. doi: 10.1111/all.12515
24. Drouet C, Désormeaux A, Robillard J, Ponard D, Bouillet L, Martin L, et al. Metallopeptidase activities in hereditary angioedema: effect of androgen prophylaxis on plasma aminopeptidase P. J Allergy Clin Immunol. (2008) 121:429–33. doi: 10.1016/j.jaci.2007.10.048
25. Lombardi AM, Sartori MT, Cabrio L, Fadin M, Zanon E, Girolami A. Severe prekallikrein (Fletcher factor) deficiency due to a compound heterozygosis (383Trp stop codon and Cys529Tyr). Thromb Haemost. (2003) 90:1040–5. doi: 10.1160/TH03–05-0275
26. Katsuda I, Maruyama F, Ezaki K, Sawamura T, Ichihara Y. A new type of plasma prekallikrein deficiency associated with homozygosity for Gly104Arg and Asn124Ser in apple domain 2 of the heavy-chain region. Eur J Haematol. (2007) 79:59–68. doi: 10.1111/j.1600-0609.2007.00871.x
27. Renné T, Sugiyama A, Gailani D, Jahnen-Dechent W, Walter U, Müller-Esterl W. Fine mapping of the H-kininogen binding site in plasma prekallikrein apple domain 2. Int Immunopharmacol. (2002) 2:1867–73. doi: 10.1016/S1567-5769(02)00170-4
28. Richards S, Aziz N, Bale S, Bick D, Das S, Gastier-Foster J, et al. Standards and guidelines for the interpretation of sequence variants: a joint consensus recommendation of the american college of medical genetics and genomics and the association for molecular pathology. Genet Med. (2015) 17:405–24. doi: 10.1038/gim.2015.30
Keywords: hereditary angioedema, Kallikrein-Kinin System, genetic variation, genotype-phenotype correlation, C1 inhibitor deficiency, F12 mutation
Citation: Veronez CL, Aabom A, Martin RP, Filippelli-Silva R, Gonçalves RF, Nicolicht P, Mendes AR, Da Silva J, Guilarte M, Grumach AS, Mansour E, Bygum A and Pesquero JB (2019) Genetic Variation of Kallikrein-Kinin System and Related Genes in Patients With Hereditary Angioedema. Front. Med. 6:28. doi: 10.3389/fmed.2019.00028
Received: 01 December 2018; Accepted: 30 January 2019;
Published: 21 February 2019.
Edited by:Giancarlo Castaman, Università degli Studi di Firenze, Italy
Reviewed by:Coen Maas, University Medical Center Utrecht, Netherlands
Christian Drouet, Université Grenoble Alpes, France
Copyright © 2019 Veronez, Aabom, Martin, Filippelli-Silva, Gonçalves, Nicolicht, Mendes, Da Silva, Guilarte, Grumach, Mansour, Bygum and Pesquero. 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: João Bosco Pesquero, email@example.com