Key genes and co-expression modules involved in asthma pathogenesis

Transcriptomic and physiological responses of Quercus acutissima and Quercus palustris to drought stress and rewatering

Establishment of oak seedlings, which is an important factor in forest restoration, is affected by drought that hampers the survival, growth, and development of seedlings. Therefore, it is necessary to understand how seedlings respond to and recover from water-shortage stress. We subjected seedlings of two oak species, Quercus acutissima and Quercus palustris, to drought stress for one month and then rewatered them for six days to observe physiological and genetic expression changes. Phenotypically, the growth of Q. acutissima was reduced and severe wilting and recovery failure were observed in Q. palustris after an increase in plant temperature. The two species differed in several physiological parameters during drought stress and recovery. Although the photosynthesis-related indicators did not change in Q. acutissima, they were decreased in Q. palustris. Moreover, during drought, content of soluble sugars was significantly increased in both species, but it recovered to original levels only in Q. acutissima. Malondialdehyde content increased in both the species during drought, but it did not recover in Q. palustris after rewatering. Among the antioxidant enzymes, only superoxide dismutase activity increased in Q. acutissima during drought, whereas activities of ascorbate peroxidase, catalase, and glutathione reductase increased in Q. palustris. Abscisic acid levels were increased and then maintained in Q. acutissima, but recovered to previous levels after rewatering in Q. palustris. RNA samples from the control, drought, recovery day 1, and recovery day 6 treatment groups were compared using transcriptome analysis. Q. acutissima exhibited 832 and 1076 differentially expressed genes (DEGs) related to drought response and recovery, respectively, whereas Q. palustris exhibited 3947 and 1587 DEGs, respectively under these conditions. Gene ontology enrichment of DEGs revealed "response to water," "apoplast," and "Protein self-association" to be common to both the species. However, in the heatmap analysis of genes related to sucrose and starch synthesis, glycolysis, antioxidants, and hormones, the two species exhibited very different transcriptome responses. Nevertheless, the levels of most DEGs returned to their pre-drought levels after rewatering. These results provide a basic foundation for understanding the physiological and genetic expression responses of oak seedlings to drought stress and recovery.

Unraveling the link between PTBP1 and severe asthma through machine learning and association rule mining method

Severe asthma is a chronic inflammatory airway disease with great therapeutic challenges. Understanding the genetic and molecular mechanisms of severe asthma may help identify therapeutic strategies for this complex condition. RNA expression data were analyzed using a combination of artificial intelligence methods to identify novel genes related to severe asthma. Through the ANOVA feature selection approach, 100 candidate genes were selected among 54,715 mRNAs in blood samples of patients with severe asthmatic and healthy groups. A deep learning model was used to validate the significance of the candidate genes. The accuracy, F1-score, AUC-ROC, and precision of the 100 genes were 83%, 0.86, 0.89, and 0.9, respectively. To discover hidden associations among selected genes, association rule mining was applied. The top 20 genes including the PTBP1, RAB11FIP3, APH1A, and MYD88 were recognized as the most frequent items among severe asthma association rules. The PTBP1 was found to be the most frequent gene associated with severe asthma among those 20 genes. PTBP1 was the gene most frequently associated with severe asthma among candidate genes. Identification of master genes involved in the initiation and development of asthma can offer novel targets for its diagnosis, prognosis, and targeted-signaling therapy.

Comprehensive analysis of the progression mechanisms of CRPC and its inhibitor discovery based on machine learning algorithms

Background: Almost all patients treated with androgen deprivation therapy (ADT) eventually develop castration-resistant prostate cancer (CRPC). Our research aims to elucidate the potential biomarkers and molecular mechanisms that underlie the transformation of primary prostate cancer into CRPC. Methods: We collected three microarray datasets (GSE32269, GSE74367, and GSE66187) from the Gene Expression Omnibus (GEO) database for CRPC. Differentially expressed genes (DEGs) in CRPC were identified for further analyses, including Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and gene set enrichment analysis (GSEA). Weighted gene coexpression network analysis (WGCNA) and two machine learning algorithms were employed to identify potential biomarkers for CRPC. The diagnostic efficiency of the selected biomarkers was evaluated based on gene expression level and receiver operating characteristic (ROC) curve analyses. We conducted virtual screening of drugs using AutoDock Vina. In vitro experiments were performed using the Cell Counting Kit-8 (CCK-8) assay to evaluate the inhibitory effects of the drugs on CRPC cell viability. Scratch and transwell invasion assays were employed to assess the effects of the drugs on the migration and invasion abilities of prostate cancer cells. Results: Overall, a total of 719 DEGs, consisting of 513 upregulated and 206 downregulated genes, were identified. The biological functional enrichment analysis indicated that DEGs were mainly enriched in pathways related to the cell cycle and metabolism. CCNA2 and CKS2 were identified as promising biomarkers using a combination of WGCNA, LASSO logistic regression, SVM-RFE, and Venn diagram analyses. These potential biomarkers were further validated and exhibited a strong predictive ability. The results of the virtual screening revealed Aprepitant and Dolutegravir as the optimal targeted drugs for CCNA2 and CKS2, respectively. In vitro experiments demonstrated that both Aprepitant and Dolutegravir exerted significant inhibitory effects on CRPC cells (p < 0.05), with Aprepitant displaying a superior inhibitory effect compared to Dolutegravir. Discussion: The expression of CCNA2 and CKS2 increases with the progression of prostate cancer, which may be one of the driving factors for the progression of prostate cancer and can serve as diagnostic biomarkers and therapeutic targets for CRPC. Additionally, Aprepitant and Dolutegravir show potential as anti-tumor drugs for CRPC.

Identifying biomarkers of differential chemotherapy response in TNBC patient-derived xenografts with a CTD/WGCNA approach

Although systemic chemotherapy remains the standard of care for TNBC, even combination chemotherapy is often ineffective. The identification of biomarkers for differential chemotherapy response would allow for the selection of responsive patients, thus maximizing efficacy and minimizing toxicities. Here, we leverage TNBC PDXs to identify biomarkers of response. To demonstrate their ability to function as a preclinical cohort, PDXs were characterized using DNA sequencing, transcriptomics, and proteomics to show consistency with clinical samples. We then developed a network-based approach (CTD/WGCNA) to identify biomarkers of response to carboplatin (MSI1, TMSB15A, ARHGDIB, GGT1, SV2A, SEC14L2, SERPINI1, ADAMTS20, DGKQ) and docetaxel (c, MAGED4, CERS1, ST8SIA2, KIF24, PARPBP). CTD/WGCNA multigene biomarkers are predictive in PDX datasets (RNAseq and Affymetrix) for both taxane- (docetaxel or paclitaxel) and platinum-based (carboplatin or cisplatin) response, thereby demonstrating cross-expression platform and cross-drug class robustness. These biomarkers were also predictive in clinical datasets, thus demonstrating translational potential.

Establishment and validation of a predictive model of preeclampsia based on transcriptional signatures of 43 genes in decidua basalis and peripheral blood

Preeclampsia (PE) has an increasing incidence worldwide, and there is no gold standard for prediction. Recent progress has shown that abnormal decidualization and impaired vascular remodeling are essential to PE pathogenesis. Therefore, it is of great significance to analyze the decidua basalis and blood changes of PE to explore new methods. Here, we performed weighted gene co-expression network analysis based on 9553 differentially expressed genes of decidua basalis data (GSE60438 includes 25 cases of PE and 23 non-cases) from Gene Expression Omnibus to screen relevant module-eigengenes (MEs). Among them, MEblue and MEgrey are the most correlated with PE, which contains 371 core genes. Subsequently, we applied the logistic least absolute shrinkage and selection operator regression, screened 43 genes most relevant to prediction from the intersections of the 371 genes and training set (GSE48424 includes 18 cases of PE and 18 non-cases) genes, and built a predictive model. The specificity and sensitivity are illustrated by receiver operating characteristic curves, and the stability was verified by two validation sets (GSE86200 includes 12 cases of PE and 48 non-cases, and GSE85307 includes 47 cases of PE and 110 non-cases). The results demonstrated that our predictive model shows good predictions, with an area under the curve of 0.991 for the training set, 0.874 and 0.986 for the validation sets. Finally, we found the 43 key marker genes in the model are closely associated with the clinically accepted predictive molecules, including FLT1, PIGF, ENG and VEGF. Therefore, this predictive model provides a potential approach for PE diagnosis and treatment.

Network-Based Analysis Reveals Novel Biomarkers in Peripheral Blood of Patients With Preeclampsia

WGCNA is a potent systems biology approach that explains the connection of gene expression based on a microarray database, which facilitates the discovery of disease therapy targets or potential biomarkers. Preeclampsia is a kind of pregnancy-induced hypertension caused by complex factors. The disease’s pathophysiology, however, remains unknown. The focus of this research is to utilize WGCNA to identify susceptible modules and genes in the peripheral blood of preeclampsia patients. Obtain the whole gene expression data of GSE48424 preeclampsia patients and normal pregnant women from NCBI’s GEO database. WGCNA is used to construct a gene co-expression network by calculating correlation coefficients between modules and phenotypic traits, screening important modules, and filtering central genes. To identify hub genes, we performed functional enrichment analysis, pathway analysis, and protein-protein interaction (PPI) network construction on key genes in critical modules. Then, the genetic data file GSE149437 and clinical peripheral blood samples were used as a validation cohort to determine the diagnostic value of these key genes. Nine gene co-expression modules were constructed through WGCNA analysis. Among them, the blue module is significantly related to preeclampsia and is related to its clinical severity. Thirty genes have been discovered by using the intersection of the genes in the blue module and the DEGs genes as the hub genes. It was found that HDC, MS4A2, and SLC18A2 scored higher in the PPI network and were identified as hub genes. These three genes were also differentially expressed in peripheral blood validation samples. Based on the above three genes, we established the prediction model of peripheral blood markers of preeclampsia and drew the nomogram and calibration curve. The ROC curves were used in the training cohort GSE48424 and the validation cohort GSE149437 to verify the predictive value of the above model. Finally, it was confirmed in the collected clinical peripheral blood samples that MS4A2 was differentially expressed in the peripheral blood of early-onset and late-onset preeclampsia, which is of great significance. This study provides a new biomarker and prediction model for preeclampsia.

Metabolomic Associations of Asthma in the Hispanic Community Health Study/Study of Latinos

Asthma disproportionally affects Hispanic and/or Latino backgrounds; however, the relation between circulating metabolites and asthma remains unclear. We conducted a cross-sectional study associating 640 individual serum metabolites, as well as twelve metabolite modules, with asthma in 3347 Hispanic/Latino background participants (514 asthmatics, 15.36%) from the Hispanic/Latino Community Health Study/Study of Latinos. Using survey logistic regression, per standard deviation (SD) increase in 1-arachidonoyl-GPA (20:4) was significantly associated with 32% high odds of asthma after accounting for clinical risk factors (p = 6.27 × 10−5), and per SD of the green module, constructed using weighted gene co-expression network, was suggestively associated with 25% high odds of asthma (p = 0.006). In the stratified analyses by sex and Hispanic and/or Latino backgrounds, the effect of 1-arachidonoyl-GPA (20:4) and the green module was predominantly observed in women (OR = 1.24 and 1.37, p < 0.001) and people of Cuban and Puerto-Rican backgrounds (OR = 1.25 and 1.27, p < 0.01). Mutations in Fatty Acid Desaturase 2 (FADS2) affected the levels of 1-arachidonoyl-GPA (20:4), and Mendelian Randomization analyses revealed that high genetically regulated 1-arachidonoyl-GPA (20:4) levels were associated with increased odds of asthma (p < 0.001). The findings reinforce a molecular basis for asthma etiology, and the potential causal effect of 1-arachidonoyl-GPA (20:4) on asthma provides an opportunity for future intervention.

Analysis of potential genetic biomarkers and molecular mechanism of smoking-related postmenopausal osteoporosis using weighted gene co-expression network analysis and machine learning

OBJECTIVES

Smoking is a significant independent risk factor for postmenopausal osteoporosis, leading to genome variations in postmenopausal smokers. This study investigates potential biomarkers and molecular mechanisms of smoking-related postmenopausal osteoporosis (SRPO).

MATERIALS AND METHODS

The GSE13850 microarray dataset was downloaded from Gene Expression Omnibus (GEO). Gene modules associated with SRPO were identified using weighted gene co-expression network analysis (WGCNA), protein-protein interaction (PPI) analysis, and pathway and functional enrichment analyses. Feature genes were selected using two machine learning methods: support vector machine-recursive feature elimination (SVM-RFE) and random forest (RF). The diagnostic efficiency of the selected genes was assessed by gene expression analysis and receiver operating characteristic curve.

RESULTS

Eight highly conserved modules were detected in the WGCNA network, and the genes in the module that was strongly correlated with SRPO were used for constructing the PPI network. A total of 113 hub genes were identified in the core network using topological network analysis. Enrichment analysis results showed that hub genes were closely associated with the regulation of RNA transcription and translation, ATPase activity, and immune-related signaling. Six genes (HNRNPC, PFDN2, PSMC5, RPS16, TCEB2, and UBE2V2) were selected as genetic biomarkers for SRPO by integrating the feature selection of SVM-RFE and RF.

CONCLUSION

The present study identified potential genetic biomarkers and provided a novel insight into the underlying molecular mechanism of SRPO.

Redox Role of ROS and Inflammation in Pulmonary Diseases

Reactive oxygen species (ROS), either derived from exogenous sources or overproduced endogenously, can disrupt the body's antioxidant defenses leading to compromised redox homeostasis. The lungs are highly susceptible to ROS-mediated damage. Oxidative stress (OS) caused by this redox imbalance leads to the pathogenesis of multiple pulmonary diseases such as asthma, chronic obstructive pulmonary disease (COPD), and acute respiratory distress syndrome (ARDS). OS causes damage to important cellular components in terms of lipid peroxidation, protein oxidation, and DNA histone modification. Inflammation further enhances ROS production inducing changes in transcriptional factors which mediate cellular stress response pathways. This deviation from normal cell function contributes to the detrimental pathological characteristics often seen in pulmonary diseases. Although antioxidant therapies are feasible approaches in alleviating OS-related lung impairment, a comprehensive understanding of the updated role of ROS in pulmonary inflammation is vital for the development of optimal treatments. In this chapter, we review the major pulmonary diseases-including COPD, asthma, ARDS, COVID-19, and lung cancer-as well as their association with ROS.

Probabilistic domain-knowledge modeling of disorder pathogenesis for dynamics forecasting of acute onset

Disease pathogenesis, a type of domain knowledge about biological mechanisms leading to diseases, has not been adequately encoded in machine-learning-based medical diagnostic models because of the inter-patient variabilities and complex dependencies of the underlying pathogenetic mechanisms. We propose 1) a novel pathogenesis probabilistic graphical model (PPGM) to quantify the dynamics underpinning patient-specific data and pathogenetic domain knowledge, 2) a Bayesian-based inference paradigm to answer the medical queries and forecast acute onsets. The PPGM model consists of two components: a Bayesian network of patient attributes and a temporal model of pathogenetic mechanisms. The model structure was reconstructed from expert knowledge elicitation, and its parameters were estimated using Variational Expectation-Maximization algorithms. We benchmarked our model with two well-established hidden Markov models (HMMs) - Input-output HMM (IO-HMM) and Switching Auto-Regressive HMM (SAR-HMM) - to evaluate the computational costs, forecasting performance, and execution time. Two case studies on Obstructive Sleep Apnea (OSA) and Paroxysmal Atrial Fibrillation (PAF) were used to validate the model. While the performance of the parameter learning step was equivalent to those of IO-HMM and SAR-HMM models, our model forecasting ability was outperforming those two models. The merits of the PPGM model are its representation capability to capture the dynamics of pathogenesis and perform medical inferences and its interpretability for physicians. The model has been used to perform medical queries and forecast the acute onset of OSA and PAF. Additional applications of the model include prognostic healthcare and preventive personalized treatments.

Stratification of Estrogen Receptor-Negative Breast Cancer Patients by Integrating the Somatic Mutations and Transcriptomic Data

Patients with estrogen receptor-negative breast cancer generally have a worse prognosis than estrogen receptor-positive patients. Nevertheless, a significant proportion of the estrogen receptor-negative cases have favorable outcomes. Identifying patients with a good prognosis, however, remains difficult, as recent studies are quite limited. The identification of molecular biomarkers is needed to better stratify patients. The significantly mutated genes may be potentially used as biomarkers to identify the subtype and to predict outcomes. To identify the biomarkers of receptor-negative breast cancer among the significantly mutated genes, we developed a workflow to screen significantly mutated genes associated with the estrogen receptor in breast cancer by a gene coexpression module. The similarity matrix was calculated with distance correlation to obtain gene modules through a weighted gene coexpression network analysis. The modules highly associated with the estrogen receptor, called important modules, were enriched for breast cancer-related pathways or disease. To screen significantly mutated genes, a new gene list was obtained through the overlap of the important module genes and the significantly mutated genes. The genes on this list can be used as biomarkers to predict survival of estrogen receptor-negative breast cancer patients. Furthermore, we selected six hub significantly mutated genes in the gene list which were also able to separate these patients. Our method provides a new and alternative method for integrating somatic gene mutations and expression data for patient stratification of estrogen receptor-negative breast cancers.

Novel Comprehensive Bioinformatics Approaches to Determine the Molecular Genetic Susceptibility Profile of Moderate and Severe Asthma

Background: Asthma is a chronic inflammatory condition linked to hyperresponsiveness in the airways. There is currently no cure available for asthma, and therapy choices are limited. Asthma is the result of the interplay between genes and the environment. The exact molecular genetic mechanism of asthma remains elusive. Aims: The aim of this study is to provide a comprehensive, detailed molecular etiology profile for the molecular factors that regulate the severity of asthma and pathogenicity using integrative bioinformatics tools. Methods: The GSE43696 omnibus gene expression dataset, which contains 50 moderate cases, 38 severe cases, and 20 healthy controls, was used to investigate differentially expressed genes (DEGs), susceptible chromosomal loci, gene networks, pathways, gene ontologies, and protein–protein interactions (PPIs) using an intensive bioinformatics pipeline. Results: The PPI network analysis yielded DEGs that contribute to interactions that differ from moderate-to-severe asthma. The combined interaction scores resulted in higher interactions for the genes STAT3, AGO2, COL1A1, CLCN6, and KSR for moderate asthma and JAK2, INSR, ERBB2, NR3C1, and PTK6 for severe asthma. Enrichment analysis (EA) demonstrated differential enrichment between moderate and severe asthma phenotypes; the ion transport regulation pathway was significantly enhanced in severe asthma phenotypes compared to that in moderate asthma phenotypes and involved PER2, GCR, IRS-2, KCNK7, KCNK6, NOX1, and SCN7A. The most enriched common pathway in both moderate and severe asthma is the development of the glucocorticoid receptor (GR) signaling pathway followed by glucocorticoid-mediated inhibition of proinflammatory and proconstrictory signaling in the airway of smooth muscle cell pathways. Gene sets were shared between severe and moderate asthma at 16 chromosome locations, including 17p13.1, 16p11.2, 17q21.31, 1p36, and 19q13.2, while 60 and 48 chromosomal locations were unique for both moderate and severe asthma, respectively. Phylogenetic analysis for DEGs showed that several genes have been intersected in phases of asthma in the same cluster of genes. This could indicate that several asthma-associated genes have a common ancestor and could be linked to the same biological function or gene family, implying the importance of these genes in the pathogenesis of asthma. Conclusion: New genetic risk factors for the development of moderate-to-severe asthma were identified in this study, and these could provide a better understanding of the molecular pathology of asthma and might provide a platform for the treatment of asthma.