Prevalence of filaggrin gene polymorphisms (exon-3) in patients with atopic dermatitis in a multiracial Brazilian population

topic dermatitis (AD) is a chronic, multifactorial disease hose clinical phenotype results from the interaction of enetic and environmental factors.1 Immune dysregulation nd skin barrier integrity determine its severity, predisosing to infections and antigen permeability.2 It is a ommon reason for dermatological consultation, especially n childhood (<12 years) representing 25.8% of dermatologcal consultations.3 The gene encoding filaggrin (FLG) is highly polymorphic nd is located in the region of the epidermal differentiation omplex (1q21.3), encoding the most important proteins nvolved in skin barrier homeostasis. FLG is the main genetic actor associated with AD and its exon-3 transcribes most of he profilaggrin protein. Skin barrier alterations are present n AD patients without filaggrin changes; however, the presnce of FLG variants leading to loss of function has been ssociated with clinical phenotypes such as persistent earlynset disease, asthma, and allergic sensitization.1,4 Marked thnic disparity has been observed in the frequency of FLG ariants leading to loss of filaggrin function in AD.5 More than 60 variants in FLG leading to loss of filaggrin unction have been identified in association with AD, with he most common among Europeans being c.1537C>T:R501X nd 2282del4:S761Cfs*36 and, in sub-Saharan Africans, .9947C>G:S3316*. Few studies have been carried out on atin-American patients with AD. The objective of the resent study was to assess the frequency of FLG variants in exon-3) in AD patients to compare Brazilian and internaional populations and explore their clinical characteristics. A cross-sectional study was carried out at the Dermaology Outpatient Clinic (FMABC; Santo André, São Paulo). ighty patients with AD (Hanifin and Rajka’s criteria) of both exes were included and examined by an experienced deratologist to assess disease severity (SCORAD, EASI) and The specimen collected for genetic analysis was obtained by swabbing the patients’ mucosa of the inner cheek and placing it in a sterile test tube (Oragene Collector OG-500®, DNA Genotek Inc., Kanata, Ontario). It was submitted to sequencing using the Sanger method. DNA extraction was performed using ethanol precipitation and a prepIT 2P reagent provided together with the Oragene kit. Polymerase chain reaction (PCR) and its sequencing analysis were performed focusing on exon-3 to identify the most common genetic variants-c.1537C>T:R501X (rs61816761) and c.2282del4:S761Cfs*36 (rs558269137) using validated primers from Thermo Fischer Scientific® (Applied Biosystems, Foster City, CA), (Hs00274028 forward: 5’CTA ACA CTG GAT CCC TGG TTC CTA 3’ and reverse 5’ CTG AGA CAG CAG AGC CAC CAA GA 3’ and Hs00395823, forward: 5’ CAG ACC TAT CTA CCG ATT GCT CGT 3’ and reverse: 5’ AAA TCA GGC ACTCGT CAC ACA CAG AA 3’). This strategy allowed the investigation of other variants in areas neighboring the target loci, but it did not cover the entire exon-3 coding region (Fig. 1). The PCR products were purified with DNA beads (Agencourt -AMPure XP-Beckman Coulter, Brea, CA). The purified samples, together with 10 L of these primers, were used for the sequencing reaction. The sequencing cycle was performed with the Big Dye Terminator v3.1 kit (Thermo Fisher Scientific). The sequencing products were submitted to capillary electrophoresis on the ABI 3500 DNA Analyzer (Applied Biosystems, Foster City, CA). The sequencing data were evaluated with Seq A (14) software (Applied Biosystems, Foster City, CA). PROVEAN (Protein Variation Effect Analyzer) v1.1 was used to predict whether a protein sequence variation caused by a missense substitution would affect protein function (available at https://provean.jcvi.org/index.php). The prevalence of each identified variant was compared to the Brazilian public database of genomic variants (ABraOm; hg38 -https://abraom.ib.usp.br/) with 1,171 samples from the population of the same region and an international database (gnomAD; v3.1.2 -https://gnomad.broadinstitute.org/) with 76,156 unrelated o collect venous blood samples for laboratory analysis and


Dear Editor,
Atopic dermatitis (AD) is a chronic, multifactorial disease whose clinical phenotype results from the interaction of genetic and environmental factors. 1 Immune dysregulation and skin barrier integrity determine its severity, predisposing to infections and antigen permeability. 2 It is a common reason for dermatological consultation, especially in childhood (<12 years) representing 25.8% of dermatological consultations. 3 The gene encoding filaggrin (FLG) is highly polymorphic and is located in the region of the epidermal differentiation complex (1q21.3), encoding the most important proteins involved in skin barrier homeostasis. FLG is the main genetic factor associated with AD and its exon-3 transcribes most of the profilaggrin protein. Skin barrier alterations are present in AD patients without filaggrin changes; however, the presence of FLG variants leading to loss of function has been associated with clinical phenotypes such as persistent earlyonset disease, asthma, and allergic sensitization. 1,4 Marked ethnic disparity has been observed in the frequency of FLG variants leading to loss of filaggrin function in AD. 5 More than 60 variants in FLG leading to loss of filaggrin function have been identified in association with AD, with the most common among Europeans being c.1537C>T:R501X and 2282del4:S761Cfs*36 and, in sub-Saharan Africans, c.9947C>G:S3316*. Few studies have been carried out on Latin-American patients with AD. The objective of the present study was to assess the frequency of FLG variants (in exon-3) in AD patients to compare Brazilian and international populations and explore their clinical characteristics.
A cross-sectional study was carried out at the Dermatology Outpatient Clinic (FMABC; Santo André, São Paulo). Eighty patients with AD (Hanifin and Rajka's criteria) of both sexes were included and examined by an experienced dermatologist to assess disease severity (SCORAD, EASI) and to collect venous blood samples for laboratory analysis and oral mucosa for genetic analysis. The participants/guardians signed a free and informed consent form.
The specimen collected for genetic analysis was obtained by swabbing the patients' mucosa of the inner cheek and placing it in a sterile test tube (Oragene Collector OG-500 ® , DNA Genotek Inc., Kanata, Ontario). It was submitted to sequencing using the Sanger method.
DNA extraction was performed using ethanol precipitation and a prepIT 2P reagent provided together with the Oragene kit.
Polymerase chain reaction (PCR) and its sequencing analysis were performed focusing on exon-3 to identify the most common genetic variants---c.1537C>T:R501X (rs61816761) and c.2282del4:S761Cfs*36 (rs558269137) using validated primers from Thermo Fischer Scientific ® (Applied Biosystems, Foster City, CA), (Hs00274028 forward: 5'CTA ACA CTG GAT CCC TGG TTC CTA 3' and reverse 5' CTG AGA CAG CAG AGC CAC CAA GA 3' and Hs00395823, forward: 5' CAG ACC TAT CTA CCG ATT GCT CGT 3' and reverse: 5' AAA TCA GGC ACTCGT CAC ACA CAG AA 3'). This strategy allowed the investigation of other variants in areas neighboring the target loci, but it did not cover the entire exon-3 coding region (Fig. 1). The PCR products were purified with DNA beads (Agencourt ---AMPure XP-Beckman Coulter, Brea, CA). The purified samples, together with 10 L of these primers, were used for the sequencing reaction. The sequencing cycle was performed with the Big Dye Terminator v3.1 kit (Thermo Fisher Scientific). The sequencing products were submitted to capillary electrophoresis on the ABI 3500 DNA Analyzer (Applied Biosystems, Foster City, CA). The sequencing data were evaluated with Seq A (14) software (Applied Biosystems, Foster City, CA). PROVEAN (Protein Variation Effect Analyzer) v1.1 was used to predict whether a protein sequence variation caused by a missense substitution would affect protein function (available at https://provean.jcvi.org/index.php).
The prevalence of each identified variant was compared to the Brazilian public database of genomic variants (ABraOm; hg38 ---https://abraom.ib.usp.br/) with 1,171 samples from the population of the same region and an international database (gnomAD; v3.1.2 --https://gnomad.broadinstitute.org/) with 76,156 unrelated individuals of different ethnicities. The statistical significance was set at p ≤ 0.001. 6 The main demographic data, clinical severity, eosinophilia, and IgE levels are shown in Table 1. Most patients had moderate and severe AD (79%), elevated IgE levels (98%), and eosinophilia (68%). The results of the genetic analysis of the AD sample and the comparison with the two population controls are shown in Table 2. Twenty-six FLG exon-3 genetic variants were detected in 60 patients (75%; 95% CI 65%---85%). Homozygosity and compound heterozygosity were not identified.
Eight synonymous variants were found in the FLG gene, four of them widely distributed in Brazil and worldwide. However, four of them (c.1665T>G:rs152285733, c.1737A>C: synonymous, c.1521G>A: synonymous, and c.1800T>C: synonymous) were observed to be more common in the AD sample than in the world (p < 0.001). The variant c.1476T>C: synonymous, rare in the Brazilian and in the international population (p < 0.001), is common among African-Americans. The synonymous c.2544T>C, considered common among AD patients and regional controls, is extremely rare in the world (p < 0.001).
Fourteen missense variants were detected, with two of them (c.1712A>G:H559R, c.1777A>T:T581S) being more common among AD patients than in regional and international controls (p < 0.001).
No FLG variant was associated with AD clinical severity, eosinophilia, or elevated serum IgE levels (p > 0.1).
These results reinforce the high polymorphism of exon-3 in FLG and its ethnic association, making it difficult to generalize genomic results in relation to AD phenotypes, especially in highly mixed populations. 5,7,8 According to the literature, null variants were expected in 14% to 42% of AD patients, and, of these, 20 were found in exon-3. In addition, variations in FLG have been reported among AD patients of different ethnicities. 8 The Brazilian population has a multiracial composition, after 500 years of miscegenation amongst individuals from Western Europe, Africa, and the Amerindian populations. 9 The c.2512C>T:R826X and 2282delCAGT:S761Cfs*36 variants in AD patients corroborate these ancestries, meanwhile the c.2544T>C: synonymous variant is characteristic of this region.
Considering the epidemiological association with disease development, isolated variants in FLG gene do not fully explain the variation in AD severity, eosinophilia, or elevated IgE levels, reinforcing the multifactorial aspects of the disease. 1,8,10 In conclusion, we authors identified null variants of the FLG gene (in exon-3) (c.1537C>T:501X, c.2282del4:S761Cfs*36 and c.2512C>T:R826X) in 10% of Brazilian patients with AD, but without an association with the main clinical characteristics of AD.