Integrated WGCNA and PPI Network to Screen Hub Genes Signatures for Infantile Hemangioma

Comparative Proteomics Identified EXOSC1 as a Target Protein of Anticancer Peptide LVTX-8 in Nasopharyngeal Carcinoma Cells

Nasopharyngeal carcinoma (NPC) is a prevalent malignancy that usually occurs among the nose and throat. Due to mild initial symptoms, most patients are diagnosed in the late stage, and the recurrence rate of tumors is high, resulting in many deaths every year. Traditional radiotherapy and chemotherapy are prone to causing drug resistance and significant side effects. Therefore, searching for new bioactive drugs including anticancer peptides is necessary and urgent. LVTX-8 is a peptide toxin synthesized from the cDNA library of the spider Lycosa vittata, which is consisting of 25 amino acids. In this study, a series of in vitro cell experiments such as cell toxicity, colony formation, and cell migration assays were performed to exam the anticancer activity of LVTX-8 in NPC cells (5-8F and CNE-2). The results suggested that LVTX-8 significantly inhibited cell proliferation and migration of NPC cells. To find the potential molecular targets for the anticancer capability of LVTX-8, high-throughput proteomic and bioinformatics analysis were conducted on NPC cells. The results identified EXOSC1 as a potential target protein with significantly differential expression levels under LVTX-8+/LVTX-8- conditions. The results in this research indicate that spider peptide toxin LVTX-8 exhibits significant anticancer activity in NPC, and EXOSC1 may serve as a target protein for its anticancer activity. These findings provide a reference for the development of new therapeutic drugs for NPC and offer new ideas for the discovery of biomarkers related to NPC diagnosis. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium (https://proteomecentral.proteomexchange.org) via the iProX partner repository with the data set identifier PXD050542.

Fucoidan suppresses proliferation and epithelial-mesenchymal transition process via Wnt/β-catenin signalling in hemangioma

Hemangioma is a common benign tumour that usually occurs on the skin of the head and neck, particularly among infants. The current clinical treatment against hemangioma is surgery excision, however, application of drug is a safer and more economical therapy for children suffering from hemangioma. As a natural sulfated polysaccharide rich in brown algae, fucoidan is widely recognized for anti-tumour bioactivity and dosage safety in humans. This study aims to demonstrate the anti-tumour effect and underlying mechanism of fucoidan against hemangioma in vivo and in vitro. We investigated the effects of fucoidan by culturing hemangioma cells in vitro and treating BALB/c mice bearing with hemangioma. At first, we measured the cell proliferation and migration ability through in vitro experiments. Then, we tested the expression of epithelial-mesenchymal transition (EMT) and Wnt/β-catenin pathway-related biomarkers by western blot and qPCR. Furthermore, we applied β-catenin-specific inhibitor, XAV939, to determine whether fucoidan suppressed EMT via the Wnt/β-catenin pathway in hemangioma cells. In vivo experiments, we applied oral gavage of fucoidan to treat EOMA-bearing mice, along with evaluating the safety and efficacy of fucoidan. We found that fucoidan remarkably inhibits the proliferation and EMT ability of hemangioma cells, which is dependent on the Wnt/β-catenin pathway. These results suggest that fucoidan exhibits tumour inhibitory effect on aggressive hemangioma via regulating the Wnt/β-catenin signalling pathway both in vitro and in vivo, providing a new potent drug candidate for treating hemangioma.

Global research trends of infantile hemangioma: A bibliometric and visualization analysis from 2000 to 2022

Background

Infantile hemangioma (IH) has received global attention, resulting in a significant volume of literature. However, there is a lack of bibliometric analyses specifically focusing on IH publications. This study aims to fill this gap by conducting a comprehensive analysis of IH publications, investigating their characteristics, contribution distribution, and developmental trends. By enhancing our understanding of IH and identifying potential research topics and collaborators, this study will contribute to the advancement of the field.

Methods

A total of 4333 articles and reviews on IH were collected from the Web of Science (WoS) database, spanning the years 2000-2022. The study encompassed a comprehensive analysis of IH publications, evaluating their quantity and quality. Additionally, we profiled publishing groups based on country, institution, author publication records, and collaboration networks. Lastly, we identified and summarized the prominent research topics.

Results

Annual publications on IH have increased over the past 20 years. The United States has the highest number of publications and the highest total number of citations. Pediatric Dermatology was the most influential journal in the IH field. The citation analysis indicated that the articles published by Léauté-Labrèze in 2008 had the highest number of citations. The articles published by North PE in 2000 and Boye E in 2001 laid a certain research foundation for this field. Concerning institutions, most of the cooperative relationships were established in the same country/region. The United States has the largest number of scientific research institutions and IH researchers, leading most of the cross-country collaboration. The University of California, San Francisco, Medical College of Wisconsin, Harvard University, and Shanghai Jiaotong University were the research centers that published the most IH-related research. Frieden IJ, Mulliken JB, and Drolet BA were the top three most influential authors. Frieden IJ, Garzon MC, and Mulliken JB were the top three authors with the most cited frequency. In addition, keywords and keyword co-occurrence networks prompted that the pathological mechanism of IH, clinical analysis, and other vascular anomalies are research hotspots. Analysis of trending topics suggests that research on IH has evolved from treatment-focused studies towards investigations of other vascular diseases and a series of clinical case studies. Currently, clinical case studies receive the most attention in the field.

Conclusions

This comprehensive bibliometric study provides a thorough analysis of post-2000 publications in the field of IH, offering insights into current research trends for the first time. The findings suggest that future investigations will continue to prioritize understanding IH mechanisms, treatment approaches, and treatment evaluation. Furthermore, the exploration of other vascular diseases and the inclusion of clinical case studies are expected to contribute to advancements in IH clinical practice. By identifying potential collaborators, partner institutions, and new research avenues, this study offers valuable guidance for future in-depth research on IH.

Genome-wide DNA methylation and transcription analysis in psoriatic epidermis

Aim: To identify DNA methylation and transcription biomarkers in the psoriatic epidermis. Materials & methods: Gene transcription and DNA methylation datasets of psoriatic epidermal tissue were obtained from the Gene Expression Omnibus. Machine learning algorithm analysis and weighted gene coexpression network analysis were carried out to screen hub genes. Results: Differentially methylated and expressed genes were identified in the psoriatic epidermis. Six hub genes were selected - GZMB, CRIP1, S100A12, ISG15, CRABP2 and VNN1 - whose transcript levels showed a significant correlation with Psoriasis Area and Severity Index scores and immune infiltration. Conclusion: Psoriatic epidermis is primarily in a hypermethylated status. Epidermis-specific hub differentially methylated and expressed genes are potential biomarkers to help judge the condition of psoriasis.

Identification of molecular subtypes of ischaemic stroke based on immune-related genes and weighted co-expression network analysis

Immune system has been reported to play a key role in the development of ischaemic stroke (IS). Nevertheless, its exact immune-related mechanism has not yet been fully revealed. Gene expression data of IS and healthy control samples was downloaded from Gene Expression Omnibus database and differentially expressed genes (DEGs) was obtained. Immune-related genes (IRGs) data was downloaded from the ImmPort database. The molecular subtypes of IS were identified based on IRGs and weighted co-expression network analysis (WGCNA). 827 DEGs and 1142 IRGs were obtained in IS. Based on 1142 IRGs, 128 IS samples were clustered into two molecular subtypes: clusterA and clusterB. Based on the WGCNA, the authors found that the blue module had the highest correlation with IS. In the blue module, 90 genes were screened as candidate genes. The top 55 genes were selected as the central nodes according to gene degree in protein-protein interactions network of all genes in blue module. Through taking overlap, nine real hub genes were obtained that might distinguish between clusterA subtype and clusterB subtype of IS. The real hub genes (IL7R, ITK, SOD1, CD3D, LEF1, FBL, MAF, DNMT1, and SLAMF1) may be associated with molecular subtypes and immune regulation of IS.

MAEL gene contributes to bovine testicular development through the m5C-mediated splicing

Knowledge of RNA molecules regulating testicular development and spermatogenesis in bulls is essential for elite bull selection and an ideal breeding program. Herein, we performed direct RNA sequencing (DRS) to explore the functional characterization of RNA molecules produced in the testicles of 9 healthy Simmental bulls at three testicular development stages (prepuberty, puberty, and postpuberty). We identified 5,043 differentially expressed genes associated with testicular weight. These genes exhibited more alternative splicing at sexual maturity, particularly alternative 3' (A3) and 5' (A5) splice sites usage and exon skipping (SE). The expression of hub genes in testicular developmental stages was also affected by both m6A and m5C RNA modifications. We found m5C-mediated splicing events significantly (p < 0.05) increased MAEL gene expression at the isoform level, likely promoting spermatogenesis. Our findings highlight the complexity of RNA processing and expression as well as the regulation of transcripts involved in testicular development and spermatogenesis.

Identification of Diagnostic Markers in Infantile Hemangiomas

Background

Infantile Hemangiomas (IHs) are common benign vascular tumors of infancy that may have serious consequences. The research on diagnostic markers for IHs is scarce.

Methods

The "limma" R package was applied to identify differentially expressed genes (DEGs) in developing IHs. Plugin ClueGO in Cytoscape software performed functional enrichment of DEGs. The Search Tool for Retrieving Interacting Genes (STRING) database was utilized to construct the PPI network. The least absolute shrinkage and selection operator (LASSO) regression model and support vector machine recursive feature elimination (SVM-RFE) analysis were used to identify diagnostic genes for IHs. The receiver operating characteristic (ROC) curve evaluated diagnostic genes' discriminatory ability. Single-gene based on Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) was conducted by Gene Set Enrichment Analysis (GSEA). The chemicals related to the diagnostic genes were excavated by the Comparative Toxicogenomics Database (CTD). Finally, the online website Network Analyst was used to predict the transcription factors targeting the diagnostic genes.

Results

A total of 205 DEGs were singled out from IHs samples of 6-, 12-, and 24-month-old infants. These genes principally participated in vasculogenesis and development-related, endothelial cell-related biological processes. Then we mined 127 interacting proteins and created a network with 127 nodes and 251 edges. Furthermore, LASSO and SVM-RRF algorithms identified five diagnostic genes, namely, TMEM2, GUCY1A2, ISL1, WARS, and STEAP4. ROC curve analysis results indicated that the diagnostic genes had a powerful ability to distinguish IHs samples from normal samples. Next, the results of GSEA for a single gene illustrated that all five diagnostic genes inhibited the "valine, leucine, and isoleucine degradation" pathway in the development of IHs. WARS, TMEM2, and STEAP4 activated the "blood vessel development" and "vasculature development" in IHs. Subsequently, inhibitors targeting TMEM2, GUCY1A2, ISL1, and STEAP4 were mined. Finally, 14 transcription factors regulating GUCY1A2, 14 transcription factors regulating STEAP4, and 26 transcription factors regulating ISL1 were predicted.

Conclusion

This study identified five diagnostic markers for IHs and further explored the mechanisms and targeting drugs, providing a basis for diagnosing and treating IHs.

Comparative transcriptome profiling and weighted gene co-expression network analysis to identify core genes in maize (Zea mays L.) silks infected by multiple fungi

Maize (Zea mays L.) is the third most popular Poaceae crop after wheat and rice and used in feed and pharmaceutical sectors. The maize silk contains bioactive components explored by traditional Chinese herbal medicine for various pharmacological activities. However, Fusarium graminearum, Fusarium verticillioides, Trichoderma atroviride, and Ustilago maydis can infect the maize, produce mycotoxins, hamper the quantity and quality of silk production, and further harm the primary consumer's health. However, the defense mechanism is not fully understood in multiple fungal infections in the silk of Z. mays. In this study, we applied bioinformatics approaches to use the publicly available transcriptome data of Z. mays silk affected by multiple fungal flora to identify core genes involved in combatting disease response. Differentially expressed genes (DEGs) were identified among intra- and inter-transcriptome data sets of control versus infected Z. mays silks. Upon further comparison between up- and downregulated genes within the control of datasets, 4,519 upregulated and 5,125 downregulated genes were found. The DEGs have been compared with genes in the modules of weighted gene co-expression network analysis to relevant specific traits towards identifying core genes. The expression pattern of transcription factors, carbohydrate-active enzymes (CAZyme), and resistance genes was analyzed. The present investigation is supportive of our findings that the gene ontology, immunity stimulus, and resistance genes are upregulated, but physical and metabolic processes such as cell wall organizations and pectin synthesis were downregulated respectively. Our results are indicative that terpene synthase TPS6 and TPS11 are involved in the defense mechanism against fungal infections in maize silk.

Single-Cell Analysis Reveals the Role of the Neuropeptide Receptor FPR2 in Monocytes in Kawasaki Disease: A Bioinformatic Study

Exploring the role of neuropeptides in the communication between monocyte subtypes facilitates an investigation of the pathogenesis of Kawasaki disease (KD). We investigated the patterns of interaction between neuropeptide-associated ligands and receptors in monocyte subpopulations in KD patients. Single-cell analysis was employed for the identification of cell subpopulations in KD patients, and monocytes were classified into 3 subpopulations: classical monocytes (CMs), intermediate monocytes (IMs), and nonclassical monocytes (NCMs). Cell-cell communication and differential analyses were used to identify ligand-receptor interactions in monocytes. Five neuropeptide-related genes (SORL1, TNF, SORT1, FPR2, and ANXA1) were involved in cell-cell interactions, wherein FPR2, a neuropeptide receptor, was significantly highly expressed in KD. Weighted gene coexpression network analysis revealed a significant correlation between the yellow module and FPR2 (p < 0.001, CC = 0.43). Using the genes in the yellow module, we constructed a PPI network to assess the possible functions of the FPR2-associated gene network. Gene set enrichment analysis showed that increased FPR2 levels may be involved in immune system regulation. FPR2 in CMs mediates the control of inflammation in KD. The findings of this study may provide a novel target for the clinical treatment of KD.

Eucommia ulmoides Oliver-Tribulus terrestris L. Drug Pair Regulates Ferroptosis by Mediating the Neurovascular-Related Ligand-Receptor Interaction Pathway- A Potential Drug Pair for Treatment Hypertension and Prevention Ischemic Stroke

Background

In this study, we used the network pharmacology approach to explore the potential disease targets of the Eucommia ulmoides Oliver (EUO)-Tribulus terrestris L. (TT) drug pair in the treatment of hypertension-associated neurovascular lesions and IS via the ferroptosis pathway.

Methods

We used the Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform to search for the key active compounds and targets of the drug pair. Based on the GeneCards database, the relevant targets for the drug pair were obtained. Then, we performed the molecular docking of the screened core active ingredients and proteins using the DAVID database and the R AutoDock Vina software. Based on the GSE22255 dataset, these screened target proteins were used to build random forest (RF) and support vector machine (SVM) models. Finally, a new IS nomogram prediction model was constructed and evaluated.

Results

There were 36 active compounds in the EUO-TT drug pair. CHRM1, NR3C1, ADRB2, and OPRD1 proteins of the neuroactive ligand-receptor interaction pathway interacted with the proteins related to the ferroptosis pathway. Molecular docking experiments identified 12 active ingredients of the drug pair that may tightly bind to those target proteins. We constructed a visual IS nomogram prediction model using four genes (CHRM1, NR3C1, ADRB2, and OPRD1). The calibration curve, DCA, and clinical impact curves all indicated that the nomogram model is clinically applicable and diagnostically capable. CHRM1, NR3C1, ADRB2, and OPRD1, the target genes of the four effective components of the EUO-TT drug pair, were considered as risk markers for IS.

Conclusions

The active ingredients of EUO-TT drug pair may act on proteins associated with the neuroactive ligand-receptor interaction pathway to regulate ferroptosis in vascular neurons cells, ultimately affecting the onset and progression of hypertension.

Comprehensive Genetic Analysis of Tuberculosis and Identification of Candidate Biomarkers

Purpose: The purpose of this study is to use the data in the GEO database to analyze, screen biomarkers that can diagnose tuberculosis, and verification of candidate biomarkers.

Materials and methods: GSE158767 dataset were used to process WGCNA analysis, differential gene analysis, Gene ontology and KEGG analysis, protein-protein network analysis and hub genes analysis. Based on our previous study, the intersect between WGCNA and differential gene analysis could be used as candidate biomarkers. Then, the enzyme-linked immunosorbent assay was used to validate candidate biomarkers, and receiver operating characteristic was used to assess diagnose ability of candidate biomarkers.

Results: A total of 412 differential genes were screened. And we obtained 105 overlapping genes between DEGs and WGCNA. GO and KEGG analysis showed that most of the differential genes were significantly enriched in innate immunity. A total of 15 hub genes were screened, and four of them were verified by Enzyme-linked immunosorbent assay. CCL5 performed well in distinguishing the healthy group from the TB group (AUC = 0.723). And CCL19 performed well in distinguishing the TB group from the ORD groups (AUC = 0.811).

Conclusion: CCL19, C1Qb, CCL5 and HLA-DMB may play important role in tuberculosis, which indicated four genes may become effective biomarkers and could be conveniently used to facilitate the individual tuberculosis diagnosis in Chinese people.

Identification of Key Signaling Pathways and Genes in Eosinophilic Asthma and Neutrophilic Asthma by Weighted Gene Co-Expression Network Analysis

Background: Asthma is a heterogeneous disease with different subtypes including eosinophilic asthma (EA) and neutrophilic asthma (NA). However, the mechanisms underlying the difference between the two subtypes are not fully understood.

Methods: Microarray datasets (GSE45111 and GSE137268) were acquired from Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) in induced sputum between EA (n = 24) and NA (n = 15) were identified by “Limma” package. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses and Gene set enrichment analysis (GSEA) were used to explore potential signaling pathways. Weighted gene co-expression network analysis (WGCNA) were performed to identify the key genes that were strongly associated with EA and NA.

Results: A total of 282 DEGs were identified in induced sputum of NA patients compared with EA patients. In GO and KEGG pathway analyses, DEGs were enriched in positive regulation of cytokine production, and cytokine-cytokine receptor interaction. The results of GSEA showed that ribosome, Parkinson’s disease, and oxidative phosphorylation were positively correlated with EA while toll-like receptor signaling pathway, primary immunodeficiency, and NOD-like receptor signaling pathway were positively correlated with NA. Using WGCNA analysis, we identified a set of genes significantly associated NA including IRFG, IRF1, STAT1, IFIH1, IFIT3, GBP1, GBP5, IFIT2, CXCL9, and CXCL11.

Conclusion: We identified potential signaling pathways and key genes involved in the pathogenesis of the asthma subsets, especially in neutrophilic asthma.

Identification of the susceptibility genes for COVID-19 in lung adenocarcinoma with global data and biological computation methods

Introduction

The risk of infection with COVID-19 is high in lung adenocarcinoma (LUAD) patients, and there is a dearth of studies on the molecular mechanism underlying the high susceptibility of LUAD patients to COVID-19 from the perspective of the global differential expression landscape.

Objectives

To fill the research void on the molecular mechanism underlying the high susceptibility of LUAD patients to COVID-19 from the perspective of the global differential expression landscape.

Methods

Herein, we identified genes, specifically the differentially expressed genes (DEGs), correlated with the susceptibility of LUAD patients to COVID-19. These were obtained by calculating standard mean deviation (SMD) values for 49 SARS-CoV-2-infected LUAD samples and 24 non-affected LUAD samples, as well as 3931 LUAD samples and 3027 non-cancer lung samples from 40 pooled RNA-seq and microarray datasets. Hub susceptibility genes significantly related to COVID-19 were further selected by weighted gene co-expression network analysis. Then, the hub genes were further analyzed via an examination of their clinical significance in multiple datasets, a correlation analysis of the immune cell infiltration level, and their interactions with the interactome sets of the A549 cell line.

Results

A total of 257 susceptibility genes were identified, and these genes were associated with RNA splicing, mitochondrial functions, and proteasomes. Ten genes, MEA1, MRPL24, PPIH, EBNA1BP2, MRTO4, RABEPK, TRMT112, PFDN2, PFDN6, and NDUFS3, were confirmed to be the hub susceptibility genes for COVID-19 in LUAD patients, and the hub susceptibility genes were significantly correlated with the infiltration of multiple immune cells.

Conclusion

In conclusion, the susceptibility genes for COVID-19 in LUAD patients discovered in this study may increase our understanding of the high risk of COVID-19 in LUAD patients.

An integrated network analysis approach to identify potential key genes, transcription factors, and microRNAs regulating human hematopoietic stem cell aging

Hematopoietic stem cells (HSCs) undergo functional deterioration with increasing age that causes loss of their self-renewal and regenerative potential. Despite various efforts, significant success in identifying molecular regulators of HSC aging has not been achieved, one prime reason being the non-availability of appropriate human HSC samples. To demonstrate the scope of integrating and re-analyzing the HSC transcriptomics data available, we used existing tools and databases to structure a sequential data analysis pipeline to predict potential candidate genes, transcription factors, and microRNAs simultaneously. This sequential approach comprises (i) collecting matched young and aged mice HSC sample datasets, (ii) identifying differentially expressed genes, (iii) identifying human homologs of differentially expressed genes, (iv) inferring gene co-expression network modules, and (v) inferring the microRNA-transcription factor-gene regulatory network. Systems-level analyses of HSC interaction networks provided various insights based on which several candidates were predicted. For example, 16 HSC aging-related candidate genes were predicted (e.g., CD38, BRCA1, AGTR1, GSTM1, etc.) from GCN analysis. Following this, the shortest path distance-based analyses of the regulatory network predicted several novel candidate miRNAs and TFs. Among these, miR-124-3p was a common regulator in candidate gene modules, while TFs MYC and SP1 were identified to regulate various candidate genes. Based on the regulatory interactions among candidate genes, TFs, and miRNAs, a potential regulation model of biological processes in each of the candidate modules was predicted, which provided systems-level insights into the molecular complexity of each module to regulate HSC aging.

Using weighted gene co-expression network analysis (WGCNA) to identify the hub genes related to hypoxic adaptation in yak (Bos grunniens)

BACKGROUND

As a mammal living at the highest altitude in the world, the yak has strong adaptability to the harsh natural environment (such as low temperature, scarce food, especially low oxygen) of Qinghai-Tibet Plateau (QTP) after a long process of natural selection.

OBJECTIVE

Here, we used Weighted Correlation Network Analysis (WGCNA), a systematic biology method, to identify hypoxic adaptation-related modules and hub genes. The research of the adaptability of yak against hypoxia is of great significance to identify the genetic characteristics and yak breeding.

METHODS

Based on the transcriptome sequencing data (PRJNA362606), the R package DESeq2 and WGCNA were conducted to analyze differentially expressed genes (DEGs) and construct the gene co-expression network. The module hub genes were identified and characterized by the correlation of gene and trait, module membership (kME). In addition, GO and KEGG enrichment analyses were used to explore the functions of hub genes.

RESULTS

Our results revealed that 1098, 1429, and 1645 DEGs were identified in muscle, spleen, and lung, respectively. Besides, a total of 13 gene co-expression modules were detected, of which two hypoxic adaptation-related modules (saddlebrown and turquoise) were found. We identified 39 and 150 hub genes in these two modules. Functional enrichment analyses showed that 12 GO terms and 18 KEGG pathways were enriched in the saddlebrown module while 85 GO terms and 22 KEGG pathways were enriched in the turquoise module. The significant pathways related to hypoxia adaptation include FoxO signaling pathway, Thermogenesis pathway, and Retrograde endocannabinoid signaling pathway, etc. CONCLUSIONS: In this study, we obtained two hypoxia-related specific modules and identified hub genes based on the connectivity by constructing a weighted gene co-expression network. Function enrichment analysis of two modules revealed mitochondrion is the most important organelle for hypoxia adaptation. Moreover, the insulin-related pathways and thermogenic-related pathways played a major role. The results of this study provide theoretical guidance for further understanding the molecular mechanism of yak adaptation to hypoxia.