Functional Annotation Of Genes Research Articles

Alternative splicing in Atlantic salmon head kidney and SHK-1 cell line during the Piscirickettsia salmonis infection: A comparative transcriptome survey

Piscirickettsia salmonis, an intracellular bacterium in salmon aquaculture, is a big challenge because it is responsible for 54.2% of Atlantic salmon mortalities. In recent years, the high relevance of Alternative Splicing (AS) as a molecular mechanism associated with infectious conditions and host-pathogen interaction processes, especially in host immune activation, has been observed. Several studies have highlighted the role of AS in the host's immune response during viral, bacterial, and endoparasite infection. In the present study, we evaluated AS transcriptome profiles during P. salmonis infection in the two most used study models, SHK-1 cell line and salmon head kidney tissue. First, the SHK-1 cell line was exposed to P. salmonis infection at 0-, 7-, and 14-days post-infection (dpi). Following, total RNA was extracted for Illumina sequencing.On the other hand, RNA-Seq datasets of Atlantic salmon head kidney infected with the same P. salmonis strayingwase used. For both study models, the highest number of differentially alternative splicing (DAS) events was observed at 7 dpi, 16,830 DAS events derived from 9213 DAS genes in SHK-1 cells, and 13,820 DAS events from 7684 DAS genes in salmon HK. Alternative first exon (AF) was the most abundant AS type in the three infection times analyzed, representing 31% in SHK-1 cells and 228.6 in salmon HK; meanwhile, mutually exclusive exon (MX) was the least abundant. Notably, functional annotation of DAS genes in SHK-1 cells infected with P. salmonis showed a high presence of genes related to nucleotide metabolism. In contrast, the salmon head kidney exhibited many GO terms associated with immune response. Our findings reported the role of AS during P. salmonis infection in Atlantic salmon. These studies would contribute to a better understanding of the molecular bases that support the pathogen-host interaction, evidencing the contribution of AS regulating the transcriptional host response.

Mechanism of procyanidin b2 in the treatment of chronic myeloid leukemia based on integrating network pharmacology and molecular docking.

To study the pharmacological mechanism of PCB2 on chronic myeloid leukemia (CML) by integrating network pharmacological methods systematically. Firstly, the potential target genes of PCB2 were predicted by the pharmacological database and analysis platform (TCMSP and Pharmmapper). Meanwhile, the relevant target genes of CML were collected from GeneCards and DisGene. Pooled data were collected to screen for common target genes. Furthermore, the above intersection genes were imported into the String website to construct a protein-protein interaction (PPI) network, and the Gene Ontology (GO) functional annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway were further analyzed. Besides, molecular docking was performed to verify the possible binding conformation between PCB2 and candidate targets. Finally, MTT and RT-PCR experiments of K562 cells were performed to verify the above results of network pharmacology. A total of 229 PCB2 target genes were retrieved, among which 186 target genes had interaction with CML. The pharmacological effects of PCB2 on CML were related to some important oncogenes and signaling pathways. The top ten core targets predicted by Network Analysis were as follows: AKT1, EGFR, ESR1, CASP3, SRC, VEGFA, HIF1A, ERBB2, MTOR, and IGF1. Molecular docking studies confirmed that hydrogen bonding was the main interaction force of PCB2 binding targets. According to the molecular docking score, the following three target proteins were most likely to bind to PCB2: VEGFA (-5.5 kcal/mol), SRC (-5.1 kcal/mol), and EGFR (-4.6 kcal/mol). After treatment of PCB2 for 24h, mRNA expression levels of VEGFA and HIF1A decreased significantly in K562 cells. Through integrating network pharmacology combined with molecular docking, the study revealed the potential mechanism of PCB2 anti-chronic myeloid leukemia.

The Resistance and Virulence Characteristics of Salmonella Enteritidis Strain Isolated from Patients with Food Poisoning Based on the Whole-Genome Sequencing and Quantitative Proteomic Analysis.

This paper explores the drug resistance, genome and proteome expression characteristics of Salmonella from a food poisoning event. A multidrug-resistant Salmonella Enteritidis strain, labeled as 27A, was isolated and identified from a food poisoning patient. Antimicrobial susceptibility testing determined the resistance of 27A strain to 14 antibiotics. Then, WGS analysis and comparative genomics analysis were performed on 27A, and the functional annotation of resistance genes, virulence genes were performed based on VFDB, ARDB, COG, CARD, GO, KEGG, and CAZY databases. Meanwhile, based on iTRAQ technology, quantitative proteomic analysis was conducted on 27A to analyze the functions and interactions of differentially expressed proteins related to bacterial resistance and pathogenicity. Strain 27A belonged to ST11 S. Enteritidis and was resistant to levofloxacin, ciprofloxacin, ampicillin, piperacillin, and ampicillin/sulbactam. There were 33 drug resistance genes, 384 virulence genes and 2 plasmid replicon, IncFIB(S) and IncFII(S), annotated by WGS. Proteomic analysis revealed significant changes in virulence and drug proteins, which were mainly involved in bacterial pathogenicity and metabolic processes. PPI prediction showed the relationship between virulence proteins and T3SS proteins, and PagN cooperated with proteins related to T3SS to jointly mediate the invasion of 27A strain on the human body. Phylogenetic analysis indicated that S. Enteritidis has potential transmission in humans, food, and animals. This study comprehensively analyzed the drug resistance and virulence phenotypes of S. Enteritidis 27A using genomic and proteomic approaches. These helps reveal the drug resistance and virulence mechanisms of S. Enteritidis, and provides important information for the source tracing and the prevention of related diseases, which lays a foundation for research on food safety, public health monitoring, and the drug resistance and pathogenicity of S. Enteritidis.

Bioinformatics screening of gene expression profile and diagnostic application of meningeal carcinoma.

The aim of this study was to research gene expression profiles and diagnostic applications of meningeal carcinoma based on bioinformatics. We used the Gene Expression Omnibus (GEO) database to obtain the GSE43290 dataset based on the expression data of normal meninges and meningiomas consisting of 51 samples divided into two groups (47 samples of meningioma tumors and four samples of normal meninges). We used the GEO2R tool to identify differentially expressed genes (DEGs) by setting the log2 fold change as greater than two and an adjusted p-value lower than 0.05. We used the database for annotation, visualization and integrated discovery (DAVID) to perform gene ontology, biological pathways and functional annotation of the DEGs. A search Tool for the Retrieval of Interacting Gene database (STRING) was used to obtain Protein-Protein Interaction (PPI) and modular networks based on the Markov clustering algorithm. Our study identified 358 significant DEGs, of which 343 were downregulated genes while 15 were upregulated. Five significant hub genes (CXCL8, AGT, CXCR4, CXCL12 and CXCL2) were associated with various biological pathways, molecular functions and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. The DEGs were enriched in biological pathways of chemokine-mediated signaling, positive regulation of leukocyte chemotaxis, second messenger-mediated signaling, induction of positive chemotaxis, CXCR chemokine receptor binding and activities of cytokines. These hub genes and pathways could be targeted in clinical research to discover new treatments for meningeal carcinoma.

Fine mapping of the major gene BhHLS1 controlling seed size in wax gourd (Benincasa hispida).

The seed size of wax gourds is an important agronomic trait; however, the associated genes have not yet been reported. In this study, we used a high-density genetic map constructed based on F8 recombinant inbred line populations derived from a cross between MY-1 (large seed) and GX-71 (small seed) strains to detect quantitative trait locis (QTLs) for seed-size-related traits in wax gourd over a two-year period. Two stable QTLs (qSL10 and qSW10) for seed length (SL) and seed width (SW) on chromosome 10 were repeatedly detected over two years (2021-2022). qSL10 had a phenotypic variation rate of 75.30% and 80.80% in 2021 and 2022, respectively. Whereas, qSW10 had a phenotypic variation rate of 66.60% and 73.80% in 2021 and 2022, respectively. Further, a single nucleotide polymorphism mutation was found to cause early termination of Bch10G006400 (BhHLS1) translation in GX-71 through sequencing analysis of candidate genes. Based on gene functional annotation and quantitative real-time PCR analyses, BhHLS1 encoded a probable N-acetyltransferase HLS1-like protein and its expression level was significantly different between parents. Therefore, BhHLS1 is a major candidate gene associated with a one-factor polymorphism regulating the SL and SW of wax gourds. Finally, based on variation in the BhHLS1 sequence, a cleaved amplified polymorphic sequence marker was developed for the molecular marker-assisted breeding of wax gourds. Overall, this study is of great significance for the genetic improvement of seed size, verification of gene functions, and cultivation of specific germplasm resources for wax gourds.

The chromosome-level genome assembly of Fraxinus americana provides insights into the evolution of Oleaceae plants

White ash (Fraxinus americana linn.) originates from the southeastern United States. It is a tall and fast-growing tree species with strong salt-alkali resistance and cold tolerance, making it an important reforestation species and widely planted worldwide. Here, we completed the chromosome-level reference genome assembly of F. americana based on Illumina, PacBio, and Hi-C reads, with a genome size of 878.98 Mb, an N50 of 3.27 Mb, and a heterozygosity rate of 0.3 %. Based on de novo prediction, transcriptome prediction, and homology-based protein prediction, we obtained 39,538 genes. Approximately 843.21 Mb of the assembly genome was composed of 37,928 annotated protein-coding genes, with a gene function annotation rate of 95.93 %. 99.94 % of the overlap clusters (877.44 Mb) were anchored to 23 chromosomes. Synteny analysis of F. americana and other Oleaceae plants showed that F. americana underwent frequent chromosome rearrangements. The amplification of the Ale transposons effectively promoted the genome size of F. americana. Compared with other Oleaceae plants, the Glutathione S-transferase (GST) gene family in the F. americana genome has undergone significant expansion, which may help F. americana cope with adverse natural environments. Furthermore, we found that key enzyme-coding gene families related to lignin biosynthesis were expanded and highly expressed in F. americana leaves. These key genes drive lignin synthesis and benefit F. americana in fast-growing, as well as resisting biotic and abiotic stress. Overall, the F. americana genome assembly provides insights into the evolution of Oleaceae plants and provides abundant resources for breeding and germplasm conservation of white ash.

Single-cell multi-omics sequencing reveals chromosome copy number inconsistency between trophectoderm and inner cell mass in human reconstituted embryos after spindle transfer.

Is the chromosome copy number of the trophectoderm (TE) of a human reconstituted embryos after spindle transfer (ST) representative of the inner cell mass (ICM)? Single-cell multi-omics sequencing revealed that ST blastocysts have a higher proportion of cell lineages exhibiting intermediate mosaicism than conventional ICSI blastocysts, and that the TE of ST blastocysts does not represent the chromosome copy number of ICM. Preimplantation genetic testing for aneuploidy (PGT-A) assumes that TE biopsies are representative of the ICM, but the TE and ICM originate from different cell lineages, and concordance between TE and ICM is not well-studied, especially in ST embryos. We recruited 30 infertile women who received treatment at our clinic and obtained 45 usable blastocysts (22 from conventional ICSI and 23 reconstituted embryos after ST). We performed single-cell multi-omics sequencing on all blastocysts to predict and verify copy number variations (CNVs) in each cell. We determined the chromosome copy number of each embryo by analysing the proportion of abnormal cells in each blastocyst. We used the Bland-Altman concordance and the Kappa test to evaluate the concordance between TE and ICM in the both groups. The study was conducted at a public tertiary hospital in China, where all the embryo operations, including oocytes retrieval, ST, and ICSI, were performed in the embryo laboratory. We utilized single-cell multi-omics sequencing technology at the Biomedical Pioneering Innovation Center, School of Life Sciences, Peking University, to analyse the blastocysts. Transcriptome sequencing was used to predict the CNV of each cell through bioinformatics analysis, and the results were validated using the DNA methylation library of each cell to confirm chromosomal normalcy. We conducted statistical analysis and graphical plotting using R 4.2.1, SPSS 27, and GraphPad Prism 9.3. Mean age of the volunteers, the blastocyst morphology, and the developmental ratewere similar in ST and ICSI groups. The blastocysts in the ST group had some additional chromosomal types that were prone to variations beyond those enriched in the blastocysts of the ICSI group. Finally, both Bland-Altman concordance test and kappa concordancetest showed good chromosomal concordance between TE and ICM in the ICSI blastocysts (kappa = 0.659, P < 0.05), but not in ST blastocysts (P = 1.000), suggesting that the TE in reconstituted embryos is not representative of ICM. Gene functional annotation (GO and KEGG analyses) suggests that there may be new or additional pathways for CNV generation in ST embryos compared to ICSI embryos. This study was mainly limited by the small sample size and the limitations of single-cell multi-omics sequencing technology. To select eligible single cells, some cells of the embryos were eliminated or not labelled, resulting in a loss of information about them. The findings of this study are innovative and exploratory. A larger sample size of human embryos (especially ST embryos) and more accurate molecular genetics techniques for detecting CNV in single cells are needed to validate our results. Our study justifies the routine clinical use of PGT-A in ICSI blastocysts, as we found that the TE is a good substitute for ICM in predicting chromosomal abnormalities. While PGT-A is not entirely accurate, our data demonstrate good clinical feasibility. This trial was able to provide correct genetic counselling to patients regarding the reliability of PGT-A. Regarding ST blastocysts, the increased mosaicism rate and the inability of the TE to represent the chromosomal copy number of the ICM are both biological characteristics that differentiate them from ICSI blastocysts. Currently, ST is not used clinically on a large scale to produce blastocysts. However, if ST becomes more widely used in the future, our study will be the first to demonstrate that the use of PGT-A in ST blastocysts may not be as accurate as PGT-A for ICSI blastocysts. This study was supported by grants from the National Key R&D Program of China (2018YFA0107601) and the National Key R&D Program of China (2018YFC1003003). The authors declare no conflict of interest. N/A.

Metabolic pathways engineering for drought or/and heat tolerance in cereals.

Drought (D) and heat (H) are the two major abiotic stresses hindering cereal crop growth and productivity, either singly or in combination (D/+H), by imposing various negative impacts on plant physiological and biochemical processes. Consequently, this decreases overall cereal crop production and impacts global food availability and human nutrition. To achieve global food and nutrition security vis-a-vis global climate change, deployment of new strategies for enhancing crop D/+H stress tolerance and higher nutritive value in cereals is imperative. This depends on first gaining a mechanistic understanding of the mechanisms underlying D/+H stress response. Meanwhile, functional genomics has revealed several stress-related genes that have been successfully used in target-gene approach to generate stress-tolerant cultivars and sustain crop productivity over the past decades. However, the fast-changing climate, coupled with the complexity and multigenic nature of D/+H tolerance suggest that single-gene/trait targeting may not suffice in improving such traits. Hence, in this review-cum-perspective, we advance that targeted multiple-gene or metabolic pathway manipulation could represent the most effective approach for improving D/+H stress tolerance. First, we highlight the impact of D/+H stress on cereal crops, and the elaborate plant physiological and molecular responses. We then discuss how key primary metabolism- and secondary metabolism-related metabolic pathways, including carbon metabolism, starch metabolism, phenylpropanoid biosynthesis, γ-aminobutyric acid (GABA) biosynthesis, and phytohormone biosynthesis and signaling can be modified using modern molecular biotechnology approaches such as CRISPR-Cas9 system and synthetic biology (Synbio) to enhance D/+H tolerance in cereal crops. Understandably, several bottlenecks hinder metabolic pathway modification, including those related to feedback regulation, gene functional annotation, complex crosstalk between pathways, and metabolomics data and spatiotemporal gene expressions analyses. Nonetheless, recent advances in molecular biotechnology, genome-editing, single-cell metabolomics, and data annotation and analysis approaches, when integrated, offer unprecedented opportunities for pathway engineering for enhancing crop D/+H stress tolerance and improved yield. Especially, Synbio-based strategies will accelerate the development of climate resilient and nutrient-dense cereals, critical for achieving global food security and combating malnutrition.

Characterization and fine mapping analysis of a major stable QTL qKnps-4A for kernel number per spike in wheat.

A major stable QTL for kernel number per spike was narrowed down to a 2.19-Mb region containing two potential candidate genes, and its effects on yield-related traits were characterized. Kernel number per spike (KNPS) in wheat is a key yield component. Dissection and characterization of major stable quantitative trait loci (QTLs) for KNPS would be of considerable value for the genetic improvement of yield potential using molecular breeding technology. We had previously reported a major stable QTL controlling KNPS, qKnps-4A. In the current study, primary fine-mapping analysis, based on the primary mapping population, located qKnps-4A to an interval of approximately 6.8-Mb from 649.0 to 655.8 Mb on chromosome 4A refering to 'Kenong 9204' genome. Further fine-mapping analysis based on a secondary mapping population narrowed qKnps-4A to an approximately 2.19-Mb interval from 653.72 to 655.91 Mb. Transcriptome sequencing, gene function annotation analysis and homologous gene related reports showed that TraesKN4A01HG38570 and TraesKN4A01HG38590 were most likely to be candidate genes of qKnps-4A. Phenotypic analysis based on paired near-isogenic lines in the target region showed that qKnps-4A increased KNPS mainly by increasing the number of central florets per spike. We also evaluated the effects of qKnps-4A on other yield-related traits. Moreover, we dissected the QTL cluster of qKnps-4A and qTkw-4A and proved that the phenotypic effects were probably due to close linkage of two or more genes rather than pleiotropic effects of a single gene. This study provides molecular marker resource for wheatmolecular breeding designed to improve yield potential, and lay the foundation for gene functional analysis of qKnps-4A.

Identification of mobile genetic elements with geNomad

Identifying and characterizing mobile genetic elements in sequencing data is essential for understanding their diversity, ecology, biotechnological applications and impact on public health. Here we introduce geNomad, a classification and annotation framework that combines information from gene content and a deep neural network to identify sequences of plasmids and viruses. geNomad uses a dataset of more than 200,000 marker protein profiles to provide functional gene annotation and taxonomic assignment of viral genomes. Using a conditional random field model, geNomad also detects proviruses integrated into host genomes with high precision. In benchmarks, geNomad achieved high classification performance for diverse plasmids and viruses (Matthews correlation coefficient of 77.8% and 95.3%, respectively), substantially outperforming other tools. Leveraging geNomad’s speed and scalability, we processed over 2.7 trillion base pairs of sequencing data, leading to the discovery of millions of viruses and plasmids that are available through the IMG/VR and IMG/PR databases. geNomad is available at https://portal.nersc.gov/genomad.

Shuangshen Ningxin capsule alleviates myocardial ischemia–reperfusion injury in miniature pigs by modulating mitophagy: network pharmacology and experiments in vivo

BackgroundMyocardial ischemia/reperfusion injury (MI/RI) is involved in a variety of pathological states for which there is no effective treatment exists. Shuangshen Ningxin (SSNX) capsule which is developed by Xiyuan Hospital, Chinese Academy of Traditional Chinese Medicine has been demonstrated to alleviate MI/RI, but its mechanism remains to be further elucidated.MethodsThe MI/RI miniature pigs model was constructed to assess the pharmacodynamics of SSNX by blocking the proximal blood flow of the left anterior descending branch of the cardiac coronary artery through an interventional balloon. The principal chemical compounds and potential targets of SSNX were screened by HPLC–MS and SwissTargetPrediction. The targets of MI/RI were identified based on Online Mendelian Inheritance in Man (OMIM) and GeneCards. Cytoscape 3.9.0 was applied to construct a protein–protein interaction (PPI) network, and Gene Ontology (GO) functional annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were performed using metascape. To further validate the mechanism of SSNX, Molecular docking, Transmission electron microscopy, and Western blot analysis were used to test the effectiveness of targets in related pathways.ResultsOur results indicated that SSNX significantly improved cardiac function, attenuated myocardial I/R injury. Through network analysis, a total of 15 active components and 201 targets were obtained from SSNX, 75 of which are potential targets for the treatment of MI/RI. KEGG and MCODE analysis showed that SSNX is involved in the mitophagy signaling pathway, and ginsenoside Rg1, ginsenoside Rb1 and ginsenoside Rb2 are key components associated with the mitophagy. Further experimental results proved that SSNX protected mitochondrial structure and function, and significantly reduced the expression of mitophagy-related proteins PTEN-induced putative kinase 1 (PINK1), Parkin, FUN14 domain containing 1 (FUNDC1) and Bcl-2/E1B-19 kDa interacting protein 3 (BNIP3) in MI/RI miniature pigs.ConclusionIn our study, the integration of network pharmacology and experiments in vivo demonstrated that SSNX interfered with MI/RI by inhibiting mitophagy.

The ester production capacity of Pichia kudriavzevii based on functional annotation of genes.

Esters were identified as the primary volatile flavor compounds in Chinese Baijiu, exerting a significant influence on its quality and aroma. This study focused on the yeast strain Pichia kudriavzevii, renowned for its high capacity to produce esters. Whole genome sequences were annotated and analyzed using the GO, KEGG, KOG, CAZy, and Pfam databases to determine the genetic basis underly the enhanced ester production capacity. Results showed that P. kudriavzevii gene function was concentrated in biosynthetic capacity, metabolic capacity, amino acid translocation capacity, glycoside hydrolysis capacity and transfer capacity. Additionally, acyltransferase and kinase were predicted as active sites contributing to P. kudriavzevii high ester production. We further compared the volatile composition differences between P. kudriavzevii and Saccharomyces cerevisiae through Headspace solid-phase microextraction-gas Chromatography-mass spectrometry (HS-SPME-GC-MS), revealing P. kudriavzevii produced 3.5 times more esters than S. cerevisiae. Overall, our findings suggest that P. kudriavzevii had potential applications in the Baijiu brewing industry.

SLA2 is a prognostic marker in HNSCC and correlates with immune cell infiltration in the tumor microenvironment

PurposeTo investigate Src-like adaptor 2 gene (SLA2) expression in head and neck squamous cell carcinoma (HNSCC), its potential prognostic value, and its effect on immune cell infiltration.MethodsThrough a variety of bioinformatics analyses, we extracted and analyzed data sets from the Cancer Genome Atlas (TCGA), Tumor Immune Estimation Resource (TIMER), and Gene Expression Profile Interaction Analysis (GEPIA) to analyze the correlation between SLA2 and the prognosis, immune checkpoint, tumor microenvironment (TME) and immune cell infiltration of HNSCC, and to explore its potential oncogenic mechanism. To further explore the potential role of SLA2 in HNSCC by Gene ontology (GO) functional annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis.ResultsSLA2 messenger ribonucleic acid (mRNA) levels were increased in HNSCC tumor tissues compared with normal tissues. In addition, we found that SLA2 may be an independent prognostic factor for HNSCC, and high SLA2 expression is associated with favorable prognosis in HNSCC. SLA2 expression was positively correlated with B cells, cluster of differentiation 8-positive T cells (CD8 + T cells), cluster of differentiation 4-positive T cells (CD4 + T cells), macrophages, neutrophil and dendritic cells infiltration. SLA2 has also been shown to co-express immune-related genes and immune checkpoints. Significant GO term analysis by Gene Set Enrichment Analysis (GSEA) indicated that genes correlated with SLA2 were located mainly in the side of membrane, receptor complex, secretory granule membrane, endocytic vesicle, membrane region, and endosome membrane, where they were involved in leukocyte cell–cell adhesion, response to interferon-gamma, and regulation of immune effector process. These related genes also served as antigen binding, cytokine receptor activity, phosphatidylinositol 3-kinase activity, peptide receptor activity, Src homology domain 3 (SH3) domain binding, and cytokine receptor binding. KEGG pathway analysis demonstrated that these genes related to SLA2 were mainly enriched in signal pathways, such as hematopoietic cell lineage, cell adhesion molecules (CAMs), natural killer cell mediated cytotoxicity, measles, and chemokine signaling pathway.ConclusionsSLA2 is increased in HNSCC, and high SLA2 expression is associated with favorable prognosis. SLA2 may affect tumor development by regulating tumor infiltrating cells in TME. SLA2 may be a potential target for immunotherapy.

Cyclin dependent kinase inhibitor 3 (CDKN3) upregulation is associated with unfavorable prognosis in clear cell renal cell carcinoma and shapes tumor immune microenvironment: A bioinformatics analysis.

Cell cycle regulatory proteins plays a pivotal role in the development and progression of many human malignancies. Identification of their biological functions as well as their prognostic utility presents an active field of research. As a continuation of the ongoing efforts to elucidate the molecular characteristics of clear cell renal cell carcinoma (ccRCC); we present a comprehensive bioinformatics study targeting the prognostic and mechanistic role of cyclin-dependent kinase inhibitor 3 (CDKN3) in ccRCC. The ccRCC cohort from the Cancer Genome Atlas Program was accessed through the UCSC Xena browser to obtain CDKN3 mRNA expression data and their corresponding clinicopathological variables. The independent prognostic signature of CDKN3 was evaluated using univariate and multivariate Cox logistic regression analysis. Gene set enrichment analysis and co-expression gene functional annotations were used to discern CDKN3-related altered molecular pathways. The tumor immune microenvironment was evaluated using TIMER 2.0 and gene expression profiling interactive analysis. CDKN3 upregulation is associated with shortened overall survival (hazard ratio [HR] = 2.325, 95% confident interval [CI]: 1.703-3.173, P < .0001) in the Cancer Genome Atlas Program ccRCC cohort. Univariate (HR: 0.426, 95% CI: 0.316-0.576, P < .001) and multivariate (HR: 0.560, 95% CI: 0.409-0.766, P < .001) Cox logistic regression analyses indicate that CDKN3 is an independent prognostic variable of the overall survival. High CDKN3 expression is associated with enrichment within the following pathways including allograph rejection, epithelial-mesenchymal transition, mitotic spindle, inflammatory response, IL-6/JAK/STAT3 signaling, spermatogenesis, TNF-α signaling via NF-kB pathway, complement activation, KRAS signaling, and INF-γ signaling. CDKN3 is also associated with significant infiltration of a wide spectrum of immune cells and correlates remarkably with immune-related genes. CDKN3 is a poor prognostic biomarker in ccRCC that alters many molecular pathways and impacts the tumor immune microenvironment.

Mechanism of Astragali Radix-Curcumae Rhizoma in treating gastric cancer based on network pharmacology and experimental verification

This study aims to investigate the mechanism of Astragali Radix-Curcumae Rhizoma(HQEZ) in the treatment of gastric cancer based on network pharmacology. Further, the SGC7901 cell model of gastric cancer was employed to validate the efficacy and key targets of the herb pair. Firstly, the CCK-8 assay was employed to evaluate the direct effect of HQEZ on the proliferation of gastric cancer SGC7901 cells. Then, network pharmacology methods were employed to investigate the active ingredients, key targets, and key signaling pathways involved in the treatment of gastric cancer with HQEZ. The results showed that HQEZ contained 18 potential active ingredients, such as quercetin, naringenin, and curcumin. The results of gene ontology(GO) functional annotation and Kyoto Encyclopedia of Genes and Genomes(KEGG) pathway enrichment suggested that the main targets of HQEZ in treating gastric cancer were involved in the regulation of protein serine/threonine kinase activity, activation of mitogen-activated protein kinase(MAPK) activity, cysteine-type endopeptidase activity, and negative regulation of protein serine/threonine kinase activity. The hypoxia-inducible factor-1(HIF-1) signaling pathway, ATP-binding cassette(ABC) transporters, cytochrome P450-mediated metabolism of xenobiotics, p53 signaling pathway, and cell apoptosis were key signaling pathways of HQEZ in treating gastric cancer. The cell experiments demonstrated that HQEZ significantly downregulated the expression of ATP-binding cassette subfamily B member 1(ABCB1), epidermal growth factor receptor(EGFR), phosphorylated serine/threonine kinase(p-AKT), hypoxia inducible factor 1 subunit alpha(HIF1A), B-cell lymphoma 2(BCL2), breast cancer susceptibility protein 1(BRCA1), DNA polymerase theta(POLH), ribonucleotide reductase M1(RRM1), and excision repair cross-complementation group 1(ERCC1), and upregulated the expression of tumor protein P53(TP53) and cysteinyl aspartate-specific proteinase(CAPS3). Finally, a multivariate COX regression model was adopted to study the relationship between gene expression and clinical information data of gastric cancer patients in the TCGA database, which demonstrated that the key targets of HQEZ were associated with the poor prognosis in gastric cancer patients. Further feature selection using the LASSO algorithm showed that EGFR, HIF1A, TP53, POLH, RRM1, and ERCC1 were closely associated with the survival of gastric can-cer patients. In conclusion, HQEZ regulates the expression of genes involved in DNA repair, survival, and apoptosis in gastric cancer cells via multiple targets and pathways, assisting the treatment of gastric cancer.

A Bacillus subtilis Strain ZJ20 with AFB1 Detoxification Ability: A Comprehensive Analysis.

As a class I carcinogen, aflatoxin can cause serious damage to various tissues and organs through oxidative stress injuries. The liver, as the target organ of AFB1, is the most seriously damaged. Biological methods are commonly used to degrade AFB1. In our study, the aflatoxin B1-degrading strain ZJ20 was screened from AFB1-contaminated feed and soil, and the degradation of AFB1 by ZJ20 was investigated. The whole genome of strain ZJ20 was analyzed, revealing the genomic complexity of strain ZJ20. The 16S rRNA analysis of strain ZJ20 showed 100% identity to Bacillus subtilis IAM 12118. Through whole gene functional annotation, it was determined that ZJ20 has high antioxidant activity and enzymatic activity; more than 100 CAZymes and 11 gene clusters are involved in the production of secondary metabolites with antimicrobial properties. In addition, B. subtilis ZJ20 was predicted to contain a cluster of genes encoding AFB1-degrading enzymes, including chitinase, laccase, lactonase, and manganese oxidase. The comprehensive analysis of B. subtilis provides a theoretical basis for the subsequent development of the biological functions of ZJ20 and the combinatorial enzyme degradation of AFB1.

Integrative miRNA-mRNA network analysis to identify crucial pathways of salinity adaptation in brain transcriptome of Labeo rohita.

Introduction: Brain being the master regulator of the physiology of animal, the current study focuses on the gene expression pattern of the brain tissue with special emphasis on regulation of growth, developmental process of an organism and cellular adaptation of Labeo rohita against unfavourable environmental conditions. Methods: RNA-seq study was performed on collected brain samples at 8ppt salt concentration and analyzed for differential gene expression, functional annotation and miRNA-mRNA regulatory network. Results: We found that 2450 genes were having significant differential up and down regulation. The study identified 20 hub genes based on maximal clique centrality algorithm. These hub genes were mainly involved in various signaling pathways, energy metabolism and ion transportation. Further, 326 up and 1214 down regulated genes were found to be targeted by 7 differentially expressed miRNAs i.e., oni-miR-10712, oni-miR-10736, ssa-miR-221-3p, ssa-miR-130d-1-5p, ssa-miR-144-5p and oni-miR-10628. Gene ontology analysis of these differentially expressed genes led to the finding that these genes were involved in signal transduction i.e., calcium, FOXO, PI3K-AKT, TGF-β, Wnt and p53 signalling pathways. Differentially expressed genes were also involved in regulation of immune response, environmental adaptation i.e., neuroactive ligand-receptor interaction, ECM-receptor interaction, cell adhesion molecules and circadian entrainment, osmoregulation and energy metabolism, which are critical for salinity adaptation. Discussion: The findings of whole transcriptomic study on brain deciphered the miRNA-mRNA interaction patterns and pathways associated with salinity adaptation of L. rohita.

Hub Gene Identification and Heat-Stress-Related Transcriptional Regulation Mechanism in Cabbage (Brassica oleracea L.)

Cabbage is a heat-stress-sensitive cool-season crop. When exposed to high temperatures, cabbage plants can experience reduced growth, wilting, leaf yellowing, and premature bolting (the formation of a flowering stem). The regulatory mechanism controlling heat stress is poorly understood in cabbage. To investigate this mechanism, physiological changes and transcriptional profiling for different heat treatment times were analyzed in this study. The results showed that superoxide dismutase (SOD), catalase (CAT), and peroxidase (POD) activities were enhanced under heat stress. In particular, the increase in SOD and POD activities after 12 h of heat treatment was greater than that after 4 h of heat treatment. With increasing heat treatment duration, the leaf CAT activity and H2O2 content decreased after the initial increase. The electrolyte leakage and malondialdehyde (MDA) content dropped significantly, while the proline content increased. Alongside that, 7007 and 5537 upregulated genes were identified in the experimental group treated with heat stress for the 4 h and 12 h treatments, respectively. We found that 10,479 DEGs were shared in the heat stress treatment, of which 1241 were associated with the heat treatment time. By integrating the expression patterns and functional annotations of genes related to heat stress, we identified 15 hub genes that respond to heat stress in cabbage. Meanwhile, we had constructed a physiological to molecular model of cabbage response to long-term heat stress. These findings provide new insights for the comprehensive analysis of cabbage response to heat stress and genetic resources for breeding new varieties of cabbage with heat tolerance via genetic engineering.

Transcriptome analysis reveals the peptide toxins diversity of Macrothele palpator venom

Spider venom is a large pharmacological repertoire of different bioactive peptide toxins. However, obtaining crude venom from some spiders is challenging. Thus, studying individual toxins through venom purification is a daunting task. In this study, we constructed the cDNA library and transcriptomic sequencing from the Macrothele palpator venom glands. Subsequently, 718 high-quality expressed sequence tags (ESTs) were identified, and grouped into three categories, including 449 toxin-like (62.53 %), 136 cellular component (18.94 %) and 133 non-matched (18.52 %) based on the gene function annotation. Additionally, 112 non-redundant toxin-like peptides were classified into 13 families (families A-M) based on their sequence homology and cysteine framework. Bioinformatics analysis revealed a high sequence similarity between families A-J and the toxins from Macrothele gigas in the NR database. In contrast, families K-M had a generally low sequence homology with known spider peptide toxins and unpredictable biological functions. Taken together, this study adds many new members to the spider toxin superfamily and provides a basis for identifying various potential biological tools in M. palpator venom.

Phylogenetic Analysis Guides Transporter Protein Deorphanization: A Case Study of the SLC25 Family of Mitochondrial Metabolite Transporters.

Homology search and phylogenetic analysis have commonly been used to annotate gene function, although they are prone to error. We hypothesize that the power of homology search in functional annotation depends on the coupling of sequence variation to functional diversification, and we herein focus on the SoLute Carrier (SLC25) family of mitochondrial metabolite transporters to survey this coupling in a family-wide manner. The SLC25 family is the largest family of mitochondrial metabolite transporters in eukaryotes that translocate ligands of different chemical properties, ranging from nucleotides, amino acids, carboxylic acids and cofactors, presenting adequate experimentally validated functional diversification in ligand transport. Here, we combine phylogenetic analysis to profile SLC25 transporters across common eukaryotic model organisms, from Saccharomyces cerevisiae, Caenorhabditis elegans, Drosophila melanogaster, Danio rerio, to Homo sapiens, and assess their sequence adaptations to the transported ligands within individual subfamilies. Using several recently studied and poorly characterized SLC25 transporters, we discuss the potentials and limitations of phylogenetic analysis in guiding functional characterization.