Organelle Crosstalk Regulators Are Regulated in Diseases, Tumors, and Regulatory T Cells: Novel Classification of Organelle Crosstalk Regulators

To examine whether the expressions of 260 organelle crosstalk regulators (OCRGs) in 16 functional groups are modulated in 23 diseases and 28 tumors, we performed extensive -omics data mining analyses and made a set of significant findings: (1) the ratios of upregulated vs. downregulated OCRGs are 1:2.8 in acute inflammations, 1:1 in metabolic diseases, 1:1.2 in autoimmune diseases, and 1:3.8 in organ failures; (2) sepsis and trauma-upregulated OCRG groups such as vesicle, mitochondrial (MT) fission, and mitophagy but not others, are termed as the cell crisis-handling OCRGs. Similarly, sepsis and trauma plus organ failures upregulated seven OCRG groups including vesicle, MT fission, mitophagy, sarcoplasmic reticulum–MT, MT fusion, autophagosome–lysosome fusion, and autophagosome/endosome–lysosome fusion, classified as the cell failure-handling OCRGs; (3) suppression of autophagosome–lysosome fusion in endothelial and epithelial cells is required for viral replications, which classify this decreased group as the viral replication-suppressed OCRGs; (4) pro-atherogenic damage-associated molecular patterns (DAMPs) such as oxidized low-density lipoprotein (oxLDL), lipopolysaccharide (LPS), oxidized-1-palmitoyl-2-arachidonoyl-sn-glycero-3-phosphocholine (oxPAPC), and interferons (IFNs) totally upregulated 33 OCRGs in endothelial cells (ECs) including vesicle, MT fission, mitophagy, MT fusion, endoplasmic reticulum (ER)–MT contact, ER– plasma membrane (PM) junction, autophagosome/endosome–lysosome fusion, sarcoplasmic reticulum–MT, autophagosome–endosome/lysosome fusion, and ER–Golgi complex (GC) interaction as the 10 EC-activation/inflammation-promoting OCRG groups; (5) the expression of OCRGs is upregulated more than downregulated in regulatory T cells (Tregs) from the lymph nodes, spleen, peripheral blood, intestine, and brown adipose tissue in comparison with that of CD4⁺CD25⁻ T effector controls; (6) toll-like receptors (TLRs), reactive oxygen species (ROS) regulator nuclear factor erythroid 2-related factor 2 (Nrf2), and inflammasome-activated regulator caspase-1 regulated the expressions of OCRGs in diseases, virus-infected cells, and pro-atherogenic DAMP-treated ECs; (7) OCRG expressions are significantly modulated in all the 28 cancer datasets, and the upregulated OCRGs are correlated with tumor immune infiltrates in some tumors; (8) tumor promoter factor IKK2 and tumor suppressor Tp53 significantly modulate the expressions of OCRGs. Our findings provide novel insights on the roles of upregulated OCRGs in the pathogenesis of inflammatory diseases and cancers, and novel pathways for the future therapeutic interventions for inflammations, sepsis, trauma, organ failures, autoimmune diseases, metabolic cardiovascular diseases (CVDs), and cancers.

Introduction

Cardiovascular diseases (CVDs), such as coronary heart disease, hypertension, stroke, and peripheral artery disease, have become the number 1 cause of death globally (1, 2). We and others recently reported that CVD stressors and risk factors such as hyperlipidemia (3, 4), hyperglycemia (5), hyperhomocysteinemia (6, 7), and chronic kidney disease (8–10) promote atherosclerosis and vascular inflammation via several mechanisms. These mechanisms include innate immune activation (11) of endothelial cells (ECs) (3, 12–15) promoting EC injury (16); Ly6Chigh inflammatory mouse monocyte and CD40⁺ human monocyte differentiation (7, 17–19); disease reprogrammed macrophages (20–22); cytokine and secretome regulation (23–30); decreased/transdifferentiated CD4⁺Foxp3⁺ regulatory T cells (Treg) (24, 31–34); and impaired vascular repairability of bone marrow-derived progenitor cells (35, 36). In addition, we recently proposed new models such as intracellular organelle dangers (37) and reactive oxygen species (ROS) as an integrated sensing system for metabolic homeostasis and alarming (38), which indicated that metabolic reprogramming and dysfunction trigger mitochondrial (MT) ROS (4, 39–42); caspase-1/inflammasome activation (8, 10) downregulated histone modification enzymes (43) and increased expressions of trained immunity pathway enzymes (22, 39, 44–47). These reports have clearly demonstrated that mitochondria play significant roles in connecting metabolic reprogramming and dysfunction to inflammation initiation and gene transcription. However, the question of how organelle crosstalk regulates the progression of various diseases remains poorly characterized.

With progress of super-multiplexed optical imaging and barcoding (48) and other system-level spectral imaging and analyses (49, 50), it has been recognized that the crosstalk between mitochondria and other organelles, such as lysosomes, peroxisomes, and the endoplasmic reticulum (ER) (51), is mediated by transcriptional programs and other signaling mechanisms (52). Organelle exchange materials including lipids, ions, and proteins at the membrane contact sites (MCSs), which may be different from membrane trafficking between organelles by vesiculotubular carriers regulated by Rab GTPase (53). Organelle interactions also likely mediate organelle dysfunction downstream of MT impairments (54). In addition, external or internal stress activates several well-orchestrated processes aimed at either restoring cellular homeostasis or committing cell death (55). These processes include unfolded protein response (UPR) in ER and mitochondria (56), autophagy, hypoxia, and MT function, which underlie the ER stress response (57), suggesting that cellular pathological responses may be mediated in organelle crosstalk mechanisms (58). Accordingly, alterations to these networks, such as impaired ER–mitochondria MCSs, have been linked to several diseases such as neurodegeneration (59, 60), CVD (61), diabetes (62), kidney diseases (63, 64), and cancers (65, 66).

Regardless of significant progress in the field, the following questions remained poorly understood: (1) how many regulators (regulatomes) participate in organelle crosstalk; (2) whether the expressions of organelle crosstalk regulatomes are modulated in various disease conditions such as acute inflammation (AI) and injuries, chronic inflammations, autoimmune diseases (ADs), metabolic diseases (MDs), and cancers; (3) whether the expressions of organelle crosstalk regulatomes are differentially modulated in various cell types in response to pathological stimuli; and (4) what are the potential mechanisms regulating expression of organelle crosstalk regulators (OCRGs). To address these important questions, we examined our novel hypothesis that pathological conditions and various DAMPs significantly modulate the expressions of organelle crosstalk regulators in disease-specific and cell type-specific manner. To test this hypothesis, we collected 260 organelle crosstalk regulators as the regulatomes that participate in 16 organelle crosstalk processes including (1) MT biogenesis, (2) MT fission, (3) MT fusion, (4) MT fission and fusion, (5) mitophagy, (6) MT protein translocation, (7) MT contact site, (8) ER (67)–MT contact (68), (9) sarcoplasmic reticulum–MT (69), (10) ER–plasma conjunctions (70, 71), (11) ER–Golgi interaction, (12) ER–endosome, (13) autophagosome–lysosome fusion, (14) autophagosome–endosome/lysosome fusion (68), (15) endosome–Golgi trafficking, and (16) vesicle (Table 1). We performed an extensive -omics database mining and determined the expressions of 260 organelle crosstalk regulators in more than 50 microarrays from more than 20 diseases and 28 cancers/tumors. Corresponding mechanisms were explored, such as the regulation of (a) (TLRs); (b) caspase-1; (c) ROS; and (d) oncogenes/tumor suppressors and made a set of significant findings. Our findings have provided novel insights on pathophysiological roles of OCRGs in various inflammations, diseases, and cancers and provide novel therapeutic targets and strategies for various inflammations, CVD, MDs, and cancers.

TABLE 1

Table 1. Two hundred sixty organelle crosstalk regulators (OCRGs) were searched from three databases such as National Center for Biotechnology Information (NCBI), Gene Set Enrichment Analysis (GSEA), and the Human Protein Atlas (HPA) website.

Methods

Expression Profile of Organelle Crosstalk Regulator Genes in Microarray Data From Patients With Various Inflammatory Diseases and Cancers

The 25 microarray datasets of AI and injuries, ADs, MDs, and CVDs (Table 2); one microarray dataset of Middle East respiratory syndrome coronavirus (MERS-CoV)-infected human microvascular ECs; one microarray dataset of severe respiratory syndrome coronavirus (SARS-CoV)-infected human airway epithelial cells; two microarray datasets of influenza virus-infected lung epithelial cells, Calu-3 cells (non-small-cell lung cancer cell line), and human umbilical vein ECs (HUVECs); one microarray of Kaposi sarcoma-associated herpes virus-infected primary human dermal ECs (Table 3); six microarrays of ECs treated by pro-atherogenic DAMPs (Table 4); six microarray datasets of Treg cells (Table 5); 15 microarrays of Treg regulator deficiency (Table 6); six microarrays of TLR deficiencies (Table 7); one microarray of ROS negative regulator nuclear factor erythroid 2-related factor 2 (Nrf2) deficiency (Table 8); and two microarrays of caspase-1 deficiency (Table 9) were collected from National Institutes of Health (NIH)–National Center for Biotechnology Information (NCBI)–Gene Expression Omnibus (GEO) databases (https://www.ncbi.nlm.nih.gov/gds/) and analyzed with an online software GEO2R (https://www.ncbi.nlm.nih.gov/geo/geo2r/). In addition, gene expression data from 28 cancers were analyzed with the Gene Expression Profiling Interactive Analysis (GEPIA2) database (http://gepia2.cancer-pku.cn/#index), in which dataset sources were from The Cancer Genome Atlas (TCGA)/Genotype-Tissue Expression (GTEx) data (Table 10). Furthermore, 19 microarray datasets of oncogene and tumor suppressor deficiencies were collected from NCBI–GEO and analyzed with GEO2R (Table 11). The differentially expressed OCRGs and their changes in all microarrays and TCGA datasets are listed in Supplementary Tables 1, 2. The original microarray experiments used different cells, which prevented us from comparing the effects of disease conditions in regulating OCRGs in the same cell types. Of note, our approach was well-justified. For example, as a common practice, we (23) and others (72) often studied gene expression in non-ideal heterogeneous peripheral blood mononuclear cell (PBMC) populations in pathophysiological conditions, which are actually composed of many cell types (also see the Discussion section).

TABLE 2

Table 2. Twenty-five microarray datasets including acute inflammations, metabolic diseases, autoimmune diseases, and organ failure diseases in the National Institutes of Health (NIH)–National Center for Biotechnology Information (NCBI)–Gene Expression Omnibus (GEO) dataset database (https://www.ncbi.nlm.nih.gov/gds/) were collected to analyze the expression changes of organelle crosstalk regulators (OCRGs).

TABLE 3

Table 3. Three microarray datasets of time course including MERS coronavirus, SARS coronavirus, avian influenza virus, one dataset of influenza virus infected human umbilical vein endothelial cells and one dataset of Kaposi sarcoma-associated herpes virus infection in the National Institutes of Health (NIH)–National Center for Biotechnology Information (NCBI)–Gene Expression Omnibus (GEO) datasets database (https://www.ncbi.nlm.nih.gov/gds/) were collected to analyze the expression changes of genes that we are interested in.

TABLE 4

Table 4. Six microarray datasets about endothelial cells in pro-atherogenic damage-associated molecular patterns (DAMPs) such as oxLDL, LPS, oxPAPC, and IFN treated conditions in the National Institutes of Health (NIH)–National Center for Biotechnology Information (NCBI)–Gene Expression Omnibus (GEO) datasets database (https://www.ncbi.nlm.nih.gov/gds/) were collected to analyze the changes of organelle interactions and vesicle-related gene.

TABLE 5

Table 5. Six microarray datasets about Treg cells in the National Institutes of Health (NIH)–National Center for Biotechnology Information (NCBI)–Gene Expression Omnibus (GEO) datasets database (https://www.ncbi.nlm.nih.gov/gds/) were collected to analyze the changes of organelle interactions and vesicle-related gene (p < 0.05, |logFC| > 1).

TABLE 6

Table 6. Fifteen microarrays of regulatory T cell regulator deficiency in the National Institutes of Health (NIH)–National Center for Biotechnology Information (NCBI)–Gene Expression Omnibus (GEO) datasets database (https://www.ncbi.nlm.nih.gov/gds/) were collected to analyze the expression changes of OCRGs (p < 0.05, |logFC| > 1).

TABLE 7

Table 7. Six microarrays of toll-like receptor deficiencies in the National Institutes of Health (NIH)–National Center for Biotechnology Information (NCBI)–Gene Expression Omnibus (GEO) datasets database (https://www.ncbi.nlm.nih.gov/gds/) were collected to analyze the expression changes of organelle crosstalk regulators (OCRGs) (p < 0.05, |logFC| > 1).

TABLE 8

Table 8. Reactive oxygen species (ROS) negative regulator Nrf2 deficiency microarray (GSE7810) in the National Institutes of Health (NIH)–National Center for Biotechnology Information (NCBI)–Gene Expression Omnibus (GEO) datasets database (https://www.ncbi.nlm.nih.gov/gds/) was collected to analyze the expression changes of organelle interactions and vesicle-related genes (p < 0.05, |logFC| > 1).

TABLE 9

Table 9. Caspase-1 deficiency microarrays in the National Institutes of Health (NIH)–National Center for Biotechnology Information (NCBI)–Gene Expression Omnibus (GEO) datasets database (https://www.ncbi.nlm.nih.gov/gds/) were collected to analyze the expression changes of organelle crosstalk regulators (OCRGs) (p < 0.05, |logFC| > 1).

TABLE 10

Table 10. Organelle crosstalk regulators (OCRGs) were examined in 28 types of cancers.

TABLE 11

Table 11. Nineteen microarrays of deficiencies of oncogene and tumor suppressor in the National Institutes of Health (NIH)–National Center for Biotechnology Information (NCBI)–Gene Expression Omnibus (GEO) datasets database (https://www.ncbi.nlm.nih.gov/gds/) were collected to analyze the expression changes of organelle crosstalk regulators (OCRGs) (p < 0.05, |logFC| > 1).

Metascape Pathway Analysis

We utilized Metascape Pathway Analysis (MPA; http://metascape.org/gp/index.html#/main/step1) (73) to characterize molecular and cellular functions related to the identified genes in our microarray analysis. Differentially expressed genes were identified and uploaded into MPA for analysis. The core and pathways analysis was used to identify molecular and cellular pathways, as we have previously reported (23).

Protein–Protein Interaction Analysis

Protein–protein interaction (PPI) networks were generated from STRING database (https://string-db.org/). Enrichment analysis results of the shared gene and top 10 connected proteins were downloaded and visualized by using Cytoscape software 3.7.2 (https://cytoscape.org/) (74).

Immune Infiltrate Analyses

TISIDB database (75) (http://cis.hku.hk/TISIDB/index.php) and Gene Set Cancer Analysis (GSCA) database (http://bioinfo.life.hust.edu.cn/GSCA/#/immune) (76) were used in the immune cell infiltrate analyses.

Statistical Analysis of Microarray Data

Twelve housekeeping genes including CHMP2A, EMC7, GPI, PSMB2, PSMB4, RAB7A, SNRPD3, VPS29, VCP, ACTB, RPL27, and OAZ1 (Supplementary Table 3 of Housekeeping Genes) in all 85 GEO datasets regardless of species were chosen for this study. The housekeeping gene list was extracted from the list provided by Eisenberg and de Jonge (77, 78). Briefly, the mean fold change (FC) of housekeeping genes between treatment and control groups varies from 0.72 to 1.55. As this variation was out of the range of FC < 0.5 or FC > 2 (|log2FC| > 1), we concluded that the datasets were of high quality. The target genes with expression changes more than 2-fold in microarrays were defined as the upregulated genes, while genes with their expression decreased more than 2-fold in microarrays were defined as downregulated genes (p < 0.05, |log2FC| > 1).

Results

The Ratios of Upregulated vs. Downregulated Organelle Crosstalk Regulators in Diseases Are Different, and Differentially Expressed Organelle Crosstalk Regulators Are Shared in Diseases

We hypothesized that pathological conditions significantly modulate the expressions of organelle crosstalk regulators in disease-specific and cell type-specific manner. To examine this hypothesis, we collected 260 organelle crosstalk regulators (regulatomic genes and OCRGs) in 16 groups (Figure 1) and were effectively detected (Table 1). Among the OCRGs, 19 genes can be classified into two groups, and one gene can be classified into three groups (Supplementary Figure 1), so that a total of 281 genes were calculated when using a donut chart to show percentages of the gene classifications. The OCRGs included (1) five genes in MT biogenesis (2%), (2) 28 genes in MT fission (10%), (3) 15 genes in MT fusion (5%), (4) seven genes in MT fission and fusion (2%), (5) 28 genes in mitophagy (10%), (6) seven genes in MT protein translocation (2%), (7) four genes in MT contact site (1%), (8) 13 genes in ER–MT contact (5%), (9) nine genes in sarcoplasmic reticulum–MT (3%), (10) 10 genes in ER–plasma conjunctions (3%), (11) six genes in ER–Golgi interaction (2%), (12) seven genes in ER–endosome (2%), (13) four genes in autophagosome–lysosome fusion (1%), (14) nine genes in autophagosome–endosome/lysosome fusion (3%), (15) four genes in endosome–Golgi trafficking (1%), and (16) 125 genes in vesicle (44%). To confirm the functional focuses of the OCRGs, we performed the MPA using the database (http://metascape.org/gp/index.html#/main/step1) (73). Our analysis showed that 260 OCRGs were indeed enriched in MT organization, MT fission, autophagy, vesicle organization, and organelle organization, which can be viewed in Supplementary Figure 2. As shown in Table 2, the expressions of 260 OCRGs were examined in four major categories of diseases, in a total of 23 types of diseases. These diseases and cell types were (a) three types of AIs such as lung injury (whole blood), septic shock (whole blood), and severe trauma (monocytes, leukocytes, and T cells); (b) nine types of MDs including obesity (adipocytes and adipose stem cells) and metabolically healthy obesity (MHO; subcutaneous adipose). Based on the criteria, patients with MHO have no metabolic syndrome (MetS) and insulin resistance (IR) (20) and are metabolically unhealthy obese (MUO; subcutaneous adipose). The difference between MUO and MHO can be found as cited (79), obese with IR (adipocytes), obese with insulin sensitivity (IS) (adipocytes), type 2 diabetes (T2D; liver, skeletal muscle, subcutaneous adipose, and visceral adipose), atherosclerosis (carotid artery plaques, macrophages, and T cells), atherosclerosis and familial combined hyperlipidemia (FCH; plaques and monocytes), and familial hypercholesterolemia (FHC; lymphoblastic cells); (c) seven types of ADs including rheumatoid arthritis (RA; macrophages), acute cutaneous lupus (ACLE; skin), chronic cutaneous lupus (CCLE; skin), psoriasis (skin), subacute cutaneous lupus (SCLE; skin), ulcerative colitis (UC; colon/rectum and PBMCs), and Crohn's disease (CD; PBMC); and (d) four types of organ failures (OFs) such as heart failure (left ventricle), hepatitis B virus-associated acute liver failure (liver), end-stage renal failure (ESRF; whole blood), and chronic kidney disease (CKD) hemodialysis (PBMC). In a total of 39 comparison datasets, we made the following findings: (i) the expressions of OCRGs were modulated except obese with IS; (ii) the expressions of OCRGs in AIs such as septic shock, trauma, ADs, and OFs except heart failure were changed more than those in MDs; (iii) among 38 datasets (except obese IS vs. lean), the downregulations of OCRGs were more than the upregulations of OCRGs in 21 datasets; (iv) the most significant downregulation of OCRGs was found in monocytes, leukocytes, and T cells in severe trauma, skeletal muscle in T2D, and whole blood and PBMC in end-stage renal disease and CKD hemodialysis; (v) the expression changes of OCRGs in T cells in severe trauma were much more than those in T cells in atherosclerosis, suggesting that OCRG expression changes were more disease-dependent rather than cell type-dependent; and (vi) adipocytes had less OCRG expression changes than other cell types.

FIGURE 1

Figure 1. Two hundred sixty organelle crosstalk regulators (OCRGs) in 16 groups were analyzed. There are 16 classification or function of these genes, including nine mitochondria (MT)-related gene populations: MT biogenesis (five genes), MT fusion (15 genes), MT fission (28 genes), MT fission and fusion (seven genes), mitophagy (28 genes), ER–MT contact site (13 genes), MT contact site (four genes), sarcoplasmic reticulum–MT contact site (nine genes), and MT protein translocation (seven genes); three endoplasmic reticulum (ER)-related populations except ER–MT contact site: ER–endosome contact site (seven genes), ER–Golgi interact genes (6), and ER–PM junctions (10 genes); endosome–Golgi trafficking (four genes); autophagosome–lysosome fusion (four genes); autophagosome–endosome/lysosome fusion (nine genes); vesicle-related (125 genes) coded proteins are enhanced by the Human Protein Atlas. Because some genes have more than one function, the final selection is 260 non-repetitive genes. The detailed gene list is in Table 1.

To determine whether among 260 OCRGs upregulated and downregulated regulators were shared in disease categories, we performed a Venn diagram analysis. As shown in Figures 2A–D, in 26 upregulated OCRGs in AIs, the upregulation of two regulators, SERPINA1 and AZU1, was shared among acute respiratory distress syndrome (ARDS), sepsis, and trauma; upregulation of CD24 and FKBP8 was shared among ARDS and trauma; and upregulation of GRN was shared between sepsis and trauma. In 46 upregulated OCRGs in MDs, upregulation of two regulators VAMP8 and SERPINA1 was shared between obese and atherosclerosis; and two upregulated regulators DNM2 and MTFP1 were shared between T2D and atherosclerosis. In 80 upregulated OCRGs in ADs, three regulators (i.e., GPRC5A, LAMP3, and MX2) were shared among RA, autoimmune skin disease (ASD), and inflammatory bowel disease (IBD); three regulators (i.e., OSBPL8, SERPINA1, and MX1) were shared between RA and IBD; seven regulators (i.e., PPARG, MYO19, MAP1LC3A, TLR9, NACC1, P4HA2, and SAMD9) were shared between RA and ASD; and four regulators (i.e., POU2F2, VMP1, ATP11A, and KDR) were shared between ASD and IBD. In 39 upregulated OCRGs in OFs, two regulators (i.e., BAX and GRN) were shared between CKD hemodialysis and hepatitis B virus liver failure. These results have demonstrated that first, the majority of the upregulated OCRGs was disease-specific; second, ADs shared more upregulated OCRGs than other diseases; and third, the ratios of upregulated OCRGs vs. downregulated OCRGs were 1:2.8 in AIs, 1:1 in MDs, 1:1.2 in ADs, and 1:3.8 in OFs, suggesting that AIs and OFs had downregulated OCRGs much more than those of upregulated in comparison with those of others. We then used a Venn diagram to analyze the overlapping OCRGs and their classifications among AIs, MDs, ADs, and OFs (Figure 2E) and listed the exclusively expressed OCRGs upregulated and downregulated in these four types of diseases in Supplementary Tables 4A,B. The results showed that a total of 24 upregulated OCRGs were shared by two or three different diseases. AIs and MDs, and AIs and ADs shared three and one OCRGs, respectively. MDs and ADs, and MDs and OFs shared five and three OCRGs, respectively. ADs and OFs shared six OCRGs. AIs, MDs, and ADs shared three OCRGs. AIs, MDs, and OFs shared one OCRG. MDs, ADs, and OFs shared one OCRG. There were 12, 17, 43, and 15 exclusively upregulated OCRGs in AIs, MDs, ADs, and OF diseases, respectively (Supplementary Table 4A). A total of 67 downregulated OCRGs were shared by two or more than two different diseases (Figure 2F). The four OCRGs including two MT fission genes (i.e., DDHD2 and MAPT) and two vesicle genes (i.e., AASS and CYB5A) were downregulated by four types of diseases. OFs and other diseases (AIs, MDs, and ADs) shared more downregulated OCRGs (12, 11, and 17, respectively). There were 22, 7, 30, and 45 downregulated OCRGs in AIs, MDs, ADs, and OFs, respectively (Supplementary Table 4B). These results showed that differentially classified OCRGs were shared in different diseases.

FIGURE 2

Figure 2. Venn diagram shows that there are several shared OCRGs in acute inflammations (AIs), metabolic diseases (MDs), autoimmune diseases (ADs), and organ failures (OFs). (A) Venn diagram shows 17 genes are shared among acute inflammatory diseases (AIs). SERPINA1 and AZU1 are shared upregulated genes; TOMM20 and ECI2 are shared downregulated genes by sepsis/ARDS, sepsis/shock, and trauma. CD24, FKBP8, and GRN are shared upregulated genes by sepsis and trauma. HSPD1, BNIP3, HSPA9, SET, DDHD2, TOMM70, KLRG1, ZNF266, AFG3L2, and CLIP4 are 10 genes that are commonly downregulated in sepsis and trauma groups. After removing five uncertain genes (upregulated in one type of disease and downregulated in another type of disease), 21 genes are upregulated and 58 genes are downregulated in AIs. (B) Venn diagram shows seven genes are shared among metabolic diseases (MDs). VAMP8 and SERPINA1 are shared upregulated genes by obese and atherosclerosis; DNM2 and MTFP1 are shared upregulated genes by type 2 diabetes (T2D) and atherosclerosis. AASS is shared a downregulated gene by three metabolic diseases, and CD24 and CYB5A are shared downregulated genes by T2D and atherosclerosis. After removal of 12 uncertain genes, 34 genes are upregulated and 34 genes are downregulated in MDs. (C) Venn diagram shows 33 genes are shared among autoimmune diseases. GPRC5A, LAMP3, and MX2 are shared upregulated genes; PPARGC1A and AASS are shared downregulated genes by rheumatoid arthritis (RA), autoimmune skin disease (ASD), and inflammatory bowel disease (IBD). Other 28 genes are common genes between two autoimmune diseases. After removal of 21 uncertain genes, 59 genes are upregulated, and 68 genes are downregulated in ADs. (D) Venn diagram shows 16 genes are shared in organ failure datasets. BAX and GRN are shared upregulated genes by hepatitis B virus-associated acute liver failure (HBV-ALF). AKAP9 is the common downregulated gene in HBV-ALF, end-stage renal failure (ESRF), and hemodialysis; NPC1, RAB7A, CLOCK, KLRG1, SEC23IP, and SNX1 are the common downregulated genes in ESRF and hemodialysis. After removal of 13 uncertain genes, 26 genes are upregulated and 99 genes are downregulated in OFs. (E) Venn diagram shows the shared and exclusive upregulated and downregulated OCRGs in these four diseases. In upregulated OCRGs, a total of 24 genes were shared by two or three different diseases. AIs and MDs shared four genes, two vesicle-related genes (RAB20 and CTSA), and two mitophagy genes (FKBP8 and ATG7). AIs and ADs shared mitophagy gene (MAP1LC3A). AIs, MDs, and OFs shared vesicle genes (GRN). MDs and OFs shared three genes autophagosome–endosome/lysosome fusion regulator (VAMP8), MT fission gene (MTFR2), and vesicle gene (POU2F2). AIs, MD, and ADs shared mitophagy (VMP1), MT fission (MX2), and vesicle (SERPINA1) genes. MDs, ADs, and OFs shared sarcoplasmic reticulum–MT gene (OSBPL8). ADs and OFs shared two mitophagy genes (ATG3 and ATG12), one MT fission gene (DCN), and three vesicle genes (CLTA, ZDHHC13, and SAMD9). MDs and ADs shared five genes including one ER–PM junctions gene (JPH3), one MT fission and fusion regulator (GDAP1), and three vesicle genes (RPTOR, ACBD5, and CCZ1). Other exclusive upregulated genes in these four diseases are listed in Supplementary Table 4A. (F) In downregulated OCRGs, a total of 67 genes were shared by two or more than two different diseases. Two MT fission genes (DDHD2 and MAPT) and two vesicle genes (AASS and CYB5A), were shared by four types of diseases (red box). ADs and OFs shared 17 OCRGs, including one ER–PM junction regulator (JPH2) and one MT contact site gene (APOOL), two MT fission regulators (DNM3 and PINK1) (also mitophagy regulator), two genes with MT fusion and mitophagy function (OMA1 and USP30), and 11 vesicle genes. AIs and OFs shared 12 OCRGs, RAB7A (ER–endosome and autophagosome–endosome/lysosome fusion), JPH3 (ER–PM junctions), VDAC1 (mitophagy, ER–MT contact, and sarcoplasmic reticulum–MT), HSPD1 (MT biogenesis), DHODH (MT fission), BNIP3 (MT fission and fusion, and mitophagy), PID1 (MT fusion), TOMM70 (MT translocation), and four vesicle genes (DNAJA3, KLRG1, SEC23IP, and STX16). MDs and OFs shared 11 OCRGs including ER–endosome genes (NPC1), ER–MT contact and ER–GC interaction regulator VAPA, ER–MT contact, ER–GC interaction and sarcoplasmic reticulum–MT regulator VAPB, two MT fission genes MYO19 and DNM1, and six vesicle genes. AIs and MDs, AIs and ADs, and MDs and ADs shared two, five, and three downregulated OCRGs, respectively. AIs, MDs, and ADs shared three downregulated OCRGs, ITPR1 (ER–MT contact, sarcoplasmic reticulum–MT), OPTN (mitophagy), and DYRK4 (vesicle gene). AIs, MDs, and OFs shared four OCRGs (EGFR, INF2, RAB30, and NOV). AIs, ADs, and OFs shared six OCRGs, VPS41 (autophagosome–endosome/lysosome fusion), two MT fission regulators DDHD1 and TMEM135, and three vesicle regulators (ECI2, CLIP4, and SNX1). Other exclusive downregulated genes in these four diseases are listed in Supplementary Table 4B. Abbreviations: AIs, acute inflammations; MDs, metabolic diseases; ADs, autoimmune diseases; OFs, organ failures; ARDS, acute respiratory distress syndrome; T2D, type 2 diabetes; HF, heart failure; HBV-ALF, hepatitis B virus-associated acute liver failure; ESRF, end-stage renal failure.

Sepsis and Trauma-Upregulated Organelle Crosstalk Regulator Groups Are Classified as the Acute Crisis-Handling Organelle Crosstalk Regulators; and Organ Failure-Upregulated Organelle Crosstalk Regulators Groups Are Classified as the Cell Failure-Handling Organelle Crosstalk Regulators

We then hypothesized that every major disease group modulates differentially 16 OCRG groups. We used the donut chart analysis. As shown in Figure 3, the 21 upregulated OCRGs in AIs were distributed in three categories including vesicle (64%), mitophagy (23%), and MT fission (14%). The 26 upregulated OCRGs in OFs were distributed in seven groups including vesicle (50%), MT fusion (12%), mitophagy (12%), MT fission (8%), sarcoplasmic reticulum–MT (8%), autophagosome–endosome/lysosome fusion (8%), and autophagosome–lysosome fusion (4%). The 34 upregulated OCRGs in MDs were distributed in 12 groups, including vesicle (43%), MT fission (11%), mitophagy (11%), ER–PM junctions (9%), and MT fusion (6%), and in the rest of the seven groups (3% of each). The 59 upregulated OCRGs in ADs were distributed in 14 groups, including vesicle (41%), MT fission (13%), mitophagy (11%), ER–PM junctions (8%), ER–MT contact (5%), and sarcoplasmic reticulum–MT (5%), and in the rest of other groups (2–3% of each). However, the downregulated OCRGs in AIs, OFs, MDs, and ADs were distributed in 15, 15, 9, and 14 groups, respectively. These results have demonstrated that first, vesicle, MT fission, and mitophagy were three top groups of organelle crosstalk regulators upregulated in all four major diseases. Since these three groups of OCRGs are only upregulated in AIs, sepsis, and trauma, we also classify those three groups including vesicle, MT fission, and mitophagy as the cell crisis-handling OCRGs, at least in the partial functions of the three groups, which were well-correlated with a previous report that mitophagy regulator PINK1 is an MT quality control gate keeper (80); second, ER crosstalk regulators were only upregulated in MD and AD groups but not in AIs and OFs; third, a few groups of OCRGs were upregulated in AIs (three groups) and OFs (seven groups), but the other two disease groups had more groups of OCRGs upregulated (MDs had 12 groups and ADs had 14 groups). Therefore, similar to the cell crisis-handling OCRGs that we defined, we also classify vesicle, MT fission, mitophagy, sarcoplasmic reticulum–MT, MT fusion, autophagosome–lysosome fusion, and autophagosome/endosome–lysosome fusion as the cellular failure-handling OCRGs, at least in the partial functions of seven groups; fourth, MT fusion regulators and sarcoplasmic reticulum–MT regulators were upregulated in MDs, ADs, and OFs; and fifth, AIs, MDs, ADs, and OFs downregulated 15, 9, 14, and 15 OCRG groups, respectively, suggesting that AIs and OFs upregulate less OCRG groups but downregulate more OCRG groups than MDs and ADs.

FIGURE 3

Figure 3. Donut chart shows the different classification/function ratios of upregulated and downregulated OCRGs in acute inflammations (AIs), metabolic diseases (MDs), autoimmune diseases (ADs), and organ failures (OFs). From the center to the outside of the donut chart, the circle represents the ratios of OCRGs in acute inflammations (AIs), organ failures (OFs), metabolic diseases (MD), autoimmune diseases (ADs), and all 260 OCRGs in turn. Groups of upregulated OCRGs in AIs (3) and OFs (7) are less than in MDs (12) and OFs (14) (left). In the 21 upregulated genes in AIs, vesicle-related genes account for 64%, and mitophagy and mitochondrial (MT) fission genes account for 14 and 23%, respectively. In the 26 upregulated genes in OFs, the top 3 high-number genes are vesicle- (50%), mitophagy- (12%), and MT fusion-related (12%) OCRGs. In the 34 upregulated genes in MDs, vesicle- (43%), MT fission- (11%), and mitophagy-related (11%) genes are the top 3 genes with high proportion. In the 59 upregulated genes of ADs, vesicle-related genes account for 41%, MT fission genes account for 13%, and mitophagy genes account for 11%. The downregulated OCRGs distributed in AIs, OFs, MDs, and ADs are 15, 15, nine, and 14 groups, respectively (right). In the 58 downregulated genes in AIs, the top 3 high-number genes are vesicle- (28%), MT fission- (14%), and mitophagy-related (13%) genes. In the 99 downregulated genes in OFs, the proportions of vesicle (56%), MT fission (12%), and ER–plasma membrane (PM) junctions (5%) genes are higher than those of others. In the 34 downregulated genes in MDs, the numbers of vesicle (45%), MT fission (13%), and endoplasmic reticulum (ER)–MT contact genes (13%) are higher. In ADs, vesicle (51%), MT fission (12%), and MT fusion genes (8%) are the top 3. These results show that different diseases have different expression patterns of OCRGs and that the classification proportion of upregulated and downregulated genes in the same type of diseases is different, which may explain the different functions of OCRGs in the occurrence and development of diseases. Abbreviations: ER, endoplasmic reticulum; GC, Golgi complex; MT, mitochondria; PM, ER–plasma membrane.

The Majority of Upregulated Pathways Are Disease Group-Specific, and Some Upregulated Pathways Are Shared by Acute Inflammations, Metabolic Diseases, Autoimmune Diseases, and Organ Failures

To determine the functions of upregulated OCRGs, we performed MPAs. As shown in Figures 4A–D, AIs upregulated the top pathways including vacuole organization, lysosomal transport, neutrophil degranulation, regulation of endocytosis, positive regulation of peptidase activity, protein transfer within lipid bilayer, MT fission, vesicle organization, response to acid chemical, lysosomal organization, and receptor-mediated endocytosis. MDs upregulated the top pathways including MT organization, organelle fusion, macroautophagy, vesicle organization, response to anesthesia, myeloid leukocyte activation, protein localization, positive regulation of transmembrane transport, and endosomal transport. ADs upregulated the top pathways including MT fission, protein localization, membrane fusion, organelle fusion, neurodegeneration with brain iron accumulation (NBIA) subtype pathway, insertion of tail-anchored proteins, calcium ion transmembrane transport, vesicle fusion, Fragile X syndrome, MT localization, response to inorganic substance, regulation of NIK/NF-kB signaling, bile acid metabolic process, cholesterol metabolism, lung fibrosis, and protein targeting. OFs upregulated the top pathways including macroautophagy, organelle fusion, autophagosome maturation, cytosolic calcium ion transport, negative regulation of neuron apoptotic process, antigen processing and presentation, and negative regulation of cellular component organization. In Figures 5A–D, AIs downregulated the top 10 pathways including mitochondrion organization, autophagy, MT fission, protein localization, MT fusion, autophagy–animal, selective autophagy, response to nutrient levels, MT calcium ion homeostasis, and frataxin complex. MDs downregulated the top pathways including MT fission, autophagy, Golgi vesicle transport, mitochondrion localization, response to light stimulus, response to cadmium ion, positive regulation of blood circulation, regulation of peptide hormone secretion, L1CAM interactions, and cellular response to unfolded protein. ADs downregulated the top pathways including mitochondrion organization, autophagy, organelle fusion, MT membrane organization, endocytosis, positive regulation of apoptotic process, establishment of protein localization to organelle, negative regulation of MT fission, vesicle organization, and lysosomal transport. OFs downregulated the top pathways including MT fission, vesicle organization, protein localization, organelle fusion, Golgi vesicle transport, regulation of vesicle-mediated transport, MT fusion, negative regulation of mitochondrion organization, lysosomal transport, ROS metabolic process, endocytosis, and response to hyperoxia.

FIGURE 4

Figure 4. Enrichment analysis results show that several pathways of upregulated OCRGs were involved in acute inflammations (AIs), metabolic diseases (MDs), autoimmune diseases (ADs), and organ failures (OFs). Gene Ontology (GO)-based enrichment of significant upregulated OCRGs in different types of diseases was analyzed using Metascape software (http://metascape.org/gp/index.html#/main/step1; PMID: 30944313). (A) There are 14 significant GO enrichment results in 21 upregulated OCRGs in acute inflammations (AIs). The top 3 GO enrichments of OCRGs are vacuole organization, lysosomal transport, and neutrophil degranulation. Except for neutrophil degranulation, relative to acute inflammatory diseases, regulation of endocytosis and receptor-mediated endocytosis are significant GO enrichment results, which show that endocytosis is increased in AIs. (B) There are nine significant GO enrichment results in 34 upregulated OCRGs in metabolic diseases (MDs). The top 3 GO enrichments of OCRGs are mitochondrion organization, organelle fusion, and macroautophagy. Some genes are enriched in myeloid leukocyte activation, which suggested that these genes are involved in chronic low-grade inflammation. (C) There are 20 significant GO enrichment results in 60 upregulated OCRGs in autoimmune diseases (ADs). Mitochondrial fission, protein localization, and macroautophagy are the top 3 GO enrichments. Additionally, calcium ion transmembrane transport and regulation of NIK/NF-kappaB signaling are the significant GO enrichment results. (D) There are seven significant GO enrichment results in 26 upregulated OCRGs in organ failures (OFs). The upregulated genes are enriched in macroautophagy organelle fusion, autophagosome maturation, cytosolic calcium ion transport, negative regulation of neuron apoptotic process, antigen processing and presentation, and negative regulation of cellular component organization. (E) There are shared and exclusive pathways in four types of diseases in upregulated OCRGs. Venn diagram shows that there were five significant pathways shared by acute inflammations (AIs), metabolic diseases (MDs), autoimmune diseases (ADs), and organ failures (OFs). One pathway vesicle organization was shared by AIs and MDs. Mitochondrial fission pathway was shared by AIs and ADs. Macroautophagy pathway was shared by MDs and OFs. Organelle fusion was shared by MDs, ADs, and OFs. (F) The exclusive upregulated pathways in AIs, MDs, ADs, and OFs were 11, five, 17, and five, respectively.

FIGURE 5

Figure 5. Enrichment analysis results show that several pathways of downregulated OCRGs were involved in acute inflammations (AIs), metabolic diseases (MDs), autoimmune diseases (ADs), and organ failures (OFs). Gene Ontology-based enrichment of significant downregulated OCRGs in different types of diseases was analyzed by using Metascape software (http://metascape.org/gp/index.html#/main/step1; PMID: 30944313). (A) In the 20 significant GO enrichment results of 58 downregulated OCRGs in acute inflammations (AIs), the top 3 of downregulated OCRGs are enriched in mitochondrion organization, autophagy, and mitochondrial fission. Mitochondrial calcium ion homeostasis is downregulated in acute inflammation diseases. (B) There are 13 significant GO enrichment results in 34 downregulated OCRGs in metabolic disease (MDs). Mitochondrial fission, autophagy, and Golgi vesicle transport are the top 3 downregulated signaling that enriched downregulated OCRGs. Additionally, lipid transport is suppressed. (C) In the 20 significant GO enrichment results of 67 downregulated OCRGs in autoimmune diseases (ADs), regulation of mitochondrion organization, organelle fusion, and autophagy are the top 3 downregulated enrichment signaling. Cytosolic calcium ion transport is suppressed, while calcium ion transmembrane transport is active in ADs (Figure 5). That is, calcium transporting between membrane-bounded organelles is active. (D) The 20 significant GO enrichment results of 99 downregulated OCRGs in organ failures (OFs). Top 3 GO enrichments of OCRGs are mitochondrial fission, vesicle organization, and protein localization. Reactive oxygen species metabolic process and metabolism of lipids were suppressed in organ failure diseases. (E) There are shared and exclusive pathways in four types of diseases in downregulated OCRGs. Venn diagram shows 13 significant pathways shared by acute inflammations (AIs), metabolic diseases (MDs), autoimmune diseases (ADs), and organ failures (OFs). Mitochondrial fragmentation involved in apoptotic process and mitochondrion organization pathways are downregulated in AIs and ADs. Autophagy signaling is decreased in AIs, MDs, and ADs. Mitochondrial fusion pathway is shared by AIs and OFs. Response to nutrient levels is downregulated in AIs, ADs, and OFs. Four pathways regulated exocytosis, lysosomal transport, vesicle organization, and organelle fusion and are downregulated in ADs and OFs. Response to light stimulus and Golgi vesicle transport are downregulated in MDs and OFs. Downregulation of R-HSA-9609507:Protein localization and mitochondrial fission are shared by AIs, MDs, and OFs. (F) The exclusive downregulated pathways in AIs, MDs, ADs, and OFs are 13, eight, 12 and 10, respectively.

A Venn diagram was used to analyze the overlapping pathways of upregulated and downregulated OCRGs among these four major types of diseases. As shown in Figure 4E, for upregulated OCRGs, vesicle organization was shared by AIs and MDs; MT fission was shared by AIs and ADs; macroautophagy was shared by MDs and OFs; organelle fusion was shared by MDs, ADs, and OFs. Figure 4F lists the exclusively enriched upregulated pathways in these four diseases. In Figure 5E, for downregulated OCRGs, MT fragmentation involved in the apoptotic process and mitochondrion organization was shared by AIs and ADs; autophagy signaling was shared by AIs, MDs, and ADs; MT fusion was shared by AIs and OFs; response to nutrient levels was shared by AIs, ADs, and OFs; regulated exocytosis, lysosomal transport, vesicle organization, and organelle fusion were shared by ADs and OFs; response to light stimulus and Golgi vesicle transport were shared by MDs and OFs; protein localization and MT fission were shared by AIs, MDs, and OFs. Figure 5F lists the exclusively enriched downregulated pathways in these four diseases. These results have demonstrated that (1) the majority of signal pathways for upregulated OCRGs are the major disease group-specific; (2) upregulated pathway organelle fusion is shared by three disease groups such as MDs, ADs, and OFs; (3) upregulated MT fission is shared by AIs and ADs; (4) upregulated vesicle organization is shared by AIs and MDs; (5) upregulated protein localization is shared by MDs and ADs; and (6) more downregulated pathways are shared than upregulated pathways among diseases.

Decreased Autophagosome–Lysosome Fusion Is Required for Viral Replications, Which Classify This Decreased Group as the Viral Replication-Suppressed Organelle Crosstalk Regulators

It has been reported that using single organelle multispectral flow cytometry identified altered energy metabolism, changes in MT size, and MT membrane potential in viral infected cells (81). We hypothesized that organelle crosstalk is modulated in cells infected by viruses. We collected 38 microarray datasets for comparison with seven groups of virus-infected cells including MERS-CoV infected human microvascular ECs (82) (0, 12, 24, 36, and 48 h post infection), SARS-CoV-infected human airway epithelium cells (83) (0, 12, 24, 36, 48, 60, 72, 84, and 96 h post infection), human influenza virus H1N1 infected human airway epithelium cells (83) (0, 6, 12, 18, 24, 36, and 48 h post infection), avian influenza virus A H7N9-infected Calu-3 human lung epithelium cells (84) (3, 7, 12, and 24 h post infection), avian influenza virus H7N7 infected Calu-3 human lung epithelium cells (84) (3, 7, 12, and 24 h post infection), avian influenza virus H5N1-infected Calu-3 human lung epithelium cells (84) (3, 7, 12, and 24 h post infection), human influenza virus H3N2-infected Calu-3 human lung epithelium cells (84) (3, 7, 12, and 24 h post infection), and H9N2-infected HUVECs (24 h). As shown in Table 3, our results indicated that (1) MERS-CoV infection in human microvascular ECs gradually modulated OCRG expressions in 12, 24, 36, and 48 h post infection with the peaks at 24 and 36 h (upregulated 41 and downregulated 44 OCRGs at 24 h and upregulated 39 and downregulated 46 OCRGs, respectively); (2) SARS-CoV infection slightly modulated OCRG expressions in human airway epithelial cells in 0, 24, 36, 48, 60, 72, 84, and 96 h post infection with the peak at 60 h post infection (upregulated 15 OCRGs and downregulated zero OCRG); (3) human influenza virus H1N1 infection significantly modulated OCRG expressions in human airway epithelial cells in 0, 6, 12, 18, 24, 36, and 48 h post infection with the peaks at 18, 24, 36, and 48 h post infection (upregulated 20, 22, 27, and 18 OCRGs and downregulated 21, 29, 30, and 16 OCRGs, respectively); (4) infections of avian influenza virus strains H7N9, H7N7, and H5N1 in human lung Calu-3 epithelium cells and strain H9N2 in HUVECs significantly modulated OCRG expressions at 24 h post infection with upregulation of 30, 38, 40, and 27 OCRGs and downregulation of 70, 109, 92, and 137 OCRGs, respectively; and (5) infection of human influenza virus strains H3N2 in human lung Calu-3 epithelium cells significantly modulated OCRG expressions at 12 and 24 h post infection with upregulation of 12 and 19 OCRGs and downregulation of 16 and 42 OCRGs, respectively.

Infections of avian influenza virus strains at 24 h significantly modulated OCRG expressions. We used a Venn diagram and MPA to examine avian influenza virus modulation of OCRGs. The result showed that 87 OCRGs were upregulated and 190 OCRGs were downregulated in influenza virus strains H1N1-, H7N9-, H7N7-, H5N1-, H3N2-, and H9N2-infected cells. MX1 and SYNPO2 were the commonly upregulated OCRGs shared by these six strains of influenza virus-infected cells; there were no commonly downregulated OCRGs by these six influenza virus-infected cells. Besides the overlapped 47 OCRGs, 40 OCRGs were upregulated (Supplementary Table 5A) and were from 12 functional groups except MT fusion, MT biogenesis, autophagosome–lysosome fusion, and endosome–Golgi complex (GC) groups; 143 OCRGs were downregulated and were from all the 16 functional groups (Figures 6A,C). The significant pathways for upregulated OCRGs included protein localization, mitochondrion organization, membrane fusion, organelle fusion, regulation of calcium ion transmembrane transporter activity, regulation of viral process, endomembrane system organization, metabolism of steroids, regulation of organelle assembly, and leukocyte chemotaxis. The top 10 pathways for downregulated OCRGs included mitochondrion organization, autophagy, MT fusion, MT transport, negative regulation of mitochondrion organization, lysosomal transport, organelle localization, endosomal transport, endocytosis, and protein localization (Supplementary Figures 3A,B).

FIGURE 6

Figure 6. More OCRGs were significantly expressed in influenza virus- and icMERS-inoculated cells. (A) In H1N1-infected human microvascular endothelial cells, H7N9-, H7N7-, H5N1-, and H3N2-infected Calu-3 cells and H9N2-infected human umbilical vein endothelial cells, MX1 and SYNPO2 were the common upregulated OCRGs shared by these six types of flu virus-treated cells; no common downregulated genes were shared by these six flu virus-treated cells. When the overlapped 47 genes were removed, 40 genes are upregulated and 143 genes are downregulated in these five influenza virus-infected cells. (B) In time course, MERS coronavirus-inoculated human microvascular endothelial cells, HSPD1, MX1, STX1A, and GTPBP2 were the common upregulated OCRGs; and OPTN, BNIP3L, VAMP8, VPS41, AASS, EPS15, and PDCD6IP were the common downregulated OCRGs. When the overlapped three genes were removed, 59 genes are upregulated and 62 genes are downregulated in MERS coronavirus-inoculated cells. (C) Donut chart shows the ratio of classification of upregulated OCRGs (A) and downregulated OCRGs (B) in virus-treated cells. Classifications of upregulated OCRGs in icMERS- and influenza virus-inoculated cells are 15 and 12, respectively. Classifications of downregulated OCRGs in icMERS- and influenza virus-inoculated cells are 11 and 16, respectively.

Since MERS-CoV infection in human microvascular ECs induced significant modulation of OCRGs, we further used a Venn diagram and MPA to examine MERS-CoV modulation of OCRGs. As shown in Figure 6B, MERS-CoV infection resulted in upregulation of 62 OCRGs and downregulation of 65 OCRGs. The majority of upregulated OCRGs in 12, 24, 36, and 48 h post infections were shared in at least two time points; four OCRGs were shared in all four time points such as HSPD1, MX1, STX1A, and GTPBP2. Seven downregulated OCRGs such as OPTN, BNIP3L, VAMP8, VPS41, AASS, EPS15, and PDCD6IP were shared in all four time points. Besides the overlapped three genes, the upregulated 59 OCRGs were from 15 functional groups, except autophagosome–lysosome fusion; and the downregulated 62 OCRGs were from 11 functional groups, except ER–PM junctions, sarcoplasmic reticulum–MT, MT biogenesis, ER–GC interaction, and endosome–GC groups (Figure 6C). As shown in Supplementary Table 5B, downregulation of five including VAMP8, VPS41, ATG14, STX17, and RILP out of nine regulators were in autophagosome–endosome/lysosome fusion; and downregulation of one out of four regulators, TIRAP, was in the autophagosome–lysosome fusion group. The top 10 pathways for upregulated OCRGs included mitochondrion organization, synaptic vesicle budding from presynaptic endocytic zone membrane, regulation of mitochondrion organization, autophagy, protein localization, response to unfolded protein, cargo recognition for clathrin-mediated endocytosis, divalent metal ion transport, MT DNA metabolic process, and tissue remodeling. The top 10 pathways for downregulated OCRGs included macroautophagy, organelle fusion, autophagosome maturation, protein localization, vacuolar transport, regulation of MT fission, membrane trafficking, autophagosome membrane docking, endocytosis, and MT transport. The Venn diagram indicated that MT transport, endocytosis, regulated exocytosis, and autophagosome maturation pathways were commonly downregulated in influenza virus- and MERS infectious clone (icMERS)-inoculated cells (Supplementary Figures 3C,D). Taken together, our results have demonstrated that (1) virus infection in ECs and lung epithelial cells significantly modulate OCRGs in all the functional groups; (2) decreased one group autophagosome–lysosome fusion and four signaling pathways including MT transport, endocytosis, regulated exocytosis, and autophagosome maturation are the significant organelle crosstalk features of viral infections. The significant modulation of OCRGs by MERS-CoV in human microvascular ECs may be the underlying mechanism for much higher (41–50%) acute kidney injuries caused by MERS-CoV than that of SARS-CoV (6.7%) and COVID-19 (3%) (85). These results suggest that increased organelle crosstalk in all 15 functional groups (except autophagosome–lysosome fusion) but decreased lysosome degradation are required for viral replication; and significant modulation of organelle crosstalk in human microvascular ECs in coronavirus family infection represented by MERS-CoV may be the important underlying mechanism for COVID-19 (caused by SARS-CoV 2)-induced cardiovascular complications (86, 87).

Organelle Crosstalk Regulators Upregulated by Pro-atherogenic Damage-Associated Molecular Patterns in Endothelial Cells Are Classified Endothelial Cell-Activation/Inflammation-Promoting Organelle Crosstalk Regulator Groups

Our previous reports showed that sterile inflammatory stimuli can cause intracellular organelle stress that occurs not only in cancer cells but also in vascular ECs (37). We reported that lysolipids are capable of transdifferentiating human aortic ECs (HAECs) into innate immune cells, including induction of potent DAMP receptors, such as CD36 molecule (14). CD36 pathways are activated by several distinct ligands, which converged on these pathways and results in inflammatory responses and endothelial dysfunction. CD36 pathway may be an underlying cause of CVDs and cerebrovascular diseases (88). Oxidized lipids upregulated by CD36 can modulate endothelial properties and may contribute to atherogenesis (89). However, how DAMP receptor signaling modulates expression of organelle crosstalking in ECs remains poorly understood. We examined the expression changes of OCRGs in ECs under stimulations of various pathogen-associated molecular patterns (PAMPs)/DAMPs (90–92). As shown in Table 3, influenza virus infection in HUVECs upregulated 27 OCRGs and downregulated 137 OCRGs. Kaposi sarcoma-associated herpes virus infection in human dermal ECs upregulated six OCRGs and downregulated 10 OCRGs, suggesting that acute inflammatory virus infection has higher OCRG expression modulation than chronic tumorigenic virus infection. Lipopolysaccharide (LPS) stimulation of human lung microvascular ECs (HLMECs) for 4, 8, and 24 h upregulated 8, 17, and 9 OCRGs, respectively, and downregulated 7, 15, and 14 OCRGs, respectively. LPS stimulation of mouse aortic ECs (MAECs) for 4 h upregulated three OCRGs and downregulated five OCRGs, which were similar to those found in HLMECs. It has been reported that Notch1 signaling is involved in regulating MT biogenesis (93). On the other hand, organelle crosstalking also regulates Notch1 signaling. For example, the endocytic trafficking of Notch receptor leads to either transportation to lysosomes for degradation via multi-vesicular bodies (MVBs) and late endosomes or recycling back to the plasma membrane for ligand binding activation (94). However, how modulation of Notch signaling regulates organelle crosstalking in ECs remains poorly characterized. Inhibition of endothelial NOTCH1 (inflammation promoting) signaling with siRNA in the absence or presence of inflammatory cytokine interleukin-1β (IL-1β) stimulation for 24 h in HUVECs resulted in no modulation or one OCRG downregulation (with IL-1β) (95). However, another report showed that reduction of NOTCH1 expression in HAECs by siRNA, in the absence of stimulations with inflammatory lipids or cytokines, increased inflammatory molecules and binding of monocytes (96). In the datasets of this paper, we found that inhibition of NOTCH1 expression with siRNA in HAECs upregulated three OCRGs and downregulated six OCRGs. In addition, interferon-α (IFNα), IFNβ and IFNγ treatment of HUVECs led to upregulation of 7, 11, and 6 OCRGs, respectively, but no downregulation of OCRGs. Of note, oscillatory shear vs. laminar shear (fibrosa), oscillatory shear vs. laminar shear (ventricularis) (97), and MAECs from atherogenic apolipoprotein E (ApoE)-deficient mice vs. MAECs from wild-type control mice had no modulation of OCRG expressions potentially due to a long-term chronic adaptation process in response to shear stress or hyperlipidemia. Pro-atherogenic DAMP stimuli oxidized low-density lipoprotein (oxLDL) (98) and proinflammatory oxidized 1-palmitoyl-2-arachidonoyl-sn-glycerol-3-phosphatidylcholine (oxPAPC) stimulation of MAECs presumably via transient receptor potential ankyrin 1 (TRPA1) (99) upregulated three, zero, and seven OCRGs and downregulated five, six, and seven OCRGs (Table 4).

As shown in Figures 7A,B (Supplementary Table 6 showed the detailed differentially expressed OCRGs), the 33 upregulated OCRGs in ECs included 10 out of 16 functional groups of OCRGs (except MT fission and fusion, MT translocation, MT biogenesis, MT contact site, autophagosome–lysosome fusion, and endosome–GC groups): (1) six MT fission regulators, DNM2, INF2, TMEM135, COX10, DDHD1, and LRRK2; (2) two mitophagy regulators, SQSTM1 (also vesicle group) and OPTN; (3) two ER–MT contact and sarcoplasmic reticulum–MT regulators, HSPA9 and ITPR1; (4) one ER–PM junction regulator, JPH2; (5) one MT fusion regulator, BAK1 (also ER–MT contact group); (6) one autophagosome–endosome/lysosome fusion regulator, STX17; and (7) the rest of the 17 vesicle regulators including BMP2, CD24, CSF2, POU2F2, ATP23, CLTC, P4HA2, VPS26A, ATP11A, CCZ1, EIF4ENIF1, RPS6KC1, GTPBP2, NOV, PICALM, PNPLA2, and RUNX1. The 36 downregulated OCRGs in ECs also included 10 out of 16 functional groups of OCRGs (except MT fusion, sarcoplasmic reticulum–MT, autophagosome–endosome/lysosome fusion, MT biogenesis, MT contact site, and endosome–GC groups). The upregulated top 5 pathways involved in DAMP-stimulated ECs included MT fission, autophagy, membrane trafficking, Golgi organization, and WNT pathway (Figure 7C). The downregulated top 5 pathways included MT fission, autophagy, synthesis of bile acids and bile salts, endocytosis, and RAB geranylgeranylation (Figure 7D). By comparison, 20 OCRGs were upregulated and were classified into five groups (Supplementary Figures 5A,B); the top pathways involved in upregulating OCRGs in LPS-stimulated ECs included MT fission, endosomal transport, receptor internalization, cytokine signaling in immune system, and ZNF410 TARGET GENES (Supplementary Figure 5C), latter four top pathways of which were different from those in DAMP-stimulated upregulation of OCRGs in ECs. The top pathways involved in downregulating OCRGs in LPS-stimulated ECs included MT fission, regulation of calcium ion transmembrane transporter activity, negative regulation of transferase activity, receptor metabolic process, and establishment of organelle localization, also latter four top pathways of which were different from those in DAMP-stimulated upregulation of OCRGs in ECs (Supplementary Figure 5D). Taken together, our results have demonstrated that (1) acute inflammatory (influenza) virus infection has higher OCRG modulation than in chronic tumorigenic (herpes) virus infection in ECs; (2) EC activation stimuli such as TLR4 stimulation by LPS and oxLDL result in modulation of OCRGs in HAECs and MAECs; (3) modulation of OCRGs by inhibition of NOTCH1 in HAEC is associated with increased inflammatory status and increased atherosclerosis in endothelial NOTCH1 heterozygous mice; (4) proinflammatory cytokines type 1 IFN (α, β, and probably ω, κ, and ε) signaling and type II IFN (γ) (100) modulated OCRG expression in ECs; and (5) oscillatory shear vs. laminar shear vs. laminar shear and ApoE KO MAECs have no modulation of OCRGs, potentially due to adaptation of ApoE KO MAECs and interplays of multiple pathways bone morphogenetic protein (BMP)-transforming growth factor-β (TGFβ), Wingless and Int-1 (WNT), Notch membrane receptors (NOTCH), hypoxia inducible factor-α (HIF1α a transcription factor), Twist family basic helix-loop-helix (BHLH) transcription factor 1 (TWIST1), and a subset of homeobox genes (HOX family) (101).

FIGURE 7

Figure 7. Classifications of OCRGs differentially expressed in endothelial cells treated with pro-atherogenic DAMPs such as oxLDL, LPS, oxPAPC, and IFNs treated. (A) When the overlapped genes in all upregulated OCRGs and in all downregulated OCRGs were removed, 33 genes are upregulated and 36 genes are downregulated exclusively in inflammatory factor-treated endothelial cells. (B) Donut chart shows the ratio of classification of upregulated OCRGs (left) and downregulated OCRGs (right) in pro-atherogenic DAMP-treated endothelial cells. There are 10 classifications in 33 upregulated and 36 downregulated OCRGs, respectively. There are no MT fission and fusion, MT translocation, MT biogenesis, MT contact site, autophagosome–lysosome fusion, and endosome–GC OCRGs in upregulated genes and not MT fusion, sarcoplasmic reticulum–MT, autophagosome–endosome/lysosome fusion, MT biogenesis, MT contact site, and endosome–GC OCRGs in downregulated genes. (C) There are 20 significant GO enrichment results in 33 upregulated OCRGs of pro-atherogenic DAMP-treated endothelial cells. Top 3 GO enrichments of OCRGs are membrane trafficking, autophagy, and mitochondrial fission. Wnt signaling pathway is upregulated, and this change suggests inflammatory factors may be via Wnt signaling pathway, promoting inflammation induced by OCRGs. (D) There are 14 significant GO enrichment results in 36 downregulated OCRGs in inflammatory factor-treated endothelial cells. Calcium ion transmembrane transporter activity signaling is downregulated, and this change suggests that inflammatory factors can change calcium ion transmembrane transporter activity via OCRGs in endothelial cells. Enrichment analysis was performed by using Metascape software (http://metascape.org/gp/index.html#/main/step1; PMID: 30944313).

The Expressions of Organelle Crosstalk Regulators Are Modulated in Regulatory T Cells in Comparison With Those of CD4⁺CD25^– T Effector Controls; Upregulated Organelle Crosstalk Regulators Are More Than Downregulated Organelle Crosstalk Regulators in Regulatory T Cells

Recent reports showed that the MT membrane potential and the proliferation inhibition function of CD4⁺Foxp3⁺ Treg from patients with myasthenia gravis are enhanced by autophagy inducing agent rapamycin and suppressed by phosphoinositide 3-kinase (PI3 kinase) and autophagy inhibitor 3-methyladenine (3-MA) (102) and that during autoimmune conditions, Treg function alterations associate with MT oxidative stress, dysfunctional mitophagy, and enhanced DNA damage and cell death (103, 104). We hypothesized that the expressions of OCRGs are modulated in CD4⁺CD25^highFOXP3⁺ Treg in comparison with those of CD4⁺CD25⁻ T effector cell controls. As shown in Table 5, Treg from eight different tissues including lymph nodes (LNs), visceral adipose tissue, brown adipose tissue (cold), brown adipose tissue (warm), spleen from mice with skeletal muscle injury for 4 days, spleen from mice with skeletal muscle injury for 14 days, and spleen from male and female mice had no downregulated OCRGs. In addition, Treg from five tissues such as warm spleen, skeletal muscle from mice with skeletal muscle injury for 4 days, skeletal muscle from mice with skeletal muscle injury for 14 days, intestine, and peripheral blood had one, four, two, one, and one downregulated OCRGs, respectively. Moreover, Treg in spleen from female mice had no upregulated OCRGs; and Treg from the rest of the 12 groups had one to 12 upregulated OCRGs (Supplementary Table 7 listed, in detail, the significant expressed OCRGs in Treg compared with conventional T cells). To examine a hypothesis that the expressions of OCRGs are required for the suppressive functions of Treg, we determined this hypothesis with Treg signature gene deficient microarray datasets. The 15 microarrays of Treg regulator deficiency datasets were collected to analyze the expression changes of OCRGs (Table 6). In the deficiency datasets of Treg signature genes GATA binding protein 3 (Gata3), histone deacetylase 9 (Hdac9), and peroxisome proliferator activated receptor gamma (Pparg) in LNs, there were no significant expression changes of OCRGs. There were seven significantly upregulated OCRGs and no downregulated OCRGs in Treg signature transcription factor B-cell lymphoma 6 (Bcl6) knockout Treg from spleen and LN. There were six significantly upregulated OCRGs and one downregulated OCRG in Foxo1 knockout Treg from the thymus, spleen, and LNs. In Dicer 1 (ribonuclease III, microRNA-maturation enzyme) knockout CD4⁺ T cells, there were five upregulated OCRGs and one downregulated OCRG. In tripartite motif containing 28 [E3 small ubiquitin-like modifier (SUMO)-protein ligase Trim28] knockout Treg, there were three upregulated OCRGs and five downregulated OCRGs. In B lymphocyte-induced maturation protein 1 (Blimp1) deficiency Treg, there were three upregulated OCRGs and six downregulated OCRGs.

As shown in Figure 8A, when combining all the upregulated and downregulated OCRGs in the Treg, besides the overlapping genes, four OCRGs downregulated in Treg included sarcoplasmic reticulum–mitochondria contact regulator OSBP-related protein 8 (Osbpl8) (105), MT fission regulators MX1 and MX2, and one vesicle regulator Hexim1 (Supplementary Table 8). There were no significant pathways of the upregulated OCRGs from MPA. The 19 upregulated OCRGs included one autophagosome–lysosome fusion regulator Tirap, one mitophagy regulator Vmp1, one MT fission gene Dcn, one MT translocation Tomm7, Vdac1 with more function (mitophagy, ER–MT contact, sarcoplasmic reticulum–MT), MT fission and fusion, mitophagy regulator Bnip3, and the other 13 vesicle regulators (Supplementary Table 8 and Figure 8B). The top signaling pathways of these upregulated OCRGs from the MPA included positive regulation of macroautophagy, regeneration, regulation of protein stability, positive regulation of apoptotic process, and leukocyte chemotaxis (Figure 8C). These results have demonstrated for the first time that Treg from various tissues and mice with injured skeletal muscle have upregulation of OCRGs compared with those of CD4⁺CD25⁻ T effector controls; and Treg slightly upregulate a few OCRGs in Treg from visceral adipose tissue, skeletal muscle, and intestine. The results of weakened Treg cells from 11 microarrays (Treg signature gene knockout datasets) showed that the upregulated OCRGs were more than downregulated OCRGs (Table 6 and Supplementary Table 9 list in detail differentially expressed OCRGs). As shown in Figure 8D, the 17 upregulated OCRGs included one autophagosome–lysosome fusion gene Tirap, one ER–endosome contact gene Npc1, two ER–MT contact genes Bcap31 and Itpr1 (also sarcoplasmic reticulum–MT group), one ER–PM junctions Stx1a, two mitophagy genes Optn and Atg7, two MT fission genes Mtfr2 and Mtfp1, and the other eight vesicle regulators (Nsdhl, Clta, Abcd3, Ctsa, Rab11fip5, Dyrk4, Mbd1, and Tyr). The eight downregulated OCRGs included two MT fission regulators Dcn and Mapt and six vesicle regulators (Ankrd6, Igf2r, Gtpbp2, Dpp7, Gprc5a, and Grn) (Supplementary Table 9 and Figure 8E). There were no significant pathways in the downregulated OCRGs from MPA. The signaling pathways of these upregulated OCRGs included regulation of autophagy, adaptive immune system, Fragile X syndrome, insulin secretion lipid transport, and organic hydroxy compound metabolic process (Figure 8F). Taken together, our results have demonstrated that (1) the expression of OCRGs are modulated in CD4⁺FOXP3⁺ Treg and Treg cells with signature genes deficiencies; and (2) upregulated OCRGs in Treg are more than downregulated OCRGs. In weakened Treg when Treg signature genes are deficient, upregulated OCRGs are more than downregulated OCRGs; and (3) positive regulation of macroautophagy, regeneration, regulation of protein stability, positive regulation of apoptotic process, and leukocyte chemotaxis signaling are upregulated in Treg; and regulation of autophagy, adaptive immune system, Fragile X syndrome, insulin secretion, lipid transport, and organic hydroxy compound metabolic process are upregulated in weakened Treg.

FIGURE 8

Figure 8. Less classification of OCRGs differentially expressed in regulatory T cells and weakened Treg cells. (A) When the overlapped three genes in all upregulated OCRGs and in all downregulated OCRGs were removed, 19 genes are upregulated and four genes are downregulated exclusively in regulatory T cells. (B) Donut chart shows the ratio of classification of upregulated OCRGs (left) and downregulated OCRGs (right) in regulatory T cells. There are eight classifications in 19 upregulated OCRGs and three classifications in four downregulated OCRGs. (C) Enrichment analysis results of 19 upregulated genes show signaling pathways of positive regulation of macroautophagy, regeneration, regulation of protein stability, positive regulation of apoptotic process, leukocyte chemotaxis, response to wounding, metabolism of lipids, transmembrane receptor protein tyrosine kinase signaling pathway, and negative regulation of transport were upregulated. (D) When the overlapped five genes in all upregulated genes and in all downregulated genes were removed, 17 genes are upregulated and eight genes are downregulated exclusively in regulatory T-cell signature gene deficiency (weakened regulatory T cells). (E) Donut chart shows the ratio of classification of upregulated OCRGs (left) and downregulated OCRGs (right) in weakened regulatory T cells. There are eight classifications (vesicle, MT fission, Mitophagy, ER–MT contact, ER–PM junctions, sarcoplasmic reticulum–MT, ER–endosome, and autophagosome–lysosome fusion) in 17 upregulated OCRGs and two classifications (vesicle and MT fission) in eight downregulated OCRGs, respectively. (F) Enrichment analysis results of 17 upregulated genes show signaling pathways of regulation of autophagy, adaptive immune system, Fragile X syndrome, insulin secretion, lipid transport, and organic hydroxy compound metabolic process.

Toll-Like Receptors, Reactive Oxygen Species Regulator Nuclear Factor Erythroid 2-Related Factor 2, and Inflammasome-Activated Caspase-1 Regulate the Expressions of Organelle Crosstalk Regulators

In order to determine the mechanisms underlying the expression changes of OCRGs in diseases, virus-infected cells, and pro-atherogenic DAMP-treated ECs, we examined expressions of OCRGs in nine microarrays including six TLRs (106), one ROS regulator Nrf2 (23), and two caspase-1-deficient microarrays. Caspase-1 is activated in the DAMP/PAMP-sensing protein complexes termed inflammasomes (3, 11, 16, 36, 39, 90, 106–110) (Tables 7–9). The results showed that TLR2, TLR4, and TLR3/7/9 deficiencies upregulated 29, seven, and nine OCRGs and downregulated 37, 12, and 12 OCRGs, respectively (Supplementary Figure 6). In addition, ROS negative regulator Nrf2 deficiency upregulated 14 and downregulated seven OCRGs. Moreover, caspase-1 deficiency upregulated seven and downregulated five OCRGs. Then we proposed two models to explore the mechanisms underlying the expression changes of OCRGs regulated by DAMP/PAMP sensing TLRs and caspase-1/inflammasomes (Figures 9A,B). A Venn diagram was used to identify the regulated OCRGs in diseases, virus infections, and cells stimulated by pro-atherogenic DAMPs (Supplementary Figures 7–11). As shown in Figure 9C and Supplementary Table 10, TLRs, caspase-1, and ROS regulators play significant roles in upregulating a total of 46 OCRGs (17.7%) from 14 groups (except for MT fission and fusion and ER–MT contact OCRGs) in different diseases, virus, and pro-atherogenic factor-treated cells. Of note, the signaling pathways in upregulating 12 OCRGs (i.e., CSF2, MX1, ABCD3, ATP2A1, APOOL, PDXDC1, MX2, CSF2, POU2F2, MCOLN1, NPC1, RUNX1, and BMP2) are shared by two groups of upstream regulators. These findings suggest that these 12 OCRGs may play important roles in pathophysiologies upregulated by DAMP/PAMP signaling.

FIGURE 9

Figure 9. Three types of regulators (TLRs, caspase-1, and ROS regulators) upregulated a total of 46 OCRGs in different diseases, and virus- and pro-atherogenic factor-treated cells. (A) Two models were used to explore the mechanisms of expression changes of OCRGs in the upstream gene-deficient microarrays. When the upstream gene A plays a positive role (promote disease) in disease (or in virus, pro-atherogenic DAMP-treated cells), OCRGs upregulated in the disease will be downregulated in deficiency of upstream gene A (upper left box). OCRGs downregulated in the disease will be upregulated in deficiency of upstream gene A (bottom right box). (B) When the upstream gene B plays a negative role (suppress disease) in disease (or in virus, pro-atherogenic DAMP-treated cells), OCRGs upregulated in the disease will be upregulated in deficiency of upstream gene B (upper right box). OCRGs downregulated in the disease will be downregulated in deficiency of upstream gene B (bottom left box). (C) Alluvial plot shows the regulated 46 OCRGs by TLR2, TLR4, TLR3/7/9, caspase-1, and ROS regulator Nrf2; and these 46 OCRGs classified 14 groups of all 16 groups (except for MT fission and fusion, and ER–MT contact OCRGs), and vesicle, MT fission, and mitophagy are the top 3 groups. TLR2 regulated 24 OCRGs and caspase-1 regulated five OCRGs. Twelve OCRGs (i.e., CSF2, MX1, ABCD3, ATP2A1, APOOL, PDXDC1, MX2, CSF2, POU2F2, MCOLN1, NPC1, RUNX1, and BMP2) were regulated by two upstream regulators; other OCRGs were regulated by one upstream regulator. The regulated OCRGs were more in autoimmune diseases and virus-infected cells than in other diseases. Alluvial plot was plotted by http://www.bioinformatics.com.cn, an online platform for data analysis and visualization. The detailed relation between OCRGs and upstream regulators is listed in Supplementary Table 10.

Organelle Crosstalk Regulators Are Significantly Modulated in Cancer; Upregulations of Organelle Crosstalk Regulators in Liver Hepatocellular Carcinoma and Lung Squamous Cell Carcinoma Are Correlated With Immune Infiltrates; Mitochondrial Import Machine Is Upregulated in Cancers and Is Associated With Immune Infiltrates

It has been reported that mitophagy pathways are intricately linked to the metabolic rewiring of cancer cells to support the high bioenergetic demand of the tumors (111); autophagy, acting as a cancer-suppressive function, is inclined to hinder metastasis by selectively downregulating critical transcription factors of the epithelial–mesenchymal transition (EMT) in the early phases (112). The molecular mechanisms underlying fusion, with either lysosomes or plasma membrane, are key determinants to maintain cell homeostasis upon stressing stimuli. The accumulation of undigested substrates leads to cancer and other diseases such as lysosomal storage disorders and age-related neurodegenerative diseases (113). We hypothesized that carcinogenesis modulates the expressions of OCRGs. To examine this hypothesis, we collected the expression data of 260 OCRGs in 28 cancer datasets in the NIH-NCI TCGA (https://www.cancer.gov/about-nci/organization/ccg/research/structural-genomics/tcga/using-tcga)/GTEx database (GTEx Portal, https://www.gtexportal.org/home/) from EGPIA2 (http://gepia2.cancer-pku.cn/#index). As shown in Table 10, all the 28 cancers from the gastrointestinal, respiratory, brain, genitourinary, hematopoietic, and digestive systems significantly modulated the expressions of OCRGs. The upregulated OCRGs range from one OCRG in lung adenocarcinoma (LUAD) to 166 OCRGs in thymoma; and downregulated OCRGs range from three in pancreatic adenocarcinoma (PAAD) to 65 in testicular germ cell tumor.

To identify the common features of OCRG modulation in cancers, we calculated the percentage (a/b%) of upregulated OCRGs (a) over total upregulated expressed genes and (b) in the cancer microarrays, and the ratios of downregulated OCRGs over total downregulated genes. We noticed that 16 cancers had upregulated OCRG percentage >1% and 12 cancers had upregulated OCRGs percentages <1%. Cancers were roughly classified into eight systems to explore the expression changes of OCRGs. The numbers of upregulated OCRGs were more than those of downregulated genes in tumors of the immune system, digestive system, and nervous system. The numbers of downregulated OCRGs were more than those of the upregulated OCRGs in tumors of the endocrine system, respiratory system, and circulation system and most tumors of the reproductive system and urinary system. In breast invasive carcinoma (BRCA) and skin cutaneous melanoma (SKCM), the numbers of upregulated OCRGs were more than those of downregulated OCRGs. In head and neck squamous cell carcinoma (HNSC), the numbers of upregulated OCRGs were less than those of downregulated OCRGs (Figures 10A,B). We then examined a hypothesis that upregulations of OCRGs are correlated with increased immune cell infiltration. To examine this hypothesis, GSCA database (http://bioinfo.life.hust.edu.cn/GSCA/#/immune) (76) was used to analyze the correlation of upregulated OCRGs with immune infiltrates. Liver hepatocellular carcinoma (LIHC) and lung squamous cell carcinoma (LUSC) were studied for the details as an example. The result in Figure 10C shows that almost 19 upregulated OCRGs in LIHC were positively correlated with B cells, natural killer (NK) cells, effector-memory T cells, naturally occurring Treg (nTreg), CD8 T cells, dendritic cells (DCs), type 1 Treg (Tr1), and type 1 T helper cells (Th1) infiltration but negatively correlated with macrophages, neutrophils, CD4 T cells, inducible Treg (iTreg), central memory T cells, monocytes, Th17, and CD4 naïve cell infiltrations. The expressions of BAX, OPTN, GRN, CTSA, CLTA, and CD24 were positively correlated with infiltration scores; and the expressions of TOMM20 were negatively correlated with infiltration scores. Almost 14 upregulated OCRGs in LUSC were positively correlated with CD8 T cells, nTreg, and effector-memory T-cell infiltrations and were negatively correlated with CD4 T cells, gamma-delta T cells, NK cells, mucosal-associated invariant T cells (MAIT), cytotoxic macrophages, NK cells, and follicular T helper (Tfh) cells infiltrations. The expressions of all 14 OCRGs were negatively correlated with immune cell infiltration scores (Figure 10D). Therefore, our findings have demonstrated that the expression modulations of OCRGs can predict the degrees of tumor immune infiltrations and serve as potential therapeutic targets. The GO (Gene Ontology) chord analysis showed that more upregulated OCRGs in LIHC and LUSC were enriched in mitochondrion organization pathway than other pathways (Figures 10E,F). In fact, the upregulated OCRGs in most tumors can be enriched in mitochondrion organization pathways (Supplementary Figure 12).

FIGURE 10

Figure 10. The ratio of upregulated and downregulated OCRGs has different dominance in tumors of different systems, and differentially expressed OCRGs are correlated with immune infiltrates in liver hepatocellular carcinoma (LIHC) and lung squamous cell carcinoma (LUSC). (A) Tumors are listed according to the ratio of upregulated and downregulated OCRGs. The results show that in pancreatic adenocarcinoma (PAAD), brain lower-grade glioma (LCC), and esophageal carcinoma (ESCA), these 15 cancers [from PAAD to uterine corpus endometrial carcinoma (UCEC)] upregulated OCRGs more than the downregulated OCRGs; in lung adenocarcinoma (LUAD), adrenocortical carcinoma (ACC), and acute myeloid leukemia (LAML), these 13 cancers (from HNSC to LUAD) downregulated OCRGs more than the upregulated genes. In order to rule out the influence of the overall expression profile on our interest OCRGs, we take the ratio of number of OCRGs and number of total differently expressed genes as the comparison criteria for the number of upregulated and downregulated genes. Result show that 16 cancers have upregulated OCRG percentage >1% and that 12 cancers have upregulated OCRGs percentages <1%. (B) Expression of OCRGs in tumors was roughly classified according to the system. The numbers of upregulated OCRGs are more than downregulated genes in tumors of immune system, digestive system, and nervous system. The numbers of downregulated OCRGs are more than upregulated OCRGs in tumors of the endocrine system, respiratory system, and circulation system and most tumors of the reproductive system and urinary system. In breast invasive carcinoma (BRCA) and skin cutaneous melanoma (SKCM), the numbers of upregulated OCRGs are more than those of downregulated OCRGs. In head and neck squamous cell carcinoma (HNSC), the numbers of upregulated OCRGs are less than those of downregulated OCRGs. (C) Almost 19 upregulated OCRGs in LIHC are positively correlated with B cell, NKT, effector memory, nTreg, CD8 T, DC, Tr1, and Th1 infiltration but negatively correlated with macrophage, neutrophil, CD4 T, iTreg, central memory, monocyte, Th17, and CD4 naïve cell infiltration. Expressions of BAX, OPTN, GRN, CTSA, CLTA, and CD24 are positively correlated with infiltration score, and expression of TOMM20 is negatively correlated with infiltration score. (D) Fourteen upregulated OCRGs in LUSC are positively correlated with CD8 T, nTreg, and effector-memory cell infiltration and negatively correlated with CD4 T, gamma-delta, NKT, MAIT, cytotoxic macrophage, NK, and Tfh cell infiltration. Expressions of all 14 OCRGs are negatively correlated with infiltration score. Correlations between expression of upregulated OCRGs and immune infiltrates were analyzed in LIHC and LUSC in Gene Set Cancer Analysis (GSCA) database (http://bioinfo.life.hust.edu.cn/GSCA/#/immune). (E) Fifteen out 19 upregulated OCRGs in LIHC enriched in seven significant pathways including regulation of protein stability, protein localization to membrane, mitochondrial fusion, TOM mitochondrial, positive regulation of apoptotic process, regulation of mitochondrion organization, and mitochondrion organization. More OCRGs enriched in TOM complex mitochondrial and mitochondrion organization (GOterm4 and GOterm7). (F) Eleven out of 14 OCRGs in LUSC enriched in five significant pathways including mitochondrial fission, regulation of T cell activation, divalent metal ion transport, mitochondrial membrane organization, and mitochondrion organization. More OCRGs (TOMM40, APOO, PARL, HSPD1, MTFP1, TOMM22, OPA1, and MTFR2) enriched mitochondrion organization pathway (GOterm5). GO (Gene Ontology) chord was generated by http://www.bioinformatics.com.cn, an online platform for data analysis and visualization.

We further examined five digestive cancers such as PAAD, stomach adenocarcinoma (STAD), rectal adenocarcinoma (READ), LIHC, and colon adenocarcinoma (COAD) among the cancer groups with >1% OCRGs and performed the Venn diagram analysis on upregulated OCRGs of these five digestive cancers. As shown in Supplementary Figure 13A, the results showed that PAAD had 80 specific OCRGs (54.8%) out of 146 upregulated OCRGs. In contrast, the other cancers had much lower specific OCRGs. For example, STAD, READ, LIHC, and COAD had three (5.6%), one (2.8%), two (11.1%), and zero (0%) specific OCRGs, respectively. Of note, one OCRG shared by all five cancer groups was translocase of outer MT membrane 40 (TOMM40). We also found that TOMM40 was upregulated in 14 out of 33 types of cancers (42.4%) (Supplementary Figure 13B). Our results are well-correlated with previous reports that TOMM40 (TOM40) and significant numbers of MT protein import machineries overexpressed in cancers such as the TOM complex (TOM20, TOM40, TOM7, and TOM70), TIM23 complex, TIM22 complex, small Tim chaperons and Mia40, and Tim 44 (111, 114). To elucidate the functional aspects of TOMM40, we searched the PPI network database in the STRING database (https://string-db.org/). As shown in Supplementary Figure 13C, the top 10 proteins functionally interacted with TOMM40 including translocase of outer MT membrane (TOM) complex members TOMM7, TOMM20, TOMM70A, and TOMM22; translocase of inner MT membrane (TIM) complex members TIMM50, TIMM44, TIMM13, and TIMM10; sorting and assembly machinery complex (SAM) member SAMM50; and CYC1 (cytochrome C1, complex III subunit 4). As shown in Supplementary Figures 13D,E, the most relative top 10 genes to TOMM40 were localized in two types of protein complexes, the translocase of outer MT membrane complexes (TOMs) and translocase of inner MT membrane complexes (TIMs), and were mainly enriched in protein targeting to mitochondria, protein transporter activity and MT protein complex in biology process, molecular function and cell component (115, 116). The enriched pathways of these 11 genes were MT protein import, Pink/Parkin-mediated mitophagy, and Ub-specific processing proteases (117). We further determined whether increased expressions of TOMM40 and 10 interaction proteins in cancers lead to better prognosis by analyzing the overall survival map of hazard ratio (HR) and disease-free survival map [relapse-free survival (RFS)] of TOMM40 and 10 interaction proteins. The result showed that overall survival of patients in most cancers were positively related to these 11 genes (red box in Supplementary Figure 13F); and RFS of patients in most cancers was positively related to upregulation of these 11 genes (red box in Supplementary Figure 13G; the blue box was the negative correlation between the expression of genes and RFS). It has been reported that MT protein transport system is modulated in cancers (118). A previous report showed that TOMM70 acts as receptor of the MT antiviral-signaling protein (MAVS) and thereby participates in the corresponding system of innate immunity against viral infections (119). It remains unknown whether this innate immune function of TOMM70 is correlated with the better prognosis of cancers upregulating TOMM40 complexes. Relations between abundance of tumor-infiltrating lymphocytes, immuno-inhibitors, immuno-stimulators, MHCs, and expression of TOMM40 were explored by using TISIDB database (http://cis.hku.hk/TISIDB/index.php) (75). The results indicated that the expressions of TOMM40 were positively correlated with abundance of act CD8, act CD4, CD56, and monocyte cells in most of the cancers and that immuno-inhibitor PVRL2; immuno-stimulators CD276, PVR, TNFRSF18, TNFRSF25, and TNFRSF4; and HLA-A, HLA-B, HLA-C, HLA-F, TAP1, TAP2, and TAPBP molecules were positively correlated with expression of TOMM40 in most of cancers (Supplementary Figure 14). We speculate that upregulation of TOMM40 may improve the prognosis of patients by promoting tumor immunity.

Tumor Promoter Factor and Tumor Suppressor Significantly Modulate the Expressions of Organelle Crosstalk Regulators; Vesicle Regulators May Attenuate or Increase Tumorigenesis When Tumor Promoter Factor or Tumor Suppressor Is Deficient

We then hypothesized that oncogenes and tumor suppressors play important roles in regulating the expressions of OCRGs in cancer cells. To test this hypothesis, we collected a total of 19 microarrays from the deficiencies of oncogenes and tumor suppressors in NIH–NCBI–GEO dataset database (https://www.ncbi.nlm.nih.gov/gds/). Since NF-kB and inflammation have been linked to Kras mutations and cancer development, inhibitor of NF-kB kinase subunit b (IKK2) depletion in lung tumor cells significantly attenuates tumor proliferation and prolongs mouse survival (120). As shown in Table 11, IKK2 KO vs. wild-type control upregulated five OCRG and downregulated five OCRGs in mouse lung tumor lines. Four additional IKK and IKK2 knockout (KO) in mouse lung tumor nodules, mouse and human leukemia cells, and human MDA-MB-231 breast cancer epithelial cells did not significantly modulate OCRGs. Deficiencies of NF-kB subunit Rela in mouse lung carcinoma cells and human MDA-MD-231 cells upregulated one and five OCRGs and downregulated two OCRGs. In addition, the Janus kinases (JAKs)/signal transducer and activator of transcription proteins (STATs) signaling pathway provide important roles in contributing to oncogenesis in hematological malignancies and solid tumors (121). JAK2 knockdown (KD) in human acute myelogenous leukemia (AML) cell line, STAT1 KD in human T cell acute lymphocytic leukemia cells, Stat3 KD in mouse mammary tumor, and STAT3 KD in human urothelial cancer cells resulted in zero to two OCRG upregulation and one to five OCRG downregulations. Moreover, tumor suppressor and transcription factor Tp53 KO in mouse liver tumor (122), mouse mammary epithelial cells, oviductal cells, neural stem cells, small intestine, and neu primary tumor led to upregulation of zero up to 23 OCRGs and downregulation of zero up to 48 OCRGs, respectively. Of note, adenomatous polyposis coli (APC) is a multi-functional tumor suppressor gene (123); and a limited role for Tp53 is found in modulating early tumor progression induced by APC loss in mouse intestine (124). APC KO led to six OCRG upregulation and three OCRG downregulation compared with that of wild-type controls, respectively. Furthermore, phosphatidylinositol 3-kinase (PI3-kinase) and phosphatase and tensin homolog deleted on chromosome 10 (PTEN) are major positive and negative regulators of PI3 kinase pathway in regulating cell growth, survival, and proliferation and are two of the most frequently mutated proteins in human cancers (125). PTEN is a tumor suppressor (126); and PTEN KD and KO resulted in zero to seven OCRG upregulation and zero to four OCRG downregulation. In all the significantly modulated OCRGs, there were 13 upregulated and 16 downregulated OCRGs in the oncogene deficiency databases examined and 47 upregulated and 50 downregulated OCRGs in tumor suppressors deficiency databases (Supplementary Figures 15A,B, Supplementary Tables 11, 12 listed, in detail, the modulated OCRGs). These results have demonstrated that in tumor conditions, upregulated oncogenes positively increased the expressions of 16 upregulated OCRGs and 13 downregulated OCRGs. In contrast, tumor suppressors deficiencies upregulated 47 OCRGs in these datasets. Upregulated and downregulated OCRGs were classified into seven groups in oncogene deficiencies and were classified into 12 and 14 groups in tumor suppressors deficiencies (Supplementary Figure 15C).

Our data in Table 11 show that IKK2 depletion upregulated five OCRGs in lung tumor cells, which were well-correlated with the finding that IKK2 depletion attenuated tumor proliferation and prolonged mouse survival. We then tested a hypothesis that IKK2 depletion upregulated OCRGs may partially overlap with LUAD-downregulated OCRGs. As shown in Supplementary Figure 16A, the Venn diagram indicated that vesicle regulator perilipin 4 (PLIN4) was overlapped between the LUAD-downregulated OCRGs and the IKK2 KO-upregulated OCRGs. As shown in Supplementary Figure 16B, the expression profile analysis using cancer and normal gene expression database GEPIA (http://gepia.cancer-pku.cn/index.html) (127) showed that the expression of PLIN4 was decreased in 23 out of 33 tumor types compared with controls. Our data in Table 8 show that tumor suppressor Tp53 KO in liver tumors upregulated 23 OCRGs and downregulated 48 OCRGs, which were associated with induction of carcinomas in Tp53 KO mice (122, 126). Similarly, we then tested a hypothesis that Tp53 KO upregulated OCRGs may partially overlap with that liver tumor (LIHC)-upregulated OCRGs; and Tp53 KO downregulated OCRGs may partially overlap with liver tumor (LIHC)-downregulated OCRGs. As shown in Supplementary Figure 16C, the Venn diagram indicated that prolyl 4-hydroxylase subunit alpha 2 (P4HA2) was shared between the Tp53 KO upregulated OCRGs and the liver tumor (LIHC)-upregulated OCRGs; and alpha-aminoadipic semialdehyde synthase (AASS) was shared between the Tp53 KO downregulated OCRGs and the liver tumor (LIHC)-downregulated OCRGs. The cancer gene expression database analysis showed that the expressions of vesicle regulator P4HA2 were increased in seven out of 33 tumor types compared with controls (Supplementary Figure 16D); and the expression of vesicle regulator AASS in seven out of 33 tumor types was decreased compared with that of controls (Supplementary Figure 16E). Taken together, our results have demonstrated that first, tumor promoter IKK2 and tumor suppressor Tp53 significantly modulate the expressions of OCRGs; second, lung tumor LUAD-downregulated vesicle regulator PLIN4 may contribute to IKK2 KO-promoted tumor formation; third, liver tumor (LIHC)-upregulated vesicle regulator P4HA2 may increase Tp53 KO-promoted liver tumorigenesis; and fourth, LIHC-downregulated vesicle regulator AASS may contribute to Tp53 suppression of liver tumorigenesis.

Discussion

Although alterations to these networks, such as impaired ER–mitochondria MCSs, have been linked to several diseases such as neurodegeneration (59, 60), CVD (61), diabetes (62), kidney diseases (63, 64), and cancers (65, 66), an important question remains poorly characterized how the expressions of organelle crosstalk regulator genes are modulated in various diseases. To examine whether the expressions of 260 OCRGs in 16 functional groups are modulated in 23 diseases and 28 tumors, we performed extensive -omics data mining analyses in a panoramic manner with the method that we pioneered in 2004 (23, 26, 29, 91, 107, 109, 128) and made a set of significant findings: (1) the ratios of upregulated vs. downregulated OCRGs are 1:2.8 in AIs, 1:1 in MDs, 1:1.1 in ADs, and 1:3.8 in OFs; and AIs and OFs upregulate less OCRG groups (three to seven out of 16) than metabolic and ADs (11 to 13 out of 16); (2) sepsis and trauma-upregulated OCRG groups including vesicle, MT fission, and mitophagy but not other groups of OCRGs; thus, the three groups of OCRGs are classified as the acute crisis-handling OCRGs; and similarly, sepsis and trauma plus OFs-upregulated seven OCRG groups including vesicle, MT fission, mitophagy, sarcoplasmic reticulum–MT, MT fusion, autophagosome–lysosome fusion, and autophagosome/endosome–lysosome fusion are classified as the cell failure-handling OCRGs; (3) the majority of upregulated pathways are disease group-specific; upregulated organelle fusion and macroautophagy are shared by MDs, ADs, and OFs; upregulated MT fission is shared by AIs, MDs, and ADs; upregulated vesicle organization is shared by AI and MDs; and upregulated protein localization is shared by MDs and ADs; (4) increased OCRG expressions in all the 15 functional groups except autophagosome–lysosome fusion but decreased autophagosome–lysosome fusion are required for viral replications, which classify this decreased group as the viral replication-suppressed OCRGs. (5) PAMPs from influenza virus and Kaposi sarcoma-associated herpes virus infections and pro-atherogenic DAMPs such as oxLDL, LPS, oxPAPC, and IFNs totally upregulated 33 OCRGs in ECs, which classify vesicle, MT fission, mitophagy, MT fusion, ER–MT contact, ER–PT junction, and autophagosome/endosome–lysosome fusion as the seven EC-activation/inflammation-promoting OCRG groups; (6) the expression of OCRGs is modulated in Treg from the LNs, spleen, peripheral blood, intestine, and brown adipose tissue in comparison with that of CD4⁺CD25⁻ T effector controls; (7) TLR2, TLR4, TLR3/7/9, caspase-1, and ROS regulator Nrf2 can upregulate 46 OCRGs from 14 groups out of all 16 groups (except for MT fission and fusion, and ER–MT contact OCRGs), suggesting that DAMP/PAMP-sensing systems play significant roles in modulating the expressions of OCRGs; (8) OCRG expressions are significantly modulated in all the 28 cancer datasets; upregulated OCRGs in LIHC and LUSC are correlated with immune cell infiltrations, suggesting that upregulation of certain OCRGs may promote immune cell infiltration; MT protein translocase TOMM40 and MT import machine are upregulated in cancers and are positively associated with better prognosis; and TOMM40 is correlated with some immune infiltration signature genes, suggesting that increase of protein transport into mitochondria may facilitate the immune cell infiltration; and (9) tumor promoter factor IKK2 and tumor suppressor Tp53 significantly modulate the expressions of OCRGs; lung tumor-downregulated vesicle regulator PLIN4 may contribute to IKK2 KO-attenuated tumor formation; liver tumor-upregulated vesicle regulator P4HA2 may increase Tp53 deficiency-promoted liver tumorigenesis; and liver tumor-downregulated vesicle regulator AASS may contribute to Tp53 suppression of liver tumorigenesis.

The original microarray experiments we analyzed used different cells, which prevented us from comparing the effects of various regulators in regulating the expressions of OCRGs in the same cell types, although our database mining method was not ideal. However, to fill in the important knowledge gaps, our approach was justified. Actually, this approach has been a common practice that we (23) and others (72, 129) often use in studying gene expression in non-ideal, heterogeneous PBMC populations in disease conditions vs. healthy conditions, and PBMCs are actually composed of many cell types, such as B cells (~15%), T cells (~70%), monocytes (~5%), and NK cells (~10%) (130).

To summarize our findings from a panoramic analysis on the expression changes of 260 OCRGs in 51 diseases and cancers presented here, we propose a novel working model to integrate these results, as follows, in Figure 11A: first, the MCSs between two organelles, at least between ER and mitochondria (Mito, MT), have been found to play several significant functions such as metabolism and Ca²⁺ transport as we reported (4); organelle communicating ROS as we recently reported (4, 38, 40–42, 131); apoptosis as we have reported (132, 133); and autophagy, mitophagy, and DAMP communication as we have reported (37). The metabolism regulatory functions between ER–MT contact sites include gluconeogenesis, glycogen synthesis and breakdown, fat storage, hormone metabolism, drug metabolism, lipid biosynthesis, and homeostasis (134); second, since our study analyzed the 260 OCRGs in 51 major diseases and cancers and many other master regulator deficient transcriptomic data, we have demonstrated for the first time that 16 functional groups of OCRGs are not equally modulated in various diseases and cancers. As shown in Figure 11B, we identified three groups of OCRGs such as vesicle, MT fission, and mitophagy as the cell crisis-handling OCRGs, at least in partial functions of the three groups; seven groups of OCRGs including three cell crisis-handling OCRGs plus sarcoplasmic reticulum–MT, MT fusion, autophagosome–lysosome fusion, and autophagosome/endosome–lysosome fusion as the cell failure-handling OCRGs; and autophagosome–lysosome fusion as virus-suppressed OCRGs. Previous reports showed that growing numbers of viruses are found to weaponize the ubiquitin modification system to suppress anti-viral type I IFN pathways (135). Similarly, we found that virus infections weaponize the suppression of autophagosome–lysosome fusion group potentially for virus replication. Since most viruses use endocytosis to enter the host cell, 11 clinically approved generic drugs are identified as potential candidates for repurposing as blockers of several potential routes for severe acute SARS-CoV-2 endocytosis (136), which needs to be further re-considered since we found that coronavirus and other virus infections in ECs and epithelial cells suppress endocytosis–lysosome pathways (Supplementary Figure 4); 10 groups of OCRGs activated by DAMPs in ECs including vesicle (137), MT fission, mitophagy, autophagosome/endosome–lysosome fusion, ER–MT contact, ER–PM junctions, sarcoplasmic reticulum–MT, MT fusion, ER–endosome, and ER–GC interaction. Treg plastic OCRGs, and oncogenes and tumor suppressors promoted vesicle regulators; and third, extracellular PAMPs/DAMPs, cytokines, and viruses act via their receptors including TLRs on the cell membrane and intracellular locations (Nod-like receptors and inflammasomes) (90) to activate OCRGs directly or indirectly by activating intracellular organelle dangers first and then OCRGs. Therefore, similar to our recently proposed new concept that ROS are an integrated system for sensing and alarming organelle metabolic homeostasis and dysfunction (38), here, we propose a new concept that the organelle crosstalk serves as both new DAMPs sensing and communicating network and new cell crisis and cell failure-handling network (37) to fulfill both physiological functions and pathological functions (Figure 11C).

FIGURE 11

Figure 11. A novel working model shows the function, classification, and regulated mechanisms of OCRGs in diseases, and DAMP- and virus-treated cells. (A) The membrane contact sites between two organelles, at least between ER and mitochondria (Mito, MT), have been found to play several significant functions from previous reports and our study, such as in metabolism (including gluconeogenesis, glycogen synthesis and breakdown, fat storage, hormone metabolism, drug metabolism, lipid biosynthesis, and homeostasis), Ca²⁺ transport, organelle communication ROS, apoptosis, autophagy/mitophagy, and DAMP communication. (B) Thirteen of 16 functional groups of OCRGs are differently modulated in various diseases and cancers. Three groups of OCRGs such as vesicle, mitochondrial fission, and mitophagy were identified as the cell crisis-handling OCRGs (red line pointed). Seven groups of OCRGs including three cell crisis-handling OCRGs plus sarcoplasmic reticulum-mito, mitochondrial fusion, autophagosome–lysosome fusion, and autophagosome/endosome–lysosome fusion were identified as the cell failure-handling OCRGs (blue line pointed). Autophagosome–lysosome fusion OCRGs can be regarded as virus infection-suppressed OCRGs (brown line pointed). Eight classification groups (i.e., vesicle, MT fission, mitophagy, sarcoplasmic reticulum–MT, autophagosome–lysosome fusion, ER–MT contact, ER–PM junctions, and ER–endosome) are Treg plastic OCRGs (purple line pointed). Ten groups of ECs activation by DAMP OCRGs include seven Treg plastic OCRGs (except autophagosome–lysosome fusion) and add MT fusion, autophagosome/endosome–lysosome fusion, and ER–GC interaction OCRGs (green line pointed). Vesicle regulators can be regarded as oncogenes and tumor suppressor-promoted OCRGs (light blue line pointed). (C) Extracellular pathogen-associated molecular patterns (PAMPs)/damage-associated molecular patterns (DAMPs), cytokines, and viruses act via their receptors including TLRs on the cell membrane and intracellular locations (Nod-like receptors and inflammasomes) to activate OCRGs directly or indirectly by activating intracellular organelle dangers first and then OCRGs. And ROS regulator NRF2 oncogenes and tumor suppressors also regulate the expression of OCRGs in diseases and tumors. A new concept that the organelle crosstalk serves as both new DAMP sensing and communicating network and new cell crisis and cell failure-handling network to fulfill both physiological functions and pathological functions. *Created with BioRender.com.

One limitation of the current study is that due to the low throughput nature of verification techniques in every laboratory including ours, we could not verify every result we identified with the analyses of high-throughput data, which are similar to all the papers with RNA-Seq, single-cell RNA-Seq, and other -omics data. We acknowledge that carefully designed in vitro and in vivo experimental models will be needed in the future to verify regulator gene deficiency-upregulated OCRGs further and the underlying mechanisms we report here. Nevertheless, our findings provide novel insights on the roles of upregulated OCRGs in the pathogenesis of inflammatory diseases and cancers, novel pathways for the future therapeutic interventions for inflammations, sepsis, trauma, OFs, ADs, metabolic CVDs, and cancers.

Data Availability Statement

The datasets generated for this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

Author Contributions

ML carried out the data gathering and data analysis and prepared the tables and figures. NW, KX, FS, EV, YS, RZ, JW, HS, WY, YL, YS, CD, LLiu, LLi, WH, JY, DP, JS, XJ, and HW aided with analysis of the data. XY supervised the experimental design, data analysis, and manuscript writing. All authors read and approved the final manuscript.

Funding

This work was partially supported by NIH grants to HW and XY.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

ML was supported by a fellowship from the School of Basic Medical Science, Shanxi Medical University.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcvm.2021.713170/full#supplementary-material

References

1. Benjamin EJ, Blaha MJ, Chiuve SE, Cushman M, Das SR, Deo R, et al. Heart disease and stroke statistics-2017 update: a report from the American Heart Association. Circulation. (2017) 135:e146–603. doi: 10.1161/CIR.0000000000000491

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Libby P, Buring JE, Badimon L, Hansson GK, Deanfield J, Bittencourt MS, et al. Atherosclerosis. Nat Rev Dis Primers. (2019) 5:56. doi: 10.1038/s41572-019-0106-z

CrossRef Full Text | Google Scholar

3. Yin Y, Li X, Sha X, Xi H, Li YF, Shao Y, et al. Early hyperlipidemia promotes endothelial activation via a caspase-1-sirtuin 1 pathway. Arterioscler Thromb Vasc Biol. (2015) 35:804–16. doi: 10.1161/ATVBAHA.115.305282

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Li X, Fang P, Li Y, Kuo YM, Andrews AJ, Nanayakkara G, et al. Mitochondrial reactive oxygen species mediate lysophosphatidylcholine-induced endothelial cell activation. arterioscler Thromb Vasc Biol. (2016) 36:1090–100. doi: 10.1161/ATVBAHA.115.306964

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Fang P, Zhang D, Cheng Z, Yan C, Jiang X, Kruger WD, et al. Hyperhomocysteinemia potentiates hyperglycemia-induced inflammatory monocyte differentiation and atherosclerosis. Diabetes. (2014) 63:4275–90. doi: 10.2337/db14-0809

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Xi H, Zhang Y, Xu Y, Yang WY, Jiang X, Sha X, et al. Caspase-1 inflammasome activation mediates homocysteine-induced pyrop-apoptosis in endothelial cells. Circ Res. (2016) 118:1525–39. doi: 10.1161/CIRCRESAHA.116.308501

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Yang J, Fang P, Yu D, Zhang L, Zhang D, Jiang X, et al. Chronic kidney disease induces inflammatory CD40+ monocyte differentiation via homocysteine elevation and DNA hypomethylation. Circ Res. (2016) 119:1226–41. doi: 10.1161/CIRCRESAHA.116.308750

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Ferrer LM, Monroy AM, Lopez-Pastrana J, Nanayakkara G, Cueto R, Li YF, et al. Caspase-1 plays a critical role in accelerating chronic kidney disease-promoted neointimal hyperplasia in the carotid artery. J Cardiovasc Transl Res. (2016) 9:135–44. doi: 10.1007/s12265-016-9683-3

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Monroy MA, Fang J, Li S, Ferrer L, Birkenbach MP, Lee IJ, et al. Chronic kidney disease alters vascular smooth muscle cell phenotype. Front Biosci. (2015) 20:784–95. doi: 10.2741/4337

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Sun Y, Johnson C, Zhou J, Wang L, Li YF, Lu Y, et al. Uremic toxins are conditional danger- or homeostasis-associated molecular patterns. Front Biosci. (2018) 23:348–87. doi: 10.2741/4595

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Shao Y, Saredy J, Yang WY, Sun Y, Lu Y, Saaoud F, et al. Vascular endothelial cells and innate immunity. Arterioscler Thromb Vasc Biol. (2020) 40:e138–52. doi: 10.1161/ATVBAHA.120.314330

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Sha X, Meng S, Li X, Xi H, Maddaloni M, Pascual DW, et al. Interleukin-35 inhibits endothelial cell activation by suppressing MAPK-AP-1 pathway. J Biol Chem. (2015) 290:19307–18. doi: 10.1074/jbc.M115.663286

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Shao Y, Cheng Z, Li X, Chernaya V, Wang H, Yang XF. Immunosuppressive/anti-inflammatory cytokines directly and indirectly inhibit endothelial dysfunction–a novel mechanism for maintaining vascular function. J Hematol Oncol. (2014) 7:80. doi: 10.1186/s13045-014-0080-6

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Li X, Wang L, Fang P, Sun Y, Jiang X, Wang H, et al. Lysophospholipids induce innate immune transdifferentiation of endothelial cells, resulting in prolonged endothelial activation. J Biol Chem. (2018) 293:11033–45. doi: 10.1074/jbc.RA118.002752

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Li A, Sun Y. Drummer Ct, Lu Y, Yu D, Zhou Y, et al. Increasing upstream chromatin long-range interactions may favor induction of circular RNAs in LysoPC-activated human aortic endothelial cells. Front Physiol. (2019) 10:433. doi: 10.3389/fphys.2019.00433

CrossRef Full Text | Google Scholar

16. Lopez-Pastrana J, Ferrer LM Li YF, Xiong X, Xi H, Cueto R, et al. Inhibition of caspase-1 activation in endothelial cells improves angiogenesis: a novel therapeutic potential for ischemia. J Biol Chem. (2015) 290:17485–94. doi: 10.1074/jbc.M115.641191

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Zhang CE, Wei W, Liu YH, Peng JH, Tian Q, Liu GP, et al. Hyperhomocysteinemia increases beta-amyloid by enhancing expression of gamma-secretase and phosphorylation of amyloid precursor protein in rat brain. Am J Pathol. (2009) 174:1481–91. doi: 10.2353/ajpath.2009.081036

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Zhang D, Fang P, Jiang X, Nelson J, Moore JK, Kruger WD, et al. Severe hyperhomocysteinemia promotes bone marrow-derived and resident inflammatory monocyte differentiation and atherosclerosis in LDLr/CBS-deficient mice. Circ Res. (2012) 111:37–49. doi: 10.1161/CIRCRESAHA.112.269472

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Fang P, Li X, Shan H, Saredy JJ, Cueto R, Xia J, et al. Ly6C(+) inflammatory monocyte differentiation partially mediates hyperhomocysteinemia-induced vascular dysfunction in type 2 diabetic db/db mice. Arterioscler Thromb Vasc Biol. (2019) 39:2097–119. doi: 10.1161/ATVBAHA.119.313138

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Johnson C, Drummer C IV, Shan H, Shao Y, Sun Y, Lu Y, et al. A novel subset of CD95(+) pro-inflammatory macrophages overcome miR155 deficiency and may serve as a switch from metabolically healthy obesity to metabolically unhealthy obesity. Front Immunol. (2020) 11:619951. doi: 10.3389/fimmu.2020.619951

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Lai B, Wang J, Fagenson A, Sun Y, Saredy J, Lu Y, et al. Twenty novel disease group-specific and 12 new shared macrophage pathways in eight groups of 34 diseases including 24 inflammatory organ diseases and 10 types of tumors. Front Immunol. (2019) 10:2612. doi: 10.3389/fimmu.2019.02612

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Drummer CT, Saaoud F, Shao Y, Sun Y, Xu K, Lu Y, et al. Trained immunity and reactivity of macrophages and endothelial cells. Arterioscler Thromb Vasc Biol. (2021) 41:1032–46. doi: 10.1161/ATVBAHA.120.315452

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Zhang R, Saredy J, Shao Y, Yao T, Liu L, Saaoud F, et al. End-stage renal disease is different from chronic kidney disease in upregulating ROS-modulated proinflammatory secretome in PBMCs - A novel multiple-hit model for disease progression. Redox Biol. (2020) 34:101460. doi: 10.1016/j.redox.2020.101460

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Zhang R, Xu K, Shao Y, Sun Y, Saredy J, Cutler E, et al. Tissue Treg secretomes and transcription factors shared with stem cells contribute to a Treg niche to maintain Treg-ness with 80% innate immune pathways, and functions of immunosuppression and tissue repair. Front Immunol. (2020) 11:632239. doi: 10.3389/fimmu.2020.632239

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Li X, Fang P, Yang WY, Wang H, Yang X. IL-35, as a newly proposed homeostasis-associated molecular pattern, plays three major functions including anti-inflammatory initiator, effector, and blocker in cardiovascular diseases. Cytokine. (2019) 122:154076. doi: 10.1016/j.cyto.2017.06.003

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Li X, Mai J, Virtue A, Yin Y, Gong R, Sha X, et al. IL-35 is a novel responsive anti-inflammatory cytokine–a new system of categorizing anti-inflammatory cytokines. PLoS ONE. (2012) 7:e33628. doi: 10.1371/journal.pone.0033628

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Li X, Shao Y, Sha X, Fang P, Kuo YM, Andrews AJ, et al. IL-35 (Interleukin-35) suppresses endothelial cell activation by inhibiting mitochondrial reactive oxygen species-mediated site-specific acetylation of H3K14 (Histone 3 Lysine 14). Arterioscler Thromb Vasc Biol. (2018) 38:599–609. doi: 10.1161/ATVBAHA.117.310626

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Mai J, Nanayakkara G, Lopez-Pastrana J, Li X, Li YF, Wang X, et al. Interleukin-17A promotes aortic endothelial cell activation via transcriptionally and post-translationally activating p38 mitogen-activated Protein Kinase (MAPK) pathway. J Biol Chem. (2016) 291:4939–54. doi: 10.1074/jbc.M115.690081

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Liu M, Saredy J, Zhang R, Shao Y, Sun Y, Yang WY, et al. Approaching inflammation paradoxes-proinflammatory cytokine blockages induce inflammatory regulators. Front Immunol. (2020) 11:554301. doi: 10.3389/fimmu.2020.554301

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Fu H, Sun Y, Shao Y, Saredy J, Cueto R, Liu L, et al. Interleukin 35 delays hindlimb ischemia-induced angiogenesis through regulating ROS-extracellular matrix but spares later regenerative angiogenesis. Front Immunol. (2020) 11:595813. doi: 10.3389/fimmu.2020.595813

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Xiong Z, Song J, Yan Y, Huang Y, Cowan A, Wang H, et al. Higher expression of Bax in regulatory T cells increases vascular inflammation. Front Biosci. (2008) 13:7143–55. doi: 10.2741/3217

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Xiong Z, Yan Y, Song J, Fang P, Yin Y, Yang Y, et al. Expression of TCTP antisense in CD25(high) regulatory T cells aggravates cuff-injured vascular inflammation. Atherosclerosis. (2009) 203:401–8. doi: 10.1016/j.atherosclerosis.2008.07.041

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Yang WY, Shao Y, Lopez-Pastrana J, Mai J, Wang H, Yang XF. Pathological conditions re-shape physiological Tregs into pathological Tregs. Burns Trauma. (2015) 3:1. doi: 10.1186/s41038-015-0001-0

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Xu K, Yang WY, Nanayakkara GK, Shao Y, Yang F, Hu W, et al. GATA3, HDAC6, and BCL6 regulate FOXP3+ Treg plasticity and determine Treg conversion into either novel antigen-presenting cell-like Treg or Th1-Treg. Front Immunol. (2018) 9:45. doi: 10.3389/fimmu.2018.00045

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Nelson J, Wu Y, Jiang X, Berretta R, Houser S, Choi E, et al. Hyperhomocysteinemia suppresses bone marrow CD34+/VEGF receptor 2+ cells and inhibits progenitor cell mobilization and homing to injured vasculature-a role of beta1-integrin in progenitor cell migration and adhesion. FASEB J. (2015) 29:3085–99. doi: 10.1096/fj.14-267989

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Li YF, Huang X, Li X, Gong R, Yin Y, Nelson J, et al. Caspase-1 mediates hyperlipidemia-weakened progenitor cell vessel repair. Front Biosci. (2016) 21:178–91. doi: 10.2741/4383

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Zeng H, Nanayakkara GK, Shao Y, Fu H, Sun Y, Cueto R, et al. DNA checkpoint and repair factors are nuclear sensors for intracellular organelle stresses-inflammations and cancers can have high genomic risks. Front Physiol. (2018) 9:516. doi: 10.3389/fphys.2018.00516

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Sun Y, Lu Y, Saredy J, Wang X, Drummer Iv C, Shao Y, et al. ROS systems are a new integrated network for sensing homeostasis and alarming stresses in organelle metabolic processes. Redox Biol. (2020) 37:101696. doi: 10.1016/j.redox.2020.101696

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Fagenson AM, Xu K, Saaoud F, Nanayakkara G, Jhala NC, Liu L, et al. Liver ischemia reperfusion injury, enhanced by trained immunity, is attenuated in caspase 1/caspase 11 double gene knockout mice. Pathogens. (2020) 9:879. doi: 10.3390/pathogens9110879

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Li X, Fang P, Mai J, Choi ET, Wang H, Yang XF. Targeting mitochondrial reactive oxygen species as novel therapy for inflammatory diseases and cancers. J Hematol Oncol. (2013) 6:19. doi: 10.1186/1756-8722-6-19

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Nanayakkara GK, Wang H, Yang X. Proton leak regulates mitochondrial reactive oxygen species generation in endothelial cell activation and inflammation - A novel concept. Arch Biochem Biophys. (2019) 662:68–74. doi: 10.1016/j.abb.2018.12.002

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Cheng J, Nanayakkara G, Shao Y, Cueto R, Wang L, Yang WY, et al. Mitochondrial proton leak plays a critical role in pathogenesis of cardiovascular diseases. Adv Exp Med Biol. (2017) 982:359–70. doi: 10.1007/978-3-319-55330-6_20

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Shao Y, Chernaya V, Johnson C, Yang WY, Cueto R, Sha X, et al. Metabolic diseases downregulate the majority of histone modification enzymes, making a few upregulated enzymes novel therapeutic targets–“sand out and gold stays”. J Cardiovasc Transl Res. (2016) 9:49–66. doi: 10.1007/s12265-015-9664-y

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Zhong C, Yang X, Feng Y, Yu J. Trained immunity: an underlying driver of inflammatory atherosclerosis. Front Immunol. (2020) 11:284. doi: 10.3389/fimmu.2020.00284

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Lu Y, Sun Y, Drummer Ct, Nanayakkara GK, Shao Y, Saaoud F, et al. Increased acetylation of H3K14 in the genomic regions that encode trained immunity enzymes in lysophosphatidylcholine-activated human aortic endothelial cells - Novel qualification markers for chronic disease risk factors and conditional DAMPs. Redox Biol. (2019) 24:101221. doi: 10.1016/j.redox.2019.101221

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Li X, Fang P, Sun Y, Shao Y, Yang WY, Jiang X, et al. Anti-inflammatory cytokines IL-35 and IL-10 block atherogenic lysophosphatidylcholine-induced, mitochondrial ROS-mediated innate immune activation, but spare innate immune memory signature in endothelial cells. Redox Biol. (2020) 28:101373. doi: 10.1016/j.redox.2019.101373

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Yang Q, Zhang R, Tang P, Sun Y, Johnson C, Saredy J, et al. Ultrasound may suppress tumor growth, inhibit inflammation, and establish tolerogenesis by remodeling innatome via pathways of ROS, immune checkpoints, cytokines, and trained immunity/tolerance. J Immunol Res. (2021) 2021:6664453. doi: 10.1155/2021/6664453

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Hu F, Zeng C, Long R, Miao Y, Wei L, Xu Q, et al. Supermultiplexed optical imaging and barcoding with engineered polyynes. Nat Methods. (2018) 15:194–200. doi: 10.1038/nmeth.4578

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Valm AM, Cohen S, Legant WR, Melunis J, Hershberg U, Wait E, et al. Applying systems-level spectral imaging and analysis to reveal the organelle interactome. Nature. (2017) 546:162–7. doi: 10.1038/nature22369

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Shai N, Yifrach E, van Roermund CWT, Cohen N, Bibi C, L IJ, et al. Systematic mapping of contact sites reveals tethers and a function for the peroxisome-mitochondria contact. Nat Commun. (2018) 9:1761. doi: 10.1038/s41467-018-03957-8

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Salvador-Gallego R, Hoyer MJ, Voeltz GK. SnapShot: functions of endoplasmic reticulum membrane contact sites. Cell. (2017) 171:1224–e1. doi: 10.1016/j.cell.2017.11.005

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Bohnert M, Schuldiner M. Stepping outside the comfort zone of membrane contact site research. Nat Rev Mol Cell Biol. (2018) 19:483–4. doi: 10.1038/s41580-018-0022-1

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Stenmark H. Rab GTPases as coordinators of vesicle traffic. Nat Rev Mol Cell Biol. (2009) 10:513–25. doi: 10.1038/nrm2728

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Diogo CV, Yambire KF, Fernandez Mosquera L, Branco FT, Raimundo N. Mitochondrial adventures at the organelle society. Biochem Biophys Res Commun. (2018) 500:87–93. doi: 10.1016/j.bbrc.2017.04.124

CrossRef Full Text | Google Scholar

55. Quiros PM, Mottis A, Auwerx J. Mitonuclear communication in homeostasis and stress. Nat Rev Mol Cell Biol. (2016) 17:213–26. doi: 10.1038/nrm.2016.23

CrossRef Full Text | Google Scholar

56. Shpilka T, Haynes CM. The mitochondrial UPR: mechanisms, physiological functions and implications in ageing. Nat Rev Mol Cell Biol. (2018) 19:109–20. doi: 10.1038/nrm.2017.110

PubMed Abstract | CrossRef Full Text | Google Scholar

57. Senft D, Ronai ZA, UPR. autophagy, and mitochondria crosstalk underlies the ER stress response. Trends Biochem Sci. (2015) 40:141–8. doi: 10.1016/j.tibs.2015.01.002

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Fransen M, Lismont C, Walton P. The peroxisome-mitochondria connection: how and why? Int J Mol Sci. (2017) 18:1126. doi: 10.3390/ijms18061126

CrossRef Full Text | Google Scholar

59. Area-Gomez E, Schon EA. Mitochondria-associated ER membranes and Alzheimer disease. Curr Opin Genet Dev. (2016) 38:90–6. doi: 10.1016/j.gde.2016.04.006

CrossRef Full Text | Google Scholar

60. Deus CM, Yambire KF, Oliveira PJ, Raimundo N. Mitochondria-lysosome crosstalk: from physiology to neurodegeneration. Trends Mol Med. (2020) 26:71–88. doi: 10.1016/j.molmed.2019.10.009

PubMed Abstract | CrossRef Full Text | Google Scholar

61. Kaasik A, Veksler V, Boehm E, Novotova M, Minajeva A, Ventura-Clapier R. Energetic crosstalk between organelles: architectural integration of energy production and utilization. Circ Res. (2001) 89:153–9. doi: 10.1161/hh1401.093440

PubMed Abstract | CrossRef Full Text | Google Scholar

62. Castano C, Novials A, Parrizas M. Exosomes and diabetes. Diabetes Metab Res Rev. (2019) 35:e3107. doi: 10.1002/dmrr.3107

CrossRef Full Text | Google Scholar

63. Inoue T, Maekawa H, Inagi R. Organelle crosstalk in the kidney. Kidney Int. (2019) 95:1318–25. doi: 10.1016/j.kint.2018.11.035

CrossRef Full Text | Google Scholar

64. Inagi R. The Implication of organelle cross talk in AKI. Nephron. (2020) 144:634–7. doi: 10.1159/000508639

PubMed Abstract | CrossRef Full Text | Google Scholar

65. Bebelman MP, Smit MJ, Pegtel DM, Baglio SR. Biogenesis and function of extracellular vesicles in cancer. Pharmacol Ther. (2018) 188:1–11. doi: 10.1016/j.pharmthera.2018.02.013

PubMed Abstract | CrossRef Full Text | Google Scholar

66. Xu J, Camfield R, Gorski SM. The interplay between exosomes and autophagy - partners in crime. J Cell Sci. (2018) 131:jcs.215210. doi: 10.1242/jcs.215210

PubMed Abstract | CrossRef Full Text | Google Scholar

67. Joshi AS, Zhang H, Prinz WA. Organelle biogenesis in the endoplasmic reticulum. Nat Cell Biol. (2017) 19:876–82. doi: 10.1038/ncb3579

CrossRef Full Text | Google Scholar

68. Soto-Heredero G, Baixauli F, Mittelbrunn M. Interorganelle communication between mitochondria and the endolysosomal system. Front Cell Dev Biol. (2017) 5:95. doi: 10.3389/fcell.2017.00095

PubMed Abstract | CrossRef Full Text | Google Scholar

69. Lopez-Crisosto C, Pennanen C, Vasquez-Trincado C, Morales PE, Bravo-Sagua R, Quest AFG, et al. Sarcoplasmic reticulum-mitochondria communication in cardiovascular pathophysiology. Nat Rev Cardiol. (2017) 14:342–60. doi: 10.1038/nrcardio.2017.23

PubMed Abstract | CrossRef Full Text | Google Scholar

70. Besprozvannaya M, Dickson E, Li H, Ginburg KS, Bers DM, Auwerx J, et al. GRAM domain proteins specialize functionally distinct ER-PM contact sites in human cells. Elife. (2018) 7:e31019. doi: 10.7554/eLife.31019

CrossRef Full Text | Google Scholar

71. Lees JA, Messa M, Sun EW, Wheeler H, Torta F, Wenk MR, et al. Lipid transport by TMEM24 at ER-plasma membrane contacts regulates pulsatile insulin secretion. Science. (2017) 355:eaah6171. doi: 10.1126/science.aah6171

PubMed Abstract | CrossRef Full Text | Google Scholar

72. Li Y, Oosting M, Smeekens SP, Jaeger M, Aguirre-Gamboa R, Le KTT, et al. A functional genomics approach to understand variation in cytokine production in humans. Cell. (2016) 167:1099–110 e14. doi: 10.1016/j.cell.2016.10.017

PubMed Abstract | CrossRef Full Text | Google Scholar

73. Zhou Y, Zhou B, Pache L, Chang M, Khodabakhshi AH, Tanaseichuk O, et al. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun. (2019) 10:1523. doi: 10.1038/s41467-019-09234-6

PubMed Abstract | CrossRef Full Text | Google Scholar

74. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. (2003) 13:2498–504. doi: 10.1101/gr.1239303

PubMed Abstract | CrossRef Full Text | Google Scholar

75. Ru B, Wong CN, Tong Y, Zhong JY, Zhong SSW, Wu WC, et al. TISIDB: an integrated repository portal for tumor-immune system interactions. Bioinformatics. (2019) 35:4200–2. doi: 10.1093/bioinformatics/btz210

PubMed Abstract | CrossRef Full Text | Google Scholar

76. Liu CJ, Hu FF, Xia MX, Han L, Zhang Q, Guo AY. GSCALite: a web server for gene set cancer analysis. Bioinformatics. (2018) 34:3771–2. doi: 10.1093/bioinformatics/bty411

PubMed Abstract | CrossRef Full Text | Google Scholar

77. Eisenberg E, Levanon EY. Human housekeeping genes, revisited. Trends Genet. (2013) 29:569–74. doi: 10.1016/j.tig.2013.05.010

CrossRef Full Text | Google Scholar

78. de Jonge HJ, Fehrmann RS, de Bont ES, Hofstra RM, Gerbens F, Kamps WA, et al. Evidence based selection of housekeeping genes. PLoS ONE. (2007) 2:e898. doi: 10.1371/journal.pone.0000898

PubMed Abstract | CrossRef Full Text | Google Scholar

79. Targher G, Byrne CD. Obesity: Metabolically healthy obesity and NAFLD. Nat Rev Gastroenterol Hepatol. (2016) 13:442–4. doi: 10.1038/nrgastro.2016.104

PubMed Abstract | CrossRef Full Text | Google Scholar

80. Leites EP, Morais VA. Mitochondrial quality control pathways: PINK1 acts as a gatekeeper. Biochem Biophys Res Commun. (2018) 500:45–50. doi: 10.1016/j.bbrc.2017.06.096

PubMed Abstract | CrossRef Full Text | Google Scholar

81. Schneider A, Kurz S, Manske K, Janas M, Heikenwalder M, Misgeld T, et al. Single organelle analysis to characterize mitochondrial function and crosstalk during viral infection. Sci Rep. (2019) 9:8492. doi: 10.1038/s41598-019-44922-9

PubMed Abstract | CrossRef Full Text | Google Scholar

82. Menachery VD, Mitchell HD, Cockrell AS, Gralinski LE, Yount BL Jr, Graham RL, et al. MERS-CoV accessory ORFs play key role for infection and pathogenesis. mBio. (2017) 8:e00665–17. doi: 10.1128/mBio.00665-17

PubMed Abstract | CrossRef Full Text | Google Scholar

83. Mitchell HD, Eisfeld AJ, Sims AC, McDermott JE, Matzke MM, Webb-Robertson BJ, et al. A network integration approach to predict conserved regulators related to pathogenicity of influenza and SARS-CoV respiratory viruses. PLoS ONE. (2013) 8:e69374. doi: 10.1371/journal.pone.0069374

PubMed Abstract | CrossRef Full Text | Google Scholar

84. Josset L, Zeng H, Kelly SM, Tumpey TM, Katze MG. Transcriptomic characterization of the novel avian-origin influenza A (H7N9) virus: specific host response and responses intermediate between avian (H5N1 and H7N7) and human (H3N2) viruses and implications for treatment options. MBio. (2014) 5:e01102–13. doi: 10.1128/mBio.01102-13

PubMed Abstract | CrossRef Full Text | Google Scholar

85. Zhu Z, Lian X, Su X, Wu W, Marraro GA, Zeng Y. From SARS and MERS to COVID-19: a brief summary and comparison of severe acute respiratory infections caused by three highly pathogenic human coronaviruses. Respir Res. (2020) 21:224. doi: 10.1186/s12931-020-01479-w

PubMed Abstract | CrossRef Full Text | Google Scholar

86. Evans PC, Rainger GE, Mason JC, Guzik TJ, Osto E, Stamataki Z, et al. Endothelial dysfunction in COVID-19: a position paper of the ESC Working Group for Atherosclerosis and Vascular Biology, and the ESC Council of Basic Cardiovascular Science. Cardiovasc Res. (2020) 116:2177–84. doi: 10.1093/cvr/cvaa230

PubMed Abstract | CrossRef Full Text | Google Scholar

87. Gustafson D, Raju S, Wu R, Ching C, Veitch S, Rathnakumar K, et al. Overcoming barriers: the endothelium as a linchpin of coronavirus disease 2019 pathogenesis? Arterioscler Thromb Vasc Biol. (2020) 40:1818–29. doi: 10.1161/ATVBAHA.120.314558

PubMed Abstract | CrossRef Full Text | Google Scholar

88. Cho S. CD36 as a therapeutic target for endothelial dysfunction in stroke. Curr Pharm Des. (2012) 18:3721–30. doi: 10.2174/138161212802002760

PubMed Abstract | CrossRef Full Text | Google Scholar

89. Silverstein RL. Oxidized Lipid Uptake by Scavenger Receptor CD36 (Cluster of Differentiation 36) modulates endothelial surface properties and may contribute to atherogenesis. Arterioscler Thromb Vasc Biol. (2018) 38:4–5. doi: 10.1161/ATVBAHA.117.310351

PubMed Abstract | CrossRef Full Text | Google Scholar

90. Yin Y, Pastrana JL Li X, Huang X, Mallilankaraman K, Choi ET, et al. Inflammasomes: sensors of metabolic stresses for vascular inflammation. Front Biosci. (2013) 18:638–49. doi: 10.2741/4127

PubMed Abstract | CrossRef Full Text | Google Scholar

91. Wang X, Li YF, Nanayakkara G, Shao Y, Liang B, Cole L, et al. Lysophospholipid receptors, as novel conditional danger receptors and homeostatic receptors modulate inflammation-novel paradigm and therapeutic potential. J Cardiovasc Transl Res. (2016) 9:343–59. doi: 10.1007/s12265-016-9700-6

PubMed Abstract | CrossRef Full Text | Google Scholar

92. Shao Y, Nanayakkara G, Cheng J, Cueto R, Yang WY, Park JY, et al. Lysophospholipids and their receptors serve as conditional DAMPs and DAMP receptors in tissue oxidative and inflammatory injury. Antioxid Redox Signal. (2018) 28:973–86. doi: 10.1089/ars.2017.7069

PubMed Abstract | CrossRef Full Text | Google Scholar

93. Chandel NS. Mitochondria as signaling organelles. BMC Biol. (2014) 12:34. doi: 10.1186/1741-7007-12-34

CrossRef Full Text | Google Scholar

94. Zhou W, He Q, Zhang C, He X, Cui Z, Liu F, et al. BLOS2 negatively regulates Notch signaling during neural and hematopoietic stem and progenitor cell development. Elife. (2016) 5:e18108. doi: 10.7554/eLife.18108

PubMed Abstract | CrossRef Full Text | Google Scholar

95. Poulsen LC, Edelmann RJ, Kruger S, Dieguez-Hurtado R, Shah A, Stav-Noraas TE, et al. Inhibition of endothelial NOTCH1 signaling attenuates inflammation by reducing cytokine-mediated histone acetylation at inflammatory enhancers. Arterioscler Thromb Vasc Biol. (2018) 38:854–69. doi: 10.1161/ATVBAHA.117.310388

PubMed Abstract | CrossRef Full Text | Google Scholar

96. Briot A, Civelek M, Seki A, Hoi K, Mack JJ, Lee SD, et al. Endothelial NOTCH1 is suppressed by circulating lipids and antagonizes inflammation during atherosclerosis. J Exp Med. (2015) 212:2147–63. doi: 10.1084/jem.20150603

PubMed Abstract | CrossRef Full Text | Google Scholar

97. Holliday CJ, Ankeny RF, Jo H, Nerem RM. Discovery of shear- and side-specific mRNAs and miRNAs in human aortic valvular endothelial cells. Am J Physiol Heart Circ Physiol. (2011) 301:H856–67. doi: 10.1152/ajpheart.00117.2011

PubMed Abstract | CrossRef Full Text | Google Scholar

98. Erbilgin A, Siemers N, Kayne P, Yang WP, Berliner J, Lusis AJ. Gene expression analyses of mouse aortic endothelium in response to atherogenic stimuli. Arterioscler Thromb Vasc Biol. (2013) 33:2509–17. doi: 10.1161/ATVBAHA.113.301989

PubMed Abstract | CrossRef Full Text | Google Scholar

99. Liu B, Tai Y, Caceres AI, Achanta S, Balakrishna S, Shao X, et al. Oxidized phospholipid OxPAPC activates TRPA1 and contributes to chronic inflammatory pain in mice. PLoS ONE. (2016) 11:e0165200. doi: 10.1371/journal.pone.0165200

PubMed Abstract | CrossRef Full Text | Google Scholar

100. Nallar SC, Kalvakolanu DV. Interferons, signal transduction pathways, and the central nervous system. J Interferon Cytokine Res. (2014) 34:559–76. doi: 10.1089/jir.2014.0021

PubMed Abstract | CrossRef Full Text | Google Scholar

101. Souilhol C, Serbanovic-Canic J, Fragiadaki M, Chico TJ, Ridger V, Roddie H, et al. Endothelial responses to shear stress in atherosclerosis: a novel role for developmental genes. Nat Rev Cardiol. (2020) 17:52–63. doi: 10.1038/s41569-019-0239-5

PubMed Abstract | CrossRef Full Text | Google Scholar

102. Wang N, Yuan J, Karim MR, Zhong P, Sun YP, Zhang HY, et al. Effects of mitophagy on regulatory T cell function in patients with myasthenia gravis. Front Neurol. (2020) 11:238. doi: 10.3389/fneur.2020.00238

PubMed Abstract | CrossRef Full Text | Google Scholar

103. Procaccini C, Matarese G. Where mitochondria meet autoimmunity: the Treg Cell link. Cell Metab. (2020) 32:507–9. doi: 10.1016/j.cmet.2020.08.006

PubMed Abstract | CrossRef Full Text | Google Scholar

104. Alissafi T, Kalafati L, Lazari M, Filia A, Kloukina I, Manifava M, et al. Mitochondrial oxidative damage underlies regulatory T cell defects in autoimmunity. Cell Metab. (2020) 32:591–604 e7. doi: 10.1016/j.cmet.2020.07.001

PubMed Abstract | CrossRef Full Text | Google Scholar

105. Beaslas O, Metso J, Nissila E, Laurila PP, Kaiharju E, Batchu KC, et al. Osbpl8 deficiency in mouse causes an elevation of high-density lipoproteins and gender-specific alterations of lipid metabolism. PLoS ONE. (2013) 8:e58856. doi: 10.1371/journal.pone.0058856

PubMed Abstract | CrossRef Full Text | Google Scholar

106. Yang XF, Yin Y, Wang H. Vascular inflammation and atherogenesis are activated via receptors for PAMPs and suppressed by regulatory T Cells. Drug Discov Today Ther Strateg. (2008) 5:125–42. doi: 10.1016/j.ddstr.2008.11.003

PubMed Abstract | CrossRef Full Text | Google Scholar

107. Yin Y, Yan Y, Jiang X, Mai J, Chen NC, Wang H, et al. Inflammasomes are differentially expressed in cardiovascular and other tissues. Int J Immunopathol Pharmacol. (2009) 22:311–22. doi: 10.1177/039463200902200208

PubMed Abstract | CrossRef Full Text | Google Scholar

108. Shen J, Yin Y, Mai J, Xiong X, Pansuria M, Liu J, et al. Caspase-1 recognizes extended cleavage sites in its natural substrates. Atherosclerosis. (2010) 210:422–9. doi: 10.1016/j.atherosclerosis.2009.12.017

PubMed Abstract | CrossRef Full Text | Google Scholar

109. Wang L, Fu H, Nanayakkara G, Li Y, Shao Y, Johnson C, et al. Novel extracellular and nuclear caspase-1 and inflammasomes propagate inflammation and regulate gene expression: a comprehensive database mining study. J Hematol Oncol. (2016) 9:122. doi: 10.1186/s13045-016-0351-5

PubMed Abstract | CrossRef Full Text | Google Scholar

110. Li YF, Nanayakkara G, Sun Y, Li X, Wang L, Cueto R, et al. Analyses of caspase-1-regulated transcriptomes in various tissues lead to identification of novel IL-1beta-, IL-18- and sirtuin-1-independent pathways. J Hematol Oncol. (2017) 10:40. doi: 10.1186/s13045-017-0406-2

PubMed Abstract | CrossRef Full Text | Google Scholar

111. Vara-Perez M, Felipe-Abrio B, Agostinis P. Mitophagy in cancer: a tale of adaptation. Cells. (2019) 8:493. doi: 10.3390/cells8050493

CrossRef Full Text | Google Scholar

112. Chen HT, Liu H, Mao MJ, Tan Y, Mo XQ, Meng XJ, et al. Crosstalk between autophagy and epithelial-mesenchymal transition and its application in cancer therapy. Mol Cancer. (2019) 18:101. doi: 10.1186/s12943-019-1030-2

PubMed Abstract | CrossRef Full Text | Google Scholar

113. Buratta S, Tancini B, Sagini K, Delo F, Chiaradia E, Urbanelli L, et al. Lysosomal exocytosis, exosome release and secretory autophagy: the autophagic- and endo-lysosomal systems go extracellular. Int J Mol Sci. (2020) 21:2576. doi: 10.3390/ijms21072576

PubMed Abstract | CrossRef Full Text | Google Scholar

114. Palmer CS, Anderson AJ, Stojanovski D. Mitochondrial protein import dysfunction: mitochondrial disease, neurodegenerative disease and cancer. FEBS Lett. (2021) 595:1107–31. doi: 10.1002/1873-3468.14022

PubMed Abstract | CrossRef Full Text | Google Scholar

115. Wiedemann N, Pfanner N. Mitochondrial machineries for protein import and assembly. Annu Rev Biochem. (2017) 86:685–714. doi: 10.1146/annurev-biochem-060815-014352

PubMed Abstract | CrossRef Full Text | Google Scholar

116. Hansen KG, Herrmann JM. Transport of proteins into mitochondria. Protein J. (2019) 38:330–42. doi: 10.1007/s10930-019-09819-6

CrossRef Full Text | Google Scholar

117. Jadiya P, Tomar D. Mitochondrial protein quality control mechanisms. Genes. (2020) 11:563. doi: 10.3390/genes11050563

PubMed Abstract | CrossRef Full Text | Google Scholar

118. Lytovchenko O, Kunji ERS. Expression and putative role of mitochondrial transport proteins in cancer. Biochim Biophys Acta. (2017) 1858:641–54. doi: 10.1016/j.bbabio.2017.03.006

PubMed Abstract | CrossRef Full Text | Google Scholar

119. Kreimendahl S, Rassow J. The mitochondrial outer membrane protein Tom70-mediator in protein traffic, membrane contact sites and innate immunity. Int J Mol Sci. (2020) 21:7262. doi: 10.3390/ijms21197262

PubMed Abstract | CrossRef Full Text | Google Scholar

120. Xia Y, Yeddula N, Leblanc M, Ke E, Zhang Y, Oldfield E, et al. Reduced cell proliferation by IKK2 depletion in a mouse lung-cancer model. Nat Cell Biol. (2015) 17:532. doi: 10.1038/ncb3145

PubMed Abstract | CrossRef Full Text | Google Scholar

121. Thomas SJ, Snowden JA, Zeidler MP, Danson SJ. The role of JAK/STAT signalling in the pathogenesis, prognosis and treatment of solid tumours. Br J Cancer. (2015) 113:365–71. doi: 10.1038/bjc.2015.233

PubMed Abstract | CrossRef Full Text | Google Scholar

122. Katz SF, Lechel A, Obenauf AC, Begus-Nahrmann Y, Kraus JM, Hoffmann EM, et al. Disruption of Trp53 in livers of mice induces formation of carcinomas with bilineal differentiation. Gastroenterology. (2012) 142:1229–39 e3. doi: 10.1053/j.gastro.2012.02.009

PubMed Abstract | CrossRef Full Text | Google Scholar

123. Aoki K, Taketo MM. Adenomatous polyposis coli (APC): a multi-functional tumor suppressor gene. J Cell Sci. (2007) 120:3327–35. doi: 10.1242/jcs.03485

PubMed Abstract | CrossRef Full Text | Google Scholar

124. Reed KR, Meniel VS, Marsh V, Cole A, Sansom OJ, Clarke AR, et al. limited role for p53 in modulating the immediate phenotype of Apc loss in the intestine. BMC Cancer. (2008) 8:162. doi: 10.1186/1471-2407-8-162

PubMed Abstract | CrossRef Full Text | Google Scholar

125. Chalhoub N, Baker SJ PTEN and the PI3-kinase pathway in cancer. Annu Rev Pathol. (2009) 4:127–50. doi: 10.1146/annurev.pathol.4.110807.092311

CrossRef Full Text | Google Scholar

126. Chatterjee N, Pazarentzos E, Mayekar MK, Gui P, Allegakoen DV, Hrustanovic G, et al. Synthetic essentiality of metabolic regulator PDHK1 in PTEN-deficient cells and cancers. Cell Rep. (2019) 28:2317–30 e8. doi: 10.1016/j.celrep.2019.07.063

PubMed Abstract | CrossRef Full Text | Google Scholar

127. Tang Z, Li C, Kang B, Gao G, Li C, Zhang Z, et al. a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. (2017) 45:W98–102. doi: 10.1093/nar/gkx247

PubMed Abstract | CrossRef Full Text | Google Scholar

128. Ng B, Yang F, Huston DP, Yan Y, Yang Y, Xiong Z, et al. Increased noncanonical splicing of autoantigen transcripts provides the structural basis for expression of untolerized epitopes. J Allergy Clin Immunol. (2004) 114:1463–70. doi: 10.1016/j.jaci.2004.09.006

PubMed Abstract | CrossRef Full Text | Google Scholar

129. Zemmour D, Zilionis R, Kiner E, Klein AM, Mathis D, Benoist C. Single-cell gene expression reveals a landscape of regulatory T cell phenotypes shaped by the TCR. Nat Immunol. (2018) 19:291–301. doi: 10.1038/s41590-018-0051-0

PubMed Abstract | CrossRef Full Text | Google Scholar

130. Corkum CP, Ings DP, Burgess C, Karwowska S, Kroll W, Michalak TI. Immune cell subsets and their gene expression profiles from human PBMC isolated by Vacutainer Cell Preparation Tube (CPT) and standard density gradient. BMC Immunol. (2015) 16:48. doi: 10.1186/s12865-015-0113-0

PubMed Abstract | CrossRef Full Text | Google Scholar

131. Li X, Fang P, Yang WY, Chan K, Lavallee M, Xu K, et al. Mitochondrial ROS, uncoupled from ATP synthesis, determine endothelial activation for both physiological recruitment of patrolling cells and pathological recruitment of inflammatory cells. Can J Physiol Pharmacol. (2017) 95:247–52. doi: 10.1139/cjpp-2016-0515

PubMed Abstract | CrossRef Full Text | Google Scholar

132. Yang XF, Weber GF, Cantor H A. novel Bcl-x isoform connected to the T cell receptor regulates apoptosis in T cells. Immunity. (1997) 7:629–39. doi: 10.1016/S1074-7613(00)80384-2

PubMed Abstract | CrossRef Full Text | Google Scholar

133. Yang Y, Yang F, Xiong Z, Yan Y, Wang X, Nishino M, et al. An N-terminal region of translationally controlled tumor protein is required for its antiapoptotic activity. Oncogene. (2005) 24:4778–88. doi: 10.1038/sj.onc.1208666

PubMed Abstract | CrossRef Full Text | Google Scholar

134. Kumar V, Maity S. ER Stress-sensor proteins and ER-mitochondrial crosstalk-signaling beyond (ER) stress response. Biomolecules. (2021) 11:173. doi: 10.3390/biom11020173

PubMed Abstract | CrossRef Full Text | Google Scholar

135. Heaton SM, Borg NA, Dixit VM. Ubiquitin in the activation and attenuation of innate antiviral immunity. J Exp Med. (2016) 213:1–13. doi: 10.1084/jem.20151531

PubMed Abstract | CrossRef Full Text | Google Scholar

136. Glebov OO. Understanding SARS-CoV-2 endocytosis for COVID-19 drug repurposing. FEBS J. (2020) 287:3664–71. doi: 10.1111/febs.15369

PubMed Abstract | CrossRef Full Text | Google Scholar

137. Yang Q, Nanayakkara GK, Drummer C, Sun Y, Johnson C, Cueto R, et al. Low-intensity ultrasound-induced anti-inflammatory effects are mediated by several new mechanisms including gene induction, immunosuppressor cell promotion, and enhancement of exosome biogenesis and docking. Front Physiol. (2017) 8:818. doi: 10.3389/fphys.2017.00818

PubMed Abstract | CrossRef Full Text | Google Scholar