In-depth systems biological evaluation of bovine alveolar macrophages suggests novel insights into molecular mechanisms underlying Mycobacterium bovis infection.

Objective: Bovine tuberculosis (bTB) is a chronic respiratory infectious disease of domestic livestock caused by intracellular Mycobacterium bovis infection, which causes ~$3 billion in annual losses to global agriculture. Providing novel tools for bTB managements requires a comprehensive understanding of the molecular regulatory mechanisms underlying the M. bovis infection. Nevertheless, a combination of different bioinformatics and systems biology methods was used in this study in order to clearly understand the molecular regulatory mechanisms of bTB, especially the immunomodulatory mechanisms of M. bovis infection. Methods: RNA-seq data were retrieved and processed from 78 (39 non-infected control vs. 39 M. bovis-infected samples) bovine alveolar macrophages (bAMs). Next, weighted gene co-expression network analysis (WGCNA) was performed to identify the co-expression modules in non-infected control bAMs as reference set. The WGCNA module preservation approach was then used to identify non-preserved modules between non-infected controls and M. bovis-infected samples (test set). Additionally, functional enrichment analysis was used to investigate the biological behavior of the non-preserved modules and to identify bTB-specific non-preserved modules. Co-expressed hub genes were identified based on module membership (MM) criteria of WGCNA in the non-preserved modules and then integrated with protein-protein interaction (PPI) networks to identify co-expressed hub genes/transcription factors (TFs) with the highest maximal clique centrality (MCC) score (hub-central genes). Results: As result, WGCNA analysis led to the identification of 21 modules in the non-infected control bAMs (reference set), among which the topological properties of 14 modules were altered in the M. bovis-infected bAMs (test set). Interestingly, 7 of the 14 non-preserved modules were directly related to the molecular mechanisms underlying the host immune response, immunosuppressive mechanisms of M. bovis, and bTB development. Moreover, among the co-expressed hub genes and TFs of the bTB-specific non-preserved modules, 260 genes/TFs had double centrality in both co-expression and PPI networks and played a crucial role in bAMs-M. bovis interactions. Some of these hub-central genes/TFs, including PSMC4, SRC, BCL2L1, VPS11, MDM2, IRF1, CDKN1A, NLRP3, TLR2, MMP9, ZAP70, LCK, TNF, CCL4, MMP1, CTLA4, ITK, IL6, IL1A, IL1B, CCL20, CD3E, NFKB1, EDN1, STAT1, TIMP1, PTGS2, TNFAIP3, BIRC3, MAPK8, VEGFA, VPS18, ICAM1, TBK1, CTSS, IL10, ACAA1, VPS33B, and HIF1A, had potential targets for inducing immunomodulatory mechanisms by M. bovis to evade the host defense response. Conclusion: The present study provides an in-depth insight into the molecular regulatory mechanisms behind M. bovis infection through biological investigation of the candidate non-preserved modules directly related to bTB development. Furthermore, several hub-central genes/TFs were identified that were significant in determining the fate of M. bovis infection and could be promising targets for developing novel anti-bTB therapies and diagnosis strategies.

Integrated Network Analysis to Identify Key Modules and Potential Hub Genes Involved in Bovine Respiratory Disease: A Systems Biology Approach.

Background: Bovine respiratory disease (BRD) is the most common disease in the beef and dairy cattle industry. BRD is a multifactorial disease resulting from the interaction between environmental stressors and infectious agents. However, the molecular mechanisms underlying BRD are not fully understood yet. Therefore, this study aimed to use a systems biology approach to systematically evaluate this disorder to better understand the molecular mechanisms responsible for BRD. Methods: Previously published RNA-seq data from whole blood of 18 healthy and 25 BRD samples were downloaded from the Gene Expression Omnibus (GEO) and then analyzed. Next, two distinct methods of weighted gene coexpression network analysis (WGCNA), i.e., module-trait relationships (MTRs) and module preservation (MP) analysis were used to identify significant highly correlated modules with clinical traits of BRD and non-preserved modules between healthy and BRD samples, respectively. After identifying respective modules by the two mentioned methods of WGCNA, functional enrichment analysis was performed to extract the modules that are biologically related to BRD. Gene coexpression networks based on the hub genes from the candidate modules were then integrated with protein-protein interaction (PPI) networks to identify hub-hub genes and potential transcription factors (TFs). Results: Four significant highly correlated modules with clinical traits of BRD as well as 29 non-preserved modules were identified by MTRs and MP methods, respectively. Among them, two significant highly correlated modules (identified by MTRs) and six nonpreserved modules (identified by MP) were biologically associated with immune response, pulmonary inflammation, and pathogenesis of BRD. After aggregation of gene coexpression networks based on the hub genes with PPI networks, a total of 307 hub-hub genes were identified in the eight candidate modules. Interestingly, most of these hub-hub genes were reported to play an important role in the immune response and BRD pathogenesis. Among the eight candidate modules, the turquoise (identified by MTRs) and purple (identified by MP) modules were highly biologically enriched in BRD. Moreover, STAT1, STAT2, STAT3, IRF7, and IRF9 TFs were suggested to play an important role in the immune system during BRD by regulating the coexpressed genes of these modules. Additionally, a gene set containing several hub-hub genes was identified in the eight candidate modules, such as TLR2, TLR4, IL10, SOCS3, GZMB, ANXA1, ANXA5, PTEN, SGK1, IFI6, ISG15, MX1, MX2, OAS2, IFIH1, DDX58, DHX58, RSAD2, IFI44, IFI44L, EIF2AK2, ISG20, IFIT5, IFITM3, OAS1Y, HERC5, and PRF1, which are potentially critical during infection with agents of bovine respiratory disease complex (BRDC). Conclusion: This study not only helps us to better understand the molecular mechanisms responsible for BRD but also suggested eight candidate modules along with several promising hub-hub genes as diagnosis biomarkers and therapeutic targets for BRD.

Differential Co-Expression Network Analysis Reveals Key Hub-High Traffic Genes as Potential Therapeutic Targets for COVID-19 Pandemic.

Background: The recent emergence of COVID-19, rapid worldwide spread, and incomplete knowledge of molecular mechanisms underlying SARS-CoV-2 infection have limited development of therapeutic strategies. Our objective was to systematically investigate molecular regulatory mechanisms of COVID-19, using a combination of high throughput RNA-sequencing-based transcriptomics and systems biology approaches. Methods: RNA-Seq data from peripheral blood mononuclear cells (PBMCs) of healthy persons, mild and severe 17 COVID-19 patients were analyzed to generate a gene expression matrix. Weighted gene co-expression network analysis (WGCNA) was used to identify co-expression modules in healthy samples as a reference set. For differential co-expression network analysis, module preservation and module-trait relationships approaches were used to identify key modules. Then, protein-protein interaction (PPI) networks, based on co-expressed hub genes, were constructed to identify hub genes/TFs with the highest information transfer (hub-high traffic genes) within candidate modules. Results: Based on differential co-expression network analysis, connectivity patterns and network density, 72% (15 of 21) of modules identified in healthy samples were altered by SARS-CoV-2 infection. Therefore, SARS-CoV-2 caused systemic perturbations in host biological gene networks. In functional enrichment analysis, among 15 non-preserved modules and two significant highly-correlated modules (identified by MTRs), 9 modules were directly related to the host immune response and COVID-19 immunopathogenesis. Intriguingly, systemic investigation of SARS-CoV-2 infection identified signaling pathways and key genes/proteins associated with COVID-19's main hallmarks, e.g., cytokine storm, respiratory distress syndrome (ARDS), acute lung injury (ALI), lymphopenia, coagulation disorders, thrombosis, and pregnancy complications, as well as comorbidities associated with COVID-19, e.g., asthma, diabetic complications, cardiovascular diseases (CVDs), liver disorders and acute kidney injury (AKI). Topological analysis with betweenness centrality (BC) identified 290 hub-high traffic genes, central in both co-expression and PPI networks. We also identified several transcriptional regulatory factors, including NFKB1, HIF1A, AHR, and TP53, with important immunoregulatory roles in SARS-CoV-2 infection. Moreover, several hub-high traffic genes, including IL6, IL1B, IL10, TNF, SOCS1, SOCS3, ICAM1, PTEN, RHOA, GDI2, SUMO1, CASP1, IRAK3, HSPA5, ADRB2, PRF1, GZMB, OASL, CCL5, HSP90AA1, HSPD1, IFNG, MAPK1, RAB5A, and TNFRSF1A had the highest rates of information transfer in 9 candidate modules and central roles in COVID-19 immunopathogenesis. Conclusion: This study provides comprehensive information on molecular mechanisms of SARS-CoV-2-host interactions and identifies several hub-high traffic genes as promising therapeutic targets for the COVID-19 pandemic.