The intrinsic substrate specificity of the human tyrosine kinome.

Phosphorylation of proteins on tyrosine (Tyr) residues evolved in metazoan organisms as a mechanism of coordinating tissue growth1. Multicellular eukaryotes typically have more than 50 distinct protein Tyr kinases that catalyse the phosphorylation of thousands of Tyr residues throughout the proteome1-3. How a given Tyr kinase can phosphorylate a specific subset of proteins at unique Tyr sites is only partially understood4-7. Here we used combinatorial peptide arrays to profile the substrate sequence specificity of all human Tyr kinases. Globally, the Tyr kinases demonstrate considerable diversity in optimal patterns of residues surrounding the site of phosphorylation, revealing the functional organization of the human Tyr kinome by substrate motif preference. Using this information, Tyr kinases that are most compatible with phosphorylating any Tyr site can be identified. Analysis of mass spectrometry phosphoproteomic datasets using this compendium of kinase specificities accurately identifies specific Tyr kinases that are dysregulated in cells after stimulation with growth factors, treatment with anti-cancer drugs or expression of oncogenic variants. Furthermore, the topology of known Tyr signalling networks naturally emerged from a comparison of the sequence specificities of the Tyr kinases and the SH2 phosphotyrosine (pTyr)-binding domains. Finally we show that the intrinsic substrate specificity of Tyr kinases has remained fundamentally unchanged from worms to humans, suggesting that the fidelity between Tyr kinases and their protein substrate sequences has been maintained across hundreds of millions of years of evolution.

Transcriptomic Maps of Colorectal Liver Metastasis: Machine Learning of Gene Activation Patterns and Epigenetic Trajectories in Support of Precision Medicine.

The molecular mechanisms of the liver metastasis of colorectal cancer (CRLM) remain poorly understood. Here, we applied machine learning and bioinformatics trajectory inference to analyze a gene expression dataset of CRLM. We studied the co-regulation patterns at the gene level, the potential paths of tumor development, their functional context, and their prognostic relevance. Our analysis confirmed the subtyping of five liver metastasis subtypes (LMS). We provide gene-marker signatures for each LMS, and a comprehensive functional characterization that considers both the hallmarks of cancer and the tumor microenvironment. The ordering of CRLMs along a pseudotime-tree revealed a continuous shift in expression programs, suggesting a developmental relationship between the subtypes. Notably, trajectory inference and personalized analysis discovered a range of epigenetic states that shape and guide metastasis progression. By constructing prognostic maps that divided the expression landscape into regions associated with favorable and unfavorable prognoses, we derived a prognostic expression score. This was associated with critical processes such as epithelial-mesenchymal transition, treatment resistance, and immune evasion. These factors were associated with responses to neoadjuvant treatment and the formation of an immuno-suppressive, mesenchymal state. Our machine learning-based molecular profiling provides an in-depth characterization of CRLM heterogeneity with possible implications for treatment and personalized diagnostics.

The Incidence and Severity of COVID-19 Infection Post Vaccination in Saudi Arabia.

Background Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of coronavirus disease 2019 (COVID-19). Presently, there is ongoing continuous research for more therapeutic options with a wide variety of vaccine availability. However, many people have worried about the vaccine's side effects. Hence, the current study was conducted to determine the prevalence of vaccinated individuals, side effects, and the rate of infectivity post vaccination including the three doses of vaccinations. Methods A cross-sectional questionnaire-based survey was conducted using Google Forms (Google, Inc., Mountain View, CA). Five hundred forty-three individuals participated and reported their status of COVID-19 infection, vaccination, and side effects. All the participants from Saudi Arabia received all the vaccine shots including the booster dose. Results Most of the Saudi nationals were fully vaccinated, and most received Pfizer vaccines for their first and second shots. Pain at the injection site was reported as the most common adverse effect followed by fever, headache, fatigue, and joint pain. Conclusion From the findings, it is concluded that most of the population of Saudi Arabia was vaccinated effectively. Pain at the injection site is identified as the primary adverse effect of vaccination. Most of the population is vaccinated with the Pfizer vaccine. Long-term side effect monitoring is recommended with large population studies to confirm the status of vaccines and adverse effects.

An atlas of substrate specificities for the human serine/threonine kinome.

Protein phosphorylation is one of the most widespread post-translational modifications in biology1,2. With advances in mass-spectrometry-based phosphoproteomics, 90,000 sites of serine and threonine phosphorylation have so far been identified, and several thousand have been associated with human diseases and biological processes3,4. For the vast majority of phosphorylation events, it is not yet known which of the more than 300 protein serine/threonine (Ser/Thr) kinases encoded in the human genome are responsible3. Here we used synthetic peptide libraries to profile the substrate sequence specificity of 303 Ser/Thr kinases, comprising more than 84% of those predicted to be active in humans. Viewed in its entirety, the substrate specificity of the kinome was substantially more diverse than expected and was driven extensively by negative selectivity. We used our kinome-wide dataset to computationally annotate and identify the kinases capable of phosphorylating every reported phosphorylation site in the human Ser/Thr phosphoproteome. For the small minority of phosphosites for which the putative protein kinases involved have been previously reported, our predictions were in excellent agreement. When this approach was applied to examine the signalling response of tissues and cell lines to hormones, growth factors, targeted inhibitors and environmental or genetic perturbations, it revealed unexpected insights into pathway complexity and compensation. Overall, these studies reveal the intrinsic substrate specificity of the human Ser/Thr kinome, illuminate cellular signalling responses and provide a resource to link phosphorylation events to biological pathways.

In Vivo Modeling of CLL Transformation to Richter Syndrome Reveals Convergent Evolutionary Paths and Therapeutic Vulnerabilities.

Transformation to aggressive disease histologies generates formidable clinical challenges across cancers, but biological insights remain few. We modeled the genetic heterogeneity of chronic lymphocytic leukemia (CLL) through multiplexed in vivo CRISPR-Cas9 B-cell editing of recurrent CLL loss-of-function drivers in mice and recapitulated the process of transformation from indolent CLL into large cell lymphoma [i.e., Richter syndrome (RS)]. Evolutionary trajectories of 64 mice carrying diverse combinatorial gene assortments revealed coselection of mutations in Trp53, Mga, and Chd2 and the dual impact of clonal Mga/Chd2 mutations on E2F/MYC and interferon signaling dysregulation. Comparative human and murine RS analyses demonstrated tonic PI3K signaling as a key feature of transformed disease, with constitutive activation of the AKT and S6 kinases, downmodulation of the PTEN phosphatase, and convergent activation of MYC/PI3K transcriptional programs underlying enhanced sensitivity to MYC/mTOR/PI3K inhibition. This robust experimental system presents a unique framework to study lymphoid biology and therapy. SIGNIFICANCE: Mouse models reflective of the genetic complexity and heterogeneity of human tumors remain few, including those able to recapitulate transformation to aggressive disease histologies. Herein, we model CLL transformation into RS through multiplexed in vivo gene editing, providing key insight into the pathophysiology and therapeutic vulnerabilities of transformed disease. This article is highlighted in the In This Issue feature, p. 101.

Activation of Notch and Myc Signaling via B-cell-Restricted Depletion of Dnmt3a Generates a Consistent Murine Model of Chronic Lymphocytic Leukemia.

Chronic lymphocytic leukemia (CLL) is characterized by disordered DNA methylation, suggesting these epigenetic changes might play a critical role in disease onset and progression. The methyltransferase DNMT3A is a key regulator of DNA methylation. Although DNMT3A somatic mutations in CLL are rare, we found that low DNMT3A expression is associated with more aggressive disease. A conditional knockout mouse model showed that homozygous depletion of Dnmt3a from B cells results in the development of CLL with 100% penetrance at a median age of onset of 5.3 months, and heterozygous Dnmt3a depletion yields a disease penetrance of 89% with a median onset at 18.5 months, confirming its role as a haploinsufficient tumor suppressor. B1a cells were confirmed as the cell of origin of disease in this model, and Dnmt3a depletion resulted in focal hypomethylation and activation of Notch and Myc signaling. Amplification of chromosome 15 containing the Myc gene was detected in all CLL mice tested, and infiltration of high-Myc-expressing CLL cells in the spleen was observed. Notably, hyperactivation of Notch and Myc signaling was exclusively observed in the Dnmt3a CLL mice, but not in three other CLL mouse models tested (Sf3b1-Atm, Ikzf3, and MDR), and Dnmt3a-depleted CLL were sensitive to pharmacologic inhibition of Notch signaling in vitro and in vivo. Consistent with these findings, human CLL samples with lower DNMT3A expression were more sensitive to Notch inhibition than those with higher DNMT3A expression. Altogether, these results suggest that Dnmt3a depletion induces CLL that is highly dependent on activation of Notch and Myc signaling. SIGNIFICANCE: Loss of DNMT3A expression is a driving event in CLL and is associated with aggressive disease, activation of Notch and Myc signaling, and enhanced sensitivity to Notch inhibition.

Phase II study of pembrolizumab in refractory esophageal cancer with correlates of response and survival.

BACKGROUND: Immune checkpoint inhibitors have revolutionized cancer treatment, but the benefits in refractory patients with esophageal cancer have been modest. Predictors of response as well as new targets for novel therapeutic combinations are needed. In this phase 2 clinical trial, we tested single-agent pembrolizumab in patients with advanced esophageal cancer, who received at least one prior line of therapy. METHODS: Pembrolizumab 200 mg every 3 weeks was tested in 49 patients with refractory esophageal cancer: 39 with adenocarcinoma and 10 with esophageal squamous cell carcinoma. Major endpoints were radiological response by Immune-related Response Evaluation Criteria In Solid Tumors and survival. Tumor samples were evaluated for programmed cell death ligand 1 (PD-L1) expression, tumor mutational burden (TMB), and immune contexture by both NanoString mRNA expression analysis and flow cytometry. Peripheral blood mononuclear cells and a panel of circulating chemokines were also analyzed. RESULTS: The overall response rate (ORR) was 8% (4 of 49 patients; 95% CI 2.3% to 19.6%). Median overall survival (OS) was 5.8 months (95% CI 4.0 to 9.5). ORR and OS were not associated with histology. For PD-L1-positive patients, ORR was 13.3% (95% CI 1.7% to 40.5%) and median OS was 7.9 months (95% CI 4.7 to 15.5). A trend toward improved OS was observed in seven patients with a TMB ≥10 mut/Mb (p=0.086). Tumors with a PD-L1 Combined Positive Score ≥1 showed enrichment of LAG3 (p=0.005) and IDO1 (p=0.04) gene expression. Baseline levels of circulating CXCL10, interleukin 2 (IL2) receptor α (IL2RA) and IL6 were associated with survival: CXCL10 favorably, (HR 0.37, p=0.002 (progression-free survival); HR 0.55, p=0.018 (OS)); IL2RA and IL6 unfavorably (HR 1.57, p=0.020 for IL6 (OS); HR 2.36, p=0.025 for IL2RA (OS)). CONCLUSIONS: Pembrolizumab monotherapy was modestly effective in refractory esophageal cancer. Circulating CXCL10 at baseline appeared to be a robust predictor of response. Other T cell exhaustion markers are upregulated in PD-L1-positive patients, suggesting that immunotherapy combinations such as anti-LAG3/programmed cell death protein 1 (PD-1) or anti-IDO1/PD-1 may be of promise in refractory esophageal cancer.

A hotspot mutation in transcription factor IKZF3 drives B cell neoplasia via transcriptional dysregulation.

Hotspot mutation of IKZF3 (IKZF3-L162R) has been identified as a putative driver of chronic lymphocytic leukemia (CLL), but its function remains unknown. Here, we demonstrate its driving role in CLL through a B cell-restricted conditional knockin mouse model. Mutant Ikzf3 alters DNA binding specificity and target selection, leading to hyperactivation of B cell receptor (BCR) signaling, overexpression of nuclear factor κB (NF-κB) target genes, and development of CLL-like disease in elderly mice with a penetrance of ~40%. Human CLL carrying either IKZF3 mutation or high IKZF3 expression was associated with overexpression of BCR/NF-κB pathway members and reduced sensitivity to BCR signaling inhibition by ibrutinib. Our results thus highlight IKZF3 oncogenic function in CLL via transcriptional dysregulation and demonstrate that this pro-survival function can be achieved by either somatic mutation or overexpression of this CLL driver. This emphasizes the need for combinatorial approaches to overcome IKZF3-mediated BCR inhibitor resistance.

Personal neoantigen vaccines induce persistent memory T cell responses and epitope spreading in patients with melanoma.

Personal neoantigen vaccines have been envisioned as an effective approach to induce, amplify and diversify antitumor T cell responses. To define the long-term effects of such a vaccine, we evaluated the clinical outcome and circulating immune responses of eight patients with surgically resected stage IIIB/C or IVM1a/b melanoma, at a median of almost 4 years after treatment with NeoVax, a long-peptide vaccine targeting up to 20 personal neoantigens per patient ( NCT01970358 ). All patients were alive and six were without evidence of active disease. We observed long-term persistence of neoantigen-specific T cell responses following vaccination, with ex vivo detection of neoantigen-specific T cells exhibiting a memory phenotype. We also found diversification of neoantigen-specific T cell clones over time, with emergence of multiple T cell receptor clonotypes exhibiting distinct functional avidities. Furthermore, we detected evidence of tumor infiltration by neoantigen-specific T cell clones after vaccination and epitope spreading, suggesting on-target vaccine-induced tumor cell killing. Personal neoantigen peptide vaccines thus induce T cell responses that persist over years and broaden the spectrum of tumor-specific cytotoxicity in patients with melanoma.

Cross-Site Concordance Evaluation of Tumor DNA and RNA Sequencing Platforms for the CIMAC-CIDC Network.

PURPOSE: Whole-exome (WES) and RNA sequencing (RNA-seq) are key components of cancer immunogenomic analyses. To evaluate the consistency of tumor WES and RNA-seq profiling platforms across different centers, the Cancer Immune Monitoring and Analysis Centers (CIMAC) and the Cancer Immunologic Data Commons (CIDC) conducted a systematic harmonization study. EXPERIMENTAL DESIGN: DNA and RNA were centrally extracted from fresh frozen and formalin-fixed paraffin-embedded non-small cell lung carcinoma tumors and distributed to three centers for WES and RNA-seq profiling. In addition, two 10-plex HapMap cell line pools with known mutations were used to evaluate the accuracy of the WES platforms. RESULTS: The WES platforms achieved high precision (> 0.98) and recall (> 0.87) on the HapMap pools when evaluated on loci using > 50× common coverage. Nonsynonymous mutations clustered by tumor sample, achieving an index of specific agreement above 0.67 among replicates, centers, and sample processing. A DV200 > 24% for RNA, as a putative presequencing RNA quality control (QC) metric, was found to be a reliable threshold for generating consistent expression readouts in RNA-seq and NanoString data. MedTIN > 30 was likewise assessed as a reliable RNA-seq QC metric, above which samples from the same tumor across replicates, centers, and sample processing runs could be robustly clustered and HLA typing, immune infiltration, and immune repertoire inference could be performed. CONCLUSIONS: The CIMAC collaborating laboratory platforms effectively generated consistent WES and RNA-seq data and enable robust cross-trial comparisons and meta-analyses of highly complex immuno-oncology biomarker data across the NCI CIMAC-CIDC Network.

Correction to: Angiosarcoma patients treated with immune checkpoint inhibitors: a case series of seven patients from a single institution.

Following publication of the original article [1], the authors have reported that the following sentence "While of the same IgG1 class as ipilimumab, preclinical data suggests this molecule may have enhanced activity against T regulatory cells".

Angiosarcoma patients treated with immune checkpoint inhibitors: a case series of seven patients from a single institution.

BACKGROUND: Angiosarcoma is an uncommon endothelial malignancy and a highly aggressive soft tissue sarcoma. Due to its infiltrative nature, successful management of localized angiosarcoma is often challenging. Systemic chemotherapy is used in the metastatic setting and occasionally in patients with high-risk localized disease in neoadjuvant or adjuvant settings. However, responses tend to be short-lived and most patients succumb to metastatic disease. Novel therapies are needed for patients with angiosarcomas. METHODS: We performed a retrospective analysis of patients with locally advanced or metastatic angiosarcoma, who were treated with checkpoint inhibitors at our institution. We collected their clinical information and outcome measurements. In one patient with achieved complete response, we analyzed circulating and infiltrating T cells within peripheral blood and tumor tissue. RESULTS: We have treated seven angiosarcoma (AS) patients with checkpoint inhibitors either in the context of clinical trials or off label [Pembrolizumab + Axitinib (NCT02636725; n = 1), AGEN1884, a CTLA-4 inhibitor (NCT02694822; n = 2), Pembrolizumab (n = 4)]. Five patients had cutaneous angiosarcoma, one primary breast angiosarcoma and one radiation-associated breast angiosarcoma. At 12 weeks, 5/7 patients (71%) had partial response of their lesions either on imaging and/or clinical exam and two (29%) had progressive disease. 6/7 patients are alive to date and, thus far, 3/7 patients (43%) have progressed (median 3.4 months)- one achieved partial response after pembrolizumab was switched to ongoing Nivolumab/Ipilimumab, one died of progressive disease at 31 weeks (primary breast angiosarcoma) and one was placed on pazopanib. One patient had a complete response (CR) following extended treatment with monotherapy AGEN1884. No patient experienced any ≥ grade 2 toxicities. CONCLUSIONS: This case series underscores the value of targeted immunotherapy in treating angiosarcoma. It also identifies genetic heterogeneity of cutaneous angiosarcomas and discusses specific genetic findings that may explain reported benefits from immunotherapy.

Mutation-Derived Neoantigens for Cancer Immunotherapy.

Mutation-derived neoantigens distinguish tumor from normal cells. T cells can sense the HLA-presented mutations, recognize tumor cells as non-self and destroy them. Therapeutically, immunotherapy antibodies can increase the virulence of the immune system by increasing T-cell cytotoxicity targeted toward neoantigens. Neoantigen vaccines act through antigen-presenting cells, such as dendritic cells, to activate patient-endogenous T cells that recognize vaccine-encoded mutations. Infusion of mutation-targeting T cells by adoptive cell therapy (ACT) directly increases the number and frequency of cytotoxic T cells recognizing and killing tumor cells. At the same time, publicly-funded consortia have profiled tumor genomes across many indications, identifying mutations in each tumor. For example, we find basal and HER2 positive tumors contain more mutated proteins and more TP53 mutations than luminal A/B breast tumors. HPV negative tumors have more mutated proteins than HPV positive head and neck tumors and in agreement with the hypothesis that HPV activity interferes with p53 activity, only 14% of the HPV positive mutations have TP53 mutations vs. 86% of the HPV negative tumors. Lung adenocarcinomas in smokers have over four times more mutated proteins relative to those in never smokers (median 248 vs. 61, respectively). With an eye toward immunotherapy applications, we review the spectrum of mutations in multiple indications, show variations in indication sub-types, and examine intra- and inter-indication prevalence of re-occurring mutation neoantigens that could be used for warehouse vaccines and ACT.

Spatiotemporal segregation of human marginal zone and memory B cell populations in lymphoid tissue.

Human memory B cells and marginal zone (MZ) B cells share common features such as the expression of CD27 and somatic mutations in their IGHV and BCL6 genes, but the relationship between them is controversial. Here, we show phenotypic progression within lymphoid tissues as MZ B cells emerge from the mature naïve B cell pool via a precursor CD27-CD45RBMEM55+ population distant from memory cells. By imaging mass cytometry, we find that MZ B cells and memory B cells occupy different microanatomical niches in organised gut lymphoid tissues. Both populations disseminate widely between distant lymphoid tissues and blood, and both diversify their IGHV repertoire in gut germinal centres (GC), but nevertheless remain largely clonally separate. MZ B cells are therefore not developmentally contiguous with or analogous to classical memory B cells despite their shared ability to transit through GC, where somatic mutations are acquired.

MicroRNA signatures discriminate between uterine and ovarian serous carcinomas.

Synchronous endometrial and ovarian malignancies occur in 5% of women presenting with endometrial cancer and 10% of patients presenting with ovarian malignancy. When a high-grade serous carcinoma concurrently involves both ovary and endometrium, pathological determination of whether they are synchronous primaries or metastatic tumors from one primary site can be challenging. MicroRNAs (miRNA) are 22-nucleotide noncoding RNAs that are aberrantly expressed in cancer cells and may inherit their cellular lineage characteristics. We explored possible differential miRNA signatures that may separate high-grade ovarian serous carcinoma from primary endometrial serous carcinoma. Forty-seven samples of histologically pure high-grade serous carcinoma of both uterine (16 case) and ovarian primaries (31 cases) were included. Expression of 384 mature miRNAs was analyzed using ABI TaqMan Low-Density Arrays technology. A random forest model was used to identify miRNAs that together could differentiate between uterine and ovarian serous carcinomas. Among 150 miRNAs detectable at various levels in the study cases, a panel of 11-miRNA signatures was identified to significantly discriminate between ovarian and uterine serous carcinoma (P < .05). A nested cross-validated convergent forest plot using 6 of the 11 miRNA signature was eventually established to classify the tumors with 91.5% accuracy. In conclusion, we have characterized a miRNA signature panel in this exploratory study that shows significant discriminatory power in separating primary ovarian high-grade serous carcinoma from its endometrial counterpart.

VDJML: a file format with tools for capturing the results of inferring immune receptor rearrangements.

BACKGROUND: The genes that produce antibodies and the immune receptors expressed on lymphocytes are not germline encoded; rather, they are somatically generated in each developing lymphocyte by a process called V(D)J recombination, which assembles specific, independent gene segments into mature composite genes. The full set of composite genes in an individual at a single point in time is referred to as the immune repertoire. V(D)J recombination is the distinguishing feature of adaptive immunity and enables effective immune responses against an essentially infinite array of antigens. Characterization of immune repertoires is critical in both basic research and clinical contexts. Recent technological advances in repertoire profiling via high-throughput sequencing have resulted in an explosion of research activity in the field. This has been accompanied by a proliferation of software tools for analysis of repertoire sequencing data. Despite the widespread use of immune repertoire profiling and analysis software, there is currently no standardized format for output files from V(D)J analysis. Researchers utilize software such as IgBLAST and IMGT/High V-QUEST to perform V(D)J analysis and infer the structure of germline rearrangements. However, each of these software tools produces results in a different file format, and can annotate the same result using different labels. These differences make it challenging for users to perform additional downstream analyses. RESULTS: To help address this problem, we propose a standardized file format for representing V(D)J analysis results. The proposed format, VDJML, provides a common standardized format for different V(D)J analysis applications to facilitate downstream processing of the results in an application-agnostic manner. The VDJML file format specification is accompanied by a support library, written in C++ and Python, for reading and writing the VDJML file format. CONCLUSIONS: The VDJML suite will allow users to streamline their V(D)J analysis and facilitate the sharing of scientific knowledge within the community. The VDJML suite and documentation are available from https://vdjserver.org/vdjml/ . We welcome participation from the community in developing the file format standard, as well as code contributions.

Protective Role of Ramipril and Candesartan against Myocardial Ischemic Reperfusion Injury: A Biochemical and Transmission Electron Microscopical Study.

The present study was designed to investigate the role of combined administration of Ramipril and Candesartan against in vitro myocardial ischemic reperfusion injury in rat. Male Wistar albino rats were divided into five groups (n = 6) and treated with saline (10 mL/kg), Ramipril (2 mg/kg), Candesartan (1 mg/kg), and the combination of both drugs, respectively 24 h before induction of global ischemia (5 min of stabilization, 9 min of global ischemia, and 12 min of reflow). Combination of Ramipril and Candesartan when compared to the monotherapy significantly increased the levels of superoxide dismutase, reduced glutathione, catalase, and nitric oxide and decreased the levels of thiobarbituric acid reactive substances. In addition, the superior protective role of combination of Ramipril and Candesartan on ischemia induced myocardial damage was further confirmed by well preserved myocardial tissue architecture in light microscopy and transmission electron microscopy analysis studies. The combination was proved to be effective in salvaging the myocardial tissue against ischemic reperfusion injury when compared to the monotherapy of individual drugs and further investigations on protective mechanism of drugs by increasing the nitric oxide level at molecular levels are needed.

Characterization of Diabetogenic CD8+ T Cells: IMMUNE THERAPY WITH METABOLIC BLOCKADE.

Type 1 diabetes mellitus is caused by the killing of insulin-producing ß cells by CD8+T cells. The disease progression, which is chronic, does not follow a course like responses to conventional antigens such as viruses, but accelerates as glucose tolerance deteriorates. To identify the unique features of the autoimmune effectors that may explain this behavior, we analyzed diabetogenic CD8+ T cells that recognize a peptide from the diabetes antigen IGRP (NRP-V7-reactive) in prediabetic NOD mice and compared them to others that shared their phenotype (CD44(+)CD62L(lo)PD-1(+)CXCR3(+)) but negative for diabetes antigen tetramers and to LCMV (lymphocytic choriomeningitis)-reactive CD8+ T cells. There was an increase in the frequency of the NRP-V7-reactive cells coinciding with the time of glucose intolerance. The T cells persisted in hyperglycemic NOD mice maintained with an insulin pellet despite destruction of ß cells. We compared gene expression in the three groups of cells compared with the other two subsets of cells, and the NRP-V7-reactive cells exhibited gene expression of memory precursor effector cells. They had reduced cellular proliferation and were less dependent on oxidative phosphorylation. When prediabetic NOD mice were treated with 2-deoxyglucose to block aerobic glycolysis, there was a reduction in the diabetes antigen versus other cells of similar phenotype and loss of lymphoid cells infiltrating the islets. In addition, treatment of NOD mice with 2-deoxyglucose resulted in improved ß cell granularity. These findings identify a link between metabolic disturbances and autoreactive T cells that promotes development of autoimmune diabetes.

Change-O: a toolkit for analyzing large-scale B cell immunoglobulin repertoire sequencing data.

UNLABELLED: Advances in high-throughput sequencing technologies now allow for large-scale characterization of B cell immunoglobulin (Ig) repertoires. The high germline and somatic diversity of the Ig repertoire presents challenges for biologically meaningful analysis, which requires specialized computational methods. We have developed a suite of utilities, Change-O, which provides tools for advanced analyses of large-scale Ig repertoire sequencing data. Change-O includes tools for determining the complete set of Ig variable region gene segment alleles carried by an individual (including novel alleles), partitioning of Ig sequences into clonal populations, creating lineage trees, inferring somatic hypermutation targeting models, measuring repertoire diversity, quantifying selection pressure, and calculating sequence chemical properties. All Change-O tools utilize a common data format, which enables the seamless integration of multiple analyses into a single workflow. AVAILABILITY AND IMPLEMENTATION: Change-O is freely available for non-commercial use and may be downloaded from http://clip.med.yale.edu/changeo. CONTACT: steven.kleinstein@yale.edu.