1
|
Kassis G, Palshikar MG, Hilchey SP, Zand MS, Thakar J. Discrete-state models identify pathway specific B cell states across diseases and infections at single-cell resolution. J Theor Biol 2024; 583:111769. [PMID: 38423206 PMCID: PMC11046450 DOI: 10.1016/j.jtbi.2024.111769] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2023] [Revised: 02/10/2024] [Accepted: 02/17/2024] [Indexed: 03/02/2024]
Abstract
Oxygen (O2) regulated pathways modulate B cell activation, migration and proliferation during infection, vaccination, and other diseases. Modeling these pathways in health and disease is critical to understand B cell states and ways to mediate them. To characterize B cells by their activation of O2 regulated pathways we develop pathway specific discrete state models using previously published single-cell RNA-sequencing (scRNA-seq) datasets from isolated B cells. Specifically, Single Cell Boolean Omics Network Invariant-Time Analysis (scBONITA) was used to infer logic gates for known pathway topologies. The simplest inferred set of logic gates that maximized the number of "OR" interactions between genes was used to simulate B cell networks involved in oxygen sensing until they reached steady network states (attractors). By focusing on the attractors that best represented sequenced cells, we identified genes critical in determining pathway specific cellular states that corresponded to diseased and healthy B cell phenotypes. Specifically, we investigate the transendothelial migration, regulation of actin cytoskeleton, HIF1A, and Citrate Cycle pathways. Our analysis revealed attractors that resembled the state of B cell exhaustion in HIV+ patients as well as attractors that promoted anerobic metabolism, angiogenesis, and tumorigenesis in breast cancer patients, which were eliminated after neoadjuvant chemotherapy (NACT). Finally, we investigated the attractors to which the Azimuth-annotated B cells mapped and found that attractors resembling B cells from HIV+ patients encompassed a significantly larger number of atypical memory B cells than HIV- attractors. Meanwhile, attractors resembling B cells from breast cancer patients post NACT encompassed a reduced number of atypical memory B cells compared to pre-NACT attractors.
Collapse
Affiliation(s)
- George Kassis
- Department of Microbiology and Immunology, University of Rochester School of Medicine and Dentistry, Rochester, USA
| | - Mukta G Palshikar
- Biophysics, Structural, and Computational Biology Program, University of Rochester School of Medicine and Dentistry, Rochester, USA
| | - Shannon P Hilchey
- Department of Medicine, Division of Nephrology, University of Rochester Medical Center, Rochester, NY, USA
| | - Martin S Zand
- Department of Medicine, Division of Nephrology, University of Rochester Medical Center, Rochester, NY, USA
| | - Juilee Thakar
- Department of Microbiology and Immunology, University of Rochester School of Medicine and Dentistry, Rochester, USA; Biophysics, Structural, and Computational Biology Program, University of Rochester School of Medicine and Dentistry, Rochester, USA; Department of Biostatistics and Computational Biology, University of Rochester School of Medicine and Dentistry, Rochester, USA; Department of Biomedical Genetics, University of Rochester School of Medicine and Dentistry, Rochester, USA.
| |
Collapse
|
2
|
Palshikar MG, Min X, Crystal A, Meng J, Hilchey SP, Zand MS, Thakar J. Executable Network Models of Integrated Multiomics Data. J Proteome Res 2023; 22:1546-1556. [PMID: 37000949 PMCID: PMC10167691 DOI: 10.1021/acs.jproteome.2c00730] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/03/2023]
Abstract
Multiomics profiling provides a holistic picture of a condition being examined and captures the complexity of signaling events, beginning from the original cause (environmental or genetic), to downstream functional changes at multiple molecular layers. Pathway enrichment analysis has been used with multiomics data sets to characterize signaling mechanisms. However, technical and biological variability between these layered data limit an integrative computational analyses. We present a Boolean network-based method, multiomics Boolean Omics Network Invariant-Time Analysis (mBONITA), to integrate omics data sets that quantify multiple molecular layers. mBONITA utilizes prior knowledge networks to perform topology-based pathway analysis. In addition, mBONITA identifies genes that are consistently modulated across molecular measurements by combining observed fold-changes and variance, with a measure of node (i.e., gene or protein) influence over signaling, and a measure of the strength of evidence for that gene across data sets. We used mBONITA to integrate multiomics data sets from RAMOS B cells treated with the immunosuppressant drug cyclosporine A under varying O2 tensions to identify pathways involved in hypoxia-mediated chemotaxis. We compare mBONITA's performance with 6 other pathway analysis methods designed for multiomics data and show that mBONITA identifies a set of pathways with evidence of modulation across all omics layers. mBONITA is freely available at https://github.com/Thakar-Lab/mBONITA.
Collapse
|
3
|
Palshikar MG, Palli R, Tyrell A, Maggirwar S, Schifitto G, Singh MV, Thakar J. Executable models of immune signaling pathways in HIV-associated atherosclerosis. NPJ Syst Biol Appl 2022; 8:35. [PMID: 36131068 PMCID: PMC9492768 DOI: 10.1038/s41540-022-00246-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2022] [Accepted: 09/01/2022] [Indexed: 11/09/2022] Open
Abstract
Atherosclerosis (AS)-associated cardiovascular disease is an important cause of mortality in an aging population of people living with HIV (PLWH). This elevated risk has been attributed to viral infection, anti-retroviral therapy, chronic inflammation, and lifestyle factors. However, the rates at which PLWH develop AS vary even after controlling for length of infection, treatment duration, and for lifestyle factors. To investigate the molecular signaling underlying this variation, we sequenced 9368 peripheral blood mononuclear cells (PBMCs) from eight PLWH, four of whom have atherosclerosis (AS+). Additionally, a publicly available dataset of PBMCs from persons before and after HIV infection was used to investigate the effect of acute HIV infection. To characterize dysregulation of pathways rather than just measuring enrichment, we developed the single-cell Boolean Omics Network Invariant Time Analysis (scBONITA) algorithm. scBONITA infers executable dynamic pathway models and performs a perturbation analysis to identify high impact genes. These dynamic models are used for pathway analysis and to map sequenced cells to characteristic signaling states (attractor analysis). scBONITA revealed that lipid signaling regulates cell migration into the vascular endothelium in AS+ PLWH. Pathways implicated included AGE-RAGE and PI3K-AKT signaling in CD8+ T cells, and glucagon and cAMP signaling pathways in monocytes. Attractor analysis with scBONITA facilitated the pathway-based characterization of cellular states in CD8+ T cells and monocytes. In this manner, we identify critical cell-type specific molecular mechanisms underlying HIV-associated atherosclerosis using a novel computational method.
Collapse
Affiliation(s)
- Mukta G Palshikar
- Biophysics, Structural, and Computational Biology Program, University of Rochester School of Medicine and Dentistry, Rochester, USA
| | - Rohith Palli
- Medical Scientist Training Program, University of Rochester School of Medicine and Dentistry, Rochester, USA
| | - Alicia Tyrell
- University of Rochester Clinical & Translational Science Institute, Rochester, USA
| | - Sanjay Maggirwar
- Department of Microbiology, Immunology and Tropical Medicine, George Washington University School of Medicine and Health Sciences, Washington, DC, USA
| | - Giovanni Schifitto
- Department of Neurology, University of Rochester School of Medicine and Dentistry, Rochester, USA
- Department of Imaging Sciences, University of Rochester School of Medicine and Dentistry, Rochester, USA
| | - Meera V Singh
- Department of Neurology, University of Rochester School of Medicine and Dentistry, Rochester, USA
- Department of Microbiology and Immunology, University of Rochester School of Medicine and Dentistry, Rochester, USA
| | - Juilee Thakar
- Biophysics, Structural, and Computational Biology Program, University of Rochester School of Medicine and Dentistry, Rochester, USA.
- Department of Microbiology and Immunology, University of Rochester School of Medicine and Dentistry, Rochester, USA.
- Department of Biostatistics and Computational Biology, University of Rochester School of Medicine and Dentistry, Rochester, USA.
- Department of Biomedical Genetics, University of Rochester School of Medicine and Dentistry, Rochester, USA.
| |
Collapse
|
4
|
Hilchey SP, Palshikar MG, Mendelson ES, Shen S, Rasam S, Emo JA, Qu J, Thakar J, Zand MS. Cyclosporine A Modulates LSP1 Protein Levels in Human B Cells to Attenuate B Cell Migration at Low O 2 Levels. Life (Basel) 2022; 12:life12081284. [PMID: 36013463 PMCID: PMC9410508 DOI: 10.3390/life12081284] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/13/2022] [Revised: 08/08/2022] [Accepted: 08/10/2022] [Indexed: 12/20/2022]
Abstract
Coordinated migration of B cells within and between secondary lymphoid tissues is required for robust antibody responses to infection or vaccination. Secondary lymphoid tissues normally expose B cells to a low O2 (hypoxic) environment. Recently, we have shown that human B cell migration is modulated by an O2-dependent molecular switch, centrally controlled by the hypoxia-induced (transcription) factor-1α (HIF1A), which can be disrupted by the immunosuppressive calcineurin inhibitor, cyclosporine A (CyA). However, the mechanisms by which low O2 environments attenuate B cell migration remain poorly defined. Proteomics analysis has linked CXCR4 chemokine receptor signaling to cytoskeletal rearrangement. We now hypothesize that the pathways linking the O2 sensing molecular switch to chemokine receptor signaling and cytoskeletal rearrangement would likely contain phosphorylation events, which are typically missed in traditional transcriptomic and/or proteomic analyses. Hence, we have performed a comprehensive phosphoproteomics analysis of human B cells treated with CyA after engagement of the chemokine receptor CXCR4 with CXCL12. Statistical analysis of the separate and synergistic effects of CyA and CXCL12 revealed 116 proteins whose abundance is driven by a synergistic interaction between CyA and CXCL12. Further, we used our previously described algorithm BONITA to reveal a critical role for Lymphocyte Specific Protein 1 (LSP1) in cytoskeletal rearrangement. LSP1 is known to modulate neutrophil migration. Validating these modeling results, we show experimentally that LSP1 levels in B cells increase with low O2 exposure, and CyA treatment results in decreased LSP1 protein levels. This correlates with the increased chemotactic activity observed after CyA treatment. Lastly, we directly link LSP1 levels to chemotactic capacity, as shRNA knock-down of LSP1 results in significantly increased B cell chemotaxis at low O2 levels. These results directly link CyA to LSP1-dependent cytoskeletal regulation, demonstrating a previously unrecognized mechanism by which CyA modulates human B cell migration. Data are available via ProteomeXchange with identifier PXD036167.
Collapse
Affiliation(s)
- Shannon P. Hilchey
- Department of Medicine, Division of Nephrology, University of Rochester Medical Center, Rochester, NY 14642, USA
| | - Mukta G. Palshikar
- Biophysics, Structural, and Computational Biology Program, University of Rochester Medical Center, Rochester, NY 14642, USA
| | - Eric S. Mendelson
- Department of Medicine, Division of Nephrology, University of Rochester Medical Center, Rochester, NY 14642, USA
| | - Shichen Shen
- Department of Pharmaceutical Sciences, State University of New York (SUNY) at Buffalo, Buffalo, NY 14203, USA
- New York State Center of Excellence in Bioinformatics & Life Sciences, State University of New York (SUNY) at Buffalo, Buffalo, NY 14203, USA
| | - Sailee Rasam
- Department of Pharmaceutical Sciences, State University of New York (SUNY) at Buffalo, Buffalo, NY 14203, USA
- New York State Center of Excellence in Bioinformatics & Life Sciences, State University of New York (SUNY) at Buffalo, Buffalo, NY 14203, USA
| | - Jason A. Emo
- Department of Medicine, Division of Nephrology, University of Rochester Medical Center, Rochester, NY 14642, USA
| | - Jun Qu
- Department of Pharmaceutical Sciences, State University of New York (SUNY) at Buffalo, Buffalo, NY 14203, USA
- New York State Center of Excellence in Bioinformatics & Life Sciences, State University of New York (SUNY) at Buffalo, Buffalo, NY 14203, USA
| | - Juilee Thakar
- Department of Microbiology and Immunology, University of Rochester Medical Center, Rochester, NY 14642, USA
- Department of Biostatistics and Computational Biology, University of Rochester Medical Center, Rochester, NY 14642, USA
| | - Martin S. Zand
- Department of Medicine, Division of Nephrology, University of Rochester Medical Center, Rochester, NY 14642, USA
- Clinical and Translational Science Institute, University of Rochester Medical Center, Rochester, NY 14642, USA
- Correspondence:
| |
Collapse
|
5
|
Palshikar MG, Hilchey SP, Zand MS, Thakar J. WikiNetworks: translating manually created biological pathways for topological analysis. Bioinformatics 2021; 38:869-871. [PMID: 34636843 PMCID: PMC8756176 DOI: 10.1093/bioinformatics/btab699] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/10/2021] [Revised: 09/28/2021] [Accepted: 10/11/2021] [Indexed: 02/03/2023] Open
Abstract
SUMMARY WikiPathways is a database of 2979 biological pathways across 31 species created using the drawing software PathVisio. Many of these pathways are not directly usable for network-based topological analyses due to differences in curation styles and drawings. We developed the WikiNetworks package to standardize and construct directed networks by combining geometric information and manual annotations from WikiPathways. WikiNetworks performs significantly better than existing tools. This enables the use of high-quality WikiPathways resource for network-based topological analysis of high-throughput data. AVAILABILITY AND IMPLEMENTATION WikiNetworks is written in Python3 and is available on github.com/Thakar-Lab/wikinetworks and on PyPI. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Mukta G Palshikar
- Biophysics, Structural, and Computational Biology Program, University of Rochester School of Medicine and Dentistry, Rochester, NY 14642, USA
| | - Shannon P Hilchey
- Division of Nephrology, University of Rochester School of Medicine and Dentistry, Rochester, NY 14642, USA
| | - Martin S Zand
- Division of Nephrology, University of Rochester School of Medicine and Dentistry, Rochester, NY 14642, USA,Clinical and Translational Science Institute, University of Rochester School of Medicine and Dentistry, Rochester, NY 14642, USA
| | | |
Collapse
|
6
|
Hilchey SP, Palshikar MG, Emo JA, Li D, Garigen J, Wang J, Mendelson ES, Cipolla V, Thakar J, Zand MS. Cyclosporine a directly affects human and mouse b cell migration in vitro by disrupting a hIF-1 αdependent, o 2 sensing, molecular switch. BMC Immunol 2020; 21:13. [PMID: 32183695 PMCID: PMC7079363 DOI: 10.1186/s12865-020-0342-8] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2019] [Accepted: 02/27/2020] [Indexed: 01/24/2023] Open
Abstract
BACKGROUND Hypoxia is a potent molecular signal for cellular metabolism, mitochondrial function, and migration. Conditions of low oxygen tension trigger regulatory cascades mediated via the highly conserved HIF-1 α post-translational modification system. In the adaptive immune response, B cells (Bc) are activated and differentiate under hypoxic conditions within lymph node germinal centers, and subsequently migrate to other compartments. During migration, they traverse through changing oxygen levels, ranging from 1-5% in the lymph node to 5-13% in the peripheral blood. Interestingly, the calcineurin inhibitor cyclosporine A is known to stimulate prolyl hydroxylase activity, resulting in HIF-1 α destabilization and may alter Bc responses directly. Over 60% of patients taking calcineurin immunosuppressant medications have hypo-gammaglobulinemia and poor vaccine responses, putting them at high risk of infection with significantly increased morbidity and mortality. RESULTS We demonstrate that O 2 tension is a previously unrecognized Bc regulatory switch, altering CXCR4 and CXCR5 chemokine receptor signaling in activated Bc through HIF-1 α expression, and controlling critical aspects of Bc migration. Our data demonstrate that calcineurin inhibition hinders this O 2 regulatory switch in primary human Bc. CONCLUSION This previously unrecognized effect of calcineurin inhibition directly on human Bc has significant and direct clinical implications.
Collapse
Affiliation(s)
- Shannon P Hilchey
- University of Rochester Medical CenterDivision of Nephrology, 601 Elmwood Ave., Rochester, 14642 NY USA
| | - Mukta G Palshikar
- University of RochesterBiophysics, Structural, and Computational Biology Program, 601 Elmwood Ave. - Box 675, Rochester, 14642 NY USA
| | - Jason A Emo
- University of Rochester Medical CenterDivision of Nephrology, 601 Elmwood Ave., Rochester, 14642 NY USA
| | - Dongmei Li
- University of RochesterClinical and Translational Science Institute, 265 Crittenden Blvd., Rochester, 14642 NY USA
| | - Jessica Garigen
- University of RochesterClinical and Translational Science Institute, 265 Crittenden Blvd., Rochester, 14642 NY USA
| | - Jiong Wang
- University of Rochester Medical CenterDivision of Nephrology, 601 Elmwood Ave., Rochester, 14642 NY USA
| | - Eric S Mendelson
- University of Rochester Medical CenterDivision of Nephrology, 601 Elmwood Ave., Rochester, 14642 NY USA
| | - Valentina Cipolla
- University of Rochester Medical CenterDivision of Nephrology, 601 Elmwood Ave., Rochester, 14642 NY USA
| | - Juilee Thakar
- University of RochesterDepartment of Microbiology and Immunology, 601 Elmwood Ave - Box 672, Rochester, 14642 NY USA
- University of RochesterDepartment of Biostatistics and Computational Biology, 265 Crittenden Blvd., Rochester, 14642 NY USA
| | - Martin S Zand
- University of Rochester Medical CenterDivision of Nephrology, 601 Elmwood Ave., Rochester, 14642 NY USA
- University of RochesterClinical and Translational Science Institute, 265 Crittenden Blvd., Rochester, 14642 NY USA
| |
Collapse
|