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

A Non-Coding Variant in SLC15A4 Modulates Enhancer Activity and Lysosomal Deacidification Linked to Lupus Susceptibility

Systemic lupus erythematosus (SLE) is a complex autoimmune disease with a strong genetic basis. Despite the identification of several single nucleotide polymorphisms (SNPs) near the SLC15A4 gene that are significantly associated with SLE across multiple populations, specific causal SNP(s) and molecular mechanisms responsible for disease susceptibility are unknown. To address this gap, we employed bioinformatics, expression quantitative trait loci (eQTLs), and 3D chromatin interaction analysis to nominate a likely functional variant, rs35907548, in an active intronic enhancer of SLC15A4. Through luciferase reporter assays followed by chromatin immunoprecipitation (ChIP)-qPCR, we observed significant allele-specific enhancer effects of rs35907548 in diverse cell lines. The rs35907548 risk allele T is associated with increased regulatory activity and target gene expression, as shown by eQTLs and chromosome conformation capture (3C)-qPCR. The latter revealed long-range chromatin interactions between the rs35907548 enhancer and the promoters of SLC15A4, GLTLD1, and an uncharacterized lncRNA. The enhancer-promoter interactions and expression effects were validated by CRISPR/Cas9 knock-out (KO) of the locus in HL60 promyeloblast cells. KO cells also displayed dramatically dysregulated endolysosomal pH regulation. Together, our data show that the rs35907548 risk allele affects multiple aspects of cellular physiology and may directly contribute to SLE.

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.

Downregulation of circETS1 disrupts Th17/Treg homeostasis by inhibiting FOXP3 transcription: A new potential biomarker in systemic lupus erythematosus

Systemic lupus erythematosus (SLE) is an autoimmune disease associated with an imbalance of T helper 17 (Th17) to regulatory T cells (Tregs). However, the underlying mechanism remains unclear. Increasing evidence suggests that circular RNAs play a crucial role in SLE. Although circETS1 was discovered 30 years ago, detailed exploration of its functions remains limited. In this study, we measured the expression levels of circETS1 in peripheral blood mononuclear cells (PBMCs) and CD4+ T cells of patients with SLE by quantitative polymerase chain reaction. The impact of circETS1 expression on the Th17/Treg balance and underlying mechanism were evaluated using double-luciferase reporter, gain-/loss-of-function, and rescue assays. Receiver operating characteristic (ROC) curve analysis was conducted to assess the diagnostic value of circETS1. Both circETS1 and FOXP3 expression were downregulated in the PBMCs and CD4+ T cells of patients with SLE (n = 28) compared with those in the cells of healthy controls (n = 20). Mechanistically, we found that circETS1 can bind directly to the microRNA miR-1205, acting as a sponge to upregulate the transcription of FOXP3, thereby maintaining the Th17/Treg balance. Notably, ROC analysis showed that the expression level of circETS1 in PBMCs had an area under the curve of 0.873 (95% confidence interval: 0.771-0.976; p < .001) for diagnosing SLE, with a sensitivity of 80.00% and a specificity of 89.29%. Finally, we found negative correlations between the level of circETS1 in PBMCs and disease severity (according to the Systemic Lupus Erythematosus Disease Activity Index) in patients with SLE (r = -0.7712, 95% CI: -0.8910 to -0.5509; p < .001). The imbalance in Th17/Treg cells in SLE may be attributed, in part, to the circETS1/miR-1205/FOXP3 pathway. CircETS1 has potential to serve as a valuable biomarker for the diagnosis and evaluation of SLE.

A Non-Coding Variant in SLC15A4 Modulates Enhancer Activity and Lysosomal Deacidification Linked to Lupus Susceptibility

Systemic lupus erythematosus (SLE) is a complex autoimmune disease with a strong genetic basis. Despite the identification of several single nucleotide polymorphisms (SNPs) near the SLC15A4 gene that are significantly associated with SLE across multiple populations, specific causal SNP(s) and molecular mechanisms responsible for disease susceptibility are unknown. To address this gap, we employed bioinformatics, expression quantitative trait loci (eQTLs), and 3D chromatin interaction analysis to nominate a likely functional variant, rs35907548, in an active intronic enhancer of SLC15A4 . Through luciferase reporter assays followed by chromatin immunoprecipitation (ChIP)-qPCR, we observed significant allele-specific enhancer effects of rs35907548 in diverse cell lines. The rs35907548 risk allele T is associated with increased regulatory activity and target gene expression, as shown by eQTLs and chromosome conformation capture (3C)-qPCR. The latter revealed long-range chromatin interactions between the rs35907548 enhancer and the promoters of SLC15A4, GLTLD1 , and an uncharacterized lncRNA. The enhancer-promoter interactions and expression effects were validated by CRISPR/Cas9 knock-out (KO) of the locus in HL60 promyeloblast cells. KO cells also displayed dramatically dysregulated endolysosomal pH regulation. Together, our data show that the rs35907548 risk allele affects multiple aspects of cellular physiology and may directly contribute to SLE.

Ets1 and IL17RA cooperate to regulate autoimmune responses and skin immunity to Staphylococcus aureus

Introduction

Ets1 is a lymphoid-enriched transcription factor that regulates B- and Tcell functions in development and disease. Mice that lack Ets1 (Ets1 KO) develop spontaneous autoimmune disease with high levels of autoantibodies. Naïve CD4 + T cells isolated from Ets1 KO mice differentiate more readily to Th17 cells that secrete IL-17, a cytokine implicated in autoimmune disease pathogenesis. To determine if increased IL-17 production contributes to the development of autoimmunity in Ets1 KO mice, we crossed Ets1 KO mice to mice lacking the IL-17 receptor A subunit (IL17RA KO) to generate double knockout (DKO) mice.

Methods

In this study, the status of the immune system of DKO and control mice was assessed utilizing ELISA, ELISpot, immunofluorescent microscopy, and flow cytometric analysis of the spleen, lymph node, skin. The transcriptome of ventral neck skin was analyzed through RNA sequencing. S. aureus clearance kinetics in in exogenously infected mice was conducted using bioluminescent S. aureus and tracked using an IVIS imaging experimental scheme.

Results

We found that the absence of IL17RA signaling did not prevent or ameliorate the autoimmune phenotype of Ets1 KO mice but rather that DKO animals exhibited worse symptoms with striking increases in activated B cells and secreted autoantibodies. This was correlated with a prominent increase in the numbers of T follicular helper (Tfh) cells. In addition to the autoimmune phenotype, DKO mice also showed signs of immunodeficiency and developed spontaneous skin lesions colonized by Staphylococcus xylosus. When DKO mice were experimentally infected with Staphylococcus aureus, they were unable to clear the bacteria, suggesting a general immunodeficiency to staphylococcal species. γδ T cells are important for the control of skin staphylococcal infections. We found that mice lacking Ets1 have a complete deficiency of the γδ T-cell subset dendritic epidermal T cells (DETCs), which are involved in skin woundhealing responses, but normal numbers of other skin γδ T cells. To determine if loss of DETC combined with impaired IL-17 signaling might promote susceptibility to staph infection, we depleted DETC from IL17RA KO mice and found that the combined loss of DETC and impaired IL-17 signaling leads to an impaired clearance of the infection.

Conclusions

Our studies suggest that loss of IL-17 signaling can result in enhanced autoimmunity in Ets1 deficient autoimmune-prone mice. In addition, defects in wound healing, such as that caused by loss of DETC, can cooperate with impaired IL-17 responses to lead to increased susceptibility to skin staph infections.

Genetic dissection of TLR9 reveals complex regulatory and cryptic proinflammatory roles in mouse lupus

In lupus, Toll-like receptor 7 (TLR7) and TLR9 mediate loss of tolerance to RNA and DNA, respectively. Yet, TLR7 promotes disease, while TLR9 protects from disease, implying differences in signaling. To dissect this 'TLR paradox', we generated two TLR9 point mutants (lacking either ligand (TLR9^K51E) or MyD88 (TLR9^P915H) binding) in lupus-prone MRL/lpr mice. Ameliorated disease of Tlr9^K51E mice compared to Tlr9^-/- controls revealed a TLR9 'scaffold' protective function that is ligand and MyD88 independent. Unexpectedly, Tlr9^P915H mice were more protected than both Tlr9^K51E and Tlr9^WT mice, suggesting that TLR9 also possesses ligand-dependent, but MyD88-independent, regulatory signaling and MyD88-mediated proinflammatory signaling. Triple-mixed bone marrow chimeras showed that TLR9-MyD88-independent regulatory roles were B cell intrinsic and restrained differentiation into pathogenic age-associated B cells and plasmablasts. These studies reveal MyD88-independent regulatory roles of TLR9, shedding light on the biology of endosomal TLRs.

WhichTF is functionally important in your open chromatin data?

We present WhichTF, a computational method to identify functionally important transcription factors (TFs) from chromatin accessibility measurements. To rank TFs, WhichTF applies an ontology-guided functional approach to compute novel enrichment by integrating accessibility measurements, high-confidence pre-computed conservation-aware TF binding sites, and putative gene-regulatory models. Comparison with prior sheer abundance-based methods reveals the unique ability of WhichTF to identify context-specific TFs with functional relevance, including NF-κB family members in lymphocytes and GATA factors in cardiac cells. To distinguish the transcriptional regulatory landscape in closely related samples, we apply differential analysis and demonstrate its utility in lymphocyte, mesoderm developmental, and disease cells. We find suggestive, under-characterized TFs, such as RUNX3 in mesoderm development and GLI1 in systemic lupus erythematosus. We also find TFs known for stress response, suggesting routine experimental caveats that warrant careful consideration. WhichTF yields biological insight into known and novel molecular mechanisms of TF-mediated transcriptional regulation in diverse contexts, including human and mouse cell types, cell fate trajectories, and disease-associated cells.

Transcription factors (TFs), a class of DNA binding proteins, regulate tissue- and cell-type-specific expression of genes. Identifying the critical TFs in a given cellular context leads to investigating molecular regulatory mechanisms in development, differentiation, and disease. Because there are more than 1,500 human TFs, experimental measurements of genome-wide occupancy across all TFs have been challenging. While computational approaches play pivotal roles, most existing methods rely on statistical enrichment, focusing either on sequence motif similarity recognized by TFs or the similarity of the genomic region of interest with the previously characterized TF occupancy profile. Here we propose WhichTF as an alternative, incorporating curated biomedical knowledge from ontology and integrating it with the high-confidence prediction of conserved TF binding sites in user-provided genomic regions of interest. We develop a new WhichTF score to rank TFs and demonstrate its applicability across human and mouse cell types, cellular differentiation trajectories, and disease-associated cells.

NFAT-dependent and -independent exhaustion circuits program maternal CD8 T cell hypofunction in pregnancy

Pregnancy is a common immunization event, but the molecular mechanisms and immunological consequences provoked by pregnancy remain largely unknown. We used mouse models and human transplant registry data to reveal that pregnancy induced exhausted CD8 T cells (Preg-TEX), which associated with prolonged allograft survival. Maternal CD8 T cells shared features of exhaustion with CD8 T cells from cancer and chronic infection, including transcriptional down-regulation of ribosomal proteins and up-regulation of TOX and inhibitory receptors. Similar to other models of T cell exhaustion, NFAT-dependent elements of the exhaustion program were induced by fetal antigen in pregnancy, whereas NFAT-independent elements did not require fetal antigen. Despite using conserved molecular circuitry, Preg-TEX cells differed from TEX cells in chronic viral infection with respect to magnitude and dependency of T cell hypofunction on NFAT-independent signals. Altogether, these data reveal the molecular mechanisms and clinical consequences of maternal CD8 T cell hypofunction and identify pregnancy as a previously unappreciated context in which T cell exhaustion may occur.

Super interactive promoters provide insight into cell type-specific regulatory networks in blood lineage cell types

Existing studies of chromatin conformation have primarily focused on potential enhancers interacting with gene promoters. By contrast, the interactivity of promoters per se, while equally critical to understanding transcriptional control, has been largely unexplored, particularly in a cell type-specific manner for blood lineage cell types. In this study, we leverage promoter capture Hi-C data across a compendium of blood lineage cell types to identify and characterize cell type-specific super-interactive promoters (SIPs). Notably, promoter-interacting regions (PIRs) of SIPs are more likely to overlap with cell type-specific ATAC-seq peaks and GWAS variants for relevant blood cell traits than PIRs of non-SIPs. Moreover, PIRs of cell-type-specific SIPs show enriched heritability of relevant blood cell trait (s), and are more enriched with GWAS variants associated with blood cell traits compared to PIRs of non-SIPs. Further, SIP genes tend to express at a higher level in the corresponding cell type. Importantly, SIP subnetworks incorporating cell-type-specific SIPs and ATAC-seq peaks help interpret GWAS variants. Examples include GWAS variants associated with platelet count near the megakaryocyte SIP gene EPHB3 and variants associated lymphocyte count near the native CD4 T-Cell SIP gene ETS1. Interestingly, around 25.7% ~ 39.6% blood cell traits GWAS variants residing in SIP PIR regions disrupt transcription factor binding motifs. Importantly, our analysis shows the potential of using promoter-centric analyses of chromatin spatial organization data to identify biologically important genes and their regulatory regions.

Toll-like Receptor Signaling Inhibitory Peptide Improves Inflammation in Animal Model and Human Systemic Lupus Erythematosus

Toll-like receptors (TLRs) play a major role in the innate immune system. Several studies have shown the regulatory effects of TLR-mediated pathways on immune and inflammatory diseases. Dysregulated functions of TLRs within the endosomal compartment, including TLR7/9 trafficking, may cause systemic lupus erythematosus (SLE). TLR signaling pathways are fine-tuned by Toll/interleukin-1 receptor (TIR) domain-containing adapters, leading to interferon (IFN)-α production. This study describes a TLR inhibitor peptide 1 (TIP1) that primarily suppresses the downstream signaling mediated by TIR domain-containing adapters in an animal model of lupus and patients with SLE. The expression of most downstream proteins of the TLR7/9/myeloid differentiation factor 88 (MyD88)/IFN regulatory factor 7 signaling was downregulated in major tissues such as the kidney, spleen, and lymph nodes of treated mice. Furthermore, the pathological analysis of the kidney tissue confirmed that TIP1 could improve inflammation in MRL/lpr mice. TIP1 treatment downregulated many downstream proteins associated with TLR signaling, such as MyD88, interleukin-1 receptor-associated kinase, tumor necrosis factor receptor-associated factor 6, and IFN-α, in the peripheral blood mononuclear cells of patients with SLE. In conclusion, our data suggest that TIP1 can serve as a potential candidate for the treatment of SLE.

Association of microRNA-34a rs2666433 (A/G) Variant with Systemic Lupus Erythematosus in Female Patients: A Case-Control Study

Several microRNAs (miRNAs) are associated with autoimmune disease susceptibility and phenotype, including systemic lupus erythematosus (SLE). We aimed to explore for the first time the role of the miRNA-34a gene (MIR34A) rs2666433A > G variant in SLE risk and severity. A total of 163 adult patients with SLE and matched controls were recruited. Real-Time allelic discrimination PCR was applied for genotyping. Correlation with disease activity and clinic-laboratory data was done. The rs2666433 variant conferred protection against SLE development under heterozygous [A/G vs. G/G; OR = 0.57, 95%CI = 0.34-0.95], homozygous [A/A vs. G/G; OR = 0.52, 95%CI = 0.29-0.94], dominant [A/G + A/A vs. GG; OR = 0.55, 95%CI = 0.35-0.88], and log-additive [OR = 0.71, 95%CI = 0.53-0.95] models. Data stratification by sex revealed a significant association with SLE development in female participants under heterozygous/homozygous models (p-interaction = 0.004). There was no clear demarcation between SLE patients carrying different genotypes regarding the disease activity index or patients stratified according to lupus nephritis. Enrichment analysis confirmed the implication of MIR34A in the SLE pathway by targeting several genes related to SLE etiopathology. In conclusion, although the MIR34A rs2666433 variant conferred protection against developing SLE disease in the study population, it showed no association with disease activity. Replication studies in other populations are warranted.

Transcription Factors in Cartilage Homeostasis and Osteoarthritis

Osteoarthritis (OA) is the most common degenerative joint disease, and it is characterized by articular cartilage loss. In part, OA is caused by aberrant anabolic and catabolic activities of the chondrocyte, the only cell type present in cartilage. These chondrocyte activities depend on the intra- and extracellular signals that the cell receives and integrates into gene expression. The key proteins for this integration are transcription factors. A large number of transcription factors exist, and a better understanding of the transcription factors activated by the various signaling pathways active during OA can help us to better understand the complex etiology of OA. In addition, establishing such a profile can help to stratify patients in different subtypes, which can be a very useful approach towards personalized therapy. In this review, we discuss crucial transcription factors for extracellular matrix metabolism, chondrocyte hypertrophy, chondrocyte senescence, and autophagy in chondrocytes. In addition, we discuss how insight into these factors can be used for treatment purposes.

Co-expression of master transcription factors determines CD4+ T cell plasticity and functions in auto-inflammatory diseases

Master CD4⁺ T cell lineage determined transcription factors are found to be dysregulated in pathogenesis of autoimmune and inflammatory diseases. CD4⁺ T cells categorized into different lineages based on their functions, cell surface markers and master transcription factors those required for expression of lineage specific cytokines. T-bet, GATA3, RORγt and Foxp3 are major transcription regulators of Th1, Th2, Th17 and Treg cells respectively. Significant progress has been made in understanding expression of lineage specific master regulators that drives CD4⁺ T cell differentiation. It is known that each CD4⁺ T cell lineage express precise determined transcription factor and due to cross regulation between these factors the CD4⁺ T cells able to maintain thier specific phenotype. However, recent studies shows that the lineage specifying transcription factors frequently co-expressed. There is an emerging area of research revealing that the co-expression of lineage-specifying transcription factors alters the potential function and flexibility of subsets of CD4⁺ T cell, this in turn favors the autoimmune pathology. Here, we discuss similarities and differences between mutually co-expressed transcription factors in CD4⁺ T cell subsets and then recapitulates on cell type specific and dynamic balance between the lineage restricted transcription factors in determining plasticity of CD4⁺ T cell subsets. Furthermore, we discuss abnormal regulation of such transcription factors that establishes a pathogenic CD4⁺ T cell phenotype in autoimmune diseases and how this understanding will provide further insight into potential therapeutic development.

A Computational Model of the Endothelial to Mesenchymal Transition

Endothelial cells (ECs) form the lining of lymph and blood vessels. Changes in tissue requirements or wounds may cause ECs to behave as tip or stalk cells. Alternatively, they may differentiate into mesenchymal cells (MCs). These processes are known as EC activation and endothelial-to-mesenchymal transition (EndMT), respectively. EndMT, Tip, and Stalk EC behaviors all require SNAI1, SNAI2, and Matrix metallopeptidase (MMP) function. However, only EndMT inhibits the expression of VE-cadherin, PECAM1, and VEGFR2, and also leads to EC detachment. Physiologically, EndMT is involved in heart valve development, while a defective EndMT regulation is involved in the physiopathology of cardiovascular malformations, congenital heart disease, systemic and organ fibrosis, pulmonary arterial hypertension, and atherosclerosis. Therefore, the control of EndMT has many promising potential applications in regenerative medicine. Despite the fact that many molecular components involved in EC activation and EndMT have been characterized, the system-level molecular mechanisms involved in this process have not been elucidated. Toward this end, hereby we present Boolean network model of the molecular involved in the regulation of EC activation and EndMT. The simulated dynamic behavior of our model reaches fixed and cyclic patterns of activation that correspond to the expected EC and MC cell types and behaviors, recovering most of the specific effects of simple gain and loss-of-function mutations as well as the conditions associated with the progression of several diseases. Therefore, our model constitutes a theoretical framework that can be used to generate hypotheses and guide experimental inquiry to comprehend the regulatory mechanisms behind EndMT. Our main findings include that both the extracellular microevironment and the pattern of molecular activity within the cell regulate EndMT. EndMT requires a lack of VEGFA and sufficient oxygen in the extracellular microenvironment as well as no FLI1 and GATA2 activity within the cell. Additionally Tip cells cannot undergo EndMT directly. Furthermore, the specific conditions that are sufficient to trigger EndMT depend on the specific pattern of molecular activation within the cell.

Quantifying functional impact of non-coding variants with multi-task Bayesian neural network

MOTIVATION

Advances in high-throughput genotyping and sequencing technologies during recent years have revealed essential roles of non-coding regions in gene regulation. Genome-wide association studies (GWAS) suggested that a large proportion of risk variants are located in non-coding regions and remain unexplained by current expression quantitative trait loci catalogs. Interpreting the causal effects of these genetic modifications is crucial but difficult owing to our limited knowledge of how regulatory elements function. Although several computational methods have been designed to prioritize regulatory variants that substantially impact human phenotypes, few of them achieve consistently high performance even when large-scale multi-omic data are integrated.

RESULTS

We propose a novel multi-task framework based on Bayesian deep neural networks, MtBNN, to quantify the deleterious impact of single nucleotide polymorphisms in non-coding genomic regions. With the high-efficiency provided by the multi-task Bayesian framework to integrate information from different sources, MtBNN is capable of extracting features from genomic sequences of large-scale chromatin-profiling data, such as chromatin accessibility and transcript factor binding affinities, and calculating the distribution of the probability that a non-coding variant disrupts regulatory activities. A series of comprehensive experiments show that MtBNN quantifies the functional impact of cis-regulatory variations with high accuracy, including expression quantitative trait locus, DNase I sensitivity quantitative trait locus and functional genetic variants located within ATAC-peaks that affect the accessibility of the corresponding peak and achieves significantly better performance than the existing methods. Moreover, MtBNN has applications in the discovery of potentially causal disease-associated single-nucleotide polymorphisms (SNPs), thus helping fine-map the GWAS SNPs.

AVAILABILITY AND IMPLEMENTATION

Code can be downloaded from https://github.com/Zoesgithub/MtBNN.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Predicting flares in patients with stable systemic lupus erythematosus

OBJECTIVES

Data on flares in Asian patients with systemic lupus erythematosus (SLE) are scarce. Here, we aim to identify the baseline predictors of flares in a cohort of Southeast Asian patients with SLE.

METHODS

Consecutive adult patients with prevalent SLE according to the 1997 ACR or 2012 SLICC criteria were enrolled and followed three-monthly. Clinical and laboratory data were collected at every visit using a standardised protocol. Flares were defined using the SELENA-SLEDAI Flare Index (SFI). Baseline predictors of flare in patients with stable disease (SLE Disease Activity Index-2K (SLEDAI-2K) of ≤ 4) were determined using Cox proportional hazards.

RESULTS

Of the 210 patients recruited, 148 (70.5%) were Chinese. The median (IQR) SLEDAI-2K at entry was 2 (0-4) and the median (IQR) disease duration was 10 (4.4-16.4) years. At baseline, 152 (72.4%) patients had stable disease. After a median (IQR) follow-up of 31.5 (24.1-36.3) months, 109 (51.9%) flared. Stable patients who flared tended to be in the lowest tertile of age (HR 3.08, 95% CI 1.72-5.48, p < 0.01), had thrombocytopenia (HR 5.01, 95% CI 1.32-18.99, p = 0.02), hypocomplementemia (HR 3.35, 95% CI 1.54-7.30, p < 0.01) and had the highest baseline prednisolone doses (HR 2.39, 95% CI 1.28-4.46, p = 0.01). Conversely, patients in the lowest tertile of disease duration tended not to flare (HR 0.41, 95% CI 0.21-0.80, p = 0.01).

CONCLUSION

Flares are common in Asian SLE patients with initial stable disease. Close monitoring is needed for patients who are younger, with longer disease duration, thrombocytopenia, hypocomplementemia, or who required a higher baseline prednisolone dose.

Ets1 Controls the Development of B Cell Autoimmune Responses in a Cell-Intrinsic Manner

Ets1 is emerging as a key transcription factor that is required to prevent autoimmunity in mice and humans. Ets1 is expressed in both B and T cells, and mice lacking Ets1 are characterized by excess B and T cell activation, leading to enhanced formation of Ab-secreting cells and high titers of autoantibodies. In humans, genome-wide association studies have detected associations of single nucleotide polymorphisms in the human ETS1 gene with autoimmune diseases, including lupus. An increased fraction of CD4⁺ T cells from Ets1^−/− mice have an activated effector-memory phenotype, and there are aberrations in differentiation that contribute to the autoimmune phenotype. In vitro studies of B cells suggest that Ets1 may have B cell–intrinsic effects as well. To confirm B cell–intrinsic roles for Ets1, we crossed CD19-Cre mice to mice with a floxed allele of Ets1. Mice with a B cell–specific deletion of Ets1 show increases in B cell activation, numbers of Ab-secreting cells, and levels of autoantibodies, despite the fact that T cells are normal. However, when compared with conventional Ets1 knockout mice, mice with B cell–specific loss of Ets1 have a significantly milder phenotype. These results demonstrate that Ets1 is required in B cells to prevent autoimmune responses but that loss of Ets1 activity in other cell types is required for maximal autoimmune phenotypes.

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.

Extrafollicular responses in humans and SLE

Chronic autoimmune diseases, and in particular Systemic Lupus Erythematosus (SLE), are endowed with a long-standing autoreactive B-cell compartment that is presumed to reactivate periodically leading to the generation of new bursts of pathogenic antibody-secreting cells (ASC). Moreover, pathogenic autoantibodies are typically characterized by a high load of somatic hypermutation and in some cases are highly stable even in the context of prolonged B-cell depletion. Long-lived, highly mutated antibodies are typically generated through T-cell-dependent germinal center (GC) reactions. Accordingly, an important role for GC reactions in the generation of pathogenic autoreactivity has been postulated in SLE. Nevertheless, pathogenic autoantibodies and autoimmune disease can be generated through B-cell extrafollicular (EF) reactions in multiple mouse models and human SLE flares are characterized by the expansion of naive-derived activated effector B cells of extrafollicular phenotype. In this review, we will discuss the properties of the EF B-cell pathway, its relationship to other effector B-cell populations, its role in autoimmune diseases, and its contribution to human SLE. Furthermore, we discuss the relationship of EF B cells with Age-Associated B cells (ABCs), a TLR-7-driven B-cell population that mediates murine autoimmune and antiviral responses.

Genome-wide DNA methylation analysis in primary antiphospholipid syndrome neutrophils

Antiphospholipid syndrome (APS) is a systemic autoimmune disease characterized by thromboembolic events and pregnancy loss. We sought to characterize the DNA methylation profile of primary APS in comparison to healthy controls and individuals with SLE. In primary APS neutrophils compared to controls, 17 hypomethylated and 25 hypermethylated CpG sites were identified. Notable hypomethylated genes included ETS1, a genetic risk locus for SLE, and PTPN2, a genetic risk locus for other autoimmune diseases. Gene ontology analysis of hypomethylated genes revealed enrichment of genes involved in pregnancy. None of the differentially methylated sites in primary APS were differentially methylated in SLE neutrophils, and there was no demethylation of interferon signature genes in primary APS as is seen in SLE. Hypomethylation within a single probe in the IFI44L promoter (cg06872964) was able to distinguish SLE from primary APS with a sensitivity of 93.3% and specificity of 80.0% at a methylation fraction of 0.329.

Longitudinal expression profiling of CD4+ and CD8+ cells in patients with active to quiescent giant cell arteritis

Background

Giant cell arteritis (GCA) is the most common form of vasculitis affecting elderly people. It is one of the few true ophthalmic emergencies but symptoms and signs are variable thereby making it a challenging disease to diagnose. A temporal artery biopsy is the gold standard to confirm GCA, but there are currently no specific biochemical markers to aid diagnosis. We aimed to identify a less invasive method to confirm the diagnosis of GCA, as well as to ascertain clinically relevant predictive biomarkers by studying the transcriptome of purified peripheral CD4+ and CD8+ T lymphocytes in patients with GCA.

Methods

We recruited 16 patients with histological evidence of GCA at the Royal Victorian Eye and Ear Hospital, Melbourne, Australia, and aimed to collect blood samples at six time points: acute phase, 2–3 weeks, 6–8 weeks, 3 months, 6 months and 12 months after clinical diagnosis. CD4+ and CD8+ T-cells were positively selected at each time point through magnetic-assisted cell sorting. RNA was extracted from all 195 collected samples for subsequent RNA sequencing. The expression profiles of patients were compared to those of 16 age-matched controls.

Results

Over the 12-month study period, polynomial modelling analyses identified 179 and 4 statistically significant transcripts with altered expression profiles (FDR < 0.05) between cases and controls in CD4+ and CD8+ populations, respectively. In CD8+ cells, two transcripts remained differentially expressed after 12 months; SGTB, associated with neuronal apoptosis, and FCGR3A, associatied with Takayasu arteritis. We detected genes that correlate with both symptoms and biochemical markers used for predicting long-term prognosis. 15 genes were shared across 3 phenotypes in CD4 and 16 across CD8 cells. In CD8, IL32 was common to 5 phenotypes including Polymyalgia Rheumatica, bilateral blindness and death within 12 months.

Conclusions

This is the first longitudinal gene expression study undertaken to identify robust transcriptomic biomarkers of GCA. Our results show cell type-specific transcript expression profiles, novel gene-phenotype associations, and uncover important biological pathways for this disease. In the acute phase, the gene-phenotype relationships we have identified could provide insight to potential disease severity and as such guide in initiating appropriate patient management.

Electronic supplementary material

The online version of this article (10.1186/s12920-018-0376-4) contains supplementary material, which is available to authorized users.

Overall Downregulation of mRNAs and Enrichment of H3K4me3 Change Near Genome-Wide Association Study Signals in Systemic Lupus Erythematosus: Cell-Specific Effects

This study was designed to define gene expression and H3K4me3 histone modifications in T cells, B cells, and monocytes in systemic lupus erythematosus (SLE). Array studies of total peripheral blood mononuclear cells have demonstrated gene expression signatures related to neutrophils, interferon, and other inflammatory pathways. It is not clear how consistent these effects are across different cell types. In this study, RNA-seq and chromatin immunoprecipitation-seq were utilized to identify gene expression patterns and H3K4me3 histone modifications related to promoter activation in SLE. Across the three cell types, there was 55% concordance for gene expression changes related to SLE. Key conserved pathways were ribosome biogenesis among upregulated genes and heat shock response among downregulated genes. ETS family transcription factors (TFs) and STAT1 were revealed as common regulators by position weight matrices. When epigenetic changes were leveraged with gene expression, the pivotal TFs ATF3 and FOS were defined with ATF3 also cross-referencing with gene expression-identified TFs. Genome-wide association study (GWAS) single nucleotide polymorphisms associated with SLE were cross-referenced with both mRNA and H3K4me3 changes in SLE. Baseline mRNA expression and H3K4me3 peak height was higher at sites that cross-referenced with GWAS signals, however, all three cell types exhibited an overall decrease in expression of GWAS-associated RNAs differentially expressed in SLE. H3K4me3 changes in SLE were also enriched in GWAS-associated sites. In summary, the SLE disease process is associated with both shared and cell-specific changes in gene expression and epigenetics. Surprisingly, GWAS-associated RNAs were overall markedly decreased across all three cell types. TF analysis identified ATF3, FOS, STAT1, and ETS family members as critical, all pathways with a recognized relationship to the SLE disease process. GWAS signals clearly mark both cell-type specific changes in SLE as well as concordant changes across all three cell types. Interpretation of single nucleotide polymorphism effects in SLE will require tissue-specific mechanistic studies and therapeutics will require mechanistic studies in multiple cell types.

Bruton's Tyrosine Kinase, a Component of B Cell Signaling Pathways, Has Multiple Roles in the Pathogenesis of Lupus

Systemic lupus erythematosus (SLE) is an autoimmune disease characterized by the loss of adaptive immune tolerance to nucleic acid-containing antigens. The resulting autoantibodies form immune complexes that promote inflammation and tissue damage. Defining the signals that drive pathogenic autoantibody production is an important step in the development of more targeted therapeutic approaches for lupus, which is currently treated primarily with non-specific immunosuppression. Here, we review the contribution of Bruton’s tyrosine kinase (Btk), a component of B and myeloid cell signaling pathways, to disease in murine lupus models. Both gain- and loss-of-function genetic studies have revealed that Btk plays multiple roles in the production of autoantibodies. These include promoting the activation, plasma cell differentiation, and class switching of autoreactive B cells. Small molecule inhibitors of Btk are effective at reducing autoantibody levels, B cell activation, and kidney damage in several lupus models. These studies suggest that Btk may promote end-organ damage both by facilitating the production of autoantibodies and by mediating the inflammatory response of myeloid cells to these immune complexes. While Btk has not been associated with SLE in GWAS studies, SLE B cells display signaling defects in components both upstream and downstream of Btk consistent with enhanced activation of Btk signaling pathways. Taken together, these observations indicate that limiting Btk activity is critical for maintaining B cell tolerance and preventing the development of autoimmune disease. Btk inhibitors, generally well-tolerated and approved to treat B cell malignancy, may thus be a useful therapeutic approach for SLE.