Objectives Rheumatoid arthritis (RA), systemic lupus erythematosus (SLE) and primary Sjögren’s syndrome (pSS) share many clinical manifestations and serological features. The aim of this study was to identify the common transcriptional profiling and composition of immune cells in peripheral blood in these autoimmune diseases (ADs). Methods We analysed bulk RNA-seq data for enrichment of biological processes, transcription factors (TFs) and deconvolution-based immune cell types from peripheral blood mononuclear cells (PBMCs) in 119 treatment-naive patients (41 RA, 38 pSS, 28 SLE and 12 polyautoimmunity) and 20 healthy controls. The single-cell RNA-seq (scRNA-seq) and flow cytometry had been performed to further define the immune cell subsets on PBMCs. Results Similar transcriptional profiles and common gene expression signatures associated with nucleosome assembly and haemostasis were identified across RA, SLE, pSS and polyautoimmunity. Distinct TF ensembles and gene regulatory network were mainly enriched in haematopoiesis. The upregulated cell-lineage-specific TFs PBX1, GATA1, TAL1 and GFI1B demonstrated a strong gene expression signature of megakaryocyte (MK) expansion. Gene expression-based cell type enrichment revealed elevated MK composition, specifically, CD41b⁺CD42b⁺ and CD41b⁺CD61⁺ MKs were expanded, further confirmed by flow cytometry in these ADs. In scRNA-seq data, MKs were defined by TFs PBX1/GATA1/TAL1 and pre-T-cell antigen receptor gene, PTCRA. Cellular heterogeneity and a distinct immune subpopulation with functional enrichment of antigen presentation were observed in MKs. Conclusions The identification of MK expansion provided new insights into the peripheral immune cell atlas across RA, SLE, pSS and polyautoimmunity. Aberrant regulation of the MK expansion might contribute to the pathogenesis of these ADs.

Objectives Familial aggregation of primary Sjögren’s syndrome (pSS), systemic lupus erythematosus (SLE), rheumatoid arthritis (RA) and co-aggregation of these autoimmune diseases (ADs) (also called familial autoimmunity) is well recognised. However, the genetic predisposition variants that explain this clustering remains poorly defined. Methods We used whole-exome sequencing on 31 families (9 pSS, 11 SLE, 6 RA and 5 mixed autoimmunity), followed by heterozygous filtering and cosegregation analysis of a family-focused approach to document rare variants predicted to be pathogenic by in silico analysis. Potential importance in immune-related processes, gene ontology, pathway enrichment and overlap analyses were performed to prioritise gene sets. Results A range from 1 to 50 rare possible pathogenic variants, including 39 variants in immune-related genes across SLE, RA and pSS families, were identified. Among this gene set, regulation of T cell activation (p=4.06×10⁻⁷) and T cell receptor (TCR) signalling pathway (p=1.73×10⁻⁶) were particularly concentrated, including PTPRC (CD45), LCK, LAT–SLP76 complex genes (THEMIS, LAT, ITK, TEC, TESPA1, PLCL1), DGKD, PRKD1, PAK2 and NFAT5, shared across 14 SLE, RA and pSS families. TCR-interactive genes P2RX7, LAG3, PTPN3 and LAX1 were also detected. Overlap analysis demonstrated that the antiviral immunity gene DUS2 variant cosegregated with SLE, RA and pSS phenotypes in an extended family, that variants in the TCR-pathway genes CD45, LCK and PRKD1 occurred independently in three mixed autoimmunity families, and that variants in CD36 and VWA8 occurred in both RA-pSS and SLE-pSS families. Conclusions Our preliminary results define common genetic characteristics linked to familial pSS, SLE and RA and highlight rare genetic variations in TCR signalling pathway genes which might provide innovative molecular targets for therapeutic interventions for those three ADs.

Rosacea is a chronic inflammatory skin disorder with high incidence rate. Although genetic predisposition to rosacea is suggested by existing evidence, the genetic basis remains largely unknown. Here we present the integrated results of whole genome sequencing (WGS) in 3 large rosacea families and whole exome sequencing (WES) in 49 additional validation families. We identify single rare deleterious variants of LRRC4, SH3PXD2A and SLC26A8 in large families, respectively. The relevance of SH3PXD2A, SLC26A8 and LRR family genes in rosacea predisposition is underscored by presence of additional variants in independent families. Gene ontology analysis suggests that these genes encode proteins taking part in neural synaptic processes and cell adhesion. In vitro functional analysis shows that mutations in LRRC4, SH3PXD2A and SLC26A8 induce the production of vasoactive neuropeptides in human neural cells. In a mouse model recapitulating a recurrent Lrrc4 mutation from human patients, we find rosacea-like skin inflammation, underpinned by excessive vasoactive intestinal peptide (VIP) release by peripheral neurons. These findings strongly support familial inheritance and neurogenic inflammation in rosacea development and provide mechanistic insight into the etiopathogenesis of the condition.

Rosacea is a common chronic inflammatory condition that mainly affects the central face. However, the molecular background of the normal central face and the transcriptional profiling and immune cell composition of rosacea lesions remain largely unknown. Here, we performed whole-skin and epidermal RNA-seq of central facial skin from healthy individuals, lesions and matched normal skin from rosacea patients. From whole-skin RNA-seq, the site-specific gene signatures for central facial skin were mainly enriched in epithelial cell differentiation, with upregulation of the activator protein-1 (AP1) transcription factor (TF). We identified the common upregulated inflammatory signatures and diminished keratinization signature for rosacea lesions. Gene ontology, pathway, TF enrichment and immunohistochemistry results suggested that STAT1 was the potential core of the critical TF networks connecting the epithelial–immune crosstalk in rosacea lesions. Epidermal RNA-seq and immunohistochemistry analysis further validated the epithelial-derived STAT1 signature in rosacea lesions. The epidermal STAT1/IRF1 signature was observed across ETR, PPR, and PhR subtypes. Immune cell composition revealed that macrophages were common in all 3 subtypes. Finally, we described subtype-specific gene signatures and immune cell composition correlated with phenotypes. These findings reveal the specific epithelial differentiation in normal central facial skin, and epithelial–immune crosstalk in lesions providing insight into an initial keratinocyte pattern in the pathogenesis of rosacea.

Objective: Idiopathic inflammatory myositis-associated interstitial lung disease (IIM-ILD) significantly increases morbidity and mortality. Lung ultrasound B-lines and Krebs von den Lungen-6 (KL-6) are identified as new sonographic and serum markers of ILD, respectively. The aim of our work was to assess the role of B-lines and KL-6 as markers of the severity of IIM-ILD. For this purpose, the correlation among B-lines score, serum KL-6 levels, high-resolution CT (HRCT) score, and pulmonary function tests were investigated in IIM-ILD patients. Methods: Thirty-eight patients with IIM-ILD underwent chest HRCT scans, lung ultrasound and pulmonary function tests (independently performed within 1 week) examination. To assess severity and extent of ILD at HRCT, the Warrick score was used. The B-lines score denoting the extension of ILD was calculated by summing the number of B-lines on a total of 50 scanning sites. Serum KL-6 levels (U/ml) was measured by chemiluminescent enzyme immunoassay. Results: A significant correlation was found between the B-lines score and serum KL-6 levels (r = 0.43, P < 0.01), and between the Warrick score and serum KL-6 levels (r = 0.45, P < 0.01). A positive correlation between B-lines score and the Warrick score (r = 0.87, P < 0.0001) was also confirmed. Both B-lines score and KL-6 levels inversely correlated to diffusion capacity for carbon monoxide (r = −0.77, P < 0.0001 and r = −0.42, P < 0.05, respectively) and total lung capacity (r = −0.73, P < 0.0001 and r = −0.36, P < 0.05, respectively). Moreover, B-lines correlated inversely with forced vital capacity (r = −0.73, P < 0.0001), forced expiratory volume in 1 s (r = −0.69, P < 0.0001). Conclusion: B-lines score and serum KL-6 levels correlate with HRCT findings and pulmonary function tests, supporting their use as measures of IIM-ILD severity.