Skip to main content

ORIGINAL RESEARCH article

Front. Med., 11 May 2023
Sec. Pulmonary Medicine
This article is part of the Research Topic Diagnosis and Treatment of Sarcoidosis View all 12 articles

Sex differences in the genetics of sarcoidosis across European and African ancestry populations

Updated
\r\nYing XiongYing Xiong1Susanna Kullberg,Susanna Kullberg1,2Lori GarmanLori Garman3Nathan PezantNathan Pezant3David EllinghausDavid Ellinghaus4Vasiliki VasilaVasiliki Vasila1Anders EklundAnders Eklund2Benjamin A. RybickiBenjamin A. Rybicki5Michael C. IannuzziMichael C. Iannuzzi6Stefan Schreiber,Stefan Schreiber4,7Joachim Müller-QuernheimJoachim Müller-Quernheim8Courtney G. MontgomeryCourtney G. Montgomery3Johan Grunewald,,Johan Grunewald1,2,9Leonid Padyukov,Leonid Padyukov9,10Natalia V. Rivera,,*Natalia V. Rivera1,9,10*
  • 1Respiratory Medicine Division, Department of Medicine Solna, Karolinska Institutet, Stockholm, Sweden
  • 2Department of Respiratory Medicine and Allergy, Theme Inflammation and Ageing, Karolinska University Hospital, Stockholm, Sweden
  • 3Genes and Human Disease, Oklahoma Medical Research Foundation, Oklahoma City, OK, United States
  • 4Institute of Clinical Molecular Biology, Christian-Albrechts-University of Kiel, Kiel, Germany
  • 5Department of Public Health Sciences, Henry Ford Health System, Detroit, MI, United States
  • 6Zucker School of Medicine, Staten Island University Hospital, Northwell/Hofstra University, Staten Island, NY, United States
  • 7Clinic for Internal Medicine I, University Hospital Schleswig-Holstein, Kiel, Germany
  • 8Department of Pneumology, University Medical Center Freiburg, Faculty of Medicine, University of Freiburg, Freiburg im Breisgau, Germany
  • 9Center for Molecular Medicine, Karolinska Institutet, Stockholm, Sweden
  • 10Division of Rheumatology, Department of Medicine Solna, Karolinska Institutet, Karolinska University Hospital, Stockholm, Sweden

Background: Sex differences in the susceptibility of sarcoidosis are unknown. The study aims to identify sex-dependent genetic variations in two clinical sarcoidosis phenotypes: Löfgren’s syndrome (LS) and non-Löfgren’s syndrome (non-LS).

Methods: A meta-analysis of genome-wide association studies was conducted on Europeans and African Americans, totaling 10,103 individuals from three population-based cohorts, Sweden (n = 3,843), Germany (n = 3,342), and the United States (n = 2,918), followed by an SNP lookup in the UK Biobank (UKB, n = 387,945). A genome-wide association study based on Immunochip data consisting of 141,000 single nucleotide polymorphisms (SNPs) was conducted in the sex groups. The association test was based on logistic regression using the additive model in LS and non-LS sex groups independently. Additionally, gene-based analysis, gene expression, expression quantitative trait loci (eQTL) mapping, and pathway analysis were performed to discover functionally relevant mechanisms related to sarcoidosis and biological sex.

Results: We identified sex-dependent genetic variations in LS and non-LS sex groups. Genetic findings in LS sex groups were explicitly located in the extended Major Histocompatibility Complex (xMHC). In non-LS, genetic differences in the sex groups were primarily located in the MHC class II subregion and ANXA11. Gene-based analysis and eQTL enrichment revealed distinct sex-specific gene expression patterns in various tissues and immune cell types. In LS sex groups, a pathway map related to antigen presentation machinery by IFN-gamma. In non-LS, pathway maps related to immune response lectin-induced complement pathway in males and related to maturation and migration of dendritic cells in skin sensitization in females were identified.

Conclusion: Our findings provide new evidence for a sex bias underlying sarcoidosis genetic architecture, particularly in clinical phenotypes LS and non-LS. Biological sex likely plays a role in disease mechanisms in sarcoidosis.

Background

Sarcoidosis is a multi-system disease with unknown etiology, characterized by non-caseating granulomas (1). In 90% of cases, sarcoidosis affects the lungs and lymphatic system (2).

Sarcoidosis is prevalent worldwide, commonly affecting individuals between 20 and 50 years of age, with males diagnosed earlier than females (3, 4). On a global scale, sarcoidosis incidence and prevalence vary based on sex, age, geography, and ethnicity (4). In the United States, the prevalence is 141 per 100,000 in African–Americans and 49.8 per 100,000 in Caucasians (5). In Nordic countries such as Sweden, the prevalence of sarcoidosis was reported between 152 and 215 per 100,000 (6), whereas in Denmark was reported at 77 per 100,000 (7).

The disease is heterogeneous, as different manifestations and clinical outcomes have been observed in clinical settings, particularly among sex groups (3, 8, 9). Interestingly, sarcoidosis is more common in females, has a late-onset, has a variable disease course depending on the affected organ, and has a higher mortality rate than men (3). Sex differences are also seen in extra-pulmonary phenotypes (10, 11), including cardiac (12, 13) and skin sarcoidosis.

Sex hormones have been shown to play a role in sarcoidosis (14, 15). For instance, sarcoidosis has a low disease activity in pregnancy or goes into remission. However, the ameliorating effect is lost as flares occur after delivery (16, 17).

In the lungs, sarcoidosis is characterized by Löfgren’s syndrome (LS) and non-Löfgren’s syndrome (LS). LS is an acute form of the disease characterized by periarticular swelling around the ankles, erythema nodosum, and bilateral hilar lymphadenopathy (18). Typically, females with LS exhibit erythema nodosum (EN), while males with LS have symptoms of bilateral ankle arthritis. Patients with LS usually have a good prognosis (19). LS often occurs in individuals of European ancestry and is relatively rare in individuals of African origin. In non-LS, sarcoidosis is characterized by a heterogeneous disease course, often with an insidious onset, unrelenting disease course, and a high risk for clinical organ impairment, predominantly pulmonary fibrosis (18, 20). Patients with fibrotic sarcoidosis have markedly decreased survival compared with the general population (21, 22).

The causes for clinical differences and disease course in sarcoidosis among sex groups are unknown. However, evidence-based factors, including genetics, epigenetics, and environmental exposures, have been postulated. Regarding genetics, few studies have addressed sex differences in the genetic variation in disease, despite most studies adjusting for sex in their analyses (23). An exception is a genome-wide admixture scan conducted in African Americans stratified by sex, identifying several genetic variants associated with sarcoidosis exclusively in females (24).

In this work, we sought to investigate the genetic associations in the sex groups and identify differences and commonalities in European and African population ancestries by a meta-analysis approach and high-density mapping.

Materials and methods

Samples and study design

The samples were from three independent sarcoidosis population-based cohorts from Sweden, Germany, and the United States. All participants provided written informed consent for the study. Study protocols of all studies had been approved by respective local institutional boards.

The Swedish cohort

The case-control study consisted of 4,133 individuals. The study was approved by the local institutional review board in Stockholm, Sweden. All participants provided written informed consent and were permitted to use their DNA for research purposes. Baseline clinical and demographic characteristics were measured when the participants were enrolled in the study.

Patients were enrolled at the time of disease investigation at the Sarcoidosis Centrum, Karolinska University Hospital Solna, Sweden. The diagnosis was established on radiographic manifestations, findings at bronchoscopy with bronchoalveolar lavage (BAL), including an elevated CD4/CD8 ratio > 3.5, and positive biopsies and in accordance with the criteria outlined by the World Association of Sarcoidosis and Other Granulomatous Disorders [WASOG, Statement on sarcoidosis (27)]. Sarcoidosis patients were further characterized into two clinical phenotypes, LS and non-LS. LS was defined by typical clinical manifestations with an acute onset of the disease, including fever, bilateral hilar lymphadenopathy on chest X-ray, bilateral ankle arthritis, and/or erythema nodosum. Non-LS was defined as a heterogeneous disease course with an insidious onset, unrelenting disease course, and a high risk for clinical organ impairment, such as developing fibrosis in the lung. Healthy controls included 3,085 individuals who were recruited from two large-scale epidemiological cohorts. Mainly, 2,025 individuals were from the Environmental Investigation of Rheumatoid arthritis (EIRA) study described in Klareskog et al. (25) and 1,060 individuals were from the Epidemiological investigation of risk factors for Multiple Sclerosis (EIMS) study described in Hedstrom et al. (26).

The German cohort

The case-control study consisted of 4,975 individuals, of which 413 were non-LS (27). The description and inclusion of these patients are described elsewhere (28, 29). Briefly, German control subjects (n = 4,498) were derived from Popgen (n = 2,485) (30), and the Heinz Nixdorf RECALL (HNR) study (n = 1,499) (31). Additionally, 304 control individuals of South German origin were recruited from the Bavarian Red Cross, and 210 control individuals were recruited from the Charité - Universitätsmedizin Berlin. The mean age and male percentage were 62.5 years (SD = 12.3) and 40% male in non-LS and 57.8 years (SD = 12.2), and 51% male in controls.

The USA African American cohort

The case-control study consisted of 1,657 individuals. The sample included 781 sarcoidosis cases characterized as non-LS and 876 healthy controls. All individuals were taken from an extensive cohort of African-American (AA) sarcoidosis patients and controls assembled from various studies. Further details are available in Refs. (3235).

Genotyping and quality control

Genotyping for sarcoidosis patients in the Swedish cohort was performed using the Illumina Immunochip platform and was performed at The SNP&SEQ Technology Platform in Uppsala University, Sweden, with Illumina Infinium assay using the Immuno Bead Chip (Immunochip version 1) as described in Rivera et al. (23). Healthy controls in the EIRA cohort (n = 2,086) were genotyped on the same platform at the Genome Institute of Singapore, as described previously (36, 37). Healthy controls in the EIMS cohort (n = 1,060) were genotyped with the same SNP array described elsewhere (26). Briefly, quality control filtering thresholds were applied using tools implemented in PLINK v1.09b software (38, 39). SNPs at call rate < 95% with minor allele frequency (MAF) > 1% and call rate < 99% with MAF ≤ 5% were filtered out. Moreover, SNPs that had Hardy–Weinberg Equilibrium (HWE) P < 1 × 10–7 (in the control group) were also excluded. Individuals with missing genotypes < 97% were also removed. Quality control (QC) resulted in 141,151 SNPs and 3,842 individuals on the Illumina Immunochip array.

Genotypes for the German cohort were quality controlled, as described in Fischer et al. (40). Genotyping for the African American cohort (US-AA) was performed at the OMRF using the Illumina HumanOmni1-Quad array for ∼1.1M variants across the genome. Of these, 121,988 were in common with the Immunochip platform and carried forward for analysis. Further details on genotyping and quality control filtering are described in Adrianto et al. (41).

Statistical analysis

Association analysis of LS and non-LS in sex groups

In each sex group, we examined the association of single nucleotide polymorphisms (SNPs) in LS and non-LS. For assessing the association in LS, we applied an additive model using logistic regression adjusted for age and four principal components (PCs). PCs were derived from principal component analysis (PCA) performed using EIGENSTRAT (42) software on a pruned genotyped data set. The pruning of genotypes was performed using the pruning function and default parameters implemented in PLINK 1.90 beta software (39). For assessing the association in non-LS, we applied an additive model using logistic regression without covariate adjustment due to the lack of age information across the cohorts. In LS and non-LS, a significance threshold was defined by genomic control-corrected P < 3.5 × 10–7 (0.05/141,151 quality-controlled SNPs) considering Bonferroni correction. A suggestive P < 5 × 10–5 was also considered to identify potential signals (4346). Furthermore, to account for high linkage disequilibrium (LD) (r2 ≥ 0.8) among associated SNPs, LD pruning was performed using PriorityPruner version 0.01.4 software1 with default parameters, which identified and prioritized SNPs in genomic loci and thus referred as tag-SNPs.

As an exploratory analysis, to assess the effect of sex on genetics, we implemented an interaction model where the disease outcome was regressed on SNP adjusted for sex and the interaction term (SNP × sex) as covariates. A significant threshold for the interaction term was set at P < 0.05. The interaction analysis was performed using the glm function implemented in PLINK 2.0 software (39).

Additionally, we also performed a heterogeneity test to assess the effects of association between males and females. Heterogeneity statistics i.e., Cochran’s Q-statistic, I2 statistic (the percentage of variability in the effect sizes), and tau2 (the between-study variance in our meta-analysis) were computed applying a random effects model and using the rma function implemented in the metafor R package (47).

Meta-analysis of significant SNPs

Meta-analysis was performed using SNP association results in non-LS. Herein, the meta-analysis was conducted using the inverse variance weighting (IVW) method implemented in the METAL software (48). For each SNP, the combined genetic effect size (defined by Beta), standard error (SE), meta-P-value (Pmeta), total variability in effect size, also known as the heterogeneity index (I2), and heterogeneity P-value were calculated. SNPs at Pmeta < 5 × 10–8 were defined as genome-wide significant. Two meta-analyses were conducted, (1) a meta-analysis considering SNP associations in the European ancestry cohorts (Sweden and Germany) and (2) a meta-analysis considering SNP associations in the multi-ethnic (Sweden, Germany, and US African American).

Besides METAL computing heterogeneity index (I2), we also performed an independent heterogeneity test applying the random effects model using the rma function implemented in the metafor R package (47). Summary of heterogeneity statistics include Cochran’s Q-statistic, I2 statistic (the percentage of variability in the effect sizes), and tau2 (the between-study variance in our meta-analysis).

Lookup of significant SNPs in the UK Biobank

Using the UKB resource, we examined SNP associations at P < 5 × 10–5 using summary statistics of a genome-wide association study (GWAS) on “doctor-diagnosed sarcoidosis” (Data-Field ID: 22133), consisting of 91,787 individuals (395 cases and 91,392 controls) in the UKBB (49). Chiefly, sex-stratified summary statistics GWAS on males (189 cases and 41,045 controls) and females (206 cases and 50,347 controls) were obtained from https://github.com/Nealelab/UK_Biobank_GWAS. GWAS analysis was conducted using an additive model and logistic regression adjusted by the first 20 PCs, age, and age2. A lookup significance threshold was set at P < 5 × 10–8.

Gene-based analysis

A gene-based analysis was conducted using SNP associations at P < 5 × 10–5 and MAGMA (50). Gene-based significance was set to P < 2.0 × 10–6 based on Bonferroni correction (0.05/24,769 genes). The number of genes was adopted from the VEGAS2 software definition, where 24,769 unique genes on the 22 autosomes were identified (51).

SNP enrichment analysis via expression quantitative trait loci (eQTLs) and pathway analysis

Expression quantitative trait loci enrichment was performed using SNP associations at P < 5 × 10–5 and gene expression data from 53 human tissues from the Genotype-Tissue Expression (GTEx) project version 8 available on the Functional Mapping and Annotation of genetic associations (FUMA) web tool (52). Additional eQTL enrichment using gene expression of lung tissue, whole blood, and immune cell types was conducted using EUGENE software (53, 54). A significant threshold for n SNP correlated with gene expression (defined as eQTL-SNP or eSNP) was set at P < 0.05.

Pathway analysis was conducted based on SNP associations at P < 5 × 10–8 to elucidate disease mechanisms using MetaCore software.

Results

The summary of phenotypes and the number of samples in the sex groups for cohorts investigated are shown in Table 1.

TABLE 1
www.frontiersin.org

Table 1. Sample size in association analysis in Sweden, Germany, and United States (A); sample size in replication lookup in the UK Biobank (B).

Association results

The genome-wide association analysis using high-density mapping Immunochip array was conducted in males and females independently in Swedish, German, and US-African American cohorts.

LS sex groups

In the Swedish cohort (Supplementary Tables 1, 2), 542 SNPs in LS males and 617 SNPs in LS females were identified at P < 3.5 × 10–7. Rs2187668 (OR = 4.715, SE = 0.1983, P = 6.02 × 10–14) located in HLA-DQA1 was the top signal identified in males, and rs9268219 (OR = 3.667, SE = 0.1586, P = 3.33 × 10–15) located in C6orf10 was the top signal identified in females. After considering linkage disequilibrium (LD) among associated signals, 77 SNPs in males and 101 SNPs in females remained, which were defined as tag-SNPs. The top results for LS sex groups (Table 2) are illustrated in the forest plot (Figure 1). Manhattan and QQ plots are shown in Supplementary Figure 1. Heterogeneity statistics are available in Supplementary Table 2X.

TABLE 2
www.frontiersin.org

Table 2. The top tag-SNPs (PGC < 3.5 × 10–7) associated with LS males and females in the Swedish cohort.

FIGURE 1
www.frontiersin.org

Figure 1. Forest plot for the 15 top common SNPs in the LS sex groups in the Swedish cohort.

Conditioning for the top SNP rs2187668 in males identified 12 SNPs at P < 5 × 10–5. The most significant SNP was rs9296079 (OR = 7.476, SE = 0.410, P = 9.51 × 10–7), located at 7.8 kb 5′ of HLA-DPB2. Similarly, conditioning for rs9268219 in females identified 8 SNPs at P < 5 × 10–5. The most significant was rs29243 (OR = 10.980, SE = 0.514, P = 6.04 × 10–6) located in GABBR1. Non-MHC SNPs at P < 5 × 10–5 were identified near LRRTM4 in males and CAST, LNPEP, TNIP1, CSMD1, B4GALNT1, SBNO2, and ABCG1 in females (Supplementary Table 3).

Non-LS sex groups

In the Swedish cohort (Supplementary Tables 4, 5), 5 SNPs in males at P < 3.5 × 10–7 and 1 SNP in females at P = 6.73 × 10–7 were identified. The top signals were rs1049550 (OR = 0.5231, SE = 0.1126, P = 1.82 × 10–8) located in ANXA11 in males and rs1964995 (OR = 0.5314, SE = 0.1198, P = 6.73 × 10–7) located at 36 kb of the 3′ of HLA-DRB5, in females. After LD assessment, 3 SNPs in males defined as tag-SNPs remained.

In the German cohort (Supplementary Tables 6, 7), rs1964995 (OR = 0.5452, SE = 0.08223, P = 3.40 × 10–11) in males and rs4502931 (OR = 0.5899, SE = 0.06549, P = 4.36 × 10–14) in females were identified as the top signals. In the US African American cohort (Supplementary Tables 8, 9), rs9271640 (OR = 1.771, SE = 0.1097, P = 3.62 × 10–7) in males and rs1964995 (OR = 0.5452, SE = 0.08223, P = 3.40 × 10–11) in females were identified as the top signals.

The top association results for each cohort are summarized in Table 3. Manhattan and QQ plots for non-LS sex groups are shown in Supplementary Figures 24.

Interaction analysis with sex

Results from the interaction analysis in the LS and non-LS sex groups in the Swedish cohort showed significant SNPs interacting with the sex variable. In LS sex groups (Supplementary Table 10), the most significant interacting SNPs were rs2853973 (Pint = 3.8 × 10–4) in males and rs1470410 (Pint = 5.44 × 10–4) in females. In non-LS sex groups (Supplementary Table 11), the most significant interacting SNPs were rs10940422 (Pint = 4.18 × 10–4) in males and rs12432418 (Pint = 4.85 × 10–5) in females.

Meta-analysis in non-LS

Non-LS sex groups

Meta-analysis in the European cohorts (Sweden and Germany) (Supplementary Tables 12, 13) identified 57 SNPs in males and 112 SNPs in females at Pmeta < 5 × 10–8. Top signal rs1964995 located 36 kb from the 3′ of HLA-DRB5 was the same in males (OR = 1.715, SE = 0.0656, Pmeta = 3.92 × 10–18) and females (OR = 1.7304, SE = 0.0598, Pmeta = 5 × 10–20). A forest plot illustrating the top findings is shown in Figure 2. Manhattan and Q-Q plots of the meta-analysis are shown in Supplementary Figure 5. Heterogeneity statistics in the SWE-GER are shown in Supplementary Table 13X.

FIGURE 2
www.frontiersin.org

Figure 2. Forest plots for the 15 top common meta-SNPs in the non-LS sex groups in the European group (Swedish-German) and the multi-ethnic group (Swedish-German-United States-African American).

TABLE 3
www.frontiersin.org

Table 3. The top tag-SNPs (PGC < 3.5 × 10–7) associated with non-LS males and females in Swedish, German, and USA-AA cohorts.

A comparison of associations at Pmeta < 5 × 10–8 showed 44 SNPs shared between non-LS males and females, 17 SNPs exclusively associated with non-LS males, and 215 SNPs exclusively associated with non-LS females. Most SNPs were in the MHC region. Non-MHC SNPs include rs694739 (OR = 1.3626, SE = 0.063, Pmeta = 9.21 × 10–07) located at 7.9 kbp from the 3′ of PRDX5 associated with non-LS males, and rs2573346 (OR = 0.7712, SE = 0.0564, Pmeta = 4.19 × 10–06) in ANXA11 associated with non-LS females.

Results from a meta-analysis in the multi-ethnic cohorts (Sweden, Germany, and US-African American) (Supplementary Tables 14, 15) identified 12 SNPs in non-LS males and 49 in non-LS females (Table 3). Top SNPs were rs1964995 (OR = 1.7736, SE = 0.0587, Pmeta = 1.54 × 10–22) located 36 kb from 3′ of HLA-DRB5 in males, and rs2395153 (Beta = 0.6565, SE = 0.0477, Pmeta = 1.11 × 10–18) located at 5.9 kb from 5′ of C6orf10 in females. The top meta-SNPs are shown in Table 4. Forest plots (Figure 2) illustrate the top common meta-SNPs in the European and multi-ancestry cohort groups. Manhattan plots of meta-analysis in multi-ethnic groups are shown in Supplementary Figure 6. Heterogeneity statistics in SWE-GER-USA are shown in Supplementary Table 15X.

TABLE 4
www.frontiersin.org

Table 4. The top significant SNPs (Pmeta < 5e-8) of GWAS meta-analysis on non-LS sex groups.

A comparison analysis at Pmeta < 5 × 10–8 showed a shared SNP between non-LS males and females, 11 SNPs exclusively associated with non-LS males, and 48 SNPs exclusively associated with non-LS females. Most associated SNPs were in the MHC region. Non-MHC SNPs were observed in non-LS females and included rs7813186, rs1049550, and rs7133604 (Data not shown).

Lookup of meta-SNPs in the UK Biobank

Genome-wide association study of “Doctor diagnosed sarcoidosis” in the UKB was used as a validation step. SNP-lookup showed several sex-specific SNPs in LS and non-LS at P < 5 × 10–5. Chiefly, 223 SNPs associated with LS males and 209 SNPs associated with LS females were validated at P < 5 × 10–8 (Supplementary Tables 16, 17). In non-LS, 11 non-LS males and 31 SNPs non-LS females in the SWE-GER cohorts were validated at P < 5 × 10–8 (Supplementary Tables 18, 19). In the multi-ethnic cohorts (SWE-GER-USA-AA), 2 SNPs in non-LS males and 8 SNPs in non-LS females were observed (Supplementary Tables 20, 21).

Gene-based analysis

The gene-based analysis revealed significant genomic loci associated with sex groups in LS and non-LS (Supplementary Tables 22, 23). The most significant genomic locus was HLA-DRA (P = 4.36 × 10–14) in LS males and FKBPL (P = 1.39 × 10–14) in LS females at a gene-based P < 2.0 × 10–6 (0.05/24,769 genes).

In the European cohorts (SWE-GER), 15 genes associated with non-LS males and 18 genes associated with non-LS females were identified, where the most significant genomic locus was HLA-DRA (P = 5.16 × 10–15) in non-LS males and HLA-DRB1 (P = 2.42 × 10–13) in non-LS females. In the multi-ethnic cohorts (SWE-GER-USA-AA), 5 genes associated with non-LS males and 8 genes associated with non-LS females were identified. The top genes were HLA-DRA (P = 1.88 × 10–14) in non-LS males and C6orf10 (P = 6.47 × 10–12) in non-LS females.

A comparison of gene-based analysis in LS and non-LS sex groups across cohorts (Supplementary Table 24) revealed differences in gene-based associations. In LS, gene-based associations showed 77 genes shared among sex groups, 11 genes exclusively associated with males, and 14 genes exclusively associated with females. In non-LS sex groups, in the European cohorts (SWE-GER), 18 genes were shared among non-LS sex groups, and 8 genes were exclusively associated with non-LS males. In the multi-ethnic cohorts (SWE-GER-USA-AA), 4 genes were shared, one gene was exclusively associated with non-LS males, and 4 were exclusively associated with non-LS females.

Figure 3 shows a circos plot depicting gene-based associations among LS and non-LS sex groups across cohorts.

FIGURE 3
www.frontiersin.org

Figure 3. Gene-based analysis in LS and non-LS sex groups across all cohorts at p < 2e-6. ls, LS in SWE; non-ls1, non-LS in SWE and GER; non-ls2, non-LS in SWE, GER, and USA-AA; shared, genes shared by sex groups; male, genes associated with males only; female, genes associated with females only.

Gene expression and expression quantitative trait loci (eQTL) enrichment

Results from tissue expression analysis using 54 tissue types in GText version 8 and FUMA revealed significant sex-specific gene expression at Bonferroni P-value (PBon < 0.05).

Löfgren’s syndrome males and females (Supplementary Figure 7), significant differential gene expression was observed in the spleen, small intestine, lung, and blood.

In the European cohorts (SWE, GER) (Supplementary Figure 8), differential gene expression was observed in the spleen, small intestine terminal Ileum, lung, and blood in non-LS males. In contrast, significant differential gene expression was observed in the brain hippocampus, and lung in non-LS females. In the multi-ethnic cohorts (SWE, GER, USA-AA) (Supplementary Figure 8), differentially expressed genes were observed in the spleen, small intestine terminal ilium, lung, and blood in non-LS males and females.

Expression quantitative trait loci enrichment using different eQTL databases across different tissues and cell types (Supplementary Table 25) revealed several significant eQTLs SNPs (eSNPs) in sex groups.

In LS sex groups, 311 eSNPs in 134 genes in males and 354 eSNPs in 140 genes in females were identified using the GText v8 database. In SWE-GER cohorts, non-LS sex groups, 32 eSNPs in 77 genes in males and 59 eSNPs in 89 genes in females were identified using the GText v8 database. In the SWE-GER-US-AA cohorts, non-LS sex groups, two eSNPs in males and 23 eSNPs in females were identified using the GText v8 database. Similar observations were also observed using eQTLGen and eQTLcatalogue databases.

Moreover, eQTL mapping using immune cell eQTLs (Supplementary Table 26) available in the EUGENE database identified significant eSNPs in relevant immune cell types. In LS sex, significant eSNPs of B cells, T cells, CD4, and CD8 T cells, granulocytes, monocytes, macrophages, neutrophils, and PBMCs were identified. In the SWE-GER-USA-AA non-LS sex groups, eQNPS of T cells, CD4 T cells, fibroblasts, monocytes, macrophages, and neutrophils were observed.

Gene enrichment and pathway analysis

In LS sex groups (Supplementary Table 27), a pathway map related to immune response induction of the antigen presentation machinery by IFN-gamma (FDR = 1.530e-8 in males; FDR = 1.926e-11 in females) were identified. In the SWE-GER non-LS sex groups (Supplementary Table 27), pathway maps related to maturation and migration of dendritic cells in skin sensitization (FDR = 2.085e-9) in males and related to immune response induction of the antigen presentation machinery by IFN-gamma (FDR = 1.998e-10) in females were identified. In the SWE-GER-USA-AA non-LS sex groups (Supplementary Table 27), pathway maps related to immune response lectin-induced complement pathway (FDR = 1.571e-5) in males and related to maturation and migration of dendritic cells in skin sensitization (FDR = 1.824e-6) in females were noted.

Discussion

The present study investigated sex differences in sarcoidosis, particularly in LS and non-LS phenotypes. This work is the first investigation for independently characterizing the genetic architecture of sarcoidosis in males and females. In brief, we analyzed genetic associations between LS and non-LS males and females conducted in European and multi-ethnic cohorts. The multi-ethnic cohorts comprised three population-based cohorts from Sweden, Germany, and the United States. The cohorts included different ancestries, i.e., white European and black African American.

For both LS and non-LS sex groups, a genetic analysis based on a genome-wide association study was conducted, followed by a meta-analysis, an SNP lookup in the UKB, gene-based analysis, differential gene expression analysis using the GText version 8 resource, eQTL tissue/cell mapping enrichment, and pathway analysis.

Löfgren’s syndrome, an acute form of sarcoidosis, is highly prevalent in Sweden (55). However, due to insufficient LS cases in the German cohort and being rare in the United States, genetic analyses in the LS sex groups were only conducted in the Swedish cohort. In contrast, genetic investigations in the non-LS sex groups were performed across all three cohorts.

Association analysis in LS sex groups identified distinct genetic associations clustering in the extended MHC (xMHC) (56), highlighting the role of genes located in this region as main component of LS susceptibility, as previously reported (23). The top LS findings were associations in HLA-DQA1 in males and C6orf10 in females. Although these genomic loci are on the same chromosome region (classical class II) and many variants in these genes are in linkage disequilibrium (LD r2 > 0.8), we observed significant differences in the genetic effects (beta estimates) among variants associated with LS males and females (Figure 1). Interestingly, after conditioning for the top MHC signals, we identified non-MHC loci worth investigating in future studies of LS. For instance, a nearby signal close to LRRTM4 associated with LS males is worth looking into, as LRRTM4 is a synapse-organizing molecule involved in protein-protein interactions that regulate the nervous and immune systems (57). Another interesting non-MHC association is LNPEP associated with LS females. LNPEP is an aminopeptidase that regulates hormones (i.e., arginine-vasopressin and oxytocin), and is involved in trimming peptide antigens for cross-presentation to T cells in autoimmune diseases (58).

In the European cohorts (SWE and GER), meta-analysis in non-LS sex groups identified mainly associations in the MHC classical class II subregion and a few non-MHC signals, chiefly in ANXA11, TMEM163, and nearby PRDX5 in both sex groups. As in LS, the differences among these signals were the genetic effects (beta estimates) in the sex groups, as illustrated in the forest plot (Figure 2). Similar findings were observed in the multi-ethnic cohorts (SWE, GER, and USA-AA).

To further understand the genetic differences, the comparison of genetic signals of LS and non-LS sex groups showed that some signals were associated with a specific sex, whereas other signals were common to both males and females; however, their genetic effects differed. These observations were then validated using the UKB resource, using the “doctors diagnosed sarcoidosis” phenotype. Based on this evidence, it is compelling to suggest that genetic differences in the sex groups could be linked to the observed differences in clinical manifestations in sarcoidosis male and female patients, as documented in various works (5961). Additionally, it is essential to highlight that ancestry also played a role in the differences in genetic effects in the sex groups, as the effects changed when multi-ethnic ancestries were considered.

Results from gene-based analysis showed patterns of genomic associations in the sex groups, such as shared and sex-specific (Figure 3). Notably, the genomic patterns of associations in the LS sex groups were more extensive since the LS susceptibility spanned across the xMHC, which harbors 252 genes and 139 pseudogenes. In the non-LS sex groups, on the contrary, the genomic patterns were limited to a few genes located in the MHC class II region. As expected, genomic patterns decreased as more ancestries were included, stressing the role of ancestry as a modifier of sarcoidosis susceptibility. Ancestry plays a significant role in the genetics of complex diseases (62, 63), which is also relevant in sarcoidosis.

Insights from tissue gene expression and eQTL enrichment showed that SNPs associated with biological sex play a role in gene expression variability across various tissues and cell types. For instance, gene expression in the spleen, small intestine, lung, and whole blood differed in LS and non-LS sex groups. Furthermore, SNP associations correlated with eQTLs of tissues and cell types, particularly immune cells eQTLs (e.g., B cells, T cells, CD4 T cells, CD8 T cells, macrophages, and monocytes), were identified. These findings underscore the functional role of genetic variants on gene expression of immune cells, thus shaping the immunopathogenesis of sarcoidosis in the sex groups. eQTL studies in complex diseases showed that eQTL SNPs act as master regulators of gene expression in several tissues and cell types (64, 65), influencing the expression of multiple genes and acting as gene regulators (66).

Expression quantitative trait loci findings in various immune cell types (Supplementary Table 26) further supports our hypothesis that sarcoidosis is an immune-regulated disease, as previously proposed (6769).

Indeed, differences in SNP associations in HLA genes (7075) may result in different immune responses and inflammatory disease course as seen in male and female patients with sarcoidosis. Thus, it is worth keeping in mind that HLA genes play a crucial role in autoimmunity and that their function has a profound sex bias on disease phenotype, as more females are affected by autoimmune (or chronic inflammatory) disorders than the counterparts. Understanding the molecular mechanisms underlying sex selection in the development of autoimmune disorders is an area of active research (76).

Compiling the evidence shown in this work, with a particular focus on the strong presence of HLA genetic associations in sarcoidosis, the involvement of various immune cells types, and the differences in immune response and clinical course in sex groups (67, 68, 77), the possibility of sarcoidosis being an autoimmune disorder cannot be ignored.

Besides HLA genes, we also identified signals in the ANXA11 gene in non-LS sex groups with different genetic effects, as denoted by rs1049550 (Figure 2). SNP associations in ANXA11 were first reported by Hofmann et al. (78) and have been linked to sarcoidosis risk, radiographic phenotype (Scadding stage IV), and to interact with HLA-DRA (79). SNPs in ANXA11 have also been reported to be associated with autoimmune diseases (80). Indeed, the functional role of ANXA11 in cell division and apoptosis (81) and the regulation of inflammatory cells (82) make ANXA11 an attractive biomarker for sarcoidosis.

Moreover, evidence for a pathway related to immune response induction of the antigen presentation machinery by IFN-gamma in LS sex groups offers an incentive to examine the role of IFN-gamma in sarcoidosis (Supplementary Table 27). IFN-gamma is a pro-inflammatory cytokine produced by immune cells and has immune effector functions on many genes involved in tissue homeostasis, immune/inflammatory responses, and tumor immunosurveillance (83). In the non-LS groups, pathway maps identified mainly were related to immune response but not specific to a sex group. Thus, more data is needed to understand the underlying mechanisms of non-LS as this sarcoidosis phenotype is heterogeneous and polygenic.

In summary, we provide new insights into a sex bias underpinning the genetic architecture of sarcoidosis. The differences observed in the genetic susceptibility in the sex groups in LS and non-LS add further information which may explain the well-recognized observations between sex groups, incidence, prevalence, age of onset, symptoms, severity of the disease, and drug reactions observed in male and female sarcoidosis patients (23, 8487).

It goes without a doubt that further genetic and omics studies ought to be conducted to broaden the characterization of the genetic architecture of sarcoidosis among sex groups. Such efforts, of course, shall consider the inclusion of multi-ancestry populations with large sample sizes and balanced sex ratios. An ongoing endeavor to meet such requirements is the MESARGEN consortium,2 an international framework aimed to include multiple ancestry populations and large sample sizes for genetics and genomics studies of sarcoidosis.

Limitations of the study

We used the Illumina Immunochip version 2 platform in our analysis—a customized design based on investigating major autoimmune and inflammatory diseases (88). Therefore, we limited our investigations to immune-mediated genomic regions, leaving out other potential loci implicated in sarcoidosis susceptibility. Despite this limitation, we offer substantial evidence that the genetic susceptibility of sarcoidosis and LS and non-LS phenotypes in males and females shared common genomic loci but also differed in their genetic architectures given the sex group. We observed these differences in gene-based analysis, eQTLs enrichment across different tissues, and pathway analysis. Further investigations with multiple ancestry populations and large sample sizes are needed to unveil discoveries.

Perspectives and significance

In this work, we revealed sex differences in the genetic architecture of sarcoidosis, particularly in clinical phenotypes LS and non-LS. Our work provides further knowledge that the genetic architectures of LS and non-LS are distinct, and that biological sex plays a role in determining the underlying genetic architecture.

Data availability statement

Access to genomic data is limited due to GPDR restrictions because the data contain personal or other sensitive information and cannot be deposited in public databases. All summary statistics from the analyses are accessible via DOI: 10.6084/m9.figshare.21785273.

Ethics statement

This study was approved by the Stockholm Ethical Review Board, Stockholm, Sweden. The patients/participants provided their written informed consent to participate in this study.

Author contributions

NR and LP designed and implemented the study. YX conducted all the genetic analyses in the Swedish cohort and a meta-analysis among all cohorts. YX and NR drafted the manuscript. SK, JG, and AE recruited patients, obtained samples, and clinically evaluated patients for phenotype characterization in the Swedish cohort. DE conducted genetic analyses in the German cohort. SS and JM-Q recruited patients, obtained samples, and clinically evaluated patients for phenotype characterization in the German cohort. BR and MI recruited patients and obtained samples. CM provided genotype data and oversaw analyses conducted by LG and NP of the US-AA cohort. All authors revised the manuscript and provided feedback.

Funding

This work was funded by the Swedish Heart-Lung Foundation awarded to NR (grant nos. 20170664, 20200505, and 20200506), SK (grant no. 20200163), and JG (grant no. 20190478); Karolinska Institutet Foundation awarded to NR (grant no. FS-2018:0007); and Swedish Research Council awarded to LP (grant no. 2018-02884) and JG (grant no. 2019-01034). This work was also funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) through grant FI 1935/1-1. The work received further infrastructure support from the DFG under Germany’s Excellence Strategy—EXC 2167-390884018 and the PopGen Biobank (Kiel, Germany; Ref Nr. 2018-032). The popgen 2.0 network (P2N) was supported by a grant from the German Federal Ministry of Education and Research (01EY1103). The US-AA sarcoidosis data and specimen collection are derived from the Trans-Omics in Precision Medicine (TOPMed) program was supported by the National Heart, Lung and Blood Institute (NHLBI). Genome Sequencing for “NHLBI TOPMed: African-American Sarcoidosis Genetics Resource” (phs001207.v3.p1) was performed at Baylor College of Medicine Human Genome Sequencing Center (3R01HL113326-04S1); Northwest Genomics Center (HHSN268201600032I); and Broad Institute Genomics Platform (HHSN268201600034I). Core support including centralized genomic read mapping and genotype calling, along with variant quality metrics and filtering were provided by the TOPMed Informatics Research Center (3R01HL-117626-02S1; contract HHSN268201800002I). Core support including phenotype harmonization, data management, sample-identity QC, and general program coordination were provided by the TOPMed Data Coordinating Center (R01HL-120393; U01HL-120393; contract HHSN268201800001I). We gratefully acknowledge the studies and participants who provided biological samples and data for TOPMed. The US-AA cohort investigators are funded via grants from the Foundation for Sarcoidosis Research (Chicago, IL) and the National Institutes of Health (R01HL113326, P30 GM110766, and U54GM104938-06). The computations and genomic data handling were enabled by resources provided by the Swedish National Infrastructure for Computing (SNIC) at the Uppsala Multidisciplinary Center for Advanced Computational Science (UPPMAX), partially funded by the Swedish Research Council through grant agreement no. 2018-05973.

Acknowledgments

We thank all patients and the personnel involved in this study.

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/fmed.2023.1132799/full#supplementary-material

Abbreviations

LS, Löfgren’s syndrome; Non-LS, non-Löfgren’s syndrome; SNPs, single nucleotide polymorphisms; eQTL, expression quantitative trait loci; UK Biobank, United Kingdom Biobank; ANXA11, Annexin A11 gene; GWAS, Genome-wide association study; US-AA, United States African Americans; PC, Principal component; PCA, principal component analysis; IVW, inverse variance weighting; OR, odds ratio; CI, confidence interval.

Footnotes

  1. ^ http://prioritypruner.sourceforge.net/
  2. ^ http://mesargen.org

References

1. Valeyre D, Prasse A, Nunes H, Uzunhan Y, Brillet P, Müller-Quernheim J. Sarcoidosis. Lancet. (2014) 383:1155–67. doi: 10.1016/S0140-6736(13)60680-7

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Sauer W, Stern B, Baughman R, Culver D, Royal W. High-risk sarcoidosis. current concepts and research imperatives. Ann Am Thorac Soc. (2017) 14(Suppl_6):S437–44. doi: 10.1513/AnnalsATS.201707-566OT

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Birnbaum A, Rifkin L. Sarcoidosis: sex-dependent variations in presentation and management. J Ophthalmol. (2014) 2014:236905. doi: 10.1155/2014/236905

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Thomeer M, Mateyo K. Sarcoidosis around the Globe. Semin Respir Crit Care Med. (2017) 38:393–403. doi: 10.1055/s-0037-1602845

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Ma X, Zhu L, Kurche J, Xiao H, Dai H, Wang C. Global and regional burden of interstitial lung disease and pulmonary sarcoidosis from 1990 to 2019: results from the Global Burden of Disease study 2019. Thorax. (2022) 77:596–605. doi: 10.1136/thoraxjnl-2020-216732

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Arkema E, Grunewald J, Kullberg S, Eklund A, Askling J. Sarcoidosis incidence and prevalence: a nationwide register-based assessment in Sweden. Eur Respir J. (2016) 48:1690–9. doi: 10.1183/13993003.00477-2016

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Sikjær M, Hilberg O, Ibsen R, Løkke A. Sarcoidosis: a nationwide registry-based study of incidence, prevalence and diagnostic work-up. Respir med. (2021) 187:106548. doi: 10.1016/j.rmed.2021.106548

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Ungprasert P, Crowson C, Matteson E. Influence of gender on epidemiology and clinical manifestations of sarcoidosis: a population-based retrospective cohort study 1976-2013. Lung. (2017) 195:87–91. doi: 10.1007/s00408-016-9952-6

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Singha A, Kirkland M, Drake W, Crouser ED. The influence of age and sex in sarcoidosis. Curr Opin Pulm Med. (2022) 28:307–13. doi: 10.1097/MCP.0000000000000882

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Lundkvist A, Kullberg S, Arkema E, Cedelund K, Eklund A, Grunewald J, et al. Differences in disease presentation between men and women with sarcoidosis: a cohort study. Respir Med. (2022) 191:106688. doi: 10.1016/j.rmed.2021.106688

PubMed Abstract | CrossRef Full Text | Google Scholar

11. James W, Koutroumpakis E, Saha B, Nathani A, Saavedra L, Yucel R, et al. Clinical features of extrapulmonary sarcoidosis without lung involvement. Chest. (2018) 154:349–56. doi: 10.1016/j.chest.2018.02.003

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Duvall C, Pavlovic N, Rosen N, Wand A, Griffin J, Okada D, et al. Sex differences in presentation and outcomes of cardiac sarcoidosis. J Cardiac Fail. (2022) 28:S16. doi: 10.1016/j.cardfail.2022.03.045

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Williamson K, Rosenbaum A, Kapa S, Blauwet L. Sex differences in cardiac sarcoidosis. J Am Coll Cardiol. (2020) 75(11_Supplement_1):999. doi: 10.1016/S0735-1097(20)31626-0

CrossRef Full Text | Google Scholar

14. Porter N, Beynon H, Randeva H. Endocrine and reproductive manifestations of sarcoidosis. QJM. (2003) 96:553–61. doi: 10.1093/qjmed/hcg103

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Singha A, Kirkland M, Drake W, Crouser ED. The influence of age and sex in sarcoidosis. Curr Opin Pulm Med. (2022) 28:307–13.

Google Scholar

16. Hadid V, Patenaude V, Oddy L, Abenhaim H. Sarcoidosis and pregnancy: obstetrical and neonatal outcomes in a population-based cohort of 7 million births. J Perinatal Med. (2015) 43:201–7. doi: 10.1515/jpm-2014-0017

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Cozier Y, Berman J, Palmer J, Boggs D, Wise L, Rosenberg L. Reproductive and hormonal factors in relation to incidence of sarcoidosis in US Black women: the black women’s health study. Am J Epidemiol. (2012) 176:635–41. doi: 10.1093/aje/kws145

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Grunewald J, Grutters J, Arkema E, Saketkoo L, Moller D, Muller-Quernheim J. Sarcoidosis. Nat Rev Dis Primers. (2019) 5:45. doi: 10.1038/s41572-019-0096-x

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Grunewald J, Eklund A. Löfgren’s syndrome. Am J Respir Crit Care Med. (2009) 179:307–12. doi: 10.1164/rccm.200807-1082OC

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Patterson K, Strek M. Pulmonary fibrosis in sarcoidosis. Clinical features and outcomes. Ann Am Thorac Soc. (2013) 10:362–70. doi: 10.1513/AnnalsATS.201303-069FR

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Nardi A, Brillet P, Letoumelin P, Girard F, Brauner M, Uzunhan Y, et al. Stage IV sarcoidosis: comparison of survival with the general population and causes of death. Eur Respir J. (2011) 38:1368–73. doi: 10.1183/09031936.00187410

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Jeny F, Uzunhan Y, Lacroix M, Gille T, Brillet P, Nardi A, et al. Predictors of mortality in fibrosing pulmonary sarcoidosis. Respir Med. (2020) 169:105997. doi: 10.1016/j.rmed.2020.105997

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Rivera N, Ronninger M, Shchetynsky K, Franke A, Nothen M, Muller-Quernheim J, et al. High-density genetic mapping identifies new susceptibility variants in sarcoidosis phenotypes and shows genomic-driven phenotypic differences. Am J Respir Crit Care Med. (2016) 193:1008–22. doi: 10.1164/rccm.201507-1372OC

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Rybicki B, Levin A, McKeigue P, Datta I, Gray-McGuire C, Colombo M, et al. A genome-wide admixture scan for ancestry-linked genes predisposing to sarcoidosis in African-Americans. Genes Immun. (2011) 12:67–77. doi: 10.1038/gene.2010.56

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Klareskog L, Stolt P, Lundberg K, Kallberg H, Bengtsson C, Grunewald J, et al. A new model for an etiology of rheumatoid arthritis: smoking may trigger HLA-DR (shared epitope)-restricted immune reactions to autoantigens modified by citrullination. Arthritis Rheum. (2006) 54:38–46. doi: 10.1002/art.21575

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Hedstrom A, Baarnhielm M, Olsson T, Alfredsson L. Tobacco smoking, but not Swedish snuff use, increases the risk of multiple sclerosis. Neurology. (2009) 73:696–701. doi: 10.1212/WNL.0b013e3181b59c40

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Statement on sarcoidosis. Joint Statement of the American Thoracic Society (ATS), the European Respiratory Society (ERS) and the World Association of Sarcoidosis and Other Granulomatous Disorders (WASOG) adopted by the ATS Board of Directors and by the ERS Executive Committee. Am J Respir Crit Care Med. (1999) 160:736–55.

PubMed Abstract | Google Scholar

28. Valentonyte R, Hampe J, Huse K, Rosenstiel P, Albrecht M, Stenzel A, et al. Sarcoidosis is associated with a truncating splice site mutation in BTNL2. Nat Genet. (2005) 37:357–64. doi: 10.1038/ng1519

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Hofmann S, Fischer A, Nothnagel M, Jacobs G, Schmid B, Wittig M, et al. Genome-wide association analysis reveals 12q13.3-q14.1 as new risk locus for sarcoidosis. Eur Respir J. (2013) 41:888–900. doi: 10.1183/09031936.00033812

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Krawczak M, Nikolaus S, von Eberstein H, Croucher P, El Mokhtari N, Schreiber S. PopGen: population-based recruitment of patients and controls for the analysis of complex genotype-phenotype relationships. Commun Genet. (2006) 9:55–61. doi: 10.1159/000090694

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Schmermund A, Mohlenkamp S, Stang A, Gronemeyer D, Seibel R, Hirche H, et al. Assessment of clinically silent atherosclerotic disease and established and novel risk factors for predicting myocardial infarction and cardiac death in healthy middle-aged subjects: rationale and design of the Heinz Nixdorf RECALL Study. Risk factors, evaluation of coronary calcium and lifestyle. Am Heart J. (2002) 144:212–8. doi: 10.1067/mhj.2002.123579

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Design of a case control etiologic study of sarcoidosis (ACCESS). ACCESS Research Group. J Clin Epidemiol. (1999) 52(12):1173–86.

Google Scholar

33. Rybicki B, Hirst K, Iyengar S, Barnard J, Judson M, Rose C, et al. A sarcoidosis genetic linkage consortium: the sarcoidosis genetic analysis (SAGA) study. Sarcoidosis Vasc Diffuse Lung Dis. (2005) 22:115–22.

PubMed Abstract | Google Scholar

34. Iannuzzi M, Maliarik M, Poisson L, Rybicki B. Sarcoidosis susceptibility and resistance HLA-DQB1 alleles in African Americans. Am J Respir Crit Care Med. (2003) 167:1225–31. doi: 10.1164/rccm.200209-1097OC

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Rasmussen A, Sevier S, Kelly J, Glenn S, Aberle T, Cooney C, et al. The lupus family registry and repository. Rheumatology. (2011) 50:47–59. doi: 10.1093/rheumatology/keq302

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Eyre S, Bowes J, Diogo D, Lee A, Barton A, Martin P, et al. High-density genetic mapping identifies new susceptibility loci for rheumatoid arthritis. Nat Genet. (2012) 44:1336–40. doi: 10.1038/ng.2462

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Padyukov L, Seielstad M, Ong R, Ding B, Ronnelid J, Seddighzadeh M, et al. A genome-wide association study suggests contrasting associations in ACPA-positive versus ACPA-negative rheumatoid arthritis. Ann Rheum Dis. (2011) 70:259–65. doi: 10.1136/ard.2009.126821

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira M, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. (2007) 81:559–75. doi: 10.1086/519795

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Chang C, Chow C, Tellier L, Vattikuti S, Purcell S, Lee J. Second-generation PLINK: rising to the challenge of larger and richer datasets. GigaScience. (2015) 4:7. doi: 10.1186/s13742-015-0047-8

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Fischer A, Ellinghaus D, Nutsua M, Hofmann S, Montgomery C, Iannuzzi M, et al. Identification of immune-relevant factors conferring sarcoidosis genetic risk. Am J Respir Crit Care Med. (2015) 192:727–36. doi: 10.1164/rccm.201503-0418OC

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Adrianto I, Lin C, Hale J, Levin A, Datta I, Parker R, et al. Genome-wide association study of African and European Americans implicates multiple shared and ethnic specific loci in sarcoidosis susceptibility. PLoS One. (2012) 7:e43907. doi: 10.1371/journal.pone.0043907

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Price A, Patterson N, Plenge R, Weinblatt M, Shadick N, Reich D. Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet. (2006) 38:904–9. doi: 10.1038/ng1847

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Steve E, John B, Dorothée D, Annette L, Anne B, Paul M, et al. High-density genetic mapping identifies new susceptibility loci for rheumatoid arthritis. Nat Genet. (2012) 44:1336.

Google Scholar

44. Beecham A, Patsopoulos N, Xifara D, Davis M, Kemppinen A, Cotsapas C, et al. Analysis of immune-related loci identifies 48 new susceptibility variants for multiple sclerosis. Nat Genet. (2013) 45:1353. doi: 10.1038/ng.2770

PubMed Abstract | CrossRef Full Text | Google Scholar

45. David E, Hansjörg B, Jorge E, Elke R, Anja M, Ingo M, et al. High-density genotyping study identifies four new susceptibility loci for atopic dermatitis. Nat Genet. (2013) 45:808.

PubMed Abstract | Google Scholar

46. Carmona F, Mackie Sarah l, Martín J, Taylor John C, Vaglio A, Eyre S, et al. A large-scale genetic analysis reveals a strong contribution of the HLA class II region to giant cell arteritis susceptibility. Am J Hum Genet. (2015) 96:565–80.

PubMed Abstract | Google Scholar

47. Viechtbauer W. Conducting meta-analyses in R with the metafor Package. J Stat Softw. (2010) 36:1–48. doi: 10.18637/jss.v036.i03

CrossRef Full Text | Google Scholar

48. Willer C, Li Y, Abecasis G. METAL: fast and efficient meta-analysis of genomewide association scans. Bioinformatics. (2010) 26:2190–1. doi: 10.1093/bioinformatics/btq340

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Sudlow C, Gallacher J, Allen N, Beral V, Burton P, Danesh J, et al. UK biobank: an open access resource for identifying the causes of a wide range of complex diseases of middle and old age. PLoS Med. (2015) 12:e1001779. doi: 10.1371/journal.pmed.1001779

PubMed Abstract | CrossRef Full Text | Google Scholar

50. De Leeuw CA, Mooij JM, Heskes T, Posthuma D. MAGMA: Generalized gene-set analysis of GWAS data. PLoS Comput Biol. (2015) 11:e1004219.

Google Scholar

51. Mishra A, Macgregor S. VEGAS2: software for more flexible gene-based testing. Twin Res Hum Genet. (2015) 18:86–91. doi: 10.1017/thg.2014.79

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Watanabe K, Taskesen E, Van Bochoven A, Posthuma D. Functional mapping and annotation of genetic associations with FUMA. Nat Commun. (2017) 8:1826. doi: 10.1038/s41467-017-01261-5

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Ferreira M, Jansen R, Willemsen G, Penninx B, Bain L, Vicente C, et al. Gene-based analysis of regulatory variants identifies 4 putative novel asthma risk genes related to nucleotide synthesis and signaling. J Allergy Clin Immunol. (2017) 139:1148–57. doi: 10.1016/j.jaci.2016.07.017

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Franke L, Jansen R. eQTL analysis in humans. Methods Mol Biol. (2009) 573:311. doi: 10.1007/978-1-60761-247-6_17

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Brown F, Modi P, Tanner LS. Lofgren syndrome. StatPearls. Treasure Island (FL): StatPearls Publishing Copyright © 2022, StatPearls Publishing LLC. (2022).

Google Scholar

56. Horton R, Wilming L, Rand V, Lovering R, Bruford E, Khodiyar V, et al. Gene map of the extended human MHC. Nat Rev Genet. (2004) 5:889–99. doi: 10.1038/nrg1489

PubMed Abstract | CrossRef Full Text | Google Scholar

57. Sinha R, Siddiqui T, Padmanabhan N, Wallin J, Zhang C, Karimi B, et al. LRRTM4: a novel regulator of presynaptic inhibition and ribbon synapse arrangements of retinal bipolar cells. Neuron. (2020) 105:1007–17.e5. doi: 10.1016/j.neuron.2019.12.028

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Paladini F, Fiorillo M, Tedeschi V, Mattorre B, Sorrentino R. The multifaceted nature of aminopeptidases ERAP1, ERAP2, and LNPEP: from evolution to disease. Front Immunol. (2020) 11:1576. doi: 10.3389/fimmu.2020.01576

PubMed Abstract | CrossRef Full Text | Google Scholar

59. Gerke A, Judson M, Cozier Y, Culver D, Koth L. Disease burden and variability in sarcoidosis. Ann Am Thoracic Soc. (2017) 14(Suppl_6):S421–8. doi: 10.1513/AnnalsATS.201707-564OT

PubMed Abstract | CrossRef Full Text | Google Scholar

60. Pandit P, Perez R, Roman J. Sex-based differences in interstitial lung disease. Am J Med Sci. (2020) 360:467–73. doi: 10.1016/j.amjms.2020.04.023

PubMed Abstract | CrossRef Full Text | Google Scholar

61. Sharp M, Psoter K, Balasubramanian A, Pulapaka A, Chen E, Brown S, et al. Heterogeneity of lung function phenotypes in sarcoidosis: role of race and sex differences. Ann Am Thoracic Soc. (2023) 20:30–7. doi: 10.1513/AnnalsATS.202204-328OC

PubMed Abstract | CrossRef Full Text | Google Scholar

62. Qiu J, Moore J, Darabos C. Studying the genetics of complex disease with ancestry-specific human phenotype networks: the case of type 2 diabetes in East Asian populations. Genet Epidemiol. (2016) 40:293–303. doi: 10.1002/gepi.21964

PubMed Abstract | CrossRef Full Text | Google Scholar

63. Dries D. Genetic ancestry, population admixture, and the genetic epidemiology of complex disease. Circulation. (2009) 2:540–3. doi: 10.1161/CIRCGENETICS.109.922898

PubMed Abstract | CrossRef Full Text | Google Scholar

64. Fairfax B, Makino S, Radhakrishnan J, Plant K, Leslie S, Dilthey A, et al. Genetics of gene expression in primary immune cells identifies cell type–specific master regulators and roles of HLA alleles. Nat Genet. (2012) 44:502–10. doi: 10.1038/ng.2205

PubMed Abstract | CrossRef Full Text | Google Scholar

65. Võsa U, Claringbould A, Westra H, Bonder M, Deelen P, Zeng B, et al. Large-scale cis- and trans-eQTL analyses identify thousands of genetic loci and polygenic scores that regulate blood gene expression. Nat Genet. (2021) 53:1300–10.

PubMed Abstract | Google Scholar

66. Tian J, Keller M, Broman A, Kendziorski C, Yandell B, Attie A, et al. The dissection of expression quantitative trait locus hotspots. Genetics. (2016) 202:1563–74. doi: 10.1534/genetics.115.183624

PubMed Abstract | CrossRef Full Text | Google Scholar

67. Kaiser Y, Eklund A, Grunewald J. Moving target: shifting the focus to pulmonary sarcoidosis as an autoimmune spectrum disorder. Eur Respir J. (2019) 54.

Google Scholar

68. Starshinova A, Malkova A, Basantsova N, Zinchenko Y, Kudryavtsev I, Ershov G, et al. Sarcoidosis as an autoimmune disease. Front Immunol. (2020) 10:2933. doi: 10.3389/fimmu.2019.02933

PubMed Abstract | CrossRef Full Text | Google Scholar

69. Song M, Manansala M, Parmar P, Ascoli C, Rubinstein I, Sweiss N. Sarcoidosis and autoimmunity. Curr Opin Pulm Med. (2021) 27:448–54. doi: 10.1097/MCP.0000000000000809

PubMed Abstract | CrossRef Full Text | Google Scholar

70. Rossman MD, Kreider ME. Is chronic beryllium disease sarcoidosis of known etiology? Sarcoidosis Vasc Diffuse Lung Dis. (2003) 20(2):104–9.

Google Scholar

71. Sharma SK, Balamurugan A, Pandey RM, Saha PK, Mehra NK. Human leukocyte antigen-DR alleles influence the clinical course of pulmonary sarcoidosis in Asian Indians. Am J Respir Cell Mol Biol. (2003) 29(2):225–31.

Google Scholar

72. Darlington P, Tallstedt L, Padyukov L, Kockum I, Cederlund K, Eklund A, et al. HLA-DRB1*alleles and symptoms associated with Heerfordt’s syndrome in sarcoidosis. Eur Respir J. (2011) 38:1151–57.

Google Scholar

73. Zhou Y, Shen L, Zhang Y, Jiang D, Li H. Human leukocyte antigen-A, -B, and -DRB1 alleles and sarcoidosis in Chinese Han subjects. Hum Immunol. (2011) 72:571–5.

Google Scholar

74. Wennerstrom A, Pietinalho A, Vauhkonen H, Lahtela L, Palikhe A, Hedman J, et al. HLA-DRB1 allele frequencies and C4 copy number variation in Finnish sarcoidosis patients and associations with disease prognosis. Hum Immunol. (2012) 73:93–100.

Google Scholar

75. van Moorsel CHM, Petrek M, Rivera NV Unravelling the genetic basis of sarcoidosis. Eur Respir Monogr. (2022) 96:41–56.

Google Scholar

76. Rubtsova K, Marrack P, Rubtsov A. Sexual dimorphism in autoimmunity. J Clin Investig. (2015) 125:2187–93. doi: 10.1172/JCI78082

PubMed Abstract | CrossRef Full Text | Google Scholar

77. Korsten P, Tampe B, Konig MF, Nikiphorou E. Sarcoidosis and autoimmune diseases: differences, similarities and overlaps. Curr Opin Pulm Med. (2018) 24(5):504–12.

Google Scholar

78. Hofmann S, Franke A, Fischer A, Jacobs G, Nothnagel M, Gaede KI. Genome-wide association study identifies ANXA11 as a new susceptibility locus for sarcoidosis. Nat Genet. (2008) 40:1103–6. doi: 10.1038/ng.198

PubMed Abstract | CrossRef Full Text | Google Scholar

79. Levin AM, Iannuzzi MC, Montgomery CG, Trudeau S, Datta I, McKeigue P, et al. Association of ANXA11 genetic variation with sarcoidosis in African Americans and European Americans. Genes Immun. (2013) 14:13–8.

Google Scholar

80. Iaccarino L, Ghirardello A, Canova M, Zen M, Bettio S, Nalotto L, et al. Anti-annexins autoantibodies: their role as biomarkers of autoimmune diseases. Autoimmun Rev. (2011) 10:553–8. doi: 10.1016/j.autrev.2011.04.007

PubMed Abstract | CrossRef Full Text | Google Scholar

81. Moss SE, Morgan RO. The annexins. Genome Biol. (2004) 5:219. doi: 10.1186/gb-2004-5-4-219

PubMed Abstract | CrossRef Full Text | Google Scholar

82. Bruschi M, Petretto A, Vaglio A, Santucci L, Candiano G, Ghiggeri G. Annexin A1 and autoimmunity: from basic science to clinical applications. Int J Mol Sci. (2018) 19:1348.

Google Scholar

83. Castro F, Cardoso A, Gonçalves R, Serre K, Oliveira M. Interferon-gamma at the crossroads of tumor immune surveillance or evasion. Front Immunol. (2018) 9:847. doi: 10.3389/fimmu.2018.00847

PubMed Abstract | CrossRef Full Text | Google Scholar

84. Dempsey O, Paterson E, Kerr K, Denison A. Sarcoidosis. BMJ. (2009) 339:b3206. doi: 10.1136/bmj.b3206

PubMed Abstract | CrossRef Full Text | Google Scholar

85. Miedema J, Kaiser Y, Broos C, Wijsenbeek M, Grunewald J, Kool M. Th17-lineage cells in pulmonary sarcoidosis and Lofgren’s syndrome: friend or foe? J Autoimmun. (2018) 87:82–96. doi: 10.1016/j.jaut.2017.12.012

PubMed Abstract | CrossRef Full Text | Google Scholar

86. Grunewald J, Eklund A. Sex-specific manifestations of Lofgren’s syndrome. Am J Respir Crit Care Med. (2007) 175:40–4. doi: 10.1164/rccm.200608-1197OC

PubMed Abstract | CrossRef Full Text | Google Scholar

87. Fischer A, Valentonyte R, Nebel A, Nothnagel M, Muller-Quernheim J, Schurmann M, et al. Female-specific association of C-C chemokine receptor 5 gene polymorphisms with Lofgren’s syndrome. J Mol Med. (2008) 86:553–61. doi: 10.1007/s00109-008-0315-5

PubMed Abstract | CrossRef Full Text | Google Scholar

88. Cortes A, Brown M. Promise and pitfalls of the immunochip. Arthritis Res Ther. (2011) 13:101. doi: 10.1186/ar3204

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: sarcoidosis, genetics, genome-wide association study (GWAS), meta-analysis, single nucleotide polymorphisms, immunogenetics and HLA

Citation: Xiong Y, Kullberg S, Garman L, Pezant N, Ellinghaus D, Vasila V, Eklund A, Rybicki BA, Iannuzzi MC, Schreiber S, Müller-Quernheim J, Montgomery CG, Grunewald J, Padyukov L and Rivera NV (2023) Sex differences in the genetics of sarcoidosis across European and African ancestry populations. Front. Med. 10:1132799. doi: 10.3389/fmed.2023.1132799

Received: 27 December 2022; Accepted: 10 April 2023;
Published: 11 May 2023.

Edited by:

Ying Zhou, Tongji University, China

Reviewed by:

Humberto Garcia-Ortiz, National Institute of Genomic Medicine (INMEGEN), Mexico
Ilias C. Papanikolaou, General Hospital of Corfu, Greece

Copyright © 2023 Xiong, Kullberg, Garman, Pezant, Ellinghaus, Vasila, Eklund, Rybicki, Iannuzzi, Schreiber, Müller-Quernheim, Montgomery, Grunewald, Padyukov and Rivera. 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: Natalia V. Rivera, natalia.rivera@ki.se

Disclaimer: 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.