Identification and prioritization of gene sets associated with schizophrenia risk by co-expression network analysis in human brain

Schizophrenia-Like Deficits and Impaired Glutamate/Gamma-aminobutyric acid Homeostasis in Zfp804a Conditional Knockout Mice

BACKGROUND AND HYPOTHESIS

Zinc finger protein 804A (ZNF804A) was the first genome-wide associated susceptibility gene for schizophrenia (SCZ) and played an essential role in the pathophysiology of SCZ by influencing neurodevelopment regulation, neurite outgrowth, synaptic plasticity, and RNA translational control; however, the exact molecular mechanism remains unclear.

STUDY DESIGN

A nervous-system-specific Zfp804a (ZNF804A murine gene) conditional knockout (cKO) mouse model was generated using clustered regularly interspaced short palindromic repeat/Cas9 technology and the Cre/loxP method.

RESULTS

Multiple and complex SCZ-like behaviors, such as anxiety, depression, and impaired cognition, were observed in Zfp804a cKO mice. Molecular biological methods and targeted metabolomics assay validated that Zfp804a cKO mice displayed altered SATB2 (a cortical superficial neuron marker) expression in the cortex; aberrant NeuN, cleaved caspase 3, and DLG4 (markers of mature neurons, apoptosis, and postsynapse, respectively) expressions in the hippocampus and a loss of glutamate (Glu)/γ-aminobutyric acid (GABA) homeostasis with abnormal GAD67 (Gad1) expression in the hippocampus. Clozapine partly ameliorated some SCZ-like behaviors, reversed the disequilibrium of the Glu/GABA ratio, and recovered the expression of GAD67 in cKO mice.

CONCLUSIONS

Zfp804a cKO mice reproducing SCZ-like pathological and behavioral phenotypes were successfully developed. A novel mechanism was determined in which Zfp804a caused Glu/GABA imbalance and reduced GAD67 expression, which was partly recovered by clozapine treatment. These findings underscore the role of altered gene expression in understanding the pathogenesis of SCZ and provide a reliable SCZ model for future therapeutic interventions and biomarker discovery.

Network-wide risk convergence in gene co-expression identifies reproducible genetic hubs of schizophrenia risk

The omnigenic model posits that genetic risk for traits with complex heritability involves cumulative effects of peripheral genes on mechanistic "core genes," suggesting that in a network of genes, those closer to clusters including core genes should have higher GWAS signals. In gene co-expression networks, we confirmed that GWAS signals accumulate in genes more connected to risk-enriched gene clusters, highlighting across-network risk convergence. This was strongest in adult psychiatric disorders, especially schizophrenia (SCZ), spanning 70% of network genes, suggestive of super-polygenic architecture. In snRNA-seq cell type networks, SCZ risk convergence was strongest in L2/L3 excitatory neurons. We prioritized genes most connected to SCZ-GWAS genes, which showed robust association to a CRISPRa measure of PGC3 regulation and were consistently identified across several brain regions. Several genes, including dopamine-associated ones, were prioritized specifically in the striatum. This strategy thus retrieves current drug targets and can be used to prioritize other potential drug targets.

Promotion of apoptosis in melanoma cells by taxifolin through the PI3K/AKT signaling pathway: Screening of natural products using WGCNA and CMAP platforms

Melanoma is a skin cancer originating from melanocytes. The global incidence rate of melanoma is rapidly increasing, posing significant public health challenges. Identifying effective therapeutic agents is crucial in addressing this growing problem. Natural products have demonstrated promising anti-tumor activity. In this study, a plant flavonoid, taxifolin, was screened using Weighted Correlation Network Analysis (WGCNA) in combination with the Connectivity Map (CMAP) platform. Taxifolin was confirmed to inhibit the proliferation, migration, and invasion ability of melanoma A375 and MV-3 cells by promoting apoptosis. Additionally, it suppressed the Epithelial-Mesenchymal Transition (EMT) process of melanoma cells. Cyber pharmacological analysis revealed that taxifolin exerts its inhibitory effect on melanoma through the PI3K/AKT signaling pathway, specifically by downregulating the protein expression of p-PI3K and p-AKT. Notably, the addition of SC-79, an activator of the PI3K/AKT signaling pathway, reversed the effects of taxifolin on cell migration and apoptosis. Furthermore, in vivo experiments demonstrated that taxifolin treatment slowed tumor growth in mice without significant toxic effects. Based on these findings, taxifolin holds promise as a potential drug for melanoma treatment.

Interaction between RNF4 and SART3 is associated with the risk of schizophrenia

The pathogenesis of schizophrenia (SCZ) is heavily influenced by genetic factors. Ring finger protein 4 (RNF4) and squamous cell carcinoma antigen recognized by T cells 3 (SART3) are thought to be involved in nervous system growth and development via oxidative stress pathways. Moreover, they have previously been linked to SCZ. Yet the role of RNF4 and SART3 in SCZ remains unclear. Here, we investigated how these two genes are involved in SCZ by studying their variants observed in patients. We first observed significantly elevated mRNA levels of RNF4 and SART3 in the peripheral blood in both first-episode (n = 30) and chronic (n = 30) SCZ patients compared to controls (n = 60). Next, we targeted-sequenced three single nucleotide polymorphisms (SNPs) in SART3 and six SNPs in RNF4 for association with SCZ using the genomic DNA extracted from peripheral blood leukocytes from SCZ participants (n = 392) and controls (n = 572). We observed a combination of SNPs that included rs1203860, rs2282765 (both in RNF4), and rs2287550 (in SART3) was associated with increased risk of SCZ, suggesting common pathogenic mechanisms between these two genes. We then conducted experiments in HEK293T cells to better understand the interaction between RNF4 and SART3. We observed that SART3 lowered the expression of RNF4 through ubiquitination and downregulated the expression of nuclear factor E2-related factor 2 (NRF2), a downstream factor of RNF4, implicating the existence of a possible shared regulatory mechanism for RNF4 and SART3. In conclusion, our study provides evidence that the interaction between RNF4 and SART3 contributes to the risk of SCZ. The findings shed light on the underlying molecular mechanisms of SCZ and may lead to the development of new therapies and interventions for this disorder.

Polygenic and transcriptional risk scores identify chronic obstructive pulmonary disease subtypes

Rationale

Genetic variants and gene expression predict risk of chronic obstructive pulmonary disease (COPD), but their effect on COPD heterogeneity is unclear.

Objectives

Define high-risk COPD subtypes using both genetics (polygenic risk score, PRS) and blood gene expression (transcriptional risk score, TRS) and assess differences in clinical and molecular characteristics.

Methods

We defined high-risk groups based on PRS and TRS quantiles by maximizing differences in protein biomarkers in a COPDGene training set and identified these groups in COPDGene and ECLIPSE test sets. We tested multivariable associations of subgroups with clinical outcomes and compared protein-protein interaction networks and drug repurposing analyses between high-risk groups.

Measurements and Main Results

We examined two high-risk omics-defined groups in non-overlapping test sets (n=1,133 NHW COPDGene, n=299 African American (AA) COPDGene, n=468 ECLIPSE). We defined "High activity" (low PRS/high TRS) and "severe risk" (high PRS/high TRS) subgroups. Participants in both subgroups had lower body-mass index (BMI), lower lung function, and alterations in metabolic, growth, and immune signaling processes compared to a low-risk (low PRS, low TRS) reference subgroup. "High activity" but not "severe risk" participants had greater prospective FEV 1 decline (COPDGene: -51 mL/year; ECLIPSE: - 40 mL/year) and their proteomic profiles were enriched in gene sets perturbed by treatment with 5-lipoxygenase inhibitors and angiotensin-converting enzyme (ACE) inhibitors.

Conclusions

Concomitant use of polygenic and transcriptional risk scores identified clinical and molecular heterogeneity amongst high-risk individuals. Proteomic and drug repurposing analysis identified subtype-specific enrichment for therapies and suggest prior drug repurposing failures may be explained by patient selection.

Dopamine signaling enriched striatal gene set predicts striatal dopamine synthesis and physiological activity in vivo

The polygenic architecture of schizophrenia implicates several molecular pathways involved in synaptic function. However, it is unclear how polygenic risk funnels through these pathways to translate into syndromic illness. Using tensor decomposition, we analyze gene co-expression in the caudate nucleus, hippocampus, and dorsolateral prefrontal cortex of post-mortem brain samples from 358 individuals. We identify a set of genes predominantly expressed in the caudate nucleus and associated with both clinical state and genetic risk for schizophrenia that shows dopaminergic selectivity. A higher polygenic risk score for schizophrenia parsed by this set of genes predicts greater dopamine synthesis in the striatum and greater striatal activation during reward anticipation. These results translate dopamine-linked genetic risk variation into in vivo neurochemical and hemodynamic phenotypes in the striatum that have long been implicated in the pathophysiology of schizophrenia.

Weighted gene co-expression network analysis reveals key biomarkers and immune infiltration characteristics for bronchial epithelial cells from asthmatic patients

BACKGROUND

Asthma ranks among the most prevalent non-communicable diseases worldwide. Previous studies have elucidated the significant role of the immune system in its pathophysiology. Nevertheless, the immune-related mechanisms underlying asthma are complex and still inadequately understood. Thus, our objective was to investigate novel key biomarkers and immune infiltration characteristics associated with asthma by employing integrated bioinformatics tools.

METHODS

In this study, we conducted a weighted gene co-expression network analysis (WGCNA) to identify key modules and genes potentially implicated in asthma. Functional annotation of these key modules and genes was carried out through gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. Additionally, we constructed a protein-protein interaction (PPI) network using the STRING database to identify 10 hub genes. Furthermore, we evaluated the relative proportion of immune cells in bronchial epithelial cell samples from 20 healthy individuals and 88 asthmatic patients using CIBERSORT. Finally, we validated the hub genes and explored their correlation with immune infiltration.

RESULTS

Furthermore, 20 gene expression modules and 10 hub genes were identified herein. Among them, complement component 3 (C3), prostaglandin I2 receptor (PTGIR), parathyroid hormone-like hormone (PTHLH), and C-X3-C motif chemokine ligand 1 (CX3CL1) were closely correlated with the infiltration of immune cells. They may be novel candidate biomarkers or therapeutic targets for asthma. Furthermore, B cells memory, and plasma cells might play an important role in immune cell infiltration after asthma.

CONCLUSIONS

C3, PTGIR, CX3CL1, and PTHLH have important clinical diagnostic values and are correlated with infiltration of multiple immune cell types in asthma. These hub genes, B cells memory, and plasma cells may become important biological targets for therapeutic asthma drug screening and drug design.

The Genomic Intersection of Oligodendrocyte Dynamics in Schizophrenia and Aging Unravels Novel Pathological Mechanisms and Therapeutic Potentials

Schizophrenia is a significant worldwide health concern, affecting over 20 million individuals and contributing to a potential reduction in life expectancy by up to 14.5 years. Despite its profound impact, the precise pathological mechanisms underlying schizophrenia continue to remain enigmatic, with previous research yielding diverse and occasionally conflicting findings. Nonetheless, one consistently observed phenomenon in brain imaging studies of schizophrenia patients is the disruption of white matter, the bundles of myelinated axons that provide connectivity and rapid signalling between brain regions. Myelin is produced by specialised glial cells known as oligodendrocytes, which have been shown to be disrupted in post-mortem analyses of schizophrenia patients. Oligodendrocytes are generated throughout life by a major population of oligodendrocyte progenitor cells (OPC), which are essential for white matter health and plasticity. Notably, a decline in a specific subpopulation of OPC has been identified as a principal factor in oligodendrocyte disruption and white matter loss in the aging brain, suggesting this may also be a factor in schizophrenia. In this review, we analysed genomic databases to pinpoint intersections between aging and schizophrenia and identify shared mechanisms of white matter disruption and cognitive dysfunction.

Integrative Pharmacology in the Treatment of Substance Use Disorders

The detrimental physical, mental, and socioeconomic effects of substance use disorders (SUDs) have been apparent to the medical community for decades. However, it has become increasingly urgent in recent years to develop novel pharmacotherapies to treat SUDs. Currently, practitioners typically rely on monotherapy. Monotherapy has been shown to be superior to no treatment at all for most substance classes. However, many randomized controlled trials (RCTs) have revealed that monotherapy leads to poorer outcomes when compared with combination treatment in all specialties of medicine. The results of RCTs suggest that monotherapy frequently fails since multiple dysregulated pathways, enzymes, neurotransmitters, and receptors are involved in the pathophysiology of SUDs. As such, research is urgently needed to determine how various neurobiological mechanisms can be targeted by novel combination treatments to create increasingly specific yet exceedingly comprehensive approaches to SUD treatment. This article aims to review the neurobiology that integrates many pathophysiologic mechanisms and discuss integrative pharmacology developments that may ultimately improve clinical outcomes for patients with SUDs. Many neurobiological mechanisms are known to be involved in SUDs including dopaminergic, nicotinic, N-methyl-D-aspartate (NMDA), and kynurenic acid (KYNA) mechanisms. Emerging evidence indicates that KYNA, a tryptophan metabolite, modulates all these major pathophysiologic mechanisms. Therefore, achieving KYNA homeostasis by harmonizing integrative pathophysiology and pharmacology could prove to be a better therapeutic approach for SUDs. We propose KYNA-NMDA-α7nAChRcentric pathophysiology, the "conductor of the orchestra," as a novel approach to treat many SUDs concurrently. KYNA-NMDA-α7nAChR pathophysiology may be the "command center" of neuropsychiatry. To date, extant RCTs have shown equivocal findings across comparison conditions, possibly because investigators targeted single pathophysiologic mechanisms, hit wrong targets in underlying pathophysiologic mechanisms, and tested inadequate monotherapy treatment. We provide examples of potential combination treatments that simultaneously target multiple pathophysiologic mechanisms in addition to KYNA. Kynurenine pathway metabolism demonstrates the greatest potential as a target for neuropsychiatric diseases. The investigational medications with the most evidence include memantine, galantamine, and N-acetylcysteine. Future RCTs are warranted with novel combination treatments for SUDs. Multicenter RCTs with integrative pharmacology offer a promising, potentially fruitful avenue to develop novel therapeutics for the treatment of SUDs.

Exploring the multifaceted role of RASGRP1 in disease: immune, neural, metabolic, and oncogenic perspectives

RAS guanyl releasing protein 1 (RASGRP1) is a guanine nucleotide exchange factor (GEF) characterized by the presence of a RAS superfamily GEF domain. It functions as a diacylglycerol (DAG)-regulated nucleotide exchange factor, specifically activating RAS through the exchange of bound GDP for GTP. Activation of RAS by RASGRP1 has a wide range of downstream effects at the cellular level. Thus, it is not surprising that many diseases are associated with RASGRP1 disorders. Here, we present an overview of the structure and function of RASGRP1, its crucial role in the development, expression, and regulation of immune cells, and its involvement in various signaling pathways. This review comprehensively explores the relationship between RASGRP1 and various diseases, elucidates the underlying molecular mechanisms of RASGRP1 in each disease, and identifies potential therapeutic targets. This study provides novel insights into the role of RASGRP1 in insulin secretion and highlights its potential as a therapeutic target for diabetes. The limitations and challenges associated with studying RASGRP1 in disease are also discussed.

Convergent and divergent genes expression profiles associated with brain-wide functional connectome dysfunction in deficit and non-deficit schizophrenia

Deficit schizophrenia (DS) is a subtype of schizophrenia characterized by the primary and persistent negative symptoms. Previous studies have identified differences in brain functions between DS and non-deficit schizophrenia (NDS) patients. However, the genetic regulation features underlying these abnormal changes are still unknown. This study aimed to detect the altered patterns of functional connectivity (FC) in DS and NDS and investigate the gene expression profiles underlying these abnormal FC. The study recruited 82 DS patients, 96 NDS patients, and 124 healthy controls (CN). Voxel-based unbiased brain-wide association study was performed to reveal altered patterns of FC in DS and NDS patients. Machine learning techniques were used to access the utility of altered FC for diseases diagnosis. Weighted gene co-expression network analysis (WGCNA) was employed to explore the associations between altered FC and gene expression of 6 donated brains. Enrichment analysis was conducted to identify the genetic profiles, and the spatio-temporal expression patterns of the key genes were further explored. Comparing to CN, 23 and 20 brain regions with altered FC were identified in DS and NDS patients. The altered FC among these regions showed significant correlations with the SDS scores and exhibited high efficiency in disease classification. WGCNA revealed associations between DS/NDS-related gene expression and altered FC. Additionally, 22 overlapped genes, including 12 positive regulation genes and 10 negative regulation genes, were found between NDS and DS. Enrichment analyses demonstrated relationships between identified genes and significant pathways related to cellular response, neuro regulation, receptor binding, and channel activity. Spatial and temporal gene expression profiles of SCN1B showed the lowest expression at the initiation of embryonic development, while DPYSL3 exhibited rapid increased in the fetal. The present study revealed different altered patterns of FC in DS and NDS patients and highlighted the potential value of FC in disease classification. The associations between gene expression and neuroimaging provided insights into specific and common genetic regulation underlying these brain functional changes in DS and NDS, suggesting a potential genetic-imaging pathogenesis of schizophrenia.

Integration of machine learning to identify diagnostic genes in leukocytes for acute myocardial infarction patients

BACKGROUND

Acute myocardial infarction (AMI) has two clinical characteristics: high missed diagnosis and dysfunction of leukocytes. Transcriptional RNA on leukocytes is closely related to the course evolution of AMI patients. We hypothesized that transcriptional RNA in leukocytes might provide potential diagnostic value for AMI. Integration machine learning (IML) was first used to explore AMI discrimination genes. The following clinical study was performed to validate the results.

METHODS

A total of four AMI microarrays (derived from the Gene Expression Omnibus) were included in bioanalysis (220 sample size). Then, the clinical validation was finished with 20 AMI and 20 stable coronary artery disease patients (SCAD). At a ratio of 5:2, GSE59867 was included in the training set, while GSE60993, GSE62646, and GSE48060 were included in the testing set. IML was explicitly proposed in this research, which is composed of six machine learning algorithms, including support vector machine (SVM), neural network (NN), random forest (RF), gradient boosting machine (GBM), decision trees (DT), and least absolute shrinkage and selection operator (LASSO). IML had two functions in this research: filtered optimized variables and predicted the categorized value. Finally, The RNA of the recruited patients was analyzed to verify the results of IML.

RESULTS

Thirty-nine differentially expressed genes (DEGs) were identified between controls and AMI individuals from the training sets. Among the thirty-nine DEGs, IML was used to process the predicted classification model and identify potential candidate genes with overall normalized weights > 1. Finally, two genes (AQP9 and SOCS3) show their diagnosis value with the area under the curve (AUC) > 0.9 in both the training and testing sets. The clinical study verified the significance of AQP9 and SOCS3. Notably, more stenotic coronary arteries or severe Killip classification indicated higher levels of these two genes, especially SOCS3. These two genes correlated with two immune cell types, monocytes and neutrophils.

CONCLUSION

AQP9 and SOCS3 in leukocytes may be conducive to identifying AMI patients with SCAD patients. AQP9 and SOCS3 are closely associated with monocytes and neutrophils, which might contribute to advancing AMI diagnosis and shed light on novel genetic markers. Multiple clinical characteristics, multicenter, and large-sample relevant trials are still needed to confirm its clinical value.

Investigating trait variability of gene co-expression network architecture in brain by controlling for genomic risk of schizophrenia

The effect of schizophrenia (SCZ) genetic risk on gene expression in brain remains elusive. A popular approach to this problem has been the application of gene co-expression network algorithms (e.g., WGCNA). To improve reliability with this method it is critical to remove unwanted sources of variance while also preserving biological signals of interest. In this WCGNA study of RNA-Seq data from postmortem prefrontal cortex (78 neurotypical donors, EUR ancestry), we tested the effects of SCZ genetic risk on co-expression networks. Specifically, we implemented a novel design in which gene expression was adjusted by linear regression models to preserve or remove variance explained by biological signal of interest (GWAS genomic scores for SCZ risk-(GS-SCZ), and genomic scores- GS of height (GS-Ht) as a negative control), while removing variance explained by covariates of non-interest. We calculated co-expression networks from adjusted expression (GS-SCZ and GS-Ht preserved or removed), and consensus between them (representative of a "background" network free of genomic scores effects). We then tested the overlap between GS-SCZ preserved modules and background networks reasoning that modules with reduced overlap would be most affected by GS-SCZ biology. Additionally, we tested these modules for convergence of SCZ risk (i.e., enrichment in PGC3 SCZ GWAS priority genes, enrichment in SCZ risk heritability and relevant biological ontologies. Our results highlight key aspects of GS-SCZ effects on brain co-expression networks, specifically: 1) preserving/removing SCZ genetic risk alters the co-expression modules; 2) biological pathways enriched in modules affected by GS-SCZ implicate processes of transcription, translation and metabolism that converge to influence synaptic transmission; 3) priority PGC3 SCZ GWAS genes and SCZ risk heritability are enriched in modules associated with GS-SCZ effects. Overall, our results indicate that gene co-expression networks that selectively integrate information about genetic risk can reveal novel combinations of biological pathways involved in schizophrenia.

Dopamine and schizophrenia from bench to bedside: Discovery of a striatal co-expression risk gene set that predicts in vivo measures of striatal function

Schizophrenia (SCZ) is characterized by a polygenic risk architecture implicating diverse molecular pathways important for synaptic function. However, how polygenic risk funnels through these pathways to translate into syndromic illness is unanswered. To evaluate biologically meaningful pathways of risk, we used tensor decomposition to characterize gene co-expression in post-mortem brain (of neurotypicals: N=154; patients with SCZ: N=84; and GTEX samples N=120) from caudate nucleus (CN), hippocampus (HP), and dorsolateral prefrontal cortex (DLPFC). We identified a CN-predominant gene set showing dopaminergic selectivity that was enriched for genes associated with clinical state and for genes associated with SCZ risk. Parsing polygenic risk score for SCZ based on this specific gene set (parsed-PRS), we found that greater pathway-specific SCZ risk predicted greater in vivo striatal dopamine synthesis capacity measured by [ 18 F]-FDOPA PET in three independent cohorts of neurotypicals and patients (total N=235) and greater fMRI striatal activation during reward anticipation in two additional independent neurotypical cohorts (total N=141). These results reveal a 'bench to bedside' translation of dopamine-linked genetic risk variation in driving in vivo striatal neurochemical and hemodynamic phenotypes that have long been implicated in the pathophysiology of SCZ.

A maturational shift in the frontal cortex synaptic transcriptional landscape underlies schizophrenia-relevant behavioural traits: A congenital rat model

Disruption of brain development early in life may underlie the neurobiology behind schizophrenia. We have reported more immature synaptic spines in the frontal cortex (FC) of adult Roman High-Avoidance (RHA-I) rats, a behavioural model displaying schizophrenia-like traits. Here, we performed a whole transcriptome analysis in the FC of 4 months old male RHA-I (n=8) and its counterpart, the Roman Low-Avoidance (RLA-I) (n=8). We identified 203 significant genes with overrepresentation of genes involved in synaptic function. Next, we performed a gene set enrichment analysis (GSEA) for genes co-expressed during neurodevelopment. Gene networks were obtained by weighted gene co-expression network analysis (WGCNA) of a transcriptomic dataset containing human FC during lifespan (n=269). Out of thirty-one functional gene networks, six were significantly enriched in the RHA-I. These were differentially regulated during infancy and enriched in biological ontologies related to myelination, synaptic function, and immune response. We validated differential gene expression in a new cohort of adolescent (<=2 months old) and young-adult (>=3 months old) RHA-I and RLA-I rats. The results confirmed overexpression of Gsn, Nt5cd1, Ppp1r1b, and Slc9a3r1 in young-adult RHA-I, while Cables1, a regulator of Cdk5 phosphorylation in actin regulation and involved in synaptic plasticity and maturation, was significantly downregulated in adolescent RHA-I. This age-related expression change was also observed for presynaptic components Snap25 and Snap29. Our results show a different maturational expression profile of synaptic components in the RHA-I strain, supporting a shift in FC maturation underlying schizophrenia-like behavioural traits and adding construct validity to this strain as a neurodevelopmental model.

The influence of tissue pH and RNA integrity number on gene expression of human postmortem brain

Background

Evaluating and controlling confounders are necessary when investigating molecular pathogenesis using human postmortem brain tissue. Particularly, tissue pH and RNA integrity number (RIN) are valuable indicators for controlling confounders. However, the influences of these indicators on the expression of each gene in postmortem brain have not been fully investigated. Therefore, we aimed to assess these effects on gene expressions of human brain samples.

Methods

We isolated total RNA from occipital lobes of 13 patients with schizophrenia and measured the RIN and tissue pH. Gene expression was analyzed and gene sets affected by tissue pH and RIN were identified. Moreover, we examined the functions of these genes by enrichment analysis and upstream regulator analysis.

Results

We identified 2,043 genes (24.7%) whose expressions were highly correlated with pH; 3,004 genes (36.3%) whose expressions were highly correlated with RIN; and 1,293 genes (15.6%) whose expressions were highly correlated with both pH and RIN. Genes commonly affected by tissue pH and RIN were highly associated with energy production and the immune system. In addition, genes uniquely affected by tissue pH were highly associated with the cell cycle, whereas those uniquely affected by RIN were highly associated with RNA processing.

Conclusion

The current study elucidated the influence of pH and RIN on gene expression profiling and identified gene sets whose expressions were affected by tissue pH or RIN. These findings would be helpful in the control of confounders for future postmortem brain studies.

Identification of schizophrenia symptom-related gene modules by postmortem brain transcriptome analysis

Schizophrenia is a multifactorial disorder, the genetic architecture of which remains unclear. Although many studies have examined the etiology of schizophrenia, the gene sets that contribute to its symptoms have not been fully investigated. In this study, we aimed to identify each gene set associated with corresponding symptoms of schizophrenia using the postmortem brains of 26 patients with schizophrenia and 51 controls. We classified genes expressed in the prefrontal cortex (analyzed by RNA-seq) into several modules by weighted gene co-expression network analysis (WGCNA) and examined the correlation between module expression and clinical characteristics. In addition, we calculated the polygenic risk score (PRS) for schizophrenia from Japanese genome-wide association studies, and investigated the association between the identified gene modules and PRS to evaluate whether genetic background affected gene expression. Finally, we conducted pathway analysis and upstream analysis using Ingenuity Pathway Analysis to clarify the functions and upstream regulators of symptom-related gene modules. As a result, three gene modules generated by WGCNA were significantly correlated with clinical characteristics, and one of these showed a significant association with PRS. Genes belonging to the transcriptional module associated with PRS significantly overlapped with signaling pathways of multiple sclerosis, neuroinflammation, and opioid use, suggesting that these pathways may also be profoundly implicated in schizophrenia. Upstream analysis indicated that genes in the detected module were profoundly regulated by lipopolysaccharides and CREB. This study identified schizophrenia symptom-related gene sets and their upstream regulators, revealing aspects of the pathophysiology of schizophrenia and identifying potential therapeutic targets.

Consensus molecular environment of schizophrenia risk genes in coexpression networks shifting across age and brain regions

Schizophrenia is a neurodevelopmental brain disorder whose genetic risk is associated with shifting clinical phenomena across the life span. We investigated the convergence of putative schizophrenia risk genes in brain coexpression networks in postmortem human prefrontal cortex (DLPFC), hippocampus, caudate nucleus, and dentate gyrus granule cells, parsed by specific age periods (total N = 833). The results support an early prefrontal involvement in the biology underlying schizophrenia and reveal a dynamic interplay of regions in which age parsing explains more variance in schizophrenia risk compared to lumping all age periods together. Across multiple data sources and publications, we identify 28 genes that are the most consistently found partners in modules enriched for schizophrenia risk genes in DLPFC; twenty-three are previously unidentified associations with schizophrenia. In iPSC-derived neurons, the relationship of these genes with schizophrenia risk genes is maintained. The genetic architecture of schizophrenia is embedded in shifting coexpression patterns across brain regions and time, potentially underwriting its shifting clinical presentation.

Identification and validation of key genes associated with atrial fibrillation in the elderly

Background

Atrial fibrillation (AF) is the most common cardiac arrhythmia and significantly increases the risk of stroke and heart failure (HF), contributing to a higher mortality rate. Increasing age is a major risk factor for AF; however, the mechanisms of how aging contributes to the occurrence and progression of AF remain unclear. This study conducted weighted gene co-expression network analysis (WGCNA) to identify key modules and hub genes and determine their potential associations with aging-related AF.

Materials and methods

WGCNA was performed using the AF dataset GSE2240 obtained from the Gene Expression Omnibus, which contained data from atrial myocardium in cardiac patients with permanent AF or sinus rhythm (SR). Hub genes were identified in clinical samples. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were also performed.

Results

Green and pink were the most critical modules associated with AF, from which nine hub genes, PTGDS, COLQ, ASTN2, VASH1, RCAN1, AMIGO2, RBP1, MFAP4, and ALDH1A1, were hypothesized to play key roles in the AF pathophysiology in elderly and seven of them have high diagnostic value. Functional enrichment analysis demonstrated that the green module was associated with the calcium, cyclic adenosine monophosphate (cAMP), and peroxisome proliferator-activated receptors (PPAR) signaling pathways, and the pink module may be associated with the transforming growth factor beta (TGF-β) signaling pathway in myocardial fibrosis.

Conclusion

We identified nine genes that may play crucial roles in the pathophysiological mechanism of aging-related AF, among which six genes were associated with AF for the first time. This study provided novel insights into the impact of aging on the occurrence and progression of AF, and identified biomarkers and potential therapeutic targets for AF.

The foundation and architecture of precision medicine in neurology and psychiatry

Neurological and psychiatric diseases have high degrees of genetic and pathophysiological heterogeneity, irrespective of clinical manifestations. Traditional medical paradigms have focused on late-stage syndromic aspects of these diseases, with little consideration of the underlying biology. Advances in disease modeling and methodological design have paved the way for the development of precision medicine (PM), an established concept in oncology with growing attention from other medical specialties. We propose a PM architecture for central nervous system diseases built on four converging pillars: multimodal biomarkers, systems medicine, digital health technologies, and data science. We discuss Alzheimer's disease (AD), an area of significant unmet medical need, as a case-in-point for the proposed framework. AD can be seen as one of the most advanced PM-oriented disease models and as a compelling catalyzer towards PM-oriented neuroscience drug development and advanced healthcare practice.

Bioinformatics analysis of ferroptosis-related gene AKR1C3 as a potential biomarker of asthma and its identification in BEAS-2B cells

Ferroptosis is a newly discovered type of cell death and has recently been shown to be associated with asthma. However, the relationship between them at the genetic level has not been elucidated via informatics analysis. In this study, bioinformatics analyses are conducted using asthma and ferroptosis datasets to identify candidate ferroptosis-related genes using the R software. Weighted gene co-expression network analysis is performed to identify co-expressed genes. Protein-protein interaction networks, the Kyoto encyclopedia of genes and genomes, and gene ontology enrichment analysis are used to identify the potential functions of the candidate genes. We experimentally validate the results of our analysis using small interfering RNAs and plasmids to silence and upregulate the expression of the candidate gene in human bronchial epithelial cells (BEAS-2B). The ferroptosis signature levels are examined. Bioinformatics analysis of the asthma dataset GDS4896 shows that the level of the aldo-keto reductase family 1 member C3 (AKR1C3) gene in the peripheral blood of patients with severe therapy-resistant asthma and controlled persistent mild asthma (MA) is significantly upregulated. The AUC values for asthma diagnosis and MA are 0.823 and 0.915, respectively. The diagnostic value of AKR1C3 is verified using the GSE64913 dataset. The gene module of AKR1C3 is evident in MA and functions through redox reactions and metabolic processes. Ferroptosis indicators are downregulated by the overexpression of AKR1C3 and upregulated by silencing AKR1C3. The ferroptosis-related gene AKR1C3 can be used as a diagnostic biomarker for asthma, particularly for MA, and regulates ferroptosis in BEAS-2B cells.

Polygenic scores for psychiatric disorders in a diverse postmortem brain tissue cohort

A new era of human postmortem tissue research has emerged thanks to the development of 'omics technologies that measure genes, proteins, and spatial parameters in unprecedented detail. Also newly possible is the ability to construct polygenic scores, individual-level metrics of genetic risk (also known as polygenic risk scores/PRS), based on genome-wide association studies, GWAS. Here, we report on clinical, educational, and brain gene expression correlates of polygenic scores in ancestrally diverse samples from the Human Brain Collection Core (HBCC). Genotypes from 1418 donors were subjected to quality control filters, imputed, and used to construct polygenic scores. Polygenic scores for schizophrenia predicted schizophrenia status in donors of European ancestry (p = 4.7 × 10^-8, 17.2%) and in donors with African ancestry (p = 1.6 × 10^-5, 10.4% of phenotypic variance explained). This pattern of higher variance explained among European ancestry samples was also observed for other psychiatric disorders (depression, bipolar disorder, substance use disorders, anxiety disorders) and for height, body mass index, and years of education. For a subset of 223 samples, gene expression from dorsolateral prefrontal cortex (DLPFC) was available through the CommonMind Consortium. In this subgroup, schizophrenia polygenic scores also predicted an aggregate gene expression score for schizophrenia (European ancestry: p = 0.0032, African ancestry: p = 0.15). Overall, polygenic scores performed as expected in ancestrally diverse samples, given historical biases toward use of European ancestry samples and variable predictive power of polygenic scores across phenotypes. The transcriptomic results reported here suggest that inherited schizophrenia genetic risk influences gene expression, even in adulthood. For future research, these and additional polygenic scores are being made available for analyses, and for selecting samples, using postmortem tissue from the Human Brain Collection Core.

Multiple-model machine learning identifies potential functional genes in dilated cardiomyopathy

Introduction

Machine learning (ML) has gained intensive popularity in various fields, such as disease diagnosis in healthcare. However, it has limitation for single algorithm to explore the diagnosing value of dilated cardiomyopathy (DCM). We aim to develop a novel overall normalized sum weight of multiple-model MLs to assess the diagnosing value in DCM.

Methods

Gene expression data were selected from previously published databases (six sets of eligible microarrays, 386 samples) with eligible criteria. Two sets of microarrays were used as training; the others were studied in the testing sets (ratio 5:1). Totally, we identified 20 differently expressed genes (DEGs) between DCM and control individuals (7 upregulated and 13 down-regulated).

Results

We developed six classification ML methods to identify potential candidate genes based on their overall weights. Three genes, serine proteinase inhibitor A3 (SERPINA3), frizzled-related proteins (FRPs) 3 (FRZB), and ficolin 3 (FCN3) were finally identified as the receiver operating characteristic (ROC). Interestingly, we found all three genes correlated considerably with plasma cells. Importantly, not only in training sets but also testing sets, the areas under the curve (AUCs) for SERPINA3, FRZB, and FCN3 were greater than 0.88. The ROC of SERPINA3 was significantly high (0.940 in training and 0.918 in testing sets), indicating it is a potentially functional gene in DCM. Especially, the plasma levels in DCM patients of SERPINA3, FCN, and FRZB were significant compared with healthy control.

Discussion

SERPINA3, FRZB, and FCN3 might be potential diagnosis targets for DCM, Further verification work could be implemented.

Gene co-expression architecture in peripheral blood in a cohort of remitted first-episode schizophrenia patients

A better understanding of schizophrenia subtypes is necessary to stratify the patients according to clinical attributes. To explore the genomic architecture of schizophrenia symptomatology, we analyzed blood co-expression modules and their association with clinical data from patients in remission after a first episode of schizophrenia. In total, 91 participants of the 2EPS project were included. Gene expression was assessed using the Clariom S Human Array. Weighted-gene co-expression network analysis (WGCNA) was applied to identify modules of co-expressed genes and to test its correlation with global functioning, clinical symptomatology, and premorbid adjustment. Among the 25 modules identified, six modules were significantly correlated with clinical data. These modules could be clustered in two groups according to their correlation with clinical data. Hub genes in each group showing overlap with risk genes for schizophrenia were enriched in biological processes related to metabolic processes, regulation of gene expression, cellular localization and protein transport, immune processes, and neurotrophin pathways. Our results indicate that modules with significant associations with clinical data showed overlap with gene sets previously identified in differential gene-expression analysis in brain, indicating that peripheral tissues could reveal pathogenic mechanisms. Hub genes involved in these modules revealed multiple signaling pathways previously related to schizophrenia, which may represent the complex interplay in the pathological mechanisms behind the disease. These genes could represent potential targets for the development of peripheral biomarkers underlying illness traits in clinical remission stages after a first episode of schizophrenia.

Lessons Learned From Parsing Genetic Risk for Schizophrenia Into Biological Pathways

The clinically heterogeneous presentation of schizophrenia is compounded by the heterogeneity of risk factors and neurobiological correlates of the disorder. Genome-wide association studies in schizophrenia have uncovered a remarkably high number of genetic variants, but the biological pathways they impact upon remain largely unidentified. Among the diverse methodological approaches employed to provide a more granular understanding of genetic risk for schizophrenia, the use of biological labels, such as gene ontologies, regulome approaches, and gene coexpression have all provided novel perspectives into how genetic risk translates into the neurobiology of schizophrenia. Here, we review the salient aspects of parsing polygenic risk for schizophrenia into biological pathways. We argue that parsed scores, compared to standard polygenic risk scores, may afford a more biologically plausible and accurate physiological modeling of the different dimensions involved in translating genetic risk into brain mechanisms, including multiple brain regions, cell types, and maturation stages. We discuss caveats, opportunities, and pitfalls inherent in the parsed risk approach.

Comprehensive analysis of omics data identifies relevant gene networks for Attention-Deficit/Hyperactivity Disorder (ADHD)

Attention-deficit/hyperactivity disorder (ADHD) is a highly prevalent neurodevelopmental disorder that results from the interaction of both genetic and environmental risk factors. Genome-wide association studies have started to identify multiple genetic risk loci associated with ADHD, however, the exact causal genes and biological mechanisms remain largely unknown. We performed a multi-step analysis to identify and characterize modules of co-expressed genes associated with ADHD using data from peripheral blood mononuclear cells of 270 ADHD cases and 279 controls. We identified seven ADHD-associated modules of co-expressed genes, some of them enriched in both genetic and epigenetic signatures for ADHD and in biological pathways relevant for psychiatric disorders, such as the regulation of gene expression, epigenetics and immune system. In addition, for some of the modules, we found evidence of potential regulatory mechanisms, including microRNAs and common genetic variants. In conclusion, our results point to promising genes and pathways for ADHD, supporting the use of peripheral blood to assess gene expression signatures in psychiatric disorders. Furthermore, they highlight that the combination of multi-omics signals provides deeper and broader insights into the biological mechanisms underlying ADHD.

Identification of immune-related key genes in the peripheral blood of ischaemic stroke patients using a weighted gene coexpression network analysis and machine learning

Background

The immune system plays a vital role in the pathological process of ischaemic stroke. However, the exact immune-related mechanism remains unclear. The current research aimed to identify immune-related key genes associated with ischaemic stroke.

Methods

CIBERSORT was utilized to reveal the immune cell infiltration pattern in ischaemic stroke patients. Meanwhile, a weighted gene coexpression network analysis (WGCNA) was utilized to identify meaningful modules significantly correlated with ischaemic stroke. The characteristic genes correlated with ischaemic stroke were identified by the following two machine learning methods: the support vector machine-recursive feature elimination (SVM-RFE) algorithm and least absolute shrinkage and selection operator (LASSO) logistic regression.

Results

The CIBERSORT results suggested that there was a decreased infiltration of naive CD4 T cells, CD8 T cells, resting mast cells and eosinophils and an increased infiltration of neutrophils, M0 macrophages and activated memory CD4 T cells in ischaemic stroke patients. Then, three significant modules (pink, brown and cyan) were identified to be significantly associated with ischaemic stroke. The gene enrichment analysis indicated that 519 genes in the above three modules were mainly involved in several inflammatory or immune-related signalling pathways and biological processes. Eight hub genes (ADM, ANXA3, CARD6, CPQ, SLC22A4, UBE2S, VIM and ZFP36) were revealed to be significantly correlated with ischaemic stroke by the LASSO logistic regression and SVM-RFE algorithm. The external validation combined with a RT‒qPCR analysis revealed that the expression levels of ADM, ANXA3, SLC22A4 and VIM were significantly increased in ischaemic stroke patients and that these key genes were positively associated with neutrophils and M0 macrophages and negatively correlated with CD8 T cells. The mean AUC value of ADM, ANXA3, SLC22A4 and VIM was 0.80, 0.87, 0.91 and 0.88 in the training set, 0.85, 0.77, 0.86 and 0.72 in the testing set and 0.87, 0.83, 0.88 and 0.91 in the validation samples, respectively.

Conclusions

These results suggest that the ADM, ANXA3, SLC22A4 and VIM genes are reliable serum markers for the diagnosis of ischaemic stroke and that immune cell infiltration plays a crucial role in the occurrence and development of ischaemic stroke.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12967-022-03562-w.

Identification of potential biomarkers and their correlation with immune infiltration cells in schizophrenia using combinative bioinformatics strategy

Many studies have identified changes in gene expression in brains of schizophrenia patients and their altered molecular processes, but the findings in different datasets were inconsistent and diverse. Here we performed the most comprehensive analysis of gene expression patterns to explore the underlying mechanisms and the potential biomarkers for early diagnosis in schizophrenia. We focused on 10 gene expression datasets in post-mortem human brain samples of schizophrenia downloaded from gene expression omnibus (GEO) database using the integrated bioinformatics analyses including robust rank aggregation (RRA) algorithm, Weighted gene co-expression network analysis (WGCNA) and CIBERSORT. Machine learning algorithm was used to construct the risk prediction model for early diagnosis of schizophrenia. We identified 15 key genes (SLC1A3, AQP4, GJA1, ALDH1L1, SOX9, SLC4A4, EGR1, NOTCH2, PVALB, ID4, ABCG2, METTL7A, ARC, F3 and EMX2) in schizophrenia by performing multiple bioinformatics analysis algorithms. Moreover, the interesting part of the study is that there is a correlation between the expression of hub genes and the immune infiltrating cells estimated by CIBERSORT. Besides, the risk prediction model was constructed by using both these genes and the immune cells with a high accuracy of 0.83 in the training set, and achieved a high AUC of 0.77 for the test set. Our study identified several potential biomarkers for diagnosis of SCZ based on multiple bioinformatics algorithms, and the constructed risk prediction model using these biomarkers achieved high accuracy. The results provide evidence for an improved understanding of the molecular mechanism of schizophrenia.

Integrative identification of immune-related key genes in atrial fibrillation using weighted gene coexpression network analysis and machine learning

Background

The immune system significantly participates in the pathologic process of atrial fibrillation (AF). However, the molecular mechanisms underlying this participation are not completely explained. The current research aimed to identify critical genes and immune cells that participate in the pathologic process of AF.

Methods

CIBERSORT was utilized to reveal the immune cell infiltration pattern in AF patients. Meanwhile, weighted gene coexpression network analysis (WGCNA) was utilized to identify meaningful modules that were significantly correlated with AF. The characteristic genes correlated with AF were identified by the least absolute shrinkage and selection operator (LASSO) logistic regression and support vector machine recursive feature elimination (SVM-RFE) algorithm.

Results

In comparison to sinus rhythm (SR) individuals, we observed that fewer activated mast cells and regulatory T cells (Tregs), as well as more gamma delta T cells, resting mast cells, and M2 macrophages, were infiltrated in AF patients. Three significant modules (pink, red, and magenta) were identified to be significantly associated with AF. Gene enrichment analysis showed that all 717 genes were associated with immunity- or inflammation-related pathways and biological processes. Four hub genes (GALNT16, HTR2B, BEX2, and RAB8A) were revealed to be significantly correlated with AF by the SVM-RFE algorithm and LASSO logistic regression. qRT–PCR results suggested that compared to the SR subjects, AF patients exhibited significantly reduced BEX2 and GALNT16 expression, as well as dramatically elevated HTR2B expression. The AUC measurement showed that the diagnostic efficiency of BEX2, HTR2B, and GALNT16 in the training set was 0.836, 0.883, and 0.893, respectively, and 0.858, 0.861, and 0.915, respectively, in the validation set.

Conclusions

Three novel genes, BEX2, HTR2B, and GALNT16, were identified by WGCNA combined with machine learning, which provides potential new therapeutic targets for the early diagnosis and prevention of AF.

A polygenic score indexing a DRD2-related co-expression network is associated with striatal dopamine function

The D2 dopamine receptor (D2R) is the primary site of the therapeutic action of antipsychotics and is involved in essential brain functions relevant to schizophrenia, such as attention, memory, motivation, and emotion processing. Moreover, the gene coding for D2R (DRD2) has been associated with schizophrenia at a genome-wide level. Recent studies have shown that a polygenic co-expression index (PCI) predicting the brain-specific expression of a network of genes co-expressed with DRD2 was associated with response to antipsychotics, brain function during working memory in patients with schizophrenia, and with the modulation of prefrontal cortex activity after pharmacological stimulation of D2 receptors. We aimed to investigate the relationship between the DRD2 gene network and in vivo striatal dopaminergic function, which is a phenotype robustly associated with psychosis and schizophrenia. To this aim, a sample of 92 healthy subjects underwent 18F-DOPA PET and was genotyped for genetic variations indexing the co-expression of the DRD2-related genetic network in order to calculate the PCI for each subject. The PCI was significantly associated with whole striatal dopamine synthesis capacity (p = 0.038). Exploratory analyses on the striatal subdivisions revealed a numerically larger effect size of the PCI on dopamine function for the associative striatum, although this was not significantly different than effects in other sub-divisions. These results are in line with a possible relationship between the DRD2-related co-expression network and schizophrenia and extend it by identifying a potential mechanism involving the regulation of dopamine synthesis. Future studies are needed to clarify the molecular mechanisms implicated in this relationship.

Identifying patterns of immune related cells and genes in the peripheral blood of acute myocardial infarction patients using a small cohort

Background

The immune system plays a vital role in the pathophysiology of acute myocardial infarction (AMI). However, the exact immune related mechanism is still unclear. This research study aimed to identify key immune-related genes involved in AMI.

Methods

CIBERSORT, a deconvolution algorithm, was used to determine the proportions of 22 subsets of immune cells in blood samples. The weighted gene co-expression network analysis (WGCNA) was used to identify key modules that are significantly associated with AMI. Then, CIBERSORT combined with WGCNA were used to identify key immune-modules. The protein–protein interaction (PPI) network was constructed and Molecular Complex Detection (MCODE) combined with cytoHubba plugins were used to identify key immune-related genes that may play an important role in the occurrence and progression of AMI.

Results

The CIBERSORT results suggested that there was a decrease in the infiltration of CD8 + T cells, gamma delta (γδ) T cells, and resting mast cells, along with an increase in the infiltration of neutrophils and M0 macrophages in AMI patients. Then, two modules (midnightblue and lightyellow) that were significantly correlated with AMI were identified, and the salmon module was found to be significantly associated with memory B cells. Gene enrichment analysis indicated that the 1,171 genes included in the salmon module are mainly involved in immune-related biological processes. MCODE analysis was used to identify four different MCODE complexes in the salmon module, while four hub genes (EEF1B2, RAC2, SPI1, and ITGAM) were found to be significantly correlated with AMI. The correlation analysis between the key genes and infiltrating immune cells showed that SPI1 and ITGAM were positively associated with neutrophils and M0 macrophages, while they were negatively associated with CD8 + T cells, γδ T cells, regulatory T cells (Tregs), and resting mast cells. The RT-qPCR validation results found that the expression of the ITGAM and SPI1 genes were significantly elevated in the AMI samples compared with the samples from healthy individuals, and the ROC curve analysis showed that ITGAM and SPI1 had a high diagnostic efficiency for the recognition of AMI.

Conclusions

Immune cell infiltration plays a crucial role in the occurrence and development of AMI. ITGAM and SPI1 are key immune-related genes that are potential novel targets for the prevention and treatment of AMI.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12967-022-03517-1.

Transcriptional Profiling Reveals Brain Region-Specific Gene Networks Regulated in Exercise in a Mouse Model of Parkinson’s Disease

Background

Exercise plays an essential role in improving motor symptoms in Parkinson’s disease (PD), but the underlying mechanism in the central nervous system remains unclear.

Methods

Motor ability was observed after 12-week treadmill exercise on a 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP)-induced mouse model of PD. RNA-sequencing on four brain regions (cerebellum, cortex, substantia nigra (SN), and striatum) from control animals, MPTP-induced PD, and MPTP-induced PD model treated with exercise for 12 weeks were performed. Transcriptional networks on the four regions were further identified by an integrative network biology approach.

Results

The 12-week treadmill exercise significantly improved the motor ability of an MPTP-induced mouse model of PD. RNA-seq analysis showed SN and striatum were remarkably different among individual region’s response to exercise in the PD model. Especially, synaptic regulation pathways about axon guidance, synapse assembly, neurogenesis, synaptogenesis, transmitter transport-related pathway, and synaptic regulation genes, including Neurod2, Rtn4rl2, and Cd5, were upregulated in SN and striatum. Lastly, immunofluorescence staining revealed that exercise rescued the loss of TH+ synapses in the striatal region in PD mice, which validates the key role of synaptic regulation pathways in exercise-induced protective effects in vivo.

Conclusion

SN and striatum are important brain regions in which critical transcriptional changes, such as in synaptic regulation pathways, occur after the exercise intervention on the PD model.

Metabolic polygenic risk scores effect on antipsychotic-induced metabolic dysregulation: A longitudinal study in a first episode psychosis cohort

OBJECTIVE

Metabolic syndrome is a health-threatening condition suffered by approximately one third of schizophrenia patients and largely attributed to antipsychotic medication. Previous evidence reports a common genetic background of psychotic and metabolic disorders. In this study, we aimed to assess the role of polygenic risk scores (PRSs) on the progression of the metabolic profile in a first-episode psychosis (FEP) cohort.

METHOD

Of the 231 FEP individuals included in the study, 192-220 participants were included in basal analysis and 118-179 in longitudinal 6-month models. Eleven psychopathologic and metabolic PRSs were constructed. Basal and longitudinal PRSs association with metabolic measurements was assessed by statistical analyses.

RESULTS

No major association of psychopathological PRSs with the metabolic progression was found. However, high risk individuals for depression and cholesterol-related PRSs reported a higher increase of cholesterol levels during the follow-up (FDR ≤ 0.023 for all analyses). Their effect was comparable to other well-established pharmacological and environmental risk factors (explaining at least 1.2% of total variance).

CONCLUSION

Our findings provide new evidence of the effects of metabolic genetic risk on the development of metabolic dysregulation. The future establishment of genetic profiling tools in clinical procedures could enable practitioners to better personalize antipsychotic treatment selection and dosage.

An opportunity for primary prevention research in psychotic disorders

An opportunity has opened for research into primary prevention of psychotic disorders, based on progress in endophenotypes, genetics, and genomics. Primary prevention requires reliable prediction of susceptibility before any symptoms are present. We studied a battery of measures where published data supports abnormalities of these measurements prior to appearance of initial psychosis symptoms. These neurobiological and behavioral measurements included cognition, eye movement tracking, Event Related Potentials, and polygenic risk scores. They generated an acceptably precise separation of healthy controls from outpatients with a psychotic disorder. METHODS: The Bipolar and Schizophrenia Network on Intermediate Phenotypes (B-SNIP) measured this battery in an ancestry-diverse series of consecutively recruited adult outpatients with a psychotic disorder and healthy controls. Participants include all genders, 16 to 50 years of age, 261 with psychotic disorders (Schizophrenia (SZ) 109, Bipolar with psychosis (BPP) 92, Schizoaffective disorder (SAD) 60), 110 healthy controls. Logistic Regression, and an extension of the Linear Mixed Model to include analysis of pairwise interactions between measures (Environmental kernel Relationship Matrices (ERM)) with multiple iterations, were performed to predict case-control status. Each regression analysis was validated with four-fold cross-validation. RESULTS AND CONCLUSIONS: Sensitivity, specificity, and Area Under the Curve of Receiver Operating Characteristic of 85%, 62%, and 86%, respectively, were obtained for both analytic methods. These prediction metrics demonstrate a promising diagnostic distinction based on premorbid risk variables. There were also statistically significant pairwise interactions between measures in the ERM model. The strong prediction metrics of both types of analytic model provide proof-of-principle for biologically-based laboratory tests as a first step toward primary prevention studies. Prospective studies of adolescents at elevated risk, vs. healthy adolescent controls, would be a next step toward development of primary prevention strategies.

Genetics and Brain Transcriptomics of Completed Suicide

OBJECTIVE

The authors sought to study the transcriptomic and genomic features of completed suicide by parsing the method chosen, to capture molecular correlates of the distinctive frame of mind of individuals who die by suicide, while reducing heterogeneity.

METHODS

The authors analyzed gene expression (RNA sequencing) from postmortem dorsolateral prefrontal cortex of patients who died by suicide with violent compared with nonviolent means, nonsuicide patients with the same psychiatric disorders, and a neurotypical group (total N=329). They then examined genomic risk scores (GRSs) for each psychiatric disorder included, and GRSs for cognition (IQ) and for suicide attempt, testing how they predict diagnosis or traits (total N=888).

RESULTS

Patients who died by suicide by violent means showed a transcriptomic pattern remarkably divergent from each of the other patient groups but less from the neurotypical group; consistently, their genomic profile of risk was relatively low for their diagnosed illness as well as for suicide attempt, and relatively high for IQ: they were more similar to the neurotypical group than to other patients. Differentially expressed genes (DEGs) associated with patients who died by suicide by violent means pointed to purinergic signaling in microglia, showing similarities to a genome-wide association study of Drosophila aggression. Weighted gene coexpression network analysis revealed that these DEGs were coexpressed in a context of mitochondrial metabolic activation unique to suicide by violent means.

CONCLUSIONS

These findings suggest that patients who die by suicide by violent means are in part biologically separable from other patients with the same diagnoses, and their behavioral outcome may be less dependent on genetic risk for conventional psychiatric disorders and be associated with an alteration of purinergic signaling and mitochondrial metabolism.

Abnormal RasGRP1 Expression in the Post-Mortem Brain and Blood Serum of Schizophrenia Patients

Schizophrenia (SCZ) is a polygenic severe mental illness. Genome-wide association studies (GWAS) have detected genomic variants associated with this psychiatric disorder and pathway analyses have indicated immune system and dopamine signaling as core components of risk in dorsolateral-prefrontal cortex (DLPFC) and hippocampus, but the mechanistic links remain unknown. The RasGRP1 gene, encoding for a guanine nucleotide exchange factor, is implicated in dopamine signaling and immune response. RasGRP1 has been identified as a candidate risk gene for SCZ and autoimmune disease, therefore representing a possible point of convergence between mechanisms involving the nervous and the immune system. Here, we investigated RasGRP1 mRNA and protein expression in post-mortem DLPFC and hippocampus of SCZ patients and healthy controls, along with RasGRP1 protein content in the serum of an independent cohort of SCZ patients and control subjects. Differences in RasGRP1 expression between SCZ patients and controls were detected both in DLPFC and peripheral blood of samples analyzed. Our results indicate RasGRP1 may mediate risk for SCZ by involving DLPFC and peripheral blood, thus encouraging further studies to explore its possible role as a biomarker of the disease and/or a target for new medication.

Electrophysiological measures from human iPSC-derived neurons are associated with schizophrenia clinical status and predict individual cognitive performance

Neurons derived from human induced pluripotent stem cells (hiPSCs) have been used to model basic cellular aspects of neuropsychiatric disorders, but the relationship between the emergent phenotypes and the clinical characteristics of donor individuals has been unclear. We analyzed RNA expression and indices of cellular function in hiPSC-derived neural progenitors and cortical neurons generated from 13 individuals with high polygenic risk scores (PRSs) for schizophrenia (SCZ) and a clinical diagnosis of SCZ, along with 15 neurotypical individuals with low PRS. We identified electrophysiological measures in the patient-derived neurons that implicated altered Na+ channel function, action potential interspike interval, and gamma-aminobutyric acid-ergic neurotransmission. Importantly, electrophysiological measures predicted cardinal clinical and cognitive features found in these SCZ patients. The identification of basic neuronal physiological properties related to core clinical characteristics of illness is a potentially critical step in generating leads for novel therapeutics.

WGCNA-Based Identification of Hub Genes and Key Pathways Involved in Nonalcoholic Fatty Liver Disease

Background

The morbidity of nonalcoholic fatty liver disease (NAFLD) has been rising, but the pathogenesis of NAFLD is still elusive. This study is aimed at determining NAFLD-related hub genes based on weighted gene coexpression network analysis (WGCNA).

Methods

GSE126848 dataset based construction of coexpression networks was performed based on WGCNA. Database for Annotation, Visualization, and Integrated Discovery (DAVID) was utilized for Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. Hub genes were identified and validated in independent datasets and mouse model.

Results

We found that the steelblue module was most significantly correlated with NAFLD. Total 15 hub genes (NDUFA9, UQCRQ, NDUFB8, COPS5, RPS17, UBL5, PSMA3, PSMA1, SF3B5, MRPL27, RPL26, PDCD5, PFDN6, SNRPD2, PSMB3) were derived from both the coexpression and PPI networks and considered “true” hub genes. Functional enrichment analysis showed that the hub genes were related to NAFLD pathway and oxidative phosphorylation. Independent dataset-based analysis and the establishment of NAFLD mouse model confirmed the involvement of two hub genes NDUFA9 and UQCRQ in the pathogenesis of NAFLD.

Conclusions

Oxidative phosphorylation and NAFLD pathway may be crucially involved in the pathogenesis of NAFLD, and two hub genes NDUFA9 and UQCRQ might be diagnostic biomarkers and therapeutic targets for NAFLD.

MiR-574-5P, miR-1827, and miR-4429 as Potential Biomarkers for Schizophrenia

Schizophrenia is a severe chronic debilitating disorder with millions of affected individuals. Diagnosis is based on clinical presentations, which are made when the progressive disease has appeared. Early diagnosis may help improve the clinical outcomes and response to treatments. Lack of a reliable molecular diagnostic invokes the identification of novel biomarkers. To elucidate the molecular basis of the disease, in this study we used two mRNA expression arrays, including GSE93987 and GSE38485, and one miRNA array, GSE54914, and meta-analysis was conducted for evaluation of mRNA expression arrays via metaDE package. Using WGCNA package, we performed network analysis for both mRNA expression arrays separately. Then, we constructed protein-protein interaction network for significant modules. Limma package was employed to analyze the miRNA array for identification of dysregulated miRNAs (DEMs). Using genes of significant modules and DEMs, a mRNA-miRNA network was constructed and hub genes and miRNAs were identified. To confirm the dysregulated genes, expression values were evaluated through available datasets including GSE62333, GSE93987, and GSE38485. The ability of the detected hub miRNAs to discriminate schizophrenia from healthy controls was evaluated by assessing the receiver-operating curve. Finally, the expression levels of genes and miRNAs were evaluated in 40 schizophrenia patients compared with healthy controls via Real-Time PCR. The results confirmed dysregulation of hsa-miR-574-5P, hsa-miR-1827, hsa-miR-4429, CREBRF, ARPP19, TGFBR2, and YWHAZ in blood samples of schizophrenia patients. In conclusion, three miRNAs including hsa-miR-574-5P, hsa-miR-1827, and hsa-miR-4429 are suggested as potential biomarkers for diagnosis of schizophrenia.

Identification of the key genes and immune infiltrating cells determined by sex differences in ischaemic stroke through co-expression network module

Stroke is one of the leading causes of patients' death and long-term disability worldwide, and ischaemic stroke (IS) accounts for nearly 80% of all strokes. Differential genes and weighted gene co-expression network analysis (WGCNA) in male and female patients with IS were compared. The authors used cell type identification by estimating relative subsets of RNA transcripts (CIBERSORT) to analyse the distribution pattern of immune subtypes between male and female patients. In this study, 141 up-regulated and 61 down-regulated genes were gathered and distributed into five modules in response to their correlation degree to clinical traits. The criterion for Gene Ontology (GO) term and Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathway indicated that detailed analysis had the potential to enhance clinical prediction and to identify the gender-related mechanism. After that, the expression levels of hub genes were measured via the quantitative real-time PCR (qRT-PCR) method. Finally, CCL20, ICAM1 and PTGS2 were identified and these may be some promising targets for sex differences in IS. Besides, the hub genes were further verified by rat experiments. Furthermore, these CIBERSORT results showed that T cells CD8 and Monocytes may be the target for the treatment of male and female patients, respectively.

Identification of Biomarkers Related to Immune Cell Infiltration with Gene Coexpression Network in Myocardial Infarction

Background

There is evidence that the immune system plays a key critical role in the pathogenesis of myocardial infarction (MI). However, the exact mechanisms associated with immunity have not been systematically uncovered.

Methods

This study used the weighted gene coexpression network analysis (WGCNA) and the CIBERSORT algorithm to analyze the MI expression data from the Gene Expression Omnibus database and then identify the module associated with immune cell infiltration. In addition, we built the coexpression network and protein-protein interactions network analysis to identify the hub genes. Furthermore, the relationship between hub genes and NK cell resting was validated by using another dataset GSE123342. Finally, receiver operating characteristic (ROC) curve analyses were used to assess the diagnostic value of verified hub genes.

Results

Monocytes and neutrophils were markedly increased, and T cell CD8, T cell CD4 naive, T cell CD4 memory resting, and NK cell resting were significantly decreased in MI groups compared with stable coronary artery disease (CAD) groups. The WGCNA results showed that the pink model had the highest correlation with the NK cell resting infiltration level. We identified 11 hub genes whose expression correlated to the NK cell resting infiltration level, among which, 7 hub genes (NKG7, TBX21, PRF1, CD247, KLRD1, FASLG, and EOMES) were successfully validated in GSE123342. And these 7 genes had diagnostic value to distinguish MI and stable CAD.

Conclusions

NKG7, TBX21, PRF1, CD247, KLRD1, FASLG, and EOMES may be a diagnostic biomarker and therapeutic target associated with NK cell resting infiltration in MI.

Step-by-Step Construction of Gene Co-Expression Network Analysis for Identifying Novel Biomarkers of Sepsis Occurrence and Progression

Background

Sepsis is the leading cause of death in critically ill patients. Although it is well known that the immune system plays a key role in sepsis, exactly how it works remains unknown.

Methods

In our study, we used weighted gene co-expression network analysis (WGCNA) to screen out the immune-related genes that may play a critical role in the process of sepsis.

Results

A total of three sepsis-related hub genes were screened for further verification. Subsequent analysis of immune subtypes suggested their potential predictive effect in the clinic.

Conclusion

Our study shows that three immune-related genes CHMP1A, MED15 and MGAT1 are important biomarkers of sepsis. The screened genes may help to distinguish normal individuals from patients with different degrees of sepsis.

Network modules linking expression and methylation in prefrontal cortex of schizophrenia

Tremendous work has demonstrated the critical roles of genetics, epigenetics as well as their interplay in brain transcriptional regulations in the pathology of schizophrenia (SZ). There is great success currently in the dissection of the genetic components underlying risk-conferring transcriptomic networks. However, the study of regulating effect of epigenetics in the etiopathogenesis of SZ still faces many challenges. In this work, we investigated DNA methylation and gene expression from the dorsolateral prefrontal cortex (DLPFC) region of schizophrenia patients and healthy controls using weighted correlation network approach. We identified and replicated two expression and two methylation modules significantly associated with SZ. Among them, one pair of expression and methylation modules were significantly overlapped in the module genes which were significantly enriched in astrocyte-associated functional pathways, and specifically expressed in astrocytes. Another two linked expression-methylation module pairs were involved ageing process with module genes mostly related to oligodendrocyte development and myelination, and specifically expressed in oligodendrocytes. Further examination of underlying quantitative trait loci (QTLs) showed significant enrichment in genetic risk of most psychiatric disorders for expression QTLs but not for methylation QTLs. These results support the coherence between methylation and gene expression at the network level, and suggest a combinatorial effect of genetics and epigenetics in regulating gene expression networks specific to glia cells in relation to SZ and ageing process.

Identification of Prognostic Genes in Neuroblastoma in Children by Weighted Gene Coexpression Network Analysis

Background

Neuroblastoma is a malignant neuroendocrine tumor from the sympathetic nervous system, the most common extracranial tumor in children. Identifying potential prognostic markers of neuroblastoma can provide clues for early diagnosis, recurrence, and treatment.

Methods

RNA sequence data and clinical features of 147 neuroblastomas were obtained from the TARGET (Therapeutically Applicable Research to Generate Effective Treatments project) database. Application weighted gene coexpression network analysis (WGCNA) was used to construct a free-scale gene coexpression network, to study the interrelationship between its potential modules and clinical features, and to identify hub genes in the module. We performed Lasso regression and Cox regression analyses to identify the three most important genes and develop a new prognostic model. Data from the GSE85047 cohort verified the predictive accuracy of the prognostic model.

Results

14 coexpression modules were constructed using WGCNA. Brown coexpression modules were found to be significantly associated with disease survival status. Multivariate Cox analysis was performed on genes from univariate Cox regression and Lasso regression analyses using the Cox proportional hazards regression model. Finally, we constructed a three-gene prognostic model: risk score = (0.003812659∗CKB) + (−0.152376975∗expDST) + (0.032032815∗expDUT). The prognosis of samples in the high-risk group was significantly poorer than that of samples in the low-risk group (P=1.225e − 06). The risk model was also regarded as an independent predictor of prognosis (HR = 1.632; 95% CI = 1.391–1.934; P < 0.001).

Conclusion

Our study constructed a neuroblastoma coexpressing gene module and identified a prognostic potential risk model for prognosis in neuroblastoma.

Identification of Potential Prognostic Biomarkers Associated With Cancerometastasis in Skin Cutaneous Melanoma

Skin cutaneous melanoma (SKCM) is a highly aggressive tumor. The mortality and drug resistance among it are high. Thus, exploring predictive biomarkers for prognosis has become a priority. We aimed to find immune cell-based biomarkers for survival prediction. Here 321 genes were differentially expressed in immune-related groups after ESTIMATE analysis and differential analysis. Two hundred nineteen of them were associated with the metastasis of SKCM via weighted gene co-expression network analysis. Twenty-six genes in this module were hub genes. Twelve of the 26 genes were related to overall survival in SKCM patients. After a multivariable Cox regression analysis, we obtained six of these genes (PLA2G2D, IKZF3, MS4A1, ZC3H12D, FCRL3, and P2RY10) that were independent prognostic signatures, and a survival model of them performed excellent predictive efficacy. The results revealed several essential genes that may act as significant prognostic factors of SKCM, which could deepen our understanding of the metastatic mechanisms and improve cancer treatment.

Identification of Macrophage Polarization-Related Genes as Biomarkers of Chronic Obstructive Pulmonary Disease Based on Bioinformatics Analyses

OBJECTIVES

Chronic obstructive pulmonary disease (COPD) is characterized by lung inflammation and remodeling. Macrophage polarization is associated with inflammation and tissue remodeling, as well as immunity. Therefore, this study attempts to investigate the diagnostic value and regulatory mechanism of macrophage polarization-related genes for COPD by bioinformatics analysis and to provide a new theoretical basis for experimental research.

METHODS

The raw gene expression profile dataset (GSE124180) was collected from the Gene Expression Omnibus (GEO) database. Next, a weighted gene coexpression network analysis (WGCNA) was conducted to screen macrophage polarization-related genes. The differentially expressed genes (DEGs) between the COPD and normal samples were generated using DESeq2 v3.11 and overlapped with the macrophage polarization-related genes. Moreover, functional annotations of overlapped genes were conducted by Database for Annotation, Visualization and Integrated Discovery (DAVID) Bioinformatics Resource. The immune-related genes were selected, and their correlation with the differential immune cells was analyzed by Pearson. Finally, receiver operating characteristic (ROC) curves were used to verify the diagnostic value of genes.

RESULTS

A total of 4922 coexpressed genes related to macrophage polarization were overlapped with the 203 DEGs between the COPD and normal samples, obtaining 25 genes related to COPD and macrophage polarization. GEM, S100B, and GZMA of them participated in the immune response, which were considered the candidate biomarkers. GEM and S100B were significantly correlated with marker genes of B cells which had a significant difference between the COPD and normal samples. Moreover, GEM was highly associated with the genes in the PI3K/Akt/GSK3β signaling pathway, regulation of actin cytoskeleton, and calcium signaling pathway based on a Pearson correlation analysis of the candidate genes and the genes in the B cell receptor signaling pathway. PPI network analysis also indicated that GEM might participate in the regulation of the PI3K/Akt/GSK3β signaling pathway. The ROC curve showed that GEM possessed an excellent accuracy in distinguishing COPD from normal samples.

CONCLUSIONS

The data provide a transcriptome-based evidence that GEM is related to COPD and macrophage polarization likely contributes to COPD diagnosis. At the same time, it is hoped that in-depth functional mining can provide new ideas for exploring the COPD pathogenesis.

Convergent and distributed effects of the 3q29 deletion on the human neural transcriptome

The 3q29 deletion (3q29Del) confers high risk for schizophrenia and other neurodevelopmental and psychiatric disorders. However, no single gene in this interval is definitively associated with disease, prompting the hypothesis that neuropsychiatric sequelae emerge upon loss of multiple functionally-connected genes. 3q29 genes are unevenly annotated and the impact of 3q29Del on the human neural transcriptome is unknown. To systematically formulate unbiased hypotheses about molecular mechanisms linking 3q29Del to neuropsychiatric illness, we conducted a systems-level network analysis of the non-pathological adult human cortical transcriptome and generated evidence-based predictions that relate 3q29 genes to novel functions and disease associations. The 21 protein-coding genes located in the interval segregated into seven clusters of highly co-expressed genes, demonstrating both convergent and distributed effects of 3q29Del across the interrogated transcriptomic landscape. Pathway analysis of these clusters indicated involvement in nervous-system functions, including synaptic signaling and organization, as well as core cellular functions, including transcriptional regulation, posttranslational modifications, chromatin remodeling, and mitochondrial metabolism. Top network-neighbors of 3q29 genes showed significant overlap with known schizophrenia, autism, and intellectual disability-risk genes, suggesting that 3q29Del biology is relevant to idiopathic disease. Leveraging "guilt by association", we propose nine 3q29 genes, including one hub gene, as prioritized drivers of neuropsychiatric risk. These results provide testable hypotheses for experimental analysis on causal drivers and mechanisms of the largest known genetic risk factor for schizophrenia and highlight the study of normal function in non-pathological postmortem tissue to further our understanding of psychiatric genetics, especially for rare syndromes like 3q29Del, where access to neural tissue from carriers is unavailable or limited.

Tumor Purity Coexpressed Genes Related to Immune Microenvironment and Clinical Outcomes of Lung Adenocarcinoma

Purpose

Lung cancer tissue includes tumor tissue, stromal cells, immune cells, and epithelial cells. These nontumor cells dilute the tumor purity in lung cancer tissues. Tumor purity plays an essential role in the immune response to lung cancer. At present, the biological processes related to the purity of lung cancer tumors remains unclear.

Methods

We measured tumor purity in 486 lung carcinoma tissues from TCGA-LUAD FPKM by using the "estimate" R package. Lung carcinoma tumor mutation burden was calculated by analyzing TCGA single nucleotide polymorphism data. The immune cell proportion was also evacuated via the CIBERSORT method. Lung carcinoma samples with P < 0.05 were considered significant. Based on the tumor purity and lung carcinoma gene matrix, we performed weighted gene coexpression network analysis (WGCNA), and the tumor purity-related module was identified. Then, we analyzed the functions of the factors involved in the module. We screened the coexpressed factors related to clinical outcome and immunophenotype. Finally, expression levels of these factors were measured at tissue and single-cell levels.

Results

A lung cancer tumor purity correlated coexpression network was determined. Five coexpressed genes (CD4, CD53, EVI2B, PLEK, and SASH3) were identified as tumor purity coexpressed genes that negatively correlated with tumor purity. Because the factors in the coexpression network often participate in similar biological processes, we found that CD4, CD53, EVI2B, PLEK, and SASH3 were most related to positive regulation of cytokine production and interleukin-2 production through functional enrichment. In a clinical phenotype analysis, we found that these five factors can be used as independent prognostic risk factors. We found that these factors were significantly negatively correlated with tumor purity and positively correlated with the immune score in the immunophenotyping analysis. Using GSEA analysis, we found that the antigen processing and presentation pathway were related to the five tumor coexpressed genes mentioned above. SASH3 and CD53 were used to conduct a prognostic model based on the interaction analysis of the Support Vector Machine and the Least Absolute Shrinkage and Selection Operator. SASH3 was verified to be related to CD8A using a single-cell analysis.

Conclusion

Tumor purity-related coexpression factors in the tumor microenvironment have essential clinical, genomic, and biological significance in lung cancer. These coexpression factors (SASH3 and CD53) can be used to classify tumor purity phenotypes and to predict clinical outcomes.

Identification of Susceptible Genes for Chronic Obstructive Pulmonary Disease with Lung Adenocarcinoma by Weighted Gene Co-Expression Network Analysis

Purpose

Chronic obstructive pulmonary disease (COPD) and lung adenocarcinoma (LUAD) are common disorders and usually co-exists. However, genetic mechanisms between COPD and LUAD are rarely reported. This study aims to identify susceptible genes of COPD with LUAD.

Methods

Using the published data of GSE106899, co-expression modules were constructed by weighted gene co-expression network analysis (WGCNA). Subsequently, top 50 genes in the most tumor-related module were identified, among which hub genes were selected and validated.

Results

Twenty co-expression modules were constructed on 13,865 genes from 62 lung tissues of COPD patients with or without LUAD, in which one module (blue) was most related to tumorigenesis. Functional enrichment analyses showed that the genes in the blue module were mainly enriched in cell cycle, DNA transcription/replication and cancer pathways, etc. Combined with protein–protein interaction network, MTA1, PKMYT1 and FZR1 genes had the most intramodular connectivity, which were regarded as the hub genes. However, only FZR1 was validated to be overexpressed in lung tissues of COPD with LUAD and cigarette smoke extract-stimulated A549 cells, a human LUAD cell line.

Conclusion

This study suggests overexpression of FZR1 may play a key role in the tumorigenesis of LUAD in patients with COPD.