Assessment of the genetic basis of rosacea by genome-wide association study

Anne Lynn S Chang, Inbar Raber, Jin Xu, Rui Li, Robert Spitale, Julia Chen, Amy K Kiefer, Chao Tian, Nicholas K Eriksson, David A Hinds, Joyce Y Tung, Anne Lynn S Chang, Inbar Raber, Jin Xu, Rui Li, Robert Spitale, Julia Chen, Amy K Kiefer, Chao Tian, Nicholas K Eriksson, David A Hinds, Joyce Y Tung

Abstract

Rosacea is a common, chronic skin disease that is currently incurable. Although environmental factors influence rosacea, the genetic basis of rosacea is not established. In this genome-wide association study, a discovery group of 22,952 individuals (2,618 rosacea cases and 20,334 controls) was analyzed, leading to identification of two significant single-nucleotide polymorphisms (SNPs) associated with rosacea, one of which replicated in a new group of 29,481 individuals (3,205 rosacea cases and 26,262 controls). The confirmed SNP, rs763035 (P=8.0 × 10(-11) discovery group; P=0.00031 replication group), is intergenic between HLA-DRA and BTNL2. Exploratory immunohistochemical analysis of HLA-DRA and BTNL2 expression in papulopustular rosacea lesions from six individuals, including one with the rs763035 variant, revealed staining in the perifollicular inflammatory infiltrate of rosacea for both proteins. In addition, three HLA alleles, all MHC class II proteins, were significantly associated with rosacea in the discovery group and confirmed in the replication group: HLA-DRB1*03:01 (P=1.0 × 10(-8) discovery group; P=4.4 × 10(-6) replication group), HLA-DQB1*02:01 (P=1.3 × 10(-8) discovery group; P=7.2 × 10(-6) replication group), and HLA-DQA1*05:01 (P=1.4 × 10(-8) discovery group; P=7.6 × 10(-6) replication group). Collectively, the gene variants identified in this study support the concept of a genetic component for rosacea, and provide candidate targets for future studies to better understand and treat rosacea.

Figures

Figure 1
Figure 1
Manhattan plot of single-nucleotide polymorphisms in the rosacea genome-wide association study of the discovery group. Single-nucleotide polymorphisms meeting a Bonferroni threshold of significance at P<5 × 10−8 (gray line) were rs763035 and rs111314066; however, only rs763035 was confirmed in the replication group (Table 2). Results are plotted as −log10-transformed P-values from the association test.
Figure 2
Figure 2
Regional association plot for rs763035, the single-nucleotide polymorphism confirmed in the replication group. The y-axis shows negative transformed −log10 of association test P-values, and x-axis shows the adjacent region on the chromosome. Colors indicate the strength of linkage disequilibrium with the index single-nucleotide polymorphism, and plot symbols indicate imputed (+) and genotyped (o) single-nucleotide polymorphisms. The human leukocyte antigen association is supported by numerous genotyped and imputed variants spanning an interval of about 1.5 Mb.
Figure 3
Figure 3
Genes associated with rs763035 are present in rosacea skin. (a) Rosacea skin from an individual of European descent and heterozygous for the minor allele at rs763035 immunostained with human leukocyte antigen-DRA antibody with a strong signal in the perifollicular inflammatory infiltrate (black and white arrow), epidermal Langerhans cells (single black arrow), and endothelial cells (double arrows). (b) Rosacea skin from the same individual as a and b immunostained with BTNL2 antibody showed mild–moderate staining in keratinocytes, perifollicular infiltrate (black and white arrow), and endothelial cells (double arrow). (c) Normal skin from the same individual stains with human leukocyte antigen-DRA antibody but lacks perifollicular infiltrate, numerous epidermal Langerhans cells, and dilated blood vessels seen in rosacea skin. (d) Normal skin from the same individual stained with BTNL2 antibody in keratinocytes but does not contain inflammatory infiltrate or dilated blood vessels. Negative controls were performed for all the tissues above by omitting either human leukocyte antigen-DRA or BTNL2 primary antibodies, and these displayed no staining (data not shown). Scale bar=100 μm.

References

    1. Abram K, Silm H, Maaroos H-I, et al. Risk factors associated with rosacea. J Eur Acad Dermatol Venereol. 2010;24:565–571.
    1. Anderson CA, Boucher G, Lees CW, et al. Meta-analysis identifies 29 additional ulcerative colitis risk loci, increasing the number of confirmed associations to 47. Nat Genet. 2011;43:246–252.
    1. Browning SR, Browning BL. Rapid and accurate haplotype phasing and missing data inference for whole genome association studies using localized haplotype clustering. Am J Hum Genet. 2007;81:1084–1097.
    1. Cribier B. Medical history of the representation of rosacea in the 19th century. J Am Acad Dermatol. 2013;69:S2–13.
    1. Dimas AS, Deutsch S, Stranger BE, et al. Common regulatory variation impacts gene expression in a cell type-dependent manner. Science. 2009;325:1246–1250.
    1. Fowler J, Jarratt M, Moore A, et al. Once-daily topical brimonidine tartrate gel 0.5% is a novel treatment for moderate to severe facial erythema of rosacea: results of two multicenter, randomized and vehicle-controlled studies. Br J Dermatol. 2012;166:633–641.
    1. Elewski BE, Draelos Z, Dreno B, et al. Rosacea—global diversity and optimized outcome: proposed international consensus from the Rosacea International Expert Group. JEADV. 2011;25:188–200.
    1. Erlich H, Valdes AM, Noble J, et al. HLA DR-DQ haplotypes and genotypes and type 1 diabetes risk: analysis of the type 1 diabetes genetics consortium families. Diabetes. 2008;57:1084–1092.
    1. Freidlin B., Zheng G, Li Z, et al. Trend tests for case-control studies of genetic markers: power, sample size, and robustness. Hum Hered. 2002;53:146–152.
    1. Gibbs JR, van der Brug MP, Hernandez DG, et al. Abundant quantitative trait loci exist for DNA methylation and gene expression in human brain. PLoS Genet. 2010;6:e1000952.
    1. Henn BM, Hon L, Macpherson JM, et al. Cryptic distant relatives are common in both isolated and cosmopolitan genetic samples. PLoS One. 2013;7:e34267.
    1. Howie B, Fuchsberger C, Stephens M, et al. Fast and accurate genotype imputation in genome-wide association studies through pre-phasing. Nat Genet. 2012;44:955–959.
    1. International Multiple Sclerosis Genetics Consortium Risk alleles for multiple sclerosis identified by a genomewide study. New Eng J Med. 2007;357:851–862.
    1. Kastelan S, Tomic M, Salopek-Rabatic J, et al. The association between the HLA system and retinopathy development in patients with type 1 diabetes mellitus. Coll Antropol. 2013;37 (Suppl 1:65–70.
    1. Krini M, Chouliaras G, Kanariou M, et al. HLA class II high-resolution genotyping in Greek children with celiac disease and impact on disease susceptibility. Pediatr Res. 2012;72:625–630.
    1. Montgomery SB, Sammeth M, Gutierrez-Arcelus M, et al. Transcriptome genetics using second generation sequencing in a Caucasian population. Nature. 2010;464:773–777.
    1. Rallis E, Korfitis C. Isotretinoin for the treatment of granulomatous rosacea: a case report and review of the literature. J Cutan Med Surg. 2012;16:438–441.
    1. Schadt EE, Molony C, Chudin E, et al. Mapping the genetic architecture of gene expression in human liver. PLoS Biol. 2008;6:e107.
    1. Sirota M, Schaub MA, Batzoglou S, et al. Autoimmune disease classification by inverse association with SNP alleles. PLoS Genet. 2009;5:e1000792.
    1. Spoendlin J, Voegel JJ, Jick SS, et al. Risk of rosacea in patients with diabetes using insulin or oral antidiabetic drugs. J Invest Dermatol. 2013;133:2790–2793.
    1. Steinhoff M, Schauber J, Leyden JJ. New insights into rosacea pathophysiology: a review of recent findings. J Am Acad Dermatol. 2013;69:S15–S26.
    1. Stranger BE, Nica AC, Forrest MS, et al. Population genomics of human gene expression. Nat Genet. 2007;39:1217–1224.
    1. Valentonyte R, Hampe J, Huse K, et al. Sarcoidosis is associated with a truncating splice site mutation in BTNL2. Nat Genet. 2005;37:357–364.
    1. Weinstock L.B., Steinhoff M. Rosacea and small intestinal bacterial overgrowth: prevalence and response to rifaximin. J Am Acad Dermatol. 2013;68:875–876.
    1. Yamasaki K, Di Nardo A, Bardan A, et al. Increased serine protease activity and cathelicidin promotes skin inflammation in rosacea. Nat Med. 2007;13:975–980.
    1. Yamasaki K, Kanada K, Macleod DT, et al. TLR2 expression is increased in rosacea and stimulates enhanced serine protease production by keratinocytes. J Invest Dermatol. 2011;131:688–697.
    1. Zeller T, Wild P, Szymczak S, et al. Genetics and beyond—the transcriptome of human monocytes and disease susceptibility. PLoS One. 2010;5:e10693.
    1. Zheng X, Shen J, Cox C, et al. HIBAG-HLA genotype imputation with attribute bagging. Pharmacogenom J. 2013;14:192–200.

Source: PubMed

Подписаться