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.

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.

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.

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.

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.

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.

WebScipio: Reconstructing alternative splice variants of eukaryotic proteins.

Accurate exon-intron structures are essential prerequisites in genomics, proteomics and for many protein family and single gene studies. We originally developed Scipio and the corresponding web service WebScipio for the reconstruction of gene structures based on protein sequences and available genome assemblies. WebScipio also allows predicting mutually exclusive spliced exons and tandemly arrayed gene duplicates. The obtained gene structures are illustrated in graphical schemes and can be analysed down to the nucleotide level. The set of eukaryotic genomes available at the WebScipio server is updated on a daily basis. The current version of the web server provides access to ∼3400 genome assembly files of >1100 sequenced eukaryotic species. Here, we have also extended the functionality by adding a module with which expressed sequence tag (EST) and cDNA data can be mapped to the reconstructed gene structure for the identification of all types of alternative splice variants. WebScipio has a user-friendly web interface, and we believe that the improved web server will provide better service to biologists interested in the gene structure corresponding to their protein of interest, including all types of alternative splice forms and tandem gene duplicates. WebScipio is freely available at http://www.webscipio.org.

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.

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.

Actomyosin organelle functions of SPIRE actin nucleators precede animal evolution.

An important question in cell biology is how cytoskeletal proteins evolved and drove the development of novel structures and functions. Here we address the origin of SPIRE actin nucleators. Mammalian SPIREs work with RAB GTPases, formin (FMN)-subgroup actin assembly proteins and class-5 myosin (MYO5) motors to transport organelles along actin filaments towards the cell membrane. However, the origin and extent of functional conservation of SPIRE among species is unknown. Our sequence searches show that SPIRE exist throughout holozoans (animals and their closest single-celled relatives), but not other eukaryotes. SPIRE from unicellular holozoans (choanoflagellate), interacts with RAB, FMN and MYO5 proteins, nucleates actin filaments and complements mammalian SPIRE function in organelle transport. Meanwhile SPIRE and MYO5 proteins colocalise to organelles in Salpingoeca rosetta choanoflagellates. Based on these observations we propose that SPIRE originated in unicellular ancestors of animals providing an actin-myosin driven exocytic transport mechanism that may have contributed to the evolution of complex multicellular animals.

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.

Predicting mutually exclusive spliced exons based on exon length, splice site and reading frame conservation, and exon sequence homology.

BACKGROUND: Alternative splicing of pre-mature RNA is an important process eukaryotes utilize to increase their repertoire of different protein products. Several types of different alternative splice forms exist including exon skipping, differential splicing of exons at their 3'- or 5'-end, intron retention, and mutually exclusive splicing. The latter term is used for clusters of internal exons that are spliced in a mutually exclusive manner. RESULTS: We have implemented an extension to the WebScipio software to search for mutually exclusive exons. Here, the search is based on the precondition that mutually exclusive exons encode regions of the same structural part of the protein product. This precondition provides restrictions to the search for candidate exons concerning their length, splice site conservation and reading frame preservation, and overall homology. Mutually exclusive exons that are not homologous and not of about the same length will not be found. Using the new algorithm, mutually exclusive exons in several example genes, a dynein heavy chain, a muscle myosin heavy chain, and Dscam were correctly identified. In addition, the algorithm was applied to the whole Drosophila melanogaster X chromosome and the results were compared to the Flybase annotation and an ab initio prediction. Clusters of mutually exclusive exons might be subsequent to each other and might encode dozens of exons. CONCLUSIONS: This is the first implementation of an automatic search for mutually exclusive exons in eukaryotes. Exons are predicted and reconstructed in the same run providing the complete gene structure for the protein query of interest. WebScipio offers high quality gene structure figures with the clusters of mutually exclusive exons colour-coded, and several analysis tools for further manual inspection. The genome scale analysis of all genes of the Drosophila melanogaster X chromosome showed that WebScipio is able to find all but two of the 28 annotated mutually exclusive spliced exons and predicts 39 new candidate exons. Thus, WebScipio should be able to identify mutually exclusive spliced exons in any query sequence from any species with a very high probability. WebScipio is freely available to academics at http://www.webscipio.org.

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.

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.

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.

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.

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.