Linear Models For Microarray Data Research Articles

A diagnostic model based on bioinformatics and machine learning to differentiate bipolar disorder from schizophrenia and major depressive disorder

Bipolar disorder (BD) showed the highest suicide rate of all psychiatric disorders, and its underlying causative genes and effective treatments remain unclear. During diagnosis, BD is often confused with schizophrenia (SC) and major depressive disorder (MDD), due to which patients may receive inadequate or inappropriate treatment, which is detrimental to their prognosis. This study aims to establish a diagnostic model to distinguish BD from SC and MDD in multiple public datasets through bioinformatics and machine learning and to provide new ideas for diagnosing BD in the future. Three brain tissue datasets containing BD, SC, and MDD were chosen from the Gene Expression Omnibus database (GEO), and two peripheral blood datasets were selected for validation. Linear Models for Microarray Data (Limma) analysis was carried out to identify differentially expressed genes (DEGs). Functional enrichment analysis and machine learning were utilized to identify. Least absolute shrinkage and selection operator (LASSO) regression was employed for identifying candidate immune-associated central genes, constructing protein-protein interaction networks (PPI), building artificial neural networks (ANN) for validation, and plotting receiver operating characteristic curve (ROC curve) for differentiating BD from SC and MDD and creating immune cell infiltration to study immune cell dysregulation in the three diseases. RBM10 was obtained as a candidate gene to distinguish BD from SC. Five candidate genes (LYPD1, HMBS, HEBP2, SETD3, and ECM2) were obtained to distinguish BD from MDD. The validation was performed by ANN, and ROC curves were plotted for diagnostic value assessment. The outcomes exhibited the prediction model to have a promising diagnostic value. In the immune infiltration analysis, Naive B, Resting NK, and Activated Mast Cells were found to be substantially different between BD and SC. Naive B and Memory B cells were prominently variant between BD and MDD. In this study, RBM10 was found as a candidate gene to distinguish BD from SC; LYPD1, HMBS, HEBP2, SETD3, and ECM2 serve as five candidate genes to distinguish BD from MDD. The results obtained from the ANN network showed that these candidate genes could perfectly distinguish BD from SC and MDD (76.923% and 81.538%, respectively).

Identification of age-specific biomarkers of spinal cord injury: A bioinformatics analysis of young and aged mice models.

Spinal cord injury (SCI) is a debilitating event that often results in long-term physical damage, disability, and a significant impact on a patient's quality of life. Past research has noted an age-dependent decline in regeneration post-SCI. This study seeks to identify potential biomarkers and enriched pathways in young and aged SCI mouse models. We retrieved the microarray data of spinal cord samples from SCI mice and control mice from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) were identified using the R software and the Linear Models for Microarray Data (limma) package. The Gene Set Enrichment Analysis (GSEA) determined enrichment differences among gene sets. The Weighted Gene Co-expression Network Analysis (WGCNA) pinpointed clinically significant modules and hub genes in SCI. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis was employed to uncover significantly related pathways of crucial genes in SCI. We found 2560 DEGs in the young SCI group, comprised of 1733 upregulated RNAs and 827 downregulated RNAs. In the aged SCI group, 3048 DEGs were revealed including 1856 upregulated and 1192 downregulated genes. The GSEA revealed 12 enriched signaling pathways in the young SCI group, such as IL6/JAK/STAT3 signaling, interferon alpha response, and interferon gamma response, and 21 signaling pathways in the aged SCI group such as IL6/JAK/STAT3 signaling, E2F targets, and angiogenesis-related pathways. The WGCNA identified the turquoise module as significantly associated with the clinical traits of both young and aged SCI, and revealed 3181 hub genes. Ultimately, 1858 significant genes in SCI were found, with associated signaling pathways including epithelial-mesenchymal transition (EMT), interferon gamma response, and KARS signaling. Our study explored the RNA expression patterns and enriched signaling pathways in young and aged SCI mice. These findings may provide potential biomarkers for targeted SCI therapy.

Identification and assessment of differentially expressed necroptosis long non-coding RNAs associated with periodontitis in human

BackgroundPeriodontitis is the most common oral disease and is closely related to immune infiltration in the periodontal microenvironment and its poor prognosis is related to the complex immune response. The progression of periodontitis is closely related to necroptosis, but there is still no systematic study of long non-coding RNA (lncRNA) associated with necroptosis for diagnosis and treatment of periodontitis.Material and methodsTranscriptome data and clinical data of periodontitis and healthy populations were obtained from the Gene Expression Omnibus (GEO) database, and necroptosis-related genes were obtained from previously published literature. FactoMineR package in R was used to perform principal component analysis (PCA) for obtaining the necroptosis-related lncRNAs. The core necroptosis-related lncRNAs were screened by the Linear Models for Microarray Data (limma) package in R, PCA principal component analysis and lasso algorithm. These lncRNAs were then used to construct a classifier for periodontitis with logistic regression. The receiver operating characteristic (ROC) curve was used to evaluate the sensitivity and specificity of the model. The CIBERSORT method and ssGSEA algorithm were used to estimate the immune infiltration and immune pathway activation of periodontitis. Spearman’s correlation analysis was used to further verify the correlation between core genes and periodontitis immune microenvironment. The expression level of core genes in human periodontal ligament cells (hPDLCs) was detected by RT-qPCR.ResultsA total of 10 core necroptosis-related lncRNAs (10-lncRNAs) were identified, including EPB41L4A-AS1, FAM30A, LINC01004, MALAT1, MIAT, OSER1-DT, PCOLCE-AS1, RNF144A-AS1, CARMN, and LINC00582. The classifier for periodontitis was successfully constructed. The Area Under the Curve (AUC) was 0.952, which suggested that the model had good predictive performance. The correlation analysis of 10-lncRNAs and periodontitis immune microenvironment showed that 10-lncRNAs had an impact on the immune infiltration of periodontitis. Notably, the RT-qPCR results showed that the expression level of the 10-lncRNAs obtained was consistent with the chip analysis results.ConclusionsThe 10-lncRNAs identified from the GEO dataset had a significant impact on the immune infiltration of periodontitis and the classifier based on 10-lncRNAs had good detection efficiency for periodontitis, which provided a new target for diagnosis and treatment of periodontitis.

Identification of a miRNA-mRNA regulatory network for post-stroke depression: a machine-learning approach.

The study aimed to explore the miRNA and mRNA biomarkers in post-stroke depression (PSD) and to develop a miRNA-mRNA regulatory network to reveal its potential pathogenesis. The transcriptomic expression profile was obtained from the GEO database using the accession numbers GSE117064 (miRNAs, stroke vs. control) and GSE76826 [mRNAs, late-onset major depressive disorder (MDD) vs. control]. Differentially expressed miRNAs (DE-miRNAs) were identified in blood samples collected from stroke patients vs. control using the Linear Models for Microarray Data (LIMMA) package, while the weighted correlation network analysis (WGCNA) revealed co-expressed gene modules correlated with the subject group. The intersection between DE-miRNAs and miRNAs identified by WGCNA was defined as stroke-related miRNAs, whose target mRNAs were stroke-related genes with the prediction based on three databases (miRDB, miRTarBase, and TargetScan). Using the GSE76826 dataset, the differentially expressed genes (DEGs) were identified. Overlapped DEGs between stroke-related genes and DEGs in late-onset MDD were retrieved, and these were potential mRNA biomarkers in PSD. With the overlapped DEGs, three machine-learning methods were employed to identify gene signatures for PSD, which were established with the intersection of gene sets identified by each algorithm. Based on the gene signatures, the upstream miRNAs were predicted, and a miRNA-mRNA network was constructed. Using the GSE117064 dataset, we retrieved a total of 667 DE-miRNAs, which included 420 upregulated and 247 downregulated ones. Meanwhile, WGCNA identified two modules (blue and brown) that were significantly correlated with the subject group. A total of 117 stroke-related miRNAs were identified with the intersection of DE-miRNAs and WGCNA-related ones. Based on the miRNA-mRNA databases, we identified a list of 2,387 stroke-related genes, among which 99 DEGs in MDD were also embedded. Based on the 99 overlapped DEGs, we identified three gene signatures (SPATA2, ZNF208, and YTHDC1) using three machine-learning classifiers. Predictions of the three mRNAs highlight four miRNAs as follows: miR-6883-5p, miR-6873-3p, miR-4776-3p, and miR-6738-3p. Subsequently, a miRNA-mRNA network was developed. The study highlighted gene signatures for PSD with three genes (SPATA2, ZNF208, and YTHDC1) and four upstream miRNAs (miR-6883-5p, miR-6873-3p, miR-4776-3p, and miR-6738-3p). These biomarkers could further our understanding of the pathogenesis of PSD.

Predictive model, miRNA-TF network, related subgroup identification and drug prediction of ischemic stroke complicated with mental disorders based on genes related to gut microbiome.

Patients with comorbid schizophrenia, depression, drug use, and multiple psychiatric diagnoses have a greater risk of carotid revascularization following stroke. The gut microbiome (GM) plays a crucial role in the attack of mental illness and IS, which may become an index for the diagnosis of IS. A genomic study of the genetic commonalities between SC and IS, as well as its mediated pathways and immune infiltration, will be conducted to determine how schizophrenia contributes to the high prevalence of IS. According to our study, this could be an indicator of ischemic stroke development. We selected two datasets of IS from the Gene Expression Omnibus (GEO), one for training and the other for the verification group. Five genes related to mental disorders and GM were extracted from Gene cards and other databases. Linear models for microarray data (Limma) analysis was utilized to identify differentially expressed genes (DEGs) and perform functional enrichment analysis. It was also used to conduct machine learning exercises such as random forest and regression to identify the best candidate for immune-related central genes. Protein-protein interaction (PPI) network and artificial neural network (ANN) were established for verification. The receiver operating characteristic (ROC) curve was drawn for the diagnosis of IS, and the diagnostic model was verified by qRT-PCR. Further immune cell infiltration analysis was performed to study the IS immune cell imbalance. We also performed consensus clustering (CC) to analyze the expression of candidate models under different subtypes. Finally, miRNA, transcription factors (TFs), and drugs related to candidate genes were collected through the Network analyst online platform. Through comprehensive analysis, a diagnostic prediction model with good effect was obtained. Both the training group (AUC 0.82, CI 0.93-0.71) and the verification group (AUC 0.81, CI 0.90-0.72) had a good phenotype in the qRT-PCR test. And in verification group 2 we validated between the two groups with and without carotid-related ischemic cerebrovascular events (AUC 0.87, CI 1-0.64). Furthermore, we investigated cytokines in both GSEA and immune infiltration and verified cytokine-related responses by flow cytometry, particularly IL-6, which played an important role in IS occurrence and progression. Therefore, we speculate that mental illness may affect the development of IS in B cells and IL-6 in T cells. MiRNA (hsa-mir-129-2-3p, has-mir-335-5p, and has-mir-16-5p) and TFs (CREB1, FOXL1), which may be related to IS, were obtained. Through comprehensive analysis, a diagnostic prediction model with good effect was obtained. Both the training group (AUC 0.82, CI 0.93-0.71) and the verification group (AUC 0.81, CI 0.90-0.72) had a good phenotype in the qRT-PCR test. And in verification group 2 we validated between the two groups with and without carotid-related ischemic cerebrovascular events (AUC 0.87, CI 1-0.64). MiRNA (hsa-mir-129-2-3p, has-mir-335-5p, and has-mir-16-5p) and TFs (CREB1, FOXL1), which may be related to IS, were obtained.

Identification of hub genes associated with osteoporosis development by comprehensive bioinformatics analysis.

Osteoporosis (OP) is a systemic bone disease caused by various factors, including, the decrease of bone density and quality, the destruction of bone microstructure, and the increase of bone fragility. It is a disease with a high incidence in a large proportion of the world's elderly population. However, osteoporosis lacks obvious symptoms and sensitive biomarkers. Therefore, it is extremely urgent to discover and identify disease-related biomarkers for early clinical diagnosis and effective intervention for osteoporosis. In our study, the Linear Models for Microarray Data (LIMMA) tool was used to screen differential expressed genes from transcriptome sequencing data of OP blood samples downloaded from the GEO database, and cluster Profiler was used for enriching analysis of differently expressed genes. In order to analyzed the relevance of gene modules, clinical symptoms, and the most related module setting genes associated with disease progression, we adapted Weighted Gene Co-expression Network Analysis (WGCNA) to screen and analyze the related pathways and relevant molecules. We used the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) database to construct protein interaction network of key modules, and Cytoscape software was used to complete network visualization and screen of core genes in the network. Various plug-in algorithms of cytoHubba were used to identify key genes of OP. Finally, correlation analysis and single-gene gene probe concentration analysis (GSEA) analysis were performed for each core gene. Results of a total of 8 key genes that were closely related to the occurrence and development of OP were screened out, which provided a brand-new idea for the clinical diagnosis and early prevention of OP. Quantitative real-time PCR (qRT-PCR) was performed for validation, the expression levels of CUL1, PTEN and STAT1 genes in the OS group were significantly higher than in the non-OS groups. Receiver operating characteristic analysis demonstrated that CUL1, PTEN and STAT1 displayed considerable diagnostic accuracy for OS.

The molecular features of chronic lung allograft dysfunction in lung transplant airway mucosa.

Many lung transplants fail due to chronic lung allograft dysfunction (CLAD). We recently showed that transbronchial biopsies (TBBs) from CLAD patients manifest severe parenchymal injury and dedifferentiation, distinct from time-dependent changes. The present study explored time-selective and CLAD-selective transcripts in mucosal biopsies from the third bronchial bifurcation (3BMBs), compared to those in TBBs. We used genome-wide microarray measurements in 324 3BMBs to identify CLAD-selective changes as well as time-dependent changes and develop a CLAD classifier. CLAD-selective transcripts were identified with linear models for microarray data (limma) and were used to build an ensemble of 12 classifiers to predict CLAD. Hazard models and random forests were then used to predict the risk of graft loss using the CLAD classifier, transcript sets associated with rejection, injury, and time. T cell-mediated rejection and donor-specific antibody were increased in CLAD 3BMBs but most had no rejection. Like TBBs, 3BMBs showed a time-dependent increase in transcripts expressed in inflammatory cells that was not associated with CLAD or survival. Also like TBBs, the CLAD-selective transcripts in 3BMBs reflected severe parenchymal injury and dedifferentiation, not inflammation or rejection. While 3BMBs and TBBs did not overlap in their top 20 CLAD-selective transcripts, many CLAD-selective transcripts were significantly increased in both for example LOXL1, an enzyme controlling matrix remodeling. In Cox models for one-year survival, the 3BMB CLAD-selective transcripts and CLAD classifier predicted graft loss and correlated with CLAD stage. Many 3BMB CLAD-selective transcripts were also increased by injury in kidney transplants and correlated with decreased kidney survival, including LOXL1. Mucosal and transbronchial biopsies from CLAD patients reveal a diffuse molecular injury and dedifferentiation state that impacts prognosis and correlates with the physiologic disturbances. CLAD state in lung transplants shares features with failing kidney transplants, indicating elements shared by the injury responses of distressed organs.

Gene Correlation Network Analysis to Identify Biomarkers of Peri-Implantitis

Background and Objectives: The histopathological and clinical conditions for transforming peri-implant mucositis into peri-implantitis (PI) are not fully clarified. We aim to uncover molecular mechanisms and new potential biomarkers of PI. Materials and Methods: Raw GSE33774 and GSE57631 datasets were obtained from the Gene Expression Omnibus (GEO) database. The linear models for microarray data (LIMMA) package in R software completes differentially expressed genes (DEGs). We conducted a weighted gene co-expression network analysis (WGCNA) on the top 25% of altered genes and identified the key modules associated with the clinical features of PI. Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed using the R software. We constructed a protein–protein interaction (PPI) network through the STRING database. After that we used Cytohubba plug-ins of Cytoscape to screen out the potential hub genes, which were subsequently verified via receiver operating characteristic (ROC) curves in another dataset, GSE178351, and revalidation of genes through the DisGeNET database. Results: We discovered 632 DEGs (570 upregulated genes and 62 downregulated genes). A total of eight modules were screened by WGCNA, among which the turquoise module was most correlated with PI. The Cytohubba plug-ins were used for filtering hub genes, which are highly linked with PI development, from the candidate genes in the protein–protein interaction (PPI) network. Conclusions: We found five key genes from PI using WGCNA. Among them, ICAM1, CXCL1, and JUN are worthy of further study of new target genes, providing the theoretical basis for further exploration of the occurrence and development mechanism of PI.

Revealing Potential Spinal Cord Injury Biomarkers and Immune Cell Infiltration Characteristics in Mice.

Spinal cord injury (SCI) is a disabling condition with significant morbidity and mortality. Currently, no effective SCI treatment exists. This study aimed to identify potential biomarkers and characterize the properties of immune cell infiltration during this pathological event. To eliminate batch effects, we concurrently analyzed two mouse SCI datasets (GSE5296, GSE47681) from the GEO database. First, we identified differentially expressed genes (DEGs) using linear models for microarray data (LIMMA) and performed functional enrichment studies on those DEGs. Next, we employed bioinformatics and machine-learning methods to identify and define the characteristic genes of SCI. Finally, we validated them using immunofluorescence and qRT-PCR. Additionally, this study assessed the inflammatory status of SCI by identifying cell types using CIBERSORT. Furthermore, we investigated the link between key markers and infiltrating immune cells. In total, we identified 561 robust DEGs. We identified Rab20 and Klf6 as SCI-specific biomarkers and demonstrated their significance using qRT-PCR in the mouse model. According to the examination of immune cell infiltration, M0, M1, and M2 macrophages, along with naive CD8, dendritic cell-activated, and CD4 Follicular T cells may have a role in the progression of SCI. Therefore, Rab20 and Klf6 could be accessible targets for diagnosing and treating SCI. Moreover, as previously stated, immune cell infiltration may significantly impact the development and progression of SCI.

Bioinformatics prediction of potential mechanisms and biomarkers underlying dilated cardiomyopathy.

BACKGROUNDHeart failure is a health burden responsible for high morbidity and mortality worldwide, and dilated cardiomyopathy (DCM) is one of the most common causes of heart failure. DCM is a disease of the heart muscle and is characterized by enlargement and dilation of at least one ventricle alongside impaired contractility with left ventricular ejection fraction < 40%. It is also associated with abnormalities in cytoskeletal proteins, mitochondrial ATP transporter, microvasculature, and fibrosis. However, the pathogenesis and potential biomarkers of DCM remain to be investigated.AIMTo investigate the candidate genes and pathways involved in DCM patients.METHODSTwo expression datasets (GSE3585 and GSE5406) were downloaded from the Gene Expression Omnibus database. The differentially expressed genes (DEGs) between the DCM patients and healthy individuals were identified using the R package “linear models for microarray data.” The pathways with common DEGs were analyzed via Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and gene set enrichment analyses. Moreover, a protein-protein interaction network (PPI) was constructed to identify the hub genes and modules. The MicroRNA Database was applied to predict the microRNAs (miRNAs) targeting the hub genes. Additionally, immune cell infiltration in DCM was analyzed using CIBERSORT.RESULTSIn total, 97 DEGs (47 upregulated and 50 downregulated) were identified. GO analysis showed that the DEGs were mainly enriched in “response to growth factor,” “extracellular matrix,” and “extracellular matrix structural constituent.” KEGG pathway analysis indicated that the DEGs were mainly enriched in “protein digestion and absorption” and “interleukin 17 (IL-17) signaling pathway.” The PPI network suggested that collagen type III alpha 1 chain (COL3A1) and COL1A2 contribute to the pathogenesis of DCM. Additionally, visualization of the interactions between miRNAs and the hub genes revealed that hsa-miR-5682 and hsa-miR-4500 interacted with both COL3A1 and COL1A2, and thus these miRNAs might play roles in DCM. Immune cell infiltration analysis revealed that DCM patients had more infiltrated plasma cells and fewer infiltrated B memory cells, T follicular helper cells, and resting dendritic cells.CONCLUSIONCOL1A2 and COL3A1 and their targeting miRNAs, hsa-miR-5682 and hsa-miR-4500, may play critical roles in the pathogenesis of DCM, which are closely related to the IL-17 signaling pathway and acute inflammatory response. These results may provide useful clues for the diagnosis and treatment of DCM.

Supervised Machine Learning Models and Protein-Protein Interaction Network Analysis of Gene Expression Profiles Induced by Omega-3 Polyunsaturated Fatty Acids

Background: Omega-3 polyunsaturated fatty acids (PUFAs), such as eicosapentaenoic (EPA) and docosahexaenoic (DHA) acids, have beneficial effects on human health, but their effect on gene expression in elderly individuals (age ≥ 65) is largely unknown. In order to examine this, the gene expression profiles were analyzed in the healthy subjects (n = 96) at baseline and after 26 weeks of supplementation with EPA+DHA to determine up-regulated and down-regulated dif-ferentially expressed genes (DEGs) triggered by PUFAs. The protein-protein interaction (PPI) networks were constructed by mapping these DEGs to a human interactome and linking them to the specific pathways. Objective: This study aimed to implement supervised machine learning models and protein-protein interaction network analysis of gene expression profiles induced by PUFAs. Methods: The transcriptional profile of GSE12375 was obtained from the Gene Expression Om-nibus database, which is based on the Affymetrix NuGO array. The probe cell intensity data were converted into the gene expression values, and the background correction was performed by the multi-array average algorithm. The LIMMA (Linear Models for Microarray Data) algo-rithm was implemented to identify relevant DEGs at baseline and after 26 weeks of supplemen-tation with a p-value &lt; 0.05. The DAVID web server was used to identify and construct the en-riched KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways. Finally, the construction of machine learning (ML) models, including logistic regression, naïve Bayes, and deep neural networks, were implemented for the analyzed DEGs associated with the specific pathways. Results: The results revealed that up-regulated DEGs were associated with neurotrophin/MAPK signaling, whereas the down-regulated DEGs were linked to cancer, acute myeloid leukemia, and long-term depression pathways. Additionally, ML approaches were able to cluster the EPA/DHA-treated and control groups by the logistic regression performing the best. Conclusion: Overall, this study highlights the pivotal changes in DEGs induced by PUFAs and provides the rationale for the implementation of ML algorithms as predictive models for this type of biomedical data.

Transcriptomic Analysis Reveals Genetic Cross-Talk between Periodontitis and Hypothyroidism.

Background Aim of this bioinformatics study based on transcriptomic analysis was to reveal the cross-talk between periodontitis (PD) and hypothyroidism (HT). Methods The gene expression datasets GSE18152 and GSE176153 of HT and GSE10334, GSE16134, and GSE173078 of PD were downloaded through the Gene Expression Omnibus (GEO) database. Differential Expression Genes (DEG) between cases and controls in each microarray were assessed by using the “limma” (linear models for microarray data) R package (|log2 fold change (FC)| >0 and P-value <0.05). To analyze the cross-talk effect between HT and PD, the intersection of DEG of HT and PD was selected. To investigate the biological function of cross-talk genes, the gene ontology (GO) functional enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were applied. Protein-Protein Interaction (PPI) network was constructed using Cytoscape software. Top 10 cross-talk genes were screened, and the expression values of these 10 genes were extracted. ROC analysis was performed by using the pROC package and GGplot2 package of R language to predict the classification accuracy. Results The overlapping DEG between HT and PD were 107 cross-talk genes. The results revealed that developmental process (P-value =1.06E-21) was the most significantly enriched biological process, followed by cell differentiation (P-value =8.49E-18) and immune system process (P-value =6.78E-11). KEGG analysis showed that Complement and coagulation cascades (P-value =2.29E-05), Hematopoietic cell lineage (P-value =2.66E-05), Phospholipase D signaling pathway (P-value =0.034367878) and Chemokine signaling pathway (P-value =0.04946333) were significantly enriched. The top 10 genes with most connections were LCE1B, LCE2B, LCE2A, LCE2C, LCE1C, LCE1F, ITGAM, C1QB, TREM2, and CD19. The AUC values of the two datasets of HT were both greater than 65% (GSE18152 = 81.42%, GSE176153 = 68.75%). AUC values of three datasets of PD were all greater than 60% (GSE10334 = 69.23%, GSE16134 = 73.72%, GSE173078 = 81.6%). Conclusions A genetic cross-talk between HT and PD was detected, whereby LCE family genes appeared to play the most important role.

Identifying and Validating of an Autophagy-Related Gene Signature for the Prediction of Early Relapse in Breast Cancer

BackgroundCompelling evidence has demonstrated the pivotal role of autophagy in the prognosis of breast cancer. Breast cancer (BC) patients with early relapse consistently exhibited worse survival.MethodsThe autophagy-related genes were derived from the Human Autophagy Database (HADb) and high-sequencing data were obtained from The Cancer Genome Atlas (TCGA). Discrepantly expressed autophagy genes (DEAGs) between early relapse and long-term survival groups were performed using the Linear Models for Microarray data (LIMMA) method. Lasso Cox regression analysis was conducted for the selection of the 4-gene autophagy-related gene signature. GSE42568 and GSE21653 databases were enrolled in this study for the external validation of the signature. Then patients were divided into high and low-risk groups based on the specific score formula. GSEA was used to discover the related signaling pathway. The Kaplan-Meier curves and the receiver operating characteristic (ROC) curves were used to evaluate the discrimination and accuracy of the 4-gene signature.ResultsA signature composed of four autophagy-related mRNA including APOL1, HSPA8, SIRT1, and TP73, was identified as significantly associated with the early relapse in BC patients. Time-dependent receiver-operating characteristic at 1 year suggested remarkable accuracy of the signature [area under the curve (AUC = 0.748)]. The risk score model based on the autophagy-related signature showed favorable predicting value in 1-, 2-, and 3-year relapse-free survival (RFS) in training and two validating cohorts. The GSEA displayed gene sets were remarkably enriched in carcinogenic activation pathways and autophagy-related pathways. The nomogram involving three variables (progesterone receptor status, T stage, and 4-gene signature) exhibited relatively good discrimination with a C-index of 0.766.ConclusionsOur study establishes an autophagy-related 4-gene signature that can effectively stratify the high-risk and low-risk BC patients for early relapse. Combined with the clinicopathological variables, the signature could significantly help oncologists tailor more efficient treatment strategies for BC patients.

Epithelial-to-Mesenchymal Transition Mediates Resistance to Maintenance Therapy with Vinflunine in Advanced Urothelial Cell Carcinoma.

Simple SummaryPlatinum-based chemotherapy is the first-line treatment for advanced urothelial cell carcinoma (aUCC). After first-line treatment, we previously showed that maintenance therapy with vinflunine improves progression-free survival. However, some patients are resistant to vinflunine and the specific mechanisms of resistance in aUCC are unclear. We analyzed the genomic landscape and the biological processes potentially related to vinflunine activity and found that epithelial-to-mesenchymal transition (EMT) plays a pivotal role as a resistance mechanism. In experiments with cell lines, curcumin reversed EMT and sensitized cells to vinflunine. We suggest that EMT mediates resistance to vinflunine and that the reversion of this process could enhance the effect of vinflunine in aUCC patients.In the phase II MAJA trial, maintenance therapy with vinflunine resulted in longer progression-free survival compared to best supportive care in advanced urothelial cell carcinoma (aUCC) patients who did not progress after first-line platinum-based chemotherapy. However, despite an initial benefit observed in some patients, unequivocal resistance appears which underlying mechanisms are presently unknown. We have performed gene expression and functional enrichment analyses to shed light on the discovery of these underlying resistance mechanisms. Differential gene expression profile of eight patients with poor outcome and nine with good outcome to vinflunine administered in the MAJA trial were analyzed. RNA was isolated from tumor tissue and gene expression was assessed by microarray. Differential expression was determined with linear models for microarray data. Gene Set Enrichment Analysis (GSEA) was used for the functional classification of the genes. In vitro functional studies were performed using UCC cell lines. Hierarchical clustering showed a differential gene expression pattern between patients with good and poor outcome to vinflunine treatment. GSEA identified epithelial-to-mesenchymal transition (EMT) as the top negatively enriched hallmark in patients with good outcome. In vitro analyses showed that the polyphenol curcumin downregulated EMT markers and sensitized UCC cells to vinflunine. We conclude that EMT mediates resistance to vinflunine and suggest that the reversion of this process could enhance the effect of vinflunine in aUCC patients.

Investigation of gene expression and genetic simultaneous control associated with erectile dysfunction and diabetes.

Diabetes can cause some diseases or abnormalities. One of the disorders caused by diabetes may be erectile dysfunction (ED). ED is sexual dysfunction characterized by the inability to establish or maintain an erect penis during sexual activity and is a common problem of men with chronic type 2 diabetes. These processes, disorders and diseases are highly influenced by the genetics of individuals. In this study, the relationship between genes and diabetes and ED has been explored by a system biology approach. For this purpose, the samples from ten control and diabetic-ED rats were collected. After a search in Gene Expression Omnibus (GEO), series with accession numberGSE2457 comprising of 5 normal and 5 diabetic-ED rats were selected. Raw CEL files of these samples were normalized with robust multi-array average (RMA) expression measure method by using the linear models for microarray data (LIMMA) R package. The extracted probe IDs were transformed into 10451 unique and validated official gene symbols. Then, differentially expressed genes (DEGs) were identified between control and normal penile mucosa by employing the LIMMA R package. DEGs were classified by utilizing KEGG to underlying pathways by Enrichr. The expression values of DEGs were used to construct a gene regulatory network (GRN), by the GENEI3 R package. To analyze the topology of constructed GRNs, betweenness centrality was calculated. Genes with higher betweenness centrality scores were then identified, through the CytoNCA. We then took the commonality of DEGs genes and high-top ranking genes from CytoNCA via a predicted interaction network using GeneMANIAas the most likely important genes in erectile dysfunction. Among the 374 DEGs studied, 146 DEGs showed up-regulation and 228 DEGs displayed down-regulation expression in diabetic-ED rats. According to the Volcano plot, the dpp4, LOC102553868, Ndufa412, Oxct1, Atp2b3 and Zfp91 gene down-regulated and Lpl, Retsat, B4galt1 and Pdk4 genes up-regulated in ED and diabetic rats. Furthermore, genes like dpp4 acted as hubs in the inferred GRN.

Identification of Key eRNAs for Spinal Cord Injury by Integrated Multinomial Bioinformatics Analysis.

Background: Spinal cord injury (SCI) is a severe neurological deficit affecting both young and older people worldwide. The potential role of key enhancer RNAs (eRNAs) in SCI remains elusive, which is a prominent challenge in the trauma repair process. This study aims to investigate the roles of key eRNAs, transcription factors (TFs), signaling pathways, and small-molecule inhibitors in SCI using multi-omics bioinformatics analysis.Methods: Microarray data of peripheral blood mononuclear cell (PBMC) samples from 27 healthy volunteers and 25 chronic-phase SCI patients were retrieved from the Gene Expression Omnibus database. Differentially expressed transcription factors (DETFs), differentially expressed enhancer RNAs (DEeRNAs), and differentially expressed target genes (DETGs) were identified using the Linear Models for Microarray Data (limma) package. Fraction of immune cells was estimated using CIBERSORT algorithm. Gene Set Variation Analysis (GSVA) was applied to identify the downstream signaling pathways. The eRNA regulatory network was constructed based on the correlation results. Connectivity Map (CMap) database was used to find potential drugs for SCI patients. The cellular communication analysis was performed to explore the molecular regulation mechanism of SCI based on single-cell RNA sequencing (scRNA-seq) data. Chromatin immunoprecipitation sequencing (ChIP-seq) and Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq) data were used to validate the key regulatory mechanisms. scRNA-seq dataset was used to validate the cell subtype localization of the key eRNAs.Results: In total, 21 DETFs, 24 DEeRNAs, and 829 DETGs were identified. A regulatory network of 13 DETFs, six DEeRNAs, seven DETGs, two hallmark pathways, two immune cells, and six immune pathways was constructed. The link of Splicing factor proline and glutamine rich (SFPQ) (TF) and vesicular overexpressed in cancer prosurvival protein 1 (VOPP1) (eRNA) (R = 0.990, p < 0.001, positive), VOPP1 (eRNA) and epidermal growth factor receptor (EGFR) (target gene) (R = 0.974, p < 0.001, positive), VOPP1, and T helper (Th) cells (R = −0.987, p < 0.001, negative), and VOPP1 and hallmark coagulation (R = 0.937, p < 0.001, positive) was selected. Trichostatin A was considered the best compound target to SCI-related eRNAs (specificity = 0.471, p < 0.001).Conclusion: VOPP1, upregulated by SFPQ, strengthened the transient expression of EGFR. Th cells and coagulation were the potential downstream pathways of VOPP1. This regulatory network and potential inhibitors provide novel diagnostic biomarkers and therapeutic targets for SCI.

The effect of air pollution on the transcriptomics of the immune response to respiratory infection

Combustion related particulate matter air pollution (PM) is associated with an increased risk of respiratory infections in adults. The exact mechanism underlying this association has not been determined. We hypothesized that increased concentrations of combustion related PM would result in dysregulation of the innate immune system. This epidemiological study includes 111 adult patients hospitalized with respiratory infections who underwent transcriptional analysis of their peripheral blood. We examined the association between gene expression at the time of hospitalization and ambient measurements of particulate air pollutants in the 28 days prior to hospitalization. For each pollutant and time lag, gene-specific linear models adjusting for infection type were fit using LIMMA (Linear Models For Microarray Data), and pathway/gene set analyses were performed using the CAMERA (Correlation Adjusted Mean Rank) program. Comparing patients with viral and/or bacterial infection, the expression patterns associated with air pollution exposure differed. Adjusting for the type of infection, increased concentrations of Delta-C (a marker of biomass smoke) and other PM were associated with upregulation of iron homeostasis and protein folding. Increased concentrations of black carbon (BC) were associated with upregulation of viral related gene pathways and downregulation of pathways related to antigen presentation. The pollutant/pathway associations differed by lag time and by type of infection. This study suggests that the effect of air pollution on the pathogenesis of respiratory infection may be pollutant, timing, and infection specific.

Drug repositioning based on gene expression data for human HER2-positive breast cancer

Human epidermal growth factor receptor 2 (HER2)-positive breast cancer represents approximately 15–30% of all invasive breast cancers. Despite the recent advances in therapeutic practices of HER2 subtype, drug resistance and tumor recurrence still have remained as major problems. Drug discovery is a long and difficult process, so the aim of this study is to find potential new application for existing therapeutic agents. Gene expression data for breast invasive carcinoma were retrieved from The Cancer Genome Atlas (TCGA) database. The normal and tumor samples were analyzed using Linear Models for Microarray Data (LIMMA) R package in order to find the differentially expressed genes (DEGs). These genes were used as entry for the library of integrated network-based cellular signatures (LINCS) L1000CDS2 software and suggested 24 repurposed drugs. According to the obtained results, some of these drugs including vorinostat, mocetinostat, alvocidib, CGP-60474, BMS-387032, AT-7519, and curcumin have significant functional similarity and structural correlation with FDA-approved breast cancer drugs. Based on the drug-target network, which consisted of the repurposed drugs and their target genes, the aforementioned drugs had the highest degrees. Moreover, the experimental approach verified curcumin as an effective therapeutic agent for HER2 positive breast cancer. Hence, our work suggested that some repurposed drugs based on gene expression data can be noticed as potential drugs for the treatment of HER2-positive breast cancer.

Clustering of genes from microarray data using hierarchical projective adaptive resonance theory: a case study of tuberculosis

We propose the hierarchical Projective Adaptive Resonance Theory (PART) algorithm for classification of gene expression data. This algorithm is realized by combing transposed quasi-supervised PART and unsupervised PART. We develop the corresponding validation statistics for each process and compare it with other clustering algorithms in a case study of tuberculosis (TB). First, we use sample-based transposed quasi-supervised PART to obtain optimal clustering results of samples distinguished by time post-infection and the representative genes for each cluster including up-regulated, down-regulated and stable genes. The up- and down-regulated genes show more than 90% similarity to the result derived from Linear Models for Microarray Data and are verified by weighted k-nearest neighbor model on TB projection. Second, we use gene-based unsupervised PART algorithm to cluster these representative genes where functional enrichment analysis is conducted in each cluster. We further confirm the main immune response of human macrophage-like THP-1 cells against TB within 2 days is type I interferon-mediated innate immunity. This study demonstrates how hierarchical PART algorithm analyzes microarray data. The sample-based quasi-supervised PART extracts representative genes and narrows down the shortlist of disease-relevant genes and gene-based unsupervised PART classifies representative genes that help to interpret immune response against TB.

ProteoSign v2: a faster and evolved user-friendly online tool for statistical analyses of differential proteomics.

Bottom-up proteomics analyses have been proved over the last years to be a powerful tool in the characterization of the proteome and are crucial for understanding cellular and organism behaviour. Through differential proteomic analysis researchers can shed light on groups of proteins or individual proteins that play key roles in certain, normal or pathological conditions. However, several tools for the analysis of such complex datasets are powerful, but hard-to-use with steep learning curves. In addition, some other tools are easy to use, but are weak in terms of analytical power. Previously, we have introduced ProteoSign, a powerful, yet user-friendly open-source online platform for protein differential expression/abundance analysis designed with the end-proteomics user in mind. Part of Proteosign's power stems from the utilization of the well-established Linear Models For Microarray Data (LIMMA) methodology. Here, we present a substantial upgrade of this computational resource, called ProteoSign v2, where we introduce major improvements, also based on user feedback. The new version offers more plot options, supports additional experimental designs, analyzes updated input datasets and performs a gene enrichment analysis of the differentially expressed proteins. We also introduce the deployment of the Docker technology and significantly increase the speed of a full analysis. ProteoSign v2 is available at http://bioinformatics.med.uoc.gr/ProteoSign.