PLD4 is a genetic determinant to systemic lupus erythematosus and involved in murine autoimmune phenotypes

Structural and mechanistic insights into disease-associated endolysosomal exonucleases PLD3 and PLD4

Endolysosomal exonucleases PLD3 and PLD4 (phospholipases D3 and D4) are associated with autoinflammatory and autoimmune diseases. We report structures of these enzymes, and the molecular basis of their catalysis. The structures reveal an intra-chain dimer topology forming a basic active site at the interface. Like other PLD superfamily members, PLD3 and PLD4 carry HxKxxxxD/E motifs and participate in phosphodiester-bond cleavage. The enzymes digest ssDNA and ssRNA in a 5'-to-3' manner and are blocked by 5'-phosphorylation. We captured structures in apo, intermediate, and product states and revealed a "link-and-release" two-step catalysis. We also unexpectedly demonstrated phosphatase activity via a covalent 3-phosphohistidine intermediate. PLD4 contains an extra hydrophobic clamp that stabilizes substrate and could affect oligonucleotide substrate preference and product release. Biochemical and structural analysis of disease-associated mutants of PLD3/4 demonstrated reduced enzyme activity or thermostability and the possible basis for disease association. Furthermore, these findings provide insight into therapeutic design.

Human and Murine Toll-like Receptor-Driven Disease in Systemic Lupus Erythematosus

The pathogenesis of systemic lupus erythematosus (SLE) is linked to the differential roles of toll-like receptors (TLRs), particularly TLR7, TLR8, and TLR9. TLR7 overexpression or gene duplication, as seen with the Y-linked autoimmune accelerator (Yaa) locus or TLR7 agonist imiquimod, correlates with increased SLE severity, and specific TLR7 polymorphisms and gain-of-function variants are associated with enhanced SLE susceptibility and severity. In addition, the X-chromosome location of TLR7 and its escape from X-chromosome inactivation provide a genetic basis for female predominance in SLE. The absence of TLR8 and TLR9 have been shown to exacerbate the detrimental effects of TLR7, leading to upregulated TLR7 activity and increased disease severity in mouse models of SLE. The regulatory functions of TLR8 and TLR9 have been proposed to involve competition for the endosomal trafficking chaperone UNC93B1. However, recent evidence implies more direct, regulatory functions of TLR9 on TLR7 activity. The association between age-associated B cells (ABCs) and autoantibody production positions these cells as potential targets for treatment in SLE, but the lack of specific markers necessitates further research for precise therapeutic intervention. Therapeutically, targeting TLRs is a promising strategy for SLE treatment, with drugs like hydroxychloroquine already in clinical use.

GWAS for systemic sclerosis identifies six novel susceptibility loci including one in the Fcγ receptor region

Here we report the largest Asian genome-wide association study (GWAS) for systemic sclerosis performed to date, based on data from Japanese subjects and comprising of 1428 cases and 112,599 controls. The lead SNP is in the FCGR/FCRL region, which shows a penetrating association in the Asian population, while a complete linkage disequilibrium SNP, rs10917688, is found in a cis-regulatory element for IRF8. IRF8 is also a significant locus in European GWAS for systemic sclerosis, but rs10917688 only shows an association in the presence of the risk allele of IRF8 in the Japanese population. Further analysis shows that rs10917688 is marked with H3K4me1 in primary B cells. A meta-analysis with a European GWAS detects 30 additional significant loci. Polygenic risk scores constructed with the effect sizes of the meta-analysis suggest the potential portability of genetic associations beyond populations. Prioritizing the top 5% of SNPs of IRF8 binding sites in B cells improves the fitting of the polygenic risk scores, underscoring the roles of B cells and IRF8 in the development of systemic sclerosis. The results also suggest that systemic sclerosis shares a common genetic architecture across populations.

Structural analysis of PLD3 reveals insights into the mechanism of lysosomal 5' exonuclease-mediated nucleic acid degradation

The phospholipase D (PLD) family is comprised of enzymes bearing phospholipase activity towards lipids or endo- and exonuclease activity towards nucleic acids. PLD3 is synthesized as a type II transmembrane protein and proteolytically cleaved in lysosomes, yielding a soluble active form. The deficiency of PLD3 leads to the slowed degradation of nucleic acids in lysosomes and chronic activation of nucleic acid-specific intracellular toll-like receptors. While the mechanism of PLD phospholipase activity has been extensively characterized, not much is known about how PLDs bind and hydrolyze nucleic acids. Here, we determined the high-resolution crystal structure of the luminal N-glycosylated domain of human PLD3 in its apo- and single-stranded DNA-bound forms. PLD3 has a typical phospholipase fold and forms homodimers with two independent catalytic centers via a newly identified dimerization interface. The structure of PLD3 in complex with an ssDNA-derived thymidine product in the catalytic center provides insights into the substrate binding mode of nucleic acids in the PLD family. Our structural data suggest a mechanism for substrate binding and nuclease activity in the PLD family and provide the structural basis to design immunomodulatory drugs targeting PLD3.

Discovery of PLD4 modulators by high-throughput screening and kinetic analysis

Phospholipase D3 (PLD3) and D4 (PLD4) are endolysosomal exonucleases of ssDNA and ssRNA that regulate innate immunity. Polymorphisms of these enzymes are correlated with numerous human diseases, including Alzheimer's, rheumatoid arthritis, and systemic sclerosis. Pharmacological modulation of these immunoregulatory proteins may yield novel immunotherapies and adjuvants. A previous study reported a high-throughput screen (N = 17,952) that discovered a PLD3-selective activator and inhibitor, as well as a nonselective inhibitor, but failed to identify selective modulators of PLD4. However, modulators selective for PLD4 are therapeutically pertinent, since recent reports have shown that regulating this protein has direct implications in cancer and autoimmune diseases. Furthermore, the high expression of PLD4 in dendritic and myeloid cells, in comparison to the broader expression of PLD3, presents the opportunity for a cell-targeted immunotherapy. Here, we describe screening of an expended diversity library (N = 45,760) with an improved platform and report the discovery of one inhibitor and three activators selective for PLD4. Furthermore, kinetic modeling and structural analysis suggest mechanistic differences in the modulation of these hits. These findings further establish the utility of this screening platform and provide a set of chemical scaffolds to guide future small-molecule development for this novel immunoregulator target.

Structural and mechanistic insights into disease-associated endolysosomal exonucleases PLD3 and PLD4

Endolysosomal exonucleases PLD3 and PLD4 (phospholipases D3 and D4) are associated with autoinflammatory and autoimmune diseases. We report structures of these enzymes, and the molecular basis of their catalysis. The structures reveal an intra-chain dimer topology forming a basic active site at the interface. Like other PLD superfamily members, PLD3 and PLD4 carry HxKxxxxD/E motifs and participate in phosphodiester-bond cleavage. The enzymes digest ssDNA and ssRNA in a 5'-to-3' manner and are blocked by 5'-phosphorylation. We captured structures in apo, intermediate, and product states and revealed a 'link-and-release' two-step catalysis. We also unexpectedly demonstrated phosphatase activity via a covalent 3' phosphistidine intermediate. PLD4 contains an extra hydrophobic clamp that stabilizes substrate and could affect oligonucleotide substrate preference and product release. Biochemical and structural analysis of disease-associated mutants of PLD3/4 demonstrated reduced enzyme activity or thermostability and the possible basis for disease association. Furthermore, these findings provide insight into therapeutic design.

Phospholipase D4 as a signature of toll-like receptor 7 or 9 signaling is expressed on blastic T-bet + B cells in systemic lupus erythematosus

BACKGROUND

In systemic lupus erythematosus (SLE), autoreactive B cells are thought to develop by-passing immune checkpoints and contribute to its pathogenesis. Toll-like receptor (TLR) 7 and 9 signaling have been implicated in their development and differentiation. Although some B cell subpopulations such as T-bet + double negative 2 (DN2) cells have been identified as autoreactive in the past few years, because the upregulated surface markers of those cells are not exclusive to them, it is still challenging to specifically target autoreactive B cells in SLE patients.

METHODS

Our preliminary expression analysis revealed that phospholipase D4 (PLD4) is exclusively expressed in plasmacytoid dendritic cells (pDCs) and B cells in peripheral blood mononuclear cells (PBMCs) samples. Monoclonal antibodies against human PLD4 were generated, and flow cytometry analyses were conducted for PBMCs from 23 healthy donors (HDs) and 40 patients with SLE. In vitro cell culture was also performed to study the conditions that induce PLD4 in B cells from HDs. Finally, recombinant antibodies were synthesized from subpopulations of PLD4 + B cells from a patient with SLE, and their antinuclear activity was measured through enzyme-linked immunosorbent assay.

RESULTS

pDCs from both groups showed comparable frequency of surface PLD4 expression. PLD4 + B cells accounted for only a few percent of HD B cells, whereas they were significantly expanded in patients with SLE (2.1% ± 0.4% vs. 10.8% ± 1.2%, P < 0.005). A subpopulation within PLD4 + B cells whose cell size was comparable to CD38 + CD43 + plasmablasts was defined as "PLD4 + blasts," and their frequencies were significantly correlated with those of plasmablasts (P < 0.005). PLD4 + blasts phenotypically overlapped with double negative 2 (DN2) cells, and, in line with this, their frequencies were significantly correlated with several clinical markers of SLE. In vitro assay using healthy PBMCs demonstrated that TLR7 or TLR9 stimulation was sufficient to induce PLD4 on the surface of the B cells. Finally, two out of three recombinant antibodies synthesized from PLD4 + blasts showed antinuclear activity.

CONCLUSION

PLD4 + B cells, especially "blastic" ones, are likely autoreactive B cells undergoing TLR stimulation. Therefore, PLD4 is a promising target marker in SLE treatment.

Disease in the Pld4thss/thss Model of Murine Lupus Requires TLR9

Phospholipase D4 (PLD4) is an endolysosomal exonuclease of ssRNA and ssDNA, rather than a phospholipase as its name suggests. Human polymorphisms in the PLD4 gene have been linked by genome-wide association studies to systemic sclerosis, rheumatoid arthritis, and systemic lupus erythematosus. However, B6.129 Pld4-/- mice develop features of a distinct disease, macrophage activation syndrome, which is reversed in mice mutated in TLR9. In this article, we compare a Pld4 null mutant identified on the BALB/c background, Pld4thss/thss, which has distinct phenotypes: short stature, thin hair, and features of systemic lupus erythematosus. All phenotypes analyzed were largely normalized in Pld4thss/thssTlr9-/- mice. Thus, Pld4thss/thss represents a rare model in which mouse lupus etiology is TLR9 dependent. Compared with PLD4-deficient B6 mice, Pld4thss/thss mice had elevated levels of serum IgG, IgG anti-dsDNA autoantibodies, BAFF, and IFN-γ and elevated B cell numbers. Overall, the data suggest that PLD4 deficiency can lead to a diverse array of rheumatological abnormalities depending upon background-modifying genes, and that these diseases of PLD4 deficiency are largely driven by TLR9 recognition of ssDNA.

A review of genetic risk in systemic lupus erythematosus

INTRODUCTION

Systemic Lupus Erythematosus (SLE) is a complex multisystem autoimmune disease with a wide range of signs and symptoms in affected individuals. The utilization of genome-wide association study (GWAS) technology has led to an explosion in the number of genetic risk factors mapped for autoimmune diseases, including SLE.

AREAS COVERED

In this review, we summarize the more recent genetic risk loci mapped in SLE, which bring the total number of loci mapped to approximately 200. We review prioritization analyses of the associated variants and experimental validation of the putative causal variants. This includes the implementation of new bioinformatic techniques to align genomic and functional data and the use of transcriptomics with single-cell RNA-sequencing, CRISPR genome editing, and Massive Parallel Reporter Assays to analyze non-coding regulatory genetics.

EXPERT OPINION

Despite progress in identifying more genetic risk loci and variant-gene pairs for SLE, understanding its pathogenesis and applying findings clinically remains challenging. The polygenic risk score (PRS) has been used as an application of SLE genetics, but with limited performance in non-EUR populations. In the next few years, advancements in proteomics, post-translational modification estimation, and whole-genome sequencing will enhance disease understanding.

Biological insights into systemic lupus erythematosus through an immune cell-specific transcriptome-wide association study

OBJECTIVE

Genome-wide association studies (GWAS) have identified >100 risk loci for systemic lupus erythematosus (SLE), but the disease genes at most loci remain unclear, hampering translation of these genetic discoveries. We aimed to prioritise genes underlying the 110 SLE loci that were identified in the latest East Asian GWAS meta-analysis.

METHODS

We built gene expression predictive models in blood B cells, CD4+ and CD8+ T cells, monocytes, natural killer cells and peripheral blood cells of 105 Japanese individuals. We performed a transcriptome-wide association study (TWAS) using data from the latest genome-wide association meta-analysis of 208 370 East Asians and searched for candidate genes using TWAS and three data-driven computational approaches.

RESULTS

TWAS identified 171 genes for SLE (p<1.0×10-5); 114 (66.7%) showed significance only in a single cell type; 127 (74.3%) were in SLE GWAS loci. TWAS identified a strong association between CD83 and SLE (p<7.7×10-8). Meta-analysis of genetic associations in the existing 208 370 East Asian and additional 1498 cases and 3330 controls found a novel single-variant association at rs72836542 (OR=1.11, p=4.5×10-9) around CD83. For the 110 SLE loci, we identified 276 gene candidates, including 104 genes at recently-identified SLE novel loci. We demonstrated in vitro that putative causal variant rs61759532 exhibited an allele-specific regulatory effect on ACAP1, and that presence of the SLE risk allele decreased ACAP1 expression.

CONCLUSIONS

Cell-level TWAS in six types of immune cells complemented SLE gene discovery and guided the identification of novel genetic associations. The gene findings shed biological insights into SLE genetic associations.

Exonucleases: Degrading DNA to Deal with Genome Damage, Cell Death, Inflammation and Cancer

Although DNA degradation might seem an unwanted event, it is essential in many cellular processes that are key to maintaining genomic stability and cell and organism homeostasis. The capacity to cut out nucleotides one at a time from the end of a DNA chain is present in enzymes called exonucleases. Exonuclease activity might come from enzymes with multiple other functions or specialized enzymes only dedicated to this function. Exonucleases are involved in central pathways of cell biology such as DNA replication, repair, and death, as well as tuning the immune response. Of note, malfunctioning of these enzymes is associated with immune disorders and cancer. In this review, we will dissect the impact of DNA degradation on the DNA damage response and its links with inflammation and cancer.

Nucleic Acid Sensing by Toll-Like Receptors in the Endosomal Compartment

Toll-like receptors (TLRs) respond to pathogen constituents, such as microbial lipids and nucleic acids (NAs). TLRs recognize NAs in endosomal compartments. Structural and functional studies have shown that recognition of NAs by TLRs depends on NA processing by RNases and DNases. DNase II-dependent DNA degradation is required for TLR9 responses to single-stranded DNAs, whereas RNase T2-dependent RNA degradation enables TLR7 and TLR8 to respond to nucleosides and oligoribonucleotides. In contrast, RNases and DNases negatively regulate TLR responses by degrading their ligands. RNase T2 negatively regulates TLR3 responses to degrading the TLR3 ligand double-stranded RNAs. Therefore, NA metabolism in the endosomal compartments affects the endosomal TLR responses. Dysregulation of NA metabolism in the endosomal compartment drives the TLR-dependent pathologies in human diseases.

Novel susceptibility loci for steroid-associated osteonecrosis of the femoral head in systemic lupus erythematosus

INTRODUCTION

Osteonecrosis of the femoral head (ONFH) involves necrosis of bone and bone marrow of the femoral head caused by ischemia with unknown etiology. Previous genetic studies on ONFH failed to produce consistent results, presumably because ONFH has various causes with different genetic backgrounds and the underlying diseases confounded the associations. Steroid-associated ONFH (S-ONFH) accounts for one-half of all ONFH, and systemic lupus erythematosus (SLE) is a representative disease underlying S-ONFH. We performed a genome-wide association study (GWAS) to identify genetic risk factors for S-ONFH in patients with SLE.

METHODS

We conducted a two-staged GWAS on 636 SLE patients with S-ONFH and 95 588 non-SLE controls. Among the novel loci identified, we determined S-ONFH specific loci by comparing allele frequencies between SLE patients without S-ONFH and non-SLE controls. We also used Korean datasets comprising 148 S-ONFH cases and 37 015 controls to assess overall significance. We evaluated the functional annotations of significant variants by in-silico analyses.

RESULTS

The Japanese GWAS identified four significant loci together with 12 known SLE susceptibility loci. The four significant variants showed comparable effect sizes on S-ONFH compared with SLE controls and non-SLE controls. Three of the four loci, MIR4293/MIR1265 (OR = 1.99, P-value = 1.1 × 10-9), TRIM49/NAALAD2 (OR = 1.65, P-value = 4.8 × 10-8) and MYO16 (OR = 3.91, P-value = 4.9 × 10-10), showed significant associations in the meta-analysis with Korean datasets. Bioinformatics analyses identified MIR4293, NAALAD2 and MYO16 as candidate causal genes. MIR4293 regulates a PPARG-related adipogenesis pathway relevant to S-ONFH.

CONCLUSIONS

We identified three novel susceptibility loci for S-ONFH in SLE.

Therapeutic effects of cationic liposomes on lupus-prone MRL/lpr mice are mediated via inhibition of TLR4-triggered B-cell activation

We previously reported that co-delivery of dihydroartemisinin and high mobility group box 1 (HMGB1) siRNAs, using cell penetrating peptide (TAT)-modified cationic liposomes (TAT-CLs-DHA/siRNA), resulted in promising activity for the treatment of inflammatory disease through TLR4 signaling pathway. In the current study, we further investigated the therapeutic effects of TAT-CLs-DHA/siRNA on lupus-prone MRL/lpr mice and explored its effects on B cell responses. In vitro, we found that TAT-CLs-DHA/siRNA suppressed the proliferation and activation of B cells through the TLR4 signaling pathway. Following parenteral administration every 4 days, TAT-CLs-DHA/siRNA significantly reduced proteinuria, glomerulonephritis, serum anti-dsDNA antibody and secretion of interleukin (IL)-6, IL-10, IL-17 and IL-21. Moreover, Western blotting showed that TAT-CLs-DHA/siRNA modulated the B-cell intrinsic pathway by downregulating expression of HMGB1, TLR4, MyD88 and NF-κB. This co-delivery system thus represents a promising treatment option for lupus nephritis, and also highlights a novel target of lupus treatment through B cell TLR4 signal pathway.

N 6-Methyladenosine Methylomic Landscape of Lung Tissues in Murine Acute Allergic Asthma

Allergic asthma is well known as a common respiratory disorder comprising an allergic inflammatory nature and excessive immune characteristic. N ⁶-methyladenosine (m6A) methylation is an RNA epigenetic modification that post-transcriptionally regulates gene expression and function by affecting the RNA fate. Currently, m6A methylation is gaining attention as a mechanism of immunoregulation. However, whether m6A methylation engages the pathological process of asthma remains uncertain. Here, we present the m6A methylomic landscape in the lung tissues of ovalbumin-induced acute asthma mice using MeRIP-seq and RNA-seq. We identified 353 hypermethylated m6A peaks within 329 messenger RNAs (mRNAs) and 150 hypomethylated m6A peaks within 143 mRNAs in the lung tissues of asthmatic mice. These differentially methylated mRNAs were found to be involved in several immune function-relevant signaling pathways. In addition, we predicted 25 RNA-binding proteins that recognize the differentially methylated peak sites by exploring public databases, and the roles of these proteins are mostly related to mRNA biogenesis and metabolism. To further investigate the expression levels of the differentially methylated genes, we performed combined analysis of the m6A methylome and transcriptome data and identified 127 hypermethylated mRNAs (107 high and 20 low expression) and 43 hypomethylated mRNAs with differential expressions (9 high and 34 low expression). Of these, there are a list of mRNAs involved in immune function and regulation. The present results highlight the essential role of m6A methylation in the pathogenesis of asthma.

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.

Identification of a risk prediction model for clinical prognosis in HER2 positive breast cancer patients

Background New biomarkers are needed to identify different clinical outcomes for HER2+ breast cancer (BC). Methods Differential genes of HER2+ BC were screened based on TCGA database. We used WGCNA to identify the genes related to the survival. Genetic Algorithm was used to structure risk prediction model. The prognostic model was validated in GSE data. Results We constructed a risk prediction model of 6 genes to identify prognosis of HER2+ BC, including CLEC9A, PLD4, PIM1, PTK2B, AKNAD1 and C15orf27. Kaplan-Meier curve showed that the model effectively distinguished the survival of HER2+ BC patients. The multivariate Cox regression suggested that the risk model was an independent predictor for HER2+ BC. Analysis related to immune showed that significant differences in immune infiltration between high- and low-risk groups classified by the prognostic model. Conclusions Our study identified a risk prediction model of 6 genes that could distinguish the prognosis of HER2+ BC.

Cleavage of DNA and RNA by PLD3 and PLD4 limits autoinflammatory triggering by multiple sensors

Phospholipase D3 (PLD3) and PLD4 polymorphisms have been associated with several important inflammatory diseases. Here, we show that PLD3 and PLD4 digest ssRNA in addition to ssDNA as reported previously. Moreover, Pld3^−/−Pld4^−/− mice accumulate small ssRNAs and develop spontaneous fatal hemophagocytic lymphohistiocytosis (HLH) characterized by inflammatory liver damage and overproduction of Interferon (IFN)-γ. Pathology is rescued in Unc93b1^3d/3dPld3^−/−Pld4^−/− mice, which lack all endosomal TLR signaling; genetic codeficiency or antibody blockade of TLR9 or TLR7 ameliorates disease less effectively, suggesting that both RNA and DNA sensing by TLRs contributes to inflammation. IFN-γ made a minor contribution to pathology. Elevated type I IFN and some other remaining perturbations in Unc93b1^3d/3dPld3^−/−Pld4^−/− mice requires STING (Tmem173). Our results show that PLD3 and PLD4 regulate both endosomal TLR and cytoplasmic/STING nucleic acid sensing pathways and have implications for the treatment of nucleic acid-driven inflammatory disease.

Loss of function polymorphisms of phospholipase D3 and D4 are associated with inflammatory diseases and their function is unclear. Here the authors show that PLD3/4 function as RNAses and deletion of these proteins in mice leads to accumulation of ssRNA which exacerbates inflammation through TLR signalling.

Azithromycin modulates Teff/Treg balance in retinal inflammation via the mTOR signaling pathway

Uveitis is one of the most common blindness-causing ocular disorders. Due to its complicated pathogenesis, the treatment of uveitis has been widely recognized as a challenge for ophthalmologists. Recently, the anti-inflammatory properties of the antibiotic Azithromycin (AZM) have been reported. However, the therapeutic effects of Azithromycin in experimental autoimmune uveitis (EAU), a representative model of human AU, have not been elucidated till date. We conducted this study to examine the therapeutic effects and potential mechanisms of Azithromycin in EAU. We observed that Azithromycin significantly attenuated retinal inflammation in EAU mice at day 14 after immunization along with a significantly decreased inflammatory cell infiltration and cytokine production in the retina. Furthermore, we observed that Azithromycin increased the number of regulatory T cells (Treg) and decreased the number of effector T cells (Teff) in both the draining lymph nodes and spleen of EAU mice. Additionally, Azithromycin suppressed the proliferation and activation of CD4 + T cells, and induced the apoptosis of CD4 + CD44 + memory T and CD4 + CXCR3 + Th1 cells. Mechanistically, we proved that Azithromycin could regulate Teff/Treg balance by inhibiting the phosphorylation of S6 ribosomal protein, a downstream target of mammalian target of rapamycin (mTOR). Together, our findings revealed that Azithromycin alleviated EAU by regulating the Teff/Treg balance through the mTOR signaling pathway, suggesting that Azithromycin could be a promising therapeutic candidate for AU.

Autoimmunity in 2019

Based on the PubMed data, we have been performing a yearly evaluation of the publications related to autoimmune diseases and immunology to ascertain the relative weight of the former in the scientific literature. It is particularly intriguing to observe that despite the numerous new avenues of immune-related mechanisms, such as cancer immunotherapy, the proportion of immunology manuscripts related to autoimmunity continues to increase and has been approaching 20% in 2019. As in the previous 13 years, we performed an arbitrary selection of the peer-reviewed articles published by the major dedicated Journals and discussed the common themes which continue to outnumber peculiarites in autoimmune diseases. The investigated areas included systemic lupus erythematosus (SLE), rheumatoid arthritis (RA), psoriatic arthritis (PsA), autoantibodies (autoAbs), and common therapeutic avenues and novel pathogenic mechanisms for autoimmune conditions. Some examples include new pathogenetic evidence which is well represented by IL21 or P2X7 receptor (P2X7R) in SLE or the application of single-cell RNA sequencing (scRNA-seq), mass cytometry, bulk RNA sequencing (RNA-seq), and flow cytometry for the analysis of different cellular populations in RA. Cumulatively and of interest to the clinicians, a large number of findings continue to underline the importance of a strict relationship between basic and clinical science to define new pathogenetic and therapeutic developments. The therapeutic pipeline in autoimmunity continues to grow and maintain a constant flow of new molecules, as well illustrated in RA and PsA, and this is most certainly derived from the new basic evidence and the high-throughput tools applied to autoimmune diseases.

Modulators of immunoregulatory exonucleases PLD3 and PLD4 identified by high-throughput screen

PLD3 and PLD4 have recently been revealed to be endosomal exonucleases that regulate the innate immune response by digesting the ligands of nucleic acid sensors. These enzymes can suppress RNA and DNA innate immune sensors like toll-like receptor 9, and PLD4-deficent mice exhibit inflammatory disease. Targeting these immunoregulatory enzymes presents an opportunity to indirectly regulate innate immune nucleic acid sensors that could yield immunotherapies, adjuvants, and nucleic acid drug stabilizers. To aid in delineating the therapeutic potential of these targets, we have developed a high-throughput fluorescence enzymatic assay to identify modulators of PLD3 and PLD4. Screening of a diversity library (N = 17952) yielded preferential inhibitors of PLD3 and PLD4 in addition to a PLD3 selective activator. The modulation models of these compounds were delineated by kinetic analysis. This work presents an inexpensive and simple method to identify modulators of these immunoregulatory exonucleases.

Lupus susceptibility region containing CDKN1B rs34330 mechanistically influences expression and function of multiple target genes, also linked to proliferation and apoptosis

OBJECTIVE

A recent genome-wide association study (GWAS) reported a significant genetic association between rs34330 of cyclin-dependent kinase inhibitor 1B (CDKN1B) and risk of systemic lupus erythematosus (SLE) in Han Chinese. This study aims to validate the reported association and elucidate the biochemical mechanisms underlying the variant's effect.

METHODS

We performed allelic association with SLE followed by meta-analysis across 11 independent cohorts (n=28,872). We applied in silico bioinformatics and experimental validation in SLE-relevant cell lines to determine the functional consequences of rs34330.

RESULTS

We replicated genetic association between SLE and rs34330 (P_meta =5.29x10^-22 , OR (95% CI)=0.84 (0.81-0.87)). Follow-up bioinformatics and eQTL analysis suggest that rs34330 is located in active chromatin and potentially regulates several target genes. Using luciferase and ChIP-qPCR, we demonstrated substantial allele-specific promoter and enhancer activity, and allele-specific binding of three histone marks (H3K27ac, H3K4me3, H3K4me1), RNA pol II, CTCF, and a critical immune transcription factor (IRF-1). Chromosome conformation capture (3C) detected long-range chromatin interactions between rs34330 and the promoters of neighboring genes APOLD1 and DDX47, and effects on CDKN1B and the other target genes were directly validated by CRISPR-based genome editing. Finally, CRISPR-dCas9-based epigenetic activation/silencing confirmed these results. Gene-edited cell lines also showed higher levels of proliferation and apoptosis.

CONCLUSION

Collectively, these findings suggest a mechanism whereby the rs34330 risk allele (C) influences the presence of histone marks, RNA pol II, and the IRF-1 transcription factor to regulate expression of several target genes linked to proliferation and apoptosis, which potentially underlie the association of rs34330 with SLE.

Meta-analysis of 208370 East Asians identifies 113 susceptibility loci for systemic lupus erythematosus

OBJECTIVE

Systemic lupus erythematosus (SLE), an autoimmune disorder, has been associated with nearly 100 susceptibility loci. Nevertheless, these loci only partially explain SLE heritability and their putative causal variants are rarely prioritised, which make challenging to elucidate disease biology. To detect new SLE loci and causal variants, we performed the largest genome-wide meta-analysis for SLE in East Asian populations.

METHODS

We newly genotyped 10 029 SLE cases and 180 167 controls and subsequently meta-analysed them jointly with 3348 SLE cases and 14 826 controls from published studies in East Asians. We further applied a Bayesian statistical approach to localise the putative causal variants for SLE associations.

RESULTS

We identified 113 genetic regions including 46 novel loci at genome-wide significance (p<5×10-8). Conditional analysis detected 233 association signals within these loci, which suggest widespread allelic heterogeneity. We detected genome-wide associations at six new missense variants. Bayesian statistical fine-mapping analysis prioritised the putative causal variants to a small set of variants (95% credible set size ≤10) for 28 association signals. We identified 110 putative causal variants with posterior probabilities ≥0.1 for 57 SLE loci, among which we prioritised 10 most likely putative causal variants (posterior probability ≥0.8). Linkage disequilibrium score regression detected genetic correlations for SLE with albumin/globulin ratio (rg=-0.242) and non-albumin protein (rg=0.238).

CONCLUSION

This study reiterates the power of large-scale genome-wide meta-analysis for novel genetic discovery. These findings shed light on genetic and biological understandings of SLE.

Dynamic landscape of immune cell-specific gene regulation in immune-mediated diseases

Genetic studies have revealed many variant loci that are associated with immune-mediated diseases. To elucidate the disease pathogenesis, it is essential to understand the function of these variants, especially under disease-associated conditions. Here, we performed a large-scale immune cell gene-expression analysis, together with whole-genome sequence analysis. Our dataset consists of 28 distinct immune cell subsets from 337 patients diagnosed with 10 categories of immune-mediated diseases and 79 healthy volunteers. Our dataset captured distinctive gene-expression profiles across immune cell types and diseases. Expression quantitative trait loci (eQTL) analysis revealed dynamic variations of eQTL effects in the context of immunological conditions, as well as cell types. These cell-type-specific and context-dependent eQTLs showed significant enrichment in immune disease-associated genetic variants, and they implicated the disease-relevant cell types, genes, and environment. This atlas deepens our understanding of the immunogenetic functions of disease-associated variants under in vivo disease conditions.

PLD3 is a neuronal lysosomal phospholipase D associated with β-amyloid plaques and cognitive function in Alzheimer's disease

Phospholipase D3 (PLD3) is a protein of unclear function that structurally resembles other members of the phospholipase D superfamily. A coding variant in this gene confers increased risk for the development of Alzheimer's disease (AD), although the magnitude of this effect has been controversial. Because of the potential significance of this obscure protein, we undertook a study to observe its distribution in normal human brain and AD-affected brain, determine whether PLD3 is relevant to memory and cognition in sporadic AD, and to evaluate its molecular function. In human neuropathological samples, PLD3 was primarily found within neurons and colocalized with lysosome markers (LAMP2, progranulin, and cathepsins D and B). This colocalization was also present in AD brain with prominent enrichment on lysosomal accumulations within dystrophic neurites surrounding β-amyloid plaques. This pattern of protein distribution was conserved in mouse brain in wild type and the 5xFAD mouse model of cerebral β-amyloidosis. We discovered PLD3 has phospholipase D activity in lysosomes. A coding variant in PLD3 reported to confer AD risk significantly reduced enzymatic activity compared to wild-type PLD3. PLD3 mRNA levels in the human pre-frontal cortex inversely correlated with β-amyloid pathology severity and rate of cognitive decline in 531 participants enrolled in the Religious Orders Study and Rush Memory and Aging Project. PLD3 levels across genetically diverse BXD mouse strains and strains crossed with 5xFAD mice correlated strongly with learning and memory performance in a fear conditioning task. In summary, this study identified a new functional mammalian phospholipase D isoform which is lysosomal and closely associated with both β-amyloid pathology and cognition.

Quantification and characterization of the 5' exonuclease activity of the lysosomal nuclease PLD3 by a novel cell-based assay

Phospholipase D3 (PLD3) and phospholipase D4 (PLD4), the most recently described lysosomal nucleases, are associated with Alzheimer’s disease, spinocerebellar ataxia, and systemic lupus erythematosus. They exhibit 5′ exonuclease activity on single-stranded DNA, hydrolyzing it at the acidic pH associated with the lysosome. However, their full cellular function is inadequately understood. To examine these enzymes, we developed a robust and automatable cell-based assay based on fluorophore- and fluorescence-quencher-coupled oligonucleotides for the quantitative determination of acidic 5′ exonuclease activity. We validated the assay under knockout and PLD-overexpression conditions and then applied it to characterize PLD3 and PLD4 biochemically. Our experiments revealed PLD3 as the principal acid 5′ exonuclease in HeLa cells, where it showed a markedly higher specific activity compared with PLD4. We further used our newly developed assay to determine the substrate specificity and inhibitory profile of PLD3 and found that proteolytic processing of PLD3 is dispensable for its hydrolytic activity. We followed the expression, proteolytic processing, and intracellular distribution of genetic PLD3 variants previously associated with Alzheimer’s disease and investigated each variant's effect on the 5′ nuclease activity of PLD3, finding that some variants lead to reduced activity, but others not. The development of a PLD3/4-specific biochemical assay will be instrumental in understanding better both nucleases and their incompletely understood roles in vitro and in vivo.

Bioinformatics analysis of the clinical value and potential mechanisms of AHNAK2 in papillary thyroid carcinoma

BACKGROUND

AHNAK2 has been recently reported as a biomarker in many cancers. However, a systematic investigation of AHNAK2 in papillary thyroid carcinoma (PTC) has not been conducted.

RESULTS

AHNAK2 is overexpressed in PTC tissues and could be an independent prognostic factor. AHNAK2 expression was significantly high in patients with advanced stage, advanced T classification, lymph node metastasis, increased BRAF mutations and decreased RAS mutations. Cell adhesion-, cell junction-, and immune-related pathways were the most frequently noted in gene set enrichment analysis. AHNAK2 expression in PTC was positively correlated with immune infiltration and negatively correlated with AHNAK2 methylation. AHNAK2 expression was significantly positively correlated with tumor progression and poor overall survival (OS) in pan-cancer patients.

CONCLUSIONS

AHNAK2 is a good biomarker for the diagnosis and prognosis of PTC. AHNAK2 may promote thyroid cancer progression through cell adhesion-, cell junction-, and immune-related pathways. Methylation may act as an upstream regulator to inhibit the expression and biological function of AHNAK2. Additionally, AHNAK2 has broad prognostic value in pan-cancer.

METHODS

Based on The Cancer Genome Atlas (TCGA) data, we screened AHNAK2-related genes through weighted gene coexpression network analysis and explored the clinical value and the potential mechanism of AHNAK2 in PTC by multiomics analysis.

Meta-analysis of genome-wide association study identifies FBN2 as a novel locus associated with systemic lupus erythematosus in Thai population

BACKGROUND

Differences in the expression of variants across ethnic groups in the systemic lupus erythematosus (SLE) patients have been well documented. However, the genetic architecture in the Thai population has not been thoroughly examined. In this study, we carried out genome-wide association study (GWAS) in the Thai population.

METHODS

Two GWAS cohorts were independently collected and genotyped: discovery dataset (487 SLE cases and 1606 healthy controls) and replication dataset (405 SLE cases and 1590 unrelated disease controls). Data were imputed to the density of the 1000 Genomes Project Phase 3. Association studies were performed based on different genetic models, and pathway enrichment analysis was further examined. In addition, the performance of disease risk estimation for individuals in Thai GWAS was assessed based on the polygenic risk score (PRS) model trained by other Asian populations.

RESULTS

Previous findings on SLE susceptible alleles were well replicated in the two GWAS. The SNPs on HLA class II (rs9270970, A>G, OR = 1.82, p value = 3.61E-26), STAT4 (rs7582694, C>G, OR = 1.57, p value = 8.21E-16), GTF2I (rs73366469, A>G, OR = 1.73, p value = 2.42E-11), and FAM167A-BLK allele (rs13277113, A>G, OR = 0.68, p value = 1.58E-09) were significantly associated with SLE in Thai population. Meta-analysis of the two GWAS identified a novel locus at the FBN2 that was specifically associated with SLE in the Thai population (rs74989671, A>G, OR = 1.54, p value = 1.61E-08). Functional analysis showed that rs74989671 resided in a peak of H3K36me3 derived from CD14+ monocytes and H3K4me1 from T lymphocytes. In addition, we showed that the PRS model trained from the Chinese population could be applied in individuals of Thai ancestry, with the area under the receiver-operator curve (AUC) achieving 0.76 for this predictor.

CONCLUSIONS

We demonstrated the genetic architecture of SLE in the Thai population and identified a novel locus associated with SLE. Also, our study suggested a potential use of the PRS model from the Chinese population to estimate the disease risk for individuals of Thai ancestry.