NFKB2 polymorphisms associate with the risk of developing rheumatoid arthritis and response to TNF inhibitors: Results from the REPAIR consortium

Polymorphisms within autophagy-related genes as susceptibility biomarkers for pancreatic cancer: A meta-analysis of three large European cohorts and functional characterization

Pancreatic ductal adenocarcinoma (PDAC) is one of the most lethal cancers with patients having unresectable or metastatic disease at diagnosis, with poor prognosis and very short survival. Given that genetic variation within autophagy-related genes influences autophagic flux and susceptibility to solid cancers, we decided to investigate whether 55,583 single nucleotide polymorphisms (SNPs) within 234 autophagy-related genes could influence the risk of developing PDAC in three large independent cohorts of European ancestry including 12,754 PDAC cases and 324,926 controls. The meta-analysis of these populations identified, for the first time, the association of the BIDrs9604789 variant with an increased risk of developing the disease (ORMeta = 1.31, p = 9.67 × 10-6). We also confirmed the association of TP63rs1515496 and TP63rs35389543 variants with PDAC risk (OR = 0.89, p = 6.27 × 10-8 and OR = 1.16, p = 2.74 × 10-5). Although it is known that BID induces autophagy and TP63 promotes cell growth, cell motility and invasion, we also found that carriers of the TP63rs1515496G allele had increased numbers of FOXP3+ Helios+ T regulatory cells and CD45RA+ T regulatory cells (p = 7.67 × 10-4 and p = 1.56 × 10-3), but also decreased levels of CD4+ T regulatory cells (p = 7.86 × 10-4). These results were in agreement with research suggesting that the TP63rs1515496 variant alters binding sites for FOXA1 and CTCF, which are transcription factors involved in modulating specific subsets of regulatory T cells. In conclusion, this study identifies BID as new susceptibility locus for PDAC and confirms previous studies suggesting that the TP63 gene is involved in the development of PDAC. This study also suggests new pathogenic mechanisms of the TP63 locus in PDAC.

Ovol1/2 loss-induced epidermal defects elicit skin immune activation and alter global metabolism

Skin epidermis constitutes the outer permeability barrier that protects the body from dehydration, heat loss, and myriad external assaults. Mechanisms that maintain barrier integrity in constantly challenged adult skin and how epidermal dysregulation shapes the local immune microenvironment and whole-body metabolism remain poorly understood. Here, we demonstrate that inducible and simultaneous ablation of transcription factor-encoding Ovol1 and Ovol2 in adult epidermis results in barrier dysregulation through impacting epithelial-mesenchymal plasticity and inflammatory gene expression. We find that aberrant skin immune activation then ensues, featuring Langerhans cell mobilization and T cell responses, and leading to elevated levels of secreted inflammatory factors in circulation. Finally, we identify failure to gain body weight and accumulate body fat as long-term consequences of epidermal-specific Ovol1/2 loss and show that these global metabolic changes along with the skin barrier/immune defects are partially rescued by immunosuppressant dexamethasone. Collectively, our study reveals key regulators of adult barrier maintenance and suggests a causal connection between epidermal dysregulation and whole-body metabolism that is in part mediated through aberrant immune activation.

Polymorphisms within Autophagy-Related Genes as Susceptibility Biomarkers for Multiple Myeloma: A Meta-Analysis of Three Large Cohorts and Functional Characterization

Multiple myeloma (MM) arises following malignant proliferation of plasma cells in the bone marrow, that secrete high amounts of specific monoclonal immunoglobulins or light chains, resulting in the massive production of unfolded or misfolded proteins. Autophagy can have a dual role in tumorigenesis, by eliminating these abnormal proteins to avoid cancer development, but also ensuring MM cell survival and promoting resistance to treatments. To date no studies have determined the impact of genetic variation in autophagy-related genes on MM risk. We performed meta-analysis of germline genetic data on 234 autophagy-related genes from three independent study populations including 13,387 subjects of European ancestry (6863 MM patients and 6524 controls) and examined correlations of statistically significant single nucleotide polymorphisms (SNPs; p < 1 × 10^-9) with immune responses in whole blood, peripheral blood mononuclear cells (PBMCs), and monocyte-derived macrophages (MDM) from a large population of healthy donors from the Human Functional Genomic Project (HFGP). We identified SNPs in six loci, CD46, IKBKE, PARK2, ULK4, ATG5, and CDKN2A associated with MM risk (p = 4.47 × 10^-4-5.79 × 10^-14). Mechanistically, we found that the ULK4_rs6599175 SNP correlated with circulating concentrations of vitamin D3 (p = 4.0 × 10^-4), whereas the IKBKE_rs17433804 SNP correlated with the number of transitional CD24⁺CD38⁺ B cells (p = 4.8 × 10^-4) and circulating serum concentrations of Monocyte Chemoattractant Protein (MCP)-2 (p = 3.6 × 10^-4). We also found that the CD46_rs1142469 SNP correlated with numbers of CD19⁺ B cells, CD19⁺CD3^- B cells, CD5⁺IgD^- cells, IgM^- cells, IgD^-IgM^- cells, and CD4^-CD8^- PBMCs (p = 4.9 × 10^-4-8.6 × 10^-4) and circulating concentrations of interleukin (IL)-20 (p = 0.00082). Finally, we observed that the CDKN2A_rs2811710 SNP correlated with levels of CD4⁺EMCD45RO⁺CD27^- cells (p = 9.3 × 10^-4). These results suggest that genetic variants within these six loci influence MM risk through the modulation of specific subsets of immune cells, as well as vitamin D3^-, MCP-2^-, and IL20-dependent pathways.

Inflammasome Genetic Variants Are Associated with Protection to Clinical Severity of COVID-19 among Patients from Rio de Janeiro, Brazil

COVID-19 has a broad spectrum of clinical manifestations, from asymptomatic or mild/moderate symptoms to severe symptoms and death. The mechanisms underlying its clinical evolution are still unclear. Upon SARS-CoV-2 infection, host factors, such as the inflammasome system, are activated by the presence of the virus inside host cells. The search for COVID-19 risk factors is of relevance for clinical management. In this study, we investigated the impact of inflammasome single-nucleotide polymorphisms (SNPs) in SARS-CoV-2-infected individuals with distinct severity profiles at clinical presentation. Patients were divided into two groups according to disease severity at clinical presentation based on the WHO Clinical Progression Scale. Group 1 included patients with mild/moderate disease (WHO < 6; n = 76), and group 2 included patients with severe/critical COVID-19 (WHO ≥ 6; n = 357). Inpatients with moderate to severe/critical profiles were recruited and followed-up at Hospital Center for COVID-19 Pandemic – National Institute of Infectology (INI)/FIOCRUZ, RJ, Brazil, from June 2020 to March 2021. Patients with mild disease were recruited at Oswaldo Cruz Institute (IOC)/FIOCRUZ, RJ, Brazil, in August 2020. Genotyping of 11 inflammasome SNPs was determined by real-time PCR. Protection and risk estimation were performed using unconditional logistic regression models. Significant differences in NLRP3 rs1539019 and CARD8 rs2043211 were observed between the two groups. Protection against disease severity was associated with the A/A genotype (OR_adj = 0.36; P = 0.032), allele A (OR_adj = 0.93; P = 0.010), or carrier-A (OR_adj = 0.45; P = 0.027) in the NLRP3 rs1539019 polymorphism; A/T genotype (OR_adj = 0.5; P = 0.045), allele T (OR_adj = 0.93; P = 0.018), or carrier-T (OR_adj = 0.48; P = 0.029) in the CARD8 rs2043211 polymorphism; and the A-C-G-C-C (OR_adj = 0.11; P = 0.018), A-C-G-C-G (OR_adj = 0.23; P = 0.003), C-C-G-C-C (OR_adj = 0.37; P = 0.021), and C-T-G-A-C (OR_adj = 0.04; P = 0.0473) in NLRP3 genetic haplotype variants. No significant associations were observed for the other polymorphisms. To the best of our knowledge, this is the first study demonstrating an association between CARD8 and NLRP3 inflammasome genetic variants and protection against COVID-19 severity, contributing to the discussion of the impact of inflammasomes on COVID-19 outcomes.

Mechanisms underlying DMARD inefficacy in difficult-to-treat rheumatoid arthritis: a narrative review with systematic literature search

Management of RA patients has significantly improved over the past decades. However, a substantial proportion of patients is difficult-to-treat (D2T), remaining symptomatic after failing biological and/or targeted synthetic DMARDs. Multiple factors can contribute to D2T RA, including treatment non-adherence, comorbidities and co-existing mimicking diseases (e.g. fibromyalgia). Additionally, currently available biological and/or targeted synthetic DMARDs may be truly ineffective ('true' refractory RA) and/or lead to unacceptable side effects. In this narrative review based on a systematic literature search, an overview of underlying (immune) mechanisms is presented. Potential scenarios are discussed including the influence of different levels of gene expression and clinical characteristics. Although the exact underlying mechanisms remain largely unknown, the heterogeneity between individual patients supports the assumption that D2T RA is a syndrome involving different pathogenic mechanisms.

Mitochondrial Dysfunction and Oxidative Stress in Rheumatoid Arthritis

Control of excessive mitochondrial oxidative stress could provide new targets for both preventive and therapeutic interventions in the treatment of chronic inflammation or any pathology that develops under an inflammatory scenario, such as rheumatoid arthritis (RA). Increasing evidence has demonstrated the role of mitochondrial alterations in autoimmune diseases mainly due to the interplay between metabolism and innate immunity, but also in the modulation of inflammatory response of resident cells, such as synoviocytes. Thus, mitochondrial dysfunction derived from several danger signals could activate tricarboxylic acid (TCA) disruption, thereby favoring a vicious cycle of oxidative/mitochondrial stress. Mitochondrial dysfunction can act through modulating innate immunity via redox-sensitive inflammatory pathways or direct activation of the inflammasome. Besides, mitochondria also have a central role in regulating cell death, which is deeply altered in RA. Additionally, multiple evidence suggests that pathological processes in RA can be shaped by epigenetic mechanisms and that in turn, mitochondria are involved in epigenetic regulation. Finally, we will discuss about the involvement of some dietary components in the onset and progression of RA.

Utility of Baseline Transcriptomic Analysis of Rheumatoid Arthritis Synovium as an Indicator for Long-Term Clinical Outcomes

Objective:

Rheumatoid arthritis is a chronic inflammatory autoimmune disease that can lead to synovial damage, persistent joint pain, and functional disability. Our objective was to evaluate baseline synovial transcriptome from early inflammatory arthritis patients (EIA) and identify pretreatment biomarkers that could potentially provide insights into long-term functional outcomes of rheumatoid arthritis (RA).

Methods

Synovial biopsies from clinically inflamed knee joints were procured from either 17 EIA patients before initiation of disease modifying anti-rheumatic drug (DMARD) therapy (DMARD-naïve EIA) using the minimally invasive closed needle biopsy technique or advanced RA patients undergoing arthroplasty. Affymetrix Human Genome U133 Plus 2 microarray platform was used to profile the synovial transcriptome. The cohort was followed clinically for a median of 12.3 years, and patient data was collected at each visit. Short-term and long-term clinical outcomes were determined by assessing RA-associated clinical parameters Statistical adjustments were made to account for asynchronous clinical visits and duration of follow up.

Results

Based on the transcriptomic analysis, we identified 5 differentially expressed genes (DEGs), including matrix metalloproteinase (MMP)-1 (fibroblast collagenase) and MMP-3 (stromelysin-1) in DMARD-naïve EIA patients, relative to advanced RA patients (q < 0.05). Dichotomous expression of MMP-1 and MMP-3 mRNA and protein was confirmed by qPCR and immunohistochemistry respectively, based on which DMARD-naïve EIA subjects were classified as MMP-high or MMP-low. Hierarchical clustering of transcriptomic data identified 947 DEGs between MMP-high and MMP-low cohorts. Co-expression and IPA analysis of DEGs in the MMP-high cohort showed an enrichment of genes that participated in metabolic or biochemical functions and intracellular immune signaling were regulated through NF-κB and β-catenin complexes and correlated with markers of systemic inflammation. Analysis of short-term clinical outcomes in MMP-high cohort showed a significant reduction in the DAS-CRP scores relative to baseline (P <0.001), whereas area under the curve analyses of modified HAQ (mHAQ) scores correlated negatively with baseline MMP-1 (R = −0.59, P = 0.03). Further, longitudinal mHAQ scores, number of swollen joints, number of DMARDs and median follow-up duration appeared to be higher in MMP-low cohort.

Conclusion

Overall, our results indicate that the gene expression profiling of synovial biopsies obtained at the DMARD-naive stage in patients with EIA categorizes them into subsets with varying degrees of inflammation and can predict the future of long-term clinical outcome.

Validation of GWAS-Identified Variants for Anti-TNF Drug Response in Rheumatoid Arthritis: A Meta-Analysis of Two Large Cohorts

We aimed to validate the association of 28 GWAS-identified genetic variants for response to TNF inhibitors (TNFi) in a discovery cohort of 1361 rheumatoid arthritis (RA) patients monitored in routine care and ascertained through the REPAIR consortium and DANBIO registry. We genotyped selected markers and evaluated their association with response to TNFi after 6 months of treatment according to the change in disease activity score 28 (ΔDAS28). Next, we confirmed the most interesting results through meta-analysis of our data with those from the DREAM cohort that included 706 RA patients treated with TNFi. The meta-analysis of the discovery cohort and DREAM registry including 2067 RA patients revealed an overall association of the LINC02549_rs7767069 SNP with a lower improvement in DAS28 that remained significant after correction for multiple testing (per-allele OR_Meta=0.83, P_Meta=0.000077; P_Het=0.61). In addition, we found that each copy of the LRRC55_rs717117G allele was significantly associated with lower improvement in DAS28 in rheumatoid factor (RF)-positive patients (per-allele OR_Meta=0.67, P=0.00058; P_Het=0.06) whereas an opposite but not significant effect was detected in RF-negative subjects (per-allele OR_Meta=1.38, P=0.10; P_Het=0.45; P_Interaction=0.00028). Interestingly, although the identified associations did not survive multiple testing correction, the meta-analysis also showed overall and RF-specific associations for the MAFB_rs6071980 and CNTN5_rs1813443 SNPs with decreased changes in DAS28 (per-allele OR_{Meta_rs6071980} = 0.85, P=0.0059; P_Het=0.63 and OR_{Meta_rs1813443_RF+}=0.81, P=0.0059; P_Het=0.69 and OR_{Meta_rs1813443_RF-}=1.00, P=0.99; P_Het=0.12; P_Interaction=0.032). Mechanistically, we found that subjects carrying the LINC02549_rs7767069T allele had significantly increased numbers of CD45RO+CD45RA+ T cells (P=0.000025) whereas carriers of the LINC02549_rs7767069T/T genotype showed significantly increased levels of soluble scavengers CD5 and CD6 in serum (P=0.00037 and P=0.00041). In addition, carriers of the LRRC55_rs717117G allele showed decreased production of IL6 after stimulation of PBMCs with B burgdorferi and E coli bacteria (P=0.00046 and P=0.00044), which suggested a reduced IL6-mediated anti-inflammatory effect of this marker to worsen the response to TNFi. In conclusion, this study confirmed the influence of the LINC02549 and LRRC55 loci to determine the response to TNFi in RA patients and suggested a weak effect of the MAFB and CNTN5 loci that need to be further investigated.

Identifying anti-TNF response biomarkers in ulcerative colitis using a diffusion-based signalling model

Motivation

Resistance to anti-TNF therapy in subgroups of ulcerative colitis (UC) patients is a major challenge and incurs significant treatment costs. Identification of patients at risk of nonresponse to anti-TNF is of major clinical importance. To date, no quantitative computational framework exists to develop a complex biomarker for the prognosis of UC treatment. Modelling patient-wise receptor to transcription factor (TF) network connectivity may enable personalized treatment.

Results

We present an approach for quantitative diffusion analysis between receptors and TFs using gene expression data. Key TFs were identified using pandaR. Network connectivities between immune-specific receptor-TF pairs were quantified using network diffusion in UC patients and controls. The patient-specific network could be considered a complex biomarker that separates anti-TNF treatment-resistant and responder patients both in the gene expression dataset used for model development and separate independent test datasets. The model was further validated in rheumatoid arthritis where it successfully discriminated resistant and responder patients to tocilizumab treatment. Our model may contribute to prognostic biomarkers that may identify treatment-resistant and responder subpopulations of UC patients.

Availability and implementation

Software is available at https://github.com/Amy3100/receptor2tfDiffusion.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

Polymorphisms within Autophagy-Related Genes Influence the Risk of Developing Colorectal Cancer: A Meta-Analysis of Four Large Cohorts

Simple Summary

We investigated the influence of autophagy-related variants in modulating colorectal cancer (CRC) risk through a meta-analysis of genome-wide association study (GWAS) data from four large European cohorts. We found that genetic variants within the DAPK2 and ATG5 loci were associated with CRC risk. This study also shed some light onto the functional mechanisms behind the observed associations and demonstrated the impact of DAPK2_rs11631973 and ATG5_rs546456 polymorphisms on the modulation of host immune responses, blood derived-cell counts and serum inflammatory protein levels, which might be involved in promoting cancer development. No effect of the DAPK2 and ATG5 polymorphisms on the autophagy flux was observed.

Abstract

The role of genetic variation in autophagy-related genes in modulating autophagy and cancer is poorly understood. Here, we comprehensively investigated the association of autophagy-related variants with colorectal cancer (CRC) risk and provide new insights about the molecular mechanisms underlying the associations. After meta-analysis of the genome-wide association study (GWAS) data from four independent European cohorts (8006 CRC cases and 7070 controls), two loci, DAPK2 (p = 2.19 × 10⁻⁵) and ATG5 (p = 6.28 × 10⁻⁴) were associated with the risk of CRC. Mechanistically, the DAPK2_rs11631973G allele was associated with IL1 β levels after the stimulation of peripheral blood mononuclear cells (PBMCs) with Staphylococcus aureus (p = 0.002), CD24 + CD38 + CD27 + IgM + B cell levels in blood (p = 0.0038) and serum levels of en-RAGE (p = 0.0068). ATG5_rs546456T allele was associated with TNF α and IL1 β levels after the stimulation of PBMCs with LPS (p = 0.0088 and p = 0.0076, respectively), CD14+CD16− cell levels in blood (p = 0.0068) and serum levels of CCL19 and cortisol (p = 0.0052 and p = 0.0074, respectively). Interestingly, no association with autophagy flux was observed. These results suggested an effect of the DAPK2 and ATG5 loci in the pathogenesis of CRC, likely through the modulation of host immune responses.

Germline Genetic Variants of Viral Entry and Innate Immunity May Influence Susceptibility to SARS-CoV-2 Infection: Toward a Polygenic Risk Score for Risk Stratification

The ongoing COVID-19 pandemic caused by the novel coronavirus, SARS-CoV-2 has affected all aspects of human society with a special focus on healthcare. Although older patients with preexisting chronic illnesses are more prone to develop severe complications, younger, healthy individuals might also exhibit serious manifestations. Previous studies directed to detect genetic susceptibility factors for earlier epidemics have provided evidence of certain protective variations. Following SARS-CoV-2 exposure, viral entry into cells followed by recognition and response by the innate immunity are key determinants of COVID-19 development. In the present review our aim was to conduct a thorough review of the literature on the role of single nucleotide polymorphisms (SNPs) as key agents affecting the viral entry of SARS-CoV-2 and innate immunity. Several SNPs within the scope of our approach were found to alter susceptibility to various bacterial and viral infections. Additionally, a multitude of studies confirmed genetic associations between the analyzed genes and autoimmune diseases, underlining the versatile immune consequences of these variants. Based on confirmed associations it is highly plausible that the SNPs affecting viral entry and innate immunity might confer altered susceptibility to SARS-CoV-2 infection and its complex clinical consequences. Anticipating several COVID-19 genomic susceptibility loci based on the ongoing genome wide association studies, our review also proposes that a well-established polygenic risk score would be able to clinically leverage the acquired knowledge.

Host immune genetic variations influence the risk of developing acute myeloid leukaemia: results from the NuCLEAR consortium

The purpose of this study was to conduct a two-stage case control association study including 654 acute myeloid leukaemia (AML) patients and 3477 controls ascertained through the NuCLEAR consortium to evaluate the effect of 27 immune-related single nucleotide polymorphisms (SNPs) on AML risk. In a pooled analysis of cohort studies, we found that carriers of the IL13_rs1295686A/A genotype had an increased risk of AML (P_Corr = 0.0144) whereas carriers of the VEGFA_rs25648T allele had a decreased risk of developing the disease (P_Corr = 0.00086). In addition, we found an association of the IL8_rs2227307 SNP with a decreased risk of developing AML that remained marginally significant after multiple testing (P_Corr = 0.072). Functional experiments suggested that the effect of the IL13_rs1295686 SNP on AML risk might be explained by its role in regulating IL1Ra secretion that modulates AML blast proliferation. Likewise, the protective effect of the IL8_rs2227307 SNP might be mediated by TLR2-mediated immune responses that affect AML blast viability, proliferation and chemorresistance. Despite the potential interest of these results, additional functional studies are still warranted to unravel the mechanisms by which these variants modulate the risk of AML. These findings suggested that IL13, VEGFA and IL8 SNPs play a role in modulating AML risk.

Biosynthesis and Chemical Diversification of Verucopeptin Leads to Structural and Functional Versatility

A synthesis program for structurally complex macrocycles is very challenging. Herein, we propose a biosynthesis pathway of the pyranylated cyclodepsipeptide verucopeptin to make enough supply and to diversify verucopeptin by genetic manipulation and one-step semisynthesis. The synthesis relies on the intrinsic reactivity of the interchangeable hemiketal pyrane and opened keto along with adjacent alkene. Biological evaluation of verucopeptin-oriented analogs delivers a potent AMP-activated protein kinase (AMPK) agonist, antibacterial agent, and selective NFκB modulator.