Radiation therapy induces immunosenescence mediated by p90RSK.

Radiation therapy (RT) to the chest increases the patients' risk of cardiovascular disease (CVD). A complete understanding of the mechanisms by which RT induces CVD could lead to specific preventive, therapeutic approaches. It is becoming evident that both genotoxic chemotherapy agents and radiation induce mitochondrial dysfunction and cellular senescence. Notably, one of the common phenotypes observed in cancer survivors is accelerated senescence, and immunosenescence is closely related to both cancer risk and CVD development. Therefore, suppression of immunosenescence can be an ideal target to prevent cancer treatment-induced CVD. However, the mechanism(s) by which cancer treatments induce immunosenescence are incompletely characterized. We isolated peripheral blood mononuclear cells (PBMCs) before and 3 months after RT from 16 thoracic cancer patients. We characterized human immune cell lineages and markers of senescence, DNA damage response (DDR), efferocytosis, and determinants of clonal hematopoiesis of indeterminant potential (CHIP), using mass cytometry (CyTOF). We found that the frequency of the B cell subtype was decreased after RT. Unsupervised clustering of the CyTOF data identified 138 functional subsets of PBMCs. Compared with baseline, RT increased TBX21 (T-bet) expression in the largest B cell subset of Ki67-/DNMT3a+naïve B cells, and T-bet expression was correlated with phosphorylation of p90RSK expression. CD38 expression was also increased in naïve B cells (CD27-) and CD8+ effector memory CD45RA T cells (TEMRA). In vitro, we found the critical role of p90RSK activation in upregulating (1) CD38+/T-bet+ memory and naïve B, and myeloid cells, (2) senescence-associated ß-gal staining, and (3) mitochondrial reactive oxygen species (ROS) after ionizing radiation (IR). These data suggest the crucial role of p90RSK activation in immunosenescence. The critical role of p90RSK activation in immune cells and T-bet induction in upregulating atherosclerosis formation has been reported. Furthermore, T-bet directly binds to the CD38 promoter region and upregulates CD38 expression. Since both T-bet and CD38 play a significant role in the process of immunosenescence, our data provide a cellular and molecular mechanism that links RT-induced p90RSK activation and the immunosenescence with T-bet and CD38 induction observed in thoracic cancer patients treated by RT and suggests that targeting the p90RSK/T-bet/CD38 pathway could play a role in preventing the radiation-associated CVD and improving cancer prognosis by inhibiting immunosenescence.

Multimodal Mapping of the Tumor and Peripheral Blood Immune Landscape in Human Pancreatic Cancer.

Pancreatic ductal adenocarcinoma (PDA) is characterized by an immune-suppressive tumor microenvironment that renders it largely refractory to immunotherapy. We implemented a multimodal analysis approach to elucidate the immune landscape in PDA. Using a combination of CyTOF, single-cell RNA sequencing, and multiplex immunohistochemistry on patient tumors, matched blood, and non-malignant samples, we uncovered a complex network of immune-suppressive cellular interactions. These experiments revealed heterogeneous expression of immune checkpoint receptors in individual patient's T cells and increased markers of CD8+ T cell dysfunction in advanced disease stage. Tumor-infiltrating CD8+ T cells had an increased proportion of cells expressing an exhausted expression profile that included upregulation of the immune checkpoint TIGIT, a finding that we validated at the protein level. Our findings point to a profound alteration of the immune landscape of tumors, and to patient-specific immune changes that should be taken into account as combination immunotherapy becomes available for pancreatic cancer.

Single-cell immune landscape of human atherosclerotic plaques.

Atherosclerosis is driven by multifaceted contributions of the immune system within the circulation and at vascular focal sites. However, specific characteristics of dysregulated immune cells within atherosclerotic lesions that lead to clinical events such as ischemic stroke or myocardial infarction are poorly understood. Here, using single-cell proteomic and transcriptomic analyses, we uncovered distinct features of both T cells and macrophages in carotid artery plaques of patients with clinically symptomatic disease (recent stroke or transient ischemic attack) compared to asymptomatic disease (no recent stroke). Plaques from symptomatic patients were characterized by a distinct subset of CD4+ T cells and by T cells that were activated and differentiated. Moreover, some T cell subsets in these plaques presented markers of T cell exhaustion. Additionally, macrophages from these plaques contained alternatively activated phenotypes, including subsets associated with plaque vulnerability. In plaques from asymptomatic patients, T cells and macrophages were activated and displayed evidence of interleukin-1ß signaling. The identification of specific features of innate and adaptive immune cells in plaques that are associated with cerebrovascular events may enable the design of more precisely tailored cardiovascular immunotherapies.

A Modified Injector and Sample Acquisition Protocol Can Improve Data Quality and Reduce Inter-Instrument Variability of the Helios Mass Cytometer.

Mass cytometry is a powerful tool for high-dimensional single cell characterization. Since the introduction of the first commercial CyTOF mass cytometer by DVS Sciences in 2009, mass cytometry technology has matured and become more widely utilized, with sequential platform upgrades designed to address specific limitations and to expand the capabilities of the platform. Fluidigm's third-generation Helios mass cytometer introduced a number of upgrades over the previous CyTOF2. One of these new features is a modified narrow bore sample injector that generates smaller ion clouds, which is expected to improve sensitivity and throughput. However, following rigorous testing, we find that the narrow-bore sample injector may have unintended negative consequences on data quality and result in lower median and higher coefficients of variation in many antibody-associated signal intensities. We describe an alternative Helios acquisition protocol using a wider bore injector, which largely mitigates these data quality issues. We directly compare these two protocols in a multisite study of 10 Helios instruments across 7 institutions and show that the modified protocol improves data quality and reduces interinstrument variability. These findings highlight and address an important source of technical variability in mass cytometry experiments that is of particular relevance in the setting of multicenter studies. © 2019 International Society for Advancement of Cytometry.

Development of a Comprehensive Antibody Staining Database Using a Standardized Analytics Pipeline.

Large-scale immune monitoring experiments (such as clinical trials) are a promising direction for biomarker discovery and responder stratification in immunotherapy. Mass cytometry is one of the tools in the immune monitoring arsenal. We propose a standardized workflow for the acquisition and analysis of large-scale mass cytometry experiments. The workflow includes two-tiered barcoding, a broad lyophilized panel, and the incorporation of a fully automated, cloud-based analysis platform. We applied the workflow to a large antibody staining screen using the LEGENDScreen kit, resulting in single-cell data for 350 antibodies over 71 profiling subsets. The screen recapitulates many known trends in the immune system and reveals potential markers for delineating MAIT cells. Additionally, we examine the effect of fixation on staining intensity and identify several markers where fixation leads to either gain or loss of signal. The standardized workflow can be seamlessly integrated into existing trials. Finally, the antibody staining data set is available as an online resource for researchers who are designing mass cytometry experiments in suspension and tissue.

Protein Barcodes Enable High-Dimensional Single-Cell CRISPR Screens.

CRISPR pools are being widely employed to identify gene functions. However, current technology, which utilizes DNA as barcodes, permits limited phenotyping and bulk-cell resolution. To enable novel screening capabilities, we developed a barcoding system operating at the protein level. We synthesized modules encoding triplet combinations of linear epitopes to generate >100 unique protein barcodes (Pro-Codes). Pro-Code-expressing vectors were introduced into cells and analyzed by CyTOF mass cytometry. Using just 14 antibodies, we detected 364 Pro-Code populations; establishing the largest set of protein-based reporters. By pairing each Pro-Code with a different CRISPR, we simultaneously analyzed multiple phenotypic markers, including phospho-signaling, on dozens of knockouts. Pro-Code/CRISPR screens found two interferon-stimulated genes, the immunoproteasome component Psmb8 and a chaperone Rtp4, are important for antigen-dependent immune editing of cancer cells and identified Socs1 as a negative regulator of Pd-l1. The Pro-Code technology enables simultaneous high-dimensional protein-level phenotyping of 100s of genes with single-cell resolution.

Comprehensive innate immune profiling of chikungunya virus infection in pediatric cases.

Chikungunya virus (CHIKV) is a mosquito-borne alphavirus that causes global epidemics of debilitating disease worldwide. To gain functional insight into the host cellular genes required for virus infection, we performed whole-blood RNA-seq, 37-plex mass cytometry of peripheral blood mononuclear cells (PBMCs), and serum cytokine measurements of acute- and convalescent-phase samples obtained from 42 children naturally infected with CHIKV Semi-supervised classification and clustering of single-cell events into 57 sub-communities of canonical leukocyte phenotypes revealed a monocyte-driven response to acute infection, with the greatest expansions in "intermediate" CD14++CD16+ monocytes and an activated subpopulation of CD14+ monocytes. Increases in acute-phase CHIKV envelope protein E2 expression were highest for monocytes and dendritic cells. Serum cytokine measurements confirmed significant acute-phase upregulation of monocyte chemoattractants. Distinct transcriptomic signatures were associated with infection timepoint, as well as convalescent-phase anti-CHIKV antibody titer, acute-phase viremia, and symptom severity. We present a multiscale network that summarizes all observed modulations across cellular and transcriptomic levels and their interactions with clinical outcomes, providing a uniquely global view of the biomolecular landscape of human CHIKV infection.

Average Overlap Frequency: A simple metric to evaluate staining quality and community identification in high dimensional mass cytometry experiments.

High dimensional cytometry now allows measurement of over 50 parameters in a single sample, and is typically visualized using sophisticated dimensionality-reducing methods and analyzed with automated clustering algorithms. While these tools facilitate the identification and presentation of key findings, it remains challenging to effectively monitor and report the staining quality of individual markers. We present the Average Overlap Frequency (AOF), a simple and efficient metric to evaluate and quantify the robustness of staining and clustering quality in high-dimensional data. We leverage the AOF to compare and determine the optimal storage conditions for stained whole blood samples prior to mass cytometry analysis. We also show that the AOF can be easily incorporated as part of automated analysis pipelines in large scale immune monitoring studies and used to flag and exclude samples with poor staining quality. We propose that the AOF may be incorporated as an essential quality control metric to better identify and report the underlying sample quality in all CyTOF and other high-dimensional cytometry experiments.

Innate Immune Landscape in Early Lung Adenocarcinoma by Paired Single-Cell Analyses.

To guide the design of immunotherapy strategies for patients with early stage lung tumors, we developed a multiscale immune profiling strategy to map the immune landscape of early lung adenocarcinoma lesions to search for tumor-driven immune changes. Utilizing a barcoding method that allows a simultaneous single-cell analysis of the tumor, non-involved lung, and blood cells, we provide a detailed immune cell atlas of early lung tumors. We show that stage I lung adenocarcinoma lesions already harbor significantly altered T cell and NK cell compartments. Moreover, we identified changes in tumor-infiltrating myeloid cell (TIM) subsets that likely compromise anti-tumor T cell immunity. Paired single-cell analyses thus offer valuable knowledge of tumor-driven immune changes, providing a powerful tool for the rational design of immune therapies. VIDEO ABSTRACT.

Host-Protozoan Interactions Protect from Mucosal Infections through Activation of the Inflammasome.

While conventional pathogenic protists have been extensively studied, there is an underappreciated constitutive protist microbiota that is an integral part of the vertebrate microbiome. The impact of these species on the host and their potential contributions to mucosal immune homeostasis remain poorly studied. Here, we show that the protozoan Tritrichomonas musculis activates the host epithelial inflammasome to induce IL-18 release. Epithelial-derived IL-18 promotes dendritic cell-driven Th1 and Th17 immunity and confers dramatic protection from mucosal bacterial infections. Along with its role as a "protistic" antibiotic, colonization with T. musculis exacerbates the development of T-cell-driven colitis and sporadic colorectal tumors. Our findings demonstrate a novel mutualistic host-protozoan interaction that increases mucosal host defenses at the cost of an increased risk of inflammatory disease.

Highly multiplexed profiling of single-cell effector functions reveals deep functional heterogeneity in response to pathogenic ligands.

Despite recent advances in single-cell genomic, transcriptional, and mass-cytometric profiling, it remains a challenge to collect highly multiplexed measurements of secreted proteins from single cells for comprehensive analysis of functional states. Herein, we combine spatial and spectral encoding with polydimethylsiloxane (PDMS) microchambers for codetection of 42 immune effector proteins secreted from single cells, representing the highest multiplexing recorded to date for a single-cell secretion assay. Using this platform to profile differentiated macrophages stimulated with lipopolysaccharide (LPS), the ligand of Toll-like receptor 4 (TLR4), reveals previously unobserved deep functional heterogeneity and varying levels of pathogenic activation. Uniquely protein profiling on the same single cells before and after LPS stimulation identified a role for macrophage inhibitory factor (MIF) to potentiate the activation of LPS-induced cytokine production. Advanced clustering analysis identified functional subsets including quiescent, polyfunctional fully activated, partially activated populations with different cytokine profiles. This population architecture is conserved throughout the cell activation process and prevails as it is extended to other TLR ligands and to primary macrophages derived from a healthy donor. This work demonstrates that the phenotypically similar cell population still exhibits a large degree of intrinsic heterogeneity at the functional and cell behavior level. This technology enables full-spectrum dissection of immune functional states in response to pathogenic or environmental stimulation, and opens opportunities to quantify deep functional heterogeneity for more comprehensive and accurate immune monitoring.

Single-cell trajectory detection uncovers progression and regulatory coordination in human B cell development.

Tissue regeneration is an orchestrated progression of cells from an immature state to a mature one, conventionally represented as distinctive cell subsets. A continuum of transitional cell states exists between these discrete stages. We combine the depth of single-cell mass cytometry and an algorithm developed to leverage this continuum by aligning single cells of a given lineage onto a unified trajectory that accurately predicts the developmental path de novo. Applied to human B cell lymphopoiesis, the algorithm (termed Wanderlust) constructed trajectories spanning from hematopoietic stem cells through to naive B cells. This trajectory revealed nascent fractions of B cell progenitors and aligned them with developmentally cued regulatory signaling including IL-7/STAT5 and cellular events such as immunoglobulin rearrangement, highlighting checkpoints across which regulatory signals are rewired paralleling changes in cellular state. This study provides a comprehensive analysis of human B lymphopoiesis, laying a foundation to apply this approach to other tissues and "corrupted" developmental processes including cancer.

viSNE enables visualization of high dimensional single-cell data and reveals phenotypic heterogeneity of leukemia.

New high-dimensional, single-cell technologies offer unprecedented resolution in the analysis of heterogeneous tissues. However, because these technologies can measure dozens of parameters simultaneously in individual cells, data interpretation can be challenging. Here we present viSNE, a tool that allows one to map high-dimensional cytometry data onto two dimensions, yet conserve the high-dimensional structure of the data. viSNE plots individual cells in a visual similar to a scatter plot, while using all pairwise distances in high dimension to determine each cell's location in the plot. We integrated mass cytometry with viSNE to map healthy and cancerous bone marrow samples. Healthy bone marrow automatically maps into a consistent shape, whereas leukemia samples map into malformed shapes that are distinct from healthy bone marrow and from each other. We also use viSNE and mass cytometry to compare leukemia diagnosis and relapse samples, and to identify a rare leukemia population reminiscent of minimal residual disease. viSNE can be applied to any multi-dimensional single-cell technology.

KinSNP software for homozygosity mapping of disease genes using SNP microarrays.

Consanguineous families affected with a recessive genetic disease caused by homozygotisation of a mutation offer a unique advantage for positional cloning of rare diseases. Homozygosity mapping of patient genotypes is a powerful technique for the identification of the genomic locus harbouring the causing mutation. This strategy relies on the observation that in these patients a large region spanning the disease locus is also homozygous with high probability. The high marker density in single nucleotide polymorphism (SNP) arrays is extremely advantageous for homozygosity mapping. We present KinSNP, a user-friendly software tool for homozygosity mapping using SNP arrays. The software searches for stretches of SNPs which are homozygous to the same allele in all ascertained sick individuals. User-specified parameters control the number of allowed genotyping 'errors' within homozygous blocks. Candidate disease regions are then reported in a detailed, coloured Excel file, along with genotypes of family members and healthy controls. An interactive genome browser has been included which shows homozygous blocks, individual genotypes, genes and further annotations along the chromosomes, with zooming and scrolling capabilities. The software has been used to identify the location of a mutated gene causing insensitivity to pain in a large Bedouin family. KinSNP is freely available from.

Enhancement of beta-sheet assembly by cooperative hydrogen bonds potential.

MOTIVATION: The roughness of energy landscapes is a major obstacle to protein structure prediction, since it forces conformational searches to spend much time struggling to escape numerous traps. Specifically, beta-sheet formation is prone to stray, since many possible combinations of hydrogen bonds are dead ends in terms of beta-sheet assembly. It has been shown that cooperative terms for backbone hydrogen bonds ease this problem by augmenting hydrogen bond patterns that are consistent with beta sheets. Here, we present a novel cooperative hydrogen-bond term that is both effective in promoting beta sheets and computationally efficient. In addition, the new term is differentiable and operates on all-atom protein models. RESULTS: Energy optimization of poly-alanine chains under the new term led to significantly more beta-sheet content than optimization under a non-cooperative term. Furthermore, the optimized structure included very few non-native patterns. AVAILABILITY: The new term is implemented within the MESHI package and is freely available at http://cs.bgu.ac.il/ approximately meshi.

Differentiable, multi-dimensional, knowledge-based energy terms for torsion angle probabilities and propensities.

Rotatable torsion angles are the major degrees of freedom in proteins. Adjacent angles are highly correlated and energy terms that rely on these correlations are intensively used in molecular modeling. However, the utility of torsion based terms is not yet fully exploited. Many of these terms do not capture the full scale of the correlations. Other terms, which rely on lookup tables, cannot be used in the context of force-driven algorithms because they are not fully differentiable. This study aims to extend the usability of torsion terms by presenting a set of high-dimensional and fully-differentiable energy terms that are derived from high-resolution structures. The set includes terms that describe backbone conformational probabilities and propensities, side-chain rotamer probabilities, and an elaborate term that couples all the torsion angles within the same residue. The terms are constructed by cubic spline interpolation with periodic boundary conditions that enable full differentiability and high computational efficiency. We show that the spline implementation does not compromise the accuracy of the original database statistics. We further show that the side-chain relevant terms are compatible with established rotamer probabilities. Despite their very local characteristics, the new terms are often able to identify native and native-like structures within decoy sets. Finally, force-based minimization of NMR structures with the new terms improves their torsion angle statistics with minor structural distortion (0.5 A RMSD on average). The new terms are freely available in the MESHI molecular modeling package. The spline coefficients are also available as a documented MATLAB file.