Identification of key candidate genes and pathways in hepatitis B virus-associated acute liver failure by bioinformatical analysis

Multi-omics analysis of a case of congenital microtia reveals aldob and oxidative stress associated with microtia etiology

BACKGROUND

Microtia is reported to be one of the most common congenital craniofacial malformations. Due to the complex etiology and the ethical barrier of embryonic study, the precise mechanisms of microtia remain unclear. Here we report a rare case of microtia with costal chondrodysplasia based on bioinformatics analysis and further verifications on other sporadic microtia patients.

RESULTS

One hundred fourteen deleterious insert and deletion (InDel) and 646 deleterious SNPs were screened out by WES, candidate genes were ranked in descending order according to their relative impact with microtia. Label-free proteomic analysis showed that proteins significantly different between the groups were related with oxidative stress and energy metabolism. By real-time PCR and immunohistochemistry, we further verified the candidate genes between other sporadic microtia and normal ear chondrocytes, which showed threonine aspartase, cadherin-13, aldolase B and adiponectin were significantly upregulated in mRNA levels but were significantly lower in protein levels. ROS detection and mitochondrial membrane potential (∆ Ψ m) detection proved that oxidative stress exists in microtia chondrocytes.

CONCLUSIONS

Our results not only spot new candidate genes by WES and label-free proteomics, but also speculate for the first time that metabolism and oxidative stress may disturb cartilage development and this might become therapeutic targets and potential biomarkers with clinical usefulness in the future.

Identifying the Mechanism of Polygoni Cuspidati Rhizoma et Radix in Treating Acute Liver Failure Based on Network Pharmacology and Molecular Docking

Materials and Methods

The potential bioactive compounds of PCRR and their targets were collected from TCMSP, TCMID, and BATMAN-TCM databases with absorption, distribution, metabolism, and excretion protocols (oral bioavailability ≥30% and drug-likeness ≥0.18). The ALF-related target genes were identified using the GeneCards and OMIM databases. A protein-protein interaction (PPI) network among these targets was constructed using the Cytoscape software to obtain the core targets. The genes associated with ALF were analyzed via Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses to identify the signaling pathways related to the therapeutic effect of PCRR in ALF.

Results

In total, 10 bioactive compounds of PCRR and 200 targets related to them were obtained, and 2913 ALF-related target genes were identified. PPI network analysis pinpointed 15 core targets, namely, TP53, AKT1, JUN, HSP90AA1, MAPK1, RELA, TNF, ESR1, IL6, MYC, MAPK14, FOS, RB1, CDKN1A, and EGFR. GO enrichment and KEGG pathway analyses revealed that the therapeutic mechanisms of PCRR in ALF are related to cell metabolism, oxidative stress, inflammation, and hepatocyte apoptosis.

Conclusion

This is the first study to explore the therapeutic mechanisms of PCRR in ALF via network pharmacology and molecular docking. This study provides a research platform with candidate ALF-related targets of PRCC for the development of therapeutics against ALF.

New Insights into the Role of the Complement System in Human Viral Diseases

The complement system (CS) is part of the human immune system, consisting of more than 30 proteins that play a vital role in the protection against various pathogens and diseases, including viral diseases. Activated via three pathways, the classical pathway (CP), the lectin pathway (LP), and the alternative pathway (AP), the complement system leads to the formation of a membrane attack complex (MAC) that disrupts the membrane of target cells, leading to cell lysis and death. Due to the increasing number of reports on its role in viral diseases, which may have implications for research on severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), this review aims to highlight significant progress in understanding and defining the role of the complement system in four groups of diseases of viral etiology: (1) respiratory diseases; (2) acute liver failure (ALF); (3) disseminated intravascular coagulation (DIC); and (4) vector-borne diseases (VBDs). Some of these diseases already present a serious global health problem, while others are a matter of concern and require the collaboration of relevant national services and scientists with the World Health Organization (WHO) to avoid their spread.

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.

Orosomucoid in liver diseases

In this editorial, the roles of orosomucoid (ORM) in the diagnoses and follow-up assessments of both nonneoplastic diseases and liver tumors are discussed with respect to the publication by Zhu et al presented in the previous issue of World Journal of Gastroenterology (2020; 26(8): 840-817). ORM, or alpha-1 acid glycoprotein (AGP), is an acute-phase protein that constitutes 1% to 3% of plasma proteins in humans and is mainly synthesized in the liver. ORM exists in serum as two variants: ORM1 and ORM2. Although the variants share 89.6% sequence identity and have similar biological properties, ORM1 constitutes the main component of serum ORM. An interesting feature of ORM is that its biological effects differ according to variations in glycosylation patterns. This variable feature makes ORM an attractive target for diagnosing and monitoring many diseases, including those of the liver. Recent findings suggest that a sharp decrease in ORM level is an important marker for HBV-associated acute liver failure (ALF), and ORM1 plays an important role in liver regeneration. In viral hepatitis, increases in both ORM and its fucosylated forms and the correlation of these increases with fibrosis progression suggest that this glycoprotein can be used with other markers as a noninvasive method in the follow-up assessment of diseases. In addition, similar findings regarding the level of the asialylated form of ORM, called asialo-AGP (AsAGP), have been reported in a follow-up assessment of fibrosis in chronic liver disease. An increase in ORM in serum has also been shown to improve hepatocellular carcinoma (HCC) diagnosis performance when combined with other markers. In addition, determination of the ORM level has been useful in the diagnosis of HCC with AFP concentrations less than 500 ng/mL. For monitoring patients with AFP-negative HCC, a unique trifucosylated tetra-antennary glycan of ORM may also be used as a new potential marker. The fact that there are very few studies investigating the expression of this glycoprotein and its variants in liver tissues constitutes a potential limitation, especially in terms of revealing all the effects of ORM on carcinogenesis and tumor behavior. Current findings indicate that ORM2 expression is decreased in tumors, and this is related to the aggressive course of the disease. Parallel to this finding, in HCC cell lines, ORM2 decreases HCC cell migration and invasion, supporting reports of its tumor suppressor role. In conclusion, the levels of ORM and its different glycosylated variants are promising additional biomarkers for identifying ALF, for monitoring fibrosis in viral hepatitis, and for diagnosing early HCC. Although there is evidence that the loss of ORM2 expression in HCC is associated with poor prognosis, further studies are needed to support these findings. Additionally, investigations of ORM expression in borderline dysplastic nodules and hepatocellular adenomas, which pose diagnostic problems in the differential diagnosis of HCC, especially in biopsy samples, may shed light on whether ORM can be used in histopathological differential diagnosis.

Potential Networks Regulated by MSCs in Acute-On-Chronic Liver Failure: Exosomal miRNAs and Intracellular Target Genes

Acute-on-chronic liver failure (ACLF) is a severe syndrome associated with high mortality. Alterations in the liver microenvironment are one of the vital causes of immune damage and liver dysfunction. Human bone marrow mesenchymal stem cells (hBMSCs) have been reported to alleviate liver injury via exosome-mediated signaling; of note, miRNAs are one of the most important cargoes in exosomes. Importantly, the miRNAs within exosomes in the hepatic microenvironment may mediate the mesenchymal stem cell (MSC)-derived regulation of liver function. This study investigated the hepatocyte exosomal miRNAs which are regulated by MSCs and the target genes which have potential in the treatment of liver failure. Briefly, ACLF was induced in mice using carbon tetrachloride and primary hepatocytes were isolated and co-cultured (or not) with MSCs under serum-free conditions. Exosomes were then collected, and the expression of exosomal miRNAs was assessed using next-generation sequencing; a comparison was performed between liver cells from healthy versus ACLF animals. Additionally, to identify the intracellular targets of exosomal miRNAs in humans, we focused on previously published data, i.e., microarray data and mass spectrometry data in liver samples from ACLF patients. The biological functions and signaling pathways associated with differentially expressed genes were predicted using gene ontology and Kyoto Encyclopedia of Genes and Genomics enrichment analyses; hub genes were also screened based on pathway analysis and the prediction of protein-protein interaction networks. Finally, we constructed the hub gene-miRNA network and performed correlation analysis and qPCR validation. Importantly, our data revealed that MSCs could regulate the miRNA content within exosomes in the hepatic microenvironment. MiR-20a-5p was down-regulated in ACLF hepatocytes and their exosomes, while the levels of chemokine C-X-C Motif Chemokine Ligand 8 (CXCL8; interleukin 8) were increased in hepatocytes. Importantly, co-culture with hBMSCs resulted in up-regulated expression of miR-20a-5p in exosomes and hepatocytes, and down-regulated expression of CXCL8 in hepatocytes. Altogether, our data suggest that the exosomal miR-20a-5p/intracellular CXCL8 axis may play an important role in the reduction of liver inflammation in ACLF in the context of MSC-based therapies and highlights CXCL8 as a potential target for alleviating liver injury.

Long non-coding RNA Gm15441 attenuates hepatic inflammasome activation in response to PPARA agonism and fasting

Exploring the molecular mechanisms that prevent inflammation during caloric restriction may yield promising therapeutic targets. During fasting, activation of the nuclear receptor peroxisome proliferator-activated receptor α (PPARα) promotes the utilization of lipids as an energy source. Herein, we show that ligand activation of PPARα directly upregulates the long non-coding RNA gene Gm15441 through PPARα binding sites within its promoter. Gm15441 expression suppresses its antisense transcript, encoding thioredoxin interacting protein (TXNIP). This, in turn, decreases TXNIP-stimulated NLR family pyrin domain containing 3 (NLRP3) inflammasome activation, caspase-1 (CASP1) cleavage, and proinflammatory interleukin 1β (IL1B) maturation. Gm15441-null mice were developed and shown to be more susceptible to NLRP3 inflammasome activation and to exhibit elevated CASP1 and IL1B cleavage in response to PPARα agonism and fasting. These findings provide evidence for a mechanism by which PPARα attenuates hepatic inflammasome activation in response to metabolic stress through induction of lncRNA Gm15441.

Identification of Hub Genes and Potential Molecular Mechanisms in Patients with HBV-Associated Acute Liver Failure

Hepatitis B virus (HBV) infection is a major cause of acute liver failure (ALF) in China, and mortality rates are high among patients who do not receive a matched liver transplant. This study aimed to determine potential mechanisms involved in HBV-ALF pathogenesis. Gene expression profiles under access numbers GSE38941 and GSE14668 were downloaded from the Gene Expression Omnibus database, including cohorts of HBV-ALF liver tissue and normal samples. Differentially expressed genes (DEGs) with false discovery rates (FDR) <0.05 and |log₂(fold change)| >1 as thresholds were screened using the Limma package. Gene modules associated with stable disease were mined using weighed gene co-expression network analysis (WGCNA). A co-expression network was constructed and DEGs were analyzed using gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. A gene-based network was constructed to explore major factors associated with disease progression. We identified 2238 overlapping DEGs as crucial gene cohorts in ALF development. Based on a WGCNA algorithm, 10 modules (modules 1-10) were obtained that ranged from 75 to 1078 genes per module. Cyclin-dependent kinase 1 (CDK1), cyclin B1 (CCNB1), and cell-division cycle protein 20 (CDC20) hub genes were screened using the co-expression network. Furthermore, 17 GO terms and 6 KEGG pathways were identified, such as cell division, immune response process, and antigen processing and presentation. Two overlapping signaling pathways that are crucial factors in HBV-ALF were screened using the Comprehensive Toxicogenomics Database (CTD). Several candidate genes including HLA-E, B2M, HLA-DPA1, and SYK were associated with HBV-ALF progression. Natural killer cell-mediated cytotoxicity and antigen presentation contributed to the progression of HBV-ALF. The HLA-E, B2M, HLA-DPA1, and SYK genes play critical roles in the pathogenesis and development of HBV-ALF.

Identification of TYMS as a promoting factor of retroperitoneal liposarcoma progression: Bioinformatics analysis and biological evidence

Retroperitoneal liposarcoma (RLPS) is one of the most common types of retroperitoneal sarcomas, and has a high recurrence rate. There is an urgent need to further explore its pathogenesis and develop more effective treatment strategies. The aim of the present study was to identify potential driver genes of RLPS through bioinformatics analysis and molecular biology to elucidate potential targets that are suitable for further analysis for the treatment of RLPS. Differentially expressed genes (DEGs) between liposarcoma and normal fatty (NF) tissues were identified based on microarray data through bioinformatics analysis, and thymidylate synthase (TYMS) was selected from the DEGs, based on high content screening (HCS). TYMS expression was evaluated in RLPS tumor tissues and cell lines. A total of 21 RLPS tissues and 10 NF frozen tissues were used for reverse transcription-quantitative PCR, and 47 RLPS formalin-fixed specimens were used for immunohistochemical analysis. The effect of TYMS downregulation on cell proliferation, apoptosis, cell cycle progression, and cell migration and invasion were evaluated using lentivirus-mediated short hairpin RNA. The underlying mechanisms of TYMS in RLPS were examined by protein microarray and verified by western blotting. A total of 855 DEGs were identified. TYMS knockdown had the most notable effect on the proliferative capacity of RLPS cells according to the HCS results. TYMS mRNA expression levels were higher in RLPS tissues compared with NF tissues (P<0.001). TYMS expression was higher in high-grade RLPS tissues compared with low-grade RLPS tissues (P=0.003). The patients with positive TYMS expression had a worse overall survival (OS) and disease-free survival (DFS) compared with the patients with negative TYMS expression (OS, P=0.024; DFS, P=0.030). The knockdown of TYMS reduced proliferation, promoted apoptosis, facilitated cell cycle progression from G₁ to S phase, and reduced cell migration and invasion of RLPS cells. Protein microarray analysis and western blotting showed that the Janus Kinase/Signal transducers and activators of transcription pathway was downregulated following TYMS knockdown. In conclusion, TYMS expression is upregulated in RLPS tissues, and downregulation of TYMS reduces RLPS progression.

Downregulation of orosomucoid 2 acts as a prognostic factor associated with cancer-promoting pathways in liver cancer

BACKGROUND

Liver cancer has a high mortality and morbidity rate throughout the world. In clinical practice, the prognosis of liver cancer patients is poor, and the complex reasons contribute to treatment failures, including fibrosis, hepatitis viral infection, drug resistance and metastasis. Thus, screening novel prognostic biomarkers is of great importance for guiding liver cancer therapy. Orosomucoid genes (ORMs) encode acute phase plasma proteins, including orosomucoid 1 (ORM1) and ORM2. Previous studies showed their upregulation upon inflammation, but the specific function of ORMs has not yet been determined, especially in the development of liver cancer.

AIM

To determine the expression of ORMs and their potential function in liver cancer.

METHODS

Analysis of the expression of ORMs in different human tissues was performed on data from the HPA RNA-seq normal tissues project. The expression ratio of ORMs was determined using the HCCDB database, including the ratio between liver cancer and other cancers, normal liver and other normal tissues, liver cancer and adjacent normal liver tissues. Analysis of ORM expression in different cancer types was performed using The Cancer Genome Atlas and TIMER database. The expression of ORMs in liver tumor tissues and adjacent normal tissues were further confirmed using Gene Expression Omnibus data, including GSE36376 and GSE14520. The 10-year overall survival (OS), progression-free survival (PFS) and relapse-free survival (RFS) rates between high and low ORM expression groups in liver cancer patients were determined using the Kaplan-Meier plotter tool. Gene Set Enrichment Analysis (GSEA) was employed to explore the ORM2-associated signaling network. Correlations between ORM2 expression and tumor purity or the infiltration level of macrophages in liver tumor tissues were determined using the TIMER database. The correlation between ORM2 gene levels, tumor-associated macrophage (TAM) markers (including CD68 and TGFβ1) and T cell immunosuppression (including CTLA4 and PD-1) in liver tumor tissues and liver GTEx was determined using the GEPIA database.

RESULTS

ORM1 and ORM2 were highly expressed in normal liver and liver tumor tissues. ORM1 and ORM2 expression was significantly decreased in liver tumor tissues compared with adjacent normal tissues, and similar results were also noted in cholangiocarcinoma, esophageal carcinoma, and lung squamous cell carcinoma. Further analysis of the Gene Expression Omnibus Database also confirmed the downregulation of ORM1 and ORM2 in liver tumors. Survival analysis showed that the high ORM2 group had better survival rates in OS, PFS and RFS. ORM1 only represented better performance in PFS, but not in OS or RFS. GSEA analysis of ORM2 from The Cancer Genome Atlas liver cancer data identified that ORM2 positively associated with the G2/M checkpoint, E2F target signaling, as well as Wnt/β-catenin and Hedgehog signaling. Moreover, apoptosis, IFN-α responses, IFN-γ responses and humoral immune responses were upregulated in the ORM2 high group. ORM2 expression was negatively correlated with the macrophage infiltration level, CD68, TGFβ1, CTLA4 and PD-1 levels.

CONCLUSION

The results showed that ORM1 and ORM2 were highly expressed specifically in liver tissues, whereas ORM1 and ORM2 were downregulated in liver tumor tissues. ORM2 is a better prognostic factor for liver cancer. Furthermore, ORM2 is closely associated with cancer-promoting pathways.

Bioinformatics analysis on multiple Gene Expression Omnibus datasets of the hepatitis B virus infection and its response to the interferon-alpha therapy

Background

Hepatitis B virus (HBV) infection is a global health problem and interferon-alpha (IFN-α) is one of the effective therapies. However, little is known about the genetic background of the HBV infection or the genetic determinants of the IFN-α treatment response. Thus, we aim to explore the possible molecular mechanisms of HBV infection and its response to the IFN-α therapy with a comprehensive bioinformatics analysis.

Methods

The Gene Expression Omnibus datasets (GSE83148, GSE84044 and GSE66698) were collected and the differentially expressed genes (DEGs), key biological processes and intersecting pathways were analyzed. The expression of the co-expressed DEGs in the clinical samples was verified by quantitative real time polymerase chain reaction (qRT-PCR).

Results

Analysis of all the 3 datasets revealed that there were eight up-regulated and one down-regulated co-expressed DEGs following the HBV infection and after IFN-α treatment. In clinical samples, the mRNA level of HKDC1, EPCAM, GSN, ZWINT and PLD3 were significantly increased, while, the mRNA level of PLEKHA2 was significantly decreased in HBV infected liver tissues compared to normal liver tissues. PI3K-Akt signaling pathway, focal adhesion, HTLV-I infection, cytokine-cytokine receptor interaction, metabolic pathways, NF-κB signaling pathway were important pathways associated with the HBV infection and the response of IFN-α treatment.

Conclusions

The co-expressed genes, common biological processes and intersecting pathways identified in the study might play an important role in HBV infection and response of IFN-α treatment. The dysregulated genes may act as novel biomarkers and therapeutic targets for HBV.

Complement System as a Target for Therapies to Control Liver Regeneration/Damage in Acute Liver Failure Induced by Viral Hepatitis

The complement system plays an important role in innate immunity inducing liver diseases as well as signaling immune cell activation in local inflammation regulating immunomodulatory effects such as liver damage and/or liver regeneration. Our aim is to evaluate the role of complement components in acute liver failure (ALF) caused by viral hepatitis, involving virus-induced ALF in human subjects using peripheral blood, samples of liver tissues, and ex vivo assays. Our findings displayed low levels of C3a in plasma samples with high frequency of C3a, C5a, and C5b/9 deposition in liver parenchyma. Meanwhile, laboratory assays using HepG2 (hepatocyte cell line) showed susceptibility to plasma samples from ALF patients impairing in vitro cell proliferation and an increase in apoptotic events submitting plasma samples to heat inactivation. In summary, our data suggest that the complement system may be involved in liver dysfunction in viral-induced acute liver failure cases using ex vivo assays. In extension to our findings, we provide insights into future studies using animal models for viral-induced ALF, as well as other associated soluble components, which need further investigation.