Author Correction: Full-length mRNA-Seq from single-cell levels of RNA and individual circulating tumor cells.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Intussusceptive Vascular Remodeling Precedes Pathological Neovascularization.

Objective- Pathological neovascularization is crucial for progression and morbidity of serious diseases such as cancer, diabetic retinopathy, and age-related macular degeneration. While mechanisms of ongoing pathological neovascularization have been extensively studied, the initiating pathological vascular remodeling (PVR) events, which precede neovascularization remains poorly understood. Here, we identify novel molecular and cellular mechanisms of preneovascular PVR, by using the adult choriocapillaris as a model. Approach and Results- Using hypoxia or forced overexpression of VEGF (vascular endothelial growth factor) in the subretinal space to induce PVR in zebrafish and rats respectively, and by analyzing choriocapillaris membranes adjacent to choroidal neovascular lesions from age-related macular degeneration patients, we show that the choriocapillaris undergo robust induction of vascular intussusception and permeability at preneovascular stages of PVR. This PVR response included endothelial cell proliferation, formation of endothelial luminal processes, extensive vesiculation and thickening of the endothelium, degradation of collagen fibers, and splitting of existing extravascular columns. RNA-sequencing established a role for endothelial tight junction disruption, cytoskeletal remodeling, vesicle- and cilium biogenesis in this process. Mechanistically, using genetic gain- and loss-of-function zebrafish models and analysis of primary human choriocapillaris endothelial cells, we determined that HIF (hypoxia-induced factor)-1α-VEGF-A-VEGFR2 signaling was important for hypoxia-induced PVR. Conclusions- Our findings reveal that PVR involving intussusception and splitting of extravascular columns, endothelial proliferation, vesiculation, fenestration, and thickening is induced before neovascularization, suggesting that identifying and targeting these processes may prevent development of advanced neovascular disease in the future. Visual Overview- An online visual overview is available for this article.

B cell alterations during BAFF inhibition with belimumab in SLE.

BACKGROUND: Systemic lupus erythematosus (SLE) is a systemic autoimmune disease, which exhibits multiple B cell abnormalities including expanded populations of memory B cells and elevated levels of autoantibodies. Belimumab is a monoclonal antibody targeting the B cell cytokine BAFF (a.k.a. BLyS), approved for the treatment of SLE. METHODS: In this prospective cohort study, B cells from peripheral blood of 23 SLE patients initiating belimumab treatment and followed longitudinally for up to three years, were assessed using mass cytometry. FINDINGS: B cells decreased during the study period, with a rapid decrease of both naïve and CD11c+CD21- B cells at the first follow-up visit, followed by a continuous reduction at subsequent follow-ups. In contrast, plasma cells and switched memory B cells remained stable throughout the study. The observed immunological changes correlated with early, but not late, clinical improvements. Moreover, high baseline B cell counts were predictive of failure to attain low disease activity. In summary, our data unveiled both rapid and gradual later therapy-associated alterations of both known and unforeseen B cell phenotypes. INTERPRETATION: Our results suggest that evaluation of B cell counts might prove useful prior to initiation of belimumab treatment and that early treatment evaluation and discontinuation might underestimate delayed clinical improvements resultant of late B cell changes.

Effect of CTLA4-Ig (abatacept) treatment on T cells and B cells in peripheral blood of patients with polymyositis and dermatomyositis.

We aimed to evaluate in vivo effects of abatacept on phenotypes of T and B cells in the circulation of myositis patients in a sub-study of the ARTEMIS trial. Twelve patients with paired frozen PBMCs before and after 6-month abatacept treatment were included in this sub-study where mass cytometry (CyTOF) was chosen as a technology to be tested for its utility in a real-life clinical immune monitoring setting. Using CyTOF, the peripheral T cell phenotypes demonstrated considerable variation over time and between individuals precluding the identification of treatment-specific changes. We therefore conclude that studies of patient cohorts displaying wide clinical heterogeneity using mass cytometry must be relatively large in order to be suited for discovery research and immune monitoring. Still, we did find some correlations with functional muscle outcome, namely positive correlations between the ratio of CD4+ T cells and CD8+ T cells (CD4/CD8) in peripheral blood samples both at baseline and after treatment with muscle endurance improvement as assessed by the functional index-2 (FI-2) test. Our data suggest that the CD4/CD8 ratio in circulation at time of active disease may be a predictor of treatment efficacy in myositis patients.

EOMES-positive CD4+ T cells are increased in PTPN22 (1858T) risk allele carriers.

The presence of the PTPN22 risk allele (1858T) is associated with several autoimmune diseases including rheumatoid arthritis (RA). Despite a number of studies exploring the function of PTPN22 in T cells, the exact impact of the PTPN22 risk allele on T-cell function in humans is still unclear. In this study, using RNA sequencing, we show that, upon TCR-activation, naïve human CD4+ T cells homozygous for the PTPN22 risk allele overexpress a set of genes including CFLAR and 4-1BB, which are important for cytotoxic T-cell differentiation. Moreover, the protein expression of the T-box transcription factor Eomesodermin (EOMES) was increased in T cells from healthy donors homozygous for the PTPN22 risk allele and correlated with a decreased number of naïve CD4+ T cells. There was no difference in the frequency of other CD4+ T-cell subsets (Th1, Th17, Tfh, Treg). Finally, an accumulation of EOMES+ CD4+ T cells was observed in synovial fluid of RA patients with a more pronounced production of Perforin-1 in PTPN22 risk allele carriers. Altogether, we propose a novel mechanism of action of PTPN22 risk allele through the generation of cytotoxic CD4+ T cells and identify EOMES+ CD4+ T cells as a relevant T-cell subset in RA pathogenesis.

Mouse Model of Alagille Syndrome and Mechanisms of Jagged1 Missense Mutations.

BACKGROUND & AIMS: Alagille syndrome is a genetic disorder characterized by cholestasis, ocular abnormalities, characteristic facial features, heart defects, and vertebral malformations. Most cases are associated with mutations in JAGGED1 (JAG1), which encodes a Notch ligand, although it is not clear how these contribute to disease development. We aimed to develop a mouse model of Alagille syndrome to elucidate these mechanisms. METHODS: Mice with a missense mutation (H268Q) in Jag1 (Jag1+/Ndr mice) were outbred to a C3H/C57bl6 background to generate a mouse model for Alagille syndrome (Jag1Ndr/Ndr mice). Liver tissues were collected at different timepoints during development, analyzed by histology, and liver organoids were cultured and analyzed. We performed transcriptome analysis of Jag1Ndr/Ndr livers and livers from patients with Alagille syndrome, cross-referenced to the Human Protein Atlas, to identify commonly dysregulated pathways and biliary markers. We used species-specific transcriptome separation and ligand-receptor interaction assays to measure Notch signaling and the ability of JAG1Ndr to bind or activate Notch receptors. We studied signaling of JAG1 and JAG1Ndr via NOTCH 1, NOTCH2, and NOTCH3 and resulting gene expression patterns in parental and NOTCH1-expressing C2C12 cell lines. RESULTS: Jag1Ndr/Ndr mice had many features of Alagille syndrome, including eye, heart, and liver defects. Bile duct differentiation, morphogenesis, and function were dysregulated in newborn Jag1Ndr/Ndr mice, with aberrations in cholangiocyte polarity, but these defects improved in adult mice. Jag1Ndr/Ndr liver organoids collapsed in culture, indicating structural instability. Whole-transcriptome sequence analyses of liver tissues from mice and patients with Alagille syndrome identified dysregulated genes encoding proteins enriched at the apical side of cholangiocytes, including CFTR and SLC5A1, as well as reduced expression of IGF1. Exposure of Notch-expressing cells to JAG1Ndr, compared with JAG1, led to hypomorphic Notch signaling, based on transcriptome analysis. JAG1-expressing cells, but not JAG1Ndr-expressing cells, bound soluble Notch1 extracellular domain, quantified by flow cytometry. However, JAG1 and JAG1Ndr cells each bound NOTCH2, and signaling from NOTCH2 signaling was reduced but not completely inhibited, in response to JAG1Ndr compared with JAG1. CONCLUSIONS: In mice, expression of a missense mutant of Jag1 (Jag1Ndr) disrupts bile duct development and recapitulates Alagille syndrome phenotypes in heart, eye, and craniofacial dysmorphology. JAG1Ndr does not bind NOTCH1, but binds NOTCH2, and elicits hypomorphic signaling. This mouse model can be used to study other features of Alagille syndrome and organ development.

Smoking and pre-existing organ damage reduce the efficacy of belimumab in systemic lupus erythematosus.

OBJECTIVES: Belimumab is the first biologic drug approved for Systemic Lupus Erythematosus (SLE). Here, we aimed to investigate the effects of belimumab on clinical and serologic outcomes, and sought to identify predictors of treatment response in three Swedish real-life settings. METHODS: Fifty-eight patients were enrolled at initiation of belimumab and followed longitudinally for up to 53months. Surveillance outcomes included the SLE Disease Activity Index 2000 (SLEDAI-2K), 100mm Visual Analogue Scales for Physician's Global Assessment (PGA), fatigue, pain and general health, and the Systemic Lupus International Collaborating Clinics/American College of Rheumatology Damage Index (SDI). Assessment of treatment response included the SLE responder index (SRI). B lymphocyte stimulator (BLyS) levels were determined using ELISA. RESULTS: SLEDAI-2K (median baseline score: 8.0; IQR: 4.0-13.8), PGA and corticosteroid use decreased during therapy, and patients reported improvements on fatigue, pain, and general health (p<0.0001 for all). SDI scores remained stable (p=0.08). Patients with baseline SDI scores >1 showed decreased probability and prolonged time to attain SRI response (HR: 0.449; 95% CI: 0.208-0.967), as did current smokers compared with non-smokers (HR: 0.103; 95% CI: 0.025-0.427). In contrast, baseline BLyS levels ≥1.2ng/mL predicted increased probability and shorter time to attain SRI response (HR: 2.566; 95% CI: 1.222-5.387). CONCLUSIONS: Disease activity and corticosteroid usage decreased, patient-reported outcomes improved, and no significant organ damage was accrued during follow-up. Smoking and organ damage predicted reduced treatment efficacy. These findings might contribute to a better selection of patients who are likely to benefit from belimumab.

Loss of CSL Unlocks a Hypoxic Response and Enhanced Tumor Growth Potential in Breast Cancer Cells.

Notch signaling is an important regulator of stem cell differentiation. All canonical Notch signaling is transmitted through the DNA-binding protein CSL, and hyperactivated Notch signaling is associated with tumor development; thus it may be anticipated that CSL deficiency should reduce tumor growth. In contrast, we report that genetic removal of CSL in breast tumor cells caused accelerated growth of xenografted tumors. Loss of CSL unleashed a hypoxic response during normoxic conditions, manifested by stabilization of the HIF1α protein and acquisition of a polyploid giant-cell, cancer stem cell-like, phenotype. At the transcriptome level, loss of CSL upregulated more than 1,750 genes and less than 3% of those genes were part of the Notch transcriptional signature. Collectively, this suggests that CSL exerts functions beyond serving as the central node in the Notch signaling cascade and reveals a role for CSL in tumorigenesis and regulation of the cellular hypoxic response.

Decoding breast cancer tissue-stroma interactions using species-specific sequencing.

INTRODUCTION: Decoding transcriptional effects of experimental tissue-tissue or cell-cell interactions is important; for example, to better understand tumor-stroma interactions after transplantation of human cells into mouse (xenografting). Transcriptome analysis of intermixed human and mouse cells has, however, frequently relied on the need to separate the two cell populations prior to transcriptome analysis, which introduces confounding effects on gene expression. METHODS: To circumvent this problem, we here describe a bioinformatics-based, genome-wide transcriptome analysis technique, which allows the human and mouse transcriptomes to be decoded from a mixed mouse and human cell population. The technique is based on a bioinformatic separation of the mouse and human transcriptomes from the initial mixed-species transcriptome resulting from sequencing an excised tumor/stroma specimen without prior cell sorting. RESULTS: Under stringent separation criteria, i.e., with a read misassignment frequency of 0.2 %, we show that 99 % of the genes can successfully be assigned to be of mouse or human origin, both in silico, in cultured cells and in vivo. We use a new species-specific sequencing technology-referred to as S(3) ("S-cube")-to provide new insights into the Notch downstream response following Notch ligand-stimulation and to explore transcriptional changes following transplantation of two different breast cancer cell lines (luminal MCF7 and basal-type MDA-MB-231) into mammary fat pad tissue in mice of different immunological status. We find that MCF7 and MDA-MB-231 respond differently to fat pad xenografting and the stromal response to transplantation of MCF7 and MDA-MB-231 cells was also distinct. CONCLUSIONS: In conclusion, the data show that the S(3) technology allows for faithful recording of transcriptomic changes when human and mouse cells are intermixed and that it can be applied to address a broad spectrum of research questions.

RNA helicase A is a downstream mediator of KIF1Bß tumor-suppressor function in neuroblastoma.

UNLABELLED: Inherited KIF1B loss-of-function mutations in neuroblastomas and pheochromocytomas implicate the kinesin KIF1B as a 1p36.2 tumor suppressor. However, the mechanism of tumor suppression is unknown. We found that KIF1B isoform ß (KIF1Bß) interacts with RNA helicase A (DHX9), causing nuclear accumulation of DHX9, followed by subsequent induction of the proapoptotic XIAP-associated factor 1 (XAF1) and, consequently, apoptosis. Pheochromocytoma and neuroblastoma arise from neural crest progenitors that compete for growth factors such as nerve growth factor (NGF) during development. KIF1Bß is required for developmental apoptosis induced by competition for NGF. We show that DHX9 is induced by and required for apoptosis stimulated by NGF deprivation. Moreover, neuroblastomas with chromosomal deletion of 1p36 exhibit loss of KIF1Bß expression and impaired DHX9 nuclear localization, implicating the loss of DHX9 nuclear activity in neuroblastoma pathogenesis. SIGNIFICANCE: KIF1Bß has neuroblastoma tumor-suppressor properties and promotes and requires nuclear-localized DHX9 for its apoptotic function by activating XAF1 expression. Loss of KIF1Bß alters subcellular localization of DHX9 and diminishes NGF dependence of sympathetic neurons, leading to reduced culling of neural progenitors, and, therefore, might predispose to tumor formation.

Full-length mRNA-Seq from single-cell levels of RNA and individual circulating tumor cells.

Genome-wide transcriptome analyses are routinely used to monitor tissue-, disease- and cell type­specific gene expression, but it has been technically challenging to generate expression profiles from single cells. Here we describe a robust mRNA-Seq protocol (Smart-Seq) that is applicable down to single cell levels. Compared with existing methods, Smart-Seq has improved read coverage across transcripts, which enhances detailed analyses of alternative transcript isoforms and identification of single-nucleotide polymorphisms. We determined the sensitivity and quantitative accuracy of Smart-Seq for single-cell transcriptomics by evaluating it on total RNA dilution series. We found that although gene expression estimates from single cells have increased noise, hundreds of differentially expressed genes could be identified using few cells per cell type. Applying Smart-Seq to circulating tumor cells from melanomas, we identified distinct gene expression patterns, including candidate biomarkers for melanoma circulating tumor cells. Our protocol will be useful for addressing fundamental biological problems requiring genome-wide transcriptome profiling in rare cells.