Graph Fourier transform for spatial omics representation and analyses of complex organs

Spatial omics technologies are capable of deciphering detailed components of complex organs or tissue in cellular and subcellular resolution. A robust, interpretable, and unbiased representation method for spatial omics is necessary to illuminate novel investigations into biological functions, whereas a mathematical theory deficiency still exists. We present SpaGFT (Spatial Graph Fourier Transform), which provides a unique analytical feature representation of spatial omics data and elucidates molecular signatures linked to critical biological processes within tissues and cells. It outperformed existing tools in spatially variable gene prediction and gene expression imputation across human/mouse Visium data. Integrating SpaGFT representation into existing machine learning frameworks can enhance up to 40% accuracy of spatial domain identification, cell type annotation, cell-to-spot alignment, and subcellular hallmark inference. SpaGFT identified immunological regions for B cell maturation in human lymph node Visium data, characterized secondary follicle variations from in-house human tonsil CODEX data, and detected extremely rare subcellular organelles such as Cajal body and Set1/COMPASS. This new method lays the groundwork for a new theoretical model in explainable AI, advancing our understanding of tissue organization and function.

T-cell phenotype including CD57+ T follicular helper cells in the tumor microenvironment correlate with a poor outcome in follicular lymphoma

T-lymphocytes are prevalent in the tumor microenvironment of follicular lymphoma (FL). However, the phenotype of T-cells may vary, and the prevalence of certain T-cell subsets may influence tumor biology and patient survival. We therefore analyzed a cohort of 82 FL patients using CyTOF to determine whether specific T-cell phenotypes were associated with distinct tumor microenvironments and patient outcome. We identified four immune subgroups with differing T-cell phenotypes and the prevalence of certain T-cell subsets was associated with patient survival. Patients with increased T cells with early differentiation stage tended to have a significantly better survival than patients with increased T-cells of late differentiation stage. Specifically, CD57+ TFH cells, with a late-stage differentiation phenotype, were significantly more abundant in FL patients who had early disease progression and therefore correlated with an inferior survival. Single cell analysis (CITE-seq) revealed that CD57+ TFH cells exhibited a substantially different transcriptome from CD57- TFH cells with upregulation of inflammatory pathways, evidence of immune exhaustion and susceptibility to apoptosis. Taken together, our results show that different tumor microenvironments among FL patients are associated with variable T-cell phenotypes and an increased prevalence of CD57+ TFH cells is associated with poor patient survival.

Abstract 3487: Molecular clusters and functional drivers of IgM monoclonal gammopathies

IgM MGUS and Waldenstrom Macroglobulinemia (WM) represent a disease spectrum with highly dissimilar therapeutic management ranging from observation to chemoimmunotherapy. The current classification relies solely on clinical features and does not explain the heterogeneity that exists within each of these conditions. To shed light on the biology that may account for the clinical differences, we used bone marrow (BM) CD19+/CD138+ sorted cells and matched BM plasma from 32 patients (pts) (7 IgM MGUS, 25 WM) and 5 healthy controls to perform the first comprehensive multiomics approach including whole exome sequencing, RNA seq, proteomics and metabolomics. Applying principal component analysis to gene expression profiling, most of WM pts clustered together, while a small subset of them grouped separately with MGUS pts, suggesting a biologic dichotomy within WM. The healthy controls formed a distinct group from most WM and MGUS. We then applied a non-negative matrix factorization consensus clustering to the gene expression data and identified three robust clusters. Cluster 1 (C1) included only WM pts, cluster 2 (C2) included both WM and MGUS pts, and cluster 3 (C3) included all controls as well as a small number of WM and MGUS pts. When mutations commonly identified in WM were analyzed, there was no difference among the three groups (excluding controls) in mutation burden of MYD88 L265P and CXCR4. Interestingly, aberrant expression of TNFAIP3 seemed a distinct feature of C1 as deletion of 6q (which encodes for TNFAIP3) and TNFAIP3 mutations were each significantly enriched in C1 (47%) compared to C2 (0%) and C3 (20%; p=0.04). Individual clusters associated with specific transcriptional signature and clinical features. While C1 displayed enrichment of cell growth, downregulation of inflammatory pathways (eg IL6 and IL8 signaling) and aggressive behavior, C2 showed increased inflammatory signaling and cell survival with indolent behavior. C3 had an intermediate feature with combined proliferative and inflammatory signatures. In accordance with the transcriptomic data, the hallmark of C1 was upregulation of proliferation-associated proteins (eg AKT, MAPK) and downregulation of inflammatory proteins (eg IL4, IL10) while the opposite was observed in C2. Once more, C3 confirmed intermediate features with combined upregulation of proliferation and inflammatory proteins. The metabolism was rewired towards fatty acid catabolism in C1, glycolysis in C2 and anabolism in C3. Accordingly, C1 showed undetectable concentration of 3-hydroxybutyric acid as opposed to C2 which had increased levels of malic and lactic acids, as end products of fatty acid oxidation and glycolysis respectively. Those metabolites had intermediate levels in C3. In conclusion, we identified three molecular clusters with distinct clinical, proteomic and metabolomic features, suggesting a potential biologic classification that may have therapeutic implications.

Citation Format: Patrizia Mondello, Jonas Paludo, Joseph Novak, Kerstin Wenzl, Shahrzad Jalali, Jordan Krull, Esteban Braggio, Surendra Dasari, Michelle Manske, Jithma Abeykoon, Asher Chanan-Khan, Robert Kyle, Morie Gertz, Zhi Zhang Yang, Anne Novak, Stephen Ansell. Molecular clusters and functional drivers of IgM monoclonal gammopathies [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 3487.

Abstract 2529: T-cell phenotype and differentiation vary in the tumor microenvironment of follicular lymphoma and are associated with patient outcome

Previous studies have found that the prevalence of subtypes of T cells may be associated with patient outcomes, but a comprehensive analysis of the immune profile in the TM of FL has not been done. In the present study, we identified groups of FL patients with discretely unique TMEs and determined whether the T-cell phenotypes in the TME differed among groups. Using a cohort of 82 FL patients with biopsy specimens collected before treatment, we defined the type of TME based on the content of major lineages (T, B, monocytes/macrophages and NK cells) determined by CyTOF analysis. Hierarchical clustering of this cohort stratified patients into 4 groups with different TME: group 1 (G1) included patients with a high percentage of monocyte/macrophages/NK cells; patients from G2 and G3 were enriched for intratumoral T and B cells, respectively. Patients with intermediate numbers of T and B cells were included in G4. CITRUS analysis revealed that T-cell clusters with phenotypes expressing KLRG1, CD57, PD-1dim and rich in TEMRA cells were significantly more abundant in G1 when compared to G2, G3 or G4. In contrast, T-cell clusters with phenotypes containing CD127, CD45RA, CCR7 or PD-1high cells were significantly less abundant in G1 when compared to other groups. When compared to G3 or G4, 2 classes of T cell clusters, all from CD4+ T cells, were significantly more (CD127+KLRG1-) or less (CD57+PD-1high) abundant in G2, respectively. Clusters with a phenotype rich in short-lived effector cells (SLEC) (from CD8+) were upregulated in G3 when compared to G4. These results suggest that patient groups with distinct TME exhibited variable T-cell phenotypes. To determine the role of T-cell differentiation in predicting patient outcome in FL, we identified T-cell subsets using tSNE plots based on the expression of T-cell maturation markers. We identified 18 subsets (S1-S18) of T cells (11 from CD4, 6 from CD8 and 1 from CD4-CD8-) in each sample. Four subsets - S4 (CD4+TN), S5 (CD4+TCM), S7 (CD8+TN) and S10 (CD8+MEPCs) - were considered to be naïve cells or cells in the early stages of differentiation. Four additional subsets - S8 (CD57+TFH), S9 (CD4+TEXH), S13 (CD8+SLECs) and S17 (CD4+PD-1+Treg) - were memory cells with expression of surface markers indicting late-stage differentiation. We found that the 4 subsets (S4: p=0.08, S5: p=0.01, S7: p=0.04 and S10: p<0.01) with an early-stage phenotype were associated with a favorable clinical outcome. In contrast, subsets (S8: p=0.02, S9: p=0.04, S13: p=0.07 and S17: p=0.06) with a late-stage differentiation phenotype had an unfavorable survival. Supporting this finding, we observed that increased numbers of CD45RA+ T cells correlated with a favorable survival. These results indicate that the differentiation stage may determine the role of T cells in predicting patient outcome in FL.

Citation Format: Zhi-Zhang Yang, Hyo Jin Kim, Hongyan Wu, Xinyi Tang, Jordan Krull, Patrizia Mondello, Jose Villasboas, Anne Novak, Stephen Ansell. T-cell phenotype and differentiation vary in the tumor microenvironment of follicular lymphoma and are associated with patient outcome [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 2529.

Predicted Cellularity using RNASeq-Based Cellular Deconvolution Differentiates Periprosthetic Joint Infection from Non-Infectious Arthroplasty Failure

Differentiation of aseptic and septic causes of arthroplasty failure is essential due to their different treatment schemes. We hypothesized that predicted cellularity profiling of sonicate fluid using the RNASeq-based the cellular deconvolution tools CIBERSORTx and ABIS-seq can differentiate between PJI and non-infected arthroplasty failure (NIAF). Profiles were created for 93 sonicate fluid samples (40 from NIAF and 53 from PJI patients) that had been subjected to RNASeq analysis. CIBERSORTx provided 22 total predicted cell types; 12 predicted cellular abundances were differentially expressed in PJI versus NIAF, including increased neutrophils, activated mast cells, and eosinophils (p≤0.0004) and decreased M0 macrophages, M2 macrophages, and Treg cells (p≤0.0001). ABIS-seq provided 17 total predicted cell types; seven predicted cellular abundances were differentially expressed in PJI versus NIAF, including increased neutrophils, Vd2 γδ T cells, and basophils (p-values of <0.0001, 0.0003, and 0.0263, respectively) and decreased nonclassical/intermediate (NC+I) monocytes, mucosal-associated invariant T (MAIT) cells, and myeloid dendritic cells (mDC) (p-values ≤0.0001). Receiver operative characteristic area under the curve (AUC) analysis identified cells types most predictive of PJI (CIBERSORTx: neutrophils [AUC = 94], activated mast cells [AUC = 93]; ABIS-seq: neutrophils [AUC = 88]) and NIAF (ABIS-seq: NC+I monocytes [AUC = 85], MAIT cells [AUC = 80], mDCs [AUC = 81]). PJI and NIAF samples were differentially clustered by principal component analysis using both tools. Overall, cellularity profiling using RNASeq-based cellular deconvolution can differentiate between PJI and NIAF sonicate fluid samples.

Work supported by NIAMS (NIH R01 AR056647). CF was supported by the Mayo Clinic Graduate School of Biomedical Sciences and the Ph.D. Training Grant in Basic Immunology (NIH R25 GM055252 24).

T-cell phenotype varies in distinct tumor microenvironments and CD57+ TFH cells are associated with disease progression and inferior survival in follicular lymphoma

The tumor microenvironment (TME) plays a crucial role in mediating the tumor immune response, thereby affecting patient outcomes in follicular lymphoma (FL). Using CyTOF, we analyzed a cohort of 82 FL patients with biopsy specimens collected before the first treatment regimen was administered. Our results showed that there was a difference in the T-cell phenotype in patients with various TMEs. Patients with TMEs rich in monocyte/macrophages/NK cells had more CD8+ cells and tended to have fewer TN, but more TEMRA and TEXH cells. Patients enriched with T cells exhibited a T-cell phenotype skewed toward to TCM and patients with B cell dominance tended to have significant higher number of TFH cells. T-cell phenotypes with more terminally differentiated, senescent and exhausted cells were highly represented in patients with less T cells when compared to patients with more T cells. We observed that CD57+TFH cells was significantly more abundant in patients who had disease progression than patients who had complete response to therapy, suggesting a role for CD57+ TFH cells in promoting malignant cell growth in FL. Consistent with this finding, we observed that increased numbers of CD57+ TFH cells correlated with an inferior survival in FL. Using CITE-seq technology, we found that CD57+ TFH cells exhibited a substantially different transcriptome from CD57− TFH cells. Genes that were differentially upregulated in CD57+ TFH cells when compared to CD57− TFH cells included genes involved in cell survival, compromised inflammatory response, and metabolic activation. Taken together, our results indicate different tumor microenvironments among patient groups that is associated with variable T-cell phenotypes.

Intratumoral expression of IL-12 from lentiviral or RNA vectors acts synergistically with TLR4 agonist (GLA) to generate anti-tumor immunological memory

Systemic interleukin-12 (IL12) anti-tumor therapy is highly potent but has had limited utility in the clinic due to severe toxicity. Here, we present two IL12-expressing vector platforms, both of which can overcome the deficiencies of previous systemic IL12 therapies: 1) an integrating lentiviral vector, and 2) a self-replicating messenger RNA formulated with polyethyleneimine. Intratumoral administration of either IL12 vector platform resulted in recruitment of immune cells, including effector T cells and dendritic cells, and the complete remission of established tumors in multiple murine models. Furthermore, concurrent intratumoral administration of the synthetic TLR4 agonist glucopyranosyl lipid A formulated in a stable emulsion (GLA-SE) induced systemic memory T cell responses that mediated complete protection against tumor rechallenge in all survivor mice (8/8 rechallenged mice), whereas only 2/6 total rechallenged mice treated with intratrumoral IL12 monotherapy rejected the rechallenge. Taken together, expression of vectorized IL12 in combination with a TLR4 agonist represents a varied approach to broaden the applicability of intratumoral immune therapies of solid tumors.

Correction to: Application of a multiphase microreactor chemostat for the determination of reaction kinetics of Staphylococcus carnosus

Unfortunately, in the "How to Cite as" section, the given and the family name of the author was incorrectly published, the correct name is Lladó Maldonado. S.

Targeting glycogen synthase kinase 3 for therapeutic benefit in lymphoma

Targeting the B-cell receptor and phosphatidylinositol 3-kinase/mTOR signaling pathways has shown meaningful, but incomplete, antitumor activity in lymphoma. Glycogen synthase kinase 3 (GSK3) α and β are 2 homologous and functionally overlapping serine/threonine kinases that phosphorylate multiple protein substrates in several key signaling pathways. To date, no agent targeting GSK3 has been approved for lymphoma therapy. We show that lymphoma cells abundantly express GSK3α and GSK3β compared with normal B and T lymphocytes at the messenger RNA and protein levels. Utilizing a new GSK3 inhibitor 9-ING-41 and by genetic deletion of GSK3α and GSK3β genes using CRISPR/CAS9 knockout, GSK3 was demonstrated to be functionally important to lymphoma cell growth and proliferation. GSK3β binds to centrosomes and microtubules, and lymphoma cells treated with 9-ING-41 become arrested in mitotic prophase, supporting the notion that GSK3β is necessary for the progression of mitosis. By analyzing recently published RNA sequencing data on 234 diffuse large B-cell lymphoma patients, we found that higher expression of GSK3α or GSK3β correlates well with shorter overall survival. These data provide rationale for testing GSK3 inhibitors in lymphoma patient trials.

Application of a multiphase microreactor chemostat for the determination of reaction kinetics of Staphylococcus carnosus

Bioreactors at the microliter scale offer a promising approach to accelerate bioprocess development. Advantages of such microbioreactors include a reduction in the use of expensive reagents. In this study, a chemostat operation mode of a cuvette-based microbubble column bioreactor made of polystyrene (working volume of 550 µL) was demonstrated. Aeration occurs through a nozzle (Ø ≤ 100 µm) and supports submerged whole-cell cultivation of Staphylococcus carnosus. Stationary concentrations of biomass and glucose were determined in the dilution rate regime ranging from 0.12 to 0.80 1/h with a glucose feed concentration of 1 g/L. For the first time, reaction kinetics of S. carnosus were estimated from data obtained from continuous cultivation. The maximal specific growth rate (µ_max = 0.824 1/h), Monod constant (K_S = 34 × 10^{- 3}g_S/L), substrate-related biomass yield coefficient (Y_X/S = 0.315 g_CDW/g_S), and maintenance coefficient (m_S = 0.0035 g_S/(g_CDW·h)) were determined. These parameters are now available for further studies in the field of synthetic biology.

A compact design for velocity-map imaging of energetic electrons and ions

We present a compact design for a velocity-map imaging spectrometer for energetic electrons and ions. The standard geometry by Eppink and Parker [Rev. Sci. Instrum. 68, 3477 (1997)] is augmented by just two extended electrodes so as to realize an additional einzel lens. In this way, for a maximum electrode voltage of 7 kV, we experimentally demonstrate imaging of electrons with energies up to 65 eV. Simulations show that energy acceptances ≲270 and ≲1200 eV with an energy resolution ΔE∕E≲5% are achievable for electrode voltages ≤20 kV when using diameters of the position-sensitive detector of 42 and 78 mm, respectively.

Application of long synthetic oligonucleotides for gene analysis: effect of probe length and stringency conditions on hybridization specificity

Two different lengths of long unique synthetic oligonucleotide probes (37- and 48-mers) specific for human major histocompatibility complex (MHC) class II beta genes were synthesized. These oligonucleotides were utilized to examine factors influencing hybridization specificity. Both probe length and stringency of washing conditions were found to be crucial factors for sequence-specific hybridization.

Pharmacokinetics and pharmacodynamics of atenolol in children

The pharmacokinetics and pharmacodynamics of intravenous atenolol were studied in 10 children during cardiac electrophysiologic studies. The intravenous pharmacokinetic data were best described by a two-compartment model and revealed the following (mean +/- SD): total body clearance, 0.15 +/- 0.06 L/hr/kg; volume of the central compartment 0.33 +/- 0.06 L/kg; volume of distribution at steady state, 0.83 +/- 0.15 L/kg; distributive elimination half-life, 0.29 +/- 0.08 hour; and terminal elimination half-life, 4.56 +/- 1.05 hours. The data suggest that children have a slightly shorter terminal elimination half-life than that of adults. Pharmacodynamic data showed a significantly (p less than 0.01) increased sinus cycle length and an increase in His to ventricle conduction time (p less than 0.05). Further studies are necessary to determine the optimal oral dose and dosing frequency of atenolol and to access the response of children to long-term treatment.

Application of synthetic oligonucleotides to detect DQ beta genes transmission within insulin-dependent diabetes families

Class II antigen genes encoded by the major histocompatibility complex region (HLA-D region) in man play an important role in susceptibility to insulin dependent diabetes mellitus (IDDM). Evidence suggests that the DQ subregion within the HLA-D region is more directly responsible for susceptibility to IDDM. Therefore, we designed a synthetic oligonucleotide specific for the DQ beta gene to further the understanding of the disease association with HLA-D region genes at the molecular level. Restriction fragment length polymorphism (RFLP) analysis was carried out using DNA isolated from nine families, each including at least two affected siblings (a total of 37 siblings). The segregation pattern of hybridizing fragments showed that: (1) for each of the DR2, DR3, and DR4 specificities, two different alleles can be identified by the DQ beta probe; (2) a 1.9 kb-Taq 1 fragment with the DR4 specificity and a 6.0 kb-Taq-1 fragment within the DR2 specificity tend to cosegregate with IDDM; (3) there was no preferential segregation of the two alleles detected within the DR3 specificity (one allele identified by a 4.7 kb-Taq 1 fragment is quite common among individuals with the DR3 specificity). The results from this study add to the evidence that certain DQ alleles appear to be more directly associated with the diabetogenic gene (or genes) in certain DR specificities.