Single-cell RNA-sequencing of virus-specific cellular immune responses in chronic hepatitis B patients.

Chronic hepatitis B (CHB) is a major global health challenge. CHB can be controlled by antivirals but a therapeutic cure is lacking. CHB is characterized by limited HBV-specific T cell reactivity and functionality and expression of inhibitory receptors. The mechanisms driving these T cell phenotypes are only partially understood. Here, we created a single-cell RNA-sequencing dataset of HBV immune responses in patients to contribute to a better understanding of the dysregulated immunity. Blood samples of a well-defined cohort of 21 CHB and 10 healthy controls, including a subset of 5 matched liver biopsies, were collected. scRNA-seq data of total immune cells (55,825) plus sorted HBV-specific (1,963), non-naive (32,773) and PD1+ T cells (96,631) was generated using the 10X Genomics platform (186,123 cells) or the full-length Smart-seq2 protocol (1,069 cells). The shared transcript count matrices of single-cells serve as a valuable resource describing transcriptional changes underlying dysfunctional HBV-related T cell responses in blood and liver tissue and offers the opportunity to identify targets or biomarkers for HBV-related immune exhaustion.

Transcriptional profiling upon T cell stimulation reveals down-regulation of inflammatory pathways in T and B cells in SLE versus Sjögren's syndrome.

Systemic lupus erythematosus (SLE) and primary Sjögren's syndrome (pSS) share clinical as well as pathogenic similarities. Although previous studies suggest various abnormalities in different immune cell compartments, dedicated cell-type specific transcriptomic signatures are often masked by patient heterogeneity. Here, we performed transcriptional profiling of isolated CD4, CD8, CD16 and CD19 lymphocytes from pSS and SLE patients upon T cell stimulation, in addition to a steady-state condition directly after blood drawing, in total comprising 581 sequencing samples. T cell stimulation, which induced a pronounced inflammatory response in all four cell types, gave rise to substantial re-modulation of lymphocyte subsets in the two autoimmune diseases compared to healthy controls, far exceeding the transcriptomic differences detected at steady-state. In particular, we detected cell-type and disease-specific down-regulation of a range of pro-inflammatory cytokine and chemokine pathways. Such differences between SLE and pSS patients are instrumental for selective immune targeting by future therapies.

Deciphering molecular and cellular ex vivo responses to bispecific antibodies PD1-TIM3 and PD1-LAG3 in human tumors.

BACKGROUND: Next-generation cancer immunotherapies are designed to broaden the therapeutic repertoire by targeting new immune checkpoints including lymphocyte-activation gene 3 (LAG-3) and T cell immunoglobulin and mucin-domain containing-3 (TIM-3). Yet, the molecular and cellular mechanisms by which either receptor functions to mediate its inhibitory effects are still poorly understood. Similarly, little is known on the differential effects of dual, compared with single, checkpoint inhibition. METHODS: We here performed in-depth characterization, including multicolor flow cytometry, single cell RNA sequencing and multiplex supernatant analysis, using tumor single cell suspensions from patients with cancer treated ex vivo with novel bispecific antibodies targeting programmed cell death protein 1 (PD-1) and TIM-3 (PD1-TIM3), PD-1 and LAG-3 (PD1-LAG3), or with anti-PD-1. RESULTS: We identified patient samples which were responsive to PD1-TIM3, PD1-LAG3 or anti-PD-1 using an in vitro approach, validated by the analysis of 659 soluble proteins and enrichment for an anti-PD-1 responder signature. We found increased abundance of an activated (HLA-DR+CD25+GranzymeB+) CD8+ T cell subset and of proliferating CD8+ T cells, in response to bispecific antibody or anti-PD-1 treatment. Bispecific antibodies, but not anti-PD-1, significantly increased the abundance of a proliferating natural killer cell subset, which exhibited enrichment for a tissue-residency signature. Key phenotypic and transcriptional changes occurred in a PD-1+CXCL13+CD4+ T cell subset, in response to all treatments, including increased interleukin-17 secretion and signaling toward plasma cells. Interestingly, LAG-3 protein upregulation was detected as a unique pharmacodynamic effect mediated by PD1-LAG3, but not by PD1-TIM3 or anti-PD-1. CONCLUSIONS: Our in vitro system reliably assessed responses to bispecific antibodies co-targeting PD-1 together with LAG-3 or TIM-3 using patients' tumor infiltrating immune cells and revealed transcriptional and phenotypic imprinting by bispecific antibody formats currently tested in early clinical trials.

Besca, a single-cell transcriptomics analysis toolkit to accelerate translational research.

Single-cell RNA sequencing (scRNA-seq) revolutionized our understanding of disease biology. The promise it presents to also transform translational research requires highly standardized and robust software workflows. Here, we present the toolkit Besca, which streamlines scRNA-seq analyses and their use to deconvolute bulk RNA-seq data according to current best practices. Beyond a standard workflow covering quality control, filtering, and clustering, two complementary Besca modules, utilizing hierarchical cell signatures and supervised machine learning, automate cell annotation and provide harmonized nomenclatures. Subsequently, the gene expression profiles can be employed to estimate cell type proportions in bulk transcriptomics data. Using multiple, diverse scRNA-seq datasets, some stemming from highly heterogeneous tumor tissue, we show how Besca aids acceleration, interoperability, reusability and interpretability of scRNA-seq data analyses, meeting crucial demands in translational research and beyond.

Critical assessment of coiled-coil predictions based on protein structure data.

Coiled-coil regions were among the first protein motifs described structurally and theoretically. The simplicity of the motif promises that coiled-coil regions can be detected with reasonable accuracy and precision in any protein sequence. Here, we re-evaluated the most commonly used coiled-coil prediction tools with respect to the most comprehensive reference data set available, the entire Protein Data Bank, down to each amino acid and its secondary structure. Apart from the 30-fold difference in minimum and maximum number of coiled coils predicted the tools strongly vary in where they predict coiled-coil regions. Accordingly, there is a high number of false predictions and missed, true coiled-coil regions. The evaluation of the binary classification metrics in comparison with naïve coin-flip models and the calculation of the Matthews correlation coefficient, the most reliable performance metric for imbalanced data sets, suggests that the tested tools' performance is close to random. This implicates that the tools' predictions have only limited informative value. Coiled-coil predictions are often used to interpret biochemical data and are part of in-silico functional genome annotation. Our results indicate that these predictions should be treated very cautiously and need to be supported and validated by experimental evidence.

Deep Phenotyping of CD11c+ B Cells in Systemic Autoimmunity and Controls.

Circulating CD11c+ B cells are a key phenomenon in certain types of autoimmunity but have also been described in the context of regular immune responses (i.e., infections, vaccination). Using mass cytometry to profile 46 different markers on individual immune cells, we systematically initially confirmed the presence of increased CD11c+ B cells in the blood of systemic lupus erythematosus (SLE) patients. Notably, significant differences in the expression of CD21, CD27, and CD38 became apparent between CD11c- and CD11c+ B cells. We observed direct correlation of the frequency of CD21-CD27- B cells and CD21-CD38- B cells with CD11c+ B cells, which were most pronounced in SLE compared to primary Sjögren's syndrome patients (pSS) and healthy donors (HD). Thus, CD11c+ B cells resided mainly within memory subsets and were enriched in CD27-IgD-, CD21-CD27-, and CD21-CD38- B cell phenotypes. CD11c+ B cells from all donor groups (SLE, pSS, and HD) showed enhanced CD69, Ki-67, CD45RO, CD45RA, and CD19 expression, whereas the membrane expression of CXCR5 and CD21 were diminished. Notably, SLE CD11c+ B cells showed enhanced expression of the checkpoint molecules CD86, PD1, PDL1, CD137, VISTA, and CTLA-4 compared to HD. The substantial increase of CD11c+ B cells with a CD21- phenotype co-expressing distinct activation and checkpoint markers, points to a quantitative increased alternate (extrafollicular) B cell activation route possibly related to abnormal immune regulation as seen under the striking inflammatory conditions of SLE which shows a characteristic PD-1/PD-L1 upregulation.

Profiling of naïve and primed human pluripotent stem cells reveals state-associated miRNAs.

Naïve human pluripotent stem cells (hPSC) resemble the embryonic epiblast at an earlier time-point in development than conventional, 'primed' hPSC. We present a comprehensive miRNA profiling of naïve-to-primed transition in hPSC, a process recapitulating aspects of early in vivo embryogenesis. We identify miR-143-3p and miR-22-3p as markers of the naïve state and miR-363-5p, several members of the miR-17 family, miR-302 family as primed markers. We uncover that miR-371-373 are highly expressed in naïve hPSC. MiR-371-373 are the human homologs of the mouse miR-290 family, which are the most highly expressed miRNAs in naïve mouse PSC. This aligns with the consensus that naïve hPSC resemble mouse naive PSC, showing that the absence of miR-371-373 in conventional hPSC is due to cell state rather than a species difference.

Large-Scale Production of Human iPSC-Derived Macrophages for Drug Screening.

Tissue-resident macrophages are key players in inflammatory processes, and their activation and functionality are crucial in health and disease. Numerous diseases are associated with alterations in homeostasis or dysregulation of the innate immune system, including allergic reactions, autoimmune diseases, and cancer. Macrophages are a prime target for drug discovery due to their major regulatory role in health and disease. Currently, the main sources of macrophages used for therapeutic compound screening are primary cells isolated from blood or tissue or immortalized or neoplastic cell lines (e.g., THP-1). Here, we describe an improved method to employ induced pluripotent stem cells (iPSCs) for the high-yield, large-scale production of cells resembling tissue-resident macrophages. For this, iPSC-derived macrophage-like cells are thoroughly characterized to confirm their cell identity and thus their suitability for drug screening purposes. These iPSC-derived macrophages show strong cellular identity with primary macrophages and recapitulate key functional characteristics, including cytokine release, phagocytosis, and chemotaxis. Furthermore, we demonstrate that genetic modifications can be readily introduced at the macrophage-like progenitor stage in order to interrogate drug target-relevant pathways. In summary, this novel method overcomes previous shortcomings with primary and leukemic cells and facilitates large-scale production of genetically modified iPSC-derived macrophages for drug screening applications.

Rates of homology directed repair of CRISPR-Cas9 induced double strand breaks are lower in naïve compared to primed human pluripotent stem cells.

Gene editing in human pluripotent stem cells (hPSC) is a powerful tool for understanding biology, for drug discovery and gene therapy. Naïve hPSC have been suggested to be superior for gene editing compared to conventional 'primed' hPSC. Using droplet digital PCR, we uncover the kinetics of Cas9-induced double strand break repair in conventional hPSC. Cut but unrepaired alleles reach their maximum after 12-24 h. Homology directed repair plateaus after 24 h, whereas repair by non-homologous end joining continues until 48 h after Cas9 introduction. Using this method, we demonstrate that the rate of homology directed repair to resolve Cas9-induced double strand breaks is 40% lower in naïve hPSC compared to conventional hPSC, correlating with, and feasibly explained by, a higher number of cells in G1 phase of the cell cycle in naïve hPSC. Therefore, naïve hPSC are less efficient for CRISPR/Cas9-mediated homology directed repair.

The Protein-Coding Human Genome: Annotating High-Hanging Fruits.

The major transcript variants of human protein-coding genes are annotated to a certain degree of accuracy combining manual curation, transcript data, and proteomics evidence. However, there is considerable disagreement on the annotation of about 2000 genes-they can be protein-coding, noncoding, or pseudogenes-and on the annotation of most of the predicted alternative transcripts. Pure transcriptome mapping approaches seem to be limited in discriminating functional expression from noise. These limitations have partially been overcome by dedicated algorithms to detect alternative spliced micro-exons and wobble splice variants. Recently, knowledge about splice mechanism and protein structure are incorporated into an algorithm to predict neighboring homologous exons, often spliced in a mutually exclusive manner. Predicted exons are evaluated by transcript data, structural compatibility, and evolutionary conservation, revealing hundreds of novel coding exons and splice mechanism re-assignments. The emerging human pan-genome is necessitating distinctive annotations incorporating differences between individuals and between populations.

Correction to: Detect tissue heterogeneity in gene expression data with BioQC.

After the publication of this work [1], a mistake was noticed in the Eq. 1. Given an m × n expression matrix with m genes and samples of n tissues, the correct definition of the Gini index for gene i is.

Targeting RNA structure in SMN2 reverses spinal muscular atrophy molecular phenotypes.

Modification of SMN2 exon 7 (E7) splicing is a validated therapeutic strategy against spinal muscular atrophy (SMA). However, a target-based approach to identify small-molecule E7 splicing modifiers has not been attempted, which could reveal novel therapies with improved mechanistic insight. Here, we chose as a target the stem-loop RNA structure TSL2, which overlaps with the 5' splicing site of E7. A small-molecule TSL2-binding compound, homocarbonyltopsentin (PK4C9), was identified that increases E7 splicing to therapeutic levels and rescues downstream molecular alterations in SMA cells. High-resolution NMR combined with molecular modelling revealed that PK4C9 binds to pentaloop conformations of TSL2 and promotes a shift to triloop conformations that display enhanced E7 splicing. Collectively, our study validates TSL2 as a target for small-molecule drug discovery in SMA, identifies a novel mechanism of action for an E7 splicing modifier, and sets a precedent for other splicing-mediated diseases where RNA structure could be similarly targeted.

Tissue-Specific Immunoregulation: A Call for Better Understanding of the "Immunostat" in the Context of Cancer.

Checkpoint inhibitor therapy has been a breakthrough in cancer research, but only some patients with cancer derive substantial benefit. Although mechanisms underlying sensitivity and resistance to checkpoint inhibitors are being elucidated, the importance of organ-specific regulation of immunity is currently underappreciated. Here, we call for a greater understanding of tissue-specific immunoregulation, namely, "tissue-specific immunostats," to make advances in treatments for cancer. A better understanding of how individual organs at baseline regulate the immune system could enable an improved precision medicine approach to cancer immunotherapy. Cancer Discov; 8(4); 395-402. ©2018 AACR.

The landscape of human mutually exclusive splicing.

Mutually exclusive splicing of exons is a mechanism of functional gene and protein diversification with pivotal roles in organismal development and diseases such as Timothy syndrome, cardiomyopathy and cancer in humans. In order to obtain a first genomewide estimate of the extent and biological role of mutually exclusive splicing in humans, we predicted and subsequently validated mutually exclusive exons (MXEs) using 515 publically available RNA-Seq datasets. Here, we provide evidence for the expression of over 855 MXEs, 42% of which represent novel exons, increasing the annotated human mutually exclusive exome more than fivefold. The data provide strong evidence for the existence of large and multi-cluster MXEs in higher vertebrates and offer new insights into MXE evolution. More than 82% of the MXE clusters are conserved in mammals, and five clusters have homologous clusters in Drosophila Finally, MXEs are significantly enriched in pathogenic mutations and their spatio-temporal expression might predict human disease pathology.

Correction: Distribution and evolution of stable single α-helices (SAH domains) in myosin motor proteins.

[This corrects the article DOI: 10.1371/journal.pone.0174639.].

Detect tissue heterogeneity in gene expression data with BioQC.

BACKGROUND: Gene expression data can be compromised by cells originating from other tissues than the target tissue of profiling. Failures in detecting such tissue heterogeneity have profound implications on data interpretation and reproducibility. A computational tool explicitly addressing the issue is warranted. RESULTS: We introduce BioQC, a R/Bioconductor software package to detect tissue heterogeneity in gene expression data. To this end BioQC implements a computationally efficient Wilcoxon-Mann-Whitney test and provides more than 150 signatures of tissue-enriched genes derived from large-scale transcriptomics studies. Simulation experiments show that BioQC is both fast and sensitive in detecting tissue heterogeneity. In a case study with whole-organ profiling data, BioQC predicted contamination events that are confirmed by quantitative RT-PCR. Applied to transcriptomics data of the Genotype-Tissue Expression (GTEx) project, BioQC reveals clustering of samples and suggests that some samples likely suffer from tissue heterogeneity. CONCLUSIONS: Our experience with gene expression data indicates a prevalence of tissue heterogeneity that often goes unnoticed. BioQC addresses the issue by integrating prior knowledge with a scalable algorithm. We propose BioQC as a first-line tool to ensure quality and reproducibility of gene expression data.

Distribution and evolution of stable single α-helices (SAH domains) in myosin motor proteins.

Stable single-alpha helices (SAHs) are versatile structural elements in many prokaryotic and eukaryotic proteins acting as semi-flexible linkers and constant force springs. This way SAH-domains function as part of the lever of many different myosins. Canonical myosin levers consist of one or several IQ-motifs to which light chains such as calmodulin bind. SAH-domains provide flexibility in length and stiffness to the myosin levers, and may be particularly suited for myosins working in crowded cellular environments. Although the function of the SAH-domains in human class-6 and class-10 myosins has well been characterised, the distribution of the SAH-domain in all myosin subfamilies and across the eukaryotic tree of life remained elusive. Here, we analysed the largest available myosin sequence dataset consisting of 7919 manually annotated myosin sequences from 938 species representing all major eukaryotic branches using the SAH-prediction algorithm of Waggawagga, a recently developed tool for the identification of SAH-domains. With this approach we identified SAH-domains in more than one third of the supposed 79 myosin subfamilies. Depending on the myosin class, the presence of SAH-domains can range from a few to almost all class members indicating complex patterns of independent and taxon-specific SAH-domain gain and loss.

Waggawagga: comparative visualization of coiled-coil predictions and detection of stable single α-helices (SAH domains).

UNLABELLED: Waggawagga is a web-based tool for the comparative visualization of coiled-coil predictions and the detection of stable single α-helices (SAH domains). Overview schemes show the predicted coiled-coil regions found in the query sequence and provide sliders, which can be used to select segments for detailed helical wheel and helical net views. A window-based score has been developed to predict SAH domains. Export to several bitmap and vector graphics formats is supported. AVAILABILITY AND IMPLEMENTATION: http://waggawagga.motorprotein.de

Spaced words and kmacs: fast alignment-free sequence comparison based on inexact word matches.

In this article, we present a user-friendly web interface for two alignment-free sequence-comparison methods that we recently developed. Most alignment-free methods rely on exact word matches to estimate pairwise similarities or distances between the input sequences. By contrast, our new algorithms are based on inexact word matches. The first of these approaches uses the relative frequencies of so-called spaced words in the input sequences, i.e. words containing 'don't care' or 'wildcard' symbols at certain pre-defined positions. Various distance measures can then be defined on sequences based on their different spaced-word composition. Our second approach defines the distance between two sequences by estimating for each position in the first sequence the length of the longest substring at this position that also occurs in the second sequence with up to k mismatches. Both approaches take a set of deoxyribonucleic acid (DNA) or protein sequences as input and return a matrix of pairwise distance values that can be used as a starting point for clustering algorithms or distance-based phylogeny reconstruction. The two alignment-free programmes are accessible through a web interface at 'Göttingen Bioinformatics Compute Server (GOBICS)': http://spaced.gobics.de http://kmacs.gobics.de and the source codes can be downloaded.

Kassiopeia: a database and web application for the analysis of mutually exclusive exomes of eukaryotes.

BACKGROUND: Alternative splicing is an important process in higher eukaryotes that allows obtaining several transcripts from one gene. A specific case of alternative splicing is mutually exclusive splicing, in which exactly one exon out of a cluster of neighbouring exons is spliced into the mature transcript. Recently, a new algorithm for the prediction of these exons has been developed based on the preconditions that the exons of the cluster have similar lengths, sequence homology, and conserved splice sites, and that they are translated in the same reading frame. DESCRIPTION: In this contribution we introduce Kassiopeia, a database and web application for the generation, storage, and presentation of genome-wide analyses of mutually exclusive exomes. Currently, Kassiopeia provides access to the mutually exclusive exomes of twelve Drosophila species, the thale cress Arabidopsis thaliana, the flatworm Caenorhabditis elegans, and human. Mutually exclusive spliced exons (MXEs) were predicted based on gene reconstructions from Scipio. Based on the standard prediction values, with which 83.5% of the annotated MXEs of Drosophila melanogaster were reconstructed, the exomes contain surprisingly more MXEs than previously supposed and identified. The user can search Kassiopeia using BLAST or browse the genes of each species optionally adjusting the parameters used for the prediction to reveal more divergent or only very similar exon candidates. CONCLUSIONS: We developed a pipeline to predict MXEs in the genomes of several model organisms and a web interface, Kassiopeia, for their visualization. For each gene Kassiopeia provides a comprehensive gene structure scheme, the sequences and predicted secondary structures of the MXEs, and, if available, further evidence for MXE candidates from cDNA/EST data, predictions of MXEs in homologous genes of closely related species, and RNA secondary structure predictions. Kassiopeia can be accessed at http://www.motorprotein.de/kassiopeia.