The Gene Ontology Resource: 20 years and still GOing strong

Integrated Bioinformatics Analysis Revealed Immune Checkpoint Genes Relevant to Type 2 Diabetes

Objective

Chronic low-grade inflammation of the pancreatic islets is the characteristic of type 2 diabetes (T2D), and some of the immune checkpoints may play important roles in the pancreatic islet inflammation. Thus, we aim to explore the immune checkpoint genes (ICGs) associated with T2D, thereby revealing the role of ICGs in the pathogenesis of T2D based on bioinformatic analyses.

Methods

Differentially expressed genes (DEGs) and immune checkpoint genes (ICGs) of islets between T2D and control group were screened from datasets of the Gene Expression Omnibus (GEO). A risk model was built based on the coefficients of ICGs calculated by ridge regression. Functional enrichment analysis and immune cell infiltration estimation were conducted. Correlations between ICGs and hub genes, T2D-related disease genes, insulin secretion genes, and beta cell function-related genes were analyzed. Finally, we conducted RT-PCR to verify the expression of these ICGs.

Results

In total, pancreatic islets from 19 cases of T2D and 84 healthy subjects were included. We identified 458 DEGs. Six significantly upregulated ICGs (CD44, CD47, HAVCR2, SIRPA, TNFSF9, and VTCN1) in T2D were screened out. These ICGs were significantly correlated with several hub genes and T2D-related genes; furthermore, they were correlated with insulin secretion and β cell function-related genes. The analysis of immune infiltration showed that the concentrations of eosinophils, T cells CD4 naive, and T cells regulatory (Tregs) were significantly higher, but CD4 memory resting T cells and monocytes were lower in islets of T2D patients. The infiltrated immune cells in T2D pancreatic islet were associated with these six ICGs. Finally, the expression levels of four ICGs were confirmed by RT-PCR, and three ICGs were validated in another independent dataset.

Conclusion

In conclusion, the identified ICGs may play an important role in T2D. Identification of these differential genes may provide new clues for the diagnosis and treatment of T2D.

Spaceflight induces changes in gene expression profiles linked to insulin and estrogen

Organismal adaptations to spaceflight have been characterized at the molecular level in model organisms, including Drosophila and C. elegans. Here, we extend molecular work to energy metabolism and sex hormone signaling in mice and humans. We found spaceflight induced changes in insulin and estrogen signaling in rodents and humans. Murine changes were most prominent in the liver, where we observed inhibition of insulin and estrogen receptor signaling with concomitant hepatic insulin resistance and steatosis. Based on the metabolic demand, metabolic pathways mediated by insulin and estrogen vary among muscles, specifically between the soleus and extensor digitorum longus. In humans, spaceflight induced changes in insulin and estrogen related genes and pathways. Pathway analysis demonstrated spaceflight induced changes in insulin resistance, estrogen signaling, stress response, and viral infection. These data strongly suggest the need for further research on the metabolic and reproductive endocrinologic effects of space travel, if we are to become a successful interplanetary species.

Analysis of the Expression and Prognostic Value of SIRTs in Hepatocellular Carcinoma

Purpose

This study contributes to the evolving understanding of the pivotal involvement of Sirtuins (SIRTs) in various human cancers, with a particular focus on elucidating their expression patterns and clinical relevance within the context of hepatocellular carcinoma (HCC). The investigation involves a comprehensive analysis of mRNA expression and prognostic implications associated with distinct SIRTs in HCC.

Patients and Methods

Initial data pertaining to SIRT expression in HCC patients were collated from publicly accessible databases. Subsequently, the expression levels of select members of the SIRT family were validated using clinicopathological specimens from HCC patients. Additionally, HCC tissue microarray was employed to scrutinize the correlation between SIRT7 expression and HCC prognosis.

Results

The findings indicated a substantial upregulation of SIRT2, SIRT3, SIRT4, SIRT6, and SIRT7 in HCC tissues. Survival analysis underscored a pronounced association between elevated mRNA levels of SIRT3, SIRT6, and SIRT7 and an adverse prognosis for HCC patients. Particularly, SIRT7 emerged as a potential independent risk factor for poor prognosis in HCC patients. Examination of the HCC tissue microarray revealed heightened expression of SIRT7 in 68 cases (54.8%) of HCC tissues. Multivariate analysis established high SIRT7 expression as an independent risk factor for diminished Disease-Free Survival (DFS) and Overall Survival (OS) in HCC patients.

Conclusion

The aberrant expression of SIRT7 presents itself may be as a novel biomarker for predicting the prognosis of HCC patients.

Identification of potential shared gene signatures between gastric cancer and type 2 diabetes: a data-driven analysis

Background

Gastric cancer (GC) and type 2 diabetes (T2D) contribute to each other, but the interaction mechanisms remain undiscovered. The goal of this research was to explore shared genes as well as crosstalk mechanisms between GC and T2D.

Methods

The Gene Expression Omnibus (GEO) database served as the source of the GC and T2D datasets. The differentially expressed genes (DEGs) and weighted gene co-expression network analysis (WGCNA) were utilized to identify representative genes. In addition, overlapping genes between the representative genes of the two diseases were used for functional enrichment analysis and protein-protein interaction (PPI) network. Next, hub genes were filtered through two machine learning algorithms. Finally, external validation was undertaken with data from the Cancer Genome Atlas (TCGA) database.

Results

A total of 292 and 541 DEGs were obtained from the GC (GSE29272) and T2D (GSE164416) datasets, respectively. In addition, 2,704 and 336 module genes were identified in GC and T2D. Following their intersection, 104 crosstalk genes were identified. Enrichment analysis indicated that "ECM-receptor interaction," "AGE-RAGE signaling pathway in diabetic complications," "aging," and "cellular response to copper ion" were mutual pathways. Through the PPI network, 10 genes were identified as candidate hub genes. Machine learning further selected BGN, VCAN, FN1, FBLN1, COL4A5, COL1A1, and COL6A3 as hub genes.

Conclusion

"ECM-receptor interaction," "AGE-RAGE signaling pathway in diabetic complications," "aging," and "cellular response to copper ion" were revealed as possible crosstalk mechanisms. BGN, VCAN, FN1, FBLN1, COL4A5, COL1A1, and COL6A3 were identified as shared genes and potential therapeutic targets for people suffering from GC and T2D.

Urine Proteomic Signatures of Mild Hypothermia Treatment in Cerebral Ischemia-Reperfusion Injury in Rats

Mild hypothermia (MH) is an effective measure to alleviate cerebral ischemia-reperfusion (I/R) injury. However, the underlying biological mechanisms remain unclear. This study set out to investigate dynamic changes in urinary proteome due to MH in rats with cerebral I/R injury and explore the neuroprotective mechanisms of MH. A Pulsinelli's four-vessel occlusion (4-VO) rat model was used to mimic global cerebral I/R injury. Liquid chromatography-tandem mass spectrometry (LC-MS/MS) was employed to profile the urinary proteome of rats with/without MH (32 °C) treatment after I/R injury. Representative differentially expressed proteins (DEPs) associated with MH were validated by western blotting in hippocampus. A total of 597 urinary proteins were identified, among which 119 demonstrated significant changes associated with MH. Gene Ontology (GO) annotation of the DEPs revealed that MH significantly enriched in endopeptidase activity, inflammatory response, aging, response to oxidative stress and reactive oxygen species, blood coagulation, and cell adhesion. Notably, changes in 12 DEPs were significantly reversed by MH treatment. Among them, 8 differential urinary proteins were previously reported to be closely associated with brain disease, including NP, FZD1, B2M, EPCR, ATRN, MB, CA1and VPS4A. Two representative proteins (FZD1, B2M) were further validated by western blotting in the hippocampus and the results were shown to be consistent with urinary proteomic analysis. Overall, this study strengthens the idea that urinary proteome can sensitively reflect pathophysiological changes in the brain, and appears to be the first study to explore the neuroprotective effects of MH by urinary proteomic analysis. FZD1 and B2M may be involved in the most fundamental molecular biological mechanisms of MH neuroprotection.

Epigenetic disruptions in the offspring hypothalamus in response to maternal infection

DNA methylation is an epigenetic modification that can alter gene expression, and the incidence can vary across developmental stages, inflammatory conditions, and sexes. The effects of viral maternal viral infection and sex on the DNA methylation patterns were studied in the hypothalamus of a pig model of immune activation during development. DNA methylation at single-base resolution in regions of high CpG density was measured on 24 individual hypothalamus samples using reduced representation bisulfite sequencing. Differential over- and under-methylated sites were identified and annotated to proximal genes and corresponding biological processes. A total of 120 sites were differentially methylated (FDR-adjusted p-value < 0.05) between maternal infection or sex groups. Among the 66 sites differentially methylated between groups exposed to inflammatory signals and control, most sites were over-methylated in the challenged group and included sites in the promoter regions of genes SIRT3 and NRBP1. Among the 54 differentially methylated sites between females and males, most sites were over-methylated in females and included sites in the promoter region of genes TNC and EIF4G1. The analysis of the genes proximal to the differentially methylated sites suggested that biological processes potentially impacted include immune response, neuron migration and ensheathment, peptide signaling, adaptive thermogenesis, and tissue development. These results suggest that translational studies should consider that the prolonged effect of maternal infection during gestation may be enacted through epigenetic regulatory mechanisms that may differ between sexes.

SNARE chaperone Sly1 directly mediates close-range vesicle tethering

The essential Golgi protein Sly1 is a member of the Sec1/mammalian Unc-18 (SM) family of SNARE chaperones. Sly1 was originally identified through remarkable gain-of-function alleles that bypass requirements for diverse vesicle tethering factors. Employing genetic analyses and chemically defined reconstitutions of ER-Golgi fusion, we discovered that a loop conserved among Sly1 family members is not only autoinhibitory but also acts as a positive effector. An amphipathic lipid packing sensor (ALPS)-like helix within the loop directly binds high-curvature membranes. Membrane binding is required for relief of Sly1 autoinhibition and also allows Sly1 to directly tether incoming vesicles to the Qa-SNARE on the target organelle. The SLY1-20 mutation bypasses requirements for diverse tethering factors but loses this ability if the tethering activity is impaired. We propose that long-range tethers, including Golgins and multisubunit tethering complexes, hand off vesicles to Sly1, which then tethers at close range to initiate trans-SNARE complex assembly and fusion in the early secretory pathway.

Identification of shared gene signatures and pathways for diagnosing osteoporosis with sarcopenia through integrated bioinformatics analysis and machine learning

BACKGROUND

Prior studies have suggested a potential relationship between osteoporosis and sarcopenia, both of which can present symptoms of compromised mobility. Additionally, fractures among the elderly are often considered a common outcome of both conditions. There is a strong correlation between fractures in the elderly population, decreased muscle mass, weakened muscle strength, heightened risk of falls, and diminished bone density. This study aimed to pinpoint crucial diagnostic candidate genes for osteoporosis patients with concomitant sarcopenia.

METHODS

Two osteoporosis datasets and one sarcopenia dataset were obtained from the Gene Expression Omnibus (GEO). Differential expression genes (DEGs) and module genes were identified using Limma and Weighted Gene Co-expression Network Analysis (WGCNA), followed by functional enrichment analysis, construction of protein-protein interaction (PPI) networks, and application of a machine learning algorithm (least absolute shrinkage and selection operator (LASSO) regression) to determine candidate hub genes for diagnosing osteoporosis combined with sarcopenia. Receiver operating characteristic (ROC) curves and column line plots were generated.

RESULTS

The merged osteoporosis dataset comprised 2067 DEGs, with 424 module genes filtered in sarcopenia. The intersection of DEGs between osteoporosis and sarcopenia module genes consisted of 60 genes, primarily enriched in viral infection. Through construction of the PPI network, 30 node genes were filtered, and after machine learning, 7 candidate hub genes were selected for column line plot construction and diagnostic value assessment. Both the column line plots and all 7 candidate hub genes exhibited high diagnostic value (area under the curve ranging from 1.00 to 0.93).

CONCLUSION

We identified 7 candidate hub genes (PDP1, ALS2CL, VLDLR, PLEKHA6, PPP1CB, MOSPD2, METTL9) and constructed column line plots for osteoporosis combined with sarcopenia. This study provides reference for potential peripheral blood diagnostic candidate genes for sarcopenia in osteoporosis patients.

Sulfur Dioxide Alleviates Organ Damage and Inflammatory Response in Cecal Ligation and Puncture-Induced Sepsis Rat

The study aimed to elucidate the mechanisms by which sulfur dioxide (SO2) alleviates organ damage during sepsis using RNA-Seq technology. A cecal ligation and puncture (CLP) sepsis model was established in rats, and the effects of SO2 treatment on organ damage were assessed through histopathological examinations. RNA-Seq was performed to analyze differentially expressed genes (DEGs), and subsequent functional annotations and enrichment analyses were conducted. The CLP model successfully induced sepsis symptoms in rats. Histopathological evaluation revealed that SO2 treatment considerably reduced tissue damage across the heart, kidney, liver, and lungs. RNA-Seq identified 950 DEGs between treated and untreated groups, with significant enrichment in genes associated with ribosomal and translational activities, amino acid metabolism, and PI3K-Akt signaling. Furthermore, gene set enrichment analysis (GSEA) showcased enrichments in pathways related to transcriptional regulation, cellular migration, proliferation, and calcium-ion binding. In conclusion, SO2 effectively mitigates multi-organ damage induced by CLP sepsis, potentially through modulating gene expression patterns related to critical biological processes and signaling pathways. These findings highlight the therapeutic promise of SO2 in managing sepsis-induced organ damage.

A complementary metaproteomic approach to interrogate microbiome cultivated from clinical colon biopsies

The human gut microbiome plays a vital role in preserving individual health and is intricately involved in essential functions. Imbalances or dysbiosis within the microbiome can significantly impact human health and are associated with many diseases. Several metaproteomics platforms are currently available to study microbial proteins within complex microbial communities. In this study, we attempted to develop an integrated pipeline to provide deeper insights into both the taxonomic and functional aspects of the cultivated human gut microbiomes derived from clinical colon biopsies. We combined a rapid peptide search by MSFragger against the Unified Human Gastrointestinal Protein database and the taxonomic and functional analyses with Unipept Desktop and MetaLab-MAG. Across seven samples, we identified and matched nearly 36,000 unique peptides to approximately 300 species and 11 phyla. Unipept Desktop provided gene ontology, InterPro entries, and enzyme commission number annotations, facilitating the identification of relevant metabolic pathways. MetaLab-MAG contributed functional annotations through Clusters of Orthologous Genes and Non-supervised Orthologous Groups categories. These results unveiled functional similarities and differences among the samples. This integrated pipeline holds the potential to provide deeper insights into the taxonomy and functions of the human gut microbiome for interrogating the intricate connections between microbiome balance and diseases.

Evolution is more repeatable in the introduction than range expansion phase of colonization

How repeatable is evolution at genomic and phenotypic scales? We studied the repeatability of evolution during 8 generations of colonization using replicated microcosm experiments with the red flour beetle, Tribolium castaneum. Based on the patterns of shared allele frequency changes that occurred in populations from the same generation or experimental location, we found adaptive evolution to be more repeatable in the introduction and establishment phases of colonization than in the spread phase, when populations expand their range. Lastly, by studying changes in allele frequencies at conserved loci, we found evidence for the theoretical prediction that range expansion reduces the efficiency of selection to purge deleterious alleles. Overall, our results increase our understanding of adaptive evolution during colonization, demonstrating that evolution can be highly repeatable while also showing that stochasticity still plays an important role.

HortGenome Search Engine, a universal genomic search engine for horticultural crops

Horticultural crops comprising fruit, vegetable, ornamental, beverage, medicinal and aromatic plants play essential roles in food security and human health, as well as landscaping. With the advances of sequencing technologies, genomes for hundreds of horticultural crops have been deciphered in recent years, providing a basis for understanding gene functions and regulatory networks and for the improvement of horticultural crops. However, these valuable genomic data are scattered in warehouses with various complex searching and displaying strategies, which increases learning and usage costs and makes comparative and functional genomic analyses across different horticultural crops very challenging. To this end, we have developed a lightweight universal search engine, HortGenome Search Engine (HSE; http://hort.moilab.net), which allows for the querying of genes, functional annotations, protein domains, homologs, and other gene-related functional information of more than 500 horticultural crops. In addition, four commonly used tools, including 'BLAST', 'Batch Query', 'Enrichment analysis', and 'Synteny Viewer' have been developed for efficient mining and analysis of these genomic data.

Combination of disulfiram and Copper-Cysteamine nanoparticles induces mitochondria damage and promotes apoptosis in endometrial cancer

Endometrial cancer (EC) stands as one of the most prevalent gynecological malignancies affecting women, with its incidence and disease-related mortality steadily on the rise. Disulfiram (DSF), an FDA-approved medication primarily used for treating alcohol addiction, has exhibited promising anti-tumor properties. Studies have revealed DSF's capacity for enhanced anti-tumor activity, particularly when combined with copper. The novel Copper-Cysteamine (CuCy) compound, Cu3Cl(SR)2 (R[bond, double bond]CH2CH2NH2), showcases photodynamic effects and demonstrates significant anti-tumor potential under various conditions, including exposure to ultraviolet light, X-ray, microwave, and ultrasound. This study delves into exploring the synergistic anti-tumor effects and underlying mechanisms by utilizing copper-cysteamine in conjunction with DSF against endometrial cancer. The investigation involved comprehensive analyses encompassing in vitro experiments utilizing Ishikawa cells, in vivo studies, and transcriptomic analyses. Remarkably, the combined administration of both compounds at a low dose of 0.5 μM exhibited pronounced efficacy in impeding tumor growth, inhibiting blood vessel formation, and stimulating cell apoptosis. Notably, experiments involving transplanted tumors in nude mice vividly demonstrated the significant in vivo anti-tumor effects of this combination treatment. Detailed examination through transmission electron microscopy unveiled compelling evidence of mitochondrial damage, cellular swelling, and rupture, indicative of apoptotic changes in morphology due to the combined treatment. Moreover, transcriptomic analysis unveiled substantial downregulation of mitochondrial-related genes at the molecular level, coupled with a significant hindrance in the DNA repair pathway. These findings strongly suggest that the combined application of CuCy and DSF induces mitochondrial impairment in Ishikawa cells, thereby fostering apoptosis and ultimately yielding potent anti-tumor effects.

The DNA-repair protein APE1 participates with hnRNPA2B1 to motif-enriched and prognostic miRNA secretion

The base excision repair (BER) Apurinic/apyrimidinic endonuclease 1 (APE1) enzyme is endowed with several non-repair activities including miRNAs processing. APE1 is overexpressed in many cancers but its causal role in the tumorigenic processes is largely unknown. We recently described that APE1 can be actively secreted by mammalian cells through exosomes. However, APE1 role in EVs or exosomes is still unknown, especially regarding a putative regulatory function on vesicular small non-coding RNAs. Through dedicated transcriptomic analysis on cellular and vesicular small RNAs of different APE1-depleted cancer cell lines, we found that miRNAs loading into EVs is a regulated process, dependent on APE1, distinctly conveying RNA subsets into vesicles. We identified APE1-dependent secreted miRNAs characterized by enriched sequence motifs and possible binding sites for APE1. In 33 out of 34 APE1-dependent-miRNA precursors, we surprisingly found EXO-motifs and proved that APE1 cooperates with hnRNPA2B1 for the EV-sorting of a subset of miRNAs, including miR-1246, through direct binding to GGAG stretches. Using TCGA-datasets, we showed that these miRNAs identify a signature with high prognostic significance in cancer. In summary, we provided evidence that the ubiquitous DNA-repair enzyme APE1 is part of the EV protein cargo with a novel post-transcriptional role for this ubiquitous DNA-repair enzyme that could explain its role in cancer progression. These findings could open new translational perspectives in cancer biology.

Defining the relationship between cellular and extracellular vesicle (EV) content in breast cancer via an integrative multi-omic analysis

Much recent research has been dedicated to exploring the utility of extracellular vesicles (EVs) as circulating disease biomarkers. Underpinning this work is the assumption that the molecular cargo of EVs directly reflects the originating cell. Few attempts have been made, however, to empirically validate this on the -omic level. To this end, we have performed an integrative multi-omic analysis of a panel of breast cancer cell lines and corresponding EVs. Whole transcriptome analysis validated that the cellular transcriptome remained stable when cultured cells are transitioned to low serum or serum-free medium for EV collection. Transcriptomic profiling of the isolated EVs indicated a positive correlation between transcript levels in cells and EVs, including disease-associated transcripts. Analysis of the EV proteome verified that HER2 protein is present in EVs, however neither the estrogen (ER) nor progesterone (PR) receptor proteins are detected regardless of cellular expression. Using multivariate analysis, we derived an EV protein signature to infer cellular patterns of ER and HER2 expression, though the ER protein could not be directly detected. Integrative analyses affirmed that the EV proteome and transcriptome captured key phenotypic hallmarks of the originating cells, supporting the potential of EVs for non-invasive monitoring of breast cancers.

Transcriptomic responses of cumulus granulosa cells to SARS-CoV-2 infection during controlled ovarian stimulation

Cumulus granulosa cells (CGCs) play a crucial role in follicular development, but so far, no research has explored the impact of SARS-CoV-2 infection on ovarian function from the perspective of CGCs. In the present study, we compared the cycle outcomes between infected and uninfected female patients undergoing controlled ovarian stimulation, performed bulk RNA-sequencing of collected CGCs, and used bioinformatic methods to explore transcriptomic changes. The results showed that women with SARS-CoV-2 infection during stimulation had significantly lower number of oocytes retrieved and follicle-oocyte index, while subsequent fertilization and embryo development were similar. CGCs were not directly infected by SARS-CoV-2, but exhibited dramatic differences in gene expression (156 up-regulated and 65 down-regulated). Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses demonstrated a high enrichment in antiviral, immune and inflammatory responses with necroptosis. In addition, the pathways related to telomere organization and double strand break repair were significantly affected by infection in gene set enrichment analysis. Further weighted gene co-expression network analysis identified a key module associated with ovarian response traits, which was mainly enriched as a decrease of leukocyte chemotaxis and migration in CGCs. For the first time, our study describes how SARS-CoV-2 infection indirectly affects CGCs at the transcriptional level, which may impair oocyte-CGC crosstalk and consequently lead to poor ovarian response during fertility treatment.

Identification of Immune-Related Gene Signature in Schizophrenia

BACKGROUND

Schizophrenia (SCZ) is a type of psychiatric disorder characterized by multiple symptoms. Our aim is to decipher the relevant mechanisms of immune-related gene signatures in SCZ.

METHODS

The SCZ dataset and its associated immunoregulatory genes were retrieved using Gene Expression Omnibus (GEO) and single-sample gene set enrichment analysis (ssGSEA). Co-expressed gene modules were determined through weighted gene correlation network analysis (WGCNA). To elucidate the functional characteristics of these clusters, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were used. Additionally, gene set enrichment analysis (GSEA) and Gene Set Variation Analysis (GSVA) were conducted to identify enriched pathways for the immune subgroups. A protein-protein interaction (PPI) network analysis was performed to identify core genes relevant to SCZ.

RESULTS

A significantly higher immune score was observed in SCZ compared to control samples. Seven distinct gene modules were identified, with genes highlighted in green selected for further analysis. Using the Cell-type Identification By Estimating Relative Subsets Of RNA Transcripts (CIBERSORT) method, degrees of immune cell adhesion and accumulation related to 22 different immune cell types were calculated. Significantly enriched bioprocesses concerning the immunoregulatory genes with differential expressions included interferon-beta, IgG binding, and response to interferon-gamma, according to GO and KEGG analyses. Eleven hub genes related to immune infiltration emerged as key players among the three top-ranked GO terms.

CONCLUSIONS

This study underscores the involvement of immunoregulatory reactions in SCZ development. Eleven immune-related genes (IFITM1 (interferon induced transmembrane protein 1), GBP1 (guanylate binding protein 1), BST2 (bone marrow stromal cell antigen 2), IFITM3 (interferon induced transmembrane protein 3), GBP2 (guanylate binding protein 2), CD44 (CD44 molecule), FCER1G (Fc epsilon receptor Ig), HLA-DRA (major histocompatibility complex, class II, DR alpha), FCGR2A (Fc gamma receptor IIa), IFI16 (interferon gamma inducible protein 16), and FCGR3B (Fc gamma receptor IIIb)) were identified as hub genes, representing potential biomarkers and therapeutic targets associated with the immune response in SCZ patients.

DNA Electrotransfer Regulates Molecular Functions in Skeletal Muscle

Background

Tissues, such as skeletal muscle, have been targeted for the delivery of plasmid DNA (pDNA) encoding vaccines and therapeutics. The application of electric pulses (electroporation or electrotransfer) increases cell membrane permeability to enhance plasmid delivery and expression. However, the molecular effects of DNA electrotransfer on the muscle tissue are poorly characterized.

Materials and Methods

Four hours after intramuscular plasmid electrotransfer, we evaluated gene expression changes by RNA sequencing. Differentially expressed genes were analyzed by gene ontology (GO) pathway enrichment analysis.

Results

GO analysis highlighted many enriched molecular functions. The terms regulated by pulse application were related to muscle stress, the cytoskeleton and inflammation. The terms regulated by pDNA injection were related to a DNA-directed response and its control. Several terms regulated by pDNA electrotransfer were similar to those regulated by pulse application. However, the terms related to pDNA injection differed, focusing on entry of the plasmid into the cells and intracellular trafficking.

Conclusion

Each muscle stimulus resulted in specific regulated molecular functions. Identifying the unique intrinsic molecular changes driven by intramuscular DNA electrotransfer will aid in the design of preventative and therapeutic gene therapies.

HMGB2 Release Promotes Pulmonary Hypertension and Predicts Severity and Mortality of Patients With Pulmonary Arterial Hypertension

BACKGROUND

Pulmonary hypertension (PH) is a progressive and life-threatening disease characterized by pulmonary vascular remodeling, which involves aberrant proliferation and apoptosis resistance of the pulmonary arterial smooth muscle cells (PASMCs), resembling the hallmark characteristics of cancer. In cancer, the HMGB2 (high-mobility group box 2) protein promotes the pro-proliferative/antiapoptotic phenotype. However, the function of HMGB2 in PH remains uninvestigated.

METHODS

Smooth muscle cell (SMC)-specific HMGB2 knockout or HMGB2-OE (HMGB2 overexpression) mice and HMGB2 silenced rats were used to establish hypoxia+Su5416 (HySu)-induced PH mouse and monocrotaline-induced PH rat models, respectively. The effects of HMGB2 and its underlying mechanisms were subsequently elucidated using RNA-sequencing and cellular and molecular biology analyses. Serum HMGB2 levels were measured in the controls and patients with pulmonary arterial (PA) hypertension.

RESULTS

HMGB2 expression was markedly increased in the PAs of patients with PA hypertension and PH rodent models and was predominantly localized in PASMCs. SMC-specific HMGB2 deficiency or silencing attenuated PH development and pulmonary vascular remodeling in hypoxia+Su5416-induced mice and monocrotaline-treated rats. SMC-specific HMGB2 overexpression aggravated hypoxia+Su5416-induced PH. HMGB2 knockdown inhibited PASMC proliferation in vitro in response to PDGF-BB (platelet-derived growth factor-BB). In contrast, HMGB2 protein stimulation caused the hyperproliferation of PASMCs. In addition, HMGB2 promoted PASMC proliferation and the development of PH by RAGE (receptor for advanced glycation end products)/FAK (focal adhesion kinase)-mediated Hippo/YAP (yes-associated protein) signaling suppression. Serum HMGB2 levels were significantly increased in patients with PA hypertension, and they correlated with disease severity, predicting worse survival.

CONCLUSIONS

Our findings indicate that targeting HMGB2 might be a novel therapeutic strategy for treating PH. Serum HMGB2 levels could serve as a novel biomarker for diagnosing PA hypertension and determining its prognosis.

Identification of Protein-Protein Interaction Associated Functions Based on Gene Ontology

Protein-protein interactions (PPIs) involve the physical or functional contact between two or more proteins. Generally, proteins that can interact with each other always have special relationships. Some previous studies have reported that gene ontology (GO) terms are related to the determination of PPIs, suggesting the special patterns on the GO terms of proteins in PPIs. In this study, we explored the special GO term patterns on human PPIs, trying to uncover the underlying functional mechanism of PPIs. The experimental validated human PPIs were retrieved from STRING database, which were termed as positive samples. Additionally, we randomly paired proteins occurring in positive samples, yielding lots of negative samples. A simple calculation was conducted to count the number of positive samples for each GO term pair, where proteins in samples were annotated by GO terms in the pair individually. The similar number for negative samples was also counted and further adjusted due to the great gap between the numbers of positive and negative samples. The difference of the above two numbers and the relative ratio compared with the number on positive samples were calculated. This ratio provided a precise evaluation of the occurrence of GO term pairs for positive samples and negative samples, indicating the latent GO term patterns for PPIs. Our analysis unveiled several nuclear biological processes, including gene transcription, cell proliferation, and nutrient metabolism, as key biological functions. Interactions between major proliferative or metabolic GO terms consistently correspond with significantly reported PPIs in recent literature.

Identification of shared potential diagnostic markers in asthma and depression through bioinformatics analysis and machine learning

BACKGROUND

There is mounting evidence that asthma might exacerbate depression. We sought to examine candidates for diagnostic genes in patients suffering from asthma and depression.

METHODS

Microarray data were downloaded from the Gene Expression Omnibus(GEO) database and used to screen for differential expressed genes(DEGs) in the SA and MDD datasets. A weighted gene co-expression network analysis(WGCNA) was used to identify the co-expression modules of SA and MDD. The least absolute shrinkage and selection operatoes(LASSO) and support vector machine(SVM) were used to determine critical biomarkers. Immune cell infiltration analysis was used to investigate the correlation between immune cell infiltration and common biomarkers of SA and MDD. Finally, validation of these analytical results was accomplished via the use of both in vivo and in vitro studies.

RESULTS

The number of DEGs that were included in the MDD dataset was 5177, whereas the asthma dataset had 1634 DEGs. The intersection of DEGs for SA and MDD included 351 genes, the strongest positive modules of SA and MDD was 119 genes, which played a function in immunity. The intersection of DEGs and modular hub genes was 54, following the analysis using machine learning algorithms,three hub genes were identified and employed to formulate a nomogram and for the evaluation of diagnostic effectiveness, which demonstrated a significant diagnostic value (area under the curve from 0.646 to 0.979). Additionally, immunocyte disorder was identified by immune infiltration. In vitro studies have revealed that STK11IP deficiency aggravated the LPS/IFN-γinduced up-regulation in M1 macrophage activation.

CONCLUSION

Asthma and MDD pathophysiology may be associated with alterations in inflammatory processes and immune pathways. Additionally, STK11IP may serve as a diagnostic marker for individuals with the two conditions.

Longitudinal analysis of epigenome-wide DNA methylation reveals novel loci associated with BMI change in East Asians

BACKGROUND

Obesity is a global public health concern linked to chronic diseases such as cardiovascular disease and type 2 diabetes (T2D). Emerging evidence suggests that epigenetic modifications, particularly DNA methylation, may contribute to obesity. However, the molecular mechanism underlying the longitudinal change of BMI has not been well-explored, especially in East Asian populations.

METHODS

This study performed a longitudinal epigenome-wide association analysis of DNA methylation to uncover novel loci associated with BMI change in 533 individuals across two Chinese cohorts with repeated DNA methylation and BMI measurements over four years.

RESULTS

We identified three novel CpG sites (cg14671384, cg25540824, and cg10848724) significantly associated with BMI change. Two of the identified CpG sites were located in regions previously associated with body shape and basal metabolic rate. Annotation of the top 20 BMI change-associated CpGs revealed strong connections to obesity and T2D. Notably, these CpGs exhibited active regulatory roles and located in genes with high expression in the liver and digestive tract, suggesting a potential regulatory pathway from genome to phenotypes of energy metabolism and absorption via DNA methylation. Cross-sectional and longitudinal EWAS comparisons indicated different mechanisms between CpGs related to BMI and BMI change.

CONCLUSION

This study enhances our understanding of the epigenetic dynamics underlying BMI change and emphasizes the value of longitudinal analyses in deciphering the complex interplay between epigenetics and obesity.

The axes of biology: a novel axes-based network embedding paradigm to decipher the functional mechanisms of the cell

Summary

Common approaches for deciphering biological networks involve network embedding algorithms. These approaches strictly focus on clustering the genes' embedding vectors and interpreting such clusters to reveal the hidden information of the networks. However, the difficulty in interpreting the genes' clusters and the limitations of the functional annotations' resources hinder the identification of the currently unknown cell's functioning mechanisms. We propose a new approach that shifts this functional exploration from the embedding vectors of genes in space to the axes of the space itself. Our methodology better disentangles biological information from the embedding space than the classic gene-centric approach. Moreover, it uncovers new data-driven functional interactions that are unregistered in the functional ontologies, but biologically coherent. Furthermore, we exploit these interactions to define new higher-level annotations that we term Axes-Specific Functional Annotations and validate them through literature curation. Finally, we leverage our methodology to discover evolutionary connections between cellular functions and the evolution of species.

Availability and implementation

Data and source code can be accessed at https://gitlab.bsc.es/sdoria/axes-of-biology.git.

A Causal Regulation Modeling Algorithm for Temporal Events with Application to Escherichia coli's Aerobic to Anaerobic Transition

Time-series experiments are crucial for understanding the transient and dynamic nature of biological phenomena. These experiments, leveraging advanced classification and clustering algorithms, allow for a deep dive into the cellular processes. However, while these approaches effectively identify patterns and trends within data, they often need to improve in elucidating the causal mechanisms behind these changes. Building on this foundation, our study introduces a novel algorithm for temporal causal signaling modeling, integrating established knowledge networks with sequential gene expression data to elucidate signal transduction pathways over time. Focusing on Escherichia coli's (E. coli) aerobic to anaerobic transition (AAT), this research marks a significant leap in understanding the organism's metabolic shifts. By applying our algorithm to a comprehensive E. coli regulatory network and a time-series microarray dataset, we constructed the cross-time point core signaling and regulatory processes of E. coli's AAT. Through gene expression analysis, we validated the primary regulatory interactions governing this process. We identified a novel regulatory scheme wherein environmentally responsive genes, soxR and oxyR, activate fur, modulating the nitrogen metabolism regulators fnr and nac. This regulatory cascade controls the stress regulators ompR and lrhA, ultimately affecting the cell motility gene flhD, unveiling a novel regulatory axis that elucidates the complex regulatory dynamics during the AAT process. Our approach, merging empirical data with prior knowledge, represents a significant advance in modeling cellular signaling processes, offering a deeper understanding of microbial physiology and its applications in biotechnology.

Tryptophan metabolites relieve intestinal Candida albicans infection by altering the gut microbiota to reduce IL-22 release from group 3 innate lymphoid cells of the colon lamina propria

Invasive candidiasis may be caused by Candida albicans (C. albicans) colonization of the intestinal tract. Preventing intestinal damage caused by Candida albicans infection and protecting intestinal barrier function have become a critical issue. Integrated analyses of the microbiome with metabolome revealed a remarkable shift of the gut microbiota and tryptophan metabolites, kynurenic acid (KynA), and indolacrylic acid (IA) in mice infected with C. albicans. The transcriptome sequencing indicated that differentially expressed genes were significantly associated with innate immune responses and inflammatory responses. The results of this study suggest that KynA and IA (KI) can alleviate intestinal damage caused by Candida albicans infection in mice by reducing intestinal permeability, increasing intestinal firmness, alleviating intestinal inflammation, and reducing the secretion of interleukin-22 (IL-22) in the 3 groups of colon innate lymphoid cells (ILC3). We performed a fecal microbiota transplantation (FMT) experiment and found that the intestinal barrier function, inflammation, and IL-22 secretion of ILC3 in the colon lamina propria of the recipient mice subjected to C. albicans infection and KI treatment were consistent with the trends of the donor mice. Our results suggest that tryptophan metabolites may directly regulate colon lamina ILC3 to promote intestinal resistance to C. albicans invasion, or indirectly regulate the ILC3 secretion of IL-22 to play a protective role in the intestinal barrier by affecting intestinal microorganisms, which may become a potential target for alleviating intestine borne C. albicans infection.

SOS1 tonoplast neo-localization and the RGG protein SALTY are important in the extreme salinity tolerance of Salicornia bigelovii

The identification of genes involved in salinity tolerance has primarily focused on model plants and crops. However, plants naturally adapted to highly saline environments offer valuable insights into tolerance to extreme salinity. Salicornia plants grow in coastal salt marshes, stimulated by NaCl. To understand this tolerance, we generated genome sequences of two Salicornia species and analyzed the transcriptomic and proteomic responses of Salicornia bigelovii to NaCl. Subcellular membrane proteomes reveal that SbiSOS1, a homolog of the well-known SALT-OVERLY-SENSITIVE 1 (SOS1) protein, appears to localize to the tonoplast, consistent with subcellular localization assays in tobacco. This neo-localized protein can pump Na+ into the vacuole, preventing toxicity in the cytosol. We further identify 11 proteins of interest, of which SbiSALTY, substantially improves yeast growth on saline media. Structural characterization using NMR identified it as an intrinsically disordered protein, localizing to the endoplasmic reticulum in planta, where it can interact with ribosomes and RNA, stabilizing or protecting them during salt stress.

Upregulation of GOLPH3 mediated by Bisphenol a promotes colorectal cancer proliferation and migration: evidence based on integrated analysis

Background

The interaction between environmental endocrine-disrupting chemicals, such as Bisphenol A (BPA), and their influence on cancer progression, particularly regarding the GOLPH3 gene in colorectal cancer, remains unclear.

Methods

We performed an integrated analysis of transcriptional profiling, clinical data, and bioinformatics analyses utilizing data from the Comparative Toxicogenomics Database and The Cancer Genome Atlas. The study employed ClueGO, Gene Set Enrichment Analysis, and Gene Set Variation Analysis for functional enrichment analysis, alongside experimental assays to examine the effects of BPA exposure on colorectal cancer cell lines, focusing on GOLPH3 expression and its implications for cancer progression.

Results

Our findings demonstrated that BPA exposure significantly promoted the progression of colorectal cancer by upregulating GOLPH3, which in turn enhanced the malignant phenotype of colorectal cancer cells. Comparative analysis revealed elevated GOLPH3 protein levels in cancerous tissues versus normal tissues, with single-cell analysis indicating widespread GOLPH3 presence across various cell types in the cancer microenvironment. GOLPH3 was also associated with multiple carcinogenic pathways, including the G2M checkpoint. Furthermore, our investigation into the colorectal cancer microenvironment and genomic mutation signature underscored the oncogenic potential of GOLPH3, exacerbated by BPA exposure.

Conclusion

This study provides novel insights into the complex interactions between BPA exposure and GOLPH3 in the context of colorectal cancer, emphasizing the need for heightened awareness and measures to mitigate BPA exposure risks. Our findings advocate for further research to validate these observations in clinical and epidemiological settings and explore potential therapeutic targets within these pathways.

Identification and Analysis of lncRNA and circRNA Related to Wheat Grain Development

The role of lncRNA and circRNA in wheat grain development is still unclear. The objectives of this study were to characterize the lncRNA and circRNA in the wheat grain development and to construct the interaction network among lncRNA, circRNA, and their target miRNA to propose a lncRNA-circRNA-miRNA module related to wheat grain development. Full transcriptome sequencing on two wheat varieties (Annong 0942 and Anke 2005) with significant differences in 1000-grain weight at 10 d (days after pollination), 20 d, and 30 d of grain development were conducted. We detected 650, 736, and 609 differentially expressed lncRNA genes, and 769, 1054, and 1062 differentially expressed circRNA genes in the grains of 10 days, 20 days and 30 days after pollination between Annong 0942 and Anke 2005, respectively. An analysis of the lncRNA-miRNA and circRNA-miRNA targeting networks reveals that circRNAs exhibit a more complex and extensive interaction network in the development of cereal grains and the formation of grain shape. Central to these interactions are tae-miR1177, tae-miR1128, and tae-miR1130b-3p. In contrast, lncRNA genes only form a singular network centered around tae-miR1133 and tae-miR5175-5p when comparing between varieties. Further analysis is conducted on the underlying genes of all target miRNAs, we identified TaNF-YB1 targeted by tae-miR1122a and TaTGW-7B targeted by miR1130a as two pivotal regulatory genes in the development of wheat grains. The quantitative real-time PCR (qRT-PCR) and dual-luciferase reporter assays confirmed the target regulatory relationships between miR1130a-TaTGW-7B and miR1122a-TaNF-YB1. We propose a network of circRNA and miRNA-mediated gene regulation in the development of wheat grains.

DExplore: An Online Tool for Detecting Differentially Expressed Genes from mRNA Microarray Experiments

Microarray experiments, a mainstay in gene expression analysis for nearly two decades, pose challenges due to their complexity. To address this, we introduce DExplore, a user-friendly web application enabling researchers to detect differentially expressed genes using data from NCBI's GEO. Developed with R, Shiny, and Bioconductor, DExplore integrates WebGestalt for functional enrichment analysis. It also provides visualization plots for enhanced result interpretation. With a Docker image for local execution, DExplore accommodates unpublished data. To illustrate its utility, we showcase two case studies on cancer cells treated with chemotherapeutic drugs. DExplore streamlines microarray data analysis, empowering molecular biologists to focus on genes of biological significance.

Suv39h-catalyzed H3K9me3 is critical for euchromatic genome organization and the maintenance of gene transcription

H3K9me3-dependent heterochromatin is critical for the silencing of repeat-rich pericentromeric regions and also has key roles in repressing lineage-inappropriate protein-coding genes in differentiation and development. Here, we investigate the molecular consequences of heterochromatin loss in cells deficient in both SUV39H1 and SUV39H2 (Suv39DKO), the major mammalian histone methyltransferase enzymes that catalyze heterochromatic H3K9me3 deposition. We reveal a paradoxical repression of protein-coding genes in Suv39DKO cells, with these differentially expressed genes principally in euchromatic (Tn5-accessible, H3K4me3- and H3K27ac-marked) rather than heterochromatic (H3K9me3-marked) or polycomb (H3K27me3-marked) regions. Examination of the three-dimensional (3D) nucleome reveals that transcriptomic dysregulation occurs in euchromatic regions close to the nuclear periphery in 3D space. Moreover, this transcriptomic dysregulation is highly correlated with altered 3D genome organization in Suv39DKO cells. Together, our results suggest that the nuclear lamina-tethering of Suv39-dependent H3K9me3 domains provides an essential scaffold to support euchromatic genome organization and the maintenance of gene transcription for healthy cellular function.

Subacromial bursa in patients with rotator cuff tear with or without adhesive capsulitis exhibits a differential transcriptome with potential molecular biomarkers

Background

Although adhesive capsulitis (AC) is a common condition, the pathological mechanisms remain understudied. The purpose of our research was to evaluate variations in gene expression across the entire genome in the subacromial bursa tissue of individuals with rotator cuff tears (RCT), with or without AC, and to explore the factors that may influence the occurrence and progression of AC.

Methods

Transcription profiles of subacromial bursa samples from 12 RCT patients, of whom 6 had also AC, were evaluated. Data were generated using RNA-seq. DESeq2 was utilized to identify the differentially expressed genes (DEGs) in both groups. In order to conduct a more in-depth examination of the DEGs, we performed Gene Ontology (GO) functional enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. A network of interactions between proteins was built, and the central genes were determined using Cytoscape. The hub genes were confirmed through qRT-PCR and immunohistochemistry.

Results

324 of the 16,251 detected genes were identified as DEGs. Analysis of GO functional enrichment showed that the DEGs were enriched in domains of biological process, molecule function and cellular component. Analysis of KEGG pathways revealed enrichment of DEGs in pathways like IL-17 signaling and ECM-receptor interaction. We verified that the association between AC and the increase in expression of the PPI network hub genes.

Conclusion

This study investigated the transcriptome differences of subacromial bursa in RCT patients with or without AC. Using bioinformatics technology, we identified the DEGs and screened out the hub genes. The research enhanced the data on gene expression profiles of DEGs in the subacromial bursa tissue of patients with RCT, offering fresh perspectives on the regulation of gene transcription.

Uncovering the shared neuro-immune-related regulatory mechanisms between spinal cord injury and osteoarthritis

Adults with spinal cord injury (SCI), a destructive neurological injury, have a significantly higher incidence of osteoarthritis (OA), a highly prevalent chronic joint disorder. This study aimed to dissect the neuroimmune-related regulatory mechanisms of SCI and OA using bioinformatics analysis. Using microarray data from the Gene Expression Omnibus database, differentially expressed genes (DEGs) were screened between SCI and sham samples and between OA and control samples. Common DEGs were used to construct a protein-protein interaction (PPI) network. Weighted gene co-expression network analysis (WGCNA) was used to mine SCI- and OA-related modules. Shared miRNAs were identified, and target genes were predicted using the Human MicroRNA Disease Database (HMDD) database. A miRNA-gene-pathway regulatory network was constructed with overlapping genes, miRNAs, and significantly enriched pathways. Finally, the expression of the identified genes and miRNAs was verified using RT-qPCR. In both the SCI and OA groups, 185 common DEGs were identified, and three hub clusters were obtained from the PPI network. WGCNA revealed three SCI-related modules and two OA-related modules. There were 43 overlapping genes between the PPI network clusters and the WGCNA network modules. Seventeen miRNAs shared between patients with SCI and OA were identified. A regulatory network consisting of five genes, six miRNAs, and six signaling pathways was constructed. Upregulation of CD44, TGFBR1, CCR5, and IGF1, while lower levels of miR-125b-5p, miR-130a-3p, miR-16-5p, miR-204-5p, and miR-204-3p in both SCI and OA were successfully verified using RT-qPCR. Our study suggests that a miRNA-gene-pathway network is implicated in the neuroimmune-related regulatory mechanisms of SCI and OA. CD44, TGFBR1, CCR5, and IGF1, and their related miRNAs (miR-125b-5p, miR-130a-3p, miR-16-5p, miR-204-5p, and miR-204-3p) may serve as promising biomarkers and candidate therapeutic targets for SCI and OA.

SLC38A3 Promotes the Proliferation and Migration of Tumor Cells and Predicts Poor Prognosis in Colorectal Cancer

Previous studies have revealed that abnormal expressions of membrane transporters were associated with colorectal cancer (CRC). We herein performed a comprehensive bioinformatics analysis to identify the key transporter protein-related genes involved in CRC and potential mechanisms. Differentially expressed transporter protein-related genes (DE-TPRGs) were identified from CRC and normal samples using The Cancer Genome Atlas database. SLC38A3 expression was validated by immunohistochemistry and RT-qPCR, and the potential mechanism was explored. A total of 63 DE-TPRGs (29 up-regulated and 34 down-regulated) were screened. Inside, ABCC2, ABCG2, SLC4A4, SLC9A3, SLC15A1, and SLC38A3 were identified as hub genes. SLC38A3 is indeed upregulated in colorectal cancer patients. Furthermore, we found that knockdown of SLC38A3 inhibited the proliferation and migration of HCT116 cells, and Hsp70 ATPase activator could rescue it. Overall, SLC38A3 is a novel potential biomarker involved in CRC progression and promotes the proliferation and migration of tumor cells by positively regulating the function of Hsp70.

Nucleobase adducts bind MR1 and stimulate MR1-restricted T cells

MR1T cells are a recently found class of T cells that recognize antigens presented by the major histocompatibility complex-I-related molecule MR1 in the absence of microbial infection. The nature of the self-antigens that stimulate MR1T cells remains unclear, hampering our understanding of their physiological role and therapeutic potential. By combining genetic, pharmacological, and biochemical approaches, we found that carbonyl stress and changes in nucleobase metabolism in target cells promote MR1T cell activation. Stimulatory compounds formed by carbonyl adducts of nucleobases were detected within MR1 molecules produced by tumor cells, and their abundance and antigenicity were enhanced by drugs that induce carbonyl accumulation. Our data reveal carbonyl-nucleobase adducts as MR1T cell antigens. Recognizing cells under carbonyl stress allows MR1T cells to monitor cellular metabolic changes with physiological and therapeutic implications.

Biclustering analysis on tree-shaped time-series single cell gene expression data of Caenorhabditis elegans

BACKGROUND

In recent years, gene clustering analysis has become a widely used tool for studying gene functions, efficiently categorizing genes with similar expression patterns to aid in identifying gene functions. Caenorhabditis elegans is commonly used in embryonic research due to its consistent cell lineage from fertilized egg to adulthood. Biologists use 4D confocal imaging to observe gene expression dynamics at the single-cell level. However, on one hand, the observed tree-shaped time-series datasets have characteristics such as non-pairwise data points between different individuals. On the other hand, the influence of cell type heterogeneity should also be considered during clustering, aiming to obtain more biologically significant clustering results.

RESULTS

A biclustering model is proposed for tree-shaped single-cell gene expression data of Caenorhabditis elegans. Detailedly, a tree-shaped piecewise polynomial function is first employed to fit non-pairwise gene expression time series data. Then, four factors are considered in the objective function, including Pearson correlation coefficients capturing gene correlations, p-values from the Kolmogorov-Smirnov test measuring the similarity between cells, as well as gene expression size and bicluster overlapping size. After that, Genetic Algorithm is utilized to optimize the function.

CONCLUSION

The results on the small-scale dataset analysis validate the feasibility and effectiveness of our model and are superior to existing classical biclustering models. Besides, gene enrichment analysis is employed to assess the results on the complete real dataset analysis, confirming that the discovered biclustering results hold significant biological relevance.

Gene Expression Profile of Benign, Intermediate, and Malignant Spitz and Spitzoid Melanocytic Lesions

Spitz and Spitzoid lesions represent one of the most challenging melanocytic neoplasms in dermatopathology. Nosologic classification has been more recently improved by the discovery of novel molecular drivers, particularly translocations. In the current study, we aimed to use an unbiased approach to explore the gene expression profile of a group of melanocytic Spitz and Spitzoid melanocytic lesions ranging from benign lesions to melanoma, including intermediate lesions such as SPARK nevi and atypical Spitz tumors/melanocytomas. Using unsupervised analysis of gene expression data, we found some distinct hierarchical clusters of lesions, including groups characterized by ALK and NTRK translocations. Few non-ALK translocated tumors demonstrated increased ALK expression, confirmed by immunohistochemistry. Spitz tumors with overlapping features of dysplastic nevi, so-called SPARK nevi, appear to have a common gene expression profile by hierarchical clustering. Finally, weighted gene correlation network analysis identified gene modules variably regulated in subtypes of these cases. Thus, gene expression profiling of Spitz and Spitzoid lesions represents a viable instrument for the characterization of these lesions.

The Arabidopsis Information Resource in 2024

Since 1999, The Arabidopsis Information Resource (www.arabidopsis.org) has been curating data about the Arabidopsis thaliana genome. Its primary focus is integrating experimental gene function information from the peer-reviewed literature and codifying it as controlled vocabulary annotations. Our goal is to produce a "gold standard" functional annotation set that reflects the current state of knowledge about the Arabidopsis genome. At the same time, the resource serves as a nexus for community-based collaborations aimed at improving data quality, access, and reuse. For the past decade, our work has been made possible by subscriptions from our global user base. This update covers our ongoing biocuration work, some of our modernization efforts that contribute to the first major infrastructure overhaul since 2011, the introduction of JBrowse2, and the resource's role in community activities such as organizing the structural reannotation of the genome. For gene function assessment, we used gene ontology annotations as a metric to evaluate: (1) what is currently known about Arabidopsis gene function and (2) the set of "unknown" genes. Currently, 74% of the proteome has been annotated to at least one gene ontology term. Of those loci, half have experimental support for at least one of the following aspects: molecular function, biological process, or cellular component. Our work sheds light on the genes for which we have not yet identified any published experimental data and have no functional annotation. Drawing attention to these unknown genes highlights knowledge gaps and potential sources of novel discoveries.

Identification and coregulation pattern analysis of long noncoding RNAs in the mouse brain after Angiostrongylus cantonensis infection

BACKGROUND

Angiostrongyliasis is a highly dangerous infectious disease. Angiostrongylus cantonensis larvae migrate to the mouse brain and cause symptoms, such as brain swelling and bleeding. Noncoding RNAs (ncRNAs) are novel targets for the control of parasitic infections. However, the role of these molecules in A. cantonensis infection has not been fully clarified.

METHODS

In total, 32 BALB/c mice were randomly divided into four groups, and the infection groups were inoculated with 40 A. cantonensis larvae by gavage. Hematoxylin and eosin (H&E) staining and RNA library construction were performed on brain tissues from infected mice. Differential expression of long noncoding RNAs (lncRNAs) and mRNAs in brain tissues was identified by high-throughput sequencing. The pathways and functions of the differentially expressed lncRNAs were determined by Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) analyses. The functions of the differentially expressed lncRNAs were further characterized by lncRNA‒microRNA (miRNA) target interactions. The potential host lncRNAs involved in larval infection of the brain were validated by quantitative real-time polymerase chain reaction (qRT‒PCR).

RESULTS

The pathological results showed that the degree of brain tissue damage increased with the duration of infection. The transcriptome results showed that 859 lncRNAs and 1895 mRNAs were differentially expressed compared with those in the control group, and several lncRNAs were highly expressed in the middle-late stages of mouse infection. GO and KEGG pathway analyses revealed that the differentially expressed target genes were enriched mainly in immune system processes and inflammatory response, among others, and several potential regulatory networks were constructed.

CONCLUSIONS

This study revealed the expression profiles of lncRNAs in the brains of mice after infection with A. cantonensis. The lncRNAs H19, F630028O10Rik, Lockd, AI662270, AU020206, and Mexis were shown to play important roles in the infection of mice with A. cantonensis infection.

Changes in Cerebrospinal Fluid Proteins across the Spectrum of Untreated and Treated Chronic HIV-1 Infection

Using the Olink Explore 1536 platform, we measured 1,463 unique proteins in 303 cerebrospinal fluid (CSF) specimens from four clinical centers that included uninfected controls and 12 groups of people living with HIV-1 infection representing the spectrum of progressive untreated and treated chronic infection. We present three initial analyses of these measurements: an overview of the CSF protein features of the sample; correlations of the CSF proteins with CSF HIV-1 RNA and neurofilament light chain protein (NfL) concentrations; and comparison of the CSF proteins in HIV-associated dementia ( HAD ) and neurosymptomatic CSF escape ( NSE ). These reveal a complex but coherent picture of CSF protein changes that includes highest concentrations of many proteins during CNS injury in the HAD and NSE groups and variable protein changes across the course of neuroasymptomatic systemic HIV-1 progression, including two common patterns, designated as lymphoid and myeloid patterns, related to the principal involvement of their underlying inflammatory cell lineages. Antiretroviral therapy reduced CSF protein perturbations, though not always to control levels. The dataset of these CSF protein measurements, along with background clinical information, is posted online. Extended studies of this unique dataset will provide more detailed characterization of the dynamic impact of HIV-1 infection on the CSF proteome across the spectrum of HIV-1 infection, and further the mechanistic understanding of HIV-1-related CNS pathobiology.

Loss of heterozygosity impacts MHC expression on the immune microenvironment in CDK12-mutated prostate cancer

BACKGROUND

In prostate cancer (PCa), well-established biomarkers such as MSI status, TMB high, and PDL1 expression serve as reliable indicators for favorable responses to immunotherapy. Recent studies have suggested a potential association between CDK12 mutations and immunotherapy response; however, the precise mechanisms through which CDK12 mutation may influence immune response remain unclear. A plausible explanation for immune evasion in this subset of CDK12-mutated PCa may be reduced MHC expression.

RESULTS

Using genomic data of CDK12-mutated PCa from 48 primary and 10 metastatic public domain samples and a retrospective cohort of 53 low-intermediate risk primary PCa, we investigated how variation in the expression of the MHC genes affected associated downstream pathways. We classified the patients based on gene expression quartiles of MHC-related genes and categorized the tumors into "High" and "Low" expression levels. CDK12-mutated tumors with higher MHC-expressed pathways were associated with the immune system and elevated PD-L1, IDO1, and TIM3 expression. Consistent with an inflamed tumor microenvironment (TME) phenotype, digital cytometric analyses identified increased CD8 + T cells, B cells, γδ T cells, and M1 Macrophages in this group. In contrast, CDK12-mutated tumors with lower MHC expression exhibited features consistent with an immune cold TME phenotype and immunoediting. Significantly, low MHC expression was also associated with chromosome 6 loss of heterozygosity (LOH) affecting the entire HLA gene cluster. These LOH events were observed in both major clonal and minor subclonal populations of tumor cells. In our retrospective study of 53 primary PCa cases from this Institute, we found a 4% (2/53) prevalence of CDK12 mutations, with the confirmation of this defect in one tumor through Sanger sequencing. In keeping with our analysis of public domain data this tumor exhibited low MHC expression at the RNA level. More extensive studies will be required to determine whether reduced HLA expression is generally associated with primary tumors or is a specific feature of CDK12 mutated PCa.

CONCLUSIONS

These data show that analysis of CDK12 alteration, in the context of MHC expression levels, and LOH status may offer improved predictive value for outcomes in this potentially actionable genomic subgroup of PCa. In addition, these findings highlight the need to explore novel therapeutic strategies to enhance MHC expression in CDK12-defective PCa to improve immunotherapy responses.

Large-language models facilitate discovery of the molecular signatures regulating sleep and activity

Sleep, locomotor and social activities are essential animal behaviors, but their reciprocal relationships and underlying mechanisms remain poorly understood. Here, we elicit information from a cutting-edge large-language model (LLM), generative pre-trained transformer (GPT) 3.5, which interprets 10.2-13.8% of Drosophila genes known to regulate the 3 behaviors. We develop an instrument for simultaneous video tracking of multiple moving objects, and conduct a genome-wide screen. We have identified 758 fly genes that regulate sleep and activities, including mre11 which regulates sleep only in the presence of conspecifics, and NELF-B which regulates sleep regardless of whether conspecifics are present. Based on LLM-reasoning, an educated signal web is modeled for understanding of potential relationships between its components, presenting comprehensive molecular signatures that control sleep, locomotor and social activities. This LLM-aided strategy may also be helpful for addressing other complex scientific questions.

Divergent selection for feed efficiency in pigs altered the duodenum transcriptomic response to feed intake and its DNA methylation profiles

Feed efficiency is a trait of interest in pigs as it contributes to lowering the ecological and economical costs of pig production. A divergent genetic selection experiment from a Large White pig population was performed for 10 generations, leading to pig lines with relatively low- (LRFI) and high- (HRFI) residual feed intake (RFI). Feeding behavior and metabolic differences have been previously reported between the two lines. We hypothesized that part of these differences could be related to differential sensing and absorption of nutrients in the proximal intestine. We investigated the duodenum transcriptome and DNA methylation profiles comparing overnight fasting with ad libitum feeding in LRFI and HRFI pigs (n = 24). We identified 1,106 differentially expressed genes between the two lines, notably affecting pathways of the transmembrane transport activity and related to mitosis or chromosome separation. The LRFI line showed a greater transcriptomic response to feed intake than the HRFI line. Feed intake affected genes from both anabolic and catabolic pathways in the pig duodenum, such as rRNA production and autophagy. Several nutrient transporter and tight junction genes were differentially expressed between lines and/or by short-term feed intake. We also identified 409 differentially methylated regions in the duodenum mucosa between the two lines, while this epigenetic mark was less affected by feeding. Our findings highlighted that the genetic selection for feed efficiency in pigs changed the transcriptome profiles of the duodenum, and notably its response to feed intake, suggesting key roles for this proximal gut segment in mechanisms underlying feed efficiency.NEW & NOTEWORTHY The duodenum is a key organ for the hunger/satiety loop and nutrient sensing. We investigated how the duodenum transcriptome and DNA methylation profiles are affected by feed intakes in pigs. We observed thousands of changes in gene expression levels between overnight-fasted and fed pigs in high-feed efficiency pig lines, but almost none in the related low-feed efficiency pig line.

Type II collagen scaffolds repair critical-sized osteochondral defects under induced conditions of osteoarthritis in rat knee joints via inhibiting TGF-β-Smad1/5/8 signaling pathway

The bidirectional relationship between osteochondral defects (OCD) and osteoarthritis (OA), with each condition exacerbating the other, makes OCD regeneration in the presence of OA challenging. Type II collagen (Col2) is important in OCD regeneration and the management of OA, but its potential applications in cartilage tissue engineering are significantly limited. This study investigated the regeneration capacity of Col2 scaffolds in critical-sized OCDs under surgically induced OA conditions and explored the underlying mechanisms that promoted OCD regeneration. Furthermore, the repair potential of Col2 scaffolds was validated in over critical-sized OCD models. After 90 days or 150 days since scaffold implantation, complete healing was observed histologically in critical-sized OCD, evidenced by the excellent integration with surrounding native tissues. The newly formed tissue biochemically resembled adjacent natural tissue and exhibited comparable biomechanical properties. The regenerated OA tissue demonstrated lower expression of genes associated with cartilage degradation than native OA tissue but comparable expression of genes related to osteochondral anabolism compared with normal tissue. Additionally, transcriptome and proteome analysis revealed the hindrance of TGF-β-Smad1/5/8 in regenerated OA tissue. In conclusion, the engrafting of Col2 scaffolds led to the successful regeneration of critical-sized OCDs under surgically induced OA conditions by inhibiting the TGF-β-Smad1/5/8 signaling pathway.

Curculigoside, a traditional Chinese medicine monomer, ameliorates oxidative stress in Alzheimer's disease mouse model via suppressing ferroptosis

Alzheimer's disease (AD) is a neurodegenerative disorder where oxidative stress, induced by ferroptosis, has been linked to neuronal damage and cognitive deficits. The objective of this study is to investigate if the potential therapeutic agent, Curculigoside (CUR), could ameliorate AD by inhibiting ferroptosis. The potential therapeutic targets, such as GPX4 and SLC7A11, were identified using weighted gene co-expression network analysis (WGCNA). Concurrently, CUR was also screened against these potential targets using various analytical methods. For the in vivo studies, intragastric administration of CUR significantly ameliorated cognitive impairment in AD model mice induced by scopolamine and okadaic acid (OA). In vitro, CUR protected neuronal cells by altering the levels of ferroptosis-related specific markers in OA and scopolamine-induced neurotoxicity. The administration of CUR through intragastric route significantly reduced the levels of AD-promoting factors (such as Aβ1-42, p-tau) and ferroptosis-promoting factors in the hippocampus and cortex of AD mice. Furthermore, CUR up-regulated the expression of GPX4 and decreased the expression of SLC7A11 in the ferroptosis signaling pathway, thereby increasing the ratio of glutathione (GSH)/oxidized glutathione (GSSG) in vivo and vitro. In conclusion, the cumulative results suggest that the natural compound CUR may serve as a promising therapeutic agent to ameliorate AD by inhibiting ferroptosis.

Disproportionate neuroanatomical effects of DCC haploinsufficiency in adolescence compared with adulthood: links to dopamine, connectivity, covariance, and gene expression brain maps in mice

BACKGROUND

Critical adolescent neural refinement is controlled by the DCC (deleted in colorectal cancer) protein, a receptor for the netrin-1 guidance cue. We sought to describe the effects of reduced DCC on neuroanatomy in the adolescent and adult mouse brain.

METHODS

We examined neuronal connectivity, structural covariance, and molecular processes in a DCC-haploinsufficient mouse model, compared with wild-type mice, using new, custom analytical tools designed to leverage publicly available databases from the Allen Institute.

RESULTS

We included 11 DCC-haploinsufficient mice and 16 wild-type littermates. Neuroanatomical effects of DCC haploinsufficiency were more severe in adolescence than adulthood and were largely restricted to the mesocorticolimbic dopamine system. The latter finding was consistent whether we identified the regions of the mesocorticolimbic dopamine system a priori or used connectivity data from the Allen Brain Atlas to determine de novo where these dopamine axons terminated. Covariance analyses found that DCC haploinsufficiency disrupted the coordinated development of the brain regions that make up the mesocorticolimbic dopamine system. Gene expression maps pointed to molecular processes involving the expression of DCC, UNC5C (encoding DCC's co-receptor), and NTN1 (encoding its ligand, netrin-1) as underlying our structural findings.

LIMITATIONS

Our study involved a single sex (males) at only 2 ages.

CONCLUSION

The neuroanatomical phenotype of DCC haploinsufficiency described in mice parallels that observed in DCC-haploinsufficient humans. It is critical to understand the DCC-haploinsufficient mouse as a clinically relevant model system.

NANOG regulate the JAK/STAT3 pathway to promote trophoblast cell migration and epithelial-mesenchymal transition (EMT) in hypertensive disorders of pregnancy (HDP) through protein interaction with CDK1

PROBLEM

Hypertensive disorders of pregnancy (HDP) are a common pregnancy disease. NANOG and Cyclin-dependent kinase 1 (CDK1) are essential for regulating the function of cell proliferation and apoptosis. However, the mechanism of action in HDP is yet unclear.

METHOD

The microarray dataset GSE6573 was downloaded from the GEO database. Emt-related gene set was downloaded from Epithelial-Mesenchymal Transition gene database 2.0 were screened differentially expressed genes by bioinformatics analysis. Pathway Commons and Scansite 4.0 databases were used to predict the interaction between proteins. Placental tissue samples were collected from HDP patients and patients with uneventful pregnancies. RT-qPCR, Western blot and immunohistochemistry were used to detect the expression of NANOG, CDK1, MMP-2, MMP-9, EMT markers and the JAK/STAT3 pathway proteins. Transfection NANOG overexpression/knockdown, and CDK1 knockdown into the human chorionic trophoblast cells (HTR-8/Svneo). CCK-8, Transwell and Wound-healing assay were used to evaluate cell proliferation, invasion and migration. CO-IP and GST pull-down assays were used to confirm the protein interaction.

RESULTS

A total obtained seven EMT-related differentially expressed genes, wherein NANOG, NODAL and LIN28A had protein interaction. In the HDP patients' tissue found that NANOG and CDK1 had lower expression. NANOG overexpression promoted HTR-8/Svneo proliferation, migration and EMT, while NANOG knockdown had the opposite effect. Further a protein interaction between STAT3 and CDK1 with NANOG. NANOG overexpression downregulated the JAK/STAT3 pathway to promote HTR-8/Svneo proliferation, migration and EMT, which was reversed by CDK1 knockdown.

CONCLUSIONS

NANOG downregulated the JAK/STAT3 pathway to promote trophoblast cell proliferation, migration and EMT through protein interaction with CDK1.

Enhancing erythromycin production in Saccharopolyspora erythraea through rational engineering and fermentation refinement: A Design-Build-Test-Learn approach

Industrial production of bioactive compounds from actinobacteria, such as erythromycin and its derivatives, faces challenges in achieving optimal yields. To this end, the Design-Build-Test-Learn (DBTL) framework, a systematic metabolic engineering approach, was employed to enhance erythromycin production in Saccharopolyspora erythraea (S. erythraea) E3 strain. A genetically modified strain, S. erythraea E3-CymRP21-dcas9-sucC (S. erythraea CS), was developed by suppressing the sucC gene using an inducible promoter and dcas9 protein. The strain exhibited improved erythromycin synthesis, attributed to enhanced precursor synthesis and increased NADPH availability. Transcriptomic and metabolomic analyses revealed altered central carbon metabolism, amino acid metabolism, energy metabolism, and co-factor/vitamin metabolism in CS. Augmented amino acid metabolism led to nitrogen depletion, potentially causing cellular autolysis during later fermentation stages. By refining the fermentation process through ammonium sulfate supplementation, erythromycin yield reached 1125.66 mg L-1, a 43.5% increase. The results demonstrate the power of the DBTL methodology in optimizing erythromycin production, shedding light on its potential for revolutionizing antibiotic manufacturing in response to the global challenge of antibiotic resistance.

Identification of Potential Neddylation-related Key Genes in Ischemic Stroke based on Machine Learning Methods

Ischemic stroke (IS) is a complex neurological disease that can lead to severe disability or even death. Understanding the molecular mechanisms involved in the occurrence and progression of IS is of great significance for developing effective treatment strategies. In this context, the role of neddylation refers to the potential impact of neddylation on various cellular processes, which may contribute to the pathogenesis and outcome of IS. First, differential analysis was conducted on the GSE16561 dataset from the GEO database to identify 350 differentially expressed genes (DEGs) between the IS and Control groups. By intersecting the differential genes with neddylation-related genes, 11 neddylation-related DEGs were obtained. Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) analyses showed that the DEGs were mainly enriched in hematopoietic cell lineage and neutrophil degranulation, while the neddylation-related DEGs were mainly enriched in apoptosis and post-translational protein modification. Further Lasso-Cox and random forest analyses were performed on the 11 neddylation-related DEGs, identifying key genes SRPK1, BIRC2, and KLHL3. Additionally, validation of the key genes was carried out using the GSE58294 dataset and clinical patients. Finally, the correlation between the key genes and ferroptosis and cuproptosis was analyzed, and a ceRNA network was constructed. Our study helps to elucidate the complex role of neddylation in the mechanism of ischemic stroke, providing potential opportunities for the development of therapeutic interventions.

Anti-inflammatory roles of type I interferon signaling in the lung

Excessive or persistent inflammation may have detrimental effects on lung structure and function. Currently, our understanding of conserved host mechanisms that control the inflammatory response remains incompletely understood. In this study, we investigated the role of type I interferon signaling in the inflammatory response against diverse clinically relevant stimuli. Using mice deficient in type I interferon signaling (IFNAR1-/-), we demonstrate that the absence of interferon signaling resulted in a robust and persistent inflammatory response against Pseudomonas aeruginosa, lipopolysaccharide, and chemotherapeutic agent bleomycin. The elevated inflammatory response in IFNAR1-/- mice was manifested as elevated myeloid cells, such as macrophages and neutrophils, in the bronchoalveolar lavage. The inflammatory cell response in the IFNAR1-/- mice persisted to 14 days and there is impaired recovery and fibrotic remodeling of the lung in IFNAR1-/- mice after bleomycin injury. In the Pseudomonas infection model, the elevated inflammatory cell response led to improved bacterial clearance in IFNAR1-/- mice, although there was similar lung injury and survival. We performed RNA sequencing of lung tissue in wild-type and IFNAR1-/- mice after LPS and bleomycin injury. Our unbiased analysis identified differentially expressed genes between IFNAR1-/- and wild-type mice, including previously unknown regulation of nucleotide-binding oligomerization domain (NOD)-like receptor signaling, retinoic acid-inducible gene-I (RIG-I) signaling, and necroptosis pathway by type I interferon signaling in both models. These data provide novel insights into the conserved anti-inflammatory mechanisms of the type I interferon signaling.NEW & NOTEWORTHY Type I interferons are known for their antiviral activities. In this study, we demonstrate a conserved anti-inflammatory role of type I interferon signaling against diverse stimuli in the lung. We show that exacerbated inflammatory response in the absence of type I interferon signaling has both acute and chronic consequences in the lung including structural changes.

Homeobox regulator Wilms Tumour 1 is displaced by androgen receptor at cis-regulatory elements in the endometrium of PCOS patients

Decidualisation, the process whereby endometrial stromal cells undergo morphological and functional transformation in preparation for trophoblast invasion, is often disrupted in women with polycystic ovary syndrome (PCOS) resulting in complications with pregnancy and/or infertility. The transcription factor Wilms tumour suppressor 1 (WT1) is a key regulator of the decidualization process, which is reduced in patients with PCOS, a complex condition characterized by increased expression of androgen receptor in endometrial cells and high presence of circulating androgens. Using genome-wide chromatin immunoprecipitation approaches on primary human endometrial stromal cells, we identify key genes regulated by WT1 during decidualization, including homeobox transcription factors which are important for regulating cell differentiation. Furthermore, we found that AR in PCOS patients binds to the same DNA regions as WT1 in samples from healthy endometrium, suggesting dysregulation of genes important to decidualisation pathways in PCOS endometrium due to competitive binding between WT1 and AR. Integrating RNA-seq and H3K4me3 and H3K27ac ChIP-seq metadata with our WT1/AR data, we identified a number of key genes involved in immune response and angiogenesis pathways that are dysregulated in PCOS patients. This is likely due to epigenetic alterations at distal enhancer regions allowing AR to recruit cofactors such as MAGEA11, and demonstrates the consequences of AR disruption of WT1 in PCOS endometrium.