Amino acid metabolism characteristics of Banxia Baizhu Tianma Decoction in realizing drug withdrawal based on transcriptomic analysis / 中国中药杂志

This study aimed to demonstrate the effect of Banxia Baizhu Tianma Decoction(BBTD) on realizing withdrawal of anti-epileptic drugs and explore the relationship between BBTD and the amino acid metabolism by transcriptomic analysis in the rat model of epilepsy induced by lithium chloride-pilocarpine. The rats with epilepsy were divided into a control group(Ctrl), an epilepsy group(Ep), a BBTD & antiepileptic drug integrative group(BADIG), and an antiepileptic drug withdrawal group(ADWG). The Ctrl and Ep were given ultrapure water by gavage for 12 weeks. The BADIG was given BBTD extract and carbamazepine solution by gavage for 12 weeks. The ADWG was given carbamazepine solution and BBTD extract by gavage for the former 6 weeks, and then only given BBTD extract for the latter 6 weeks. The therapeutic effect was evaluated by behavioral observation, electroencephalogram(EEG), and hippocampal neuronal morphological changes. High-throughput sequencing was used to obtain amino acid metabolism-related differen-tial genes in the hippocampus, and the mRNA expression in the hippocampus of each group was verified by real-time quantitative polymerase chain reaction(RT-qPCR). The hub genes were screened out through protein-protein interaction(PPI) network, and Gene Ontology(GO) functional enrichment analysis and Kyoto Encyclopedia of Genes and Genomes(KEGG) pathway enrichment analysis were performed. Two ceRNA networks, namely circRNA-miRNA-mRNA and lncRNA-miRNA-mRNA, were constructed for ADWG vs BADIG. The experimental results showed that compared with those in Ep, rats in ADWG were significantly improved in the behavioral observation, EEG, and hippocampal neuronal impairment. Thirty-four amino acid metabolism-related differential genes were obtained by transcriptomic analysis, and the sequencing results were confirmed by RT-qPCR. Eight hub genes were obtained through PPI network, involving several biological processes, molecular functions, and signal pathways related to amino acid metabolism. Finally, the circRNA-miRNA-mRNA ternary transcription network of 17 circRNA, 5 miRNA, and 2 mRNA, and a lncRNA-miRNA-mRNA ternary network of 10 lncRNA, 5 miRNA, and 2 mRNA were constructed in ADWG vs BADIG. In conclusion, BBTD can effectively achieve the withdrawal of antiepileptic drugs, which may be related to the transcriptomic regulation of amino acid metabolism.

Analysis of Significant Genes and Pathways in Esophageal Cancer Based on Gene Expression Omnibus Database / 中国医学科学杂志(英文版)

Objective To screen antigen targets for immunotherapy by analyzing over-expressed genes, and to identify significant pathways and molecular mechanisms in esophageal cancer by using bioinformatic methods such as enrichment analysis, protein-protein interaction (PPI) network, and survival analysis based on the Gene Expression Omnibus (GEO) database.Methods By screening with highly expressed genes, we mainly analyzed proteins MUC13 and EPCAM with transmembrane domain and antigen epitope from TMHMM and IEDB websites. Significant genes and pathways associated with the pathogenesis of esophageal cancer were identified using enrichment analysis, PPI network, and survival analysis. Several software and platforms including Prism 8, R language, Cytoscape, DAVID, STRING, and GEPIA platform were used in the search and/or figure creation.Results Genes MUC13 and EPCAM were over-expressed with several antigen epitopes in esophageal squamous cell carcinoma (ESCC) tissue. Enrichment analysis revealed that the process of keratinization was focused and a series of genes were related with the development of esophageal cancer. Four genes including ALDH3A1, C2, SLC6A1,and ZBTB7C were screened with significant P value of survival curve.Conclusions Genes MUC13 and EPCAM may be promising antigen targets or biomarkers for esophageal cancer. Keratinization may greatly impact the pathogenesis of esophageal cancer. Genes ALDH3A1, C2, SLC6A1,and ZBTB7C may play important roles in the development of esophageal cancer.

Cross - species regulation and underlying mechanisms of parasite - derived non-coding RNAs: a review / 中国血吸虫病防治杂志

Parasite-derived non-coding RNAs (ncRNAs) not only contribute to life activities of parasites, and microRNA (miRNA), long non-coding RNA (lncRNA), and circular RNA (circRNA) may generate a competitive endogenous RNA (ceRNA) regulatory network with host miRNAs and mRNAs via extracellular vesicles, thereby participating in infection and pathogenic processes. This article presents an overview of characterizing ncRNAs derived from parasites and the cross-species regulatory role of parasite-derived ncRNAs in host gene expression and its underlying mechanisms.

Screening of Genes Co-Associated with Sudden Infant Death Syndrome and Infectious Sudden Death in Infancy and Bioinformatics Analysis of Their Regulatory Networks / 法医学杂志

OBJECTIVES@#The common differentially expressed mRNAs in brain, heart and liver tissues of deceased sudden infant death syndrome (SIDS) and infectious sudden death in infancy (ISDI) confirmed by autopsy was screened by bioinformatics to explore the common molecular markers and pathogenesis of SIDS and ISDI.@*METHODS@#The datasets of GSE70422 and GSE136992 were downloaded, the limma of R software was used to screen differentially expressed mRNA in different tissue samples of SIDS and ISDI decedents for overlapping analysis. The clusterProfiler of R software was used to conduct gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis. The protein-protein interaction (PPI) network was constructed by STRING database, while the hub gene was screened by cytoHubba plug-in.@*RESULTS@#Compared with the control group, there were 19 significant differentially expressed genes in the tissue samples of SIDS and ISDI decedents, among which 16 in the heart tissue and 3 in the liver tissue, and the astrotactin 1 (ASTN1) gene expression difference in the heart tissue was most significant. The PPI network identified Ras homolog family member A (RHOA), integrin subunit alpha 1 (ITGA1), and H2B clustered histone 5 (H2BC5) were hub genes. The analysis of GO and KEGG showed that differentially expressed genes were enriched in the molecular pathways of actin cytoskeleton regulation, focal adhesion and response to mycophenolic acid.@*CONCLUSIONS@#ASTN1, RHOA and ITGA1 may participate in the development of SIDS and ISDI. The enrichment of differentially expressed genes in immune and inflammatory pathways suggests a common molecular regulatory mechanism between SIDS and ISDI. These findings are expected to provide new biomarkers for molecular anatomy and forensic identification of SIDS and ISDI.

Exploratory research on the probable shared molecular mechanism and transcription factors between chronic periodontitis and chronic obstructive pulmonary disease / 华西口腔医学杂志

OBJECTIVES@#To investigate possible cross-talk genes, associated pathways, and transcription factors between chronic periodontitis (CP) and chronic obstructive pulmonary disease (COPD).@*METHODS@#The gene expression profiles of CP (GSE10334 and GSE16134) and COPD (GSE76925) were downloaded from the GEO database. Differential expression and functional clustering analyses were performed. The protein‑protein interaction (PPI) network was constructed. The core cross-talk genes were filtered using four topological analysis algorithms and modular segmentation. Then, functional clustering analysis was performed again.@*RESULTS@#GSE10334 detected 164 differentially expressed genes (DEGs) (119 upregulated and 45 downregulated). GSE16134 identified 208 DEGs (154 upregulated and 54 downregulated). GSE76925 identified 1 408 DEGs (557 upregulated and 851 downregulated). The PPI network included 21 nodes and 20 edges. The final screening included seven cross-talk genes: CD79A, FCRLA, CD19, IRF4, CD27, SELL, and CXCL13. Relevant pathways included primary immunodeficiency, the B-cell receptor signaling pathway, and cytokine-cytokine receptor interaction.@*CONCLUSIONS@#This study indicates the probability of shared pathophysiology between CP and COPD, and their cross-talk genes, associated pathways, and transcription factors may offer novel concepts for future mechanistic investigations.

Identification and validation of novel biomarkers for cold-dampness syndrome of rheumatoid arthritis based on integration of multiple bioinformatics methods / 中国中药杂志

This study aims to identify the novel biomarkers of cold-dampness syndrome(RA-Cold) of rheumatoid arthritis(RA) by gene set enrichment analysis(GSEA), weighted gene correlation network analysis(WGCNA), and clinical validation. Firstly, transcriptome sequencing was carried out for the whole blood samples from RA-Cold patients, RA patients with other traditional Chinese medicine(TCM) syndromes, and healthy volunteers. The differentially expressed gene(DEG) sets of RA-Cold were screened by comparison with the RA patients with other TCM syndromes and healthy volunteers. Then, GSEA and WGCNA were carried out to screen the key DEGs as candidate biomarkers for RA-Cold. Experimentally, the expression levels of the candidate biomarkers were determined by RT-qPCR for an independent clinical cohort(not less than 10 cases/group), and the clinical efficacy of the candidates was assessed using the receiver operating characteristic(ROC) curve. The results showed that 3 601 DEGs associated with RA-Cold were obtained, including 106 up-regulated genes and 3 495 down-regulated genes. The DEGs of RA-Cold were mainly enriched in the pathways associated with inflammation-immunity regulation, hormone regulation, substance and energy metabolism, cell function regulation, and synovial pannus formation. GSEA and WGCNA showed that recombinant proteasome 26S subunit, ATPase 2(PSMC2), which ranked in the top 50% in terms of coefficient of variation, representativeness of pathway, and biological modules, was a candidate biomarker of RA-Cold. Furthermore, the validation results based on the clinical independent sample set showed that the F1 value, specificity, accuracy, and precision of PSMC2 for RA-Cold were 70.3%, 61.9%, 64.5%, and 81.3%, respectively, and the area under the curve(AUC) value was 0.96. In summary, this study employed the "GSEA-WGCNA-validation" integrated strategy to identify novel biomarkers of RA-Cold, which helped to improve the TCM clinical diagnosis and treatment of core syndromes in RA and provided an experimental basis for TCM syndrome differentiation.

Identification of biomarkers in laryngeal cancer by weighted gene co-expression network analysis / 中南大学学报(医学版)

OBJECTIVES@#Laryngeal cancer (LC) is a globally prevalent and highly lethal tumor. Despite extensive efforts, the underlying mechanisms of LC remain inadequately understood. This study aims to conduct an innovative bioinformatic analysis to identify hub genes that could potentially serve as biomarkers or therapeutic targets in LC.@*METHODS@#We acquired a dataset consisting of 117 LC patient samples, 16 746 LC gene RNA sequencing data points, and 9 clinical features from the Cancer Genome Atlas (TCGA) database in the United States. We employed weighted gene co-expression network analysis (WGCNA) to construct multiple co-expression gene modules. Subsequently, we assessed the correlations between these co-expression modules and clinical features to validate their associations. We also explored the interplay between modules to identify pivotal genes within disease pathways. Finally, we used the Kaplan-Meier plotter to validate the correlation between enriched genes and LC prognosis.@*RESULTS@#WGCNA analysis led to the creation of a total of 16 co-expression gene modules related to LC. Four of these modules (designated as the yellow, magenta, black, and brown modules) exhibited significant correlations with 3 clinical features: The age of initial pathological diagnosis, cancer status, and pathological N stage. Specifically, the yellow and magenta gene modules displayed negative correlations with the age of pathological diagnosis (r=-0.23, P<0.05; r=-0.33, P<0.05), while the black and brown gene modules demonstrated negative associations with cancer status (r=-0.39, P<0.05; r=-0.50, P<0.05). The brown gene module displayed a positive correlation with pathological N stage. Gene Ontology (GO) enrichment analysis identified 77 items, and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis identified 30 related signaling pathways, including the calcium signaling pathway, cytokine-cytokine receptor interaction, neuro active ligand-receptor interaction, and regulation of lipolysis in adipocytes, etc. Consequently, central genes within these modules that were significantly linked to the overall survival rate of LC patients were identified. Central genes included CHRNB4, FOXL2, KCNG1, LOC440173, ADAMTS15, BMP2, FAP, and KIAA1644.@*CONCLUSIONS@#This study, utilizing WGCNA and subsequent validation, pinpointed 8 genes with potential as gene biomarkers for LC. These findings offer valuable references for the clinical diagnosis, prognosis, and treatment of LC.

Construction of predictive ceRNA network and identification of the patterns of immune cells infiltrated in Graves ' ophthalmopathy / 中南大学学报(医学版)

OBJECTIVES@#Graves' ophthalmopathy (GO) is a multifactorial disease, and the mechanism of non coding RNA interactions and inflammatory cell infiltration patterns are not fully understood. This study aims to construct a competing endogenous RNA (ceRNA) network for this disease and clarify the infiltration patterns of inflammatory cells in orbital tissue to further explore the pathogenesis of GO.@*METHODS@#The differentially expressed genes were identified using the GEO2R analysis tool. The Kyoto encyclopedia of genes and genomes (KEGG) and gene ontology analysis were used to analyze differential genes. RNA interaction relationships were extracted from the RNA interactome database. Protein-protein interactions were identified using the STRING database and were visualized using Cytoscape. StarBase, miRcode, and DIANA-LncBase Experimental v.2 were used to construct ceRNA networks together with their interacted non-coding RNA. The CIBERSORT algorithm was used to detect the patterns of infiltrating immune cells in GO using R software.@*RESULTS@#A total of 114 differentially expressed genes for GO and 121 pathways were detected using both the KEGG and gene ontology enrichment analysis. Four hub genes (SRSF6, DDX5, HNRNPC,and HNRNPM) were extracted from protein-protein interaction using cytoHubba in Cytoscape, 104 nodes and 142 edges were extracted, and a ceRNA network was identified (MALAT1-MIR21-DDX5). The results of immune cell analysis showed that in GO, the proportions of CD8+ T cells and CD4+ memory resting T cells were upregulated and downregulated, respectively. The proportion of CD4 memory resting T cells was positively correlated with the expression of MALAT1, MIR21, and DDX5.@*CONCLUSIONS@#This study has constructed a ceRNA regulatory network (MALAT1-MIR21-DDX5) in GO orbital tissue, clarifying the downregulation of the proportion of CD4+ stationary memory T cells and their positive regulatory relationship with ceRNA components, further revealing the pathogenesis of GO.

Construction of endogenous RNA regulatory network for colorectal cancer based on bioinformatics / 中南大学学报(医学版)

OBJECTIVES@#The high morbidity and mortality of colorectal cancer (CRC) have posed great threats to human health. Circular RNA (circRNA) and microRNA (miRNA), acting as competing endogenous RNAs (ceRNAs), have been found to play vital roles in carcinogenesis. This paper aims to construct a circRNA/miRNA/mRNA regulatory network so as to explore the molecular mechanism of CRC.@*METHODS@#The sequencing data of circRNA from CRC were obtained from Gene Expression Omnibus (GEO). The differential circRNA was screened and its structure was identified by Cancer-specific CircRNA Database (CSCD); the sequencing data of miRNA and messenger RNA (mRNAs) were downloaded from The Cancer Genome Atlas (TCGA) database and the differentially expressed genes were screened; the corresponding miRNA of differential circRNAs were predicted by CircInteractome database; DIANA, Miranda, PicTar, and TargetScan databases were used to predict the target genes of different miRNAs; the target genes from Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) were enriched by R language; String database combined with Cytoscape 3.7.2 software was used to construct protein-protein interaction (PPI) network and hub genes were screened; the expressions of mRNAs in the Top10 hub genes were verified in CRC. The network diagrams of circRNAs/miRNAs/mRNAs and circRNAs/miRNAs/Top10 hub mRNAs were constructed by Cytoscape3.7.2. Real-time PCR was used to examine the expression levels of hsa_circRNA_0065173, hsa-mir-450b, hsa-mir-582, adenylate cyclase 5 (ADCY5), muscarinic acetylcholine receptor M2 (CHRM2), cannabinoid receptor 1 (CNR1), and lysophosphatidic acid receptor 1 (LPAR1) in the CRC tissues and the adjacent normal tissues.@*RESULTS@#A total of 14 differential circRNAs were identified, and 8 were found in CSCD; 34 miRNAs targeted by circRNAs were obtained. The PPI network was constructed, and the Top10 hub genes were identified, which were CHRM2, melanin concentrating hormone receptor 2 (MCHR2), G-protein gamma 3 subunit (GNG3), neuropeptide Y receptor Y1 (NPY1R), CNR1, LPAR1, ADCY5, adenylate cyclase 2 (ADCY2), gamma 7 (GNG7) and chemokine 12 (CXCL12), respectively. The expressions of Top 10 hub genes were also verified, and the results showed that the Top 10 hub genes were down-regulated in CRC; the constructed network diagram showed that hsa_circRNA_0065173 may regulate ADCY5, CHRM2, and Hsa-mir-450b by modulating hsa-mir-450b and hsa-mir-582. CNR1 and LPAR1 genes might serve as potentially relevant targets for the treatment of CRC. Real-time PCR results showed that the expression levels of hsa_circRNA_0065173, ADCY5, CHRM2, CNR1 and LPAR1 in the CRC tissues were significantly reduced compared with the adjacent normal tissues (all P<0.05); the expression levels of hsa-mir-450b and hsa-miR-582 were significantly increased (both P<0.05).@*CONCLUSIONS@#In this study, a potential circRNAs/miRNAs/mRNAs network is successfully constructed, which provides a new insight for CRC development mechanism through ceRNA mediated by circRNAs.

Bioinformatics screening and analysis of key differentially expressed genes characteristics in nonalcoholic fatty liver disease / 中华肝脏病杂志

Objective: To screen and analyze the key differentially expressed genes characteristics in nonalcoholic fatty liver disease (NAFLD) with bioinformatics method. Methods: NAFLD-related expression matrix GSE89632 was downloaded from the GEO database. Limma package was used to screen differentially expressed genes (DEGs) in healthy, steatosis (SS), and nonalcoholic steatohepatitis (NASH) samples. WGCNA was used to analyze the output gene module. The intersection of module genes and differential genes was used to determine the differential genes characteristic, and then GO function and KEGG signaling pathway enrichment analysis were performed. The protein-protein interaction network (PPI) was constructed using the online website STRING and Cytoscape software, and the key (Hub) genes were screened. Finally, R software was used to analyze the receiver operating characteristic curve (ROC) of the Hub gene. Results: 92 differentially expressed genes characteristic were obtained through screening, which were mainly enriched in inflammatory response-related functions of "lipopolysaccharide response and molecular response of bacterial origin", as well as cancer signaling pathways of "proteoglycan in cancer" and "T-cell leukemia virus infection-related". 10 hub genes (FOS, CXCL8, SERPINE1, CYR61, THBS1, FOSL1, CCL2, MYC, SOCS3 and ATF3) had good diagnostic value. Conclusion: The differentially expressed hub genes among the 10 NAFLD disease-related characteristics obtained with bioinformatics analysis may become a diagnostic and prognostic marker and potential therapeutic target for NAFLD. However, further basic and clinical studies are needed to validate.

Research advances on the role of competing endogenous RNAs in wound healing / 中华烧伤杂志

Wound healing, as one of the important public health issues, has been a worldwide problem. Due to the unique biological wound environment, wound healing is a very complex process with current treatments requiring long cycles, being poorly effective, and bringing high economic burden to patients. An increasing number of studies have shown that non-coding RNAs (ncRNAs) play important roles in wound healing process. The competing endogenous RNAs (ceRNAs) hypothesis in recent years is a new proposal on the inter-regulation of RNAs, which suggests a "mode of communication" between different RNAs. ceRNA regulatory network (ceRNET) combines the functions of protein-coding mRNA with ncRNA (e.g., microRNA, long non-coding RNA, pseudogenes, and circular RNA). Recent studies have shown that ceRNAs play important roles in wound healing, which may provide new effective therapeutic targets for wound healing. This paper starting with ceRNET systematically reviewed the research progress on the effects of various ceRNAs in wound healing and the future research challenges, with the aim to deeply explore the molecular mechanisms and clinical significance of ceRNAs in the process of wound healing.

Identification of onco-miRNAs in hepatocellular carcinoma and analysis of their regulatory network / 南方医科大学学报

OBJECTIVE@#To construct the regulatory network of survival-related onco-miRNAs and their target genes in hepatocellular carcinoma (HCC) and verify the interactions between the key miRNAs and their targets.@*METHODS@#We screened survival-related miRNAs in HCC in OncomiR and Oncolnc databases, predicted their target genes using miRNet, and conducted survival and expression analysis using GEPIA2 and Ualcan, respectively. The miRNA-target gene co-expression analysis was performed and the miRNA-target network was constructed. Enrichment analysis was performed in Enrichr and protein-protein interaction analysis in STRING database. We tested the effects of transfection with the mimic or inhibitor of hsa-miR-1226-3p or hsa-miR-221-5p on proliferation of HepG2 cells using CCK8 assay and examined the changes in the expressions of the target genes using RT-qPCR. The effect of transfection with hsa-miR-221-5p mimic or inhibitor on protein expressions of the target genes was examined using Western blotting in. A dual luciferase reporter assay was used to test the interaction between hsa-miR-221-5p and its potential target gene GCDH. We further examined the effect of transfection with hsa-miR-221-5p mimic and pEGFP N1-GCDH, alone or in combination, on proliferation, migration and invasion of HepG2 cells.@*RESULTS@#We identified 223 survival-related miRNAs in HCC from OncomiR and 146 miRNAs from Oncolnc with an intersection of 131 miRNAs, and 48 miRNAs were identified as onco-miRNAs in HCC after survival and expression analysis. Twenty-seven eligible target genes were identified after miRNA-mRNA co-expression analysis. The constructed miRNA-target gene network consisted of 25 miRNAs and 27 target genes. The most enriched term was fatty acid metabolism for the target genes. In HepG2 cells, transfection with the mimic or inhibitor of hsa-miR-1226-3p or hsa-miR-221-5p caused significant changes of the mRNA and protein levels of their respective target genes (P < 0.05). The results of dual luciferase reporter assay confirmed the targeting relationship between hsa-miR-221-5p and GCDH gene (P < 0.05). Transfection with hsa-miR-221-5p mimic significantly suppressed the proliferation, migration and invasion of HepG2 cells, but this effect was obviously relieved by co-transformation with pEGFP N1-GCDH (P < 0.05).@*CONCLUSION@#Fatty acid metabolism might be one of the most crucial pathways that mediate the effect of the oncomiRNAs in HCC, and the hsa-miR-221-5p/GCDH axis is an important molecular mechanism for HCC progression.

Therapeutic cell engineering: designing programmable synthetic genetic circuits in mammalian cells

Cell therapy approaches that employ engineered mammalian cells for on-demand production of therapeutic agents in the patient's body are moving beyond proof-of-concept in translational medicine. The therapeutic cells can be customized to sense user-defined signals, process them, and respond in a programmable and predictable way. In this paper, we introduce the available tools and strategies employed to design therapeutic cells. Then, various approaches to control cell behaviors, including open-loop and closed-loop systems, are discussed. We also highlight therapeutic applications of engineered cells for early diagnosis and treatment of various diseases in the clinic and in experimental disease models. Finally, we consider emerging technologies such as digital devices and their potential for incorporation into future cell-based therapies.

Characterization of chromatin accessibility in psoriasis / 医学前沿

The pathological hallmarks of psoriasis involve alterations in T cell genes associated with transcriptional levels, which are determined by chromatin accessibility. However, to what extent these alterations in T cell transcriptional levels recapitulate the epigenetic features of psoriasis remains unknown. Here, we systematically profiled chromatin accessibility on Th1, Th2, Th1-17, Th17, and Treg cells and found that chromatin remodeling contributes significantly to the pathogenesis of the disease. The chromatin remodeling tendency of different subtypes of Th cells were relatively consistent. Next, we profiled chromatin accessibility and transcriptional dynamics on memory Th/Treg cells. In the memory Th cells, 803 increased and 545 decreased chromatin-accessible regions were identified. In the memory Treg cells, 713 increased and 1206 decreased chromatin-accessible regions were identified. A total of 54 and 53 genes were differentially expressed in the peaks associated with the memory Th and Treg cells. FOSL1, SPI1, ATF3, NFKB1, RUNX, ETV4, ERG, FLI1, and ETC1 were identified as regulators in the development of psoriasis. The transcriptional regulatory network showed that NFKB1 and RELA were highly connected and central to the network. NFKB1 regulated the genes of CCL3, CXCL2, and IL1RN. Our results provided candidate transcription factors and a foundational framework of the regulomes of the disease.

Quantifying the state of cell differentiation based on the gene networks entropy / 生物工程学报

Studies of cellular dynamic processes have shown that cells undergo state changes during dynamic processes, controlled mainly by the expression of genes within the cell. With the development of high-throughput sequencing technologies, the availability of large amounts of gene expression data enables the acquisition of true gene expression information of cells at the single-cell level. However, most existing research methods require the use of information beyond gene expression, thus introducing additional complexity and uncertainty. In addition, the prevalence of dropout events hampers the study of cellular dynamics. To this end, we propose an approach named gene interaction network entropy (GINE) to quantify the state of cell differentiation as a means of studying cellular dynamics. Specifically, by constructing a cell-specific network based on the association between genes through the stability of the network, and defining the GINE, the unstable gene expression data is converted into a relatively stable GINE. This method has no additional complexity or uncertainty, and at the same time circumvents the effects of dropout events to a certain extent, allowing for a more reliable characterization of biological processes such as cell fate. This method was applied to study two single-cell RNA-seq datasets, head and neck squamous cell carcinoma and chronic myeloid leukaemia. The GINE method not only effectively distinguishes malignant cells from benign cells and differentiates between different periods of differentiation, but also effectively reflects the disease efficacy process, demonstrating the potential of using GINE to study cellular dynamics. The method aims to explore the dynamic information at the level of single cell disorganization and thus to study the dynamics of biological system processes. The results of this study may provide scientific recommendations for research on cell differentiation, tracking cancer development, and the process of disease response to drugs.

Screening Biomarkers of Sudden Coronary Death Based on mRNA Expression Profile of Rat Myocardial Tissues / 法医学杂志

OBJECTIVES@#To explore the differential expression of messenger RNA (mRNA) in myocardial tissues of rats with sudden coronary death (SCD), and to provide ideas for the forensic identification of SCD.@*METHODS@#The rat SCD model was established, and the transcriptome sequencing was performed by next-generation sequencing technology. Differentially expressed genes (DEGs) in myocardial tissues of SCD rats were screened by using the R package limma. A protein-protein interaction (PPI) network was constructed by using the STRING database and Cytoscape 3.8.2 on DEG, and hub genes were screened based on cytoHubba plug-in. Finally, the R package clusterProfiler was used to analyze the biological function and signal pathway enrichment of the selected DEG.@*RESULTS@#A total of 177 DEGs were associated with SCD and were mainly involved in the renin-angiotensin system and PI3K-Akt signaling pathway. The genes including angiotensinogen (AGT), complement component 4a (C4a), Fos proto-oncogene (FOS) and others played key roles in the development of SCD.@*CONCLUSIONS@#Genes such as AGT, C4a, FOS and other genes are expected to be potential biomarkers for forensic identification of SCD. The study based on mRNA expression profile can provide a reference for forensic identification of SCD.

Integration-based co-expression network analysis to investigate tumor-associated modules across three cancer types / 生物工程学报

In case/control gene expression data, differential expression (DE) represents changes in gene expression levels across various biological conditions, whereas differential co-expression (DC) represents an alteration of correlation coefficients between gene pairs. Both DC and DE genes have been studied extensively in human diseases. However, effective approaches for integrating DC-DE analyses are lacking. Here, we report a novel analytical framework named DC&DEmodule for integrating DC and DE analyses and combining information from multiple case/control expression datasets to identify disease-related gene co-expression modules. This includes activated modules (gaining co-expression and up-regulated in disease) and dysfunctional modules (losing co-expression and down-regulated in disease). By applying this framework to microarray data associated with liver, gastric and colon cancer, we identified two, five and two activated modules and five, five and one dysfunctional module(s), respectively. Compared with the other methods, pathway enrichment analysis demonstrated the superior sensitivity of our method in detecting both known cancer-related pathways and those not previously reported. Moreover, we identified 17, 69, and 11 module hub genes that were activated in three cancers, which included 53 known and three novel cancer prognostic markers. Random forest classifiers trained by the hub genes showed an average of 93% accuracy in differentiating tumor and adjacent normal samples in the TCGA and GEO database. Comparison of the three cancers provided new insights into common and tissue-specific cancer mechanisms. A series of evaluations demonstrated the framework is capable of integrating the rapidly accumulated expression data and facilitating the discovery of dysregulated processes.

Met-Controlled Allosteric Module of Neural Generation as A New Therapeutic Target in Rodent Brain Ischemia / 中国结合医学杂志

OBJECTIVE@#To investigate a Met-controlled allosteric module (AM) of neural generation as a potential therapeutic target for brain ischemia.@*METHODS@#We selected Markov clustering algorithm (MCL) to mine functional modules in the related target networks. According to the topological similarity, one functional module was predicted in the modules of baicalin (BA), jasminoidin (JA), cholic acid (CA), compared with I/R model modules. This functional module included three genes: Inppl1, Met and Dapk3 (IMD). By gene ontology enrichment analysis, biological process related to this functional module was obtained. This functional module participated in generation of neurons. Western blotting was applied to present the compound-dependent regulation of IMD. Co-immunoprecipitation was used to reveal the relationship among the three members. We used IF to determine the number of newborn neurons between compound treatment group and ischemia/reperfusion group. The expressions of vascular endothelial growth factor (VEGF) and matrix metalloproteinase 9 (MMP-9) were supposed to show the changing circumstances for neural generation under cerebral ischemia.@*RESULTS@#Significant reduction in infarction volume and pathological changes were shown in the compound treatment groups compared with the I/R model group (P<0.05). Three nodes in one novel module of IMD were found to exert diverse compound-dependent ischemic-specific excitatory regulatory activities. An anti-ischemic excitatory allosteric module (AM@*CONCLUSIONS@#AM

Changes in gene expression of adipose tissue CD14 / 中南大学学报(医学版)

OBJECTIVES@#To study the gene expression of adipose tissue CD14@*METHODS@#The data of GSE54350 were obtained from the public database of gene expression profiling. The data were pre-processed by Network Analyst, String 11.0, Cytoscape 3.7.1, and other analytical software. The differentially expressed genes were analyzed by gene ontology biological function and kyoto encycopedia of genes and genomes (KEGG) signaling pathway to establish differential gene protein interaction network, transcription factor-gene regulatory network, microRNA-gene regulatory network, environmental factors-gene regulatory network, and other interaction systems.@*RESULTS@#The gene expression pattern of CD14@*CONCLUSIONS@#The gene expression of adipose tissue CD14

Differential mRNA expression in C57BL/6 mice with bleomycin-induced pulmonary fibrosis and its association with LncRNA co-expression network / 南方医科大学学报

OBJECTIVE@#To study the changes in mRNA and long non-coding RNA (lncRNA) expression profiles in a mouse model of bleomycin-induced lung fibrosis and identify lung fibrosis-related mRNA for coding-noncoding coexpression (CNC) bioinformatics analysis of the differential lncRNAs.@*METHODS@#Lung fibrosis was induced by intratracheal injection of bleomycin in 10 C57BL/6 mice and another 10 mice with intratracheal injection of saline served as the control group. Lung tissues were harvested from the mice at 14 days after the injections and lung fibrosis was assessed using Masson and HE staining. LncRNA chip technology was used to screen the differentially expressed mRNAs and lncRNAs in mice with lung fibrosis, and GO and KEGG pathway analyses of the differential mRNAs were performed using NCBI database and UCSC database to identify possible fibrosis-related mRNAs, which were validated by qRT-PCR to construct a coding and non-coding co- expression network with the differential lncRNAs.@*RESULTS@#Compared with the control mice, the mice with intratracheal injection of bleomycin showed obvious lung fibrosis. The results of gene chip analysis showed that 127 mRNAs were upregulated and 184 mRNAs were down-regulated in the model group as compared with the control group. GO and pathway analysis suggested that the differentially expressed genes participated mainly in immune response, cell differentiation, and cytoskeletons; the involved signal pathways were associated mainly with cytokine and cytokine receptor interaction and chemokine signal transduction. Bioinformatics analysis identified a significant coexpression network between the fibrosisrelated mRNA and the differentially expressed lncRNA.@*CONCLUSIONS@#In mice with lung fibrosis, the differential expressions of fibrosis-related mRNAs in the lung tissues are closely correlated with the co- expressions of a large number of differential lncRNAs, which points to a new direction for investigation of the pathogenesis of pulmonary fibrosis.