Comprehensive molecular characterization and identification of prognostic signature in stomach adenocarcinoma on the basis of energy-metabolism-related genes

BACKGROUND Stomach adenocarcinoma(STAD)is a leading cause of cancer deaths,but its molecular and prognostic characteristics has never been fully illustrated.AIM To describe a molecular evaluation of primary STAD and develop new therapies and identify promising prognostic signatures.METHODS We describe a comprehensive molecular evaluation of primary STAD based on comprehensive analysis of energy-metabolism-related gene(EMRG)expression profiles.RESULTS On the basis of 86 EMRGs that were significantly associated to patients’progression-free survival(PFS),we propose a molecular classification dividing gastric cancer into two subtypes:Cluster 1,most of which are young patients and display more immune and stromal cell components in tumor microenvironment and lower tumor priority;and Cluster 2,which show early stages and better PFS.Moreover,we construct a 6-gene signature that can classify the prognostic risk of patients after a three-phase training test and validation process.Compared with patients with low-risk score,patients with high-risk score had shorter overall survival.Furthermore,calibration and DCA analysis plots indicate the excellent predictive performance of the 6-gene signature,and which present higher robustness and clinical usability compared with three previous reported prognostic gene signatures.According to gene set enrichment analysis,gene sets related to the high-risk group were participated in the ECM receptor interaction and hedgehog signaling pathway.CONCLUSION Identification of the EMRG-based molecular subtypes and prognostic gene model provides a roadmap for patient stratification and trials of targeted therapies.

In Silico analysis and linking of metabolism-related genes with the immune landscape in head and neck squamous cell carcinoma

Metabolic reprogramming and immunologic suppression are two critical characteristics promoting the progression of head and neck squamous cell carcinoma(HNSCC).The integrative analysis of all the metabolismrelated genes(MRGs)in HNSCC is lacking and the interaction between the metabolism and the immune characteristics also requires more exploration to uncover the potential mechanisms.Therefore,this study was designed to establish a prognostic signature based on all the MRGs in HNSCC.Genes of HNSCC samples were available from the TCGA and GEO databases while the MRGs were retrieved from a previous study.Ultimately 4 prognostic MRGs were selected to construct a model possessing robust prognostic value and accuracy in TCGA cohorts.The favorable reproducibility of this model was confirmed in validation cohorts from GEO databases.The risk score calculated by this model was an independent prognostic factor that further classified these HNSCC patients into high-/low-risk groups.GSEA analyses and somatic mutations indicated the low-risk group could activate several anti-tumor pathways and possessed lower TP53 mutation.The results of ESTIMATE,single-sample GSEA,CIBERSORT,and some immune-related molecules analyses suggested the low-risk group exhibited lower metabolic activities and higher immune characteristics.The Spearman correlation test implied most metabolic pathways with tumor-promoting function were negatively correlated with the immune activity,indicating a plausible approach of combining the anti-metabolism and the immunotherapy drugs in the high-risk group to enhance therapeutic effects than applied separately.In conclusion,this prognostic signature linking MRGs with the immune landscape could promote the individualized treatment for HNSCC patients.

Prognostic Biomarkers for Hepatocellular Carcinoma Based on Serine and Glycine Metabolism-related Genes

Background and Aims:Targeted therapy and immunotherapy have emerged as treatment options for hepatocellular carcinoma(HCC)in recent years.The significance of serine and glycine metabolism in various cancers is widely acknowledged.This study aims to investigate their correlation with the prognosis and tumor immune microenvironment(TIME)of HCC.Methods:Based on the public database,different subtypes were identified by cluster analysis,and the prognostic model was constructed through regression analysis.The gene expression omnibus(GEO)data set was used as the validation set to verify the performance of the model.The survival curve evaluated prognostic ability.CIBERSORT was used to evaluate the level of immune cell infiltration,and maftools analyzed the mutations.DsigDB screened small molecule compounds related to prognostic genes.Results:HCC was found to have two distinct subtypes.Subsequently,we constructed a risk score prognostic model through regression analysis based on serine and glycine metabolismrelated genes(SGMGs).A nomogram was constructed based on risk scores and other clinical factors.HCC patients with a higher risk score showed a poor prognosis,and there were significant differences in immune cell infiltration between the high-and low-risk groups.In addition,three potential drugs associated with prognostic genes,streptozocin,norfloxacin,and hydrocotarnine,were identified.Conclusions:This study investigated the expression patterns of SGMGs and their relationship with tumor characteristics,resulting in the development of a novel model for predicting the prognosis of HCC patients.The study provides a reference for clinical prognosis prediction and treatment of HCC patients.

Expression Analysis of Proline Metabolism-related Genes From Halophyte Arabis stelleri under Osmotic Stress Conditions

Arabis stelleri var. japonica evidenced stronger osmotic stress tolerance than Arabidopsis thaliana. Using an A. thaliana microarray chip, we determined changes in the expression of approximately 2 800 genes between A. stelleri plants treated with 0.2 M mannitol versus mock-treated plants. The most significant changes in the gene expression patterns were in genes defining cellular components or in genes associated with the endomembrane system, stimulus response, stress response, chemical stimulus response, and defense response. The expression patterns of three de novo proline biosynthesis enzymes were evaluated in A. stelleri var. japonica seedlings treated with 0.2 M mannitol, 0.2 M sorbitol, and 0.2 M NaCI. The expression of At-pyrroline-5-carboxylate synthetase was not affected by NaCI stress but was similarly induced by mannitol and sorbitol. The proline dehydrogenase gene, which is known to be repressed by dehydration stress and induced by free L-proline, was induced at an early stage by mannitol treatment, but the level of proline dehydrogenase was increased later by treatment with both mannitol and NaCI. The level of free L-proline accumulation increased progressively in response to treatments with mannitol, sorbitol, and NaCI. Mannitol induced L-proline accumulation more rapidly than NaCI or sorbitol. These findings demonstrate that the osmotic tolerance of the novel halophyte, Arabis stelleri, is associated with the accumulation of L-proline.

Dysregulation of genes involved in the long-chain fatty acid transport in pancreatic ductal adenocarcinoma

BACKGROUND Pancreatic ductal adenocarcinoma(PDAC)is an aggressive lethal malignancy with limited options for treatment and a 5-year survival rate of 11%in the United States.As for other types of tumors,such as colorectal cancer,aberrant de novo lipid synthesis and reprogrammed lipid metabolism have been suggested to be associated with PDAC development and progression.AIM To identify the possible involvement of lipid metabolism in PDAC by analyzing in tumoral and non-tumoral tissues the expression level of the most relevant genes involved in the long-chain fatty acid(FA)import into cell.METHODS A gene expression analysis of FASN,CD36,SLC27A1,SLC27A2,SLC27A3,SLC27A4,SLC27A5,ACSL1,and ACSL3 was performed by qRT-PCR in 24 tumoral PDAC tissues and 11 samples from non-tumoral pancreatic tissues obtained via fine needle aspiration or via surgical resection.The genes were considered significantly dysregulated between the groups when the p value was<0.05 and the fold change(FC)was≤0.5 and≥2.RESULTS We found that three FA transporters and two long-chain acyl-CoA synthetases genes were significantly upregulated in the PDAC tissue compared to the non-tumoral tissue:SLC27A2(FC=5.66;P=0.033),SLC27A3(FC=2.68;P=0.040),SLC27A4(FC=3.13;P=0.033),ACSL1(FC=4.10;P<0.001),and ACSL3(FC=2.67;P=0.012).We further investigated any possible association between the levels of the analyzed mRNAs and the specific characteristics of the tumors,including the anatomic location,the lymph node involvement,and the presence of metastasis.A significant difference in the expression of SLC27A3(FC=3.28;P=0.040)was found comparing patients with and without lymph nodes involvement with an overexpression of this transcript in 17 patients presenting tumoral cells in the lymph nodes.CONCLUSION Despite the low number of patients analyzed,these preliminary results seem to be promising.Addressing lipid metabolism through a broad strategy could be a beneficial way to treat this malignancy.Future in vitro and in vivo studies on these genes may offer important insights into the mechanisms linking PDAC with the long-chain FA import pathway.

PHOSPHO1 Serves as a Key Metabolism-Related Biomarker in the Tumorigenesis of Diffuse Large B-cell Lymphoma

Objective:Diffuse large B-cell lymphoma(DLBCL)is an aggressive type of non-Hodgkin lymphoma.Due to its genetic heterogeneity and abnormal metabolism,many DLBCL patients have a poor prognosis.This study investigated the key metabolism-related genes and potential mechanisms.Methods:Differentially expressed genes,differentially expressed transcription factors(TFs),and differentially expressed metabolism-related genes(DEMRGs)of glucose and lipid metabolic processes were identified using the edgeR package.Key DEMRGs were screened by Lasso regression,and a prediction model was constructed.The cell type identification by estimating relative subsets of RNA transcripts algorithm was utilized to assess the fraction of immune cells,and Gene Set Enrichment Analysis was used to determine immune-related pathways.A regulatory network was constructed with significant co-expression interactions among TFs,DEMRGs,immune cells/pathways,and hallmark pathways.Results:A total of 1551 DEMRGs were identified.A prognostic model with a high applicability(area under the curve=0.921)was constructed with 13 DEMRGs.Tumorigenesis of DLBCL was highly related to the neutrophil count.Four DEMRGs(PRXL2AB,CCN1,DECR2 and PHOSPHO1)with 32 TF-DEMRG,36 DEMRG-pathway,14 DEMRG-immune-cell,9 DEMRG-immune-gene-set,and 67 DEMRG-protein-chip interactions were used to construct the regulatory network.Conclusion:We provided a prognostic prediction model based on 13 DEMRGs for DLBCL.We found that phosphatase,orphan 1(PHOSPHO1)is positively regulated by regulatory factor X5(RFX5)and mediates MYC proto-oncogene(MYC)targeting the V2 pathway and neutrophils.

RPLP0/TBP are the most stable reference genes for human dental pulp stem cells under osteogenic differentiation

BACKGROUND Validation of the reference gene(RG)stability during experimental analyses is essential for correct quantitative real-time polymerase chain reaction(RT-qPCR)data normalisation.Commonly,in an unreliable way,several studies use genes involved in essential cellular functions[glyceraldehyde-3-phosphate dehydro-genase(GAPDH),18S rRNA,andβ-actin]without paying attention to whether they are suitable for such experimental conditions or the reason for choosing such genes.Furthermore,such studies use only one gene when Minimum Information for Publication of Quantitative Real-Time PCR Experiments guidelines recom-mend two or more genes.It impacts the credibility of these studies and causes dis-tortions in the gene expression findings.For tissue engineering,the accuracy of gene expression drives the best experimental or therapeutical approaches.We cultivated DPSCs under two conditions:Undifferentiated and osteogenic dif-ferentiation,both for 35 d.We evaluated the gene expression of 10 candidates for RGs[ribosomal protein,large,P0(RPLP0),TATA-binding protein(TBP),GAPDH,actin beta(ACTB),tubulin(TUB),aminolevulinic acid synthase 1(ALAS1),tyro-sine 3-monooxygenase/tryptophan 5-monooxygenase activation protein,zeta(YWHAZ),eukaryotic translational elongation factor 1 alpha(EF1a),succinate dehydrogenase complex,subunit A,flavoprotein(SDHA),and beta-2-micro-globulin(B2M)]every 7 d(1,7,14,21,28,and 35 d)by RT-qPCR.The data were analysed by the four main algorithms,ΔCt method,geNorm,NormFinder,and BestKeeper and ranked by the RefFinder method.We subdivided the samples into eight subgroups.RESULTS All of the data sets from clonogenic and osteogenic samples were analysed using the RefFinder algorithm.The final ranking showed RPLP0/TBP as the two most stable RGs and TUB/B2M as the two least stable RGs.Either theΔCt method or NormFinder analysis showed TBP/RPLP0 as the two most stable genes.However,geNorm analysis showed RPLP0/EF1αin the first place.These algorithms’two least stable RGs were B2M/GAPDH.For BestKeeper,ALAS1 was ranked as the most stable RG,and SDHA as the least stable RG.The pair RPLP0/TBP was detected in most subgroups as the most stable RGs,following the RefFinfer ranking.CONCLUSION For the first time,we show that RPLP0/TBP are the most stable RGs,whereas TUB/B2M are unstable RGs for long-term osteogenic differentiation of human DPSCs in traditional monolayers.

Identification of hub genes associated with Helicobacter pylori infection and type 2 diabetes mellitus:A pilot bioinformatics study

BACKGROUND Helicobacter pylori(H.pylori)infection is related to various extragastric diseases including type 2 diabetes mellitus(T2DM).However,the possible mechanisms connecting H.pylori infection and T2DM remain unknown.AIM To explore potential molecular connections between H.pylori infection and T2DM.METHODS We extracted gene expression arrays from three online datasets(GSE60427,GSE27411 and GSE115601).Differentially expressed genes(DEGs)commonly present in patients with H.pylori infection and T2DM were identified.Hub genes were validated using human gastric biopsy samples.Correlations between hub genes and immune cell infiltration,miRNAs,and transcription factors(TFs)were further analyzed.RESULTS A total of 67 DEGs were commonly presented in patients with H.pylori infection and T2DM.Five significantly upregulated hub genes,including TLR4,ITGAM,C5AR1,FCER1G,and FCGR2A,were finally identified,all of which are closely related to immune cell infiltration.The gene-miRNA analysis detected 13 miRNAs with at least two gene cross-links.TF-gene interaction networks showed that TLR4 was coregulated by 26 TFs,the largest number of TFs among the 5 hub genes.CONCLUSION We identified five hub genes that may have molecular connections between H.pylori infection and T2DM.This study provides new insights into the pathogenesis of H.pylori-induced onset of T2DM.

Mining elite loci and candidate genes for root morphology-related traits at the seedling stage by genome-wide association studies in upland cotton(Gossypium hirsutum L.)

Root system architecture plays an essential role in water and nutrient acquisition in plants,and it is significantly involved in plant adaptations to various environmental stresses.In this study,a panel of 242 cotton accessions was collected to investigate six root morphological traits at the seedling stage,including main root length(MRL),root fresh weight(RFW),total root length(TRL),root surface area(RSA),root volume(RV),and root average diameter(AvgD).The correlation analysis of the six root morphological traits revealed strong positive correlations of TRL with RSA,as well as RV with RSA and AvgD,whereas a significant negative correlation was found between TRL and AvgD.Subsequently,a genome-wide association study(GWAS)was performed using the root phenotypic and genotypic data reported previously for the 242 accessions using 56,010 single nucleotide polymorphisms(SNPs)from the CottonSNP80K array.A total of 41 quantitative trait loci(QTLs)were identified,including nine for MRL,six for RFW,nine for TRL,12 for RSA,12 for RV and two for AvgD.Among them,eight QTLs were repeatedly detected in two or more traits.Integrating these results with a transcriptome analysis,we identified 17 candidate genes with high transcript values of transcripts per million(TPM)≥30 in the roots.Furthermore,we functionally verified the candidate gene GH_D05G2106,which encodes a WPP domain protein 2in root development.A virus-induced gene silencing(VIGS)assay showed that knocking down GH_D05G2106significantly inhibited root development in cotton,indicating its positive role in root system architecture formation.Collectively,these results provide a theoretical basis and candidate genes for future studies on cotton root developmental biology and root-related cotton breeding.

PbrARF4 contributes to calyx shedding of fruitlets in ‘Dangshan Suli’ pear by partly regulating the expression of abscission genes

Fruitlet calyx shedding in pear plants is apparently regulated via numerous pathways that involve both environmental triggers and phytohormones cues such as auxin. In this study, we found at 10 days after full bloom (DAFB) higher levels of indoleacetic acid (IAA) and tryptophan (Trp) in calyx persistence fruitlet (CPF) than calyx shedding fruitlet (CSF) ofDanshan Suli’ pear (Pyrus bretschneideri Rhed.). Consisting with this, the activity of indolealdehyde oxidase (IAAIdO), which promotes IAA synthesis, was remarkably increased, and that of peroxidase(POD), which degrades IAA, dropped markedly in CPF but not in CSF. Further, qRT-PCR results revealed that most of 31 PbrARFs (encoding auxin response factors) in Pyrus bretschneideri were highly expressed in CPF, whereas PbrARF4, PbrARF24 and PbrARF26 were significantly downregulated in CPF vis-a-vis CSF. Phylogenetic analysis revealed that 6 PbrARFs clustered in the group III, where PbrARF4 showed the closest affinity with AtARF1 that promotes organ abscission, indicating a putative role of PbrARF4 in mediating the process of calyx shedding in pear. In fact, the ectopic overexpression of PbrARF4 in Solanum lycopersicum resulted in an earlier-formed and deeper abscission layer (AL) in the transgenic plants, whose calyxes were more prone to wilt at the mature red stage (MR) compared with the control plants (wild-type). More importantly, expression levels of the abscission genes SILS and Sl Cel2 in transgenic plants overexpressing PbrARF4 were significantly upregulated in comparation with the WT, whereas those of Sl BI and Sl TAPG2 were considerably inhibited. Further, PbrJOINTLESS and PbrIDA,the two genes related to calyx shedding in pear, were up-regulated more in CSF than CPF. The findings contribute to a better understanding of PbrARFs involved in fruitlet calyx shedding of pear, which could prove beneficial to improving the quality of pear fruit.

Regulatory potential of soil available carbon,nitrogen,and functional genes on N_(2)O emissions in two upland plantation systems

Dynamic nitrification and denitrification processes are affected by changes in soil redox conditions,and they play a vital role in regulating soil N_(2)O emissions in rice-based cultivation.It is imperative to understand the influences of different upland crop planting systems on soil N_(2)O emissions.In this study,we focused on two representative rotation systems in Central China:rapeseed–rice(RR)and wheat–rice(WR).We examined the biotic and abiotic processes underlying the impacts of these upland plantings on soil N_(2)O emissions.The results revealed that during the rapeseed-cultivated seasons in the RR rotation system,the average N_(2)O emissions were 1.24±0.20 and 0.81±0.11 kg N ha^(–1)for the first and second seasons,respectively.These values were comparable to the N_(2)O emissions observed during the first and second wheat-cultivated seasons in the WR rotation system(0.98±0.25 and 0.70±0.04 kg N ha^(–1),respectively).This suggests that upland cultivation has minimal impacts on soil N_(2)O emissions in the two rotation systems.Strong positive correlations were found between N_(2)O fluxes and soil ammonium(NH_(4)^(+)),nitrate(NO_(3)^(–)),microbial biomass nitrogen(MBN),and the ratio of soil dissolved organic carbon(DOC)to NO_(3)^(–)in both RR and WR rotation systems.Moreover,the presence of the AOA-amoA and nirK genes were positively associated with soil N_(2)O fluxes in the RR and WR systems,respectively.This implies that these genes may have different potential roles in facilitating microbial N_(2)O production in various upland plantation models.By using a structural equation model,we found that soil moisture,mineral N,MBN,and the AOA-amoA gene accounted for over 50%of the effects on N_(2)O emissions in the RR rotation system.In the WR rotation system,soil moisture,mineral N,MBN,and the AOA-amoA and nirK genes had a combined impact of over 70%on N_(2)O emissions.These findings demonstrate the interactive effects of functional genes and soil factors,including soil physical characteristics,available carbon and nitrogen,and their ratio,on soil N_(2)O emissions during upland cultivation seasons under rice-upland rotations.

Knockdown of the atypical protein kinase genes GhABC1K2-A05 and GhABC1K12-A07 make cotton more sensitive to salt and PEG stress

Activity of bc1 complex kinase(ABC1K)is an atypical protein kinase(aPK)that plays a crucial role in plant mitochondrial and plastid stress responses,but little is known about the responses of ABC1Ks to stress in cotton(Gossypium spp.).Here,we identified 40 ABC1Ks in upland cotton(Gossypium hirsutum L.)and found that the Gh ABC1Ks were unevenly distributed across 17 chromosomes.The GhABC1K family members included 35 paralogous gene pairs and were expanded by segmental duplication.The GhABC1K promoter sequences contained diverse cis-acting regulatory elements relevant to hormone or stress responses.The qRT-PCR results revealed that most Gh ABC1Ks were upregulated by exposure to different stresses.Gh ABC1K2-A05 and Gh ABC1K12-A07 expression levels were upregulated by at least three stress treatments.These genes were further functionally characterized by virus-induced gene silencing(VIGS).Compared with the controls,the Gh ABC1K2-A05-and Gh ABC1K12-A07-silenced cotton lines exhibited higher malondialdehyde(MDA)contents,lower catalase(CAT),peroxidase(POD)and superoxide dismutase(SOD)activities and reduced chlorophyll and soluble sugar contents under NaCl and PEG stress.In addition,the expression levels of six stress marker genes(Gh DREB2A,Gh SOS1,Gh CIPK6,Gh SOS2,Gh WRKY33,and Gh RD29A)were significantly downregulated after stress in the Gh ABC1K2-A05-and Gh ABC1K12-A07-silenced lines.The results indicate that knockdown of Gh ABC1K2-A05 and Gh ABC1K12-A07 make cotton more sensitive to salt and PEG stress.These findings can provide valuable information for intensive studies of Gh ABC1Ks in the responses and resistance of cotton to abiotic stresses.

Predictive model using four ferroptosis-related genes accurately predicts gastric cancer prognosis

BACKGROUND Gastric cancer(GC)is a common malignancy of the digestive system.According to global 2018 cancer data,GC has the fifth-highest incidence and the thirdhighest fatality rate among malignant tumors.More than 60%of GC are linked to infection with Helicobacter pylori(H.pylori),a gram-negative,active,microaerophilic,and helical bacterium.This parasite induces GC by producing toxic factors,such as cytotoxin-related gene A,vacuolar cytotoxin A,and outer membrane proteins.Ferroptosis,or iron-dependent programmed cell death,has been linked to GC,although there has been little research on the link between H.pylori infection-related GC and ferroptosis.AIM To identify coregulated differentially expressed genes among ferroptosis-related genes(FRGs)in GC patients and develop a ferroptosis-related prognostic model with discrimination ability.METHODS Gene expression profiles of GC patients and those with H.pylori-associated GC were obtained from The Cancer Genome Atlas and Gene Expression Omnibus(GEO)databases.The FRGs were acquired from the FerrDb database.A ferroptosis-related gene prognostic index(FRGPI)was created using least absolute shrinkage and selection operator–Cox regression.The predictive ability of the FRGPI was validated in the GEO cohort.Finally,we verified the expression of the hub genes and the activity of the ferroptosis inducer FIN56 in GC cell lines and tissues.RESULTS Four hub genes were identified(NOX4,MTCH1,GABARAPL2,and SLC2A3)and shown to accurately predict GC and H.pylori-associated GC.The FRGPI based on the hub genes could independently predict GC patient survival;GC patients in the high-risk group had considerably worse overall survival than did those in the low-risk group.The FRGPI was a significant predictor of GC prognosis and was strongly correlated with disease progression.Moreover,the gene expression levels of common immune checkpoint proteins dramatically increased in the highrisk subgroup of the FRGPI cohort.The hub genes were also confirmed to be highly overexpressed in GC cell lines and tissues and were found to be primarily localized at the cell membrane.The ferroptosis inducer FIN56 inhibited GC cell proliferation in a dose-dependent manner.CONCLUSION In this study,we developed a predictive model based on four FRGs that can accurately predict the prognosis of GC patients and the efficacy of immunotherapy in this population.

Key genes and regulatory networks for diabetic retinopathy based on hypoxia-related genes:a bioinformatics analysis

AIM:To prevent neovascularization in diabetic retinopathy(DR)patients and partially control disease progression.METHODS:Hypoxia-related differentially expressed genes(DEGs)were identified from the GSE60436 and GSE102485 datasets,followed by gene ontology(GO)functional annotation and Kyoto Encyclopedia of Genes and Genomes(KEGG)pathway enrichment analysis.Potential candidate drugs were screened using the CMap database.Subsequently,a protein-protein interaction(PPI)network was constructed to identify hypoxia-related hub genes.A nomogram was generated using the rms R package,and the correlation of hub genes was analyzed using the Hmisc R package.The clinical significance of hub genes was validated by comparing their expression levels between disease and normal groups and constructing receiver operating characteristic curve(ROC)curves.Finally,a hypoxia-related miRNA-transcription factor(TF)-Hub gene network was constructed using the NetworkAnalyst online tool.RESULTS:Totally 48 hypoxia-related DEGs and screened 10 potential candidate drugs with interaction relationships to upregulated hypoxia-related genes were identified,such as ruxolitinib,meprylcaine,and deferiprone.In addition,8 hub genes were also identified:glycogen phosphorylase muscle associated(PYGM),glyceraldehyde-3-phosphate dehydrogenase spermatogenic(GAPDHS),enolase 3(ENO3),aldolase fructose-bisphosphate C(ALDOC),phosphoglucomutase 2(PGM2),enolase 2(ENO2),phosphoglycerate mutase 2(PGAM2),and 6-phosphofructo-2-kinase/fructose-2,6-biphosphatase 3(PFKFB3).Based on hub gene predictions,the miRNA-TF-Hub gene network revealed complex interactions between 163 miRNAs,77 TFs,and hub genes.The results of ROC showed that the except for GAPDHS,the area under curve(AUC)values of the other 7 hub genes were greater than 0.758,indicating their favorable diagnostic performance.CONCLUSION:PYGM,GAPDHS,ENO3,ALDOC,PGM2,ENO2,PGAM2,and PFKFB3 are hub genes in DR,and hypoxia-related hub genes exhibited favorable diagnostic performance.

GST family genes in jujube actively respond to phytoplasma infection

Jujube witches’broom(JWB)caused by phytoplasma has a severely negative effect on multiple metabolisms in jujube.The GST gene family in plants participates in the regulation of a variety of biotic and abiotic stresses.This study aims to identify and reveal the changes in the jujube GST gene family in response to phytoplasma infection.Here,70 ZjGSTs were identified in the jujube genome and divided into 8 classes.Among them,the Tau-class,including 44 genes,was the largest.Phylogenetic analysis indicated that Tau-class genes were highly conserved among species,such as Arabidopsis,cotton,chickpea,and rice.Through chromosome location analysis,37.1%of genes were clustered,and 8 of 9 gene clusters were composed of Tau class members.Through RT-PCR,qRT-PCR and enzyme activity detection,the results showed that the expression of half(20/40)of the tested ZjGSTs was inhibited by phytoplasma infection in field and tissue culture conditions,and GST activity was also significantly reduced.In the resistant and susceptible varieties under phytoplasma infection,ZjGSTU49-ZjGSTU54 in the cluster IV showed opposite expression patterns,which may be due to functional divergence during evolution.Some upregulated genes(ZjGSTU45,ZjGSTU49,ZjGSTU59,and ZjGSTU70)might be involved in the process of jujube against JWB.The yeast two-hybrid results showed that all 6 Tauclass proteins tested could form homodimers or heterodimers.Overall,the comprehensive analysis of the jujube GST gene family revealed that ZjGSTs responded actively to phytoplasma infection.Furthermore,some screened genes(ZjGSTU24,ZjGSTU49-52,ZjGSTU70,and ZjDHAR10)will contribute to further functional studies of jujube-phytoplasma interactions.

Key Genes Involved in the Beneficial Mechanism of Hyperbaric Oxygen for Glioblastoma and Predictive Indicators of Hyperbaric Oxygen Prolonging Survival in Glioblastoma Patients

Objective The prognosis of glioblastoma is poor,and therapy-resistance is largely attributed to intratumor hypoxia.Hyperbaric oxygen(HBO)effectively alleviates hypoxia.However,the sole role of HBO in glioblastoma remains controversial.We previously reported that HBO can promote apoptosis,shorten protrusions,and delay growth of glioblastoma,but the molecular mechanism is unclear.We aimed to test candidate genes in HBO-exposed glioblastoma cells and to analyze their correlation with the survival of glioblastoma patients.Methods Glioblastoma cell lines exposed to repetitive HBO or normobaric air(NBA)were collected for RNA isolation and microarray data analysis.GO analysis,KEGG pathway analysis and survival analysis of the differentially expressed genes(DEGs)were performed.Results HBO not only inhibited hypoxia-inducing genes including CA9,FGF11,PPFIA4,TCAF2 and SLC2A12,but also regulated vascularization by downregulating the expression of COL1A1,COL8A1,COL12A1,RHOJ and FILIP1L,ultimately attenuated hypoxic microenvironment of glioblastoma.HBO attenuated inflammatory microenvironment by reducing the expression of NLRP2,CARD8,MYD88 and CD180.HBO prevented metastasis by downregulating the expression of NTM,CXCL12,CXCL13,CXCR4,CXCR5,CDC42,IGFBP3,IGFBP5,GPC6,MMP19,ADAMTS1,EFEMP1,PTBP3,NF1 and PDCD1.HBO upregulated the expression of BAK1,PPIF,DDIT3,TP53I11 and FAS,whereas downregulated the expression of MDM4 and SIVA1,thus promoting apoptosis.HBO upregulated the expression of CDC25A,MCM2,PCNA,RFC33,DSCC1 and CDC14A,whereas downregulated the expression of ASNS,CDK6,CDKN1B,PTBP3 and MAD2L1,thus inhibiting cell cycle progression.Among these DEGs,17 indicator-genes of HBO prolonging survival were detected.Conclusions HBO is beneficial for glioblastoma.Glioblastoma patients with these predictive indicators may prolong survival with HBO therapy.These potential therapeutic targets especially COL1A1,ADAMTS1 and PTBP3 deserve further validation.

Influence of N, P, Fe Nutrients Availability on Nitrogen Metabolism-Relevant Genes Expression in Skeletonema marinoi

Nitrogen(N), Phosphorus(P), and Iron(Fe) are essential elements for cellular structure and metabolism. In addition to dissolved inorganic nitrogen(DIN), phytoplankton is able to utilize dissolved organic nitrogen(DON). There is general consensus that both bacteria and higher plants nitrogen metabolism is affected by phosphate availability; this was also found to be true in coccolithophorid. Iron affects the structure and function of ecosystems through its effects on nitrogen metabolism. However, it is unclear how these nutrients affect Skeletonema marinoi's nitrogen metabolism. Here, using RT-qPCR, we investigate effects of N, P, and Fe on S. marinoi's nitrogen metabolism and nitrate reductase activity. These results illuminate that in S. marinoi, various nutrients have direct regulation on these genes expression at the molecular level. The varying degree of responses for these genes expression with differing N sources may act to increase the efficiency of nutrient capture when nitrate is limited. Suitable gene expression occurs at a N/P ratio of 16, which represents the atomic N/P ratio of phytoplankton cells and N/P concentrations in ocean; thus, nitrogen metabolism gene expression should be regulated by the existing N/P ratios in the phytoplankton's internal and external environment. Fe concentration has a direct and significant effect on nitrogen metabolism by regulating gene expression and nitrate reductase activity. Gene expression profiles identified in S. marinoi provide a foundation for understanding molecular mechanisms behind diatom nitrogen metabolism with changing N, P, and Fe nutrients allowing a basic understanding of how diatom growth is affected by nutrient utilization.

Fate and Behavior of Tetracycline Resistance Genes in Activated Carbon Adsorption

The accessibility of tetracycline resistance gene (tetG) into the pores of activated carbon (AC), as well as the impact of the pore size distribution (PSD) of AC on the uptake capacity of tetG, were investigated using eight types of AC (four coal-based and four wood-based). AC showed the capability to admit tetG and the average reduction of tetG for coal-based and wood-based ACs at the AC dose of 1 g·L-1 was 3.12 log and 3.65 log, respectively. The uptake kinetic analysis showed that the uptake of the gene followed the pseudo-second-order kinetics reaction, and the uptake rate constant for the coal-based and wood-based ACs was in the range of 5.97 × 10-12 - 4.64 × 10-9 and 7.02 × 10-11 - 1.59 × 10-8 copies·mg-1·min-1, respectively. The uptake capacity analysis by fitting the obtained experiment data with the Freundlich isotherm model indicated that the uptake constant (KF) values were 1.71 × 103 - 8.00 × 109 (copies·g-1)1-1/n for coal-based ACs and 7.00 × 108 - 3.00 × 1010 (copies·g-1)1-1/n for wood-based ones. In addition, the correlation analysis between KF values and pore volume as well as pore surface at different pore size regions of ACs showed that relatively higher positive correlation was found for pores of 50 - 100 Å, suggesting ACs with more pores in this size region can uptake more tetG. The findings of this study are valuable as reference for optimizing the adsorption process regarding antibiotic resistance-related concerns in drinking water treatment.

Comprehensive analysis of the potential pathogenesis of COVID-19 infection and liver cancer

BACKGROUND A growing number of clinical examples suggest that coronavirus disease 2019(COVID-19)appears to have an impact on the treatment of patients with liver cancer compared to the normal population,and the prevalence of COVID-19 is significantly higher in patients with liver cancer.However,this mechanism of action has not been clarified.Gene sets for COVID-19(GSE180226)and liver cancer(GSE87630)were obtained from the Gene Expression Omnibus database.After identifying the common differentially expressed genes(DEGs)of COVID-19 and liver cancer,functional enrichment analysis,protein-protein interaction network construction and scree-ning and analysis of hub genes were performed.Subsequently,the validation of the differential expression of hub genes in the disease was performed and the regulatory network of transcription factors and hub genes was constructed.RESULTS Of 518 common DEGs were obtained by screening for functional analysis.Fifteen hub genes including aurora kinase B,cyclin B2,cell division cycle 20,cell division cycle associated 8,nucleolar and spindle associated protein 1,etc.,were further identified from DEGs using the“cytoHubba”plugin.Functional enrichment analysis of hub genes showed that these hub genes are associated with P53 signalling pathway regulation,cell cycle and other functions,and they may serve as potential molecular markers for COVID-19 and liver cancer.Finally,we selected 10 of the hub genes for in vitro expression validation in liver cancer cells.CONCLUSION Our study reveals a common pathogenesis of liver cancer and COVID-19.These common pathways and key genes may provide new ideas for further mechanistic studies.

Effects of Heterologously Overexpressing PIP5K-Family Genes in Arabidopsis on Inflorescence Development

Castor is one of the top 10 oil crops in the world and has extremely valuable uses.Castor inflorescences directly affect yield,so the study of inflorescence development is very important in increasing castor yield.Our previous studies have shown that the PIP5K gene family(PIP5Ks)is associated with inflorescence development.In this study,to determine the function of each PIP5K gene in castor,a female Lm-type castor line,aLmAB2,was used to determine the relative expression levels of the PIP5Ks in castor inflorescences.Six PIP5K genes were heterologously overexpressed in Arabidopsis thaliana,the relative expression of each gene and the effect on plants was determined in A.thaliana,and the relationships among the PIP5Ks in castor were inferred.The expression levels of the PIP5Ks in the female Lm-type castor line aLmAB2 were analyzed.The relative expression levels of the PIP5K9 and PIP5K11 genes were high(p<0.05)in isofemale inflorescences,and those of PIP5K1,PIP5K2,PIP5K6,and PIP5K8 were high(p<0.05)in female inflorescences but low(p<0.05)in bisexual inflorescences.The PIP5Ks were heterologously overexpressed in A.thaliana,and T3-generation plants with stable genetic resistance,i.e.,AT-PIP5K^(+)plants(AT-PIP5K1^(+),AT-PIP5K2^(+),AT-PIP5K6^(+),AT-PIP5K8^(+),AT-PIP5K9^(+),and ATPIP5K11^(+) plants),were obtained.Biological tests of the AT-PIP5K+plants showed that the growth of the main stem was significantly delayed in AT-PIP5K+plants compared with Columbia wild-type(WT)A.thaliana plants;the PIP5K1 and PIP5K2 genes promoted lateral stem growth and flower and silique development;and the PIP5K6,PIP5K8,PIP5K9 and PIP5K11 genes inhibited lateral stem growth and flower and silique development.The correlations among PIP5Ks in castor suggest that there may be a synergistic relationship among PIP5K1,PIP5K2,and PIP5K6 in castor inflorescences,and PIP5K8,PIP5K9,and PIP5K11 are complementary to the other three genes.