High levels of SOX5 decrease proliferative capacity of human B cells, but permit plasmablast differentiation

Enhanced interpretation of immune cell phenotype and function through a rhesus macaque single-cell atlas

Single-cell RNA sequencing (scRNA-seq) allows cell classification using genome-wide transcriptional state; however, high-dimensional transcriptomic profiles, and the unsupervised analyses employed to interpret them, provide a systematically different view of biology than well-established functional/lineage definitions of immunocytes. Understanding these differences and limits is essential for accurate interpretation of these rich data. We present the Rhesus Immune Reference Atlas (RIRA), the first immune-focused macaque single-cell multi-tissue atlas. We contrasted transcriptional profiles against immune lineages, using surface protein and marker genes as ground truth. While the pattern of clustering can align with cell type, this is not always true. Especially within T and natural killer (NK) cells, many functionally distinct subsets lack defining markers, and strong shared expression programs, such as cytotoxicity, result in systematic intermingling by unsupervised clustering. We identified gene programs with high discriminatory/diagnostic value, including multi-gene signatures that model T/NK cell maturation. Directly measuring these diagnostic programs complements unsupervised analyses.

Single cell transcriptional analysis of human adenoids identifies molecular features of airway microfold cells

The nasal, oropharyngeal, and bronchial mucosa are primary contact points for airborne pathogens like Mycobacterium tuberculosis (Mtb), SARS-CoV-2, and influenza virus. While mucosal surfaces can function as both entry points and barriers to infection, mucosa-associated lymphoid tissues (MALT) facilitate early immune responses to mucosal antigens. MALT contains a variety of specialized epithelial cells, including a rare cell type called a microfold cell (M cell) that functions to transport apical antigens to basolateral antigen-presenting cells, a crucial step in the initiation of mucosal immunity. M cells have been extensively characterized in the gastrointestinal (GI) tract in murine and human models. However, the precise development and functions of human airway M cells is unknown. Here, using single-nucleus RNA sequencing (snRNA-seq), we generated an atlas of cells from the human adenoid and identified 16 unique cell types representing basal, club, hillock, and hematopoietic lineages, defined their developmental trajectories, and determined cell-cell relationships. Using trajectory analysis, we found that human airway M cells develop from progenitor club cells and express a gene signature distinct from intestinal M cells. Surprisingly, we also identified a heretofore unknown epithelial cell type demonstrating a robust interferon-stimulated gene signature. Our analysis of human adenoid cells enhances our understanding of mucosal immune responses and the role of M cells in airway immunity. This work also provides a resource for understanding early interactions of pathogens with airway mucosa and a platform for development of mucosal vaccines.

Single-cell RNA-seq mapping of chicken peripheral blood leukocytes

BACKGROUND

Single-cell transcriptomics provides means to study cell populations at the level of individual cells. In leukocyte biology this approach could potentially aid the identification of subpopulations and functions without the need to develop species-specific reagents. The present study aimed to evaluate single-cell RNA-seq as a tool for identification of chicken peripheral blood leukocytes. For this purpose, purified and thrombocyte depleted leukocytes from 4 clinically healthy hens were subjected to single-cell 3' RNA-seq. Bioinformatic analysis of data comprised unsupervised clustering of the cells, and annotation of clusters based on expression profiles. Immunofluorescence phenotyping of the cell preparations used was also performed.

RESULTS

Computational analysis identified 31 initial cell clusters and based on expression of defined marker genes 28 cluster were identified as comprising mainly B-cells, T-cells, monocytes, thrombocytes and red blood cells. Of the remaining clusters, two were putatively identified as basophils and eosinophils, and one as proliferating cells of mixed origin. In depth analysis on gene expression profiles within and between the initial cell clusters allowed further identification of cell identity and possible functions for some of them. For example, analysis of the group of monocyte clusters revealed subclusters comprising heterophils, as well as putative monocyte subtypes. Also, novel aspects of TCRγ/δ + T-cell subpopulations could be inferred such as evidence of at least two subtypes based on e.g., different expression of transcription factors MAF, SOX13 and GATA3. Moreover, a novel subpopulation of chicken peripheral B-cells with high SOX5 expression was identified. An overall good correlation between mRNA and cell surface phenotypic cell identification was shown.

CONCLUSIONS

Taken together, we were able to identify and infer functional aspects of both previously well known as well as novel chicken leukocyte populations although some cell types. e.g., T-cell subtypes, proved more challenging to decipher. Although this methodology to some extent is limited by incomplete annotation of the chicken genome, it definitively has benefits in chicken immunology by expanding the options to distinguish identity and functions of immune cells also without access to species specific reagents.

Genome-wide CRISPR activation screening in senescent cells reveals SOX5 as a driver and therapeutic target of rejuvenation

Our understanding of the molecular basis for cellular senescence remains incomplete, limiting the development of strategies to ameliorate age-related pathologies by preventing stem cell senescence. Here, we performed a genome-wide CRISPR activation (CRISPRa) screening using a human mesenchymal precursor cell (hMPC) model of the progeroid syndrome. We evaluated targets whose activation antagonizes cellular senescence, among which SOX5 outperformed as a top hit. Through decoding the epigenomic landscapes remodeled by overexpressing SOX5, we uncovered its role in resetting the transcription network for geroprotective genes, including HMGB2. Mechanistically, SOX5 binding elevated the enhancer activity of HMGB2 with increased levels of H3K27ac and H3K4me1, raising HMGB2 expression so as to promote rejuvenation. Furthermore, gene therapy with lentiviruses carrying SOX5 or HMGB2 rejuvenated cartilage and alleviated osteoarthritis in aged mice. Our study generated a comprehensive list of rejuvenators, pinpointing SOX5 as a potent driver for rejuvenation both in vitro and in vivo.

Double-negative-2 B cells are the major synovial plasma cell precursor in rheumatoid arthritis

B cells are key pathogenic drivers of chronic inflammation in rheumatoid arthritis (RA). There is limited understanding of the relationship between synovial B cell subsets and pathogenic antibody secreting cells (ASCs). This knowledge is crucial for the development of more targeted B-cell depleting therapies. While CD11c+ double-negative 2 (DN2) B cells have been suggested as an ASC precursor in lupus, to date there is no proven link between the two subsets in RA. We have used both single-cell gene expression and BCR sequencing to study synovial B cells from patients with established RA, in addition to flow cytometry of circulating B cells. To better understand the differentiation patterns within the diseased tissue, a combination of RNA-based trajectory inference and clonal lineage analysis of BCR relationships were used. Both forms of analysis indicated that DN2 B cells serve as a major precursors to synovial ASCs. This study advances our understanding of B cells in RA and reveals the origin of pathogenic ASCs in the RA synovium. Given the significant role of DN2 B cells as a progenitor to pathogenic B cells in RA, it is important to conduct additional research to investigate the origins of DN2 B cells in RA and explore their potential as therapeutic targets in place of the less specific pan-B cells depletion therapies currently in use.

Circular RNA_0003800 exacerbates IL-1β-induced chondrocyte injury via miR-197-3p/SOX5 axis

BACKGROUND

Osteoarthritis (OA) is a serious degenerative disease of articular cartilage, which has a great impact on the quality of life of patients. Circular RNA (circRNA) plays an important role in OA progression. Our study aims to explore the role and mechanism of circ_0003800 in OA.

METHODS

Circ_0003800, microRNA-197-3p (miR-197-3p) and SRY-box transcription factor 5 (SOX5) contents were measured by quantitative real-time polymerase chain reaction (qRT-PCR) and western blot. Cell Counting Kit-8 (CCK8), 5-ethynyl-2'-deoxyuridine (EdU), flow cytometry, western blot and enzyme-linked immunosorbent assay (ELISA) were deployed to evaluate cell proliferation, apoptosis, extracellular matrix (ECM) degradation, inflammatory response and oxidative stress. Interaction of miR-197-3p and circ_0003800 or SOX5 was evidenced by dual-luciferase reporter system, RNA immunoprecipitation (RIP) and RNA pull down assays.

RESULTS

OA tissues and model cells had higher abundance of circ_0003800 and SOX5, while miR-197-3p content was lower. Functionally, circ_0003800 knockdown alleviated IL-1β-mediated injury in C28/I2 cells. Mechanistically, circ_0003800 could sponge miR-197-3p, and miR-197-3p could target SOX5. Besides, in-miR-197-3p reversed the suppressive effect of circ_0003800 downregulation on IL-1β-induced C28/I2 cell injury, and SOX5 overexpression could also diminish the inhibitory effect of miR-197-3p on IL-1β-induced C28/I2 cell injury.

CONCLUSION

Circ_0003800 exacerbates IL-1β-induced chondrocyte injury via miR-197-3p/SOX5 axis.

TLR7 gain-of-function genetic variation causes human lupus

Although circumstantial evidence supports enhanced Toll-like receptor 7 (TLR7) signalling as a mechanism of human systemic autoimmune disease1-7, evidence of lupus-causing TLR7 gene variants is lacking. Here we describe human systemic lupus erythematosus caused by a TLR7 gain-of-function variant. TLR7 is a sensor of viral RNA8,9 and binds to guanosine10-12. We identified a de novo, previously undescribed missense TLR7Y264H variant in a child with severe lupus and additional variants in other patients with lupus. The TLR7Y264H variant selectively increased sensing of guanosine and 2',3'-cGMP10-12, and was sufficient to cause lupus when introduced into mice. We show that enhanced TLR7 signalling drives aberrant survival of B cell receptor (BCR)-activated B cells, and in a cell-intrinsic manner, accumulation of CD11c+ age-associated B cells and germinal centre B cells. Follicular and extrafollicular helper T cells were also increased but these phenotypes were cell-extrinsic. Deficiency of MyD88 (an adaptor protein downstream of TLR7) rescued autoimmunity, aberrant B cell survival, and all cellular and serological phenotypes. Despite prominent spontaneous germinal-centre formation in Tlr7Y264H mice, autoimmunity was not ameliorated by germinal-centre deficiency, suggesting an extrafollicular origin of pathogenic B cells. We establish the importance of TLR7 and guanosine-containing self-ligands for human lupus pathogenesis, which paves the way for therapeutic TLR7 or MyD88 inhibition.

Molecular Mechanisms Driving IL-10- Producing B Cells Functions: STAT3 and c-MAF as Underestimated Central Key Regulators?

Regulatory B cells (Bregs) have been highlighted in very different pathology settings including autoimmune diseases, allergy, graft rejection, and cancer. Improving tools for the characterization of Bregs has become the main objective especially in humans. Transitional, mature B cells and plasma cells can differentiate into IL-10 producing Bregs in both mice and humans, suggesting that Bregs are not derived from unique precursors but may arise from different competent progenitors at unrestricted development stages. Moreover, in addition to IL-10 production, regulatory B cells used a broad range of suppressing mechanisms to modulate the immune response. Although Bregs have been consistently described in the literature, only a few reports described the molecular aspects that control the acquisition of the regulatory function. In this manuscript, we detailed the latest reports describing the control of IL-10, TGFβ, and GZMB production in different Breg subsets at the molecular level. We focused on the understanding of the role of the transcription factors STAT3 and c-MAF in controlling IL-10 production in murine and human B cells and how these factors may represent an important crossroad of several key drivers of the Breg response. Finally, we provided original data supporting the evidence that MAF is expressed in human IL-10- producing plasmablast and could be induced in vitro following different stimulation cocktails. At steady state, we reported that MAF is expressed in specific human B-cell tonsillar subsets including the IgD⁺ CD27⁺ unswitched population, germinal center cells and plasmablast.

Maturation trajectories and transcriptional landscape of plasmablasts and autoreactive B cells in COVID-19

In parasite and viral infections, aberrant B cell responses can suppress germinal center reactions thereby blunting long-lived memory and may provoke immunopathology including autoimmunity. Using COVID-19 as model, we set out to identify serological, cellular, and transcriptomic imprints of pathological responses linked to autoreactive B cells at single-cell resolution. We show that excessive plasmablast expansions are prognostically adverse and correlate with autoantibody production but do not hinder the formation of neutralizing antibodies. Although plasmablasts followed interleukin-4 (IL-4) and BAFF-driven developmental trajectories, were polyclonal, and not enriched in autoreactive B cells, we identified two memory populations (CD80⁺/ISG15⁺ and CD11c⁺/SOX5⁺/T-bet^+/−) with immunogenetic and transcriptional signs of autoreactivity that may be the cellular source of autoantibodies in COVID-19 and that may persist beyond recovery. Immunomodulatory interventions discouraging such adverse responses may be useful in selected patients to shift the balance from autoreactivity toward long-term memory.

•

Plasmablast expansions correlate with disease severity and autoantibodies in COVID-19

•

Patients with high plasmablast levels exhibit IGHV4-34 skewing

•

Autoreactive BCRs are enriched in atypical memory, not plasmablast populations

Uncovering transcriptional dark matter via gene annotation independent single-cell RNA sequencing analysis

Conventional scRNA-seq expression analyses rely on the availability of a high quality genome annotation. Yet, as we show here with scRNA-seq experiments and analyses spanning human, mouse, chicken, mole rat, lemur and sea urchin, genome annotations are often incomplete, in particular for organisms that are not routinely studied. To overcome this hurdle, we created a scRNA-seq analysis routine that recovers biologically relevant transcriptional activity beyond the scope of the best available genome annotation by performing scRNA-seq analysis on any region in the genome for which transcriptional products are detected. Our tool generates a single-cell expression matrix for all transcriptionally active regions (TARs), performs single-cell TAR expression analysis to identify biologically significant TARs, and then annotates TARs using gene homology analysis. This procedure uses single-cell expression analyses as a filter to direct annotation efforts to biologically significant transcripts and thereby uncovers biology to which scRNA-seq would otherwise be in the dark.

Conventional single-cell RNA sequencing analysis rely on genome annotations that may be incomplete or inaccurate especially for understudied organisms. Here the authors present a bioinformatic tool that leverages single-cell data to uncover biologically relevant transcripts beyond the best available genome annotation.

ORN: Inferring patient-specific dysregulation status of pathway modules in cancer with OR-gate Network

Pathway level understanding of cancer plays a key role in precision oncology. However, the current amount of high-throughput data cannot support the elucidation of full pathway topology. In this study, instead of directly learning the pathway network, we adapted the probabilistic OR gate to model the modular structure of pathways and regulon. The resulting model, OR-gate Network (ORN), can simultaneously infer pathway modules of somatic alterations, patient-specific pathway dysregulation status, and downstream regulon. In a trained ORN, the differentially expressed genes (DEGs) in each tumour can be explained by somatic mutations perturbing a pathway module. Furthermore, the ORN handles one of the most important properties of pathway perturbation in tumours, the mutual exclusivity. We have applied the ORN to lower-grade glioma (LGG) samples and liver hepatocellular carcinoma (LIHC) samples in TCGA and breast cancer samples from METABRIC. Both datasets have shown abnormal pathway activities related to immune response and cell cycles. In LGG samples, ORN identified pathway modules closely related to glioma development and revealed two pathways closely related to patient survival. We had similar results with LIHC samples. Additional results from the METABRIC datasets showed that ORN could characterize critical mechanisms of cancer and connect them to less studied somatic mutations (e.g., BAP1, MIR604, MICAL3, and telomere activities), which may generate novel hypothesis for targeted therapy.

A genome-wide association study for harness racing success in the Norwegian-Swedish coldblooded trotter reveals genes for learning and energy metabolism

Background

Although harness racing is of high economic importance to the global equine industry, significant genomic resources have yet to be applied to mapping harness racing success. To identify genomic regions associated with harness racing success, the current study performs genome-wide association analyses with three racing performance traits in the Norwegian-Swedish Coldblooded Trotter using the 670 K Axiom Equine Genotyping Array.

Results

Following quality control, 613 horses and 359,635 SNPs were retained for further analysis. After strict Bonferroni correction, nine genome-wide significant SNPs were identified for career earnings. No genome-wide significant SNPs were identified for number of gallops or best km time. However, four suggestive genome-wide significant SNPs were identified for number of gallops, while 19 were identified for best km time. Multiple genes related to intelligence, energy metabolism, and immune function were identified as potential candidate genes for harness racing success.

Conclusions

Apart from the physiological requirements needed for a harness racing horse to be successful, the results of the current study also advocate learning ability and memory as important elements for harness racing success. Further exploration into the mental capacity required for a horse to achieve racing success is likely warranted.

Electronic supplementary material

The online version of this article (10.1186/s12863-018-0670-3) contains supplementary material, which is available to authorized users.

Transcription Factor SOX5 Promotes the Migration and Invasion of Fibroblast-Like Synoviocytes in Part by Regulating MMP-9 Expression in Collagen-Induced Arthritis

Objectives

Fibroblast-like synoviocytes (FLS) exhibit a unique aggressive phenotype in rheumatoid arthritis (RA). Increased FLS migration and subsequent invasion of the extracellular matrix are essential to joint destruction in RA. Our previous research reported that transcription factor SOX5 was highly expressed in RA-FLS. Here, the effects of SOX5 in RA-FLS migration and invasion will be investigated.

Methods

The migration and invasion of RA-FLS were evaluated using a transwell chamber assay. The expression of several potential SOX5-targeted genes, including matrix metalloproteinases (MMP-1, 2, 3 and 9), chemokines (CCL4, CCL2, CCR5 and CCR2), and pro-inflammatory cytokines (TNF-α and IL-6), were examined in RA-FLS using SOX5 gain- and loss-of-function study. The molecular mechanisms of SOX5-mediated MMP-9 expressions were assayed by luciferase reporter gene and chromatin immunoprecipitation (ChIP) studies. The in vivo effect of SOX5 on FLS migration and invasion was examined using collagen-induced arthritis (CIA) in DBA/1J mice.

Results

Knockdown SOX5 decreased lamellipodium formation, migration, and invasion of RA-FLS. The expression of MMP-9 was the only gene tested to be concomitantly affected by silencing or overexpressing SOX5. ChIP assay revealed that SOX5 was bound to the MMP-9 promoter in RA-FLS. The overexpression of SOX5 markedly enhanced the MMP-9 promoter activity, and specific deletion of a putative SOX5-binding site in MMP-9 promoter diminished this promoter-driven transcription in FLS. Locally knocked down SOX5 inhibited MMP-9 expression in the joint tissue and reduced pannus migration and invasion into the cartilage in CIA mice.

Conclusion

SOX5 plays a novel role in mediating migration and invasion of FLS in part by regulating MMP-9 expression in RA.

High SYK Expression Drives Constitutive Activation of CD21low B Cells

Human CD21^low B cells present with an activated phenotype and accumulate in distinct disorders connected with chronic immune stimulation. Signaling studies had revealed an increased basal phosphorylation of spleen tyrosine kinase (SYK) and phospholipase Cγ2. Additional BCR stimulation of these constitutively active cells, however, led to reduced activation of these signaling molecules and subsequently NF-κB and Ca²⁺ activation. In this article, we demonstrate that high SYK expression is a common feature of CD21^low B cells independent of the underlying disorder, and that this high expression is sufficient to drive constitutive phosphorylation of SYK and its immediate targets Bruton's tyrosine kinase and phospholipase Cγ2. Inhibition of SYK activity eliminated features of the constitutive activation in these cells and partly restored BCR signaling. High SYK expression is especially induced by CpG or CD40L in combination with IL-21, but not BCR stimulation, suggesting the importance of the immune-stimulatory context for the induction of this B cell phenotype. In summary, high SYK expression is a common feature of human CD21^low B cells and presumably results from chronic activation in inflammatory environments present in a subgroup of patients with heterogeneous disorders like chronic infection, autoimmunity, and immunodeficiency. High SYK expression by itself drives the constitutive activation observed in these B cells, which in turn may contribute to the hyporesponsiveness upon BCR stimulation. Given the high prevalence of autoreactive clones among CD21^low B cells in autoimmune disorders, the dominant role of SYK in CD21^low B cells may provide a new option for therapeutic interventions in patients with expanded CD21^low B cells and humoral autoimmunity.

Memory B cells: total recall

Immunological memory is a cornerstone of adaptive immune responses in higher vertebrates. The remarkable ability to generate memory cells following Ag exposure, in the context of natural infection or immunization, provides long-lived protection against infectious diseases, often for the hosts' lifetime. Indeed, the generation of memory B cells and long-lived plasma cells underpins the success of most vaccines. The concept of immunological memory is not new-it was first proposed nearly 2500 years ago. While our understanding of the complexities of humoral and cell-mediated memory continues to evolve, important aspects of this process remain unresolved. Here, we will provide an overview of recent advances in B-cell memory in mice and humans, and in health and disease.

Modulation of IL-6 induced RANKL expression in arthritic synovium by a transcription factor SOX5

Receptor activator of nuclear factor κB ligand (RANKL) is critically involved in bone erosion of rheumatoid arthritis (RA). We previously reported association between younger age at onset of RA and a RANKL promoter SNP that conferred an elevated promoter activity via binding to a transcription factor SOX5. Here we study the regulation of SOX5 levels in relation to RANKL expression in RA synovial fibroblasts (SF) and the development of bone erosion in the collagen-induced arthritis (CIA) mouse. Our data indicated SOX5 levels were higher in synovium and synovial fluid from RA compared to osteoarthritis patients. Pro-inflammatory cytokines upregulated SOX5 and RANKL expression in both primary RA SF and the rheumatoid synovial fibroblast cell line, MH7A. Overexpression of SOX5 resulted in significantly increased RANKL levels, while knockdown of SOX5 resulted in diminished IL-6 mediated RANKL upregulation in MH7A cells. Chromatin immunoprecipitation (ChIP) showed approximately 3-fold enrichment of RANKL-specific DNA in anti-SOX5 immunoprecipitate in IL-6 treated MH7A cells as compared to untreated cells. Locally silencing SOX5 gene significantly diminished RANKL positive cells and bone erosion in CIA mice. These findings suggest SOX5 is an important regulator of IL-6-induced RANKL expression in RA SF.