Identification of a Systemic Lupus Erythematosus Risk Locus Spanning ATG16L2, FCHSD2, and P2RY2 in Koreans

Type I interferon associated epistasis may contribute to early disease-onset and high disease activity in juvenile-onset lupus

Pathologic type I interferon (T1IFN) expression is a key feature in systemic lupus erythematosus (SLE) that associates with disease activity. When compared to adult-onset disease, juvenile-onset (j)SLE is characterized by increased disease activity and damage, which likely relates to increased genetic burden. To identify T1IFN-associated gene polymorphisms (TLR7, IRAK1, miR-3142/miR-146a, IRF5, IRF7, IFIH1, IRF8, TYK2, STAT4), identify long-range linkage disequilibrium and gene:gene interrelations, 319 jSLE patients were genotyped using panel sequencing. Coupling phenotypic quantitative trait loci (QTL) analysis identified 10 jSLE QTL that associated with young age at onset (<12 years; IRAK1 [rs1059702], TLR7 [rs3853839], IFIH1 [rs11891191, rs1990760, rs3747517], STAT4 [rs3021866], TYK2 [rs280501], IRF8 [rs1568391, rs6638]), global disease activity (SLEDAI-2 K >10; IFIH1 [rs1990760], STAT4 [rs3021866], IRF8 [rs903202, rs1568391, rs6638]), and mucocutaneous involvement (TLR7 [rs3853839], IFIH1 [rs11891191, rs1990760]). This study suggests T1IFN-associated polymorphisms and gene:gene interrelations in jSLE. Genotyping of jSLE patients may allow for individualized treatment and care.

Non-coding DNA variants for risk in lupus

Systemic Lupus Erythematosus (SLE) is a multifactorial autoimmune disease that arises from a dynamic interplay between genetics and environmental triggers. The advent of sophisticated genomics technology has catalyzed a shift in our understanding of disease etiology, spotlighting the pivotal role of non-coding DNA variants in SLE pathogenesis. In this review, we present a comprehensive examination of the non-coding variants associated with SLE, shedding light on their role in influencing disease risk and progression. We discuss the latest methodological advancements that have been instrumental in the identification and functional characterization of these genomic elements, with a special focus on the transformative power of CRISPR-based gene-editing technologies. Additionally, the review probes into the therapeutic opportunities that arise from modulating non-coding regions associated with SLE. Through an exploration of the complex network of non-coding DNA, this review aspires to decode the genetic puzzle of SLE and set the stage for groundbreaking gene-based therapeutic interventions and the advancement of precision medicine strategies tailored to SLE management.

Expression quantitative trait loci analysis in rheumatoid arthritis identifies tissue specific variants associated with severity and outcome

OBJECTIVE

Genome-wide association studies have successfully identified more than 100 loci associated with susceptibility to rheumatoid arthritis (RA). However, our understanding of the functional effects of genetic variants in causing RA and their effects on disease severity and response to treatment remains limited.

METHODS

In this study, we conducted expression quantitative trait locus (eQTL) analysis to dissect the link between genetic variants and gene expression comparing the disease tissue against blood using RNA-Sequencing of synovial biopsies (n=85) and blood samples (n=51) from treatment-naïve patients with RA from the Pathobiology of Early Arthritis Cohort.

RESULTS

This identified 898 eQTL genes in synovium and genes loci in blood, with 232 genes in common to both synovium and blood, although notably many eQTL were tissue specific. Examining the HLA region, we uncovered a specific eQTL at HLA-DPB2 with the critical triad of single-nucleotide polymorphisms (SNPs) rs3128921 driving synovial HLA-DPB2 expression, and both rs3128921 and HLA-DPB2 gene expression correlating with clinical severity and increasing probability of the lympho-myeloid pathotype.

CONCLUSIONS

This analysis highlights the need to explore functional consequences of genetic associations in disease tissue. HLA-DPB2 SNP rs3128921 could potentially be used to stratify patients to more aggressive treatment immediately at diagnosis.

An update on the roles of transcription factor Ets1 in autoimmune diseases

Transcription factors are crucial to regulate gene expression in immune cells and in other cell types. In lymphocytes, there are a large number of different transcription factors that are known to contribute to cell differentiation and the balance between quiescence and activation. One such transcription factor is E26 oncogene homolog 1 (Ets1). Ets1 expression is high in quiescent B and T lymphocytes and its levels are decreased upon activation. The human ETS1 gene has been identified as a susceptibility locus for many autoimmune and inflammatory diseases. In accord with this, gene knockout of Ets1 in mice leads to development of a lupus-like autoimmune disease, with enhanced activation and differentiation of both B cells and T cells. Prior reviews have summarized functional roles for Ets1 based on studies of Ets1 knockout mice. In recent years, numerous additional studies have been published that further validate ETS1 as a susceptibility locus for human diseases where immune dysregulation plays a causative role. In this update, new information that further links Ets1 to human autoimmune diseases is organized and collated to serve as a resource. This update also describes recent studies that seek to understand molecularly how Ets1 regulates immune cell activation, either using human cells and tissues or mouse models. This resource is expected to be useful to investigators seeking to understand how Ets1 may regulate the human immune response, particularly in terms of its roles in autoimmunity and inflammation. This article is categorized under: Immune System Diseases > Genetics/Genomics/Epigenetics Immune System Diseases > Molecular and Cellular Physiology.

HLA-DRB1*14:54 Is Associated with Pulmonary Alveolar Proteinosis: A Retrospective Real-World Audit

BACKGROUND

Pulmonary alveolar proteinosis (PAP) is a rare pulmonary disease characterized by abnormal accumulation of pulmonary surfactant lipids in alveoli or terminal bronchioles, leading to increased infection risk and progressive respiratory failure. Approximately more than 90% of all cases are autoimmune PAP (aPAP). Since one of the predisposing factors has been identified as genes located within the major-histocompatibility-complex region, an investigation of human leukocyte antigen (HLA) alleles associated with the risk of aPAP is warranted.

METHODS

We retrospectively studied 60 patients pathologically diagnosed with PAP from 2019 to 2022. Patients were divided into the aPAP group or secondary PAP (sPAP) group according to their clinical information. Qualified DNA was extracted from the paraffin-embedded tissue of 28 patients, and the PCR-sequence-based typing method was used for HLA-DRB1 genotyping.

RESULTS

A similar HLA-DRB1 allele profile (including the HLA-DRB1*08:03) between the aPAP group and sPAP group was revealed, except that HLA-DRB1*14:54, which has never been reported in aPAP patients, was only detected in the aPAP group rather than the sPAP group (19.4% vs. 0.0%, p = 0.030). Under inhaled granulocyte-macrophage colony-stimulating factor therapy, more clinical remission was observed in HLA-DRB1*14:54 carriers rather than in HLA-DRB1*08:03 carriers (80.0% vs. 57.1%).

CONCLUSIONS

Our real-world study revealed for the first time that a population with HLA-DRB1*14:54 was subject to aPAP, and HLA-DRB1*14:54 might imply a response in aPAP patients to inhaled granulocyte-macrophage colony-stimulating factor in aPAP patients.

Genetics of SLE: mechanistic insights from monogenic disease and disease-associated variants

The past few years have provided important insights into the genetic architecture of systemic autoimmunity through aggregation of findings from genome-wide association studies (GWAS) and whole-exome or whole-genome sequencing studies. In the prototypic systemic autoimmune disease systemic lupus erythematosus (SLE), monogenic disease accounts for a small fraction of cases but has been instrumental in the elucidation of disease mechanisms. Defects in the clearance or digestion of extracellular or intracellular DNA or RNA lead to increased sensing of nucleic acids, which can break B cell tolerance and induce the production of type I interferons leading to tissue damage. Current data suggest that multiple GWAS SLE risk alleles act in concert with rare functional variants to promote SLE development. Moreover, introduction of orthologous variant alleles into mice has revealed that pathogenic X-linked dominant and recessive SLE can be caused by novel variants in TLR7 and SAT1, respectively. Such bespoke models of disease help to unravel pathogenic pathways and can be used to test targeted therapies. Cell type-specific expression data revealed that most GWAS SLE risk genes are highly expressed in age-associated B cells (ABCs), which supports the view that ABCs produce lupus autoantibodies and contribute to end-organ damage by persisting in inflamed tissues, including the kidneys. ABCs have thus emerged as key targets of promising precision therapeutics.

Identification of novel hypermethylated or hypomethylated CpG sites and genes associated with anthracycline-induced cardiomyopathy

Anthracycline-induced cardiomyopathy is a leading cause of late morbidity in childhood cancer survivors. Aberrant DNA methylation plays a role in de novo cardiovascular disease. Epigenetic processes could play a role in anthracycline-induced cardiomyopathy but remain unstudied. We sought to examine if genome-wide differential methylation at 'CpG' sites in peripheral blood DNA is associated with anthracycline-induced cardiomyopathy. This report used participants from a matched case-control study; 52 non-Hispanic White, anthracycline-exposed childhood cancer survivors with cardiomyopathy were matched 1:1 with 52 survivors with no cardiomyopathy. Paired ChAMP (Chip Analysis Methylation Pipeline) with integrated reference-based deconvolution of adult peripheral blood DNA methylation was used to analyze data from Illumina HumanMethylation EPIC BeadChip arrays. An epigenome-wide association study (EWAS) was performed, and the model was adjusted for GrimAge, sex, interaction terms of age at enrollment, chest radiation, age at diagnosis squared, and cardiovascular risk factors (CVRFs: diabetes, hypertension, dyslipidemia). Prioritized genes were functionally validated by gene knockout in human induced pluripotent stem cell cardiomyocytes (hiPSC-CMs) using CRISPR/Cas9 technology. DNA-methylation EPIC array analyses identified 32 differentially methylated probes (DMP: 15 hyper-methylated and 17 hypo-methylated probes) that overlap with 23 genes and 9 intergenic regions. Three hundred and fifty-four differential methylated regions (DMRs) were also identified. Several of these genes are associated with cardiac dysfunction. Knockout of genes EXO6CB, FCHSD2, NIPAL2, and SYNPO2 in hiPSC-CMs increased sensitivity to doxorubicin. In addition, EWAS analysis identified hypo-methylation of probe 'cg15939386' in gene RORA to be significantly associated with anthracycline-induced cardiomyopathy. In this genome-wide DNA methylation profile study, we observed significant differences in DNA methylation at the CpG level between anthracycline-exposed childhood cancer survivors with and without cardiomyopathy, implicating differential DNA methylation of certain genes could play a role in pathogenesis of anthracycline-induced cardiomyopathy.

Molecular pathways identified from single nucleotide polymorphisms demonstrate mechanistic differences in systemic lupus erythematosus patients of Asian and European ancestry

Systemic lupus erythematosus (SLE) is a multi-organ autoimmune disorder with a prominent genetic component. Individuals of Asian-Ancestry (AsA) disproportionately experience more severe SLE compared to individuals of European-Ancestry (EA), including increased renal involvement and tissue damage. However, the mechanisms underlying elevated severity in the AsA population remain unclear. Here, we utilized available gene expression data and genotype data based on all non-HLA SNP associations in EA and AsA SLE patients detected using the Immunochip genotyping array. We identified 2778 ancestry-specific and 327 trans-ancestry SLE-risk polymorphisms. Genetic associations were examined using connectivity mapping and gene signatures based on predicted biological pathways and were used to interrogate gene expression datasets. SLE-associated pathways in AsA patients included elevated oxidative stress, altered metabolism and mitochondrial dysfunction, whereas SLE-associated pathways in EA patients included a robust interferon response (type I and II) related to enhanced cytosolic nucleic acid sensing and signaling. An independent dataset derived from summary genome-wide association data in an AsA cohort was interrogated and identified similar molecular pathways. Finally, gene expression data from AsA SLE patients corroborated the molecular pathways predicted by SNP associations. Identifying ancestry-related molecular pathways predicted by genetic SLE risk may help to disentangle the population differences in clinical severity that impact AsA and EA individuals with SLE.

The role of lysosomes in metabolic and autoimmune diseases

Lysosomes are catabolic organelles that contribute to the degradation of intracellular constituents through autophagy and of extracellular components through endocytosis, phagocytosis and macropinocytosis. They also have roles in secretory mechanisms, the generation of extracellular vesicles and certain cell death pathways. These functions make lysosomes central organelles in cell homeostasis, metabolic regulation and responses to environment changes including nutrient stresses, endoplasmic reticulum stress and defects in proteostasis. Lysosomes also have important roles in inflammation, antigen presentation and the maintenance of long-lived immune cells. Their functions are tightly regulated by transcriptional modulation via TFEB and TFE3, as well as by major signalling pathways that lead to activation of mTORC1 and mTORC2, lysosome motility and fusion with other compartments. Lysosome dysfunction and alterations in autophagy processes have been identified in a wide variety of diseases, including autoimmune, metabolic and kidney diseases. Deregulation of autophagy can contribute to inflammation, and lysosomal defects in immune cells and/or kidney cells have been reported in inflammatory and autoimmune pathologies with kidney involvement. Defects in lysosomal activity have also been identified in several pathologies with disturbances in proteostasis, including autoimmune and metabolic diseases such as Parkinson disease, diabetes mellitus and lysosomal storage diseases. Targeting lysosomes is therefore a potential therapeutic strategy to regulate inflammation and metabolism in a variety of pathologies.

Mendelian randomization and pathway analysis demonstrate shared genetic associations between lupus and coronary artery disease

Coronary artery disease (CAD) is a leading cause of death in patients with systemic lupus erythematosus (SLE). Despite clinical evidence supporting an association between SLE and CAD, pleiotropy-adjusted genetic association studies are limited and focus on only a few common risk loci. Here, we identify a net positive causal estimate of SLE-associated non-HLA SNPs on CAD by traditional Mendelian randomization (MR) approaches. Pathway analysis using SNP-to-gene mapping followed by unsupervised clustering based on protein-protein interactions (PPIs) identifies biological networks composed of positive and negative causal sets of genes. In addition, we confirm the casual effects of specific SNP-to-gene modules on CAD using only SNP mapping to each PPI-defined functional gene set as instrumental variables. This PPI-based MR approach elucidates various molecular pathways with causal implications between SLE and CAD and identifies biological pathways likely causative of both pathologies, revealing known and novel therapeutic interventions for managing CAD in SLE.

The role of ATG16L2 in autophagy and disease

Macroautophagy/autophagy, a fundamental cell process for nutrient recycling and defense against pathogens (termed xenophagy), is crucial to human health. ATG16L2 (autophagy related 16 like 2) is an autophagic protein and a paralog of ATG16L1. Both proteins are implicated in similar diseases such as cancer and other chronic diseases; however, most autophagy studies to date have primarily focused on the function of ATG16L1, with ATG16L2 remaining uncharacterized and understudied. Overexpression of ATG16L2 has been reported in various cancers including colorectal, gastric, and prostate carcinomas, whereas altered methylation of ATG16L2 has been associated with lung cancer formation and poorer response to therapy in leukemia. In addition, ATG16L2 polymorphisms have been implicated in a range of other diseases including inflammatory bowel diseases and neurodegenerative disorders. Despite this likely role in human health, the function of this enigmatic protein in autophagy remains unknown. Here, we review current studies on ATG16L2 and collate evidence that suggests that this protein is a potential modulator of autophagy as well as the implications this has on pathogenesis.Abbreviations: ATG5: autophagy related 5; ATG12: autophagy related 12; ATG16L1: autophagy related 16 like 1; ATG16L2: autophagy related 16 like 2; CD: Crohn disease; IBD: inflammatory bowel diseases; IRGM: immunity related GTPase M; MAP1LC3/LC3: microtubule associated protein 1 light chain 3; PE: phosphatidylethanolamine; RB1CC1: RB1 inducible coiled-coil 1; SLE: systemic lupus erythematosus; WIPI2B: WD repeat domain, phosphoinositide interacting 2B.

Systematic molecular profiling of acute leukemia cancer stem cells allows identification of druggable targets

Acute myeloid leukemia (AML) is one of the most prevalent and acute blood cancers with a poor prognosis and low overall survival rate, especially in the elderly. Although several new AML markers and drug targets have been recently identified, the rate of long-term cancer eradication has not improved significantly due to the presence and drug resistance of AML cancer stem cells (CSCs). Here we develop a novel computational pipeline to analyze the transcriptomic profiles of AML cancer (stem) cells and identify novel candidate AML CSC markers and drug targets. In our novel pipeline we apply a top-down meta-analysis strategy to integrate The Cancer Genome Atlas data with CSC datasets to infer cell stemness features. As a result, a set of genes termed the "AML key CSC genes" along with all the available drugs/compounds that could target them were identified. Overall, our novel computational pipeline could retrieve known cancer drugs (Carfilzomib) and predicted novel drugs such as Zonisamide, Amitriptyline, and their targets amongst the top ranked drugs and drug targets for targeting AML. Additionally, the pipeline applied in this study could be used for the identification of CSC-specific markers, drivers and their respective targeting drugs in other cancer types.

Early DNA methylation changes in children developing beta cell autoimmunity at a young age

AIMS/HYPOTHESIS

Type 1 diabetes is a chronic autoimmune disease of complex aetiology, including a potential role for epigenetic regulation. Previous epigenomic studies focused mainly on clinically diagnosed individuals. The aim of the study was to assess early DNA methylation changes associated with type 1 diabetes already before the diagnosis or even before the appearance of autoantibodies.

METHODS

Reduced representation bisulphite sequencing (RRBS) was applied to study DNA methylation in purified CD4⁺ T cell, CD8⁺ T cell and CD4^-CD8^- cell fractions of 226 peripheral blood mononuclear cell samples longitudinally collected from seven type 1 diabetes-specific autoantibody-positive individuals and control individuals matched for age, sex, HLA risk and place of birth. We also explored correlations between DNA methylation and gene expression using RNA sequencing data from the same samples. Technical validation of RRBS results was performed using pyrosequencing.

RESULTS

We identified 79, 56 and 45 differentially methylated regions in CD4⁺ T cells, CD8⁺ T cells and CD4^-CD8^- cell fractions, respectively, between type 1 diabetes-specific autoantibody-positive individuals and control participants. The analysis of pre-seroconversion samples identified DNA methylation signatures at the very early stage of disease, including differential methylation at the promoter of IRF5 in CD4⁺ T cells. Further, we validated RRBS results using pyrosequencing at the following CpG sites: chr19:18118304 in the promoter of ARRDC2; chr21:47307815 in the intron of PCBP3; and chr14:81128398 in the intergenic region near TRAF3 in CD4⁺ T cells.

CONCLUSIONS/INTERPRETATION

These preliminary results provide novel insights into cell type-specific differential epigenetic regulation of genes, which may contribute to type 1 diabetes pathogenesis at the very early stage of disease development. Should these findings be validated, they may serve as a potential signature useful for disease prediction and management.

Nucleic Acid-Sensing and Interferon-Inducible Pathways Show Differential Methylation in MZ Twins Discordant for Lupus and Overexpression in Independent Lupus Samples: Implications for Pathogenic Mechanism and Drug Targeting

Systemic lupus erythematosus (SLE) is a chronic, multisystem, autoimmune inflammatory disease with genomic and non-genomic contributions to risk. We hypothesize that epigenetic factors are a significant contributor to SLE risk and may be informative for identifying pathogenic mechanisms and therapeutic targets. To test this hypothesis while controlling for genetic background, we performed an epigenome-wide analysis of DNA methylation in genomic DNA from whole blood in three pairs of female monozygotic (MZ) twins of European ancestry, discordant for SLE. Results were replicated on the same array in four cell types from a set of four Danish female MZ twin pairs discordant for SLE. Genes implicated by the epigenetic analyses were then evaluated in 10 independent SLE gene expression datasets from the Gene Expression Omnibus (GEO). There were 59 differentially methylated loci between unaffected and affected MZ twins in whole blood, including 11 novel loci. All but two of these loci were hypomethylated in the SLE twins relative to the unaffected twins. The genes harboring these hypomethylated loci exhibited increased expression in multiple independent datasets of SLE patients. This pattern was largely consistent regardless of disease activity, cell type, or renal tissue type. The genes proximal to CpGs exhibiting differential methylation (DM) in the SLE-discordant MZ twins and exhibiting differential expression (DE) in independent SLE GEO cohorts (DM-DE genes) clustered into two pathways: the nucleic acid-sensing pathway and the type I interferon pathway. The DM-DE genes were also informatically queried for potential gene–drug interactions, yielding a list of 41 drugs including a known SLE therapy. The DM-DE genes delineate two important biologic pathways that are not only reflective of the heterogeneity of SLE but may also correlate with distinct IFN responses that depend on the source, type, and location of nucleic acid molecules and the activated receptors in individual patients. Cell- and tissue-specific analyses will be critical to the understanding of genetic factors dysregulating the nucleic acid-sensing and IFN pathways and whether these factors could be appropriate targets for therapeutic intervention.

Multi-ancestral GWAS identifies shared and Asian-specific loci for SLE and links type III interferon signaling and lysosomal function to the disease

OBJECTIVES

Systemic lupus erythematosus (SLE) is a prototype autoimmune disease with differences in prevalence and severity among ancestral groups. This study aims to identify novel genetic components either shared or distinct between Asian and European populations.

METHODS

Both trans-ancestral and ancestry-specific meta-analyses of genome-wide association studies (GWAS) for SLE were performed, involving 30,604 participants of European, Chinese or Thai origin. Using public epigenomic data and expression quantitative trait loci, fine-mapping analyses were conducted to identify putative causal variants and genes for the newly identified loci. Performance of polygenic risk scores (PRS) for the Thai cohort was evaluated comparing different training data.

RESULTS

We identified ten novel SLE susceptibility loci, four of which were found by Asian-specific meta-analyses. A 1bp deletion upstream of IFNLR1 was found associated with SLE, with the risk allele correlated with increased expression of IFNLR1. This gene encodes interferon lambda receptor 1, pointing to the role of type III interferon signaling in SLE. An intronic variant in SLC29A3 was found associated with SLE only in Asians. The putative risk variant may modulate SLC29A3 expression in a monocyte-specific manner. SLC29A3 encodes a lysosomal nucleoside transporter, and subsequent analyses suggest reduced lysosomal function and phagocytosis might be the mechanism underlying this association. In addition, trans-ancestral meta-analysis was proved to be valuable in risk prediction for individuals without ancestry-matched data.

CONCLUSION

Multi-ancestral GWAS identified both shared and Asian-specific loci for SLE, and functional annotation pointed to the involvement of increased type III interferon signaling and reduced lysosomal function in SLE.

Recent advances in understanding the genetic basis of systemic lupus erythematosus

Systemic lupus erythematosus (SLE) is a polygenic chronic autoimmune disease leading to multiple organ damage. A large heritability of up to 66% is estimated in SLE, with roughly 180 reported susceptibility loci that have been identified mostly by genome-wide association studies (GWASs) and account for approximately 30% of genetic heritability. A vast majority of risk variants reside in non-coding regions, which makes it quite challenging to interpret their functional implications in the SLE-affected immune system, suggesting the importance of understanding cell type-specific epigenetic regulation around SLE GWAS variants. The latest genetic studies have been highly fruitful as several dozens of SLE loci were newly discovered in the last few years and many loci have come to be understood in systemic approaches integrating GWAS signals with other biological resources. In this review, we summarize SLE-associated genetic variants in both the major histocompatibility complex (MHC) and non-MHC loci, examining polygenetic risk scores for SLE and their associations with clinical features. Finally, variant-driven pathogenetic functions underlying genetic associations are described, coupled with discussion about challenges and future directions in genetic studies on SLE.

Prognostic Impact of Genetic Variants of MECP2 and TIRAP on Clinical Outcomes of Systemic Lupus Erythematosus with and without Nephritis

Systemic lupus erythematosus (SLE) is a chronic autoimmune illness with a growing prevalence in many populations. Few studies have examined genetic predisposition to SLE, so we aimed to examine the clinical impact of the genetic polymorphisms MECP2 rs2734647and TIRAP rs8177374 on the outcomes and therapeutic precision of SLE with and without nephritis. This study included 110 SLE patients-divided into 63 with lupus nephritis (LN), and 47 without nephritis-and 100 controls. Laboratory measurements including CRP, ESR, ACR, CBC, anti-ds-DNA, vitamin A, C3, and C4 were carried out, along with genotyping of MECP2 rs2734647and TIRAP rs8177374 by real-time PCR and sequencing. Treg %, vitamin A, C3, and C4 were lower, whereas Th17 % was higher, in patients vs. controls (p < 0.001). The T allele of MECP2 rs2734647 was higher in LN than in non-nephritis and control subjects. Moreover, the T allele of TIRAP rs8177374 was higher in LN than in non-nephritis and control subjects. The MECP2 and TIRAP genes could play a role in predisposition to SLE, and can also predict disease progress to nephritis, helping to personalize medicine.

RNA sequencing of blood in coronary artery disease: involvement of regulatory T cell imbalance

Background

Cardiovascular disease had a global prevalence of 523 million cases and 18.6 million deaths in 2019. The current standard for diagnosing coronary artery disease (CAD) is coronary angiography. Surprisingly, despite well-established clinical indications, up to 40% of the one million invasive cardiac catheterizations return a result of ‘no blockage’. The present studies employed RNA sequencing of whole blood to identify an RNA signature in patients with angiographically confirmed CAD.

Methods

Whole blood RNA was depleted of ribosomal RNA (rRNA) and analyzed by single-molecule sequencing of RNA (RNAseq) to identify transcripts associated with CAD (TRACs) in a discovery group of 96 patients presenting for elective coronary catheterization. The resulting transcript counts were compared between groups to identify differentially expressed genes (DEGs).

Results

Surprisingly, 98% of DEGs/TRACs were down-regulated ~ 1.7-fold in patients with mild to severe CAD (> 20% stenosis). The TRACs were independent of comorbid risk factors for CAD, such as sex, hypertension, and smoking. Bioinformatic analysis identified an enrichment in transcripts such as FoxP1, ICOSLG, IKZF4/Eos, SMYD3, TRIM28, and TCF3/E2A that are likely markers of regulatory T cells (Treg), consistent with known reductions in Tregs in CAD. A validation cohort of 80 patients confirmed the overall pattern (92% down-regulation) and supported many of the Treg-related changes. TRACs were enriched for transcripts associated with stress granules, which sequester RNAs, and ciliary and synaptic transcripts, possibly consistent with changes in the immune synapse of developing T cells.

Conclusions

These studies identify a novel mRNA signature of a Treg-like defect in CAD patients and provides a blueprint for a diagnostic test for CAD. The pattern of changes is consistent with stress-related changes in the maturation of T and Treg cells, possibly due to changes in the immune synapse.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12920-021-01062-2.

Lack of WDFY4 Aggravates Ovalbumin-Induced Asthma via Enhanced Th2 Cell Differentiation

BACKGROUND

Asthma is a chronic inflammatory airway disease, and Th2 cells play an important role in asthma. WDFY4 (WDFY family member 4) is a susceptibility gene in several autoimmune diseases.

OBJECTIVE

In this study, the roles of WDFY4 in Th2 cell differentiation and Th2-dependent asthma were investigated.

METHODS

Naïve CD4+ T cells were isolated from wild-type and WDFY4-deficient mice and induced to differentiate in vitro. Subsequently, a mouse model of asthma was established by sensitization with ovalbumin.

RESULTS

Study results showed that WDFY4 deficiency could promote the differentiation of Th2 cells and the production of Th2 cytokines. WDFY4-deficient asthmatic mice showed higher levels of Th2 cytokines in the lungs and bronchoalveolar lavage fluid than wild-type mice. Moreover, infiltration of inflammatory cells, hyperplasia of goblet cells, production of mucus, and deposition of collagen were enhanced in WDFY4-deficient asthmatic mice.

CONCLUSIONS

Our study demonstrates the pivotal role of WDFY4 in the pathogenesis of asthma and in Th2 cell differentiation.

Understanding the genetics of systemic lupus erythematosus using Bayesian statistics and gene network analysis

The publication of genetic epidemiology meta-analyses has increased rapidly, but it has been suggested that many of the statistically significant results are false positive. In addition, most such meta-analyses have been redundant, duplicate, and erroneous, leading to research waste. In addition, since most claimed candidate gene associations were false-positives, correctly interpreting the published results is important. In this review, we emphasize the importance of interpreting the results of genetic epidemiology meta-analyses using Bayesian statistics and gene network analysis, which could be applied in other diseases.

Genetic determinants of risk in autoimmune pulmonary alveolar proteinosis

Pulmonary alveolar proteinosis (PAP) is a devastating lung disease caused by abnormal surfactant homeostasis, with a prevalence of 6-7 cases per million population worldwide. While mutations causing hereditary PAP have been reported, the genetic basis contributing to autoimmune PAP (aPAP) has not been thoroughly investigated. Here, we conducted a genome-wide association study of aPAP in 198 patients and 395 control participants of Japanese ancestry. The common genetic variant, rs138024423 at 6p21, in the major-histocompatibility-complex (MHC) region was significantly associated with disease risk (Odds ratio [OR] = 5.2; P = 2.4 × 10-12). HLA fine-mapping revealed that the common HLA class II allele, HLA-DRB1*08:03, strongly drove this signal (OR = 4.8; P = 4.8 × 10-12), followed by an additional independent risk allele at HLA-DPβ1 amino acid position 8 (OR = 0.28; P = 3.4 × 10-7). HLA-DRB1*08:03 was also associated with an increased level of anti-GM-CSF antibody, a key driver of the disease (β = 0.32; P = 0.035). Our study demonstrated a heritable component of aPAP, suggesting an underlying genetic predisposition toward an abnormal antibody production.

Trans-Ancestral Fine-Mapping and Epigenetic Annotation as Tools to Delineate Functionally Relevant Risk Alleles at IKZF1 and IKZF3 in Systemic Lupus Erythematosus

Background: Prioritizing tag-SNPs carried on extended risk haplotypes at susceptibility loci for common disease is a challenge. Methods: We utilized trans-ancestral exclusion mapping to reduce risk haplotypes at IKZF1 and IKZF3 identified in multiple ancestries from SLE GWAS and ImmunoChip datasets. We characterized functional annotation data across each risk haplotype from publicly available datasets including ENCODE, RoadMap Consortium, PC Hi-C data from 3D genome browser, NESDR NTR conditional eQTL database, GeneCards Genehancers and TF (transcription factor) binding sites from Haploregv4. Results: We refined the 60 kb associated haplotype upstream of IKZF1 to just 12 tag-SNPs tagging a 47.7 kb core risk haplotype. There was preferential enrichment of DNAse I hypersensitivity and H3K27ac modification across the 3′ end of the risk haplotype, with four tag-SNPs sharing allele-specific TF binding sites with promoter variants, which are eQTLs for IKZF1 in whole blood. At IKZF3, we refined a core risk haplotype of 101 kb (27 tag-SNPs) from an initial extended haplotype of 194 kb (282 tag-SNPs), which had widespread DNAse I hypersensitivity, H3K27ac modification and multiple allele-specific TF binding sites. Dimerization of Fox family TFs bound at the 3′ and promoter of IKZF3 may stabilize chromatin looping across the locus. Conclusions: We combined trans-ancestral exclusion mapping and epigenetic annotation to identify variants at both IKZF1 and IKZF3 with the highest likelihood of biological relevance. The approach will be of strong interest to other complex trait geneticists seeking to attribute biological relevance to risk alleles on extended risk haplotypes in their disease of interest.

Analysis of Trans-Ancestral SLE Risk Loci Identifies Unique Biologic Networks and Drug Targets in African and European Ancestries

Systemic lupus erythematosus (SLE) is a multi-organ autoimmune disorder with a prominent genetic component. Individuals of African ancestry (AA) experience the disease more severely and with an increased co-morbidity burden compared to European ancestry (EA) populations. We hypothesize that the disparities in disease prevalence, activity, and response to standard medications between AA and EA populations is partially conferred by genomic influences on biological pathways. To address this, we applied a comprehensive approach to identify all genes predicted from SNP-associated risk loci detected with the Immunochip. By combining genes predicted via eQTL analysis, as well as those predicted from base-pair changes in intergenic enhancer sites, coding-region variants, and SNP-gene proximity, we were able to identify 1,731 potential ancestry-specific and trans-ancestry genetic drivers of SLE. Gene associations were linked to upstream and downstream regulators using connectivity mapping, and predicted biological pathways were mined for candidate drug targets. Examination of trans-ancestral pathways reflect the well-defined role for interferons in SLE and revealed pathways associated with tissue repair and remodeling. EA-dominant genetic drivers were more often associated with innate immune and myeloid cell function pathways, whereas AA-dominant pathways mirror clinical findings in AA subjects, suggesting disease progression is driven by aberrant B cell activity accompanied by ER stress and metabolic dysfunction. Finally, potential ancestry-specific and non-specific drug candidates were identified. The integration of all SLE SNP-predicted genes into functional pathways revealed critical molecular pathways representative of each population, underscoring the influence of ancestry on disease mechanism and also providing key insight for therapeutic selection.

Structurally distinct endocytic pathways for B cell receptors in B lymphocytes

B lymphocytes play a critical role in adaptive immunity. On antigen binding, B cell receptors (BCR) cluster on the plasma membrane and are internalized by endocytosis. In this process, B cells capture diverse antigens in various contexts and concentrations. However, it is unclear whether the mechanism of BCR endocytosis changes in response to these factors. Here, we studied the mechanism of soluble antigen-induced BCR clustering and internalization in a cultured human B cell line using correlative superresolution fluorescence and platinum replica electron microscopy. First, by visualizing nanoscale BCR clusters, we provide direct evidence that BCR cluster size increases with F(ab’)2 concentration. Next, we show that the physical mechanism of internalization switches in response to BCR cluster size. At low concentrations of antigen, B cells internalize small BCR clusters by classical clathrin-mediated endocytosis. At high antigen concentrations, when cluster size increases beyond the size of a single clathrin-coated pit, B cells retrieve receptor clusters using large invaginations of the plasma membrane capped with clathrin. At these sites, we observed early and sustained recruitment of actin and an actin polymerizing protein FCHSD2. We further show that actin recruitment is required for the efficient generation of these novel endocytic carriers and for their capture into the cytosol. We propose that in B cells, the mechanism of endocytosis switches to accommodate large receptor clusters formed when cells encounter high concentrations of soluble antigen. This mechanism is regulated by the organization and dynamics of the cortical actin cytoskeleton.

In Vivo Remodeling of Altered Autophagy-Lysosomal Pathway by a Phosphopeptide in Lupus

The phosphopeptide P140/Lupuzor, which improves the course of lupus disease in mice and patients, targets chaperone-mediated autophagy (CMA), a selective form of autophagy that is abnormally upregulated in lupus-prone MRL/lpr mice. Administered intravenously to diseased mice, P140 reduces the expression level of two major protein players of CMA, LAMP2A and HSPA8, and inhibits CMA in vitro in a cell line that stably expresses a CMA reporter. Here, we aimed to demonstrate that P140 also affects CMA in vivo and to unravel the precise cellular mechanism of how P140 interacts with the CMA process. MRL/lpr mice and CBA/J mice used as control received P140 or control peptides intravenously. Lysosome-enriched fractions of spleen or liver were prepared to examine lysosomal function. Highly purified lysosomes were further isolated and left to incubate with the CMA substrate to study at which cellular step P140 interacts with the CMA process. The data show that P140 effectively regulates CMA in vivo in MRL/lpr mice at the step of substrate lysosomal uptake and restores some alterations of defective lysosomes. For the first time, it is demonstrated that by occluding the intralysosome uptake of CMA substrates, a therapeutic molecule can attenuate excessive CMA activity in a pathological pro-inflammatory context and protect against hyperinflammation. This recovery effect of P140 on hyperactivated CMA is not only important for lupus therapy but potentially also for treating other (auto)inflammatory diseases, including neurologic and metabolic disorders, where CMA modulation would be highly beneficial.

Analysis of clinical and laboratory characteristics and pathology of lupus nephritis-based on 710 renal biopsies in China

OBJECTIVES

This study aimed to analyse the clinical and laboratory characteristics of different pathologic classifications of lupus nephritis (LN) patients in terms of age at systemic lupus erythematosus (SLE) diagnosis and nephritis onset.

METHOD

Clinical, laboratory, and pathological data of 710 LN patients diagnosed by renal biopsy at our institution between 2000 and 2018 were retrospectively analysed. Patients were divided into the different pathological classification groups; childhood-, adult- and elderly-onset SLE groups and early- and late-onset LN groups.

RESULTS

Class IV occurred most frequently and had the lowest complement C3 level. There was an obvious increase in active index in class IV and class V + IV. Patients with class VI showed some clinical characteristics similar to end-stage renal disease. Patients with proliferative nephritis were younger at SLE diagnosis and had higher blood pressure, higher frequency of proteinuria and urinary erythrocyte and lower haemoglobin and complement C3. Pathologic classification between childhood-, adult- and elderly-onset SLE patients or between early- and late-onset LN patients was not significantly different. Elderly-onset SLE patients had the highest chronic index (CI), IgA, IgG and Sjögren's syndrome A antibodies and Sjögren's syndrome B antibodies rates, whereas late-onset LN patients showed significantly higher CI, haemoglobin, complement C3 and C4 but lower uric acid, IgM and IgG.

CONCLUSIONS

LN patients present with different clinical and laboratory characteristics according to pathological classification, age at SLE diagnosis and nephritis onset. These results might be valuable for estimating the pathology and guiding treatment and prognosis. Key Points • Patients with proliferative nephritis have more severe immune disorders, worse renal function and stronger inflammatory state. • The elderly-onset SLE patients showed a poorer condition. • The late-onset LN patients might have a more stable status.

The association of the UHRF1BP1 gene with systemic lupus erythematosus was replicated in a Han Chinese population from mainland China

Single-nucleotide polymorphisms (SNPs) in the UHRF gene have been shown to be associated with systemic lupus erythematosus (SLE) in European and Hong Kong Chinese, but statistically significant evidence for association has not been found in a mainland Han Chinese population. Therefore, we selected SNP rs13205210 located in UHRF1BP1 as a candidate association from our previously published genome-wide association study (GWAS) data of SLE (1,047 cases and 1,205 controls from a mainland Han Chinese population) to explore the association between the UHRF1BP1 gene and SLE. We conducted a large-scale replication study in an additional independent sample of 3,509 cases and 8,246 controls from a mainland Han Chinese population. Real-time PCR was used to determine gene expression differences in peripheral blood mononuclear cells (PBMCs) from cases and controls. As a result, we replicated the association between the UHRF1BP1 gene and SLE (rs13205210, missense, Pmeta  = 2.26E-17, odds ratio = 1.41) by a meta-analysis of our previous GWAS and this replication study involving a total of 4,556 cases and 9,451 controls. The UHRF1BP1 mRNA expression level in PBMCs was significantly decreased in patients with SLE compared with that in healthy controls. SNP rs13205210 exhibited an expression quantitative trait loci effect on the UHRF1BP1 gene in PBMCs from patients. In conclusion, this study not only suggests that the UHRF1BP1 gene was associated with SLE in a mainland Han Chinese population, but also implied that it might be a common genetic factor contributing to SLE susceptibility in multiple populations.

An epigenome-wide association study of sex-specific chronological ageing

BACKGROUND

Advanced age is associated with cognitive and physical decline and is a major risk factor for a multitude of disorders. There is also a gap in life expectancy between males and females. DNA methylation differences have been shown to be associated with both age and sex. Here, we investigate age-by-sex differences in blood-based DNA methylation in an unrelated cohort of 2586 individuals between the ages of 18 and 87 years, with replication in a further 4450 individuals between the ages of 18 and 93 years.

METHODS

Linear regression models were applied, with stringent genome-wide significance thresholds (p < 3.6 × 10-8) used in both the discovery and replication data. A second, highly conservative mixed linear model method that better controls the false-positive rate was also applied, using the same genome-wide significance thresholds.

RESULTS

Using the linear regression method, 52 autosomal and 597 X-linked CpG sites, mapping to 251 unique genes, replicated with concordant effect size directions in the age-by-sex interaction analysis. The site with the greatest difference mapped to GAGE10, an X-linked gene. Here, DNA methylation levels remained stable across the male adult age range (DNA methylation by age r = 0.02) but decreased across female adult age range (DNA methylation by age r = - 0.61). One site (cg23722529) with a significant age-by-sex interaction also had a quantitative trait locus (rs17321482) that is a genome-wide significant variant for prostate cancer. The mixed linear model method identified 11 CpG sites associated with the age-by-sex interaction.

CONCLUSION

The majority of differences in age-associated DNA methylation trajectories between sexes are present on the X chromosome. Several of these differences occur within genes that have been implicated in sexually dimorphic traits.

Lysosomes as a therapeutic target

Lysosomes are membrane-bound organelles with roles in processes involved in degrading and recycling cellular waste, cellular signalling and energy metabolism. Defects in genes encoding lysosomal proteins cause lysosomal storage disorders, in which enzyme replacement therapy has proved successful. Growing evidence also implicates roles for lysosomal dysfunction in more common diseases including inflammatory and autoimmune disorders, neurodegenerative diseases, cancer and metabolic disorders. With a focus on lysosomal dysfunction in autoimmune disorders and neurodegenerative diseases - including lupus, rheumatoid arthritis, multiple sclerosis, Alzheimer disease and Parkinson disease - this Review critically analyses progress and opportunities for therapeutically targeting lysosomal proteins and processes, particularly with small molecules and peptide drugs.

Distinct Tissue-Specific Roles for the Disease-Associated Autophagy Genes ATG16L2 and ATG16L1

The clear role of autophagy in human inflammatory diseases such as Crohn disease was first identified by genome-wide association studies and subsequently dissected in multiple mechanistic studies. ATG16L1 has been particularly well studied in knockout and hypomorph settings as well as models recapitulating the Crohn disease-associated T300A polymorphism. Interestingly, ATG16L1 has a single homolog, ATG16L2, which is independently implicated in diseases, including Crohn disease and systemic lupus erythematosus. However, the contribution of ATG16L2 to canonical autophagy pathways and other cellular functions is poorly understood. To better understand its role, we generated and analyzed the first, to our knowledge, ATG16L2 knockout mouse. Our results show that ATG16L1 and ATG16L2 contribute very distinctly to autophagy and cellular ontogeny in myeloid, lymphoid, and epithelial lineages. Dysregulation of any of these lineages could contribute to complex diseases like Crohn disease and systemic lupus erythematosus, highlighting the value of examining cell-specific effects. We also identify a novel genetic interaction between ATG16L2 and epithelial ATG16L1. These findings are discussed in the context of how these genes may contribute distinctly to human disease.

Update on the Genetics of Systemic Lupus Erythematosus: Genome-Wide Association Studies and Beyond

Systemic lupus erythematosus (SLE) is a chronic autoimmune disease of complex etiology that primarily affects women of childbearing age. The development of SLE is attributed to the breach of immunological tolerance and the interaction between SLE-susceptibility genes and various environmental factors, resulting in the production of pathogenic autoantibodies. Working in concert with the innate and adaptive arms of the immune system, lupus-related autoantibodies mediate immune-complex deposition in various tissues and organs, leading to acute and chronic inflammation and consequent end-organ damage. Over the past two decades or so, the impact of genetic susceptibility on the development of SLE has been well demonstrated in a number of large-scale genetic association studies which have uncovered a large fraction of genetic heritability of SLE by recognizing about a hundred SLE-susceptibility loci. Integration of genetic variant data with various omics data such as transcriptomic and epigenomic data potentially provides a unique opportunity to further understand the roles of SLE risk variants in regulating the molecular phenotypes by various disease-relevant cell types and in shaping the immune systems with high inter-individual variances in disease susceptibility. In this review, the catalogue of SLE susceptibility loci will be updated, and biological signatures implicated by the SLE-risk variants will be critically discussed. It is optimistically hoped that identification of SLE risk variants will enable the prognostic and therapeutic biomarker armamentarium of SLE to be strengthened, a major leap towards precision medicine in the management of the condition.

SNP rs4937333 in the miRNA-5003-Binding Site of the ETS1 3'-UTR Decreases ETS1 Expression

Mutations in and reduced expression of the ETS1 gene may be associated with systemic lupus erythematosus (SLE). Here, we report a replication study to investigate associations of eight ETS1 single-nucleotide polymorphisms in the 3′-untranslated region (3′-UTR) with SLE and their regulation of ETS1 expression in a study population. We found that the rs4937333 T allele was associated with a significantly increased risk of SLE (odds ratio: 1.800, 95% confidence interval: 1.02–3.157, P = 0.040) and with dramatically reduced levels of ETS1 in B cells from SLE subjects. Functionally, the rs4937333 T allele alters the binding affinity between miR-5003 and its ETS1 3′-UTR target, thus enhancing suppression of ETS1 expression. Furthermore, immunoglobulin M-secreting plasmacytes were significantly reduced among B cells with the rs4937333 C allele versus the T allele according to FACS and ELISA. Additionally, miR-5003 expression was higher in B cells than in T cells from SLE patients, and a negative correlation between miR-5003 and ETS1 was found, especially in B cells with the T allele. These findings suggest that the rs4937333 T allele is a risk factor for susceptibility to SLE in the studied population. The rs4937333 T allele may enhance the binding of miR-5003 to ETS1, which probably promotes the involvement of ETS1 in the differentiation of B cells into plasmacytes.

Amino acid signatures of HLA Class-I and II molecules are strongly associated with SLE susceptibility and autoantibody production in Eastern Asians

Human leukocyte antigen (HLA) is a key genetic factor conferring risk of systemic lupus erythematosus (SLE), but precise independent localization of HLA effects is extremely challenging. As a result, the contribution of specific HLA alleles and amino-acid residues to the overall risk of SLE and to risk of specific autoantibodies are far from completely understood. Here, we dissected (a) overall SLE association signals across HLA, (b) HLA-peptide interaction, and (c) residue-autoantibody association. Classical alleles, SNPs, and amino-acid residues of eight HLA genes were imputed across 4,915 SLE cases and 13,513 controls from Eastern Asia. We performed association followed by conditional analysis across HLA, assessing both overall SLE risk and risk of autoantibody production. DR15 alleles HLA-DRB1*15:01 (P = 1.4x10^-27, odds ratio (OR) = 1.57) and HLA-DQB1*06:02 (P = 7.4x10^-23, OR = 1.55) formed the most significant haplotype (OR = 2.33). Conditioned protein-residue signals were stronger than allele signals and mapped predominantly to HLA-DRB1 residue 13 (P = 2.2x10^-75) and its proxy position 11 (P = 1.1x10^-67), followed by HLA-DRB1-37 (P = 4.5x10^-24). After conditioning on HLA-DRB1, novel associations at HLA-A-70 (P = 1.4x10^-8), HLA-DPB1-35 (P = 9.0x10^-16), HLA-DQB1-37 (P = 2.7x10^-14), and HLA-B-9 (P = 6.5x10^-15) emerged. Together, these seven residues increased the proportion of explained heritability due to HLA to 2.6%. Risk residues for both overall disease and hallmark autoantibodies (i.e., nRNP: DRB1-11, P = 2.0x10^-14; DRB1-13, P = 2.9x10^-13; DRB1-30, P = 3.9x10^-14) localized to the peptide-binding groove of HLA-DRB1. Enrichment for specific amino-acid characteristics in the peptide-binding groove correlated with overall SLE risk and with autoantibody presence. Risk residues were in primarily negatively charged side-chains, in contrast with rheumatoid arthritis. We identified novel SLE signals in HLA Class I loci (HLA-A, HLA-B), and localized primary Class II signals to five residues in HLA-DRB1, HLA-DPB1, and HLA-DQB1. These findings provide insights about the mechanisms by which the risk residues interact with each other to produce autoantibodies and are involved in SLE pathophysiology.

The Human leukocyte antigen (HLA) region is a key genetic factor conferring risk of systemic lupus erythematosus (SLE). In spite of multiple SLE association signals identified in the HLA region, only amino-acid residues within HLA-DRB1 have been specifically described previously. In this study, we performed an imputation-based analysis on individuals with East Asian ancestry, and characterized SLE risk within the HLA region for all involved independent genes (HLA-DRB1, HLA-DPB1, HLA-DQB1, HLA-A, and HLA-B). Furthermore, we identified a characteristic SLE risk residue signature as well as a pattern of specific nRNP and Ro/La autoantibody residues located in the peptide-binding grooves, suggesting their key involvement in autoantibody production.

Cross-disease analysis of Alzheimer's disease and type-2 Diabetes highlights the role of autophagy in the pathophysiology of two highly comorbid diseases

Evidence is accumulating that the main chronic diseases of aging Alzheimer's disease (AD) and type-2 diabetes mellitus (T2DM) share common pathophysiological mechanisms. This study aimed at applying systems biology approaches to increase the knowledge of the shared molecular pathways underpinnings of AD and T2DM. We analysed transcriptomic data of post-mortem AD and T2DM human brains to obtain disease signatures of AD and T2DM and combined them with protein-protein interaction information to construct two disease-specific networks. The overlapping AD/T2DM network proteins were then used to extract the most representative Gene Ontology biological process terms. The expression of genes identified as relevant was studied in two AD models, 3xTg-AD and ApoE3/ApoE4 targeted replacement mice. The present transcriptomic data analysis revealed a principal role for autophagy in the molecular basis of both AD and T2DM. Our experimental validation in mouse AD models confirmed the role of autophagy-related genes. Among modulated genes, Cyclin-Dependent Kinase Inhibitor 1B, Autophagy Related 16-Like 2, and insulin were highlighted. In conclusion, the present investigation revealed autophagy as the central dys-regulated pathway in highly co-morbid diseases such as AD and T2DM allowing the identification of specific genes potentially involved in disease pathophysiology which could become novel targets for therapeutic intervention.

Exploring Impact of Rare Variation in Systemic Lupus Erythematosus by a Genome Wide Imputation Approach

The importance of low frequency and rare variation in complex disease genetics is difficult to estimate in patient populations. Genome-wide association studies are therefore, underpowered to detect rare variation. We have used a combined approach of genome-wide-based imputation with a highly stringent sequence kernel association (SKAT) test and a case-control burden test. We identified 98 candidate genes containing rare variation that in aggregate show association with SLE many of which have recognized immunological function, but also function and expression related to relevant tissues such as the joints, skin, blood or central nervous system. In addition we also find that there is a significant enrichment of genes annotated for disease-causing mutations in the OMIM database, suggesting that in complex diseases such as SLE, such mutations may be involved in subtle or combined phenotypes or could accelerate specific organ abnormalities found in the disease. We here provide an important resource of candidate genes for SLE.

Elucidating the environmental risk factors for rheumatic diseases: An umbrella review of meta-analyses

AIMS

Although rheumatic diseases constitute a leading cause of disability, the environmental risk factors for these diseases are not clarified. In the present study, we aim to systematically appraise the epidemiological credibility of the environmental risk factors for rheumatic diseases.

METHODS

We systematically searched PubMed to capture meta-analyses of observational studies on environmental risk factors for the most prevalent rheumatic diseases. For each association, we estimated the summary effect size estimate, the 95% confidence and prediction intervals, and the I² metric. We further examined the presence of small-study effects and excess significance bias.

RESULTS

Overall, we identified 30 eligible papers describing 42 associations. Thirty-three associations were statistically significant at P < 0.05, whereas 13 of them were statistically significant at P < 1 × 10^-6 . Thirty-two associations had large or very large between-study heterogeneity. In 12 associations, evidence of small-study effects and/or excess significance bias was found. Six risk factors (nine associations) presented convincing or highly suggestive evidence of association: smoking and pack-years of smoking for rheumatoid arthritis; BMI (per 5 kg/m² increase) for gout and hip osteoarthritis; alcohol consumption for gout; BMI (overweight vs lean, obese vs lean), knee injury and participation in heavy work for knee osteoarthritis.

CONCLUSION

Our umbrella review indicated that a narrow range of risk factors has been examined for rheumatic diseases. Current evidence strongly supports that smoking, obesity, alcohol consumption, knee injury, and work activities are associated with risk for at least one rheumatic disease.

Genetic variants in systemic lupus erythematosus susceptibility loci, XKR6 and GLT1D1 are associated with childhood-onset SLE in a Korean cohort

Impact of genetic variants on the age of systemic lupus erythematosus (SLE) onset was not fully understood. We investigated a cumulative effect of SLE-risk variants on the age of SLE onset and scanned genome-wide SNPs to search for new risk loci of childhood-onset SLE (cSLE). We analyzed 781 Korean single-center SLE subjects who previously genotyped by both Immunochip and genome-wide SNP arrays. Individual genetic risk scores (GRS) from well-validated SLE susceptibility loci were calculated and tested for their association with cSLE (<16 years at onset). Single-variant association tests were performed using a multivariable logistic regression adjusting for population stratification. GRS from SLE susceptibility loci was significantly higher in cSLE than aSLE (p = 1.23 × 10⁻³). Two SNPs, rs7460469 in XKR6 (p = 1.26 × 10⁻⁸, OR = 5.58) and rs7300146 in GLT1D1 p = 1.49 × 10⁻⁸, OR = 2.85), showed the most significant associations with cSLE. The model consisting of GRS of SLE and two newly identified loci showed an area under curve (AUC) of 0.71 in a receiver operating characteristics (ROC) curve for prediction of cSLE. In conclusion, cSLE is associated with a high cumulative SLE-risk effect and two novel SNPs rs7460469 and rs7300146, providing the first predictive model for cSLE in Koreans.

A plausibly causal functional lupus-associated risk variant in the STAT1-STAT4 locus

Systemic lupus erythematosus (SLE or lupus) (OMIM: 152700) is a chronic autoimmune disease with debilitating inflammation that affects multiple organ systems. The STAT1-STAT4 locus is one of the first and most highly replicated genetic loci associated with lupus risk. We performed a fine-mapping study to identify plausible causal variants within the STAT1-STAT4 locus associated with increased lupus disease risk. Using complementary frequentist and Bayesian approaches in trans-ancestral Discovery and Replication cohorts, we found one variant whose association with lupus risk is supported across ancestries in both the Discovery and Replication cohorts: rs11889341. In B cell lines from patients with lupus and healthy controls, the lupus risk allele of rs11889341 was associated with increased STAT1 expression. We demonstrated that the transcription factor HMGA1, a member of the HMG transcription factor family with an AT-hook DNA-binding domain, has enriched binding to the risk allele compared with the non-risk allele of rs11889341. We identified a genotype-dependent repressive element in the DNA within the intron of STAT4 surrounding rs11889341. Consistent with expression quantitative trait locus (eQTL) analysis, the lupus risk allele of rs11889341 decreased the activity of this putative repressor. Altogether, we present a plausible molecular mechanism for increased lupus risk at the STAT1-STAT4 locus in which the risk allele of rs11889341, the most probable causal variant, leads to elevated STAT1 expression in B cells due to decreased repressor activity mediated by increased binding of HMGA1.

IRF5 is elevated in childhood-onset SLE and regulated by histone acetyltransferase and histone deacetylase inhibitors

Interferon regulatory factor 5 (IRF5) plays a critical role in the induction of type I interferon, proinflammatory cytokines and chemokines, and participates in the pathogenesis of autoimmune diseases such as systemic lupus erythematosus (SLE). However, the relationship between IRF5 and childhood-onset SLE remains elusive. In the present study, we demonstrated that levels of mRNA expression of IRF5, IFN-α, and Sp1 were significantly increased in childhood-onset SLE, as seen on quantitative real-time PCR, and the expression of Sp1 and IFN-α was positively correlated with IRF5. In addition to being used as antitumor drugs, a number of histone deacetylase inhibitors (HDACi) display potent anti-inflammatory properties; however, their effects on IRF5 expression remain unclear. In this study, we identified that HDACi trichostatin A (TSA) and histone acetyltransferase (HAT)-p300 downregulated IRF5 promoter activity, mRNA expression, and protein level, whereas the HAT-p300/CBP-associated factor had no effect. Moreover, TSA inhibited the production of TNF-α and IL-6 in differentiated THP-1cells. Furthermore, chromatin immunoprecipitation assays revealed that TSA inhibited DNA binding of Sp1, RNA polymerase II, HDAC3, and p300 to the core promoter region of IRF5. Our results suggest that HDACi may have therapeutic potential in patients with autoimmune diseases such as SLE through repression of IRF5 expression.

Potentially Functional Variants of ATG16L2 Predict Radiation Pneumonitis and Outcomes in Patients with Non-Small Cell Lung Cancer after Definitive Radiotherapy

INTRODUCTION

Autophagy not only plays an important role in the progression of cancer but is also involved in tissue inflammatory response. However, few published studies have investigated associations between functional genetic variants of autophagy-related genes and radiation pneumonitis (RP) as well as clinical outcomes in patients with NSCLC after definitive radiotherapy.

METHODS

We genotyped nine potentially functional single-nucleotide polymorphisms (SNPs) in four autophagy-related genes (autophagy related 2B gene [ATG2B], autophagy related 10 gene [ATG10], autophagy related 12 gene [ATG12], and autophagy related 16 like 2 gene [ATG16L2]) in 393 North American patients with NSCLC treated by definitive radiotherapy and assessed their associations with RP, local recurrence-free survival (LRFS), progression-free survival (PFS), and overall survival (OS) in multivariable Cox proportional hazard regression analyses.

RESULTS

We found that patients with the ATG16L2 rs10898880 CC variant genotype had a better LRFS, PFS, and OS (adjusted hazard ratio = 0.59, 0.64, and 0.64; 95% confidence interval: 0.45-0.79, 0.48-0.84, and 0.48-0.86; p = 0.0004, 0.002, and 0.003, respectively), but a greater risk for development of severe RP (adjusted hazard ratio = 1.80, 95% confidence interval: 1.04-3.12, p = 0.037) than did patients with AA/AC genotypes. Further functional analyses suggested that the ATG16L2 rs10898880 C variant allele modulated expression of the ATG16L2 gene.

CONCLUSION

This is the first report that one potentially functional SNP rs10898880 in ATG16L2 may be a predictor of RP, LRFS, PFS, and OS in patients with NSCLC after definitive radiotherapy. Additional larger, prospective studies are needed to confirm these findings.

Updates in Lupus Genetics

PURPOSE OF REVIEW

Our understanding on genetic basis of SLE has been advanced through genome-wide association studies. We review recent progress in lupus genetics with a focus on SLE-associated loci that have been functionally characterized, and discuss the potential for clinical translation of genetics data.

RECENT FINDINGS

Over 100 loci have been confirmed to show robust association with SLE and many share with other immune-mediated diseases. Although causative variants captured at these established loci are limited, they guide biological studies of gene targets for functional characterization which highlight the importance of aberrant recognition of self-nucleic acid, type I interferon overproduction, and defective immune cell signaling underlying the pathogenesis of SLE. Increasing examples illustrate a predictive value of genetic findings in susceptibility/prognosis prediction, clinical classification, and pharmacological implication. Genetic findings provide a foundation for better understanding of disease pathogenic mechanisms and opportunities for target selection in lupus drug development.

Genetics of systemic lupus erythematosus and Sjögren's syndrome: an update

PURPOSE OF REVIEW

To describe the recent studies on the genetics of systemic lupus erythematosus (SLE) and Sjögren's syndrome.

RECENT FINDINGS

We overview the most recent findings on the genetic susceptibility of the diseases and provide information on their genetic similarities and differences.

SUMMARY

SLE and Sjögren's syndrome are two closely related systemic autoimmune diseases that share multiple clinical and molecular aspects, including a significant number of susceptibility genes. Several genome-wide association studies were recently published in different populations that provide a better picture of their molecular mechanisms. It is becoming clear that their genetic architecture is quite well established, but more information is required on expression quantitative trait loci, epigenetic genome-wide analyses, gene × gene interactions and the role of rare variants.

Multiple variants in 5q31.1 are associated with systemic lupus erythematosus susceptibility and subphenotypes in the Han Chinese population

BACKGROUND

A previous study provided evidence for a genetic association between PPP2CA on 5q31.1 and systemic lupus erythematosus (SLE) across multi-ancestral cohorts, but failed to find significant evidence for an association in the Han Chinese population.

OBJECTIVES

To explore the association between this locus and SLE using data from our previously published genome-wide association study (GWAS).

METHODS

Single-nucleotide polymorphisms (SNPs) rs7726414 and rs244689 (near TCF7 and PPP2CA in 5q31.1) were selected as candidate independent associations from a large-scale study in a Han Chinese population consisting of 1047 cases and 1205 controls. Subsequently, 3509 cases and 8246 controls were genotyped in two further replication studies. We then investigated the SNPs' associations with SLE subphenotypes and gene expression in peripheral blood mononuclear cells.

RESULTS

Highly significant associations with SLE in the Han Chinese population were detected for SNPs rs7726414 and rs244689 by combining the genotype data from our previous GWAS and two independent replication cohorts. Further conditional analyses indicated that these two SNPs contribute to disease susceptibility independently. A significant association with SLE, age at diagnosis < 20 years, was found for rs7726414 (P = 0·001). The expression levels of TCF7 and PPP2CA messenger RNA in patients with SLE were significantly decreased compared with those in healthy controls.

CONCLUSIONS

This study found evidence for multiple associations with SLE in 5q31.1 at genome-wide levels of significance for the first time in a Han Chinese population, in a combined genotype dataset. These findings suggest that variants in the 5q31.1 locus not only provide novel insights into the genetic architecture of SLE, but also contribute to the complex subphenotypes of SLE.

Confirmation of five novel susceptibility loci for systemic lupus erythematosus (SLE) and integrated network analysis of 82 SLE susceptibility loci

We recently identified ten novel SLE susceptibility loci in Asians and uncovered several additional suggestive loci requiring further validation. This study aimed to replicate five of these suggestive loci in a Han Chinese cohort from Hong Kong, followed by meta-analysis (11,656 cases and 23,968 controls) on previously reported Asian and European populations, and to perform bioinformatic analyses on all 82 reported SLE loci to identify shared regulatory signatures. We performed a battery of analyses for these five loci, as well as joint analyses on all 82 SLE loci. All five loci passed genome-wide significance: MYNN (rs10936599, Pmeta = 1.92 × 10-13, OR = 1.14), ATG16L2 (rs11235604, Pmeta = 8.87 × 10 -12, OR = 0.78), CCL22 (rs223881, Pmeta = 5.87 × 10-16, OR = 0.87), ANKS1A (rs2762340, Pmeta = 4.93 × 10-15, OR = 0.87) and RNASEH2C (rs1308020, Pmeta = 2.96 × 10-19, OR = 0.84) and co-located with annotated gene regulatory elements. The novel loci share genetic signatures with other reported SLE loci, including effects on gene expression, transcription factor binding, and epigenetic characteristics. Most (56%) of the correlated (r2 > 0.8) SNPs from the 82 SLE loci were implicated in differential expression (9.81 × 10-198 < P < 5 × 10-3) of cis-genes. Transcription factor binding sites for p53, MEF2A and E2F1 were significantly (P < 0.05) over-represented in SLE loci, consistent with apoptosis playing a critical role in SLE. Enrichment analysis revealed common pathways, gene ontology, protein domains, and cell type-specific expression. In summary, we provide evidence of five novel SLE susceptibility loci. Integrated bioinformatics using all 82 loci revealed that SLE susceptibility loci share many gene regulatory features, suggestive of conserved mechanisms of SLE etiopathogenesis.

Recapitulation of Candidate Systemic Lupus Erythematosus-Associated Variants in Koreans

Systemic lupus erythematosus (SLE) is a chronic autoimmune disease that affects multiple organ systems. Although the etiology of SLE remains unclear, it is widely accepted that genetic factors could be involved in its pathogenesis. A number of genome-wide association studies (GWASs) have identified novel single-nucleotide polymorphisms (SNPs) associated with the risk of SLE in diverse populations. However, not all the SNP candidates identified from non-Asian populations have been validated in Koreans. In this study, we aimed to replicate the SNPs that were recently discovered in the GWAS; these SNPs have not been validated in Koreans or have only been replicated in Koreans with an insufficient sample size to conclude any association. For this, we selected five SNPs (rs1801274 in FCGR2A and rs2286672 in PLD2, rs887369 in CXorf21, rs9782955 in LYST, and rs3794060 in NADSYN1). Through the replication study with 656 cases and 622 controls, rs1801274 in FCGR2A was found to be significantly associated with SLE in Koreans (odds ratio, 1.26, 95% confidence interval, 1.06 to 1.50; p = 0.01 in allelic model). This association was also significant in two other models (dominant and recessive). The other four SNPs did not show a significant association. Our data support that FCGR polymorphisms play important roles in the susceptibility to SLE in diverse populations, including Koreans.

Genome-Wide Association Study in an Amerindian Ancestry Population Reveals Novel Systemic Lupus Erythematosus Risk Loci and the Role of European Admixture

OBJECTIVE

Systemic lupus erythematosus (SLE) is a chronic autoimmune disease with a strong genetic component. We undertook the present work to perform the first genome-wide association study on individuals from the Americas who are enriched for Native American heritage.

METHODS

We analyzed 3,710 individuals from the US and 4 countries of Latin America who were diagnosed as having SLE, and healthy controls. Samples were genotyped with HumanOmni1 BeadChip. Data on out-of-study controls genotyped with HumanOmni2.5 were also included. Statistical analyses were performed using SNPtest and SNPGWA. Data were adjusted for genomic control and false discovery rate. Imputation was performed using Impute2 and, for classic HLA alleles, HiBag. Odds ratios (ORs) and 95% confidence intervals (95% CIs) were calculated.

RESULTS

The IRF5-TNPO3 region showed the strongest association and largest OR for SLE (rs10488631: genomic control-adjusted P [Pgcadj ] = 2.61 × 10(-29), OR 2.12 [95% CI 1.88-2.39]), followed by HLA class II on the DQA2-DQB1 loci (rs9275572: Pgcadj  = 1.11 × 10(-16), OR 1.62 [95% CI 1.46-1.80] and rs9271366: Pgcadj  = 6.46 × 10(-12), OR 2.06 [95% CI 1.71-2.50]). Other known SLE loci found to be associated in this population were ITGAM, STAT4, TNIP1, NCF2, and IRAK1. We identified a novel locus on 10q24.33 (rs4917385: Pgcadj  = 1.39 × 10(-8)) with an expression quantitative trait locus (eQTL) effect (Peqtl  = 8.0 × 10(-37) at USMG5/miR1307), and several new suggestive loci. SLE risk loci previously identified in Europeans and Asians were corroborated. Local ancestry estimation showed that the HLA allele risk contribution is of European ancestral origin. Imputation of HLA alleles suggested that autochthonous Native American haplotypes provide protection against development of SLE.

CONCLUSION

Our results demonstrate that studying admixed populations provides new insights in the delineation of the genetic architecture that underlies autoimmune and complex diseases.

The Role of the Transcription Factor Ets1 in Lupus and Other Autoimmune Diseases

Systemic lupus erythematosus (SLE) is a chronic autoimmune disease characterized by excess B- and T-cell activation, the development of autoantibodies against self-antigens including nuclear antigens, and immune complex deposition in target organs, which triggers an inflammatory response and tissue damage. The genetic and environmental factors that contribute to the development of SLE have been studied extensively in both humans and mouse models of the disease. One of the important genetic contributions to SLE development is an alteration in the expression of the transcription factor Ets1, which regulates the functional differentiation of lymphocytes. Here, we review the genetic, biochemical, and immunological studies that have linked low levels of Ets1 to aberrant lymphocyte differentiation and to the pathogenesis of SLE.