Comparative proteomic analysis across the developmental stages of the Eimeria tenella.

Eimeria tenella is the main pathogen responsible for coccidiosis in chickens. The life cycle of E. tenella is, arguably, the least complex of all Coccidia, with only one host. However, it presents different developmental stages, either in the environment or in the host and either intracellular or extracellular. Its signaling and metabolic pathways change with its different developmental stages. Until now, little is known about the developmental regulation and transformation mechanisms of its life cycle. In this study, protein profiles from the five developmental stages, including unsporulated oocysts (USO), partially sporulated (7 h) oocysts (SO7h), sporulated oocysts (SO), sporozoites (S) and second-generation merozoites (M2), were harvested using the label-free quantitative proteomics approach. Then the differentially expressed proteins (DEPs) for these stages were identified. A total of 314, 432, 689, and 665 DEPs were identified from the comparison of SO7h vs USO, SO vs SO7h, S vs SO, and M2 vs S, respectively. By conducting weighted gene coexpression network analysis (WGCNA), six modules were dissected. Proteins in blue and brown modules were calculated to be significantly positively correlated with the E. tenella developmental stages of sporozoites (S) and second-generation merozoites (M2), respectively. In addition, hub proteins with high intra-module degree were identified. Gene Ontology (GO) and Kyoto Encyclopedia of Gene and Genomes (KEGG) pathway enrichment analyses revealed that hub proteins in blue modules were involved in electron transport chain and oxidative phosphorylation. Hub proteins in the brown module were involved in RNA splicing. These findings provide new clues and ideas to enhance our fundamental understanding of the molecular mechanisms underlying parasite development.

Transcriptional landscape of cholangiocarcinoma revealed by weighted gene coexpression network analysis.

Cholangiocarcinoma (CCA) is a type of cancer with limited treatment options and a poor prognosis. Although some important genes and pathways associated with CCA have been identified, the relationship between coexpression and phenotype in CCA at the systems level remains unclear. In this study, the relationships underlying the molecular and clinical characteristics of CCA were investigated by employing weighted gene coexpression network analysis (WGCNA). The gene expression profiles and clinical features of 36 patients with CCA were analyzed to identify differentially expressed genes (DEGs). Subsequently, the coexpression of DEGs was determined by using the WGCNA method to investigate the correlations between pairs of genes. Network modules that were significantly correlated with clinical traits were identified. In total, 1478 mRNAs were found to be aberrantly expressed in CCA. Seven coexpression modules that significantly correlated with clinical characteristics were identified and assigned representative colors. Among the 7 modules, the green and blue modules were significantly related to tumor differentiation. Seventy-eight hub genes that were correlated with tumor differentiation were found in the green and blue modules. Survival analysis showed that 17 hub genes were prognostic biomarkers for CCA patients. In addition, we found five new targets (ISM1, SULT1B1, KIFC1, AURKB and CCNB1) that have not been studied in the context of CCA and verified their differential expression in CCA through experiments. Our results not only promote our understanding of the relationship between the transcriptome and clinical data in CCA but will also guide the development of targeted molecular therapy for CCA.

Identification of functional pathways and potential genes associated with interferon signaling during human adenovirus type 7 infection by weighted gene coexpression network analysis.

Human adenovirus type 7 (HAdV-7) can cause severe pneumonia and complications in children. However, the mechanism of pathogenesis and the genes involved remain largely unknown. We collected HAdV-7-infected and mock-infected A549 cells at 24, 48, and 72 hours postinfection (hpi) for RNA sequencing (RNA-Seq) and identified potential genes and functional pathways associated with HAdV-7 infection using weighted gene coexpression network analysis (WGCNA). Based on bioinformatics analysis, 12 coexpression modules were constructed by WGCNA, with the blue, tan, and brown modules significantly positively correlated with adenovirus infection at 24, 48, and 72 hpi, respectively. Functional enrichment analysis indicated that the blue module was mainly enriched in DNA replication and viral processes, the tan module was largely enriched in metabolic pathways and regulation of superoxide radical removal, and the brown module was predominantly enriched in regulation of cell death. qPCR was used to determine transcript abundance of some identified hub genes, and the results were consistent with those from RNA-Seq. Comprehensively analyzing hub genes and differentially expressed genes in the GSE68004 dataset, we identified SOCS3, OASL, ISG15, and IFIT1 as potential candidate genes for use as biomarkers or drug targets in HAdV-7 infection. We propose a multi-target inhibition of the interferon signaling mechanism to explain the association of HAdV-7 infection with the severity of clinical consequences. This study has allowed us to construct a framework of coexpression gene modules in A549 cells infected with HAdV-7, thus providing a basis for identifying potential genes and pathways involved in adenovirus infection and for investigating the pathogenesis of adenovirus-associated diseases.

An integrated bioinformatics approach to identify key biomarkers in the tubulointerstitium of patients with focal segmental glomerulosclerosis and construction of mRNA-miRNA-lncRNA/circRNA networks.

OBJECTIVE: The purpose of this study was to identify potential biomarkers in the tubulointerstitium of focal segmental glomerulosclerosis (FSGS) and comprehensively analyze its mRNA-miRNA-lncRNA/circRNA network. METHODS: The expression data (GSE108112 and GSE200818) were downloaded from the Gene Expression Omnibus database (https://www.ncbi.nlm.nih.gov/geo/). Identification and enrichment analysis of differentially expressed genes (DEGs) were performed. the PPI networks of the DEGs were constructed and classified using the Cytoscape molecular complex detection (MCODE) plugin. Weighted gene coexpression network analysis (WGCNA) was used to identify critical gene modules. Least absolute shrinkage and selection operator regression analysis were used to screen for key biomarkers of the tubulointerstitium in FSGS, and the receiver operating characteristic curve was used to determine their diagnostic accuracy. The screening results were verified by quantitative real-time-PCR (qRT-PCR) and Western blot. The transcription factors (TFs) affecting the hub genes were identified by Cytoscape iRegulon. The mRNA-miRNA-lncRNA/circRNA network for identifying potential biomarkers was based on the starBase database. RESULTS: A total of 535 DEGs were identified. MCODE obtained eight modules. The green module of WGCNA had the greatest association with the tubulointerstitium in FSGS. PPARG coactivator 1 alpha (PPARGC1A) was screened as a potential tubulointerstitial biomarker for FSGS and verified by qRT-PCR and Western blot. The TFs FOXO4 and FOXO1 had a regulatory effect on PPARGC1A. The ceRNA network yielded 17 miRNAs, 32 lncRNAs, and 50 circRNAs. CONCLUSIONS: PPARGC1A may be a potential biomarker in the tubulointerstitium of FSGS. The ceRNA network contributes to the comprehensive elucidation of the mechanisms of tubulointerstitial lesions in FSGS.

Identification of key biomarkers of the glomerulus in focal segmental glomerulosclerosis and their relationship with immune cell infiltration based on WGCNA and the LASSO algorithm.

OBJECTIVE: The aim of our study was to identify key biomarkers of glomeruli in focal glomerulosclerosis (FSGS) and analyze their relationship with the infiltration of immune cells. METHODS: The expression profiles (GSE108109 and GSE200828) were obtained from the GEO database. The differentially expressed genes (DEGs) were filtered and analyzed by gene set enrichment analysis (GSEA). MCODE module was constructed. Weighted gene coexpression network analysis (WGCNA) was performed to obtain the core gene modules. Least absolute shrinkage and selection operator (LASSO) regression was applied to identify key genes. ROC curves were employed to explore their diagnostic accuracy. Transcription factor prediction of the key biomarkers was performed using the Cytoscape plugin IRegulon. The analysis of the infiltration of 28 immune cells and their correlation with the key biomarkers were performed. RESULTS: A total of 1474 DEGs were identified. Their functions were mostly related to immune-related diseases and signaling pathways. MCODE identified five modules. The turquoise module of WGCNA had significant relevance to the glomerulus in FSGS. TGFB1 and NOTCH1 were identified as potential key glomerular biomarkers in FSGS. Eighteen transcription factors were obtained from the two hub genes. Immune infiltration showed significant correlations with T cells. The results of immune cell infiltration and their relationship with key biomarkers implied that NOTCH1 and TGFB1 were enhanced in immune-related pathways. CONCLUSION: TGFB1 and NOTCH1 may be strongly correlated with the pathogenesis of the glomerulus in FSGS and are new candidate key biomarkers. T-cell infiltration plays an essential role in the FSGS lesion process.

The Landscapes of Gluten Regulatory Network in Elite Wheat Cultivars Contrasting in Gluten Strength.

Yangmai-13 (YM13) is a wheat cultivar with weak gluten fractions. In contrast, Zhenmai-168 (ZM168) is an elite wheat cultivar known for its strong gluten fractions and has been widely used in a number of breeding programs. However, the genetic mechanisms underlying the gluten signatures of ZM168 remain largely unclear. To address this, we combined RNA-seq and PacBio full-length sequencing technology to unveil the potential mechanisms of ZM168 grain quality. A total of 44,709 transcripts were identified in Y13N (YM13 treated with nitrogen) and 51,942 transcripts in Z168N (ZM168 treated with nitrogen), including 28,016 and 28,626 novel isoforms in Y13N and Z168N, respectively. Five hundred and eighty-four differential alternative splicing (AS) events and 491 long noncoding RNAs (lncRNAs) were discovered. Incorporating the sodium-dodecyl-sulfate (SDS) sedimentation volume (SSV) trait, both weighted gene coexpression network analysis (WGCNA) and multiscale embedded gene coexpression network analysis (MEGENA) were employed for network construction and prediction of key drivers. Fifteen new candidates have emerged in association with SSV, including 4 transcription factors (TFs) and 11 transcripts that partake in the post-translational modification pathway. The transcriptome atlas provides new perspectives on wheat grain quality and would be beneficial for developing promising strategies for breeding programs.

Transcriptome-wide analyses of adipose tissue in outbred rats reveal genetic regulatory mechanisms relevant for human obesity.

Transcriptomic analysis in metabolically active tissues allows a systems genetics approach to identify causal genes and networks involved in metabolic disease. Outbred heterogeneous stock (HS) rats are used for genetic mapping of complex traits, but to-date, a systems genetics analysis of metabolic tissues has not been done. We investigated whether adiposity-associated genes and gene coexpression networks in outbred heterogeneous stock (HS) rats overlap those found in humans. We analyzed RNAseq data from adipose tissue of 415 male HS rats, correlated these transcripts with body weight (BW) and compared transcriptome signatures to two human cohorts: the "African American Genetics of Metabolism and Expression" and "Metabolic Syndrome in Men." We used weighted gene coexpression network analysis to identify adiposity-associated gene networks and mediation analysis to identify genes under genetic control whose expression drives adiposity. We identified 554 orthologous "consensus genes" whose expression correlates with BW in the rat and with body mass index (BMI) in both human cohorts. Consensus genes fell within eight coexpressed networks and were enriched for genes involved in immune system function, cell growth, extracellular matrix organization, and lipid metabolic processes. We identified 19 consensus genes for which genetic variation may influence BW via their expression, including those involved in lipolysis (e.g., Hcar1), inflammation (e.g., Rgs1), adipogenesis (e.g., Tmem120b), or no previously known role in obesity (e.g., St14 and Ms4a6a). Strong concordance between HS rat and human BW/BMI associated transcripts demonstrates translational utility of the rat model, while identification of novel genes expands our knowledge of the genetics underlying obesity.

Complexity untwined: The structure and function of cucumber (Cucumis sativus L.) shoot phloem.

Cucurbit phloem is complex, with large sieve tubes on both sides of the xylem (bicollateral phloem), and extrafascicular elements that form an intricate web linking the rest of the vasculature. Little is known of the physical interconnections between these networks or their functional specialization, largely because the extrafascicular phloem strands branch and turn at irregular angles. Here, export in the phloem from specific regions of the lamina of cucumber (Cucumis sativus L.) was mapped using carboxyfluorescein and 14 C as mobile tracers. We also mapped vascular architecture by conventional microscopy and X-ray computed tomography using optimized whole-tissue staining procedures. Differential gene expression in the internal (IP) and external phloem (EP) was analyzed by laser-capture microdissection followed by RNA-sequencing. The vascular bundles of the lamina form a nexus at the petiole junction, emerging in a predictable pattern, each bundle conducting photoassimilate from a specific region of the blade. The vascular bundles of the stem interconnect at the node, facilitating lateral transport around the stem. Elements of the extrafascicular phloem traverse the stem and petiole obliquely, joining the IP and EP of adjacent bundles. Using pairwise comparisons and weighted gene coexpression network analysis, we found differences in gene expression patterns between the petiole and stem and between IP and EP, and we identified hub genes of tissue-specific modules. Genes related to transport were expressed primarily in the EP while those involved in cell differentiation and development as well as amino acid transport and metabolism were expressed mainly in the IP.

Integrated weighted gene coexpression network analysis identifies Frizzled 2 (FZD2) as a key gene in invasive malignant pleomorphic adenoma.

BACKGROUND: Invasive malignant pleomorphic adenoma (IMPA) is a highly malignant neoplasm of the oral salivary glands with a poor prognosis and a considerable risk of recurrence. Many disease-causing genes of IMPA have been identified in recent decades (e.g., P53, PCNA and HMGA2), but many of these genes remain to be explored. Weighted gene coexpression network analysis (WGCNA) is a newly emerged algorithm that can cluster genes and form modules based on similar gene expression patterns. This study constructed a gene coexpression network of IMPA via WGCNA and then carried out multifaceted analysis to identify novel disease-causing genes. METHODS: RNA sequencing (RNA-seq) was performed for 10 pairs of IMPA and normal tissues to acquire the gene expression profiles. Differentially expressed genes (DEGs) were screened out with the cutoff criteria of |log2 Fold change (FC)|> 1 and adjusted p value < 0.05. Then, WGCNA was applied to systematically identify the hidden diagnostic hub genes of IMPA. RESULTS: In this research, a total of 1970 DEGs were screened out in IMPA tissues, including 1056 upregulated DEGs and 914 downregulated DEGs. Functional enrichment analysis was performed for identified DEGs and revealed an enrichment of tumor-associated GO terms and KEGG pathways. We used WGCNA to identify gene module most relevant with the histological grade of IMPA. The gene FZD2 was then recognized as the hub gene of the selected module with the highest module membership (MM) value and intramodule connectivity in protein-protein interaction (PPI) network. According to immunohistochemistry (IHC) staining, the expression level of FZD2 was higher in low-grade IMPA than in high-grade IMPA. CONCLUSION: FZD2 shows an expression dynamic that is negatively correlated with the clinical malignancy of IMPA and it plays a central role in the transcription network of IMPA. Thus, FZD2 serves as a promising histological indicator for the precise prediction of IMPA histological stages.

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.

The potential role of chemotaxis and the complement system in the formation and progression of thoracic aortic aneurysms inferred from the weighted gene coexpression network analysis.

BACKGROUND: Thoracic aortic aneurysm (TAA) can be life-threatening due to the progressive weakening and dilatation of the aortic wall. Once the aortic wall has ruptured, no effective pharmaceutical therapies are available. However, studies on TAA at the gene expression level are limited. Our study aimed to identify the driver genes and critical pathways of TAA through gene coexpression networks. METHODS: We analyzed the genetic data of TAA patients from a public database by weighted gene coexpression network analysis (WGCNA). Modules with clinical significance were identified, and the differentially expressed genes (DEGs) were intersected with the genes in these modules. Gene Ontology and pathway enrichment analyses were performed. Finally, hub genes that might be driving factors of TAA were identified. Furthermore, we evaluated the diagnostic accuracy of these genes and analyzed the composition of immune cells using the CIBERSORT algorithm. RESULTS: We identified 256 DEGs and two modules with clinical significance. The immune response, including leukocyte adhesion, mononuclear cell proliferation and T cell activation, was identified by functional enrichment analysis. CX3CR1, C3, and C3AR1 were the top 3 hub genes in the module correlated with TAA, and the areas under the curve (AUCs) by receiver operating characteristic (ROC) analysis of all the hub genes exceeded 0.7. Finally, we found that the proportions of infiltrating immune cells in TAA and normal tissues were different, especially in terms of macrophages and natural killer (NK) cells. CONCLUSION: Chemotaxis and the complement system were identified as crucial pathways in TAA, and macrophages with interactive immune cells may regulate this pathological process.

ESR1 as a recurrence-related gene in intrahepatic cholangiocarcinoma: a weighted gene coexpression network analysis.

BACKGROUND: Intrahepatic cholangiocarcinoma (iCCA) is the second most common malignant hepatic tumor and has a high postoperative recurrence rate and a poor prognosis. The key roles of most tumor recurrence-associated molecules in iCCA remain unclear. This study aimed to explore hub genes related to the postsurgical recurrence of iCCA. METHOD: Differentially expressed genes (DEGs) between iCCA samples and normal liver samples were screened from The Cancer Genome Atlas (TCGA) database and used to construct a weighted gene coexpression network. Module-trait correlations were calculated to identify the key module related to recurrence in iCCA patients. Genes in the key module were subjected to functional enrichment analysis, and candidate hub genes were filtered through coexpression and protein-protein interaction (PPI) network analysis. Validation studies were conducted to detect the "real" hub gene. Furthermore, the biological functions and the underlying mechanism of the real hub gene in iCCA tumorigenesis and progression were determined via in vitro experiments. RESULTS: A total of 1019 DEGs were filtered and used to construct four coexpression modules. The red module, which showed the highest correlations with the recurrence status, family history, and day to death of patients, was identified as the key module. Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses demonstrated that genes in the red module were enriched in genes and pathways related to tumorigenesis and tumor progression. We performed validation studies and identified estrogen receptor 1 (ESR1), which significantly impacted the prognosis of iCCA patients, as the real hub gene related to the recurrence of iCCA. The in vitro experiments demonstrated that ESR1 overexpression significantly suppressed cell proliferation, migration, and invasion, whereas ESR1 knockdown elicited opposite effects. Further investigation into the mechanism demonstrated that ESR1 acts as a tumor suppressor by inhibiting the JAK/STAT3 signaling pathway. CONCLUSIONS: ESR1 was identified as the real hub gene related to the recurrence of iCCA that plays a critical tumor suppressor role in iCCA progression. ESR1 significantly impacts the prognosis of iCCA patients and markedly suppresses cholangiocarcinoma cell proliferation, migration and invasion by inhibiting JAK/STAT3 signaling pathway.

Low expression of CHRDL1 and SPARCL1 predicts poor prognosis of lung adenocarcinoma based on comprehensive analysis and immunohistochemical validation.

PURPOSE: Exploring the molecular mechanisms of lung adenocarcinoma (LUAD) is beneficial for developing new therapeutic strategies and predicting prognosis. This study was performed to select core genes related to LUAD and to analyze their prognostic value. METHODS: Microarray datasets from the GEO (GSE75037) and TCGA-LUAD datasets were analyzed to identify differentially coexpressed genes in LUAD using weighted gene coexpression network analysis (WGCNA) and differential gene expression analysis. Functional enrichment analysis was conducted, and a protein-protein interaction (PPI) network was established. Subsequently, hub genes were identified using the CytoHubba plug-in. Overall survival (OS) analyses of hub genes were performed. The Clinical Proteomic Tumor Analysis Consortium (CPTAC) and the Human Protein Atlas (THPA) databases were used to validate our findings. Gene set enrichment analysis (GSEA) of survival-related hub genes were conducted. Immunohistochemistry (IHC) was carried out to validate our findings. RESULTS: We identified 486 differentially coexpressed genes. Functional enrichment analysis suggested these genes were primarily enriched in the regulation of epithelial cell proliferation, collagen-containing extracellular matrix, transforming growth factor beta binding, and signaling pathways regulating the pluripotency of stem cells. Ten hub genes were detected using the maximal clique centrality (MCC) algorithm, and four genes were closely associated with OS. The CPTAC and THPA databases revealed that CHRDL1 and SPARCL1 were downregulated at the mRNA and protein expression levels in LUAD, whereas SPP1 was upregulated. GSEA demonstrated that DNA-dependent DNA replication and catalytic activity acting on RNA were correlated with CHRDL1 and SPARCL1 expression, respectively. The IHC results suggested that CHRDL1 and SPARCL1 were significantly downregulated in LUAD. CONCLUSIONS: Our study revealed that survival-related hub genes closely correlated with the initiation and progression of LUAD. Furthermore, CHRDL1 and SPARCL1 are potential therapeutic and prognostic indicators of LUAD.

Screening of characteristic genes in ulcerative colitis by integrating gene expression profiles.

BACKGROUND: This study aimed to screen the feature modules and characteristic genes related to ulcerative colitis (UC) and construct a support vector machine (SVM) classifier to distinguish UC patients. METHODS: Four datasets that contained UC and control samples were obtained from the Gene Expression Omnibus database. Differentially expressed genes (DEGs) with consistency were screened via the MetaDE method. The weighted gene coexpression network (WGCNA) was used to distinguish significant modules based on the four datasets. The protein-protein interaction network was established based on intersection genes. Enrichment analysis of Gene Ontology (GO) biological processes (BPs) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment were established based on DAVID. An SVM combined with recursive feature elimination was also applied to construct a disease classifier for the disease diagnosis of UC patients. The efficacy of the SVM classifier was evaluated through receiver operating characteristic curves. RESULTS: Twelve highly preserved modules were obtained using the WGCNA, and 2009 DEGs with significant consistency were selected using the MetaDE method. Sixteen significantly related GO BPs and 12 KEGG pathways were obtained, such as cytokine-cytokine receptor interaction, cell adhesion molecules, and leukocyte transendothelial migration. Subsequently, 41 genes were used to construct an SVM classifier, such as CXCL1, CCR2, IL1B, and IL1A. The area under the curve (AUC) was 0.999 in the training dataset, whereas the AUC was 0.886, 0.790, and 0.819 in the validation set (GSE65114, GSE37283, and GSE36807, respectively). CONCLUSIONS: An SVM classifier based on feature genes might correctly identify healthy people or UC patients.

Identification of an immune-based mRNA-lncRNA signature for overall survival in cervical squamous cell carcinoma.

Aim: To better predict the survival of cervical squamous cell carcinoma (CESC) patients, we aimed to construct a signature according to different immune infiltration. Methods: We downloaded the RNA sequences of CESC patients from the Cancer Genome Atlas database. By using single-sample gene set enrichment analysis, we separated the samples into high- and low-immunity groups. Then we separated the samples into training and testing datasets and performed the following analyses: univariate, least absolute shrinkage and selection operator analysis, multivariate Cox regression analyses and weighted gene coexpression network analysis using R software. Gene ontology and Kyoto Encyclopedia of Genes and Genomes studies were performed using the Database for Annotation, Visualization and Integrated Discovery website. Results & conclusion: We finally identified a signature with three mRNAs and two lncRNAs: ADGRG5, HSH2D, ZMAT4, RBAKDN and LINC00200. In short, our study constructed an mRNA-lncRNA signature related to immune infiltration to better predict the survival of CESC patients.

Lay abstract Cervical squamous cell carcinoma is a prevalent cancer type among females. Our study was to construct a signature which can predict the overall survival time in cervical squamous cell carcinoma patients. We performed some analysis of genetic expression patterns using a public database of genetic data, and successfully constructed the model which holds a good prediction value.

Transcriptome analyses and weighted gene coexpression network analysis reveal key pathways and genes involved in the rapid cold resistance of the Chinese white wax scale insect.

The Chinese white wax scale insect, Ericerus pela, is an important resource insect in China. The rapid response of E. pela to decreasing temperatures plays key roles in the population distribution. In this study, we analyzed the gene expression of E. pela treated with low temperature using transcriptome analyses and weighted gene coexpression network analysis (WGCNA). The results showed that the cold resistance of E. pela involved changes in the expression of many genes. The genes were mainly involved in alcohol formation activity, lipid metabolism, membrane and structure, and oxidoreductase activity. According to the WGCNA results, some pathways related to cold resistance were found in the genes in the modules, such as cytoskeleton proteins, cytoskeleton protein pathway, biosynthesis of unsaturated fatty acids, glycerophospholipid metabolism, ether lipid metabolism, and thermogenesis. Some of the hub genes were nonspecific lipid-transfer proteins, DnaJ homolog subfamily C member 13, paramyosin, tropomodulin, and tubulin beta chain. In particular, the hub genes of the tan module included the heat shock protein (hsp) 10, hsp 60, hsp 70, and hsp 90 genes. Thirty-five antifreeze protein (afp) genes were identified according to the annotation results. Three afp genes were further identified among the hub genes. Six of these genes were selected for heterogeneous protein expression. One of them was expressed successfully. The thermal hysteresis activity (THA) analyses showed that the THA was 1.73°C. These results showed that the cytoskeleton, lipid metabolism, thermogenesis, HSPs and AFPs may play important roles in the cold resistance of E. pela.

Identification of crucial genes involved in pathogenesis of regional weakening of the aortic wall.

BACKGROUND: The diameter of the abdominal aortic aneurysm (AAA) is the most commonly used parameter for the prediction of occurrence of AAA rupture. However, the most vulnerable region of the aortic wall may be different from the most dilated region of AAA under pressure. The present study is the first to use weighted gene coexpression network analysis (WGCNA) to detect the coexpressed genes that result in regional weakening of the aortic wall. METHODS: The GSE165470 raw microarray dataset was used in the present study. Differentially expressed genes (DEGs) were filtered using the "limma" R package. DEGs were assessed by Gene Ontology biological process (GO-BP) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses. WGCNA was used to construct the coexpression networks in the samples with regional weakening of the AAA wall and in the control group to detect the gene modules. The hub genes were defined in the significant functional modules, and a hub differentially expressed gene (hDEG) coexpression network was constructed with the highest confidence based on protein-protein interactions (PPIs). Molecular compound detection (MCODE) was used to identify crucial genes in the hDEG coexpression network. Crucial genes in the hDEG coexpression network were validated using the GSE7084 and GSE57691 microarray gene expression datasets. RESULT: A total of 350 DEGs were identified, including 62 upregulated and 288 downregulated DEGs. The pathways were involved in immune responses, vascular smooth muscle contraction and cell-matrix adhesion of DEGs in the samples with regional weakening in AAA. Antiquewhite3 was the most significant module and was used to identify downregulated hDEGs based on the result of the most significant modules negatively related to the trait of weakened aneurysm walls. Seven crucial genes were identified and validated: ACTG2, CALD1, LMOD1, MYH11, MYL9, MYLK, and TPM2. These crucial genes were associated with the mechanisms of AAA progression. CONCLUSION: We identified crucial genes that may play a significant role in weakening of the AAA wall and may be potential targets for medical therapies and diagnostic biomarkers. Further studies are required to more comprehensively elucidate the functions of crucial genes in the pathogenesis of regional weakening in AAA.

Network-based analysis reveals novel gene signatures in the peripheral blood of patients with sporadic nonsyndromic thoracic aortic aneurysm.

Thoracic aortic aneurysm (TAA), a serious cardiovascular disease that causes morbidity and mortality worldwide. At present, few biomarkers can accurately diagnose the appearance of TAA before dissection or rupture. Our research has the intention to investigate the developing applicable biomarkers for TAA promising clinically diagnostic biomarkers or probable regulatory targets for TAA. In our research, we built correlation networks utilizing the expression profile of peripheral blood mononuclear cell obtained from a public microarray data set (GSE9106). Furthermore, we chose the turquoise module, which has the strongest significance with TAA and was further analyzed. Fourteen genes that overlapped with differentially expressed proteins in the medial aortic layer were obtained. Subsequently, we verified the results applying quantitative polymerase chain reaction (Q-PCR) to our clinical specimen. In general, the Q-PCR results coincide with the majority of the expression profile. Fascinatingly, a notable change occurred in CLU, DES, MYH10, and FBLN5. In summary, using weighted gene coexpression analysis, our study indicates that CLU, DES, MYH10, and FBLN5 were identified and validated to be related to TAA and might be candidate biomarkers or therapeutic targets for TAA.

Identification of potential proteases for abdominal aortic aneurysm by weighted gene coexpression network analysis.

Proteases are involved in the degradation of the extracellular matrix (ECM), which contributes to the formation of abdominal aortic aneurysm (AAA). To identify new disease targets in addition to the results of previous microarray studies, we performed next-generation sequencing (NGS) of the whole transcriptome of Angiotensin II-treated ApoE-/- male mice (n = 4) and control mice (n = 4) to obtain differentially expressed genes (DEGs). Identified DEGs of proteases were analyzed using weighted gene coexpression network analysis (WGCNA). RT-qPCR was conducted to validate the differential expression of selected hub genes. We found that 43 DEGs were correlated with the expression of the protease profile, and most were clustered in immune response module. Among 26 hub genes, we found that Mmp16 and Mmp17 were significantly downregulated in AAA mice, while Ctsa, Ctsc, and Ctsw were upregulated. Our functional annotation analysis of genes coexpressed with the five hub genes indicated that Ctsw and Mmp17 were involved in T cell regulation and Cell adhesion molecule pathway, respectively, and that both were involved in general regulation of the cell cycle and gene expression. Overall, our data suggest that these ectopic genes are potentially crucial to AAA formation and may act as biomarkers for the diagnosis of AAA.

A Selective P2Y Purinergic Receptor Agonist 2-MesADP Enhances Locomotor Recovery after Acute Spinal Cord Injury.

INTRODUCTION: Spinal cord injury (SCI) causes most severe motor and sensory dysfunctions. In Chinese traditional medicine, the agonist of a purinergic receptor is believed to have a positive effect on SCIs, and 2-Methylthio-adenosine-5'-diphosphate (2-MesADP) is a selective agonist of the P2Y purinergic receptor. METHODS: To investigate its therapeutic function and molecular mechanism in SCI, transcriptome analysis associated with weighted gene co-expression network analysis (WGCNA) was carried out at various time points after T9 crush injury. RESULTS: 2-MesADP demonstrated recovery of limb motor function at the 6 weeks after injury, accompanied by neuronal regeneration and axon remyelination at 2 and 6 weeks. Furthermore, gene profiling revealed alternated gene expression with the treatment of 2-MesADP. These genes were assigned to a total of 38 modules, followed by gene ontology analysis; of these, 18 represented neuronal apoptosis and regeneration, immune response, synaptic transmission, cell cycle, and angiogenesis. In the neuronal apoptosis and regeneration module, Nefh, NeuroD6, and Dcx in the 2-MesADP group were noticed due to their interesting expression pattern. The gene expression patterns of Mag, Mog, and Cnp, which played key roles in myelination, were significantly changed with the treatment of 2-MesADP. Wnt signal pathway was the most important pathway in 2-MesADP treatment for acute SCI. CONCLUSION: 2-MesADP enhanced locomotor recovery in mouse SCI by altering the expression of neuronal apoptosis and remyelination-related genes and Wnt signaling pathways.