Admixture mapping analysis reveals differential genetic ancestry associated with Chagas disease susceptibility in the Colombian population.

Chagas disease is an infection caused by the parasite Trypanosoma cruzi, endemic in Latino America. Leveraging the three-way admixture between Native American (AMR), European (EUR) and African (AFR) populations in Latin Americans, we aimed to better understand the genetic basis of Chagas disease by performing an admixture mapping study in a Colombian population. A two-stage study was conducted, and subjects were classified as seropositive and seronegative for T. cruzi. In stage 1, global and local ancestries were estimated using reference data from the 1000 Genomes Project (1KGP), and local ancestry associations were performed by logistic regression models. The AMR ancestry showed a protective association with Chagas disease within the major histocompatibility complex region [Odds ratio (OR) = 0.74, 95% confidence interval (CI) = 0.66-0.83, lowest P-value = 4.53 × 10-8]. The fine mapping assessment on imputed genotypes combining data from stage 1 and 2 from an independent Colombian cohort, revealed nominally associated variants in high linkage disequilibrium with the top signal (rs2032134, OR = 0.93, 95% CI = 0.90-0.97, P-value = 3.54 × 10-4) in the previously associated locus. To assess ancestry-specific adaptive signals, a selective sweep scan in an AMR reference population from 1KGP together with an in silico functional analysis highlighted the Tripartite Motif family and the human leukocyte antigen genes, with crucial role in the immune response against pathogens. Furthermore, these analyses emphasized the macrophages, neutrophils and eosinophils, as key players in the defense against T. cruzi. This first admixture mapping study in Chagas disease provided novel insights underlying the host immune response in the pathogenesis of this neglected disease.

Non-classical circulating monocytes expressing high levels of microsomal prostaglandin E2 synthase-1 tag an aberrant IFN-response in systemic sclerosis.

Systemic sclerosis (SSc) is a complex disease that affects the connective tissue, causing fibrosis. SSc patients show altered immune cell composition and activation in the peripheral blood (PB). PB monocytes (Mos) are recruited into tissues where they differentiate into macrophages, which are directly involved in fibrosis. To understand the role of CD14+ PB Mos in SSc, a single-cell transcriptome analysis (scRNA-seq) was conducted on 8 SSc patients and 8 controls. Using unsupervised clustering methods, CD14+ cells were assigned to 11 clusters, which added granularity to the known monocyte subsets: classical (cMos), intermediate (iMos) and non-classical Mos (ncMos) or type 2 dendritic cells. NcMos were significantly overrepresented in SSc patients and showed an active IFN-signature and increased expression levels of PTGES, in addition to monocyte motility and adhesion markers. We identified a SSc-related cluster of IRF7+ STAT1+ iMos with an aberrant IFN-response. Finally, a depletion of M2 polarised cMos in SSc was observed. Our results highlighted the potential of PB Mos as biomarkers for SSc and provided new possibilities for putative drug targets for modulating the innate immune response in SSc.

FLT3 functional low-frequency variant rs76428106-C is associated with susceptibility to systemic sclerosis.

OBJECTIVES: rs76428106-C, a low frequency polymorphism that affects the splicing of the FLT3 gene, has recently been associated with several seropositive autoimmune diseases. Here, we aimed to evaluate the potential implication of rs76428106-C in the susceptibility to systemic sclerosis (SSc). METHODS: We analysed a total of 26 598 European ancestry individuals, 9063 SSc and 17 535 healthy controls, to test the association between FLT3 rs76428106-C and SSc and its different subphenotypes. Genotype data of rs76428106 were obtained by imputation of already available genome-wide association study data and analysed by logistic regression analysis. RESULTS: In accordance with that observed in other autoimmune disorders, the FLT3 rs76428106-C allele was significantly increased [P-value = 2.03 × 10-3, odds ratio (OR) = 1.34] in SSc patients compared with healthy controls. A similar risk effect was found when the main SSc clinical and serological subgroups were compared with controls. When comparing SSc patients with and without digital ulcers (DU), the rs76428106-C frequency was significantly increased in DU-positive SSc patients in comparison with DU-negative patients (P-value = 0.036, OR = 2.16). CONCLUSION: This study is the first to report an association between rs76428176-C and SSc. Our results support the role of FLT3 as a relevant gene in seropositive immune-mediated diseases and a potential biomarker for SSc microangiopathy.

Contribution of Telomere Length to Systemic Sclerosis Onset: A Mendelian Randomization Study.

Although previous studies have suggested a relationship between telomere shortening and systemic sclerosis (SSc), the association between these two traits remains poorly understood. The objective of this study was to assess the causal relationship between telomere length in leukocytes (LTL) and SSc using the two-sample Mendelian randomization approach, with the genome-wide association study data for both LTL and SSc. The results of inverse-variance weighted regression (OR = 0.716 [95% CI 0.528-0.970], p = 0.031) and the Mendelian randomization pleiotropy residual sum and outlier method (OR = 0.716 [95% CI 0.563-0.911], p = 0.035) indicate an association between telomere length and SSc. Specifically, longer genetically predicted LTL is associated with a reduced risk of SSc. Sensitivity tests highlight the significant roles of the variants rs10936599 and rs2736100 annotated to the TERC and TERT genes, respectively. Our findings suggest an influence of telomere length in leukocytes on the development of SSc.

A Genome-Wide Association Study Identifies Novel Susceptibility loci in Chronic Chagas Cardiomyopathy.

BACKGROUND: Chagas disease is an infectious disease caused by the parasite Trypanosoma cruzi and is endemic from Latin American countries. The goal of our study was to identify novel genetic loci associated with chronic Chagas cardiomyopathy development in Chagas disease patients from different Latin American populations. METHODS: We performed a cross-sectional, nested case-control study including 3 sample collections from Colombia, Argentina, and Bolivia. Samples were genotyped to conduct a genome-wide association study (GWAS). These results were meta-analyzed with summary statistic data from Brazil, gathering a total of 3413 Chagas disease patients. To identify the functional impact of the associated variant and its proxies, we performed an in silico analysis of this region. RESULTS: The meta-analysis revealed a novel genome-wide statistically significant association with chronic Chagas cardiomyopathy development in rs2458298 (OR = 0.90, 95%CI = 0.87-0.94, P-value = 3.27 × 10-08), nearby the SAC3D1 gene. In addition, further in silico analyses displayed functional relationships between the associated variant and the SNX15, BAFT2, and FERMT3 genes, related to cardiovascular traits. CONCLUSIONS: Our findings support the role of the host genetic factors in the susceptibility to the development of the chronic cardiac form of this neglected disease.

Comprehensive analysis of the major histocompatibility complex in systemic sclerosis identifies differential HLA associations by clinical and serological subtypes.

OBJECTIVE: The greatest genetic effect reported for systemic sclerosis (SSc) lies in the major histocompatibility complex (MHC) locus. Leveraging the largest SSc genome-wide association study, we aimed to fine-map this region to identify novel human leucocyte antigen (HLA) genetic variants associated with SSc susceptibility and its main clinical and serological subtypes. METHODS: 9095 patients with SSc and 17 584 controls genome-wide genotyped were used to impute and test single-nucleotide polymorphisms (SNPs) across the MHC, classical HLA alleles and their composite amino acid residues. Additionally, patients were stratified according to their clinical and serological status, namely, limited cutaneous systemic sclerosis (lcSSc), diffuse cutaneous systemic sclerosis (dcSSc), anticentromere (ACA), antitopoisomerase (ATA) and anti-RNApolIII autoantibodies (ARA). RESULTS: Sequential conditional analyses showed nine SNPs, nine classical alleles and seven amino acids that modelled the observed associations with SSc. This confirmed previously reported associations with HLA-DRB1*11:04 and HLA-DPB1*13:01, and revealed a novel association of HLA-B*08:01. Stratified analyses showed specific associations of HLA-DQA1*02:01 with lcSSc, and an exclusive association of HLA-DQA1*05:01 with dcSSc. Similarly, private associations were detected in HLA-DRB1*08:01 and confirmed the previously reported association of HLA-DRB1*07:01 with ACA-positive patients, as opposed to the HLA-DPA1*02:01 and HLA-DQB1*03:01 alleles associated with ATA presentation. CONCLUSIONS: This study confirms the contribution of HLA class II and reveals a novel association of HLA class I with SSc, suggesting novel pathways of disease pathogenesis. Furthermore, we describe specific HLA associations with SSc clinical and serological subtypes that could serve as biomarkers of disease severity and progression.

Genomic Risk Score impact on susceptibility to systemic sclerosis.

OBJECTIVES: Genomic Risk Scores (GRS) successfully demonstrated the ability of genetics to identify those individuals at high risk for complex traits including immune-mediated inflammatory diseases (IMIDs). We aimed to test the performance of GRS in the prediction of risk for systemic sclerosis (SSc) for the first time. METHODS: Allelic effects were obtained from the largest SSc Genome-Wide Association Study (GWAS) to date (9 095 SSc and 17 584 healthy controls with European ancestry). The best-fitting GRS was identified under the additive model in an independent cohort that comprised 400 patients with SSc and 571 controls. Additionally, GRS for clinical subtypes (limited cutaneous SSc and diffuse cutaneous SSc) and serological subtypes (anti-topoisomerase positive (ATA+) and anti-centromere positive (ACA+)) were generated. We combined the estimated GRS with demographic and immunological parameters in a multivariate generalised linear model. RESULTS: The best-fitting SSc GRS included 33 single nucleotide polymorphisms (SNPs) and discriminated between patients with SSc and controls (area under the receiver operating characteristic (ROC) curve (AUC)=0.673). Moreover, the GRS differentiated between SSc and other IMIDs, such as rheumatoid arthritis and Sjögren's syndrome. Finally, the combination of GRS with age and immune cell counts significantly increased the performance of the model (AUC=0.787). While the SSc GRS was not able to discriminate between ATA+ and ACA+ patients (AUC<0.5), the serological subtype GRS, which was based on the allelic effects observed for the comparison between ACA+ and ATA+ patients, reached an AUC=0.693. CONCLUSIONS: GRS was successfully implemented in SSc. The model discriminated between patients with SSc and controls or other IMIDs, confirming the potential of GRS to support early and differential diagnosis for SSc.

Polymorphisms in CEP68 gene associated with risk of immediate selective reactions to non-steroidal anti-inflammatory drugs.

Non-steroidal anti-inflammatory drugs (NSAIDs) are the main triggers of drug hypersensitivity reactions. Such reactions can be pharmacologically or immunologically mediated, but in both cases individual susceptibility can be influenced by genetic factors. Polymorphisms in centrosomal protein of 68 kDa (CEP68) have been associated with pharmacologically mediated NSAIDs reactions. Here, we evaluated this gene in immunologically mediated single-NSAID-induced urticaria/angioedema or anaphylaxis (SNIUAA) by analyzing 52 single nucleotide polymorphisms in CEP68 in 176 patients and 363 NSAIDs-tolerant controls. Two intronic variants (rs2241160 and rs2241161) were significantly associated with an increased risk of SNIUAA, suggesting CEP68 to be a key player in both types of NSAIDs hypersensitivity. However, we found no overlap with genetic variants previously associated with pharmacologically mediated hypersensitivity, pointing to a complex role for this gene and its potential use in the development of biomarkers of clinical utility to diagnose patients at risk of these reactions and to differentiate entities.

Genome-wide meta-analysis reveals shared new loci in systemic seropositive rheumatic diseases.

OBJECTIVE: Immune-mediated inflammatory diseases (IMIDs) are heterogeneous and complex conditions with overlapping clinical symptoms and elevated familial aggregation, which suggests the existence of a shared genetic component. In order to identify this genetic background in a systematic fashion, we performed the first cross-disease genome-wide meta-analysis in systemic seropositive rheumatic diseases, namely, systemic sclerosis, systemic lupus erythematosus, rheumatoid arthritis and idiopathic inflammatory myopathies. METHODS: We meta-analysed ~6.5 million single nucleotide polymorphisms in 11 678 cases and 19 704 non-affected controls of European descent populations. The functional roles of the associated variants were interrogated using publicly available databases. RESULTS: Our analysis revealed five shared genome-wide significant independent loci that had not been previously associated with these diseases: NAB1, KPNA4-ARL14, DGQK, LIMK1 and PRR12. All of these loci are related with immune processes such as interferon and epidermal growth factor signalling, response to methotrexate, cytoskeleton dynamics and coagulation cascade. Remarkably, several of the associated loci are known key players in autoimmunity, which supports the validity of our results. All the associated variants showed significant functional enrichment in DNase hypersensitivity sites, chromatin states and histone marks in relevant immune cells, including shared expression quantitative trait loci. Additionally, our results were significantly enriched in drugs that are being tested for the treatment of the diseases under study. CONCLUSIONS: We have identified shared new risk loci with functional value across diseases and pinpoint new potential candidate loci that could be further investigated. Our results highlight the potential of drug repositioning among related systemic seropositive rheumatic IMIDs.

Towards a Better Classification and Novel Therapies Based on the Genetics of Systemic Sclerosis.

PURPOSE OF THE REVIEW: Nowadays, important advances have occurred in our understanding of the pathogenesis of systemic sclerosis (SSc), which is a rare immune-mediated inflammatory disease (IMID) characterized by vascular damage, immune imbalance, and fibrosis. Its etiology remains unknown; nevertheless, both environmental and genetic factors play a major role in the disease. This review will focus on the main advances made in the field of genetics of SSc. RECENT FINDINGS: The assessment of how interindividual genetic variability affects disease onset and progression has enhanced our knowledge of disease biology, and this will eventually translate in the development of new diagnostic and therapeutic tools, which is the final goal of personalized medicine. We will provide an overview of the most relevant achievements in the genetics of SSc, its shared genetics among IMIDs with special attention on drug repurposing, current challenges for the functional characterization of risk variants, and future directions.

Correction to: Towards a Better Classification and Novel Therapies Based on the Genetics of Systemic Sclerosis.

The original version of this article unfortunately contained a mistake. The legend of Fig. 1 was incorrect.

Genetic Variants of Thymic Stromal Lymphopoietin in Nonsteroidal Anti-Inflammatory Drug-Induced Urticaria/Angioedema.

BACKGROUND: Nonsteroidal anti-inflammatory drugs (NSAIDs) are the most frequent agents involved in hypersensitivity drug reactions, with NSAID-induced urticaria and/or angioedema (NIUA) being the most common entity. Mast cells are key players in NIUA and are activated by thymic stromal lymphopoietin (TSLP). This cytokine functions through recognition by its receptor, composed of IL7Rα (interleukin-7 receptor alpha) and TSLPR (TSLP receptor). These genes have been previously associated with other inflammatory diseases. METHODS: We assessed the genetic association between single nucleotide polymorphisms (SNPs) in TSLP, IL7R and TSLPR and NIUA in Spanish individuals, using genotyped and imputed data. A total of 369 unrelated NIUA patients and 580 NSAID-tolerant control subjects were included, and 6 SNPs in TSLP, 6 in IL7R and 3 in TSLPR were genotyped. Further variants were imputed using Mach and the 1,000 Genomes Project (Phase 3) data. Association testing and statistical analyses were performed with Mach2dat and R. RESULTS: A total of 139 SNPs were tested for association following quality control. Two SNPs in TSLP (rs1816678 and rs764917) showed a nominal association (p = 0.033 and 0.024, respectively) with NIUA, although these results were not statistically significant after correcting for multiple comparisons. CONCLUSIONS: Although TSLP, IL7R and TSLPR are important genes involved in the development of the inflammatory response, we found no significant genetic association with NIUA in our population for common SNPs in these genes.

Common variants of NFE2L2 gene predisposes to acute respiratory distress syndrome in patients with severe sepsis.

INTRODUCTION: The purpose of this study was to investigate whether common variants across the nuclear factor erythroid 2-like 2 (NFE2L2) gene contribute to the development of the acute respiratory distress syndrome (ARDS) in patients with severe sepsis. NFE2L2 is involved in the response to oxidative stress, and it has been shown to be associated with the development of ARDS in trauma patients. METHODS: We performed a case-control study of 321 patients fulfilling international criteria for severe sepsis and ARDS who were admitted to a Spanish network of post-surgical and critical care units, as well as 871 population-based controls. Six tagging single-nucleotide polymorphisms (SNPs) of NFE2L2 were genotyped, and, after further imputation of additional 34 SNPs, association testing with ARDS susceptibility was conducted using logistic regression analysis. RESULTS: After multiple testing adjustments, our analysis revealed 10 non-coding SNPs in tight linkage disequilibrium (0.75 ≤ r (2) ≤ 1) that were associated with ARDS susceptibility as a single association signal. One of those SNPs (rs672961) was previously associated with trauma-induced ARDS and modified the promoter activity of the NFE2L2 gene, showing an odds ratio of 1.93 per T allele (95 % confidence interval, 1.17-3.18; p = 0.0089). CONCLUSIONS: Our findings support the involvement of NFE2L2 gene variants in ARDS susceptibility and reinforce further exploration of the role of oxidant stress response as a risk factor for ARDS in critically ill patients.

Genome-Wide Association Study Identifies the First Germline Genetic Variant Associated With Erdheim-Chester Disease.

OBJECTIVE: Erdheim-Chester disease (ECD) is rare histiocytosis with a wide range of clinical manifestations. Somatic mutations are key to the pathogenesis of the disease; however, the relationship between germline genetic variants and ECD has not been examined so far. The present study aims to explore the inherited genetic component of ECD by performing the first genome-wide association study. METHODS: After quality controls, a cohort of 255 patients with ECD and 7,471 healthy donors was included in this study. Afterward, a logistic regression followed by in silico functional annotation was performed. RESULTS: A signal at the 18q12.3 genomic region was identified as a new susceptibility locus for ECD (P = 2.75 × 10-11 ; Odds Ratio = 2.09). This association was annotated to the SETBP1 gene, which is involved in clonal haematopoiesis. Functional annotation of this region and of the identified suggestive signals revealed additional genes that could be potentially involved in the pathogenesis of the disease. CONCLUSION: Overall, this work demonstrates that germline genetic variants can impact on the development of ECD and suggests new pathways with a potential pathogenic role.

Functional promoter variants in sphingosine 1-phosphate receptor 3 associate with susceptibility to sepsis-associated acute respiratory distress syndrome.

The genetic mechanisms underlying the susceptibility to acute respiratory distress syndrome (ARDS) are poorly understood. We previously demonstrated that sphingosine 1-phosphate (S1P) and the S1P receptor S1PR3 are intimately involved in lung inflammatory responses and vascular barrier regulation. Furthermore, plasma S1PR3 protein levels were shown to serve as a biomarker of severity in critically ill ARDS patients. This study explores the contribution of single nucleotide polymorphisms (SNPs) of the S1PR3 gene to sepsis-associated ARDS. S1PR3 SNPs were identified by sequencing the entire gene and tagging SNPs selected for case-control association analysis in African- and ED samples from Chicago, with independent replication in a European case-control study of Spanish individuals. Electrophoretic mobility shift assays, luciferase activity assays, and protein immunoassays were utilized to assess the functionality of associated SNPs. A total of 80 variants, including 29 novel SNPs, were identified. Because of limited sample size, conclusive findings could not be drawn in African-descent ARDS subjects; however, significant associations were found for two promoter SNPs (rs7022797 -1899T/G; rs11137480 -1785G/C), across two ED samples supporting the association of alleles -1899G and -1785C with decreased risk for sepsis-associated ARDS. In addition, these alleles significantly reduced transcription factor binding to the S1PR3 promoter; reduced S1PR3 promoter activity, a response particularly striking after TNF-α challenge; and were associated with lower plasma S1PR3 protein levels in ARDS patients. These highly functional studies support S1PR3 as a novel ARDS candidate gene and a potential target for individualized therapy.

Complement component C4 structural variation and quantitative traits contribute to sex-biased vulnerability in systemic sclerosis.

Copy number (CN) polymorphisms of complement C4 play distinct roles in many conditions, including immune-mediated diseases. We investigated the association of C4 CN with systemic sclerosis (SSc) risk. Imputed total C4, C4A, C4B, and HERV-K CN were analyzed in 26,633 individuals and validated in an independent cohort. Our results showed that higher C4 CN confers protection to SSc, and deviations from CN parity of C4A and C4B augmented risk. The protection contributed per copy of C4A and C4B differed by sex. Stronger protection was afforded by C4A in men and by C4B in women. C4 CN correlated well with its gene expression and serum protein levels, and less C4 was detected for both in SSc patients. Conditioned analysis suggests that C4 genetics strongly contributes to the SSc association within the major histocompatibility complex locus and highlights classical alleles and amino acid variants of HLA-DRB1 and HLA-DPB1 as C4-independent signals.

The Effect of Body Fat Distribution on Systemic Sclerosis.

Obesity contributes to a chronic proinflammatory state, which is a known risk factor to develop immune-mediated diseases. However, its role in systemic sclerosis (SSc) remains to be elucidated. Therefore, we conducted a two-sample mendelian randomization (2SMR) study to analyze the effect of three body fat distribution parameters in SSc. As instrumental variables, we used the allele effects described for single nucleotide polymorphisms (SNPs) in different genome-wide association studies (GWAS) for SSc, body mass index (BMI), waist-to-hip ratio (WHR) and WHR adjusted for BMI (WHRadjBMI). We performed local (pHESS) and genome-wide (LDSC) genetic correlation analyses between each of the traits and SSc and we applied several Mendelian randomization (MR) methods (i.e., random effects inverse-variance weight, MR-Egger regression, MR pleiotropy residual sum and outlier method and a multivariable model). Our results show no genetic correlation or causal relationship between any of these traits and SSc. Nevertheless, we observed a negative causal association between WHRadjBMI and SSc, which might be due to the effect of gastrointestinal complications suffered by the majority of SSc patients. In conclusion, reverse causality might be an especially difficult confounding factor to define the effect of obesity in the onset of SSc.