Mechanisms of brain damage caused by inorganic fluoride using proteomics-based techniques / 环境与职业医学

Background Chronic excessive exposure to fluoride can cause damage to the central nervous system and a certain degree of learning and memory impairment. However, the associated mechanism is not yet clear and further exploration is needed. Objective Using 4D unlabelled quantitative proteomics techniques to explore differentially expressed proteins and their potential mechanisms of action in chronic excessive fluoride exposure induced brain injury. Methods Twenty-four SPF-grade adult SD rats, half male and half male, were selected and divided into a control group and a fluoride group by random number table method, with 12 rats in each group. Among them, the control group drank tap water (fluorine content<1 mg·L−1), the fluoride group drank sodium fluoride solution (fluorine content 10 mg·L−1), and both groups were fed with ordinary mouse feed (fluoride content<0.6 mg·kg−1). After 180 d of feeding, the SD rats were weighed, and then part of the brain tissue was sampled for pathological examination by hematoxylin-eosin (HE) staining and Nissl staining. The rest of the brain tissue was frozen and stored at −80 ℃. Three brain tissue samples from each group were randomly selected for proteomics detection. Differentially expressed proteins were screened and subcellular localization analysis was performed, followed by Gene Ontology (GO) function analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, cluster analysis, and protein-protein interaction analysis. Finally, Western blotting was used to detect the expression levels of key proteins extracted from the brain tissue samples. Results After 180 d of feeding, the average weight of the rats in the fluoride group was significantly lower than that in the control group (P<0.05). The brain tissue stained with HE showed no significant morphological changes in the cerebral cortex of the fluoride treated rats, and neuron loss, irregular arrangement of neurons, eosinophilic changes, and cell body pyknosis were observed in the hippocampus. The Nissl staining results showed that the staining of neurons in the cerebral cortex and hippocampus of rats exposed to fluoride decreased (Nissl bodies decreased). The proteomics results showed that a total of 6927 proteins were identified. After screening, 206 differentially expressed proteins were obtained between the control group and the fluoride group, including 96 up-regulated proteins and 110 down-regulated proteins. The differential proteins were mainly located in cytoplasm (30.6%), nucleus (27.2%), mitochondria (13.6%), plasma membrane (13.6%), and extracellular domain (11.7%). The GO analysis results showed that differentially expressed proteins mainly participated in biological processes such as iron ion transport, regulation of dopamine neuron differentiation, and negative regulation of respiratory burst in inflammatory response, exercised molecular functions such as ferrous binding, iron oxidase activity, and cytokine activity, and were located in the smooth endoplasmic reticulum membrane, fixed components of the membrane, chloride channel complexes, and other cellular components. The KEGG significantly enriched pathways included biosynthesis of secondary metabolites, carbon metabolism, and microbial metabolism in diverse environments. The results of differential protein-protein interaction analysis showed that the highest connectivity was found in glucose-6-phosphate isomerase (Gpi). The expression level of Gpi in the brain tissue of the rats in the fluoride group was lower than that in the control group by Western blotting (P<0.05). Conclusion Multiple differentially expressed proteins are present in the brain tissue of rats with chronic fluorosis, and their functions are related to biosynthesis of secondary metabolites, carbon metabolism, and microbial metabolism in diverse environments; Gpi may be involved in cerebral neurological damage caused by chronic overdose fluoride exposure.

Identification of overlay differentially expressed genes in both rats and goats with blast lung injury through comparative transcriptomics / 中华创伤杂志(英文版)

PURPOSE@#To identify the potential target genes of blast lung injury (BLI) for the diagnosis and treatment.@*METHODS@#This is an experimental study. The BLI models in rats and goats were established by conducting a fuel-air explosive power test in an unobstructed environment, which was subsequently validated through hematoxylin-eosin staining. Transcriptome sequencing was performed on lung tissues from both goats and rats. Differentially expressed genes were identified using the criteria of q ≤ 0.05 and |log2 fold change| ≥ 1. Following that, enrichment analyses were conducted for gene ontology and the Kyoto Encyclopedia of Genes and Genomes pathways. The potential target genes were further confirmed through quantitative real-time polymerase chain reaction and enzyme linked immunosorbent assay.@*RESULTS@#Observations through microscopy unveiled the presence of reddish edema fluid, erythrocytes, and instances of focal or patchy bleeding within the alveolar cavity. Transcriptome sequencing analysis identified a total of 83 differentially expressed genes in both rats and goats. Notably, 49 genes exhibited a consistent expression pattern, with 38 genes displaying up-regulation and 11 genes demonstrating down-regulation. Enrichment analysis highlighted the potential involvement of the interleukin-17 signaling pathway and vascular smooth muscle contraction pathway in the underlying mechanism of BLI. Furthermore, the experimental findings in both goats and rats demonstrated a strong association between BLI and several key genes, including anterior gradient 2, ankyrin repeat domain 65, bactericidal/permeability-increasing fold containing family A member 1, bactericidal/permeability-increasing fold containing family B member 1, and keratin 4, which exhibited up-regulation.@*CONCLUSIONS@#Anterior gradient 2, ankyrin repeat domain 65, bactericidal/permeability-increasing fold containing family A member 1, bactericidal/permeability-increasing fold containing family B member 1, and keratin 4 hold potential as target genes for the prognosis, diagnosis, and treatment of BLI.

Integrated aqueous humor ceRNA and miRNA–TF–mRNA network analysis reveals potential molecular mechanisms governing primary open-angle glaucoma pathogenesis

To conduct an integrated bioinformatics analysis of extant aqueous humor (AH) gene expression datasets in order to identify key genes and the regulatory mechanism governing primary open?angle glaucoma (POAG) progression. Methods: Two datasets (GSE101727 and GSE105269) were downloaded from the Gene Expression Omnibus, and the messenger RNAs (mRNAs), microRNAs (miRNAs), and long noncoding RNAs (lncRNAs) were identified between controls and POAG patients. Differentially expressed (DE) mRNAs and DElncRNAs were then subjected to pathway enrichment analyses, after which a protein–protein interaction (PPI) network was generated. This network was then expanded to establish lncRNA–miRNA–mRNA and miRNA–transcription factor (TF)–mRNA networks. Results: The GSE101727 dataset was used to identify 2746 DElncRNAs and 2208 DEmRNAs, while the GSE105269 dataset was used to identify 45 DEmiRNAs. We ultimately constructed a competing endogenous RNA (ceRNA) network incorporating 47 lncRNAs, six miRNAs, and 17 mRNAs. The proteins encoded by these 17 hub mRNAs were found to be significantly enriched for activities that may be linked to POAG pathogenesis. In addition, we generated a miRNA–TF–mRNA regulatory network containing two miRNAs (miR?135a?5p and miR?139?5p), five TFs (TGIF2, TCF3, FOS, and so on), and five mRNAs (SHISA7, ST6GAL2, TXNIP, and so on). Conclusion: The SHISA7, ST6GAL2, TXNIP, FOS, and DCBLD2 genes may be viable therapeutic targets for the prevention or treatment of POAG and are regulated by the TFs (TGIF2, HNF1A, TCF3, and FOS)

Procine recombinant NK-lysin inhibits hepatocellular carcinoma metastasis by downregulating FKBP3 and inhibiting oxidative phosphorylation and glycolysis: a proteomic analysis / 南方医科大学学报

OBJECTIVE@#To investigate the potential mechanisms that mediate the inhibitory effect of porcine recombinant NKlysin (prNK-lysin) against liver cancer cell metastasis.@*METHODS@#HPLC-tandem mass spectrometry was used to identify the differentially expressed proteins in prNK-lysin-treated hepatocellular carcinoma SMMOL/LC-7721 cells in comparison with the control and PBS-treated cells. GO functional annotation and KEGG pathway analysis of the differentially expressed proteins were performed using GO and KEGG databases. RT-qPCR was used to determine the mRNA expression levels of polypeptide-N-acetylgalactosaminotransferase 13 (GALNT13), transmembrane protein 51 (TMEM51) and FKBP prolyl isomerase 3 (FKBP3) in the cells, and the protein expression of FKBP3 was verified using Western blotting.@*RESULTS@#Proteomic analysis identified 1989 differentially expressed proteins in prNK-lysin-treated cells compared with the control cells, and 2753 compared with PBS-treated cells. Fifteen proteins were differentially expressed between PBS-treated and the control cells, and 1909 were differentially expressed in prNK- lysin group compared with both PBS and control groups. These differentially expressed proteins were involved mainly in the viral process, translational initiation and RNA binding and were enriched mainly in ribosome, protein process in endoplasmic reticulum, and RNA transport pathways. RT-qPCR showed that compared with the control group, prNK-lysin treatment significantly increased the mRNA expressions of GALNT13 (P < 0.05) and TMEM51 (P < 0.01) and lowered FKBP3 mRNA expression (P < 0.05). Western blotting also showed a significantly decreased expression of FKBP3 protein in prNK-lysin-treated cells (P < 0.001).@*CONCLUSION@#Treatment with prNK-lysin causes significant changes in protein expression profile of SMMOL/LC-7721 cells and inhibits hepatocellular carcinoma metastasis by downregulating FKBP3 protein and affecting the cellular oxidative phosphorylation and glycolysis pathways.

Comparative analysis of differentially expressed genes for biosynthesis of active ingredients in fruits of different cultivars of Lycium barbarum L. based on transcriptome sequencing / 生物工程学报

To explore the differentially expressed genes (DEGs) related to biosynthesis of active ingredients in wolfberry fruits of different varieties of Lycium barbarum L. and reveal the molecular mechanism of the differences of active ingredients, we utilized Illumina NovaSeq 6000 high-throughput sequencing technology to conduct transcriptome sequencing on the fruits of 'Ningqi No.1' and 'Ningqi No.7' during the green fruit stage, color turning stage and maturity stage. Subsequently, we compared the profiles of related gene expression in the fruits of the two varieties at different development stages. The results showed that a total of 811 818 178 clean reads were obtained, resulting in 121.76 Gb of valid data. There were 2 827, 2 552 and 2 311 DEGs obtained during the green fruit stage, color turning stage and maturity stage of 'Ningqi No. 1' and 'Ningqi No. 7', respectively, among which 2 153, 2 050 and 1 825 genes were annotated in six databases, including gene ontology (GO), Kyoto encyclopedia of genes and genomes (KEGG) and clusters of orthologous groups of proteins (KOG). In GO database, 1 307, 865 and 624 DEGs of green fruit stage, color turning stage and maturity stage were found to be enriched in biological processes, cell components and molecular functions, respectively. In the KEGG database, the DEGs at three developmental stages were mainly concentrated in metabolic pathways, biosynthesis of secondary metabolites and plant-pathogen interaction. In KOG database, 1 775, 1 751 and 1 541 DEGs were annotated at three developmental stages, respectively. Searching the annotated genes against the PubMed database revealed 18, 26 and 24 DEGs related to the synthesis of active ingredients were mined at the green fruit stage, color turning stage and maturity stage, respectively. These genes are involved in carotenoid, flavonoid, terpenoid, alkaloid, vitamin metabolic pathways, etc. Seven DEGs were verified by RT-qPCR, which showed consistent results with transcriptome sequencing. This study provides preliminary evidences for the differences in the content of active ingredients in different Lycium barbarum L. varieties from the transcriptional level. These evidences may facilitate further exploring the key genes for active ingredients biosynthesis in Lycium barbarum L. and analyzing their expression regulation mechanism.

Gene expression profile in human cervical epithelial cells infected with Chlamydia trachomatis / 中华微生物学和免疫学杂志

Objective:To compare gene expression profiles in normal human cervical epithelial cells (HcerEpic) before and after Chlamydia trachomatis ( Ct) infection. Methods:HcerEpic cells that were pretreated with DEAE-D were infected with Ct serotype E standard strain and then cultured for 44 h. Uninfected HcerEpic cells were used as the control group. Total RNA was extracted from the cells in each group and reverse transcribed to construct a cDNA library. Differences in gene expression profiles between the two groups were analyzed by high-throughput sequencing and the representative genes were selected for verification by qPCR. Results:A total of 23 997 genes were detected, including 125 differentially expressed genes. Among the 125 genes, 119 were up-regulated and six were down-regulated. GO analysis showed that the differentially expressed genes were enriched in several biological processes including defense response to virus, typeⅠinterferon signaling pathway and cellular responses to typeⅠinterferons. KEGG enrichment analysis showed the differentially expressed genes were mainly enriched in the pathways related to virus infections, such as influenza A virus, herpes simplex virus, EB virus and HPV, and NOD-like receptor pathway.Conclusions:There were significant differences in transcriptome profiles of HcerEpic cells before and after Ct infection. The differentially expressed genes were mainly enriched in the interferon pathway, which was closely related to the antiviral processes in cells. qPCR verified the differentially expressed genes and the genes closely related to the interferon pathway, such as ISG15, IFIT2, OASL and UBE2L6.

Search for potential target genes in dexamethasone-induced open-angle glaucoma by bioinformatics / 中华实验眼科杂志

Objective:To predict potential target genes in dexamethasone-induced open-angle glaucoma via bioinformatics technology.Methods:The GEO datasets GSE16643, GSE37474, and GSE124114 were used to analyze the differentially expressed genes by GEO2R.Gene Set Enrichment Analysis (GSEA) was performed on the differentially expressed genes between GSE37474 and GSE124114.Intersection of the three datasets were displayed by Venn diagram.The annotation and enrichment analysis of the intersection genes were performed through Gene Ontology (GO)/Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis and then were compared with normal tissue in GTEx Portal database.The corresponding protein interaction network was obtained by STRING.Finally, the candidate genes were searched for their transcription factors in UCSC and JASPAR.Primary human trabecular cells were divided into dexamethasone group and control group, treated with 2 ml 500 nmol/L dexamethasone and the same amount of ethanol, respectively.The expression of BDKRB1 and TAGLN in trabecular cells was detected by Western blot.Results:Differential genes between GSE37474 and GSE124114 datasets enriched in complement and coagulation cascade by GSEA.There were 89 intersecting genes of the three datasets.These genes mainly regulated the formation of extracellular matrix by GO analysis.The gene with the highest enrichment score and collagen-containing extracellular matrix was found to be associated with fibroblasts in GTEx Portal database.ACTA2, MYL9, TAGLN, and LMOD1 were closely related in STRING protein-protein interaction network.Transcription factor SP1 in UCSC and JASPAR according to related genes, BDKRB1, NID1, MFGE8 and TAGLN.The relative expression levels of BDKRB1 and TAGLN proteins were 1.32±0.14 and 0.44±0.09 in dexamethasone group, respectively, which were significantly higher than 1.00±0.00 and 0.20±0.10 in the control group, respectively ( t=-3.61, 2.89; both at P<0.05). Conclusions:Bioinformatics analysis showed that transcription factor SP1 may play a role in human trabecular meshwork cells to myofibroblasts transition after dexamethasone treatment.

Key Prediction Genes of Nasopharyngeal Carcinoma:Screening Based on Systematic Bioinformatics and Validation by Cell Experiments / 中国医学科学院学报

Objective To screen out the potential prediction genes for nasopharyngeal carcinoma(NPC)from the gene microarray data of NPC samples and then verify the genes by cell experiments.Methods The NPC dataset was downloaded from Gene Expression Omnibus,and limma package was employed to screen out the differentially expressed genes.Weighted correlation network analysis package was used for weighted gene co-expression network analysis,and Venn diagram was drawn to find the common genes.The gene ontology annotation and Kyoto encyclopedia of genes and genomes pathway enrichment were then performed for the common genes.The biomarkers for NPC were further explored by protein-protein interaction network,LASSO regression,and non-parametric tests.Real-time quantitative PCR and Western blotting were employed to determine the mRNA and protein levels of key predictors of NPC,so as to verify the screening results.Results There were 622 up-regulated genes and 351 down-regulated genes in the GSE12452 dataset.A total of 116 common genes were obtained by limma analysis and weighted gene co-expression network analysis.The common genes were mainly involved in the biological processes of cell proliferation and regulation and regulation of intercellular adhesion.They were mainly enriched in Rap1,Ras,and tumor necrosis factor signaling pathways.Six key genes were screened out,encoding angiopoietin-2(ANGPT2),dual oxidase 2(DUOX2),coagulation factor Ⅲ(F3),interleukin-15(IL-15),lipocalin-2,and retinoic acid receptor-related orphan receptor B(RORB).Real-time quantitative PCR and Western blotting showed that the NPC cells had up-regulated mRNA and protein levels of ANGPT2 and IL-15 and down-regulated mRNA and protein levels of DUOX2,F3,and RORB,which was consistent with the results predicted by bioinformatics.Conclusion ANGPT2,DUOX2,F3,IL-15 and RORB are potential predictive molecular markers and therapeutic targets for NPC,which may be involved in Rap1,Ras,tumor necrosis factor and other signaling pathways.

Integrated bioinformatics approach reveals methylation-regulated differentially expressed genes in obesity

ABSTRACT Objective: To identify DNA methylation and gene expression profiles involved in obesity by implementing an integrated bioinformatics approach. Materials and methods: Gene expression (GSE94752, GSE55200, and GSE48964) and DNA methylation (GSE67024 and GSE111632) datasets were obtained from the GEO database. Differentially expressed genes (DEGs) and differentially methylated genes (DMGs) in subcutaneous adipose tissue of patients with obesity were identified using GEO2R. Methylation-regulated DEGs (MeDEGs) were identified by overlapping DEGs and DMGs. The protein-protein interaction (PPI) network was constructed with the STRING database and analyzed using Cytoscape. Functional modules and hub-bottleneck genes were identified by using MCODE and CytoHubba plugins. Functional enrichment analyses were performed based on Gene Ontology terms and KEGG pathways. To prioritize and identify candidate genes for obesity, MeDEGs were compared with obesity-related genes available at the DisGeNET database. Results: A total of 54 MeDEGs were identified after overlapping the lists of significant 274 DEGs and 11,556 DMGs. Of these, 25 were hypermethylated-low expression genes and 29 were hypomethylated-high expression genes. The PPI network showed three hub-bottleneck genes (PTGS2, TNFAIP3, and FBXL20) and one functional module. The 54 MeDEGs were mainly involved in the regulation of fibroblast growth factor production, the molecular function of arachidonic acid, and ubiquitin-protein transferase activity. Data collected from DisGeNET showed that 11 of the 54 MeDEGs were involved in obesity. Conclusion: This study identifies new MeDEGs involved in obesity and assessed their related pathways and functions. These results data may provide a deeper understanding of methylation-mediated regulatory mechanisms of obesity.

Expression analysis of transcription factors in sugarcane during cold stress / Análise de expressão de fatores de transcrição em cana-de-açúcar sob estresse por frio

Transcription factors (TF) are a wide class of genes in plants, and these can regulate the expression of other genes in response to various environmental stresses (biotic and abiotic). In the current study, transcription factor activity in sugarcane was examined during cold stress. Initially, RNA transcript reads of two sugarcane cultivars (ROC22 and GT08-1108) under cold stress were downloaded from SRA NCBI database. The reads were aligned into a reference genome and the differential expression analyses were performed with the R/Bioconductor edgeR package. Based on our analyses in the ROC22 cultivar, 963 TF genes were significantly upregulated under cold stress among a total of 5649 upregulated genes, while 293 TF genes were downregulated among a total of 3,289 downregulated genes. In the GT08-1108 cultivar, 974 TF genes were identified among 5,649 upregulated genes and 283 TF genes were found among 3,289 downregulated genes. Most transcription factors were annotated with GO categories related to protein binding, transcription factor binding, DNA-sequence-specific binding, transcription factor complex, transcription factor activity in RNA polymerase II, the activity of nucleic acid binding transcription factor, transcription corepressor activity, sequence-specific regulatory region, the activity of transcription factor of RNA polymerase II, transcription factor cofactor activity, transcription factor activity from plastid promoter, transcription factor activity from RNA polymerase I promoter, polymerase II and RNA polymerase III. The findings of above results will help to identify differentially expressed transcription factors during cold stress. It also provides a comprehensive analysis of the regulation of the transcription activity of many genes. Therefore, this study provides the molecular basis for improving cold tolerance in sugarcane and other economically important grasses.

Fatores de transcrição (FT) são uma ampla classe de genes em plantas e podem regular a expressão de outros genes em resposta a vários estresses ambientais (estresses bióticos e abióticos). No presente estudo, a atividade do fator de transcrição na cana-de-açúcar foi examinada durante o estresse pelo frio. Inicialmente, as leituras de transcrição de RNA de duas cultivares de cana-de-açúcar (ROC22 e GT08-1108) sob estresse frio foram baixadas do banco de dados SRA NCBI. As leituras foram alinhadas em um genoma de referência e as análises de expressão diferencial foram realizadas com o pacote R / Bioconductor edgeR. Com base em nossas análises no cultivar ROC22, 963 genes TF foram significativamente regulados positivamente sob estresse pelo frio entre um total de 5.649 genes regulados positivamente, enquanto 293 genes TF foram regulados negativamente entre um total de 3.289 genes regulados negativamente. No cultivar GT08-1108, 974 genes TF foram identificados entre 5.649 genes regulados positivamente e 283 genes TF foram encontrados entre 3.289 genes regulados negativamente. Os fatores de transcrição, em sua maioria, foram anotados com categorias GO relacionadas à ligação de proteína, ligação de fator de transcrição, ligação específica de sequência de DNA, complexo de fator de transcrição, atividade de fator de transcrição em RNA polimerase II, atividade de fator de transcrição de ligação de ácido nucleico, atividade de corepressor de transcrição, sequência específica da região reguladora, atividade do fator de transcrição da RNA polimerase II, atividade do cofator do fator de transcrição, atividade do fator de transcrição do promotor do plastídio, atividade do fator de transcrição do promotor da RNA polimerase I, polimerase II e RNA polimerase III. As descobertas dos resultados acima ajudarão a identificar fatores de transcrição expressos diferencialmente durante o estresse pelo frio. Ele também fornece uma análise abrangente da regulação da atividade [...].

Expression analysis of transcription factors in sugarcane during cold stress

Abstract Transcription factors (TF) are a wide class of genes in plants, and these can regulate the expression of other genes in response to various environmental stresses (biotic and abiotic). In the current study, transcription factor activity in sugarcane was examined during cold stress. Initially, RNA transcript reads of two sugarcane cultivars (ROC22 and GT08-1108) under cold stress were downloaded from SRA NCBI database. The reads were aligned into a reference genome and the differential expression analyses were performed with the R/Bioconductor edgeR package. Based on our analyses in the ROC22 cultivar, 963 TF genes were significantly upregulated under cold stress among a total of 5649 upregulated genes, while 293 TF genes were downregulated among a total of 3,289 downregulated genes. In the GT08-1108 cultivar, 974 TF genes were identified among 5,649 upregulated genes and 283 TF genes were found among 3,289 downregulated genes. Most transcription factors were annotated with GO categories related to protein binding, transcription factor binding, DNA-sequence-specific binding, transcription factor complex, transcription factor activity in RNA polymerase II, the activity of nucleic acid binding transcription factor, transcription corepressor activity, sequence-specific regulatory region, the activity of transcription factor of RNA polymerase II, transcription factor cofactor activity, transcription factor activity from plastid promoter, transcription factor activity from RNA polymerase I promoter, polymerase II and RNA polymerase III. The findings of above results will help to identify differentially expressed transcription factors during cold stress. It also provides a comprehensive analysis of the regulation of the transcription activity of many genes. Therefore, this study provides the molecular basis for improving cold tolerance in sugarcane and other economically important grasses.

Resumo Fatores de transcrição (FT) são uma ampla classe de genes em plantas e podem regular a expressão de outros genes em resposta a vários estresses ambientais (estresses bióticos e abióticos). No presente estudo, a atividade do fator de transcrição na cana-de-açúcar foi examinada durante o estresse pelo frio. Inicialmente, as leituras de transcrição de RNA de duas cultivares de cana-de-açúcar (ROC22 e GT08-1108) sob estresse frio foram baixadas do banco de dados SRA NCBI. As leituras foram alinhadas em um genoma de referência e as análises de expressão diferencial foram realizadas com o pacote R / Bioconductor edgeR. Com base em nossas análises no cultivar ROC22, 963 genes TF foram significativamente regulados positivamente sob estresse pelo frio entre um total de 5.649 genes regulados positivamente, enquanto 293 genes TF foram regulados negativamente entre um total de 3.289 genes regulados negativamente. No cultivar GT08-1108, 974 genes TF foram identificados entre 5.649 genes regulados positivamente e 283 genes TF foram encontrados entre 3.289 genes regulados negativamente. Os fatores de transcrição, em sua maioria, foram anotados com categorias GO relacionadas à ligação de proteína, ligação de fator de transcrição, ligação específica de sequência de DNA, complexo de fator de transcrição, atividade de fator de transcrição em RNA polimerase II, atividade de fator de transcrição de ligação de ácido nucleico, atividade de corepressor de transcrição, sequência específica da região reguladora, atividade do fator de transcrição da RNA polimerase II, atividade do cofator do fator de transcrição, atividade do fator de transcrição do promotor do plastídio, atividade do fator de transcrição do promotor da RNA polimerase I, polimerase II e RNA polimerase III. As descobertas dos resultados acima ajudarão a identificar fatores de transcrição expressos diferencialmente durante o estresse pelo frio. Ele também fornece uma análise abrangente da regulação da atividade de transcrição de muitos genes. Portanto, este estudo fornece base molecular para melhorar a tolerância ao frio em cana-de-açúcar e outras gramíneas economicamente importantes.

Expression analysis of transcription factors in sugarcane during cold stress / Análise de expressão de fatores de transcrição em cana-de-açúcar sob estresse por frio

Abstract Transcription factors (TF) are a wide class of genes in plants, and these can regulate the expression of other genes in response to various environmental stresses (biotic and abiotic). In the current study, transcription factor activity in sugarcane was examined during cold stress. Initially, RNA transcript reads of two sugarcane cultivars (ROC22 and GT08-1108) under cold stress were downloaded from SRA NCBI database. The reads were aligned into a reference genome and the differential expression analyses were performed with the R/Bioconductor edgeR package. Based on our analyses in the ROC22 cultivar, 963 TF genes were significantly upregulated under cold stress among a total of 5649 upregulated genes, while 293 TF genes were downregulated among a total of 3,289 downregulated genes. In the GT08-1108 cultivar, 974 TF genes were identified among 5,649 upregulated genes and 283 TF genes were found among 3,289 downregulated genes. Most transcription factors were annotated with GO categories related to protein binding, transcription factor binding, DNA-sequence-specific binding, transcription factor complex, transcription factor activity in RNA polymerase II, the activity of nucleic acid binding transcription factor, transcription corepressor activity, sequence-specific regulatory region, the activity of transcription factor of RNA polymerase II, transcription factor cofactor activity, transcription factor activity from plastid promoter, transcription factor activity from RNA polymerase I promoter, polymerase II and RNA polymerase III. The findings of above results will help to identify differentially expressed transcription factors during cold stress. It also provides a comprehensive analysis of the regulation of the transcription activity of many genes. Therefore, this study provides the molecular basis for improving cold tolerance in sugarcane and other economically important grasses.

Resumo Fatores de transcrição (FT) são uma ampla classe de genes em plantas e podem regular a expressão de outros genes em resposta a vários estresses ambientais (estresses bióticos e abióticos). No presente estudo, a atividade do fator de transcrição na cana-de-açúcar foi examinada durante o estresse pelo frio. Inicialmente, as leituras de transcrição de RNA de duas cultivares de cana-de-açúcar (ROC22 e GT08-1108) sob estresse frio foram baixadas do banco de dados SRA NCBI. As leituras foram alinhadas em um genoma de referência e as análises de expressão diferencial foram realizadas com o pacote R / Bioconductor edgeR. Com base em nossas análises no cultivar ROC22, 963 genes TF foram significativamente regulados positivamente sob estresse pelo frio entre um total de 5.649 genes regulados positivamente, enquanto 293 genes TF foram regulados negativamente entre um total de 3.289 genes regulados negativamente. No cultivar GT08-1108, 974 genes TF foram identificados entre 5.649 genes regulados positivamente e 283 genes TF foram encontrados entre 3.289 genes regulados negativamente. Os fatores de transcrição, em sua maioria, foram anotados com categorias GO relacionadas à ligação de proteína, ligação de fator de transcrição, ligação específica de sequência de DNA, complexo de fator de transcrição, atividade de fator de transcrição em RNA polimerase II, atividade de fator de transcrição de ligação de ácido nucleico, atividade de corepressor de transcrição, sequência específica da região reguladora, atividade do fator de transcrição da RNA polimerase II, atividade do cofator do fator de transcrição, atividade do fator de transcrição do promotor do plastídio, atividade do fator de transcrição do promotor da RNA polimerase I, polimerase II e RNA polimerase III. As descobertas dos resultados acima ajudarão a identificar fatores de transcrição expressos diferencialmente durante o estresse pelo frio. Ele também fornece uma análise abrangente da regulação da atividade de transcrição de muitos genes. Portanto, este estudo fornece base molecular para melhorar a tolerância ao frio em cana-de-açúcar e outras gramíneas economicamente importantes.

Analysis of key genes and targeted protection methods affecting the survival of human islets / 器官移植

Objective To identify the key genes and targeted protection methods affecting the survival of human islets. Methods Using bioinformatics method, the gene expression profile (GSE53454) was selected through screening and comparison from Gene Expression Omnibus(GEO) database. GEO2R tool was employed to screen the differentially expressed gene(DEG) between the human islets exposed (exposure group) and non-exposed (non-exposure group) to interleukin (IL)-1β and interferon (IFN)-γ for 24, 48 and 72 h, respectively. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed by DAVID. Protein-protein interaction (PPI) network was constructed by STRING and Cytoscape apps. Results A total of 69 up-regulated DEGs and 2 down-regulated DEGs were identified. GO analysis showed that during the biological process, DEGs were enriched in the aspects of virus defense and inflammatory response. In cellular components, DEGs were significantly enriched in extracellular space, outside plasma membrane and extracellular regions. Regarding molecular functions, DEGs were significantly enriched in chemokine activity and cytokine activity. KEGG analysis revealed that DEGs were mainly enriched in multiple signaling pathways, such as cytokine-cytokine receptor interaction, virus protein-cytokine and cytokine-receptor interaction, etc. Ten key genes (STAT1, CXCL10, IRF1, IL6, CXCL9, CCL5, CXCL11, ISG15, CD274, IFIT3) with high connectivity were selected by STRING analysis, all of which were significantly up-regulated in human islets exposed to IL-1β and IFN-γ. Six genes (STAT1, CXCL10, CXCL9, CXCL11, CCL5, IL6) were screened by KEGG enrichment analysis, mainly in Toll-like receptor signaling pathway. Conclusions STAT1, CXCL10, CXCL9, CXCL11, CCL5 and IL6 are the key genes affecting the survival of human islets, which are mainly enriched in Toll-like receptor signaling pathway and act as important targets for islet protection.

Mechanism of Huashi Baidu Prescription in Treating Influenza Viral Pneumonia Based on Transcriptome / 中国实验方剂学杂志

ObjectiveTo clarify the therapeutic effect of Huashi Baidu prescription on pneumonia in mice caused by influenza A （H1N1） virus and explore its mechanism based on the transcriptome. MethodA mouse influenza viral pneumonia model was built by intranasal infection with influenza A virus， and mice were continuously administered the drug for five days， so as to investigate the general condition， lung index， viral load， pathological morphology of lung tissue， survival time， and prolongation rate of survival time of mice and clarify the therapeutic effect of Huashi Baidu prescription on influenza viral pneumonia. Transcriptome technology was used to detect the differentially expressed genes in the lung tissue of mice in the model group and the normal group， as well as the Huashi Baidu prescription group and the model group， and the potential core target of the Huashi Baidu prescription for the treatment of influenza viral pneumonia was screened. Real-time fluorescence quantitative polymerase chain reaction （Real-time PCR） was used to verify the effect of Huashi Baidu prescription on the mRNA expression level of core target genes. ResultCompared with the normal group， the lung index and viral load in the lung tissue of the model group were significantly increased （P<0.05， P<0.01）. Compared with the model group， the high-dose group of Huashi Baidu prescription significantly prolonged the survival time of mice infected with influenza A virus （P<0.05） and significantly reduced the lung index value of mice （P<0.05） and the viral load of lung tissue. The high-dose， medium-dose， and low-dose groups of Huashi Baidu prescription could significantly reduce lung tissue inflammation， blood stasis， swelling， and other pathological changes in mice （P<0.05， P<0.01）. Transcriptome analysis of lung tissue showed that core genes were mainly enriched in the nuclear transcription factor-κB （NF-κB） signaling pathway， interleukin-17 （IL-17） signaling pathway， cytokine-cytokine receptor interaction， and other pathways after the intervention of Huashi Baidu prescription. TRAF6， NFKBIA， CCL2， CCL7， and CXCL2 were the top five node genes with combined score values. Real-time PCR validation showed that Huashi Baidu prescription significantly downregulated the mRNA expression of key genes TRAF6 and NFKBIA in the NF-κB signaling pathway， as well as chemokines CCL2， CCL7， and CXCL2 （P<0.05， P<0.01）. ConclusionHuashi Baidu prescription has a therapeutic effect on influenza viral pneumonia， possibly by inhibiting the expression of key nodes TRAF6 and NFKBIA in the NF-κB signaling pathway and that of chemokines CCL2， CCL7， and CXCL2， reducing the recruitment of inflammatory cells and viral load， and exerting anti-influenza viral pneumonia effects.

Identification of potential hub genes involved in alcohol use disorder via bioinformatics analysis / 四川精神卫生

BackgroundAlcohol use disorder （AUD） is a type of chronic relapsing brain disorder. Genetic factors play an important role in the pathogenesis of AUD. And screening for molecular markers of AUD is of great significance for further elucidating the pathogenesis of the disease， discovering novel therapeutic targets and preventing relapse. ObjectiveTo explore relevant hub genes and potential signal pathways associated with the development of AUD through bioinformatics analysis， and to provide a new direction for the prevention and treatment of AUD. MethodsThe GSE161986 dataset was acquired from the Gene Expression Omnibus （GEO） database. The limma package in R was utilized to identify differentially expressed genes （DEGs）. Gene set enrichment analysis （GSEA） was carried out using the Database for Annotation， Visualization and Integrated Discovery （DAVID）. A protein–protein interaction （PPI） network of DEGs was assessed using the STRING database and visualized by Cytoscape software. Finally， hub genes were validated in GSE44456 dataset. ResultsA total of 114 DEGs were identified. GSEA revealed that these genes were mainly involved in the regulation of signal transduction， protein binding， membrane trafficking and MAPK signaling pathway. PPI network and validation study indicated that GAD1， TIMP1 and CD44 were potential hub genes involved in AUD. ConclusionAberrant expression of GAD1 and TIMP1 as well as MAPK signaling pathway may play a key role in the pathogenesis of AUD， and may serve as potential molecular targets for the diagnosis and treatment of AUD. ［Funded by "Flying Project" of Shanghai Mental Health Center （number， 2022-FX-01）］

Identification of M1 macrophage-related genes in rejection after kidney transplantation based on weighted gene co-expression network analysis / 器官移植

Objective To identify M1 macrophage-related genes in rejection after kidney transplantation and construct a risk prediction model for renal allograft survival. Methods GSE36059 and GSE21374 datasets after kidney transplantation were downloaded from Gene Expression Omnibus (GEO) database. GSE36059 dataset included the samples from the recipients with rejection and stable allografts. Using this dataset, weighted gene co-expression network analysis (WGCNA) and differential analysis were conducted to screen the M1 macrophage-related differentially expressed gene (M1-DEG). Then, GSE21374 dataset (including the follow-up data of graft loss) was divided into the training set and validation set according to a ratio of 7∶3. In the training set, a multivariate Cox's model was constructed using the variables screened by least absolute shrinkage and selection operator (LASSO), and the ability of this model to predict allograft survival was evaluated. CIBERSORT was employed to analyze the differences of infiltrated immune cells between the high-risk group and low-risk group, and the distribution of human leukocyte antigen (HLA)-related genes was analyzed between two groups. Gene set enrichment analysis (GSEA) was used to further clarify the biological process and pathway enrichment in the high-risk group. Finally, the database was employed to predict the microRNA (miRNA) interacting with the prognostic genes. Results In the GSE36059 dataset, 14 M1-DEG were screened. In the GSE21374 dataset, Toll-like receptor 8 (TLR8), Fc gamma receptor 1B (FCGR1B), BCL2 related protein A1 (BCL2A1), cathepsin S (CTSS), guanylate binding protein 2(GBP2) and caspase recruitment domain family member 16 (CARD16) were screened by LASSO-Cox regression analysis, and a multivariate Cox's model was constructed based on these 6 M1-DEG. The area under curve (AUC) of receiver operating characteristic of this model for predicting the 1- and 3-year graft survival was 0.918 and 0.877 in the training set, and 0.765 and 0.736 in the validation set, respectively. Immune cell infiltration analysis showed that the infiltration of rest and activated CD4+ memory T cells, γδT cells and M1 macrophages were increased in the high-risk group (all P < 0.05). The expression level of HLA I gene was up-regulated in the high-risk group. GSEA analysis suggested that immune response and graft rejection were enriched in the high-risk group. CTSS interacted with 8 miRNA, BCL2A1 and GBP2 interacted with 3 miRNA, and FCGR1B interacted with 1 miRNA. Conclusions The prognostic risk model based on 6 M1-DEG has high performance in predicting graft survival, which may provide evidence for early interventions for high-risk recipients.

Effect of Flavanomarein from Coreopsis tinctoria Ethanol Extract in Small Intestinal Organoids in Insulin-resistant Mice Based on Transcriptome Sequencing / 中国实验方剂学杂志

ObjectiveTo investigate the effects of flavanomarein on the transcriptome of small intestinal organoids in insulin-resistant mice. MethodFirstly， small intestinal organoids of C57BL/6J and db/db mice were established. Ki-67 and E-cadherin expression was determined by immunofluorescence. Small intestinal organoids were divided into the following three groups： C57BL/6J mouse small intestinal organoids as the normal control group， db/db mouse small intestinal organoids as the model group （IR group）， and db/db mouse small intestinal organoids treated with flavanomarein as the administration group （FM group）. Western blot was used to detect the expression of glucagon-like peptide-1（GLP-1） protein on the small intestinal organoids of the three groups. Finally， transcriptome sequencing was performed on samples from the three groups. ResultOn the 6th day of small intestine organoids culture， a cyclic structure was formed around the lumen， and a small intestine organoids culture model was preliminarily established. Immunofluorescence detection showed that ki-67 and E-cadherin were expressed in small intestinal organoids. Western blot results showed that the expression of GLP-1 protein was increased by flavanomarein. In the results of differential expressed gene （DEG） screening， there were 1 862 DEGs in the IR group as compared with the normal control group， and 2 282 DEGs in the FM group as compared with the IR group. Through protein-protein interaction（PPI） network analysis of the DEGs of the two groups， 10 Hub genes， including Nr1i3， Cyp2c44， Ugt2b1， Gsta1， Gstm2， Ptgs1， Gstm4， Cyp2c38， Cyp4a32， and Gpx3， were obtained. These genes were highly expressed in the normal control group， and their expression was reduced in the IR group. After the intervention of flavanomarein， the expression of the above genes was reversed. ConclusionFlavanomarein may play its role in improving insulin resistance by reversing the expression levels of 10 Hub genes， including Nr1i3， Cyp2c44， Ugt2b1， Gsta1， Gstm2， Ptgs1， Gstm4， Cyp2c38， Cyp4a32， and Gpx3.

Mechanism of Huangjing Qianshi Decoction in treatment of prediabetic mice based on transcriptome sequencing / 中国中药杂志

Based on transcriptome sequencing technology, the mouse model of prediabetes treated with Huangjing Qianshi Decoction was sequenced to explore the possible mechanism of treating prediabetes. First of all, transcriptome sequencing was performed on the normal BKS-DB mouse group, the prediabetic model group, and the Huangjing Qianshi Decoction treatment group(treatment group) to obtain differentially expressed genes in the skeletal muscle samples of mice. The serum biochemical indexes were detected in each group to screen out the core genes of Huangjing Qianshi Decoction in prediabetes. Gene Ontology(GO) database and Kyoto Encyclopedia of Genes and Genomes(KEGG) database were used to conduct signaling pathway enrichment analysis of differentially expressed genes, and real-time quantitative polymerase chain reaction(RT-qPCR) was used to verify them. The results showed that the levels of fasting blood glucose(FBG), fasting insulin(FINS), insulin resistance index(HOMA-IR), total cholesterol(TC), triglycerides(TG), and low-density lipoprotein cholesterol(LDL-C) in the mouse model were significantly decreased after treatment with Huangjing Qianshi Decoction. In the results of differential gene screening, there were 1 666 differentially expressed genes in the model group as compared with the normal group, and there were 971 differentially expressed genes in the treatment group as compared with the model group. Among them, interleukin-6(IL-6) and NR3C2 genes, which were closely related to the regulation of insulin resis-tance function, were significantly up-regulated between the model group and the normal group, and vascular endothelial growth factor A(VEGFA) genes were significantly down-regulated between the model group and the normal group. However, the expression results of IL-6, NR3C2, and VEGFA genes were adverse between the treatment group and the model group. GO functional enrichment analysis found that the biological process annotation mainly focused on cell synthesis, cycle, and metabolism; cell component annotation mainly focused on organelles and internal components; and molecular function annotation mainly focused on binding molecular functions. KEGG pathway enrichment analysis found that it involved the protein tyrosine kinase 6(PTK6) pathway, CD28-dependent phosphoinositide 3-kinase/protein kinase B(PI3K/AKT) pathway, p53 pathway, etc. Therefore, Huangjing Qianshi Decoction can improve the state of prediabetes, and the mechanism may be related to cell cycle and apoptosis, PI3K/AKT pathway, p53 pathway, and other biological pathways regulated by IL-6, NR3C2, and VEGFA.

Identification of hub genes in airway epithelial cells of asthma patients by WGCNA and PPI network

Bronchial asthma is a common chronic disease of airway inflammation, high mucus secretion and airway hyper responsiveness. The pathogenetic mechanisms of asthma remain unclear. In this study, we aimed at identifying genes playing an import role in disease-related pathways in airway epithelial cells of asthma patients. Microarray data GSE41861 of asthma airway epithelial cells was used to screen differentially expressed genes (DEGs) through GEO2R analysis. The weighted gene co-expression network analysis (WGCNA) was performed to identify gene co-expression network modules in bronchial asthma. The DAVID database was then used to perform functional and pathway enrichment analysis of these DEGs. In addition, we have conducted protein-protein interaction (PPI) network of DEGs by STRING, and eventually found key genes and significant modules. A total of 315 DEGs (111 up-regulated and 204 down-regulated) were identified between severe asthma and healthy individual, which were mainly involved in pathways of cilium assembly, cilium morphogenesis, axon guidance, positive regulation of fat cell differentiation, and positive regulation of cell substrate adhesion. A total of 60 genes in the black module and green module were considered to be correlated with the severity of asthma. Combining PPI network, several key genes were identified, such as BP2RY14, PTGS1, SLC18A2, SIGLEC6, RGS13, CPA3, and HPGDS. Our findings revealed several genes that may be involved in the process of development of bronchial asthma and potentially be candidate targets for diagnosis or therapy of bronchial asthma.

Prediction of protein-protein interaction networks and druggable genes associated with parkinson’s disease-Supplementary data

Parkinson’s disease (PD) affects about 2-3% of the global population over 65 years of age and hence, it is the second most common neurodegenerative disorder in the world. This study explored the key genes and miRNA involved in PD. Microarray dataset (accession number GSE19587) comprising of two regions of medulla: dorsal motor nucleus of vagus (DMNV) and inferior olivary nucleus (ION) was downloaded from Gene Expression Omnibus (GEO) database. A total of 697 DEGs from ION (605 up-regulated genes and 92 down-regulated genes) and 663 DEGs from DMNV (638 up-regulated genes and 25 down-regulated genes) were screened. These DEGs were found to be enriched in 46 (DMNV) and 24 (ION) pathways common in DAVID and Comparative Toxicogenomics Database. In PPI network analysis, IGF1 and CD44 were identified as hub genes in DMNV whereas for ION, the hub genes identified were CSF2 and CD44. In TF-miRNA-target gene networks, an aggregate of 11 transcription factors and 46 miRNA were observed to influence the target genes. In drug-gene interaction studies, CYP3A5 and ESR1 had higher connective degrees and hence, they might be novel druggable targets for Parkinson’s disease.