Gene variants associated with acne vulgaris presentation and severity: a systematic review and meta-analysis

Background Multiple factors have been attributed to acne vulgaris predisposition and individual variations in the severity of skin symptoms, and genetics stood out as one of the major factors. Methods We performed a systematic review on the genes and their variants that have been investigated for association with acne presentation and severity. A random-effect meta-analysis using the allele model (minor allele vs. major allele) was also conducted to provide an overall estimation of risk effects of frequently reported gene variants. This included a subset data of 982 acne cases and 846 controls extracted from our existing GWAS database on various allergic and skin diseases among Singapore Chinese. Results Systematic review of 51 articles covering Asians and Caucasians found 60 genes/loci and their 100 variants implicated in acne; majority of them were in the intron, coding region/missense, and promoter regions. The commonly studied candidate genes/gene families include tumor necrosis factor (TNF), and the interleukin (IL) and cytochrome P450 (CYP) gene families. Our meta-analysis showed that most of the analyzed gene variants exhibited insignificant pooled odds ratio (pOR) and significant heterogeneity between studies. Nevertheless, we found that TNF rs1800629 A allele carriers and CYP17A1 rs743572 T allele carriers had significantly reduced mild acne risk [pOR: 0.60; 95% Confidence Interval (CI): 0.33–0.86] and severe acne risk (pOR: 0.59; 95% CI: 0.40–0.79), respectively, across populations. Overall, FST (follistatin) rs629725 A allele poses a significantly modest increased risk for acne presentation (pOR: 1.19, 95% CI: 1.14, 1.23), but neither TIMP2 (TIMP metallopeptidase inhibitor 2) rs8179090 nor CYP1A1 rs4646903 (pOR: 0.96, 95% CI: 0.80–1.12; pOR: 0.95, 95% CI: 0.83, 1.08), respectively. We discovered 15 novel SNPs in the 3′ UTR region of the Toll-like Receptor 4 gene (TLR4) associated with acne presentation. Conclusions This systematic review and meta-analysis suggest that genes influencing inflammatory responses, specifically TNF, and genes influencing the function and activity of sebaceous glands, specifically CYP17A1 and FST, have potential risk variants for acne presentation and severity across populations. Understanding the genetic susceptibility factors and biological pathways involved in the pathogenesis of acne will help us to gain insights into developing effective acne treatments. Supplementary Information The online version contains supplementary material available at 10.1186/s12920-021-00953-8.

of non-inflammatory closed and open comedones, inflammatory pustules and papules, as well as residual pathology like nodules and cysts [1,2]. Chronic acne inflammatory symptoms like scars, erythema and hyperpigmentation, lead to psycho-social consequences such as depression, anxiety [3] and unemployment [4].
The etiology of acne is a complex interplay between androgen-induced sebum production, follicular keratinization, inflammation, and colonization of pilosebaceous follicles by Cutibacterium acnes (formerly Propionibacterium acnes) [1]. Acne is a multifactorial disease, and we have recently reviewed that the epidemiological risk factors for acne and acne severity include demographics, genetics/hormonal, dietary habits and lifestyle factors [5]. Of note, earlier twin [6] and family studies [7,8] have established that acne susceptibility has a strong genetic component. Multiple case-control familial studies and twin studies involving numerous ethnicities showed strong heritability with estimates upwards of 78% [6][7][8], while concordance between monozygotic twins compared with dizygotic twins was higher for both acne presentation and severity [8]. Recent candidate gene studies and genome-wide association studies (GWAS) further revealed that the genes and loci associated with acne presentation and severity influence the function and activity of sebaceous glands or immune and inflammatory responses (reviewed in 9,10). Two such gene susceptibility loci were found in Han Chinese [9], while up to 15 were found in Caucasians [10][11][12].
Reviews on the genes and gene variants associated with acne presentation and severity are scarce [13][14][15]. Therefore, we hereby performed a systematic review on that, and also performed a random-effect meta-analysis on selected risk gene variants, which included a subset data of 982 acne cases and 846 controls extracted from our existing GWAS database on various allergic and skin diseases among Singapore Chinese. Lastly, we discussed the biological pathways implicated by these genes in the pathogenesis of acne.

Literature search
The articles reviewed were retrieved from searches conducted on the Web of Science database on 19 October 2020. In the first search, the search criteria used were: document type 'article' , search terms 'acne' in the topic and 'polymorphism' in the title. Additional searches using the following criteria-document type 'article' and either the search terms 'acne' and 'predisposition' in the topic or search terms 'acne risk' and 'gene' in the topic-were also conducted. The searches aimed to garner articles about polymorphisms associated with acne presentation (presence/absence of acne vulgaris, severe acne, or teenage acne) and acne severity (mild, moderate, severe grades), thus more general search terms were chosen. Monogenic disorders and syndromes associated with (severe) acne (such as Apert syndrome, Frank-ter Haar syndrome and Winchester syndrome), and genes and gene variants associated with co-morbidities of acne such as endocrine system syndromes were not included in the systematic review as they have been reviewed recently [15]. Quality of included studies was evaluated using the JBI Critical Appraisal Tool Checklist [16] containing eight criteria. All papers selected for inclusion in the systematic review (that is-those that meet the inclusion criteria described above) were subjected to rigorous appraisal by two critical appraisers.

Criteria for meta-analysis
Studies were included in the meta-analysis if they satisfied the following inclusion criteria: 1. Gene variants that were involved in at least two case-control studies that evaluated acne presentation and/or acne severity; 2. Genotype/allele frequencies were available for cases and controls; 3. The distribution of genotypes in the control group fulfils the Hardy-Weinberg equilibrium (HWE); 4. Studies that have provided an estimation of effect size, such as the odds ratio (OR) with corresponding 95% confidence intervals (CI); or if not reported in article, derived from statistical analysis using genotype/allele frequencies or from communication with corresponding authors; 5. When publications involved the overlapping data sets, only the study with the largest number of participants was included. The supporting PRISMA checklist is available as Additional file 1: Table S1.

Data extraction and retrieval
The following data were extracted: the name of the first author, year of publication, country of origin, ethnicity, genotype frequencies in acne cases and controls or acne severity grades, presence/absence of association, OR, CI and p values for genotypes and alleles. Additional information about the gene variants were retrieved from National Centre for Biotechnology Information (NCBI) dbSNP or European Bioinformatics Institute (EBI) Ensembl databases: chromosomal location, common name of variant, rs number (if available), most severe consequence. Availability of expression quantitative trait loci (eQTL) data for the SNPs and the tissue types was sourced from the GTex Portal (available at https:// gtexp ortal. org/ home/). Gene ontology analysis/network analysis was performed for the list of genes using the online Database for Annotation, Visualization and Integrated Discovery (DAVID) v6.8 software (available at https:// david. ncifc rf. gov/), with the highest classification stringency and other default settings for functional annotation clustering.
A subset of genotyping data from our existing GWAS database on various allergic and skin diseases among Singapore Chinese were extracted for 982 acne cases and 846 controls (full details in Additional file 2). After exclusion, a total of 4517 SNPs were analyzed and logistic regression analysis was conducted using the plink software to investigate the association between the minor allele counts and acne presentation.

Statistical analysis for meta-analysis
To perform the random-effect meta-analysis, we extracted the OR and 95% CI reported from each study of interest. Where not reported, OR and CI were calculated with binary logistic regression test using IBM ® SPSS ® Statistics software (IBM Inc., NY). For consistency, we estimated the association based on the allele model (minor allele vs. major allele). These study findings were combined using the random effect model and the pooled OR and 95% CI were computed using the Stata/MP ® version 16.0 statistical software (StataCorp LLC, TX). We used a chi-square-based test to examine any heterogeneity presented in the pooled risk estimate, with the inconsistency index (I 2 ) also computed; I 2 ≥ 50%, p < 0.05 was considered statistically significant. To evaluate whether the association showed any ethnicity-specific or severity-specific effects, we analyzed the data for separate subgroups defined by ethnicity or acne severity grade. Sensitivity analysis was performed by excluding individual studies and recalculating the results in order to assess the stability of the results. It was performed only for SNPs with ≥ 4 studies as we have set the criteria that the meta-analysis should only be performed for ≥ 3 studies. Publication bias was assessed using Begg's funnel plots and Egger's test.

Overview of systematic review and meta-analysis of included studies
The flow chart that displays the study selection process is shown in Fig. 1. The three searches identified 64, 45 and 115 results respectively; however, there were overlaps in the results obtained from the three searches. After excluding articles describing other diseases, types of acne and acne treatments; articles with poorly specified study design and methods; repeated articles and other irrelevant articles, 46 unique articles remained. Five crossreferencing articles were also included, thus a total of 51 articles were included in the systematic review. Excluding 26 articles that had gene variants that were involved in less than two case-control studies (excluding ours) or had OR and CI values that could not be derived or not reported by authors, 25 articles were included in the meta-analysis.
Sixty genes and their 100 variants investigated in the reviewed articles and their associations with acne presentation and acne severity are summarized in Table 1. The most common region where the variants are located is the intron, followed by the coding region-missense, promoter, intergene, coding region-synonymous, 3′ UTR and non-coding exon region. Populations covered in the studies included Asians, Europeans and North Americans, but not Africans, South Americans and Oceanians. DAVID gene ontology analysis revealed 15 functional annotation clusters (Additional file 3: Figure S2). Most of the genes identified in the reviewed studies and from DAVID analysis can be categorized into two major biological function groups, namely 'immune and inflammatory responses' and 'sebaceous gland function and activity' (Table 1). Genes and gene variants from different populations studies and their differing association findings will be first described, with meta-analysis results and in vitro findings investigating the biological functional consequences of these SNPs included as well, where available.
Sensitivity analysis showed that the corresponding pooled ORs were not significantly altered by sequential omission of individual studies (data not shown). They remained significant or insignificant even after sequential omission of individual studies. The results of sensitivity analysis indicated the stability of our results. Begg's funnel plot and Egger's test showed that there was statistically significant publication bias for the meta-analysis of TNF rs1800629 acne presentation, mild acne and severe acne; SELL rs7531806; PPARG rs1801282; and CYP17A1 rs743572 acne presentation, moderate acne and severe acne. However, the shape of the funnel plots for the rest of the SNPs analyzed seemed symmetrical and Egger's test p value > 0.05, indicating that there were no evidences for obvious publication bias (Additional file 3: Figure S3).

Genes and gene variants involved in immune and inflammatory responses Tumor necrosis factor (TNF)
Some of the genes implicated in acne presentation and acne severity are involved in immune and inflammatory responses. The pro-inflammatory factor, tumor necrosis factor (TNF) is the most frequently studied gene in the reviewed articles. In particular, the studies investigated variations at positions − 238, − 308, − 857, − 863 and − 1031 localized in the promoter region of TNF, meaning that they may affect the gene expression positively or negatively, conferring a protective or detrimental effect on acne.
The − 308 SNP (rs1800629) is the most common among all and is a potential risk factor among Caucasians, as shown by two meta-analyses [37,62], but the conclusions are different in other populations. Several studies found a significant association with acne alone [24,36,37] or both acne and acne severity [33], but other studies found insignificant associations with acne severity alone [24,35,36], acne alone [34] or both acne and acne severity [22]. Furthermore, the results of one study suggested the possibility of gender differences in the association of − 308 SNP with acne as a significant association was observed for females, but not for males or the combined  Petridis et al. [11]  Petridis et al. [11] LGR6 1q32. Petridis et al. [11]  Zhang et al. [12] population [35]. Meta-analysis of these and our study for acne presentation showed that the pooled OR was 0.97 (95% CI: 0.82-1.11) (Fig. 2a) (Fig. 2b). No significant heterogeneity was found for all severity grades (I 2 < 50%, p > 0.05), suggesting that rs1800629 A allele carriers had significantly reduced mild acne risk compared with G allele carriers across all populations. Studies that investigated − 238 SNP (rs361525) showed mixed results. The − 238 SNP was significantly associated with acne and acne severity among Pakistanis [33], while no association with acne or with both acne and acne severity was reported among Caucasians in Poland [34] and Hungary/Romania [35]. Meta-analysis of these and our study showed that this SNP was not significantly associated with acne, as the pooled OR was 1.24 (95% CI: 0.85-1.63), with no significant heterogeneity observed (I 2 = 0%, p = 0.585; Fig. 2c). Lastly, TNF − 238, − 308 and − 376 haplotypes did not show association with acne risk [37].

Genes involved in both inflammatory and immune responses and sebaceous gland function and activity
The − 857 SNP (rs1799724) was significantly associated with acne and acne severity among Caucasians [12,35], while − 863 and − 1031 SNPs were not [34]. In the case of rs1799724, the major C allele exhibited a positive association with acne, whereas the minor T allele seemed to have a protective effect [35]; consistent with studies on other chronic inflammatory diseases [63]. The C/T base change generates a novel transcription factor binding site (OCT-1) at the promoter region of TNF, immediately next to a pre-existing NF-κB binding site, resulting in an altered TNF regulation in response to various stimuli [64]. This suggests the general importance of this SNP in regulating TNF expression.
In addition, two studies which investigated the M196R variation in TNFR2 (Tumor Necrosis Factor Receptor 2) found a significant association with acne presentation and severity among Han Chinese [32], but not among Caucasians [12]. Several studies found that the IL1A − 889 C/T SNP (rs1800587) was associated with acne among Asian and European populations [18][19][20] while + 4845 G/T SNP only showed association among Caucasians in Hungary/ Romania [21], but not in USA [12]. Meta-analysis of rs1800587 studies and our study showed that this SNP was significantly associated with acne, as the pooled OR was 1.34 (95% CI: 1.14-1.54), with significant heterogeneity observed (I 2 = 77.3%, p = 0.004; Fig. 3a). Meanwhile, IL1A + 4845 G/T SNP causes an Ala114Ser substitution close to the proteolytic cleavage site of nuclear pre-IL-α (117/118) to mature IL-α, and might result in enhanced cleavage when the rare T allele is present [65,66]. Pre-IL-1α has a predominantly nuclear localization, whereas mature IL-α exhibits a cytoplasmic localization and can be secreted in response to the appropriate signals [67]. Therefore, carriers of the T allele may have a higher risk of more severe acne symptoms that result from enhanced inflammatory reactions by increased IL-α secretion [13,22].

Interleukins (IL) and their associated antagonists and receptors
In contrast, IL1B − 511 C/T, IL4 − 590 T/C, IL10 − 1082 A/G SNPs and IL1RN variable number tandem repeats (VNTR) all showed no significant association with acne and/or acne severity in the studied populations [21][22][23][24]. Furthermore, a study in Saudi Arabia found that the IL4R Q551R SNP was significantly associated with acne but not acne severity [23].
Two studies also found a significant association between IL6 − 572 G/C variants and acne presentation [19,60]. Meta-analysis of these and our study showed (See figure on next page.) Fig. 2 Meta-analysis of TNF SNPs associated with acne presentation and acne severity grades. a Subgroup meta-analysis of association of rs1800629 with acne presentation among Asians and Caucasians. b Subgroup meta-analysis of association of rs1800629 with acne severity among different grades. c Meta-analysis of association of rs361525 with acne presentation. Only gene variants that were investigated in at least two other previous studies were included in the meta-analysis. Analysis was performed under allele model (minor allele vs. major allele)., i.e. A versus G for both rs1800629 and rs361525 Heng et al. BMC Med Genomics (2021) Fig. 3b). Lastly, IL8 − 251 T/A SNP was found to be significantly associated with acne and acne severity among Pakistanis [61], but not among Polish [18]; meta-analysis of these and our study showed that this SNP was not significantly associated with acne, as the pooled OR was 1.03 (95% CI: 0.91-1.15), with no significant heterogeneity observed (I 2 = 61.2%, p = 0.076; Fig. 3c). Lastly, Ehm et al. [26] found that IL17A rs2275913, IL17F rs763780, IL17RB rs1043261 and IL23R rs11209026 all had no significant association with both severe acne and teenage acne.

Toll-like receptor (TLR) family
Another gene family implicated is the toll-like receptor (TLR) gene family. None of the TLR4 variants (rs4986790 and rs4986791) investigated in the reviewed articles were found to be associated with acne risk [12,31]. Both of these SNPs change amino acids in the ligand recognition part of the receptor, but most in vitro functional studies suggest that the LPS-induced cytokine response of the Asp299Gly/Thr399Ile haplotype does not differ from that of the wild-type cytokine response [reviewed in [68]]. Nevertheless, there was a study which found increased tumor necrosis factor-α (TNF-α) cytokine response of the African-exclusive Asp299Gly/wild-type haplotype [69]. Using a subset of genotyping data of 982 acne cases and 846 controls extracted from our existing GWAS database (Additional file 2), we found that out of all the gene variants in this systematic review (Table 1), only TLR4 (Additional file 2: Table S2.1) was significantly associated with acne presentation. Specifically, within 2 kb of this gene, 15 SNPs located in the 3′ UTR region (Additional file 2: Fig. S1), were significantly associated with acne presentation (p = 0.021222 for all SNPs; Additional file 2: Table S2.2).
Furthermore, no significant association with acne was found for two TLR2 variants + 2179 C/T and + 2258 G/A among Caucasians [12,31]. However, the latter SNP was significantly associated with acne and acne severity among Han Chinese [32]. The possible explanation in the discrepancy in the findings is that a highly homologous duplicated segment of exon 3 (harboring the SNP) exists 23 kb upstream of the TLR2 locus [70]. In vitro data showed that the Arg to Gly change may impair TLR2mediated immune signaling and the expression of various downstream target genes in response to microbial ligands [71,72], suggesting that if this SNP may be a functional one in populations where it is truly present with reasonable frequency.

Genes and gene variants involved in sebaceous gland function and activity
Some genes implicated in acne and severe acne risk may influence the function and activity of the sebaceous gland. Two frequently studied gene families that have been implicated are the cytochrome P450 (CYP) gene family and the 3-beta hydroxysteroid dehydrogenase/isomerase (HSD3B) gene family.
For CYP17A1, a Chinese study found that the − 34 T/C SNP was not significantly associated with both acne in the overall population and acne severity in females, but was significantly associated with acne severity in males, suggesting gender differences in acne severity genetic predisposition [45]. Another CYP17A1 SNP rs743572 was investigated in three studies-two studies reported significant association with acne and acne severity [43,44], while another did not [12]. Meta-analysis of these and our study for acne presentation showed that this SNP was not significantly associated with acne, as the pooled OR was 1.03 (95% CI: 0.91-1.15) (Fig. 5a), while a subgroup analysis revealed an OR of 0.95 (95% CI: 0.68-1.23), 1.08 (95% CI: 0.76-1.40), and 0.59 (95% CI: 0.40-0.79) among mild, moderate and severe grades, respectively (Fig. 5b). Significant heterogeneity was observed across these reported findings, suggesting these analyzed findings had different study outcomes across each other (all I 2 ≥ 50%, p < 0.05). These heterogeneities also correspond to the publication bias detected for this Only gene variants that were investigated in at least two other previous studies were included in the meta-analysis. Analysis was performed under allele model (minor allele vs. major allele), i.e. A versus G for a, T vs. A for b, C versus G for c, G versus A for d, A versus G for e, 192/192 versus non192/non192 for f, and G versus C for g. OR and CI values for a, d and e for Petridis et al. [11] were derived from the authors' meta-analysis. OR and CI values for a, d for Navarini et al. [10] were derived from personal communication with the authors, while for e was derived from [11] Fig. 5 Meta-analysis of CYP SNPs associated with acne presentation and acne severity grades. a Meta-analysis of association of CYP17A1 rs743572 with acne presentation. b Subgroup meta-analysis of association of rs743572 with acne severity among different grades. c Meta-analysis of association of CYP1A1 rs4646903 with acne presentation. Only gene variants that were investigated in at least two other previous studies were included in the meta-analysis. Analysis was performed under allele model (minor allele vs. major allele)., i.e. T versus C for rs743572 and C versus T for rs464903 (See figure on next page.) SNP. Although this SNP may affect CYP17 regulation by creating a putative novel SP1 promoter-binding site [74], in vitro EMSA experiments have so far not showed any evidence on this [75].

HSD3B family
Three members of the 3-β HSD family were investigated in the reviewed studies. A Chinese study found that HSD3B1 (Hydroxy-Delta-5-Steroid Dehydrogenase, 3 Beta-and Steroid Delta-Isomerase 1) rs6428829 and AAT haplotype, and HSD17B3 (Hydroxysteroid 17-Beta Dehydrogenase 3) H8 haplotype were significantly associated with acne risk [53]. In contrast, a study in Egypt found that HSD11B1 rs12086634 variants were associated with acne risk while the relationship between HSD11B1 rs846910 variants and acne risk was unclear [52].

Other genes
Other frequently studied genes thought to influence sebaceous gland function and activity include the androgen receptor (AR), peroxisome proliferator activated receptor gamma (PPARG) and insulin like growth factor 1 (IGF-1) genes. For AR, a Chinese study observed that GGN (Gly) repeat VNTR was not significantly associated with acne risk [39]. However, for the CAG (Gln) repeat VNTR, a significant association with acne risk was observed in one study [39] while other studies showed that it was significantly associated with acne and acne severity in males, but not in females [40] and vice versa [41]. This suggests that there may be gender differences in the association between AR CAG repeat VNTR and both acne and acne severity. Biologically, these VNTRs cause variation in the lengths of polyGln and polyGly stretch in the N-terminal domain of AR receptor. The variations in repeat lengths in both VNTRs have been shown to be associated with subtle modulations of AR expression, resulting in the modified transcriptional activity of various downstream targets. The transcriptional activity of AR is inversely correlated with number of the polyGln repeat length; shorter alleles exhibit greatest activity [76][77][78]. Furthermore, in vitro results suggest that acne patients carrying fewer CAG repeats may exhibit a higher AR mRNA and protein expression, leading to higher sensitivity to androgens than in control individuals [79].
For PPARG , a study in Iran found that rs3856806 and rs1801282 SNPs were not associated with acne and acne severity in the overall population [56]. However, the study found that rs1801282 was associated with acne risk in a subset of patients who developed acne at the age of 20 or older. Another study observed that the Pro/ Ala genotype of rs1801282 appeared to be protective for acne and more severe acne [55]. However, our metaanalysis revealed that the rs1801282 G allele do not pose any increased/decreased risk for acne [pooled OR: 1.03, 95% CI: 0.76, 1.29; significant heterogeneity observed (I 2 = 66.6%, p = 0.050); Fig. 4g]. This heterogeneity also corresponds to the publication bias detected for this SNP. Similarly, with reference to IGF1, the length of CA repeats was found to be significantly associated with acne risk only [49], and the risk of both acne and more severe forms of acne [50,51]. However, our meta-analysis revealed that the IGF1 VNTR 192/192 genotype do not pose any increased/decreased risk for acne [pooled OR: 1.02, 95% CI: 0.08, 1.97; no significant heterogeneity observed (I 2 = 0%, p = 0.439); Fig. 4f ].
In addition, the following variants were significantly associated with acne risk: ADH7 (alcohol dehydrogenase 7) rs1154469 [38], SRD5A2 (steroid 5 alpha-reductase 2) TA repeat VNTR [57] and VDR (vitamin D receptor) rs731236 and rs7975232 [58]. However, no significant association emerged between SRD5A2 V89L SNP and acne risk [57]. All these genes have been implicated in sebaceous gland function and activity. Figure 6 summarizes the acne-implicated genes and their potential biological functions in the pathogenesis of acne. These biological functions are also part of the top DAVID functional annotation clusters, as analyzed earlier. Notably, some of these genes have overlapping functions. The following parts will further discuss how these acne-implicated genes play plausible roles in the biological pathways influencing inflammation, sebaceous gland function and activity.

Immune and inflammatory responses
A contributing factor to the pathogenesis of acne is C. acnes, a species of bacteria commonly found in pilosebaceous units that is thought to promote the development of acne lesions [80]. Studies suggest that C. acnes can activate the NLRP3 inflammasome pathway, eventually resulting in elevated IL-1β levels and contributing to the pathogenesis of acne [81]. In addition, TGFβ may not only have roles in sebocyte lipid synthesis, but may also regulate immune responses against C. acnes. Furthermore, C. acnes can be recognized by pattern recognition receptors such as TLR2 and TLR4, which mount immune responses against the bacteria [80], thus influencing acne development. Upon recognizing C. acnes, TLRs can activate cellular components such as TNF and IL-1A, which in turn activate chemokines such as IL-8, IL-10 and downstream components that eventually promote inflammation [13,24]. IL-1A activity may also be modulated by IL1RN [14].
Furthermore, TNF can interact with its receptors, such as TNFR2, activating intracellular signaling cascades. These cascades can activate MAPKs such as MAPK11 that upregulate the expression of pro-inflammatory mediators, and also activate matrix MMPs such as MMP-2 to bring about tissue remodeling [82]. MMP activity is regulated by TIMPs; altering the balance of MMP and TIMP activity may contribute to acne development [27]. Similarly, IL-6 (which could be stimulated by IL17A, IL17F and ACE) is reported play a role in modulating host defense mechanisms [14] and inflammation [83].
Pro-inflammatory factors TNF and IL-6 also have roles in modulating RETN transcription and translation. Studies demonstrate a bidirectional relationship between resistin and TNF and IL-6: TNF and IL-6 not only modulate resistin expression but their activity is in turn modulated by resistin levels [83]. Variation in any of these genes can affect their expression, function or activity, which may predispose individuals to acne or severe acne. EGF/EGFR signaling plays an important role in the regulation of inflammation by inhibiting proinflammatory cytokines including IL-1α, IL-8 and TNF-α in C. acnes stimulated epidermal keratinocytes [84].
While the specific roles of IL-4 and IL-4R in the pathogenesis of acne are unclear, these proteins are known to play a role in inflammatory responses. Thus it is hypothesized that IL-4 and IL-4R are involved in the inflammatory processes associated with acne [23]. In addition, studies have shown that adipokines secreted by adipocytes may play a role in inflammation associated with acne [25]. The adipokine omentin may thus influence inflammatory processes that contribute to the pathogenesis of acne, such as by preventing inflammation caused by TNF [25]. Similarly, L-selectin is expressed in certain classes of immune cells and plays a role in immune cell invasion of inflamed tissue [9]. Evidence suggests that SELL, which encodes L-selectin, may contribute to inflammatory processes associated with more severe acne [14]. In addition, while the role of SEMA4B in acne pathogenesis is uncertain, it is known to be involved in immune responses and downregulates the secretion of IL-4 and IL-6 in basophils [85]. Thus, it is possible that SEMA4B may also play a role in acnerelated immune processes.
Finally, two recent bioinformatics studies using gene expression data of the lesional and nonlesional skin of acne patients from the NCBI Gene Expression Omnibus found 12 top upregulated differentially expressed genes (DEGs) involved in immune and inflammatory responses, namely FPR2, ITGB2, CXCL8, C3AR1, CXCL1, FCER1G, LILRB2, PTPRC, SAA1, CCR2, ICAM1, and FPR1 [86,87]. Specifically, the putative pathways involved in acne pathogenesis include chemokine signaling, cytokinecytokine receptor interaction, and Fc gamma R-mediated phagocytosis [86,87]. Of note, neither of the variants in these genes have been investigated in candidate gene studies nor were picked up by GWAS for association with acne presentation and severity, thus warrants future investigation.

Sebaceous gland function and activity
Sebaceous glands, which are involved in the biosynthesis and secretion of sebum, are present on almost the entire body surface, and particularly on the face [88]. Changes in the activity of these glands may lead to changes in skin physiology and have been implicated in the pathogenesis of acne. Increased sebaceous gland activity can result in higher sebum levels, which may hamper normal cellular processes in the skin and eventually lead to the development of acne lesions [88].
Some of the genes found to be associated with acne and acne severity in this review may have roles in the sebaceous gland function and activity, and thus variations in these genes may affect cell signaling and functioning in sebaceous glands, thus contributing to acne development or exacerbating acne. One cell-signaling pathway in sebaceous gland cells (sebocytes) is the VDR pathway. Sebocytes have been reported to express VDR, and studies have shown that VDR activation affects sebocyte growth and division, lipid production and levels of ILs, such as IL-6 and IL-8 [89]. Thus, variants in VDR, IL6 and IL8 may affect sebaceous gland function and contribute to acne development or exacerbation.
Two other pathways, the AR and PPAR signaling pathways, may also regulate sebocyte function and activity. Sebocytes have important functions in androgen metabolism, and high androgen levels can upregulate the activity of sebocytes. The effects of androgens on the sebaceous gland are thought to be mediated via AR and PPAR signaling pathways [88]. AR activity may, in turn, be regulated by DDB2. The protein encoded by DDB2 facilitates interactions between AR and a complex responsible for targeting AR for ubiquitination and degradation [14]. In addition, variations in genes involved in androgen synthesis can affect androgen levels and thus influence sebocyte activity. The gene SRD5A2 [57], some members of the CYP gene family, such as CYP17A1 [13], and some members of the HSD gene family encode enzymes that have roles in androgen synthesis [53]. For example, the enzymes encoded by HSD3B1 and HSD17B3 are involved in the biosynthesis of testosterone [53].
In addition, variations in genes such as IGF1 and TGFB are suggested to influence the production of lipids in sebocytes, and may contribute to acne development [14]. Notably, dietary variables such as intake of fish, fruits and vegetables were suggested to influence acne and severe acne risk by modulating IGF-1 levels, showing a possible interaction between genetic factors and modifiable risk factors [90]. In addition, the activity of the TGFβ signaling pathway is regulated by the transcription factor OVOL1 and the protein follistatin (encoded by FST), thus variations in these genes may influence sebocyte lipid synthesis and thus play a part in acne development and exacerbation [14]. Navarini et al. [10] also revealed that transcript levels of TGFB2 and OVOL1 were significantly decreased in fresh inflammatory acne papules compared to normal skin. MUC1 is another gene expressed in sebaceous glands and it encodes for a glycoprotein that may be involved in inhibiting the adhesion of bacteria onto the skin, including the acne-causing bacteria C. acnes [54]. Variations in MUC1 may influence the anti-bacterial properties of the glycoprotein, thus influencing acne development.
Furthermore, reports have suggested that some genes implicated in the pathogenesis of acne may be involved in the development of sebaceous glands [91]. One of these genes is WNT10A. WNT10A activity can influence sebum production in the sebaceous gland [11] and may also regulate the differentiation of progenitors to form sebocytes [91]. Thus, changes in Wnt signaling may affect the normal development and functioning of the sebaceous gland and in turn, the risk of acne or severe acne. Wnt signaling is regulated by LGR6, and high LGR6 expression can inhibit sebaceous progenitor differentiation [91]. SPECC1L is also reported to be a component of the Wnt signaling pathway [92] and may be implicated in acne development. Another gene involved in the development of sebaceous glands is LAMC2, which facilitates the migration of sebaceous gland progenitor cells to the site of the gland [91]. Variations in these genes can affect sebaceous gland development, contributing to acne development and exacerbation.
Recently, ADH7 rs1154469 has been identified as a novel susceptibility locus for severe acne among Han Chinese [38]. The most active enzyme encoded by this gene is retinol dehydrogenase; thus it may participate in the synthesis of retinoic acid. Retinoic acid is a retinoid which acts to normalize desquamation within the sebaceous follicles that leads to obstruction of the pilosebaceous canal as seen in acne, by reducing keratinocyte proliferation and promoting differentiation [93].
Notably, environments with high sebum levels favor the growth of C. acnes and promote the colonization of pilosebaceous units by C. acnes [94]. This suggests a link between the two implicated pathways: changes in the sebaceous gland that lead to increased sebum production may promote C. acnes infection, triggering inflammatory and immune responses that contribute to acne pathogenesis.

Other genes and potential pathways implicated in acne
The rest of the genes identified in the reviewed articles were not implicated in inflammation, immune responses or sebaceous gland function and activity. The genes and SNPs significantly associated with severe acne found by Petridis et al. [11] were: BCL11A (BAF chromatin remodeling complex subunit BCL11A) rs2901000, FGF2 (fibroblast growth factor 2) rs4487353, GLI2 (GLI family zinc finger 2) rs1092479, LOC105378977 rs158639, PINX1 (PIN2 (TERF1) interacting telomerase inhibitor 1) rs28570522 and SUGCT (succinyl-CoA:glutarate-CoA transferase) rs7809981; and by Yang et al. [38]: F13A1 (Coagulation Factor XIII A Chain) rs435048. Furthermore, a few SNPs that were not located within a gene were investigated in the reviewed studies. The SNPs rs2727365, rs4133274 and rs13248513 SNPs showed a significant association with severe acne risk [11,12]; the mechanism(s) by which variations at these polymorphisms influence severe acne risk are unclear.
It is also unclear how PINX1, which has functions in telomere length maintenance and chromosome stability [95]; F13A1, which encodes coagulation factor XIII, the last zymogen to become activated in the blood coagulation cascade and a transglutaminase enzyme [96]; BCL11A, a transcription factor involved in hematopoietic development [97]; FGF2, which is involved in processes that contribute to scar formation and wound healing [11]; and SUGCT , which encodes an enzyme that converts glutarate to glutaryl-CoA [98] may influence the pathogenesis of acne.
Genes implicated in acne presentation and severity may influence the pathogenesis of acne through mechanisms other than altering sebaceous gland function or inflammatory and immune responses. For instance, the transcription factor GLI2 regulates the expression of sonic hedgehog (Shh), a component of signaling pathways that influence the development of hair follicles [99]. A previous study reported that several genes involved in hair follicle development were implicated in severe acne risk, suggesting that this pathway may contribute to the pathogenesis of acne [11].

Conclusions
The reviewed papers reported on genes and their variants associated with the risk of acne presentation and severity. However, the results of different studies show inconsistencies, with some studies reporting that a particular variant increases the risk, while others did not. Out of the 51 studies reviewed, only five (including ours) were GWAS while others where candidate gene studies. Hence, we acknowledge the poor reproducibility of association studies, which could be attributed to factors such as insufficient availability of genetic markers, inadequate handling of population structure, lack of statistical power due to low sample size, improper control of multiple testing and extensive publication bias [reviewed in [100]]. Furthermore, it is also possible that rather than particular SNPs, it is the haplotype of multiple SNPs or changes in gene expression or function that is associated with acne presentation or severity. Rather than considering the impact of individual variants on acne presentation or severity, it may be more helpful to consider the gene as a functional unit and to understand the effects of altered gene function on acne presentation and severity.
In conclusion, this review summarizes the literature on candidate genes implicated in the risk of acne presentation and severity and possible mechanisms by which these genes may affect acne pathogenesis. Notably, a large majority of the candidate genes identified are suggested to have roles in the function and activity of sebaceous glands or immune and inflammatory responses-in line with the literature that describes acne as a chronic inflammatory disease of the pilosebaceous unit. Understanding the genetic susceptibility factors and biological pathways involved in the pathogenesis of acne will help us to gain insights into developing effective acne treatments.