Diabetic cardiomyopathy (DCM) is a common cardiovascular complication of diabetes, which may threaten the quality of life and shorten life expectancy in the diabetic population. However, the molecular mechanisms underlying the diabetes cardiomyopathy are not fully elucidated. We analyzed two datasets from Gene Expression Omnibus (GEO). Differentially expressed and weighted gene correlation network analysis (WGCNA) was used to screen key genes and molecules. Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis, and protein-protein interaction (PPI) network analysis were constructed to identify hub genes. The diagnostic value of the hub gene was evaluated using the receiver operating characteristic (ROC). Quantitative real-time PCR (RT-qPCR) was used to validate the hub genes. A total of 13 differentially co-expressed modules were selected by WGCNA and differential expression analysis. KEGG and GO analysis showed these DEGs were mainly enriched in lipid metabolism and myocardial hypertrophy pathway, cytomembrane, and mitochondrion. As a result, six genes were identified as hub genes. Finally, five genes (Pdk4, Lipe, Serpine1, Igf1r, and Bcl2l1) were found significantly changed in both the validation dataset and experimental mice with DCM. In conclusion, the present study identified five genes that may help provide novel targets for diagnosing and treating DCM.

The coronavirus disease 2019 (COVID-19) pandemic has had a significant impact globally, resulting in a higher death toll and persistent health issues for survivors, particularly those with pre-existing medical conditions. Numerous studies have demonstrated a strong correlation between catastrophic COVID-19 results and diabetes. To gain deeper insights, we analysed the transcriptome dataset from COVID-19 and diabetic peripheral neuropathic patients. Using the R programming language, differentially expressed genes (DEGs) were identified and classified based on up and down regulations. The overlaps of DEGs were then explored between these groups. Functional annotation of those common DEGs was performed using Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), Bio-Planet, Reactome, and Wiki pathways. A protein-protein interaction (PPI) network was created with bioinformatics tools to understand molecular interactions. Through topological analysis of the PPI network, we determined hub gene modules and explored gene regulatory networks (GRN). Furthermore, the study extended to suggesting potential drug molecules for the identified mutual DEG based on the comprehensive analysis. These approaches may contribute to understanding the molecular intricacies of COVID-19 in diabetic peripheral neuropathy patients through insights into potential therapeutic interventions.

UNASSIGNED: Kidney transplantation (KT) is the best treatment for end-stage renal disease. Although long and short-term survival rates for the graft have improved significantly with the development of immunosuppressants, acute rejection (AR) remains a major risk factor attacking the graft and patients. The innate immune response plays an important role in rejection. Therefore, our objective is to determine the biomarkers of congenital immunity associated with AR after KT and provide support for future research.
UNASSIGNED: A differential expression genes (DEGs) analysis was performed based on the dataset GSE174020 from the NCBI gene Expression Synthesis Database (GEO) and then combined with the GSE5099 M1 macrophage-related gene identified in the Molecular Signatures Database. We then identified genes in DEGs associated with M1 macrophages defined as DEM1Gs and performed gene ontology (GO) and Kyoto Encyclopedia of Genomes (KEGG) enrichment analysis. Cibersort was used to analyze the immune cell infiltration during AR. At the same time, we used the protein-protein interaction (PPI) network and Cytoscape software to determine the key genes. Dataset, GSE14328 derived from pediatric patients, GSE138043 and GSE9493 derived from adult patients, were used to verify Hub genes. Additional verification was the rat KT model, which was used to perform HE staining, immunohistochemical staining, and Western Blot. Hub genes were searched in the HPA database to confirm their expression. Finally, we construct the interaction network of transcription factor (TF)-Hub genes and miRNA-Hub genes.
UNASSIGNED: Compared to the normal group, 366 genes were upregulated, and 423 genes were downregulated in the AR group. Then, 106 genes related to M1 macrophages were found among these genes. GO and KEGG enrichment analysis showed that these genes are mainly involved in cytokine binding, antigen binding, NK cell-mediated cytotoxicity, activation of immune receptors and immune response, and activation of the inflammatory NF-κB signaling pathway. Two Hub genes, namely CCR7 and CD48, were identified by PPI and Cytoscape analysis. They have been verified in external validation sets, originated from both pediatric patients and adult patients, and animal experiments. In the HPA database, CCR7 and CD48 are mainly expressed in T cells, B cells, macrophages, and tissues where these immune cells are distributed. In addition to immunoinfiltration, CD4+T, CD8+T, NK cells, NKT cells, and monocytes increased significantly in the AR group, which was highly consistent with the results of Hub gene screening. Finally, we predicted that 19 TFs and 32 miRNAs might interact with the Hub gene.
UNASSIGNED: Through a comprehensive bioinformatic analysis, our findings may provide predictive and therapeutic targets for AR after KT.

Bone non-union is a common fracture complication that can severely impact patient outcomes, yet its mechanism is not fully understood. This study used differential analysis and weighted co-expression network analysis (WGCNA) to identify susceptibility modules and hub genes associated with fracture healing. Two datasets, GSE125289 and GSE213891, were downloaded from the GEO website, and differentially expressed miRNAs and genes were analysed and used to construct the WGCNA network. Gene ontology (GO) analysis of the differentially expressed genes showed enrichment in cytokine and inflammatory factor secretion, phagocytosis, and trans-Golgi network regulation pathways. Using bioinformatic site prediction and crossover gene search, miR-29b-3p was identified as a regulator of LIN7A expression that may negatively affect fracture healing. Potential miRNA-mRNA interactions in the bone non-union mechanism were explored, and miRNA-29-3p and LIN7A were identified as biomarkers of skeletal non-union. The expression of miRNA-29b-3p and LIN7A was verified in blood samples from patients with fracture non-union using qRT-PCR and ELISA. Overall, this study identified characteristic modules and key genes associated with fracture non-union and provided insight into its molecular mechanisms. Downregulated miRNA-29b-3p was found to downregulate LIN7A protein expression, which may affect the healing process after fracture in patients with bone non-union. These findings may serve as a prognostic biomarker and potential therapeutic target for bone non-union.

BACKGROUND: A widely used approach for extracting information from gene expression data employs the construction of a gene co-expression network and the subsequent computational detection of gene clusters, called modules. WGCNA and related methods are the de facto standard for module detection. The purpose of this work is to investigate the applicability of more sophisticated algorithms toward the design of an alternative method with enhanced potential for extracting biologically meaningful modules.
RESULTS: We present self-learning gene clustering pipeline (SGCP), a spectral method for detecting modules in gene co-expression networks. SGCP incorporates multiple features that differentiate it from previous work, including a novel step that leverages gene ontology (GO) information in a self-leaning step. Compared with widely used existing frameworks on 12 real gene expression datasets, we show that SGCP yields modules with higher GO enrichment. Moreover, SGCP assigns highest statistical importance to GO terms that are mostly different from those reported by the baselines.
CONCLUSIONS: Existing frameworks for discovering clusters of genes in gene co-expression networks are based on relatively simple algorithmic components. SGCP relies on newer algorithmic techniques that enable the computation of highly enriched modules with distinctive characteristics, thus contributing a novel alternative tool for gene co-expression analysis.

BACKGROUND: Dual RNA sequencing is a powerful tool that enables a comprehensive understanding of the molecular dynamics underlying plant-microbe interactions. RNA sequencing (RNA-seq) poses technical hurdles in the transcriptional analysis of plant-bacterial interactions, especially in bacterial transcriptomics, owing to the presence of abundant ribosomal RNA (rRNA), which potentially limits the coverage of essential transcripts. Therefore, to achieve cost-effective and comprehensive sequencing of the bacterial transcriptome, it is imperative to devise efficient methods for eliminating rRNA and enhancing the proportion of bacterial mRNA. In this study, we modified a strand-specific dual RNA-seq method with the goal of enriching the proportion of bacterial mRNA in the bacteria-infected plant samples. The enriched method involved the sequential separation of plant mRNA by poly A selection and rRNA removal for bacterial mRNA enrichment followed by strand specific RNA-seq library preparation steps. We assessed the efficiency of the enriched method in comparison to the conventional method by employing various plant-bacterial interactions, including both host and non-host resistance interactions with pathogenic bacteria, as well as an interaction with a beneficial rhizosphere associated bacteria using pepper and tomato plants respectively.
RESULTS: In all cases of plant-bacterial interactions examined, an increase in mapping efficiency was observed with the enriched method although it produced a lower read count. Especially in the compatible interaction with Xanthmonas campestris pv. Vesicatoria race 3 (Xcv3), the enriched method enhanced the mapping ratio of Xcv3-infected pepper samples to its own genome (15.09%; 1.45-fold increase) and the CDS (8.92%; 1.49-fold increase). The enriched method consistently displayed a greater number of differentially expressed genes (DEGs) than the conventional RNA-seq method at all fold change threshold levels investigated, notably during the early stages of Xcv3 infection in peppers. The Gene Ontology (GO) enrichment analysis revealed that the DEGs were predominantly enriched in proteolysis, kinase, serine type endopeptidase and heme binding activities.
CONCLUSIONS: The enriched method demonstrated in this study will serve as a suitable alternative to the existing RNA-seq method to enrich bacterial mRNA and provide novel insights into the intricate transcriptomic alterations within the plant-bacterial interplay.

Automated protein function prediction is a crucial and widely studied problem in bioinformatics. Computationally, protein function is a multilabel classification problem where only positive samples are defined and there is a large number of unlabeled annotations. Most existing methods rely on the assumption that the unlabeled set of protein function annotations are negatives, inducing the false negative issue, where potential positive samples are trained as negatives. We introduce a novel approach named PU-GO, wherein we address function prediction as a positive-unlabeled ranking problem. We apply empirical risk minimization, i.e. we minimize the classification risk of a classifier where class priors are obtained from the Gene Ontology hierarchical structure. We show that our approach is more robust than other state-of-the-art methods on similarity-based and time-based benchmark datasets.
METHODS: Data and code are available at https://github.com/bio-ontology-research-group/PU-GO.

BACKGROUND: Wilms tumor (WT) is the most common pediatric embryonal tumor. Improving patient outcomes requires advances in understanding and targeting the multiple genes and cellular control pathways, but its pathogenesis is currently not well-researched. We aimed to identify the potential molecular biological mechanism of WT and develop new prognostic markers and molecular targets by comparing gene expression profiles of Wilms tumors and fetal normal kidneys.
METHODS: Differential gene expression analysis was performed on Wilms tumor transcriptomic data from the GEO and TARGET databases. For biological functional analysis, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment were utilized. Out of 24 hub genes identified, nine were found to be prognostic-related through univariate Cox regression analysis. These nine genes underwent LASSO regression analysis to enhance the predictive capability of the model. The key hub genes were validated in the GSE73209 datasets, and cell function experiments were conducted to identify the genes\' functions in WiT-49 cells.
RESULTS: The enrichment analysis revealed that DEGs were significantly involved in the regulation of angiogenesis and regulation of cell differentiation. 24 DEGs were identified through PPI networks and the MCODE algorithm, and 9 of 24 genes were related to WT patients\' prognosis. EMCN and CCNA1 were identified as key hub genes, and related to the progression of WT. Functionally, over-expression of EMCN and CCNA1 knockdown inhibited cell viability, proliferation, migration, and invasion of Wilms tumor cells.
CONCLUSIONS: EMCN and CCNA1 were identified as key prognostic markers in Wilms tumor, suggesting their potential as therapeutic targets. Differential gene expression and enrichment analyses indicate significant roles in angiogenesis and cell differentiation.

BACKGROUND: Long noncoding RNAs (lncRNAs) are noncoding RNA transcripts greater than 200 nucleotides in length and are known to play a role in regulating the transcription of genes involved in vital cellular functions. We hypothesized the disease process in dysferlinopathy is linked to an aberrant expression of lncRNAs and messenger RNAs (mRNAs).
OBJECTIVE: In this study, we compared the lncRNA and mRNA expression profiles between wild-type and dysferlin-deficient murine myoblasts (C2C12 cells).
METHODS: LncRNA and mRNA expression profiling were performed using a microarray. Several lncRNAs with differential expression were validated using quantitative real-time polymerase chain reaction. Gene Ontology (GO) analysis was performed to understand the functional role of the differentially expressed mRNAs. Further bioinformatics analysis was used to explore the potential function, lncRNA-mRNA correlation, and potential targets of the differentially expressed lncRNAs.
RESULTS: We found 3195 lncRNAs and 1966 mRNAs that were differentially expressed. The chromosomal distribution of the differentially expressed lncRNAs and mRNAs was unequal, with chromosome 2 having the highest number of lncRNAs and chromosome 7 having the highest number of mRNAs that were differentially expressed. Pathway analysis of the differentially expressed genes indicated the involvement of several signaling pathways including PI3K-Akt, Hippo, and pathways regulating the pluripotency of stem cells. The differentially expressed genes were also enriched for the GO terms, developmental process and muscle system process. Network analysis identified 8 statistically significant (P<.05) network objects from the upregulated lncRNAs and 3 statistically significant network objects from the downregulated lncRNAs.
CONCLUSIONS: Our results thus far imply that dysferlinopathy is associated with an aberrant expression of multiple lncRNAs, many of which may have a specific function in the disease process. GO terms and network analysis suggest a muscle-specific role for these lncRNAs. To elucidate the specific roles of these abnormally expressed noncoding RNAs, further studies engineering their expression are required.

BACKGROUND: Tuberculosis is a serious life-threatening disease among the top global health challenges and rapid and effective diagnostic biomarkers are vital for early diagnosis especially given the increasing prevalence of multidrug resistance.
METHODS: Two human whole blood microarray datasets, GSE42826 and GSE42830 were retrieved from publicly available gene expression omnibus (GEO) database. Deregulated genes (DEGs) were identified using GEO2R online tool and Gene Ontology (GO), protein-protein interaction (PPI) network analysis was performed using Metascape and STRING databases. Significant genes (n = 8) were identified using T-test/ANOVA and Molecular Complex Detection (MCODE) score ≥10, which was validated in GSE34608 dataset. The diagnostic potential of three biomarkers was assessed using Area Under Curve (AUC) of Receiver Operating Characteristic (ROC) plot. The transcriptional levels of these genes were also examined in a separate dataset GSE31348, to monitor the patterns of variation during tuberculosis treatment.
RESULTS: A total of 62 common DEGs (57 upregulated, 7 downregulated genes) were identified in two discovery datasets. GO functions and pathway enrichment analysis shed light on the functional roles of these DEGs in immune response and type-II interferon signaling. The genes in Module-1 (n = 18) were linked to innate immune response, interferon-gamma signaling. The common genes (n = 8) were validated in GSE34608 dataset, that corroborates the results obtained from discovery sets. The gene expression levels demonstrated responsiveness to Mtb infection during anti-TB therapy in GSE31348 dataset. In GSE34608 dataset, the expression levels of three specific genes, GBP5, IFITM3, and EPSTI1, emerged as potential diagnostic makers. In combination, these genes scored remarkable diagnostic performance with 100% sensitivity and 89% specificity, resulting in an impressive Area Under Curve (AUC) of 0.958. However, GBP5 alone showed the highest AUC of 0.986 with 100% sensitivity and 89% specificity.
CONCLUSIONS: The study presents valuable insights into the critical gene network perturbed during tuberculosis. These genes are determinants for assessing the effectiveness of an anti-TB response and distinguishing between active TB and healthy individuals. GBP5, IFITM3 and EPSTI1 emerged as candidate core genes in TB and holds potential as novel molecular targets for the development of interventions in the treatment of TB.