Sex-specific effects of injury and beta-adrenergic activation on metabolic and inflammatory mediators in a murine model of post-traumatic osteoarthritis.

OBJECTIVE: Metabolic processes are intricately linked to the resolution of innate inflammation and tissue repair, two critical steps for treating post-traumatic osteoarthritis (PTOA). Based on lipolytic and immunoregulatory actions of norepinephrine, we hypothesized that intra-articular ß-adrenergic receptor (ßAR) stimulation would suppress PTOA-associated inflammation in the infrapatellar fat pad (IFP) and synovium. DESIGN: We used the ßAR agonist isoproterenol to perturb intra-articular metabolism 3.5 weeks after applying a non-invasive single-load compression injury to knees of 12-week-old male and female mice. We examined the acute effects of intra-articular isoproterenol treatment relative to saline on IFP histology, multiplex gene expression of synovium-IFP tissue, synovial fluid metabolomics, and mechanical allodynia. RESULTS: Injured knees developed PTOA pathology characterized by heterotopic ossification, articular cartilage loss, and IFP atrophy and fibrosis. Isoproterenol suppressed the upregulation of pro-fibrotic genes and downregulated the expression of adipose genes and pro-inflammatory genes (Adam17, Cd14, Icam1, Csf1r, and Casp1) in injured joints of female (but not male) mice. Analysis of published single-cell RNA-seq data identified elevated catecholamine-associated gene expression in resident-like synovial-IFP macrophages after injury. Injury substantially altered synovial fluid metabolites by increasing amino acids, peptides, sphingolipids, phospholipids, bile acids, and dicarboxylic acids, but these changes were not appreciably altered by isoproterenol. Intra-articular injection of either isoproterenol or saline increased mechanical allodynia in female mice, whereas neither substance affected male mice. CONCLUSIONS: Acute ßAR activation altered synovial-IFP transcription in a sex and injury-dependent manner, suggesting that women with PTOA may be more sensitive than men to treatments targeting sympathetic neural signaling pathways.

Novel HLA associations with outcomes of Mycobacterium tuberculosis exposure and sarcoidosis in individuals of African ancestry using nearest-neighbor feature selection.

Tuberculosis and sarcoidosis are inflammatory diseases characterized by granulomas that may occur in any organ but are often found in the lung. The panoply of classical human leukocyte antigen (HLA) alleles associated with occurrence and/or severity of both diseases varies considerably across studies. This heterogeneity of results, due to variation in factors like ancestry and disease subphenotype, as well as the use of simple modeling strategies to elucidate likely complex relationships, has made conclusions about underlying commonalities difficult. Here we perform HLA association analyses in individuals of African ancestry, using a greater resolution to include subphenotypes of disease and employing more comprehensive analytical techniques. Using a novel application of nearest-neighbor feature selection to score allelic importance, we investigated HLA allele association with Mycobacterium tuberculosis exposure outcomes in the first analysis of both latent Mycobacterium tuberculosis infection and active disease compared with those who, despite long-term exposure to active index cases, have neither positive diagnostic tests nor display clinical symptoms. We also compared persistent to resolved sarcoidosis. This led to the identification of novel HLA associations and evidence of main effects and interaction effects. We found strikingly similar main effects and interaction effects at HLA-DRB1, -DQB1, and -DPB1 in those resistant to tuberculosis (either latent or active) and persistent sarcoidosis.

Multiple Correspondence Analysis and HLA-Associations of Organ Involvement in a Large Cohort of African-American and European-American Patients with Sarcoidosis.

Sarcoidosis is a systemic granulomatous disease with predominant pulmonary involvement and vast heterogeneity of clinical manifestations and disease outcomes. African American (AA) patients suffer greater morbidity and mortality. Using Multiple Correspondence Analysis, we identified seven clusters of organ involvement in European American (EA; n = 385) patients which were similar to those previously described in a Pan-European (GenPhenReSa) and a Spanish cohort (SARCOGEAS). In contrast, AA (n = 987) had six, less well-defined and overlapping clusters with little similarity to the cluster identified in the EA cohort evaluated at the same U.S. institutions. Association of cluster membership with two-digit HLA-DRB1 alleles demonstrated ancestry-specific patterns of association and replicated known HLA effects.These results further support the notion that genetically influenced immune risk profiles, which differ based on ancestry, play a role in phenotypic heterogeneity. Dissecting such risk profiles will move us closer to personalized medicine for this complex disease.

Extended methods for gene-environment-wide interaction scans in studies of admixed individuals with varying degrees of relationships.

The etiology of many complex diseases involves both environmental exposures and inherited genetic predisposition as well as interactions between them. Gene-environment-wide interaction studies (GEWIS) provide a means to identify the interactions between genetic variation and environmental exposures that underlie disease risk. However, current GEWIS methods lack the capability to adjust for the potentially complex correlations in studies with varying degrees of relationships (both known and unknown) among individuals in admixed populations. We developed novel generalized estimating equation (GEE) based methods-GEE-adaptive and GEE-joint-to account for phenotypic correlations due to kinship while accounting for covariates, including, measures of genome-wide ancestry. In simulation studies of admixed individuals, both methods controlled family-wise error rates, an advantage over the case-only approach. They demonstrated higher power than traditional case-control methods across a wide range of underlying alternative hypotheses, especially where both marginal and interaction effects were present. We applied the proposed method to conduct a GEWIS of a known sarcoidosis risk factor (insecticide exposure) and risk of sarcoidosis in African Americans and identified two novel loci with suggestive evidence of G × E interaction.

Recent advances in sarcoidosis genomics: epigenetics, gene expression, and gene by environment (G × E) interaction studies.

PURPOSE OF REVIEW: We aim to review the most recent findings in genomics of sarcoidosis and highlight the gaps in the field. RECENT FINDINGS: Original explorations of sarcoidosis subphenotypes, including cases associated with the World Trade Center and ocular sarcoidosis, have identified novel risk loci. Innovative gene--environment interaction studies utilizing modern analytical techniques have discovered risk loci associated with smoking and insecticide exposure. The application of whole-exome sequencing has identified genetic variants associated with persistent sarcoidosis and rare functional variations. A single epigenomics study has provided background knowledge of DNA methylation mechanisms in comparison with gene expression data. The application of machine-learning techniques has suggested new drug repositioning for the treatment of sarcoidosis. Several gene expression studies have identified prominent inflammatory pathways enriched in the affected tissue. SUMMARY: Certainly, sarcoidosis research has recently advanced in the exploration of disease subphenotypes, utilizing novel analytical techniques, and including measures of clinical variation. Nevertheless, large-scale and diverse cohorts investigated with advanced sequencing methods, such as whole-genome and single-cell RNA sequencing, epigenomics, and meta-analysis coupled with cutting-edge analytic approaches, when employed, will broaden and translate genomics findings into clinical applications, and ultimately open venues for personalized medicine.

ARID3a gene profiles are strongly associated with human interferon alpha production.

Type I interferons (IFN) causes inflammatory responses to pathogens, and can be elevated in autoimmune diseases such as systemic lupus erythematosus (SLE). We previously reported unexpected associations of increased numbers of B lymphocytes expressing the DNA-binding protein ARID3a with both IFN alpha (IFNα) expression and increased disease activity in SLE. Here, we determined that IFNα producing low density neutrophils (LDNs) and plasmacytoid dendritic cells (pDCs) from SLE patients exhibit strong associations between ARID3a protein expression and IFNα production. Moreover, SLE disease activity indices correlate most strongly with percentages of ARID3a+ LDNs, but were also associated, less significantly, with IFNα expression in LDNs and pDCs. Hierarchical clustering and transcriptome analyses of LDNs and pDCs revealed SLE patients with low ARID3a expression cluster with healthy controls and identified gene profiles associated with increased proportions of ARID3a- and IFNα-expressing cells of each type. These data identify ARID3a as a potential transcription regulator of IFNα-related inflammatory responses and other pathways important for SLE disease activity.

Downregulation of E Protein Activity Augments an ILC2 Differentiation Program in the Thymus.

Innate lymphoid cells (ILCs) are important regulators in various immune responses. The current paradigm states that all newly made ILCs originate from common lymphoid progenitors in the bone marrow. Id2, an inhibitor of E protein transcription factors, is indispensable for ILC differentiation. Unexpectedly, we found that ectopically expressing Id1 or deleting two E protein genes in the thymus drastically increased ILC2 counts in the thymus and other organs where ILC2 normally reside. Further evidence suggests a thymic origin of these mutant ILC2s. The mutant mice exhibit augmented spontaneous infiltration of eosinophils and heightened responses to papain in the lung and increased ability to expulse the helminth parasite, Nippostrongylus brasiliensis These results prompt the questions of whether the thymus naturally has the capacity to produce ILC2s and whether E proteins restrain such a potential. The abundance of ILC2s in Id1 transgenic mice also offers a unique opportunity for testing the biological functions of ILC2s.

The IRF5-TNPO3 association with systemic lupus erythematosus has two components that other autoimmune disorders variably share.

Exploiting genotyping, DNA sequencing, imputation and trans-ancestral mapping, we used Bayesian and frequentist approaches to model the IRF5-TNPO3 locus association, now implicated in two immunotherapies and seven autoimmune diseases. Specifically, in systemic lupus erythematosus (SLE), we resolved separate associations in the IRF5 promoter (all ancestries) and with an extended European haplotype. We captured 3230 IRF5-TNPO3 high-quality, common variants across 5 ethnicities in 8395 SLE cases and 7367 controls. The genetic effect from the IRF5 promoter can be explained by any one of four variants in 5.7 kb (P-valuemeta = 6 × 10(-49); OR = 1.38-1.97). The second genetic effect spanned an 85.5-kb, 24-variant haplotype that included the genes IRF5 and TNPO3 (P-valuesEU = 10(-27)-10(-32), OR = 1.7-1.81). Many variants at the IRF5 locus with previously assigned biological function are not members of either final credible set of potential causal variants identified herein. In addition to the known biologically functional variants, we demonstrated that the risk allele of rs4728142, a variant in the promoter among the lowest frequentist probability and highest Bayesian posterior probability, was correlated with IRF5 expression and differentially binds the transcription factor ZBTB3. Our analytical strategy provides a novel framework for future studies aimed at dissecting etiological genetic effects. Finally, both SLE elements of the statistical model appear to operate in Sjögren's syndrome and systemic sclerosis whereas only the IRF5-TNPO3 gene-spanning haplotype is associated with primary biliary cirrhosis, demonstrating the nuance of similarity and difference in autoimmune disease risk mechanisms at IRF5-TNPO3.

Identification of Immune-Relevant Factors Conferring Sarcoidosis Genetic Risk.

RATIONALE: Genetic variation plays a significant role in the etiology of sarcoidosis. However, only a small fraction of its heritability has been explained so far. OBJECTIVES: To define further genetic risk loci for sarcoidosis, we used the Immunochip for a candidate gene association study of immune-associated loci. METHODS: Altogether the study population comprised over 19,000 individuals. In a two-stage design, 1,726 German sarcoidosis cases and 5,482 control subjects were genotyped for 128,705 single-nucleotide polymorphisms using the Illumina Immunochip for the screening step. The remaining 3,955 cases, 7,514 control subjects, and 684 parents of affected offspring were used for validation and replication of 44 candidate and two established risk single-nucleotide polymorphisms. MEASUREMENTS AND MAIN RESULTS: Four novel susceptibility loci were identified with genome-wide significance in the European case-control populations, located on chromosomes 12q24.12 (rs653178; ATXN2/SH2B3), 5q33.3 (rs4921492; IL12B), 4q24 (rs223498; MANBA/NFKB1), and 2q33.2 (rs6748088; FAM117B). We further defined three independent association signals in the HLA region with genome-wide significance, peaking in the BTNL2 promoter region (rs5007259), at HLA-B (rs4143332/HLA-B*0801) and at HLA-DPB1 (rs9277542), and found another novel independent signal near IL23R (rs12069782) on chromosome 1p31.3. CONCLUSIONS: Functional predictions and protein network analyses suggest a prominent role of the drug-targetable IL23/Th17 signaling pathway in the genetic etiology of sarcoidosis. Our findings reveal a substantial genetic overlap of sarcoidosis with diverse immune-mediated inflammatory disorders, which could be of relevance for the clinical application of modern therapeutics.

Association of HLA-DRB1 with Sarcoidosis Susceptibility and Progression in African Americans.

HLA-DRB1 is a sarcoidosis risk gene, and the *03:01 allele is strongly associated with disease resolution in European sarcoidosis cases. Whereas the HLA-DRB1 variation is associated with sarcoidosis susceptibility in African Americans, DRB1 risk alleles are not as well defined, and associations with disease resolution have not been studied. Associations between genotyped and imputed HLA-DRB1 alleles and disease susceptibility/resolution were evaluated in a sample of 1,277 African-American patients with sarcoidosis and 1,467 control subjects. In silico binding assays were performed to assess the functional significance of the associated alleles. Increased disease susceptibility was associated with the HLA-DRB1 alleles *12:01 (odds ratio [OR], 2.11; 95% confidence interval [CI], 1.65-2.69; P = 3.2 × 10(-9)) and *11:01 (OR, 1.69; 95% CI, 1.42-2.01; P = 3.0 × 10(-9)). The strongest protective association was found with *03:01 (OR, 0.56; 95% CI, 0.44-0.73; P = 1.0 × 10(-5)). The African-derived allele *03:02 was associated with decreased risk of persistent radiographic disease (OR, 0.52; 95% CI, 0.37-0.72; P = 1.3 × 10(-4)), a finding consistent across the three component studies comprising the analytic sample. The DRB1*03:01 association with disease persistence was dependent upon local ancestry, with carriers of at least one European allele at DRB1 at a decreased risk of persistent disease (OR, 0.36; 95% CI, 0.14-0.94; P = 0.037). Results of in silico binding analyses showed that DRB1*03:01 consistently demonstrated the highest binding affinities for six bacterial peptides previously found in sarcoidosis granulomas, whereas *12:01 displayed the lowest binding affinities. This study has identified DRB1*03:01 and *03:02 as novel alleles associated with disease susceptibility and course in African Americans. Further investigation of DRB1*03 alleles may uncover immunologic factors that favor sarcoidosis protection and resolution among African Americans.

Efficient generalized least squares method for mixed population and family-based samples in genome-wide association studies.

Genome-wide association studies (GWAS) that draw samples from multiple studies with a mixture of relationship structures are becoming more common. Analytical methods exist for using mixed-sample data, but few methods have been proposed for the analysis of genotype-by-environment (G×E) interactions. Using GWAS data from a study of sarcoidosis susceptibility genes in related and unrelated African Americans, we explored the current analytic options for genotype association testing in studies using both unrelated and family-based designs. We propose a novel method-generalized least squares (GLX)-to estimate both SNP and G×E interaction effects for categorical environmental covariates and compared this method to generalized estimating equations (GEE), logistic regression, the Cochran-Armitage trend test, and the WQLS and MQLS methods. We used simulation to demonstrate that the GLX method reduces type I error under a variety of pedigree structures. We also demonstrate its superior power to detect SNP effects while offering computational advantages and comparable power to detect G×E interactions versus GEE. Using this method, we found two novel SNPs that demonstrate a significant genome-wide interaction with insecticide exposure-rs10499003 and rs7745248, located in the intronic and 3' UTR regions of the FUT9 gene on chromosome 6q16.1.

Encore: Genetic Association Interaction Network centrality pipeline and application to SLE exome data.

Open source tools are needed to facilitate the construction, analysis, and visualization of gene-gene interaction networks for sequencing data. To address this need, we present Encore, an open source network analysis pipeline for genome-wide association studies and rare variant data. Encore constructs Genetic Association Interaction Networks or epistasis networks using two optional approaches: our previous information-theory method or a generalized linear model approach. Additionally, Encore includes multiple data filtering options, including Random Forest/Random Jungle for main effect enrichment and Evaporative Cooling and Relief-F filters for enrichment of interaction effects. Encore implements SNPrank network centrality for identifying susceptibility hubs (nodes containing a large amount of disease susceptibility information through the combination of multivariate main effects and multiple gene-gene interactions in the network), and it provides appropriate files for interactive visualization of a network using tools from our online Galaxy instance. We implemented these algorithms in C++ using OpenMP for shared-memory parallel analysis on a server or desktop. To demonstrate Encore's utility in analysis of genetic sequencing data, we present an analysis of exome resequencing data from healthy individuals and those with Systemic Lupus Erythematous (SLE). Our results verify the importance of the previously associated SLE genes HLA-DRB and NCF2, and these two genes had the highest gene-gene interaction degrees among the susceptibility hubs. An additional 14 genes previously associated with SLE emerged in our epistasis network model of the exome data, and three novel candidate genes, ST8SIA4, CMTM4, and C2CD4B, were implicated in the model. In summary, we present a comprehensive tool for epistasis network analysis and the first such analysis of exome data from a genetic study of SLE.

Comparison of the American-European Consensus Group Sjogren's syndrome classification criteria to newly proposed American College of Rheumatology criteria in a large, carefully characterised sicca cohort.

OBJECTIVE: To compare the performance of the American-European Consensus Group (AECG) and the newly proposed American College of Rheumatology (ACR) classification criteria for Sjögren's Syndrome (SS) in a well-characterised sicca cohort, given ongoing efforts to resolve discrepancies and weaknesses in the systems. METHODS: In a multidisciplinary clinic for the evaluation of sicca, we assessed features of salivary and lacrimal gland dysfunction and autoimmunity as defined by tests of both AECG and ACR criteria in 646 participants. Global gene expression profiles were compared in a subset of 180 participants. RESULTS: Application of the AECG and ACR criteria resulted in classification of 279 and 268 participants with SS, respectively. Both criteria were met by 244 participants (81%). In 26 of the 35 AECG+/ACR participants, the minor salivary gland biopsy focal score was ≥1 (74%), while nine had positive anti-Ro/La (26%). There were 24 AECG-/ACR+ who met ACR criteria mainly due to differences in the scoring of corneal staining. All patients with SS, regardless of classification, had similar gene expression profiles, which were distinct from the healthy controls. CONCLUSIONS: The two sets of classification criteria yield concordant results in the majority of cases and gene expression profiling suggests that patients meeting either set of criteria are more similar to other SS participants than to healthy controls. Thus, there is no clear evidence for increased value of the new ACR criteria over the old AECG criteria from the clinical or biological perspective. It is our contention, supported by this report, that improvements in diagnostic acumen will require a more fundamental understanding of the pathogenic mechanisms than is at present available.

Performance of HLA allele prediction methods in African Americans for class II genes HLA-DRB1, -DQB1, and -DPB1.

BACKGROUND: The expense of human leukocyte antigen (HLA) allele genotyping has motivated the development of imputation methods that use dense single nucleotide polymorphism (SNP) genotype data and the region's haplotype structure, but the performance of these methods in admixed populations (such as African Americans) has not been adequately evaluated. We compared genotype-based-derived from both genome-wide genotyping and targeted sequencing-imputation results to existing allele data for HLA-DRB1, -DQB1, and -DPB1. RESULTS: In European Americans, the newly-developed HLA Genotype Imputation with Attribute Bagging (HIBAG) method outperformed HLA*IMP:02. In African Americans, HLA*IMP:02 performed marginally better than HIBAG pre-built models, but HIBAG models constructed using a portion of our African American sample with both SNP genotyping and four-digit HLA class II allele typing had consistently higher accuracy than HLA*IMP:02. However, HIBAG was significantly less accurate in individuals heterozygous for local ancestry (p ≤0.04). Accuracy improved in models with equal numbers of African and European chromosomes. Variants added by targeted sequencing and SNP imputation further improved both imputation accuracy and the proportion of high quality calls. CONCLUSION: Combining the HIBAG approach with local ancestry and dense variant data can produce highly-accurate HLA class II allele imputation in African Americans.

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

[This corrects the article DOI: 10.3389/fmed.2023.1132799.].

Association of two independent functional risk haplotypes in TNIP1 with systemic lupus erythematosus.

OBJECTIVE: Systemic lupus erythematosus (SLE) is an autoimmune disease characterized by autoantibody production and altered type I interferon expression. Genetic surveys and genome-wide association studies have identified >30 SLE susceptibility genes. One of these genes, TNIP1, encodes the ABIN1 protein. ABIN1 functions in the immune system by restricting NF-κB signaling. The present study was undertaken to investigate the genetic factors that influence association with SLE in genes that regulate the NF-κB pathway. METHODS: We analyzed a dense set of genetic markers spanning TNIP1 and TAX1BP1, as well as the TNIP1 homolog TNIP2, in case-control populations of diverse ethnic origins. TNIP1, TNIP2, and TAX1BP1 were fine-mapped in a total of 8,372 SLE cases and 7,492 healthy controls from European-ancestry, African American, Hispanic, East Asian, and African American Gullah populations. Levels of TNIP1 messenger RNA (mRNA) and ABIN1 protein in Epstein-Barr virus-transformed human B cell lines were analyzed by quantitative reverse transcription-polymerase chain reaction and Western blotting, respectively. RESULTS: We found significant associations between SLE and genetic variants within TNIP1, but not in TNIP2 or TAX1BP1. After resequencing and imputation, we identified 2 independent risk haplotypes within TNIP1 in individuals of European ancestry that were also present in African American and Hispanic populations. Levels of TNIP1 mRNA and ABIN1 protein were reduced among subjects with these haplotypes, suggesting that they harbor hypomorphic functional variants that influence susceptibility to SLE by restricting ABIN1 expression. CONCLUSION: Our results confirm the association signals between SLE and TNIP1 variants in multiple populations and provide new insight into the mechanism by which TNIP1 variants may contribute to SLE pathogenesis.

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

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.

Identification of Environmental Exposures Associated with Risk of Sarcoidosis in African Americans.

Rationale: Sarcoidosis is a racially disparate granulomatous disease likely caused by environmental exposures, genes, and their interactions. Despite increased risk in African Americans, few environmental risk factor studies in this susceptible population exist. Objectives: To identify environmental exposures associated with the risk of sarcoidosis in African Americans and those that differ in effect by self-identified race and genetic ancestry. Methods: The study sample comprised 2,096 African Americans (1,205 with and 891 without sarcoidosis) compiled from three component studies. Unsupervised clustering and multiple correspondence analyses were used to identify underlying clusters of environmental exposures. Mixed-effects logistic regression was used to evaluate the association of these exposure clusters and the 51 single-component exposures with risk of sarcoidosis. A comparison case-control sample of 762 European Americans (388 with and 374 without sarcoidosis) was used to assess heterogeneity in exposure risk by race. Results: Seven exposure clusters were identified, five of which were associated with risk. The exposure cluster with the strongest risk association was composed of metals (P < 0.001), and within this cluster, exposure to aluminum had the highest risk (odds ratio, 3.30; 95% confidence interval [95% CI], 2.23-4.09; P < 0.001). This effect also differed by race (P < 0.001), with European Americans having no significant association with exposure (odds ratio, 0.86; 95% CI, 0.56-1.33). Within African Americans, the increased risk was dependent on genetic African ancestry (P = 0.047). Conclusions: Our findings support African Americans having sarcoidosis environmental exposure risk profiles that differ from those of European Americans. These differences may underlie racially disparate incidence rates that are partially explained by genetic variation differing by African ancestry.

Low Gut Microbial Diversity Augments Estrogen-driven Pulmonary Fibrosis in Female-Predominant Interstitial Lung Disease.

Although profibrotic cytokines such as IL-17A and TGF-ß1 have been implicated in interstitial lung disease (ILD) pathogenesis, interactions between gut dysbiosis, gonadotrophic hormones and molecular mediators of profibrotic cytokine expression, such as phosphorylation of STAT3, have not been defined. Here we show by chromatin immunoprecipitation sequencing (ChIP-seq) analysis of primary human CD4+ T cells that regions within the STAT3 locus are significantly enriched for binding by the transcription factor estrogen receptor alpha (ERa). Using the murine model of bleomycin-induced pulmonary fibrosis, we found significantly increased regulatory T cells compared to Th17 cells in the female lung. Genetic absence of ESR1 or ovariectomy in mice significantly increased pSTAT3 and IL-17A expression in pulmonary CD4+ T cells, which was reduced after repletion of female hormones. Remarkably, there was no significant reduction in lung fibrosis under either condition, suggesting that factors outside of ovarian hormones also contribute. Assessment of lung fibrosis among menstruating females in different rearing environments revealed that environments favoring gut dysbiosis augment fibrosis. Furthermore, hormone repletion following ovariectomy further augmented lung fibrosis, suggesting pathologic interactions between gonadal hormones and gut microbiota on lung fibrosis severity. Analysis in female sarcoidosis patients revealed a significant reduction in pSTAT3 and IL-17A levels and a concomitant increase in TGF-ß1 levels in CD4+ T cells, compared to male sarcoidosis patients. These studies reveal that estrogen is profibrotic in females and that gut dysbiosis in menstruating females augments lung fibrosis severity, supporting a critical interaction between gonadal hormones and gut flora in lung fibrosis pathogenesis.