RT Journal Article
SR Electronic
T1 Executable models of pathways built using single-cell RNA seq data reveal immune signaling dysregulations in people living with HIV and atherosclerosis
JF medRxiv
FD Cold Spring Harbor Laboratory Press
SP 2022.03.07.22271522
DO 10.1101/2022.03.07.22271522
A1 Palshikar, Mukta G.
A1 Palli, Rohith
A1 Tyrell, Alicia
A1 Maggirwar, Sanjay
A1 Schifitto, Giovanni
A1 Singh, Meera V.
A1 Thakar, Juilee
YR 2022
UL http://medrxiv.org/content/early/2022/03/07/2022.03.07.22271522.abstract
AB Background Atherosclerosis (AS)-associated cardiovascular disease is an important cause of mortality in an aging population of people living with HIV (PLWH). This elevated risk of atherosclerosis has been attributed to viral infection, prolonged usage of anti-retroviral therapy, and subsequent chronic inflammation.Methods To investigate dysregulated immune signaling in PLWH with and without AS, we sequenced 9368 peripheral blood mononuclear cells (PBMCs) from 8 PLWH, 4 of whom also had atherosclerosis (AS+). To develop executable models of signaling pathways that drive cellular states in HIV-associated atherosclerosis, we developed the single-cell Boolean Omics Network Invariant Time Analysis (scBONITA) algorithm. ScBONITA (a) uses single-cell RNA sequencing data to infer Boolean rules for topologically characterized networks, (b) prioritizes genes based on their impact on signaling, (c) performs pathway analysis, and (d) maps sequenced cells to characteristic signaling states. We used scBONITA to identify dysregulated pathways in different cell-types from AS+ PLWH and AS-PLWH. To compare our findings with pathways associated with HIV infection, we used scBONITA to analyze a publicly available dataset of PBMCs from subjects before and after HIV infection. Additionally, the executable Boolean models characterized by scBONITA were used to analyze observed cellular states corresponding to the steady states of signaling pathwaysResults We identified an increased subpopulation of CD8+ T cells and a decreased subpopulation of monocytes in AS+ PLWH. Dynamic modeling of signaling pathways and pathway analysis with scBONITA provided a new perspective on the mechanisms of HIV-associated atherosclerosis. Lipid metabolism and cell migration pathways are induced by AS rather than by HIV infection. These pathways included AGE-RAGE and PI3K-AKT signaling in CD8+ T cells, and glucagon and cAMP signaling pathways in monocytes. Further analysis of other cell subpopulations suggests that the highly interconnected PI3K-AKT signaling pathway drives cell migratory state in response to dyslipidemia. scBONITA attractor analysis mapped cells to pathway-specific signaling states that correspond to distinct cellular states.Conclusions Dynamic network modeling and pathway analysis with scBONITA indicates that dysregulated lipid signaling regulates cell migration into the vascular endothelium in AS+ PLWH. Attractor analysis with scBONITA facilitated pathway-based characterization of cellular states that are not apparent in gene expression analyses.Competing Interest StatementThe authors have declared no competing interest.Funding StatementThe study was supported by U.S. National Institutes of Health. MGP is supported by R01 AI134058. RP is supported by T32 GM07356. JT was supported by UM1 AI069511, P30 AI078498, R01 AI134058 and R21 AI136668. SBM, GS and AT were supported by R01 HL123346, SBM and MVS are supported by R01 HL128155, R01 NS066801. The University of Rochester Center for AIDS Research (UR-CFAR; P30 AI078498) provided support and core facilities.Author DeclarationsI confirm all relevant ethical guidelines have been followed, and any necessary IRB and/or ethics committee approvals have been obtained.YesThe details of the IRB/oversight body that provided approval or exemption for the research described are given below:All methods were carried out in accordance with University of Rochester guidelines and regulations, and all experimental and study protocols were approved by the University of Rochester Institutional Review Board (#RSRB00063845). Informed consent was obtained from all subjects.I confirm that all necessary patient/participant consent has been obtained and the appropriate institutional forms have been archived, and that any patient/participant/sample identifiers included were not known to anyone (e.g., hospital staff, patients or participants themselves) outside the research group so cannot be used to identify individuals.YesI understand that all clinical trials and any other prospective interventional studies must be registered with an ICMJE-approved registry, such as ClinicalTrials.gov. I confirm that any such study reported in the manuscript has been registered and the trial registration ID is provided (note: if posting a prospective study registered retrospectively, please provide a statement in the trial ID field explaining why the study was not registered in advance).YesI have followed all appropriate research reporting guidelines and uploaded the relevant EQUATOR Network research reporting checklist(s) and other pertinent material as supplementary files, if applicable.YesThe HIV/AS scRNA-seq dataset presented in this manuscript has been deposited in Single Cell Portal and will be made public upon publication of the manuscript (accession number pending). Due to the sensitive nature of HIV data, we have not made the raw data public; however, all results presented in this manuscript may be recapitulated from the data in Single Cell Portal. We also analyzed a previously published dataset that is freely accessible at https://singlecell.broadinstitute.org/single_cell/study/SCP256. All source code and instructions for the scBONITA package is available at github.com/Thakar-Lab/scBONITA.