Insufficient PD-1 expression during active autoimmune responses: a deep single-cell proteomics analysis in inflammatory arthritis.

Objectives: Programmed cell death protein-1 (PD-1) maintains peripheral immune tolerance by preventing T cell continuous activation. Aiming to understand the extent of PD-1 expression in inflammatory arthritis beyond its involvement with T cells, we assess its presence on various circulating single cells. Methods: Mass cytometry analysis of patients with active seropositive/seronegative rheumatoid (RA; n=9/8) and psoriatic (PsA; n=9) arthritis versus healthy controls (HC; n=13), re-evaluating patients after 3 months of anti-rheumatic treatment. Results: PD-1 was expressed in all leukocyte subpopulations, with the highest PD-1+ cell frequencies in eosinophils (59-73%) and T cells (50-60%), and the lowest in natural-killer cells (1-3%). PD-1+ cell frequencies and PD-1 median expression were comparable between patient subgroups and HC, in the majority of cell subsets. Exceptions included increases in certain T cell/B cell subsets of seropositive RA and specific monocyte subsets and dendritic cells of PsA; an expanded PD-1+CD4+CD45RA+CD27+CD28+ T subset, denoting exhausted T cells, was common across patient subgroups. Strikingly, significant inverse correlations between individual biomarkers of systemic inflammation (ESR and/or serum CRP) and PD-1+ cell frequencies and/or median expression were evident in several innate and adaptive immunity cell subsets of RA and PsA patients. Furthermore, all inverse correlations noted in individuals with active arthritis were no longer discernible in those who attained remission/low disease activity post-treatment. Conclusion: PD-1 expression may be insufficient, relative to the magnitude of the concomitant systemic inflammatory response on distinct leukocyte subsets, varying between RA and PsA. Our results point to the potential therapeutic benefits of pharmacological PD-1 activation, to rebalance the autoimmune response and reduce inflammation.

Single-cell analysis of MYD88L265P and MYD88WT Waldenström macroglobulinemia patients.

Waldenström macroglobulinemia (WM) is characterized by the expansion of clonal lymphoplasmacytic cells; the MYD88L265P somatic mutation is found in >90% of patients, but malignant B cells may still display intra-clonal heterogeneity. To assess clonal heterogeneity in WM, we generated and performed single-cell RNA sequencing of CD19+ sorted cells from five patients with MYD88 L265P and two patients with MYD88 WT genotype as well as two healthy donors. We identified distinct transcriptional patterns in the clonal subpopulations not only between the two genetically distinct WM subgroups but also among MYD88 L265P patients, which affected the B cell composition in the different subgroups. Comparison of clonal and normal/polyclonal B cells within each patient sample enabled the identification of patient-specific transcriptional changes. We identified gene signatures active in a subset of MYD88L265P patients, while other signatures were active in MYD88 WT patients. Finally, gene expression analysis showed common transcriptional features between patients compared to the healthy control but also differentially expressed genes between MYD88 L265P and MYD88 WT patients involved in distinct pathways, including NFκΒ, BCL2, and BTK. Overall, our data highlight the intra-tumor clonal heterogeneity in WM with potential prognostic and therapeutic implications.

Distinct innate and adaptive immunity phenotypic profile at the circulating single-cell level in Psoriatic Arthritis.

Mass cytometry was employed to investigate 47 circulating leukocyte subsets in patients with active psoriatic arthritis (PsA, n = 16) compared to healthy controls (n = 13), seropositive (RF and/or anti-CCP, n = 12) and seronegative (n = 9) RA patients. Comparing PsA to controls, different cell frequencies were found in both innate and adaptive immunity cell subsets, as well as in cells bridging innate and adaptive immunity. In some T cell subsets increased costimulatory molecules' expression in PsA, was also noted.No changes were observed in patients who remained disease-active after 3 months of treatment, in contrast to those who achieved remission/low-disease activity. Comparing PsA to seropositive RA, elevated frequencies of naïve and activated CD8+ T cells, B cells, MAIT/iNKT and ILCs were found, while the opposite was the case for terminal effector, senescent, and Th2-like cells. Strikingly, the composition of the leukocyte pool in PsA was comparable to seronegative RA, providing evidence for the pathogenetic similarities between these two entities.

SARS-CoV-2 mRNA Dual Immunization Induces Innate Transcriptional Signatures, Establishes T-Cell Memory and Coordinates the Recall Response.

BACKGROUND: mRNA vaccines have played a crucial role in controlling the SARS-CoV-2 global pandemic. However, the immunological mechanisms involved in the induction, magnitude and longevity of mRNA-vaccine-induced protective immunity are still unclear. METHODS: In our study, we used whole-RNA sequencing along with detailed immunophenotyping of antigen-specific T cells and humoral RBD-specific response to dual immunization with the Pfizer-BioNTech mRNA vaccine (BNT162b2) and correlated them with response to an additional dose, administered 10 months later, in order to comprehensively profile the immune response of healthy volunteers to BNT162b2. RESULTS: Primary dual immunization induced upregulation of the Type I interferon pathway and generated spike protein (S)-specific IFN-γ+ and TNF-α+ CD4 T cells, S-specific memory CD4 T cells, and RBD-specific antibodies against SARS-CoV-2. S-specific CD4 T cells induced by the primary series correlated with the RBD-specific antibody titers to a third dose. CONCLUSIONS: This study demonstrates the induction of both innate and adaptive immunity in response to the BNT162b2 mRNA vaccine in a coordinated manner and identifies the central role of primarily induced CD4+ T cells as a predictive biomarker of the magnitude of anamnestic immune response.

Whole blood transcriptome identifies interferon-regulated genes as key drivers in thrombotic primary antiphospholipid syndrome.

OBJECTIVE: Pathogenesis of antiphospholipid syndrome (APS) isn't fully elucidated. We aimed to identify gene signatures characterizing thrombotic primary APS (thrPAPS) and subgroups at high risk for worse outcomes. METHODS: We performed whole blood next-generation RNA-sequencing in 62 patients with thrPAPS and 29 age-/sex-matched healthy controls (HCs), followed by differential gene expression analysis (DGEA) and enrichment analysis. We trained models on transcriptomics data using machine learning. RESULTS: DGEA of 12.306 genes revealed 34 deregulated genes in thrPAPS versus HCs; 33 were upregulated by at least 2-fold, and 14/33 were type I and II interferon-regulated genes (IRGs) as determined by interferome database. Machine learning applied to deregulated genes returned 79% accuracy to discriminate thrPAPS from HCs, which increased to 82% when only the most informative IRGs were analyzed. Comparison of thrPAPS subgroups versus HCs showed an increased presence of IRGs among upregulated genes in venous thrombosis (21/23, 91%), triple-antiphospholipid antibody (aPL) positive (30/50, 60%), and recurrent thrombosis (19/42, 45%) subgroups. Enrichment analysis of upregulated genes in triple-aPL positive patients revealed terms related to 'type I interferon signaling pathway' and 'innate immune response'. DGEA among thrPAPS subgroups revealed upregulated genes, including IRGs, in patients with venous versus arterial thrombosis (n = 11, 9 IRGs), triple-aPL versus non-triple aPL (n = 10, 9 IRGs), and recurrent versus non-recurrent thrombosis (n = 10, 3 IRGs). CONCLUSION: Upregulated IRGs may better discriminate thrPAPS from HCs than all deregulated genes in peripheral blood. Taken together with DGEA data, IRGs are highly expressed in thrPAPS and high-risk subgroups of triple-aPL and recurrent thrombosis, with potential treatment implications.

Deregulated DNA damage response network in Behcet's disease.

Behcet's disease (BD) is a chronic, relapsing systemic vasculitis of unknown etiology. Since the DNA repair enzyme NEIL1 has been identified as one of the two genetic risk factors for BD by whole exome study, we examined the potential involvement of the DNA damage response (DDR) network in BD. Peripheral blood mononuclear cells from 26 patients and 26 age-/sex-matched healthy controls were studied. Endogenous DNA damage levels were increased in active BD patients compared to controls or patients in remission. In parallel, BD patients had defective nucleotide excision repair capacity. RNA-sequencing revealed reduced expression of NEIL1 that negatively correlated with DNA damage accumulation. On the other hand, expression of genes involved in senescence and senescence-associated secretory phenotype positively correlated with individual endogenous DNA damage levels. We conclude that deregulated DDR contributes to the proinflammatory environment in BD.

Protocol for unbiased, consolidated variant calling from whole exome sequencing data.

Whole Exome Sequencing (WES) is used for querying DNA variants using the protein coding parts of genomes (exomes). However, WES analysis can be challenging because of the complexity of the data. Here, we describe a consolidated protocol for unbiased WES analysis. The protocol uses three variant callers (HaplotypeCaller, FreeBayes, and DeepVariant), which have different underlying models. We provide detailed execution steps, as well as basic variant filtering, annotation, visualization, and consolidation aspects.

Lung tumor MHCII immunity depends on in situ antigen presentation by fibroblasts.

A key unknown of the functional space in tumor immunity is whether CD4 T cells depend on intratumoral MHCII cancer antigen recognition. MHCII-expressing, antigen-presenting cancer-associated fibroblasts (apCAFs) have been found in breast and pancreatic tumors and are considered to be immunosuppressive. This analysis shows that antigen-presenting fibroblasts are frequent in human lung non-small cell carcinomas, where they seem to actively promote rather than suppress MHCII immunity. Lung apCAFs directly activated the TCRs of effector CD4 T cells and at the same time produced C1q, which acted on T cell C1qbp to rescue them from apoptosis. Fibroblast-specific MHCII or C1q deletion impaired CD4 T cell immunity and accelerated tumor growth, while inducing C1qbp in adoptively transferred CD4 T cells expanded their numbers and reduced tumors. Collectively, we have characterized in the lungs a subset of antigen-presenting fibroblasts with tumor-suppressive properties and propose that cancer immunotherapies might be strongly dependent on in situ MHCII antigen presentation.

Blood Transcriptomes of Anti-SARS-CoV-2 Antibody-Positive Healthy Individuals Who Experienced Asymptomatic Versus Clinical Infection.

The reasons behind the clinical variability of SARS-CoV-2 infection, ranging from asymptomatic infection to lethal disease, are still unclear. We performed genome-wide transcriptional whole-blood RNA sequencing, bioinformatics analysis and PCR validation to test the hypothesis that immune response-related gene signatures reflecting baseline may differ between healthy individuals, with an equally robust antibody response, who experienced an entirely asymptomatic (n=17) versus clinical SARS-CoV-2 infection (n=15) in the past months (mean of 14 weeks). Among 12.789 protein-coding genes analysed, we identified six and nine genes with significantly decreased or increased expression, respectively, in those with prior asymptomatic infection relatively to those with clinical infection. All six genes with decreased expression (IFIT3, IFI44L, RSAD2, FOLR3, PI3, ALOX15), are involved in innate immune response while the first two are interferon-induced proteins. Among genes with increased expression six are involved in immune response (GZMH, CLEC1B, CLEC12A), viral mRNA translation (GCAT), energy metabolism (CACNA2D2) and oxidative stress response (ENC1). Notably, 8/15 differentially expressed genes are regulated by interferons. Our results suggest that subtle differences at baseline expression of innate immunity-related genes may be associated with an asymptomatic disease course in SARS-CoV-2 infection. Whether a certain gene signature predicts, or not, those who will develop a more efficient immune response upon exposure to SARS-CoV-2, with implications for prioritization for vaccination, warrant further study.

Interferon (IFN)-stimulated gene 15: A novel biomarker for lymphoma development in Sjögren's syndrome.

OBJECTIVE: We investigated whether interferon (IFN) induced genes could serve as biomarkers for the detection of lymphoma development among patients with Sjögren's syndrome (SS). METHODS: Total RNA was extracted from 98 labial minor salivary glands (LMSG) biopsies of SS patients [61 not complicated by lymphoma (SS-nL) and 37 complicated by Non-Hodgkin Lymphoma (NHL) (SS-L)] and 67 matched peripheral blood (PB) samples, as well as from 30 LMSG biopsies and 17 matched PB derived from sicca controls (SC). RNA sequencing was performed in LMSG biopsies of high and low risk SS patients for lymphoma development and SC. Expression analysis of type I (MX-1, IFIT-1, IFI44 and ISG-15) and type II IFN induced (CXCL9/MIG-1, GBP-1) genes was performed by real time PCR. RESULTS: ISG-15 transcript levels were significantly higher in SS-L patients compared to SS-nL patients in both LMSG tissues and PB specimens. Additionally, MIG-1 was found to display higher expression values in LMSG tissues, but not in PB derived from SS-L patients compared to the SS-nL group. A coordinate expression in PB/LMSG of type I IFN (ISG-15, MX-1 and IFI44), but not type II IFN induced genes was also observed. CONCLUSION: ISG-15 gene expression was able to distinguish SS-nL and SS-L at both periphery and tissue level and therefore could represent a novel biomarker for lymphoma development among SS patients. PB and LSMG seem to share a common transcriptional profile of type I IFN pathway.

Molecular Epidemiology of SARS-CoV-2 in Greece Reveals Low Rates of Onward Virus Transmission after Lifting of Travel Restrictions Based on Risk Assessment during Summer 2020.

The novel coronavirus severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spread rapidly during the first months of 2020 and continues to expand in multiple areas across the globe. Molecular epidemiology has provided an added value to traditional public health tools by identifying SARS-CoV-2 clusters or providing evidence that clusters based on virus sequences and contact tracing are highly concordant. Our aim was to infer the levels of virus importation and to estimate the impact of public health measures related to travel restrictions to local transmission in Greece. Our phylogenetic and phylogeographic analyses included 389 full-genome SARS-CoV-2 sequences collected during the first 7 months of the pandemic in Greece and a random collection in five replicates of 3,000 sequences sampled globally, as well as the best hits to our data set identified by BLAST. Phylogenetic trees were reconstructed by the maximum likelihood method, and the putative source of SARS-CoV-2 infections was inferred by phylogeographic analysis. Phylogenetic analyses revealed the presence of 89 genetically distinct viruses identified as independent introductions into Greece. The proportion of imported strains was 41%, 11.5%, and 8.8% during the three periods of sampling, namely, March (no travel restrictions), April to June (strict travel restrictions), and July to September (lifting of travel restrictions based on thorough risk assessment), respectively. The results of phylogeographic analysis were confirmed by a Bayesian approach. Our findings reveal low levels of onward transmission from imported cases during summer and underscore the importance of targeted public health measures that can increase the safety of international travel during a pandemic. IMPORTANCE Our study based on current state-of-the-art molecular epidemiology methods suggests that virus screening and public health measures after the lifting of travel restrictions prevented SARS-CoV-2 onward transmission from imported cases during summer 2020 in Greece. These findings provide important data on the efficacy of targeted public health measures and have important implications regarding the safety of international travel during a pandemic. Our results can provide a roadmap about prevention policy in the future regarding the reopening of borders in the presence of differences in vaccination coverage, the circulation of the virus, and the presence of newly emergent variants across the globe.

The role of A-to-I RNA editing in infections by RNA viruses: Possible implications for SARS-CoV-2 infection.

RNA editing is a fundamental biological process with 2 major forms, namely adenosine-to-inosine (A-to-I, recognized as A-to-G) and cytosine-to-uracil (C-to-U) deamination, mediated by ADAR and APOBEC enzyme families, respectively. A-to-I RNA editing has been shown to directly affect the genome/transcriptome of RNA viruses with significant repercussions for viral protein synthesis, proliferation and infectivity, while it also affects recognition of double-stranded RNAs by cytosolic receptors controlling the host innate immune response. Recent evidence suggests that RNA editing may be present in SARS-CoV-2 genome/transcriptome. The majority of mapped mutations in SARS-CoV-2 genome are A-to-G/U-to-C(opposite strand) and C-to-U/G-to-A(opposite strand) substitutions comprising potential ADAR-/APOBEC-mediated deamination events. A single nucleotide substitution can have dramatic effects on SARS-CoV-2 infectivity as shown by the D614G(A-to-G) substitution in the spike protein. Future studies utilizing serial sampling from patients with COVID-19 are warranted to delineate whether RNA editing affects viral replication and/or the host immune response to SARS-CoV-2.

Distinct transcriptional profile of blood mononuclear cells in Behçet's disease: insights into the central role of neutrophil chemotaxis.

OBJECTIVES: Both innate and adaptive immune responses are reportedly increased in Behçet's disease (BD), a chronic, relapsing systemic vasculitis lying at the intersection between autoinflammation and autoimmunity. To further study pathophysiologic molecular mechanisms operating in BD, we searched for transcriptome-wide changes in blood mononuclear cells from these patients. METHODS: We performed 3' mRNA next-generation sequencing-based genome-wide transcriptional profiling followed by analysis of differential expression signatures, Kyoto Encyclopedia of Genes and Genomes pathways, GO biological processes and transcription factor signatures. RESULTS: Differential expression analysis clustered the transcriptomes of 13 patients and one healthy subject separately from those of 10 healthy age/gender-matched controls and one patient. Among the total of 17 591 expressed protein-coding genes, 209 and 31 genes were significantly upregulated and downregulated, respectively, in BD vs controls by at least 2-fold. The most upregulated genes comprised an abundance of CC- and CXC-chemokines. Remarkably, the 5 out of top 10 upregulated biological processes involved leucocyte recruitment to peripheral tissues, especially for neutrophils. Moreover, NF-kB, TNF and IL-1 signalling pathways were prominently enhanced in BD, while transcription factor activity analysis suggested that the NF-kB p65/RELA subunit action underlies the observed differences in the BD transcriptome. CONCLUSION: This RNA-sequencing analysis in peripheral blood mononuclear cells derived from patients with BD does not support a major pathogenetic role for adaptive immunity-driven mechanisms, but clearly points to the action of aberrant innate immune responses with a central role played by upregulated neutrophil chemotaxis.

Learning Pathway Dynamics from Single-Cell Proteomic Data: A Comparative Study.

Single-cell platforms provide statistically large samples of snapshot observations capable of resolving intrercellular heterogeneity. Currently, there is a growing literature on algorithms that exploit this attribute in order to infer the trajectory of biological mechanisms, such as cell proliferation and differentiation. Despite the efforts, the trajectory inference methodology has not yet been used for addressing the challenging problem of learning the dynamics of protein signaling systems. In this work, we assess this prospect by testing the performance of this class of algorithms on four proteomic temporal datasets. To evaluate the learning quality, we design new general-purpose evaluation metrics that are able to quantify performance on (i) the biological meaning of the output, (ii) the consistency of the inferred trajectory, (iii) the algorithm robustness, (iv) the correlation of the learning output with the initial dataset, and (v) the roughness of the cell parameter levels though the inferred trajectory. We show that experimental time alone is insufficient to provide knowledge about the order of proteins during signal transduction. Accordingly, we show that the inferred trajectories provide richer information about the underlying dynamics. We learn that established methods tested on high-dimensional data with small sample size, slow dynamics, and complex structures (e.g. bifurcations) cannot always work in the signaling setting. Among the methods we evaluate, Scorpius and a newly introduced approach that combines Diffusion Maps and Principal Curves were found to perform adequately in recovering the progression of signal transduction although their performance on some metrics varies from one dataset to another. The novel metrics we devise highlight that it is difficult to conclude, which one method is universally applicable for the task. Arguably, there are still many challenges and open problems to resolve. © 2020 The Authors. Cytometry Part A published by Wiley Periodicals, Inc. on behalf of International Society for Advancement of Cytometry.

Phosphorylation mapping of Laminin ß1-chain: Kinases in association with active sites.

Laminins are a major constituent of the extracellular matrix (ECM). Laminin-111, the most extensively studied laminin isoform, consists of the α1, the ß1 and the γ1 chain, and is involved in many cellular processes, like adhesion, migration and differentiation. Given the regulatory role of phosphorylation in protein function, it is important to identify the phosphorylation sites of human laminin ß1-chain sequence (LAMB1). Therefore, we computationally predicted all possible phosphorylation sites in LAMB1. For the first time, we identified the possibly responsible kinases for already in vitro experimentally observed phosphorylated residues in LAMB1. All known functional (active) sites of LAMB1, were recorded after an extensive literature search and combined with the experimentally observed and our predicted phosphorylated residues. This generated a detailed phosphorylation map of LAMB1. Five kinases (PKA, PKC, CKII, CKI and GPCR1) were indicated important, while the role of PKA, PKC and CKII, kinases known for ectophosphorylation activity, was highlighted. The activity of PKA and PKC was associated with the active site RIQNLLKITNLRIKFVKLHTLGDNLLDS. Also, predicted phosphorylations inside two amyloidogenic (DSITKYFQMSLE, VILQHSAADIAR) and two anti-cancerous (YIGSR and PDSGR) sites suggested a possible role in the development of the corresponding diseases.

Phosphorylation mapping of laminin α1-chain: Kinases in association with active sites.

Laminin-111 is a trimeric glycoprotein of the extracellular matrix (ECM) that holds a significant role in cell adhesion, migration and differentiation. Laminin-111 is the most studied laminin isoform, composed of three chains; α1, ß1 and γ1. Phosphorylation is the most common eukaryotic post - translational modification and has regulatory effect on protein function. Using bioinformatic tools we computationally predicted all the possible phosphorylation sites on human laminin α1-chain sequence (LAMA1) according to kinases binding motifs. Thus, we predicted, for the first time, the possibly responsible kinases for fifteen of the nineteen already published experimentally observed phosphorylated residues in LAMA1. Searching the literature extensively, we recorded all the known functional sites (active sites) in LAMA1. We combined the experimentally observed and predicted phosphorylated residues as well as the active sites in LAMA1, generating an analytic phosphorylation map of human laminin α1-chain, which is useful for further analysis. Our results indicated fourteen kinases that might be important for the phosphorylation of human laminin α1-chain, out of which three kinases with reported ecto-phosphorylation activity (PKA, PKC and CKII) were suggested to have a more significant role. Six cancer associated-active sites were correlated with kinases, three out which were correlated with only the above ecto - kinases.