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.

Screening biomarkers for spinal cord injury using weighted gene co-expression network analysis and machine learning

Immune changes and inflammatory responses have been identified as central events in the pathological process of spinal co rd injury.They can greatly affect nerve regeneration and functional recovery.However,there is still limited understanding of the peripheral immune inflammato ry response in spinal cord inju ry.In this study.we obtained microRNA expression profiles from the peripheral blood of patients with spinal co rd injury using high-throughput sequencing.We also obtained the mRNA expression profile of spinal cord injury patients from the Gene Expression Omnibus(GEO)database(GSE151371).We identified 54 differentially expressed microRNAs and 1656 diffe rentially expressed genes using bioinformatics approaches.Functional enrichment analysis revealed that various common immune and inflammation-related signaling pathways,such as neutrophil extracellular trap formation pathway,T cell receptor signaling pathway,and nuclear factor-κB signal pathway,we re abnormally activated or inhibited in spinal cord inju ry patient samples.We applied an integrated strategy that combines weighted gene co-expression network analysis,LASSO logistic regression,and SVM-RFE algorithm and identified three biomarke rs associated with spinal cord injury:ANO10,BST1,and ZFP36L2.We verified the expression levels and diagnostic perfo rmance of these three genes in the original training dataset and clinical samples through the receiver operating characteristic curve.Quantitative polymerase chain reaction results showed that ANO20 and BST1 mRNA levels were increased and ZFP36L2 mRNA was decreased in the peripheral blood of spinal cord injury patients.We also constructed a small RNA-mRNA interaction network using Cytoscape.Additionally,we evaluated the proportion of 22 types of immune cells in the peripheral blood of spinal co rd injury patients using the CIBERSORT tool.The proportions of naive B cells,plasma cells,monocytes,and neutrophils were increased while the proportions of memory B cells,CD8^(+)T cells,resting natural killer cells,resting dendritic cells,and eosinophils were markedly decreased in spinal cord injury patients increased compared with healthy subjects,and ANO10,BST1 and ZFP26L2we re closely related to the proportion of certain immune cell types.The findings from this study provide new directions for the development of treatment strategies related to immune inflammation in spinal co rd inju ry and suggest that ANO10,BST2,and ZFP36L2 are potential biomarkers for spinal cord injury.The study was registe red in the Chinese Clinical Trial Registry(registration No.ChiCTR2200066985,December 12,2022).

Identify the signature genes for diagnose of uveal melanoma by weight gene co-expression network analysis

AIM: To identify and understand the relationship between co-expression pattern and clinic traits in uveal melanoma, weighted gene co-expression network analysis(WGCNA) is applied to investigate the gene expression levels and patient clinic features. Uveal melanoma is the most common primary eye tumor in adults. Although many studies have identified some important genes and pathways that were relevant to progress of uveal melanoma, the relationship between co-expression and clinic traits in systems level of uveal melanoma is unclear yet. We employ WGCNA to investigate the relationship underlying molecular and phenotype in this study.METHODS: Gene expression profile of uveal melanoma and patient clinic traits were collected from the Gene Expression Omnibus(GEO) database. The gene co-expression is calculated by WGCNA that is the R package software. The package is used to analyze the correlation between pairs of expression levels of genes.The function of the genes were annotated by gene ontology(GO).RESULTS: In this study, we identified four co-expression modules significantly correlated with clinictraits. Module blue positively correlated with radiotherapy treatment. Module purple positively correlates with tumor location(sclera) and negatively correlates with patient age. Module red positively correlates with sclera and negatively correlates with thickness of tumor. Module black positively correlates with the largest tumor diameter(LTD). Additionally, we identified the hug gene(top connectivity with other genes) in each module. The hub gene RPS15 A, PTGDS, CD53 and MSI2 might play a vital role in progress of uveal melanoma.CONCLUSION: From WGCNA analysis and hub gene calculation, we identified RPS15 A, PTGDS, CD53 and MSI2 might be target or diagnosis for uveal melanoma.

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 key genes involved in axon regeneration and Wallerian degeneration by weighted gene co-expression network analysis

Peripheral nerve injury repair requires a certain degree of cooperation between axon regeneration and Wallerian degeneration.Therefore,investigating how axon regeneration and degeneration work together to repair peripheral nerve injury may uncover the molecular mechanisms and signal cascades underlying peripheral nerve repair and provide potential strategies for improving the low axon regeneration capacity of the central nervous system.In this study,we applied weighted gene co-expression network analysis to identify differentially expressed genes in proximal and distal sciatic nerve segments from rats with sciatic nerve injury.We identified 31 and 15 co-expression modules from the proximal and distal sciatic nerve segments,respectively.Functional enrichment analysis revealed that the differentially expressed genes in proximal modules promoted regeneration,while the differentially expressed genes in distal modules promoted neurodegeneration.Next,we constructed hub gene networks for selected modules and identified a key hub gene,Kif22,which was up-regulated in both nerve segments.In vitro experiments confirmed that Kif22 knockdown inhibited proliferation and migration of Schwann cells by modulating the activity of the extracellular signal-regulated kinase signaling pathway.Collectively,our findings provide a comparative framework of gene modules that are co-expressed in injured proximal and distal sciatic nerve segments,and identify Kif22 as a potential therapeutic target for promoting peripheral nerve injury repair via Schwann cell proliferation and migration.All animal experiments were approved by the Institutional Animal Ethics Committee of Nantong University,China(approval No.S20210322-008)on March 22,2021.

Identification of key genes regulating the synthesis of quercetin derivatives in Rosa roxburghii through integrated transcriptomics and metabolomics

Rosa roxburghii fruit is rich in flavonoids, but little is known about their biosynthetic pathways. In this study, we employed transcriptomics and metabolomics to study changes related to the flavonoids at five different stages of R. roxburghii fruit development. Flavonoids and the genes related to their biosynthesis were found to undergo significant changes in abundance across different developmental stages, and numerous quercetin derivatives were identified. We found three gene expression modules that were significantly associated with the abundances of the different flavonoids in R. roxburghii and identified three structural UDP-glycosyltransferase genes directly involved in the synthesis of quercetin derivatives within these modules. In addition, we found that RrBEH4, RrLBD1 and RrPIF8could significantly increase the expression of downstream quercetin derivative biosynthesis genes. Taken together,these results provide new insights into the metabolism of flavonoids and the accumulation of quercetin derivatives in R. roxburghii.

Identification of M2 macrophage-related genes for establishing a prognostic model in pancreatic cancer: FCGR3A as key gene

Background:Pancreatic ductal adenocarcinoma(PDAC)has a rich and complex tumor immune microenvironment(TIME).M2 macrophages are among the most extensively infiltrated immune cells in the TIME and are necessary for the growth and migration of cancers.However,the mechanisms and targets mediating M2 macrophage infiltration in pancreatic cancer remain elusive.Methods:The M2 macrophage infiltration score of patients was assessed using the xCell algorithm.Using weighted gene co-expression network analysis(WGCNA),module genes associated with M2 macrophages were identified,and a predictive model was designed.The variations in immunological cell patterns,cancer mutations,and enrichment pathways between the cohorts with the high-and low-risk were examined.Additionally,the expression of FCGR3A and RNASE2,as well as their association with M2 macrophages were evaluated using the HPA,TNMplot,and GEPIA2 databases and verified by tissue immunofluorescence staining.Moreover,in vitro cell experiments were conducted,where FCGR3A was knocked down in pancreatic cancer cells using siRNA to analyze its effects on M2 macrophage infiltration,tumor proliferation,and metastasis.Results:The prognosis of patients in high-risk and low-risk groups was successfully distinguished using a prognostic risk score model of M2 macrophage-related genes(p=0.024).Between the high-and low-risk cohorts,there have been notable variations in immune cell infiltration patterns,tumor mutations,and biological functions.The risk score was linked to the manifestation of prevalent immunological checkpoints,immunological scores,and stroma values(all p<0.05).In vitro experiments and tissue immunofluorescence staining revealed that FCGR3A can promote the infiltration or polarization of M2 macrophages and enhance tumor proliferation and migration.Conclusions:In this study,an M2 macrophage-related pancreatic cancer risk score model was established,and found that FCGR3A was correlated with tumor formation,metastasis,and M2 macrophage infiltration.

Functional investigation and two-sample Mendelian randomization study of primary biliary cholangitis hub genes

BACKGROUND The identification of specific gene expression patterns is crucial for understanding the mechanisms underlying primary biliary cholangitis(PBC)and finding relevant biomarkers for diagnosis and therapeutic evaluation.AIM To determine PBC-associated hub genes and assess their clinical utility for disease prediction.METHODS PBC expression data were obtained from the Gene Expression Omnibus database.Overlapping genes from differential expression analysis and weighted gene coexpression network analysis(WGCNA)were identified as key genes for PBC.Kyoto Encyclopedia of Genes and Genomes and Gene Ontology analyses were performed to explore the potential roles of key genes.Hub genes were identified in protein-protein interaction(PPI)networks using the Degree algorithm in Cytoscape software.The relationship between hub genes and immune cells was investigated.Finally,a Mendelian randomization study was conducted to determine the causal effects of hub genes on PBC.RESULTS We identified 71 overlapping key genes using differential expression analysis and WGCNA.These genes were primarily enriched in pathways related to cytokinecytokine receptor interaction,and Th1,Th2,and Th17 cell differentiation.We utilized Cytoscape software and identified five hub genes(CD247,IL10,CCL5,CCL3,and STAT3)in PPI networks.These hub genes showed a strong correlation with immune cell infiltration in PBC.However,inverse variance weighting analysis did not indicate the causal effects of hub genes on PBC risk.CONCLUSION Hub genes can potentially serve as valuable biomarkers for PBC prediction and treatment,thereby offering significant clinical utility.

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.

Identification of the key genes and mechanisms associated with transcatheter arterial chemoembolisation refractoriness in hepatocellular carcinoma

BACKGROUND Transcatheter arterial embolisation(TACE)is the primary treatment for intermediate-stage hepatocellular carcinoma(HCC)patients while some HCC cases have shown resistance to TACE.AIM To investigate the key genes and potential mechanisms correlated with TACE refractoriness in HCC.METHODS The microarray datasets of TACE-treated HCC tissues,HCC and non-HCC tissues were collected by searching multiple public databases.The respective differentially expressed genes(DEGs)were attained via limma R package.Weighted gene co-expression network analysis was employed for identifying the significant modules related to TACE non-response.TACE refractoriness-related genes were obtained by intersecting up-regulated TACE-associated and HCC-associated DEGs together with the genes in significant modules related to TACE nonresponse.The key genes expression in the above two pairs of samples was compared respectively via Wilcoxon tests and standard mean differences model.The prognostic value of the key genes was evaluated by Kaplan-Meier curve.Multivariate analysis was utilised to investigate the independent prognostic factor in key genes.Single-cell RNA(scRNA)sequencing analysis was conducted to explore the cell types in HCC.TACE refractoriness-related genes activity was calculated via AUCell packages.The CellChat R package was used for the investigation of the cell–cell communication between the identified cell types.RESULTS HCC tissues of TACE non-responders(n=66)and TACE responders(n=81),HCC(n=3941)and non-HCC(n=3443)tissues were obtained.The five key genes,DLG associated protein 5(DLGAP5),Kinesin family member 20A(KIF20A),Assembly factor for spindle microtubules(ASPM),Kinesin family member 11(KIF11)and TPX2 microtubule nucleation factor(TPX2)in TACE refractoriness-related genes,were identified.The five key genes were all up-regulated in the TACE non-responders group and the HCC group.High expression of the five key genes predicted poor prognosis in HCC.Among the key genes,TPX2 was an independent prognostic factor.Four cell types,hepatocytes,embryonic stem cells,T cells and B cells,were identified in the HCC tissues.The TACE refractoriness-related genes expressed primarily in hepatocytes and embryonic stem cells.Hepatocytes,as the providers of ligands,had the strongest interaction with embryonic stem cells that provided receptors.CONCLUSION Five key genes(DLGAP5,KIF20A,ASPM,KIF11 and TPX2)were identified as promoting refractory TACE.Hepatocytes and embryonic stem cells were likely to boost TACE refractoriness.

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.

Identification of Differentially Expressed Genes in Grape Skin at Veraison and Maturity and Construction of Co-expression Network

The ripening process of grape is an important stage during grape growth and development. During this process, color of grape skin is the most obvious change. The molecular mechanism for the ripening of grape(a non-climacteric fruit, which ripens without ethylene and respiration bursts) is still unclear. Although numerous studies have been done on the changes in the contents of metabolites during grape ripening, the differentially expressed genes at veraison and maturity stages have not been systematically analyzed. In this study, 1 524 genes that are significantly differentially expressed in grape(Pinot Noir) skin at veraison and maturity stages were identified, and a co-expression network of these genes was built. Some of the eight co-expression modules we identified may be closely related to the synthesis or metabolism of anthocyanins, sugar acids, and other flavor substances. The transcription factor families WRKY, b ZIP, HSF and WOX may play an important role in the regulation of anthocyanin synthesis or metabolism. The results provide a foundation for further study of the molecular mechanism of grape ripening.

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.

Identification and characterization of genes related to m^(6)A modification in kiwifruit using RNA-seq and ATAC-seq

N6-methyladenosine(m^(6)A)RNA modification is a conserved mechanism that regulates the fate of RNA across eukaryotic organisms.Despite its significance,a comprehensive analysis of m^(6)A-related genes in non-model plants,such as kiwifruit,is lacking.Here,we identified 36 m^(6)A-related genes in the kiwifruit genome according to homology and phylogenetic inference.We performed bioinformatics and evolutionary analyses of the writer,eraser,and reader families of m^(6)A modification.Reanalysis of public RNA-seq data collected from samples under various biotic and abiotic stresses indicated that most m^(6)A-related genes were remarkably expressed under different conditions.Through construction of gene co-expression networks,we found significant correlations between several m^(6)A-related genes and transcription factors(TFs)as well as receptor-like genes during the development and ripening of kiwifruit.Furthermore,we performed ATAC-seq assays on diverse kiwifruit tissues to investigate the regulatory mechanisms of m^(6)A-related genes.We identified 10 common open chromatin regions that were present in at least two tissues,and these regions might serve as potential binding sites for MADS protein,C2H2 protein,and other predicted TFs.Our study offers comprehensive insights into the gene family of m^(6)A-related components in kiwifruit,which will lay foundation for exploring mechanisms of post-transcriptional regulation involved in development and adaptation of kiwifruit.

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.

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.