Expression and prognostic significance of cell division cycle associated protein 5 in pancreatic cancer tissues / 肿瘤研究与临床

Objective:To analyze the expression of cell division cycle associated protein 5 (CDCA5) in pancreatic cancer tissues and its correlation with prognosis based on the bioinformatics.Methods:The RNA sequencing data (HTSeq-FPKM) and corresponding clinical information of 168 pancreatic cancer samples from January to December 2021 were downloaded from the Cancer Genome Atlas (TCGA) database, and the data of 179 pancreatic patients from January to December 2021 were downloaded from the GEPIA2 database, and 171 normal pancreatic tissues from TCGA and GTEx databases were simultaneously integrated. The relative expression level of CDCA5 mRNA in pancreatic cancer patients in GEPIA2 database and its relationship with overall survival (OS) and disease-free survival (DFS) were explored. Combined with the clinical data of the patients, univariate and multivariate Cox regression model analysis was used to analyze the factors influencing the OS of pancreatic cancer patients. Gene set enrichment analysis (GSEA) was performed to investigate the possibly involved signal pathways of CDCA5 in pancreatic cancer.Results:In the GEPIA2 database, the relative expression level of CDCA5 mRNA in pancreatic cancer tissues was higher than that in normal pancreatic tissues, and the difference was statistically significant ( P < 0.05). The pancreatic cancer patients were divided into the high CDCA5 mRNA expression group (89 cases) and the low CDCA5 mRNA expression group (89 cases) according to the median of relative expression level of CDCA5 mRNA (the case equal to the median value was not subgrouped). Survival analysis showed that patients with high CDCA5 mRNA expression had shorter OS ( P = 0.024) and DFS ( P = 0.025) compared with those with low CDCA5 mRNA expression. Multivariate Cox analysis showed that in TCGA database, N staging ( HR = 2.15, 95% CI 1.24-3.72, P = 0.006) and CDCA5 expression ( HR = 1.71, 95% CI 1.23-2.38, P = 0.001) were independent influencing factors of OS for pancreatic cancer patients. The results of GSEA enrichment analysis indicated that high CDCA5 mRNA expression was enriched in 13 biological pathways [all P < 0.05, false discovery rate (FDR) < 0.005] including cell cycle, DNA replication, homologous recombination, pyrimidine metabolism, mismatch repair, pentose phosphate pathway, glycolysis gluconeogenesis and p53. The expression of CDCA5 mRNA was positively correlated with the expressions of HK2, PKM, PGK1, ALDOA, EN01 and LDHA (all P < 0.05). Conclusions:CDCA5 is highly expressed in pancreatic cancer tissues and is associated with poor prognosis of patients, and it can be used as a prognostic marker for pancreatic cancer.

Expression of hsa-miR-422a in hypertrophic scars: a bioinformatics analysis / 中华医学美学美容杂志

Objective:To evaluate the expression level of hsa-miR-422a in hypertrophic scars and to identify the target genes of hsa-miR-422a along with their biological functions using bioinformatics approaches.Methods:From June 2020 to December 2020, tissue samples of 3 hypertrophic scar and 3 normal skin were collected from patients (3 males, 3 females, aged 20-42 years) in Department of Plastic and Reconstructive Surgery, Shanghai Ninth People′s Hospital, Shanghai Jiaotong University School of Medicine. Primary fibroblasts were isolated and cultured. Real-time quantitative PCR was performed to quantify the expression of hsa-miR-422a. To construct a ceRNA network, starbase and Target Scandata bases were utilized to predict genes as well as long noncoding RNAs (lncRNAs) that may sponge hsa-miR-422a. GO and KEGG pathway enrichment analyses were conducted on the target genes of hsa-miR-422a; protein-protein interaction (PPI) networks were constructed to identify the hub genes whose functions were predicted by functional enrichment analyses. The expression of hub genes was validated through real-time quantitative PCR in hypertrophic scars.Results:The expression of hsa-miR-422a was significantly lower in the hypertrophic scar tissue samples and fibroblasts compared to that in the normal skin ( P<0.05). 133 target genes as well as 1033 lncRNAs were predicted by starBase and TargetScandata bases and used to construct an hsa-miR-422a-centered ceRNA network. PPI networks of the target genes revealed 10 hub genes, including MAPK1, GRB2, and IGF1R, which were discovered to be related to protein serine/threonine/tyrosine kinase activity, ubiquitin protein ligase binding, fibroblast growth factor receptor signaling pathway, muscle cell proliferation, and many others; besides, they may be involved in FoxO, mTOR, Toll-like receptor, Ras, MAPK, PI3K-Akt signaling pathways and signaling pathways regulating pluripotency of stem cells. Three hub genes (MAPK1, GRB2, and IGF1R) were significantly upregulated in hypertrophic scars ( P<0.05). Conclusions:hsa-miR-422a is significantly downregulated in the hypertrophic scars and may target hub genes such as MAPK1 in ceRNA networks, ultimately modulating hypertrophic scar formation.

Screening of key targets for diabetic cardiomyopathy based on transcriptomics and proteomics / 中国药理学通报

Aim To explore a potential new target for the prevention and treatment of diabetic cardiomyopathy ( DCM) in mice. Methods The myocardial proteomics of normal and diabetic mice was studied. The GEO database GSE161931 dataset was analyzed using R language with P < 0.05 and I log

Construction of a Prognostic Model of Multiple Myeloma Based on Metabolism-Related Genes / 中国实验血液学杂志

OBJECTIVE@#To screen the prognostic biomarkers of metabolic genes in patients with multiple myeloma (MM), and construct a prognostic model of metabolic genes.@*METHODS@#The histological database related to MM patients was searched. Data from MM patients and healthy controls with complete clinical information were selected for analysis.The second generation sequencing data and clinical information of bone marrow tissue of MM patients and healthy controls were collected from human protein atlas (HPA) and multiple myeloma research foundation (MMRF) databases. The gene set of metabolism-related pathways was extracted from Molecular Signatures Database (MSigDB) by Perl language. The biomarkers related to MM metabolism were screened by difference analysis, univariate Cox risk regression analysis and LASSO regression analysis, and the risk prognostic model and Nomogram were constructed. Risk curve and survival curve were used to verify the grouping effect of the model. Gene set enrichment analysis (GSEA) was used to study the difference of biological pathway enrichment between high risk group and low risk group. Multivariate Cox risk regression analysis was used to verify the independent prognostic ability of risk score.@*RESULTS@#A total of 8 mRNAs which were significantly related to the survival and prognosis of MM patients were obtained (P<0.01). As molecular markers, MM patients could be divided into high-risk group and low-risk group. Survival curve and risk curve showed that the overall survival time of patients in the low-risk group was significantly better than that in the high risk group (P<0.001). GSEA results showed that signal pathways related to basic metabolism, cell differentiation and cell cycle were significantly enriched in the high-risk group, while ribosome and N polysaccharide biosynthesis signaling pathway were more enriched in the low-risk group. Multivariate Cox regression analysis showed that the risk score composed of the eight metabolism-related genes could be used as an independent risk factor for the prognosis of MM patients, and receiver operating characteristic curve (ROC) showed that the molecular signatures of metabolism-related genes had the best predictive effect.@*CONCLUSION@#Metabolism-related pathways play an important role in the pathogenesis and prognosis of patients with MM. The clinical significance of the risk assessment model for patients with MM constructed based on eight metabolism-related core genes needs to be confirmed by further clinical studies.

Research progress on genetic control of host susceptibility to tuberculosis / 浙江大学学报·医学版

The "Lübeck disaster", twins studies, adoptees studies, and other epidemiological observational studies have shown that host genetic factors play a significant role in determining the host susceptibility to Mycobacterium tuberculosis infection and pathogenesis of tuberculosis. From linkage analyses to genome-wide association studies, it has been discovered that human leucocyte antigen (HLA) genes as well as non-HLA genes (such as SLC11A1, VDR, ASAP1 as well as genes encoding cytokines and pattern recognition receptors) are associated with tuberculosis susceptibility. To provide ideas for subsequent studies about risk prediction of MTB infection and the diagnosis and treatment of tuberculosis, we review the research progress on tuberculosis susceptibility related genes in recent years, focusing on the correlation of HLA genes and non-HLA genes with the pathogenesis of tuberculosis. We also report the results of an enrichment analysis of the genes mentioned in the article. Most of these genes appear to be involved in the regulation of immune system and inflammation， and are also closely related to autoimmune diseases.

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.

Progress in stimuli-responsive molecularly imprinted polymers for enrichment analysis of trace components in traditional Chinese medicine / 中国中药杂志

Molecularly imprinted polymers demonstrate outstanding performance in the research on trace ingredients because of their high selectivity. Stimuli-responsive molecularly imprinted polymers(STR-MIPs) with the introduction of different responsive groups on the basis of traditionally imprinted materials can undergo reversible transformations when exposed to external stimuli such as temperature, magnetism, pH or light. Such responsiveness, combined with the specific recognition, endows STR-MIPs with excellent perfor-mance in trace component studies. Traditional Chinese medicine(TCM) contains complex components with trace content, and thus STR-MIPs have broad application prospects in the enrichment analysis of trace components in TCM. This paper elaborates on the application of STR-MIPs in the enrichment analysis of trace components in TCM from the perspectives of different stimuli, summarized relevant research achievements in the recent five years to broaden the application fields of molecular imprinting, and proposed a few opi-nions about their future development.

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.

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.

Bioinformatics Analysis on Key Genes and Immune Infiltration of Osteosarcoma / 中国医学科学院学报

Objective To screen the potential key genes of osteosarcoma by bioinformatics methods and analyze their immune infiltration patterns. Methods The gene expression profiles GSE16088 and GSE12865 associated with osteosarcoma were obtained from the Gene Expression Omnibus(GEO),and the differentially expressed genes(DEGs)related to osteosarcoma were screened by bioinformatics tools.Gene Ontology(GO)annotation,Kyoto Encyclopedia of Genes and Genomes(KEGG)pathway enrichment,and analysis of immune cell infiltration were then carried out for the DEGs.The potential Hub genes of osteosarcoma were identified by protein-protein interaction network,and the expression of Hub genes in osteosarcoma and normal tissue samples was verified via the Cancer Genome Atlas(TCGA). Results A total of 108 DEGs were screened out.GO annotation and KEGG pathway enrichment revealed that the DEGs were mainly involved in integrin binding,extracellular matrix (ECM) structural components,ECM receptor interactions,and phosphatidylinositol 3-kinase/protein kinase B(PI3K/Akt)signaling pathway.Macrophages were the predominant infiltrating immune cells in osteosarcoma.Secreted phosphoprotein 1(SPP1),matrix metallopeptidase 2(MMP2),lysyl oxidase(LOX),collagen type V alpha(II)chain(COL5A2),and melanoma cell adhesion molecule(MCAM)presented differential expression between osteosarcoma and normal tissue samples(all P<0.05). Conclusions SPP1,MMP2,LOX,COL5A2,and MCAM are all up-regulated in osteosarcoma,which may serve as potential biomarkers of osteosarcoma.Macrophages are the key infiltrating immune cells in osteosarcoma,which may provide new perspectives for the treatment of osteosarcoma.

Differential gene sequencing alignment analysis of hyperplastic stenosis in murine arteriovenous fistula / 中华肾脏病杂志

Objective:To establish a mouse model of intra-jugular arteriovenous fistula (AVF) to screen differentially expressed genes in the process of intimal stenosis of AVF for investigating the abnormal expression signaling pathways and the mechanisms.Methods:Forty-six male C57BL/6 mice were randomly divided into AVF group ( n=23) and sham-operated group ( n=23). The AVF group underwent internal jugular arteriovenous fistuloplasty, and the sham-operated group separated the right external jugular vein and common carotid artery and then sutured the incision. The whole-genome sequences of mice with AVF stenosis were determined by transcriptomic reversible chain terminator and synthetic sequencing. The microarray data set was established, and the Benjamini & Hochberg method of gene microarray data analysis was applied to screen the differentially expressed genes. The differentially expressed genes were screened by R-language enrichment analysis. Then, gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) were performed. The subcellular localization of the differentially expressed genes was performed by BUSCA software. The protein network interaction of differentially expressed genes was analyzed by using STRING database and Cytoscape software. Results:In the AVF group, 21 mice were successfully modeled and 2 mice failed. Therefore, there were 21 mice in the AVF group and only 21 mice in the sham-operated group. This mouse internal jugular AVF model was innovated using the continuous-interrupted suture method, which improved the success rate of modeling this model. The differential gene sequencing analysis showed that there were 2 514 differentially expressed genes in the AVF process, including 1 323 up-regulated genes and 1 191 down-regulated genes. GO functional enrichment analysis showed that the differential genes were mainly enriched in metabolic process, activation, redox, mitochondria and so on. KEGG pathway enrichment analysis showed that the differential genes were enriched in metabolism, energy substance synthesis, diabetes, oxidative stress and so on. Statistical analysis of subcellular localization showed that the differences were mainly in mitochondrial proteins (24.24%), cytoplasmic proteins (17.51%), nuclear proteins (13.13%), cell membrane proteins (11.45%), and extracellular proteins (10.77%).Conclusions:Mitochondrial oxidative stress injury may be involved in the pathological damage process of endothelial proliferation stenosis in the AVF.

Macrophage differential genes and their functions in silicosis based on single cell RNA sequencing and spatial transcriptome sequencing / 环境与职业医学

Background The rise of single cell RNA sequencing (scRNA-seq) and spatial transcriptome sequencing technologies has allowed for intensive study of lung diseases, but both have been poorly studied in silicosis. Objective To explore differentially expressed genes DEGs in silicosis macrophages by scRNA-seq combined with spatial transcriptome sequencing and analyze the potential diagnostic genes. Methods Male C57BL/6 mice (5-6 weeks old, 22-30 g) were randomly divided into 4 groups: normal saline (NS) group for 7 d, NS group for 56 d, SiO2 group for 7 d, and SiO2 group for 56 d, with 1 mouse in each group. A silicosis model was constructed by tracheal drip injection of SiO2 suspension (0.2 g·kg−1, 50 g·cm−2), and the control mice were given the same volume of NS. The right lung was removed for scRNA-seq and the left lung for spatial transcriptome sequencing on day 7 and day 56, respectively. Cell populations were captured using principal component analysis techniques and dimensionality reduction of uniform manifold approximation and projection. The Find Markers function in R language was applied to analyze the DEGs changes of macrophages in two groups of lung tissues, and the corresponding DEGs were subjected to Gene Ontology enrichment analysis and Kyoto Encyclopedia of Genes and Genomes signaling pathway analysis, while STRING and CytoHubba plug-ins of Cytoscape software were applied to protein-protein interaction network analysis to screen out key (Hub) genes. Spatial transcriptome sequencing was used to explore the original location of Hub genes on lung tissue sections and their mapping in lung macrophages. Finally, the correlation of Hub gene expression levels in lung tissues of silicosis patients and mouse silicosis models was verified, the diagnostic efficacy of Hub gene using subject operating characteristic curves (ROC). In vitro experiments by applying cell viability assay were conducted to verify the changes in viability of mouse macrophages (RAW264.7) under SiO2 stimulation. Results The scRNA-seq revealed a total of 20 clusters captured and defined. The results of scRNA-seq and spatial transcriptome sequencing showed an increased number of macrophages in the lung tissue of the SiO2 group compared to the NS group and clustered in the focal areas. Among the 97 macrophage DEGs screened out, 75 were up-regulated genes, and mainly enriched in chemotaxis and migration of neutrophils, chemokine receptor binding, tumor necrosis factor signaling pathway, cytokine-cytokine receptor interaction pathway, and interleukin-17 signaling pathway; and 22 were down-regulated genes, and mainly enriched in late endosomes, peroxisome proliferator-activated receptors signaling pathway, and alcoholic liver disease signaling pathway. A total of 2 core modules and 3 Hub genes were screened out, including Ccl2, Ccl7, and Ptgs2. The scRNA-seq showed that they were expressed at elevated levels in the SiO2 group compared to the NS group and clustered in additional macrophages, and the spatial transcriptome sequencing showed that they clustered in inflammatory areas with nodular lesions. The CCL7 and PTGS2 expressions were increased in the lung tissue of SiO2 patients compared with the healthy subjects, and the areas under the working curve of the subjects were 0.850 and 0.786, respectively. The viability of RAW264.7 cells was enhanced under SiO2 stimulation at 3 h, 6 h, and 12 h compared to those without the stimulation (P<0.05). Conclusion Bioinformatics screening have identified 3 Hub genes (Ccl2, Ccl7, and Ptgs2)and 2 potential diagnostic genes (CCL7 and PTGS2) in the lung tissue of silicosis mice, which may be potential molecular markers of early-stage silicosis with implications for the development and prognosis of silicosis.

Bioinformatics analysis of gene expression profile of peripheral blood mononuclear cells in female patients with osteoarthritis / 西安交通大学学报(医学版)

【Objective】 Through bioinformatics methods to analyze the differences in the gene expression profiles of peripheral blood mononuclear cells （PBMCs） between middle-aged and elderly women and normal people， so as to explore the diagnosis and treatment targets of OA. 【Methods】 We downloaded the GSE48556 data set from GEO databases. We utilized the R language to screen out the differentially expressed genes （DEGs） between OA and NC. By gene set enrichment analysis （GSEA）， we obtained the target gene subset. The GO and KEGG pathways of the target gene subset were analyzed by DAVID. We applied STRING and Cytoscape software to construct PPI network. The module analysis was performed by the Mcode and centiscape plug-in， and the key genes were screened out by Cytohubba. 【Results】 By GSEA analysis and P.adjust 0.2， a total of 292 target genes were screened， consisting of 81 upregulated genes and 211 downregulated genes. The GO enrichment analysis of all target genes mainly focused on the biological functions， such as “regulation of NIK/NF-κB”， “monocytes”， “proliferation”， “regulation of apoptosis signaling pathway”， “TNF-mediated signaling pathway”， “regulation of Wnt signaling pathway”， “regulation of MAP kinase activity”， and “regulation of autophagy”. KEGG was mainly enriched in four pathways： cytotoxicity of natural killer cell mediation， TNF signaling pathway， MAPK signaling pathway， and apoptosis. We employed PPI network and related plug-ins to screen out eight core genes highly related to OA inflammation and apoptosis， namely， MAPK1， IL10， PTGS2， IL18， GSK3B， NFKBIA， TNFRSF1A， and EGR1. 【Conclusion】 Bioinformatics analysis revealed that the differences in PBMCs gene expressions between OA and NC were concentrated in the biological events of apoptosis and inflammation， making blood expression profile an effective breakthrough for monitoring OA target markers.

Identification of key genes and pathways associated with dilated cardiomyopathy development by bioinformatics analysis / 西安交通大学学报(医学版)

【Objective】 To explore the potential biomarkers and related enrichment pathways of dilated cardiomyopathy （DCM） by bioinformatics methods. 【Methods】 The data sets related to DCM in GEO database were searched， and microarray data sets GSE42955 and GSE1869 of human cardiomyocytes were extracted. Then， the differentially expressed genes （DEGS） were analyzed using R language， and the protein-protein interaction network was analyzed to identify the core genes and core modules of differential expression. The gene ontology database （GO） enrichment and Kyoto Gene and Genome Encyclopedia （KEGG） pathway analyses were performed. The data sets related to DCM in ArrayExpress database were searched， and the human cardiomyocytes microarray data set E-TABM-480 was extracted to verify the expressions of core genes and modules. 【Results】 We identified 10 DEGS， namely， DZIP3， FBXO32， BTBD6， FBXL5， ASB8， COMMD1， LTN1， FBXO21， RCHY1 and ARIH2， and the core DEG was DZIP3. After GO and KEGG analyses， the GO and KEGG of the above DEGS were mainly related to the ubiquitin-proteasome system. 【Conclusion】 Bioinformatics analysis shows that the ubiquitin-proteasome system plays an important role in the pathogenesis of DCM， and the mechanism remains to be further studied.

Exploring the association between de novo mutations and non-syndromic cleft lip with or without palate based on whole exome sequencing of case-parent trios / 北京大学学报(医学版)

OBJECTIVE@#To explore the association between de novo mutations (DNM) and non-syndromic cleft lip with or without palate (NSCL/P) using case-parent trio design.@*METHODS@#Whole-exome sequencing was conducted for twenty-two NSCL/P trios and Genome Analysis ToolKit (GATK) was used to identify DNM by comparing the alleles of the cases and their parents. Information of predictable functions was annotated to the locus with SnpEff. Enrichment analysis for DNM was conducted to test the difference between the actual number and the expected number of DNM, and to explore whether there were genes with more DNM than expected. NSCL/P-related genes indicated by previous studies with solid evidence were selected by literature reviewing. Protein-protein interactions analysis was conducted among the genes with protein-altering DNM and NSCL/P-related genes. R package "denovolyzeR" was used for the enrichment analysis (Bonferroni correction: P=0.05/n, n is the number of genes in the whole genome range). Protein-protein interactions among genes with DNM and genes with solid evidence on the risk factors of NSCL/P were predicted depending on the information provided by STRING database.@*RESULTS@#A total of 339 908 SNPs were qualified for the subsequent analysis after quality control. The number of high confident DNM identified by GATK was 345. Among those DNM, forty-four DNM were missense mutations, one DNM was nonsense mutation, two DNM were splicing site mutations, twenty DNM were synonymous mutations and others were located in intron or intergenic regions. The results of enrichment analysis showed that the number of protein-altering DNM on the exome regions was larger than expected (P < 0.05), and five genes (KRTCAP2, HMCN2, ANKRD36C, ADGRL2 and DIPK2A) had more DNM than expected (P < 0.05/(2×19 618)). Protein-protein interaction analysis was conducted among forty-six genes with protein-altering DNM and thirteen genes associated with NSCL/P selected by literature reviewing. Six pairs of interactions occurred between the genes with DNM and known NSCL/P-related genes. The score measuring the confidence level of the predicted interaction between RGPD4 and SUMO1 was 0.868, which was higher than the scores for other pairs of genes.@*CONCLUSION@#Our study provided novel insights into the development of NSCL/P and demonstrated that functional analyses of genes carrying DNM were warranted to understand the genetic architecture of complex diseases.

Clinical implication of xenotropic and polytropic retrovirus receptor 1 in papillary thyroid carcinoma / 浙江大学学报·医学版

To investigate the expression of xenotropic and polytropic retrovirus receptor 1 （） in papillary thyroid cancer （PTC） and its clinical implication. The HPA and UALCAN databases were used to explore the expression of XPR1 in PTC and normal tissues. The cBioPortal database was used to obtain the clinical data of PTC patients and gene expression profile. The correlation of expression with gender，age，sub-types，T stage，N stage，M stage and clinical stage of patients were analyzed. Cox regression was conducted to analysis the factors affecting the prognosis of PTC patients. The mutation of was assessed through cBioPortal database. GO and KEGG analyses were used to explore the related biological pathway of involved in PTC. HPA database analysis showed that XPR1 was highly expressed in PTC tissue compared with normal tissues. UALCAN analysis displayed that expression was significantly higher in PTC tissue compared with normal tissues （0.05）. Cox regression analysis showed that was an independent prognostic factor of PTC patients （=2.894，<0.05）. The cBioPortal database indicated that the mutation appeared in 6% PTC patients; the mutation type mainly was missense and the mutation point was located at the E615K. Enrichment analysis indicated that might affect the PTC progression through involvement in metabolic pathway. is highly expressed in PTC tissues，which is associated with the prognosis of patients. Metabolic pathway associated with might play an important role in PTC progression，indicating that might be a novel biomarker for diagnosis and treatment of PTC.

Study on active components and targets of Ginseng in Prevention and treatment of PTSD based on Network Pharmacology / 中国药理学通报

Aim To obtain the active components and targets of ginseng in the prevention and treatment of traumatic stress disorder(PTSD) through the method of network pharmacology. Methods The active components and target information of ginseng with medicinal value were obtained by TCMSP research platform, and the gene information closely related to the pathogenesis of PTSD was obtained by searching GeneCard and OMIM database. The two were matched to obtain the medicinal components and target genes of ginseng in the prevention and treatment of PTSD. The drug-dis- ease-target network diagram was drawn by R and Perl computer languages, and the target genes were analyzed by PPI network analysis, gene ontology ( GO ) and signal transduction pathway ( KEGG) enrichment analysis. Results According to the general pharmacological research methods of traditional Chinese medi cine, the screening parameters of active components were set, and nine kinds of high value medicinal ingredients of Panax ginseng were obtained. There was a drug-target relationship between the nine medicinal components and sixteen target genes related to PTSD disease. Through PPI, GO and KEGG analysis, it was found that the target genes were mainly enriched in physiological functions such as neurotransmitters, syn-aptic plasticity, ion channels and so on. Conclusions Ginseng has the pharmacological effect of preventing and treating PTSD, which may play a role in regulating the metabolism and receptor activity of monoamine neurotransmitters.

Pathogenesis and therapeutic targets of inflammatory cardiomyopathy based on bioinformatics analysis / 西安交通大学学报(医学版)

【Objective】 To make bioinformatics analysis of inflammatory cardiomyopathy so as to screen out hub genes related to etiology and therapeutic targets. 【Methods】 Differential expression analysis of inflammatory cardiomyopathy gene chip data from Gene Expression Omnibus （GEO） Database was carried out via GEO2R tool. Protein-protein interaction（PPI）network and hub genes identification were realized by String database and CytoHubba. GO and KEGG enrichment analysis for functional annotation and pathway analysis of hub genes were conducted by R language. Web-based enrichment analysis platform Enrichr and Drug Signatures database were applied to screen out candidate drugs targeting hub genes for inflammatory cardiomyopathy. 【Results】 The 149 DEGs were statistically significant， among which 44 were upregulated and 105 were downregulated. To identify hub genes， PPI network consisting of 37 nodes and 116 edges was constructed， and 16 hub genes were NDUFB7， POLR2L， NDUFS7， UQCR11， NDUFA13， NDUFA2， PHPT1， NDUFB10， UBA52， ATP5D， NDUFA3， COX6B1， POLR2J， COX4I2， AURKAIP1 and MRPL41. Hub genes were enriched to 113 different GO terms， and the most significant terms were mitochondrial ATP synthesis coupled electron transport， respiratory electron transport chain， oxidative phosphorylation， respiratory chain， mitochondrial inner membrane， NADH dehydrogenase activity and oxidoreductase activity. DEGs were enriched to 13 different signal pathways， including oxidative phosphorylation， non-alcoholic fatty liver disease， diabetic cardiomyopathy， and cardiac muscle contraction. We screened out candidate drugs targeting hub genes， namely， metformin hydrochloride， clindamycin， and hydralazine. 【Conclusion】 Hub genes screened out by decoding the expression profiles are convolved in the etiology and mechanism of inflammatory cardiomyopathy， which might serve as latent therapeutic targets and benefit patients with inflammatory cardiomyopathy.

Expression analysis of miRNA profiles in apheresis platelets at the end of storage under riboflavin and ultraviolet-Blight treatment / 中国输血杂志

【Objective】 To investigate the changes of platelet microRNA (miRNA) expression profiles on d1 and d5 during storage with riboflavin and ultraviolet-B (UVB) light (VB2-PRT) treatment, and to explore the molecular mechanism of miRNAs involved in the regulation of platelet storage lesion (PSL) under VB2-PRT treatment. 【Methods】 20 apheresis platelet concentrates (5mL / sample) were collected from voluntary blood donors. After mixing and shaking, the samples was treated with riboflavin (final concentration 50 μmol/L) and 6.24J/mL UVB light for 8min, then split into two aliquots and agitated stored at (22±2) ℃. The concentrates were sampled (5mL) on d1 and d5, respectively, and sequenced by DNA nanoball (DNB) sequencing technology. The differentially expressed miRNAs between the two groups (at different storage periods) were screened by DEGseq and MA-plot analysis software. The miRNAs, reached more than 2 times different expression between groups, were considered significant different(P<0.01). The miRanda and TargetScan softwares were used to predict the target genes. Gene Ontology (GO) function enrichment analysis and Kyoto Encyclopedia of genes and genomes (KEGG) pathway enrichment analysis were performed on the target genes of significant differentially expressed miRNAs. The expression of miRNAs was verified by real-time fluorescence quantitative PCR (qRT-PCR). 【Results】 miRNA expression profile: compared with d1 platelets, there were 590 miRNAs with significantly different expression (P< 0.01) in d5 group, including 255 up-regulated miRNAs (such as miR-99b, miR-7) and 335 down regulated miRNAs (such as miR-451a, miR-19b). Among the 272 known miRNAs, 112 were up-regulated and 160 were down regulated. There were 318 new miRNAs sequences. The enriched GO terms of target genes of differentially expressed miRNAs in d5 and d1 groups included cell components, organelles, cell membranes and other cellular structures, molecular functions such as adhesion, catalysis, molecular conversion, transportation, transcription factor and receptor activity, and biological processes such as cell processing, metabolism, biological regulation and stress. The corresponding pathways in the top 10 of KEGG enrichment were mainly secretion, glucose metabolism, signal transduction, membrane transport, translation, environmental adaptation and other signal pathways. The six randomly selected differentially expressed miRNAs verified by qRT-PCR was consistent with those of DNB sequencing. 【Conclusion】 The expression profiles of platelets miRNAs has changed significantly between d1- and d5-storage under VB2-PRT treatment. Functional prediction suggests that these miRNAs might be involved in the regulation of platelet PSL underVB2-PRT treatment.

Bioinformatic analysis between epithelioid and mixed uveal melanoma / 国际眼科杂志(Guoji Yanke Zazhi)

@#AIM: To explore the differentially expressed genes and crucial genes between epithelioid and mixed uveal melanoma(UM)based on bioinformatics analysis.<p>METHODS: Microarray datasets GSE22138 was extracted from gene expression omnibus database(GEO). The differentially expressed genes(DEGs)were screened out between epithelioid and mixed UM, and functional enrichment analysis were performed with DAVID database. STRING and cytoscape was applied to explore the protein-protein interaction(PPI)network and hub genes. Subsequently, cBioPortal was applied to explore the network of the hub genes, and GEPIA was adopted to study the survival analysis of hub genes.<p>RESULTS: Overall, 241 DEGs including 125 upregulated and 116 down regulated genes were identified. The DEGs mainly enriched in cell adhesion, response to drug and Positive regulation of endothelial cell proliferation. A total of 10 hub genes were identified. Survival analysis revealed the hub genes was associated with the prognosis of UM.<p>CONCLUSION: DEGs and hub genes identified by Bioinformatics analysis in the present study would be beneficial to understand mechanism and biological characteristics between epithelioid and mixed UM.