Safety and tolerability of AMG 330 in adults with relapsed/refractory AML: a phase 1a dose-escalation study.

AMG 330, a bispecific T-cell engager (BiTE®) that binds CD33 and CD3 on T cells facilitates T-cell-mediated cytotoxicity against CD33+ cells. This first-in-human, open-label, dose-escalation study evaluated the safety, pharmacokinetics, pharmacodynamics, and preliminary efficacy of AMG 330 in adults with relapsed/refractory acute myeloid leukemia (R/R AML). Amongst 77 patients treated with AMG 330 (0.5 µg/day-1.6 mg/day) on 14-day or 28-day cycles, maximum tolerated dose was not reached; median duration of treatment was 29 days. The most frequent treatment-related adverse events were cytokine release syndrome (CRS; 78%) and rash (30%); 10% of patients experienced grade 3/4 CRS. CRS was mitigated with stepwise dosing of AMG 330, prophylactic dexamethasone, and early treatment with tocilizumab. Among 60 evaluable patients, eight achieved complete remission or morphologic leukemia-free state; of the 52 non-responders, 37% had ≥50% reduction in AML bone marrow blasts. AMG 330 is a promising CD33-targeted therapeutic strategy for R/R AML.

Single agent subcutaneous blinatumomab for advanced acute lymphoblastic leukemia.

Blinatumomab is a BiTE® (bispecific T-cell engager) molecule that redirects CD3+ T-cells to engage and lyse CD19+ target cells. Here we demonstrate that subcutaneous (SC) blinatumomab can provide high efficacy and greater convenience of administration. In the expansion phase of a multi-institutional phase 1b trial (ClinicalTrials.gov, NCT04521231), heavily pretreated adults with relapsed/refractory B-cell acute lymphoblastic leukemia (R/R B-ALL) received SC blinatumomab at two doses: (1) 250 µg once daily (QD) for week 1 and 500 µg three times weekly (TIW) thereafter (250 µg/500 µg) or (2) 500 µg QD for week 1 and 1000 µg TIW thereafter (500 µg/1000 µg). The primary endpoint was complete remission/complete remission with partial hematologic recovery (CR/CRh) within two cycles. At the data cutoff of September 15, 2023, 29 patients were treated: 14 at the 250 µg/500 µg dose and 13 at 500 µg/1000 µg dose. Data from two ineligible patients were excluded. At the end of two cycles, 12 of 14 patients (85.7%) from the 250 µg/500 µg dose achieved CR/CRh of which nine patients (75.0%) were negative for measurable residual disease (MRD; <10-4 leukemic blasts). At the 500 µg/1000 µg dose, 12 of 13 patients (92.3%) achieved CR/CRh; all 12 patients (100.0%) were MRD-negative. No treatment-related grade 4 cytokine release syndrome (CRS) or neurologic events (NEs) were reported. SC injections were well tolerated and all treatment-related grade 3 CRS and NEs responded to standard-of-care management, interruption, or discontinuation. Treatment with SC blinatumomab resulted in high efficacy, with high MRD-negativity rates and acceptable safety profile in heavily pretreated adults with R/R B-ALL.

Intratumoral follicular regulatory T cells curtail anti-PD-1 treatment efficacy.

Immune-checkpoint blockade (ICB) has shown remarkable clinical success in boosting antitumor immunity. However, the breadth of its cellular targets and specific mode of action remain elusive. We find that tumor-infiltrating follicular regulatory T (TFR) cells are prevalent in tumor tissues of several cancer types. They are primarily located within tertiary lymphoid structures and exhibit superior suppressive capacity and in vivo persistence as compared with regulatory T cells, with which they share a clonal and developmental relationship. In syngeneic tumor models, anti-PD-1 treatment increases the number of tumor-infiltrating TFR cells. Both TFR cell deficiency and the depletion of TFR cells with anti-CTLA-4 before anti-PD-1 treatment improve tumor control in mice. Notably, in a cohort of 271 patients with melanoma, treatment with anti-CTLA-4 followed by anti-PD-1 at progression was associated with better a survival outcome than monotherapy with anti-PD-1 or anti-CTLA-4, anti-PD-1 followed by anti-CTLA-4 at progression or concomitant combination therapy.

Multi-cell type gene coexpression network analysis reveals coordinated interferon response and cross-cell type correlations in systemic lupus erythematosus.

Systemic lupus erythematosus (SLE) is an incurable autoimmune disease disproportionately affecting women. A major obstacle in finding targeted therapies for SLE is its remarkable heterogeneity in clinical manifestations as well as in the involvement of distinct cell types. To identify cell-specific targets as well as cross-correlation relationships among expression programs of different cell types, we here analyze six major circulating immune cell types from SLE patient blood. Our results show that presence of an interferon response signature stratifies patients into two distinct groups (IFNneg vs. IFNpos). Comparing these two groups using differential gene expression and differential gene coexpression analysis, we prioritize a relatively small list of genes from classical monocytes including two known immune modulators: TNFSF13B/BAFF (target of belimumab, an approved therapeutic for SLE) and IL1RN (the basis of anakinra, a therapeutic for rheumatoid arthritis). We then develop a multi-cell type extension of the weighted gene coexpression network analysis (WGCNA) framework, termed mWGCNA. Applying mWGCNA to RNA-seq data from six sorted immune cell populations (15 SLE, 10 healthy donors), we identify a coexpression module with interferon-stimulated genes (ISGs) among all cell types and a cross-cell type correlation linking expression of specific T helper cell markers to B cell response as well as to TNFSF13B expression from myeloid cells, all of which in turn correlates with disease severity of IFNpos patients. Our results demonstrate the power of a hypothesis-free and data-driven approach to discover drug targets and to reveal novel cross-correlation across cell types in SLE with implications for other autoimmune diseases.

Bioinformatics drives discovery in Biomedicine.

Bioinformatics has evolved from providing basic solutions, such as sequence alignment, structure predictions, and phylogenetic analysis to an independent data-driven field. The unprecedented growth of genomic technologies and the enormous data have opened an avenue for bioinformaticians (Bioinformatics professionals) never been seen before in the history of mankind. The novel opportunity also requires creative solutions that need skills to deal with noisy, unstructured information to offer valuable biological insights. Currently, we are seeing only the tip of an iceberg and the future will revolve around big data sets in all forms of biological research. The emerging challenge is to unfold the hidden iceberg of data.

Human Eosinophils Express a Distinct Gene Expression Program in Response to IL-3 Compared with Common ß-Chain Cytokines IL-5 and GM-CSF.

Despite recent advances in asthma management with anti-IL-5 therapies, many patients have eosinophilic asthma that remains poorly controlled. IL-3 shares a common ß subunit receptor with both IL-5 and GM-CSF but, through α-subunit-specific properties, uniquely influences eosinophil biology and may serve as a potential therapeutic target. We aimed to globally characterize the transcriptomic profiles of GM-CSF, IL-3, and IL-5 stimulation on human circulating eosinophils and identify differences in gene expression using advanced statistical modeling. Human eosinophils were isolated from the peripheral blood of healthy volunteers and stimulated with either GM-CSF, IL-3, or IL-5 for 48 h. RNA was then extracted and bulk sequencing performed. DESeq analysis identified differentially expressed genes and weighted gene coexpression network analysis independently defined modules of genes that are highly coexpressed. GM-CSF, IL-3, and IL-5 commonly upregulated 252 genes and downregulated 553 genes, producing a proinflammatory and survival phenotype that was predominantly mediated through TWEAK signaling. IL-3 stimulation yielded the most numbers of differentially expressed genes that were also highly coexpressed (n = 119). These genes were enriched in pathways involving JAK/STAT signaling. GM-CSF and IL-5 stimulation demonstrated redundancy in eosinophil gene expression. In conclusion, IL-3 produces a distinct eosinophil gene expression program among the ß-chain receptor cytokines. IL-3-upregulated genes may provide a foundation for research into therapeutics for patients with eosinophilic asthma who do not respond to anti-IL-5 therapies.

Single-cell transcriptomic analysis of tissue-resident memory T cells in human lung cancer.

High numbers of tissue-resident memory T (TRM) cells are associated with better clinical outcomes in cancer patients. However, the molecular characteristics that drive their efficient immune response to tumors are poorly understood. Here, single-cell and bulk transcriptomic analysis of TRM and non-TRM cells present in tumor and normal lung tissue from patients with lung cancer revealed that PD-1-expressing TRM cells in tumors were clonally expanded and enriched for transcripts linked to cell proliferation and cytotoxicity when compared with PD-1-expressing non-TRM cells. This feature was more prominent in the TRM cell subset coexpressing PD-1 and TIM-3, and it was validated by functional assays ex vivo and also reflected in their chromatin accessibility profile. This PD-1+TIM-3+ TRM cell subset was enriched in responders to PD-1 inhibitors and in tumors with a greater magnitude of CTL responses. These data highlight that not all CTLs expressing PD-1 are dysfunctional; on the contrary, TRM cells with PD-1 expression were enriched for features suggestive of superior functionality.

Reduced expression of phosphatase PTPN2 promotes pathogenic conversion of Tregs in autoimmunity.

Genetic variants at the PTPN2 locus, which encodes the tyrosine phosphatase PTPN2, cause reduced gene expression and are linked to rheumatoid arthritis (RA) and other autoimmune diseases. PTPN2 inhibits signaling through the T cell and cytokine receptors, and loss of PTPN2 promotes T cell expansion and CD4- and CD8-driven autoimmunity. However, it remains unknown whether loss of PTPN2 in FoxP3+ regulatory T cells (Tregs) plays a role in autoimmunity. Here we aimed to model human autoimmune-predisposing PTPN2 variants, the presence of which results in a partial loss of PTPN2 expression, in mouse models of RA. We identified that reduced expression of Ptpn2 enhanced the severity of autoimmune arthritis in the T cell-dependent SKG mouse model and demonstrated that this phenotype was mediated through a Treg-intrinsic mechanism. Mechanistically, we found that through dephosphorylation of STAT3, PTPN2 inhibits IL-6-driven pathogenic loss of FoxP3 after Tregs have acquired RORγt expression, at a stage when chromatin accessibility for STAT3-targeted IL-17-associated transcription factors is maximized. We conclude that PTPN2 promotes FoxP3 stability in mouse RORγt+ Tregs and that loss of function of PTPN2 in Tregs contributes to the association between PTPN2 and autoimmunity.

Accurate prediction of personalized olfactory perception from large-scale chemoinformatic features.

Background: The olfactory stimulus-percept problem has been studied for more than a century, yet it is still hard to precisely predict the odor given the large-scale chemoinformatic features of an odorant molecule. A major challenge is that the perceived qualities vary greatly among individuals due to different genetic and cultural backgrounds. Moreover, the combinatorial interactions between multiple odorant receptors and diverse molecules significantly complicate the olfaction prediction. Many attempts have been made to establish structure-odor relationships for intensity and pleasantness, but no models are available to predict the personalized multi-odor attributes of molecules. In this study, we describe our winning algorithm for predicting individual and population perceptual responses to various odorants in the DREAM Olfaction Prediction Challenge. Results: We find that random forest model consisting of multiple decision trees is well suited to this prediction problem, given the large feature spaces and high variability of perceptual ratings among individuals. Integrating both population and individual perceptions into our model effectively reduces the influence of noise and outliers. By analyzing the importance of each chemical feature, we find that a small set of low- and nondegenerative features is sufficient for accurate prediction. Conclusions: Our random forest model successfully predicts personalized odor attributes of structurally diverse molecules. This model together with the top discriminative features has the potential to extend our understanding of olfactory perception mechanisms and provide an alternative for rational odorant design.

Comparison of spinocerebellar ataxia type 3 mouse models identifies early gain-of-function, cell-autonomous transcriptional changes in oligodendrocytes.

Spinocerebellar ataxia type 3 (SCA3) is a neurodegenerative disorder caused by a polyglutamine-encoding CAG repeat expansion in the ATXN3 gene. This expansion leads to misfolding and aggregation of mutant ataxin-3 (ATXN3) and degeneration of select brain regions. A key unanswered question in SCA3 and other polyglutamine diseases is the extent to which neurodegeneration is mediated through gain-of-function versus loss-of-function. To address this question in SCA3, we performed transcriptional profiling on the brainstem, a highly vulnerable brain region in SCA3, in a series of mouse models with varying degrees of ATXN3 expression and aggregation. We include two SCA3 knock-in mouse models: our previously published model that erroneously harbors a tandem duplicate of the CAG repeat-containing exon, and a corrected model, introduced here. Both models exhibit dose-dependent neuronal accumulation and aggregation of mutant ATXN3, but do not exhibit a behavioral phenotype. We identified a molecular signature that correlates with ATXN3 neuronal aggregation yet is primarily linked to oligodendrocytes, highlighting early white matter dysfunction in SCA3. Two robustly elevated oligodendrocyte transcripts, Acy3 and Tnfrsf13c, were confirmed as elevated at the protein level in SCA3 human disease brainstem. To determine if mutant ATXN3 acts on oligodendrocytes cell-autonomously, we manipulated the repeat expansion in the variant SCA3 knock-in mouse by cell-type specific Cre/LoxP recombination. Changes in oligodendrocyte transcripts are driven cell-autonomously and occur independent of neuronal ATXN3 aggregation. Our findings support a primary toxic gain of function mechanism and highlight a previously unrecognized role for oligodendrocyte dysfunction in SCA3 disease pathogenesis.

GATA3 Abundance Is a Critical Determinant of T Cell Receptor ß Allelic Exclusion.

Allelic exclusion describes the essential immunological process by which feedback repression of sequential DNA rearrangements ensures that only one autosome expresses a functional T or B cell receptor. In wild-type mammals, approximately 60% of cells have recombined the DNA of one T cell receptor ß (TCRß) V-to-DJ-joined allele in a functional configuration, while the second allele has recombined only the DJ sequences; the other 40% of cells have recombined the V to the DJ segments on both alleles, with only one of the two alleles predicting a functional TCRß protein. Here we report that the transgenic overexpression of GATA3 leads predominantly to biallelic TCRß gene (Tcrb) recombination. We also found that wild-type immature thymocytes can be separated into distinct populations based on intracellular GATA3 expression and that GATA3LO cells had almost exclusively recombined only one Tcrb locus (that predicted a functional receptor sequence), while GATA3HI cells had uniformly recombined both Tcrb alleles (one predicting a functional and the other predicting a nonfunctional rearrangement). These data show that GATA3 abundance regulates the recombination propensity at the Tcrb locus and provide new mechanistic insight into the historic immunological conundrum for how Tcrb allelic exclusion is mediated.

miRmine: a database of human miRNA expression profiles.

MOTIVATION: MicroRNAs (miRNAs) are small non-coding RNAs that are involved in post-transcriptional regulation of gene expression. In this high-throughput sequencing era, a tremendous amount of RNA-seq data is accumulating, and full utilization of publicly available miRNA data is an important challenge. These data are useful to determine expression values for each miRNA, but quantification pipelines are in a primitive stage and still evolving; there are many factors that affect expression values significantly. RESULTS: We used 304 high-quality microRNA sequencing (miRNA-seq) datasets from NCBI-SRA and calculated expression profiles for different tissues and cell-lines. In each miRNA-seq dataset, we found an average of more than 500 miRNAs with higher than 5x coverage, and we explored the top five highly expressed miRNAs in each tissue and cell-line. This user-friendly miRmine database has options to retrieve expression profiles of single or multiple miRNAs for a specific tissue or cell-line, either normal or with disease information. Results can be displayed in multiple interactive, graphical and downloadable formats. AVAILABILITY AND IMPLEMENTATION: http://guanlab.ccmb.med.umich.edu/mirmine. CONTACT: bharatpa@umich.edu. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

COMPASS: A computational model to predict changes in MMSE scores 24-months after initial assessment of Alzheimer's disease.

We present COMPASS, a COmputational Model to Predict the development of Alzheimer's diSease Spectrum, to model Alzheimer's disease (AD) progression. This was the best-performing method in recent crowdsourcing benchmark study, DREAM Alzheimer's Disease Big Data challenge to predict changes in Mini-Mental State Examination (MMSE) scores over 24-months using standardized data. In the present study, we conducted three additional analyses beyond the DREAM challenge question to improve the clinical contribution of our approach, including: (1) adding pre-validated baseline cognitive composite scores of ADNI-MEM and ADNI-EF, (2) identifying subjects with significant declines in MMSE scores, and (3) incorporating SNPs of top 10 genes connected to APOE identified from functional-relationship network. For (1) above, we significantly improved predictive accuracy, especially for the Mild Cognitive Impairment (MCI) group. For (2), we achieved an area under ROC of 0.814 in predicting significant MMSE decline: our model has 100% precision at 5% recall, and 91% accuracy at 10% recall. For (3), "genetic only" model has Pearson's correlation of 0.15 to predict progression in the MCI group. Even though addition of this limited genetic model to COMPASS did not improve prediction of progression of MCI group, the predictive ability of SNP information extended beyond well-known APOE allele.

Crowdsourced assessment of common genetic contribution to predicting anti-TNF treatment response in rheumatoid arthritis.

Rheumatoid arthritis (RA) affects millions world-wide. While anti-TNF treatment is widely used to reduce disease progression, treatment fails in ∼one-third of patients. No biomarker currently exists that identifies non-responders before treatment. A rigorous community-based assessment of the utility of SNP data for predicting anti-TNF treatment efficacy in RA patients was performed in the context of a DREAM Challenge (http://www.synapse.org/RA_Challenge). An open challenge framework enabled the comparative evaluation of predictions developed by 73 research groups using the most comprehensive available data and covering a wide range of state-of-the-art modelling methodologies. Despite a significant genetic heritability estimate of treatment non-response trait (h(2)=0.18, P value=0.02), no significant genetic contribution to prediction accuracy is observed. Results formally confirm the expectations of the rheumatology community that SNP information does not significantly improve predictive performance relative to standard clinical traits, thereby justifying a refocusing of future efforts on collection of other data.

Genome-Wide Functional Annotation of Human Protein-Coding Splice Variants Using Multiple Instance Learning.

The vast majority of human multiexon genes undergo alternative splicing and produce a variety of splice variant transcripts and proteins, which can perform different functions. These protein-coding splice variants (PCSVs) greatly increase the functional diversity of proteins. Most functional annotation algorithms have been developed at the gene level; the lack of isoform-level gold standards is an important intellectual limitation for currently available machine learning algorithms. The accumulation of a large amount of RNA-seq data in the public domain greatly increases our ability to examine the functional annotation of genes at isoform level. In the present study, we used a multiple instance learning (MIL)-based approach for predicting the function of PCSVs. We used transcript-level expression values and gene-level functional associations from the Gene Ontology database. A support vector machine (SVM)-based 5-fold cross-validation technique was applied. Comparatively, genes with multiple PCSVs performed better than single PCSV genes, and performance also improved when more examples were available to train the models. We demonstrated our predictions using literature evidence of ADAM15, LMNA/C, and DMXL2 genes. All predictions have been implemented in a web resource called "IsoFunc", which is freely available for the global scientific community through http://guanlab.ccmb.med.umich.edu/isofunc .

Algorithms for modeling global and context-specific functional relationship networks.

Functional genomics has enormous potential to facilitate our understanding of normal and disease-specific physiology. In the past decade, intensive research efforts have been focused on modeling functional relationship networks, which summarize the probability of gene co-functionality relationships. Such modeling can be based on either expression data only or heterogeneous data integration. Numerous methods have been deployed to infer the functional relationship networks, while most of them target the global (non-context-specific) functional relationship networks. However, it is expected that functional relationships consistently reprogram under different tissues or biological processes. Thus, advanced methods have been developed targeting tissue-specific or developmental stage-specific networks. This article brings together the state-of-the-art functional relationship network modeling methods, emphasizes the need for heterogeneous genomic data integration and context-specific network modeling and outlines future directions for functional relationship networks.

Computational Inferences of the Functions of Alternative/Noncanonical Splice Isoforms Specific to HER2+/ER-/PR- Breast Cancers, a Chromosome 17 C-HPP Study.

This study was conducted as a part of the Chromosome-Centric Human Proteome Project (C-HPP) of the Human Proteome Organization. The main objective is to identify and evaluate functionality of a set of specific noncanonical isoforms expressed in HER2-neu positive, estrogen receptor negative (ER-), and progesterone receptor negative (PR-) breast cancers (HER2+/ER-/PR- BC), an aggressive subtype of breast cancers that cause significant morbidity and mortality. We identified 11 alternative splice isoforms that were differentially expressed in HER2+/ER-/PR- BC compared to normal mammary, triple negative breast cancer and triple positive breast cancer tissues (HER2+/ER+/PR+). We used a stringent criterion that differentially expressed noncanonical isoforms (adjusted p value < 0.05) and have to be expressed in all replicates of HER2+/ER-/PR- BC samples, and the trend in differential expression (up or down) is the same in all comparisons. Of the 11 protein isoforms, six were overexpressed in HER2+/ER-/PR- BC. We explored possible functional roles of these six proteins using several complementary computational tools. Biological processes including cell cycle events and glycolysis were linked to four of these proteins. For example, glycolysis was the top ranking functional process for DMXL2 isoform 3, with a fold change of 27 compared to just two for the canonical protein. No previous reports link DMXL2 with any metabolic processes; the canonical protein is known to participate in signaling pathways. Our results clearly indicate distinct functions for the six overexpressed alternative splice isoforms, and these functions could be specific to HER2+/ER-/PR- tumor progression. Further detailed analysis is warranted as these proteins could be explored as potential biomarkers and therapeutic targets for HER2+/ER-/PR- BC patients.

Functional Networks of Highest-Connected Splice Isoforms: From The Chromosome 17 Human Proteome Project.

Alternative splicing allows a single gene to produce multiple transcript-level splice isoforms from which the translated proteins may show differences in their expression and function. Identifying the major functional or canonical isoform is important for understanding gene and protein functions. Identification and characterization of splice isoforms is a stated goal of the HUPO Human Proteome Project and of neXtProt. Multiple efforts have catalogued splice isoforms as "dominant", "principal", or "major" isoforms based on expression or evolutionary traits. In contrast, we recently proposed highest connected isoforms (HCIs) as a new class of canonical isoforms that have the strongest interactions in a functional network and revealed their significantly higher (differential) transcript-level expression compared to nonhighest connected isoforms (NCIs) regardless of tissues/cell lines in the mouse. HCIs and their expression behavior in the human remain unexplored. Here we identified HCIs for 6157 multi-isoform genes using a human isoform network that we constructed by integrating a large compendium of heterogeneous genomic data. We present examples for pairs of transcript isoforms of ABCC3, RBM34, ERBB2, and ANXA7. We found that functional networks of isoforms of the same gene can show large differences. Interestingly, differential expression between HCIs and NCIs was also observed in the human on an independent set of 940 RNA-seq samples across multiple tissues, including heart, kidney, and liver. Using proteomic data from normal human retina and placenta, we showed that HCIs are a promising indicator of expressed protein isoforms exemplified by NUDFB6 and M6PR. Furthermore, we found that a significant percentage (20%, p = 0.0003) of human and mouse HCIs are homologues, suggesting their conservation between species. Our identified HCIs expand the repertoire of canonical isoforms and are expected to facilitate studying main protein products, understanding gene regulation, and possibly evolution. The network is available through our web server as a rich resource for investigating isoform functional relationships (http://guanlab.ccmb.med.umich.edu/hisonet). All MS/MS data were available at ProteomeXchange Web site (http://www.proteomexchange.org) through their identifiers (retina: PXD001242, placenta: PXD000754).

MI-PVT: A Tool for Visualizing the Chromosome-Centric Human Proteome.

We have developed the web-based Michigan Proteome Visualization Tool (MI-PVT) to visualize and compare protein expression and isoform-level function across human chromosomes and tissues (http://guanlab.ccmb.med.umich.edu/mipvt). As proof of principle, we have populated the tool with Human Proteome Map (HPM) data. We were able to observe many biologically interesting features. From the vantage point of our chromosome 17 team, for example, we found more than 300 proteins from chromosome 17 expressed in each of the 30 tissues and cell types studied, with the highest number of expressed proteins being 685 in testis. Comparisons of expression levels across tissues showed low numbers of proteins expressed in esophagus, but esophagus had 12 cytoskeletal proteins coded on chromosome 17 with very high expression (>1000 spectral counts). This customized MI-PVT should be helpful for biologists to browse and study specific proteins and protein data sets across tissues and chromosomes. Users can upload any data of interest in MI-PVT for visualization. Our aim is to integrate extensive mass-spectrometric proteomic data into the tool to facilitate finding chromosome-centric protein expression and correlation across tissues.