Gene Ontology Terms Research Articles

Dietary Sugar and Saturated Fat Consumption Associated with the Gastrointestinal Microbiome during Pregnancy

BackgroundGrowing evidence supports changes in the gastrointestinal microbiome over the course of pregnancy may have an impact on the short- and long-term health of both the mother and the child. ObjectiveOur objective was to explore the association of diet quality, as measured by the Healthy Eating Index (HEI), with the composition and gene ontology (GO) representation of microbial function in the maternal gastrointestinal microbiome during pregnancy. MethodsWe conducted a retrospective, observational analysis of n = 185 pregnant participants in the Pregnancy Eating Attributes Study. Maternal dietary intake was assessed in the first trimester using the automated self-administered 24-h recall method, from which the HEI 2015 was calculated. Rectal swabs were obtained in the second trimester and sequenced using the NovaSeq 6000 system shotgun platform. We used unsupervised clustering to identify microbial enterotypes representative of maternal taxa and GO functional term composition. Multivariable linear models were used to identify associations between taxa, functional terms, and food components while controlling for relevant covariates. Multinomial regression was then used to predict enterotype membership based on a participant’s HEI food component score. ResultsThose in the high diet quality tertile had a lower early pregnancy BMI [mean (M) = 23.48 kg/m2, SD = 3.38] compared with the middle (M = 27.35, SD = 6.01) and low (M = 27.49, SD = 6.99) diet quality tertiles (P < 0.01). There were no statistically significant associations between the HEI components or the total HEI score and the 4 alpha diversity measures. Differences in taxa and GO term enterotypes were found in participants with, but not limited to, a higher saturated fat component score (β = 1.35, P = 0.01), added sugar HEI component (β = 0.07, P < 0.001), and higher total dairy score (β = 1.58, P = 0.01). ConclusionsSpecific dietary components are associated with microbial composition and function in the second trimester of pregnancy. These findings provide a foundation for future testable hypotheses.

Analysis of the Long Non-Coding and Messenger RNA Expression Profiles in the Skin Tissue of Super Merino and Small-Tailed Han Sheep.

Wool quality and yield are two important economic livestock traits. However, there are relatively few molecular studies on lncRNA for improving sheep wool, so these require further exploration. In this study, we examined skin tissue from the upper scapula of Super Merino (SM) and Small-Tailed Han (STH) sheep during the growing period. The apparent difference was verified via histological examination. High-throughput RNA sequencing identified differentially expressed (DE) long non-coding (lncRNAs) and messenger RNAs (mRNAs). The target gene of DE lncRNA and DE genes were enrichment analyzed using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG). A Reverse Transcription quantitative Polymerase Chain Reaction (RT-qPCR) was used to verify randomly selected DE lncRNAs and mRNAs. Finally, the DE, RAC2, WNT11, and FZD2 genes, which were enriched in the Wnt signaling pathway, were detected via immunohistochemistry. The results showed that a total of 20,888 lncRNAs and 31,579 mRNAs were identified in the skin tissues of the two sheep species. Among these, 56 lncRNAs and 616 mRNAs were differentially expressed. Through qRT-PCR, the trends in the randomly selected DE genes' expression were confirmed to be aligned with the RNA-seq results. GO and KEGG enrichment analysis showed that DE lncRNA target genes were enriched in GO terms as represented by epidermal and skin development and keratin filature and in KEGG terms as represented by PI3K-Akt, Ras, MAPK, and Wnt signaling pathways, which were related to hair follicle growth and development. Finally, immunohistochemistry staining results indicated that RAC2, WNT11, and FZD2 were expressed in dermal papilla (DP). The lncRNAs MSTRG.9225.1 and MSTRG.98769.1 may indirectly participate in the regulation of hair follicle growth, development, and fiber traits by regulating their respective target genes, LOC114113396(KRTAP15-1), FGF1, and IGF1. In addition, MSTRG.84658.1 may regulate the Wnt signaling pathway involved in the development of sheep hair follicles by targeting RAC2. This study provides a theoretical reference for improving sheep breeding in the future and lays a foundation for further research on the effects of MSTRG.84658.1 and the target gene RAC2 on dermal papilla cells (DPC).

Multi-omics approach for understanding the response of Bacteroides fragilis to carbapenems

BackgroundThe prevalence of Bacteroides fragilis isolates resistant to first-line beta-lactam drugs is increasing, resulting in reduced treatment efficacy. Investigating the bacterial transcriptome and proteome can uncover links between bacterial genes and resistance mechanisms. In this study, we experimentally assessed in vitro the transcriptional and proteomic profiles of B. fragilis exposed to SICs of meropenem, an effective antimicrobial agent, collected from patients with intra-abdominal diseases at Astana City Hospital, Kazakhstan. MethodsB. fragilis was cultured in brain heart infusion broth and sub-cultured every 48 h for 8 days in media with and without meropenem. Total RNA was extracted from the liquid cultures using a commercial RNeasy mini kit, and strand-specific RNA sequencing (RNA-seq) was performed on the DNBSEQ platform. Raw RNA-seq data were retrieved from BioProject No. PRJNA531645 and uploaded to the NCBI Sequence Read Archive (accession no. SRX22081155). Proteins of B. fragilis were extracted and separated using sodium dodecyl sulphate–polyacrylamide gel electrophoresis, followed by analysis of the eluted peptides using liquid chromatography–tandem mass spectrometry. Cluster analysis utilised the Database for Annotation, Visualisation, and Integrated Discovery. ResultsThe subinhibitory concentration (SIC) of meropenem was determined to be 0.5 μg/L (minimum inhibitory concentration: 1). Mapping of reads to the reference genome identified 2477 expressed genes in all B. fragilis BFR KZ01 samples. Ten differentially expressed genes (DEGs) were common across comparison groups during and post-antibiotic exposure (wMEM vs. MEM2 and MEM2 vs. rMEM8); however, no substantially enriched Gene Ontology terms were identified. The cluster analysis highlighted a significant enrichment cluster (W-0560 oxidoreductase) of DEGs following antibiotic withdrawal. In total, 859 B. fragilis proteins were identified, with the expressions of three proteins, 3-oxoacyl-[acyl carrier protein] reductase, acetyl-CoA carboxylase biotin carboxylase subunit, and beta-ketoacyl-ACP synthase III, being upregulated in the enriched protein folding category. Notably, chaperone proteins such as FKBP-type peptidyl-prolyl cis-trans isomerases (involved in the cis-trans isomerisation of prolyl peptide bonds) and GroES (a co-chaperone functioning with GroEL) were also identified. ConclusionsUnder the influence of low doses of antibiotics defense mechanisms are activated which contribute to the emergence of resistance. These results provide insight into the response of B. fragilis to meropenem exposure, mainly at the SIC, contributing to the understanding bacterial survival strategies under stress conditions.

Unveiling circRNA-mediated ceRNA networks in ischemic stroke by integrative analysis of multi-source gene expression profiling

BackgroundCircular RNAs (circRNAs) are emerging as potential therapeutic targets for ischemic stroke (IS) due to their regulatory roles in inflammation and apoptosis. This study aimed to develop a comprehensive and robust IS-specific competing endogenous RNA (ceRNA) network to facilitate the identification of novel diagnostic and therapeutic targets. MethodsWe integrated expression data from 15 IS studies using the Rank-In algorithm to minimize batch effects. Differentially expressed circRNAs, miRNAs, and mRNAs were identified by comparing IS and control samples. Functional enrichment analysis of differentially expressed circRNA host genes revealed significantly enriched pathways and Gene Ontology (GO) terms relevant to IS pathogenesis. We further predicted miRNA-circRNA and mRNA-miRNA interactions, enabling the construction of a comprehensive ceRNA network to identify circRNA-related genes with diagnostic potential for IS. ResultsIntegrated analysis revealed 199 differentially expressed circRNAs, 103 miRNAs, and 1736 mRNAs in IS patients. Functional enrichment analysis implicated these molecules in relevant pathways like the neurotrophin signaling pathway and p53 signaling pathway. The constructed circRNA-miRNA-mRNA regulatory network provided insights into potential mechanisms underlying IS. Three circRNA-related genes (RGS2, CDK5R1, and NSF) displayed promising diagnostic potential for IS when combined. ConclusionsWe successfully constructed a robust and informative IS-specific ceRNA network by integrating data from diverse sources. This network identified differentially expressed RNAs and revealed enriched pathways potentially involved in IS pathogenesis. Notably, our analysis identified CDK5R1, RGS2, and NSF as potential diagnostic biomarkers for IS. This study sheds light on a circRNA-mediated regulatory network with potential diagnostic and therapeutic implications for ischemic stroke.

A genome-wide methylation analysis of Chinese Han patients with chronic insomnia disorder.

As the most common sleep disorder, chronic insomnia disorder (CID) has become a global health burden to the public. However, it remains unclear about the pathogenesis of this disease. Epigenetic changes may provide important insights into the gene-environment interaction in CID. Therefore, this study was conducted to investigate the DNA methylation pattern in CID and reveal the epigenetic mechanism of this disease. In this study, whole blood DNA was extracted from 8 CID patients (the CID group) and 8 healthy controls (the control group), respectively. Besides, genome-wide DNA methylation was detected by Illumina Human Methylation 850K Beadchip. Moreover, the sleep quality and insomnia severity were evaluated by the Pittsburgh Sleep Quality Index (PSQI) and Insomnia Severity Index (ISI), respectively. A total of 369 differentially methylated positions (DMPs) and 23 differentially methylated regions (DMRs) were identified between the CID and control groups. LHX6 was identified as the most important differentially methylated gene (DMG). The Gene Ontology (GO) analysis results corroborated that DMPs were significantly enriched in 105 GO terms, including cell signaling, homogenous cell adhesion of plasma membrane adhesion molecules, nervous system development, cell adhesion, and calcium ion binding. In addition, it was demonstrated that DMPs were significantly enriched in the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, including the hippo signaling pathway, Ras signaling pathway, and vitamin B6 metabolism. The DMR-related GO analysis results revealed the positive regulation of protein kinase activities. DNA methylation plays a critical role in the development of CID, and LHX6 is validated to be an important DMG.

Cashmere cyclic growth affected by different photoperiods alters DNA methylation patterns

Whole-genome bisulfite sequencing was used to study epigenetic differences in hair follicles from Inner Mongolia Albas cashmere goats during long and short photoperiods. Results showed earlier anagen entry in short photoperiod. Genome-wide cytosine methylation was 3.23%, with CG, CHG, and CHH methylation rates at 48.00%, 0.19%, and 0.17%, respectively. Methylation was higher for CG than for CHG and CHH, and higher for introns than exons and downstream regions. Both groups had similar levels of CG, CHG, and CHH methylation, but CHH methylation levels differed slightly between photoperiods. A total of 9,700 differentially methylated regions (DMRs) and 4,187 genes (DMGs) were identified. Gene Ontology (GO) and KEGG analyses revealed significant enrichment in 67 GO terms and 41 signaling pathways related to circadian rhythms, bringing new insights in the role of DNA methylation in cashmere growth and development.

Altered transcriptomes, cell type proportions, and dendritic spine morphology in hippocampus of suicide decedents

BackgroundSuicide is a manner of death resulting from complex environmental and genetic risks that affect millions of people globally. Both structural and functional studies identified the hippocampus as one of the vulnerable brain regions contributing to suicide risk. MethodsWe have identified the hippocampal tissue transcriptomes, gene ontology, cell type proportions, and dendritic spine morphology in controls (n = 28) and suicide decedents (n = 22). In addition, the transcriptomic signature in iPSC-derived neuronal precursor cells (NPCs) and neurons were also investigated in controls (n = 2) and suicide decedents (n = 2). ResultsThe hippocampal tissue transcriptomic data revealed that NPAS4 gene expression was downregulated while ALDH1A2, NAAA, and MLXIPL gene expressions were upregulated in hippocampal tissue of suicide decedents. The gene ontology identified 29 significant pathways including NPAS4-associated gene ontology terms “excitatory post-synaptic potential”, “regulation of postsynaptic membrane potential” and “long-term memory” indicating alteration of glutamatergic synapses in the hippocampus of suicide decedents. The cell type deconvolution identified decreased excitatory neuron proportion and an increased inhibitory neuron proportion providing evidence of excitation/inhibition imbalance in the hippocampus of suicide decedents. In addition, suicide decedents had increased dendric spine density in the hippocampus, due to an increase of thin (relatively unstable) dendritic spines, compared to controls. The transcriptomes of iPSC-derived hippocampal-like NPCs and neurons revealed 31 and 33 differentially expressed genes in NPC and neurons, respectively, of suicide decedents. ConclusionsOur findings will provide new insights into the hippocampal neuropathology of suicide.

Upregulated solute-carrier family genes in the hippocampus of schizophrenia can be rescued by antipsychotic medications

Background and hypothesisOur previous studies have found that functional changes in the hippocampal circuit from dentate gyrus (DG) to cornu ammonis 3 and 1 (CA3, CA1) are highly associated with schizophrenia (SZ). However, no studies have explored the genetic expression across the three and two human hippocampal subfields (DG-CA3-CA1 and CA3-CA1) between subjects with SZ and healthy controls (CT). Study designWe matched cohorts between CT (n = 13) and SZ (n = 13). Among SZ, 6 subjects were on antipsychotics (AP) while 7 were off AP. We combined RNA-seq data from all three and two hippocampal subfields and performed differentially expressed gene analyses across DG-CA3-CA1 and CA3-CA1 affected by either SZ or AP. Study resultsWe found that differentially expressed genes (DEGs) from effects of SZ and AP across DG-CA3-CA1 and CA3-CA1 were highly associated with gene ontology terms related to hormonal and immune signaling, cellular mitosis and apoptosis, ion and amino acid transports, and protein modification and degradation. Additionally, we found that multiple genes related to solute-carrier family and immune signaling were significantly upregulated across DG-CA3-CA1 and CA3-CA1 in patients with SZ relative to CT, and AP consistently and robustly repressed the expression of these upregulated genes in the DG-CA3-CA1 and CA3-CA1 from subjects with SZ. ConclusionsTogether, these data suggest that the upregulated solute-carrier family genes in the hippocampus might have important roles in the pathophysiology of SZ, and that AP may reduce the symptoms of psychosis in SZ via rescuing the solute-carrier gene expression.

Randomized Crossover Clinical Trial of Nicotinamide Riboside and Coenzyme Q10 on Metabolic Health and Mitochondrial Bioenergetics in CKD.

Mitochondria-driven oxidative/redox stress and inflammation play a major role in chronic kidney disease (CKD) pathophysiology. Compounds targeting mitochondrial metabolism may improve mitochondrial function, inflammation, and redox stress; however, there is limited evidence of their efficacy in CKD. We conducted a randomized, double-blind, placebo-controlled crossover trial comparing the effects of 1200 mg/day of coenzyme Q10 (CoQ10) or 1000 mg/day of nicotinamide riboside (NR) supplementation to placebo in 25 people with moderate-to-severe CKD (eGFR <60mL/min/1.73 m2). We assessed changes in the blood transcriptome using 3'-Tag-Seq gene expression profiling and changes in pre-specified secondary outcomes of inflammatory and oxidative stress biomarkers. For a subsample of participants (n=14), we assessed lymphocyte and monocyte bioenergetics using an extracellular flux analyzer. The (mean±SD) age, eGFR, and BMI of the participants were 61±11 years, 37±9 mL/min/1.73m2, and 28±5 kg/m2 respectively. Of the participants, 16% had diabetes and 40% were female. Compared to placebo, NR-mediated transcriptomic changes were enriched in gene ontology (GO) terms associated with carbohydrate/lipid metabolism and immune signaling while, CoQ10 changes were enriched in immune/stress response and lipid metabolism GO terms. NR increased plasma IL-2 (estimated difference, 0.32, 95% CI of 0.14 to 0.49 pg/mL), and CoQ10 decreased both IL-13 (estimated difference, -0.12, 95% CI of -0.24 to -0.01 pg/mL) and CRP (estimated difference, -0.11, 95% CI of -0.22 to 0.00 mg/dL) compared to placebo. Both NR and CoQ10 reduced 5 series F2-Isoprostanes (estimated difference, -0.16 and -0.11 pg/mL, respectively; P<0.05 for both). NR, but not CoQ10, increased the bioenergetic health index (BHI) (estimated difference, 0.29, 95% CI of 0.06 to 0.53) and spare respiratory capacity (estimated difference, 3.52, 95% CI of 0.04 to 7 pmol/min/10,000 cells) in monocytes. Six weeks of NR and CoQ10 improved in oxidative stress, inflammation, and cell bioenergetics in persons with moderate to severe CKD.

Exploring the Biomarkers and Potential Mechanisms of Botulinum Toxin Type A in the Treatment of Microglial Inflammatory Activation through P2X7 Receptors based on Transcriptome Sequencing.

This study aims to explore the potential mechanism by which Botulinum toxin type A (BoNT/ A) inhibits microglial inflammatory activation through P2X7 receptors (P2X7R). BoNT/A is a promising analgesic drug, and previous studies have established that it alleviates Neuropathic Pain (NP) by inhibiting microglial inflammatory activation. This study examined the biomarkers and potential mechanisms by which BoNT/A relieves neuropathic pain by mediating microglial P2X7R and analyzing transcriptome sequencing data from mouse BV-2 microglial cells. The P2X7R agonist Bz-ATP was used to induce microglial inflammatory activation, whilst RNAseq technology was used to explore the biomarkers and potential mechanisms through which BoNT/A suppresses microglial inflammation. RNA sequencing was performed on three BV-2 cell samples treated with a P2X7R specific activator (Bz-ATP) and three BV-2 cell samples pre-treated with BoNT/A. Only data that successfully passed quality control measures were included in subsequent analysis. Initially, Differentially Expressed Genes (DEGs) were identified from BoNT/A and control samples, followed by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. Biomarkers were then identified by constructing a Protein- Protein Interaction (PPI) network and utilizing the CytoHubba plug-in in Cytoscape software. Lastly, enrichment analysis and regulatory network analysis were performed to elucidate the potential mechanism of BoNT/A in the treatment of NP. 93 DEGs related to the "cell component size regulation" GO term and enriched in the "axon guidance" KEGG pathway were identified. Subsequently, 6 biomarkers were identified, namely PTPRF, CHDH, CKM, Ky, Sema3b, and Sema3f, which were enriched in pathways related to biosynthesis and metabolism, disease progression, signal transduction, and organelle function, including the "ribosome" and "Wnt signaling pathway." Finally, a competing endogenous RNA (ceRNAs) network was constructed from 6 mRNAs, 66 miRNAs, and 31 lncRNAs, forming a complex relationship network. Six genes (PTPRF, Sema3b, Sema3f, CHDH, CKM, and Ky) were identified as biomarkers of microglial inflammatory activation following BoNT/A treatment. This finding may provide a valuable reference for the relief and treatment of neuropathic pain.

Comparative profiling of whole-cell and exosome samples reveals protein signatures that stratify breast cancer subtypes

Identifying novel breast cancer biomarkers will improve patient stratification, enhance therapeutic outcomes, and help develop non-invasive diagnostics. We compared the proteomic profiles of whole-cell and exosomal samples of representative breast cancer cell subtypes to evaluate the potential of extracellular vesicles as non-invasive disease biomarkers in liquid biopsies. Overall, differentially-expressed proteins in whole-cell and exosome samples (which included markers for invasion, metastasis, angiogenesis, and drug resistance) effectively discriminated subtypes; furthermore, our results confirmed that the proteomic profile of exosomes reflects breast cancer cell-of-origin, which underscores their potential as disease biomarkers. Our study will contribute to identifying biomarkers that support breast cancer patient stratification and developing novel therapeutic strategies. We include an open, interactive web tool to explore the data as a molecular resource that can explain the role of these protein signatures in breast cancer classification.Graphical (A) We quantified proteomic profiles of breast cancer cells (BCCs) and breast cancer-derived exosomes (BCDEs) samples from four breast cancer cell lines representative of common breast cancer subtypes – Her2-positive (Her2, MDA-MB-453), Luminal A (LA, MCF7), Luminal B (LB, ZR-75), and triple-negative (TN, MDA-MB-231). (B) We independently performed four comparisons for BCCs and BCDEs, comparing each cell line against the remaining three. (C) We identified differentially-expressed proteins (DEPs) for BCCs and BCDEs, defined protein signatures, and functionally analyzed resulting networks. Through pathway inference analysis (PIA), we identified Kyoto Encyclopedia of genes and genomes (KEGG) subpaths and biological Gene Ontology (GO) terms that displayed differential activation. (D) We validated our proteomic signature using the Cancer Genome Atlas (TCGA) and the Cancer Proteome Atlas (TCPA) databases, verified that differentially-activated pathways in BCDEs caused a corresponding response in receptor cells, and confirmed that the BCDE proteomic signature reflects their cell-of-origin and identifies candidate disease biomarkers in liquid biopsies

MicroRNA expression profiles reveal wool development and fineness regulation in Gansu alpine fine-wool sheep

The development of wool has a complex regulatory mechanism both influenced by genetic and environmental factors. MicroRNAs (miRNA) were involved in various biological processes of animals, and may play an important role in the regulation of wool development. In this study, we comprehensively analyzed and identified the histological parameters of hair follicles, as well as the miRNAs, target genes, pathways, and Gene Ontology terms related to wool fineness regulation and wool growth and development using HE staining and RNA-Seqs methods. Both coarse (group C, mean fiber diameter (MFD) = 22.26 ± 0.69 μm, n = 6) and fine (group F, MFD = 16.91 ± 0.29 μm, n = 6) of Gansu alpine fine-wool sheep with different wool fineness were used in this study. The results showed that the primary follicle diameter and secondary wool fiber diameter in group C were significantly higher than those in group F (P < 0.05). And the number of primary and secondary hair follicles in group C was significantly lower than that in group F (P < 0.05). Furthermore, a total of 67 DE miRNAs and 290 potential DE miRNAs target genes were screened in the skin tissues of sheep from groups F and C, and some potential target genes related to wool fineness regulation were screened, such as CDH2, KRT82, FOXN1, LOC101106296, KRT20, MCOLN3, KRT71, and TERT. These genes were closely related to Glutathione metabolism, epidermal cell differentiation, keratinization, and regulation of hair cycle. Moreover, the regulatory network of miRNAs-mRNAs suggested that miRNAs (miR-129-x, novel m0079-3p, miR-2484-z, novel m0025-5P, etc.) may play a key role in the wool development and wool fineness regulation of Gansu alpine fine-wool sheep. In summary, this study expands the existing miRNAs database and provides new information for studying the regulation of wool development in Gansu alpine fine wool sheep.

Bioinformatics analysis to disclose shared molecular mechanisms between type-2 diabetes and clear-cell renal-cell carcinoma, and therapeutic indications

Type 2 diabetes (T2D) and Clear-cell renal cell carcinoma (ccRCC) are both complicated diseases which incidence rates gradually increasing. Population based studies show that severity of ccRCC might be associated with T2D. However, so far, no researcher yet investigated about the molecular mechanisms of their association. This study explored T2D and ccRCC causing shared key genes (sKGs) from multiple transcriptomics profiles to investigate their common pathogenetic processes and associated drug molecules. We identified 259 shared differentially expressed genes (sDEGs) that can separate both T2D and ccRCC patients from control samples. Local correlation analysis based on the expressions of sDEGs indicated significant association between T2D and ccRCC. Then ten sDEGs (CDC42, SCARB1, GOT2, CXCL8, FN1, IL1B, JUN, TLR2, TLR4, and VIM) were selected as the sKGs through the protein–protein interaction (PPI) network analysis. These sKGs were found significantly associated with different CpG sites of DNA methylation that might be the cause of ccRCC. The sKGs-set enrichment analysis with Gene Ontology (GO) terms and KEGG pathways revealed some crucial shared molecular functions, biological process, cellular components and KEGG pathways that might be associated with development of both T2D and ccRCC. The regulatory network analysis of sKGs identified six post-transcriptional regulators (hsa-mir-93-5p, hsa-mir-203a-3p, hsa-mir-204-5p, hsa-mir-335-5p, hsa-mir-26b-5p, and hsa-mir-1-3p) and five transcriptional regulators (YY1, FOXL1, FOXC1, NR2F1 and GATA2) of sKGs. Finally, sKGs-guided top-ranked three repurposable drug molecules (Digoxin, Imatinib, and Dovitinib) were recommended as the common treatment for both T2D and ccRCC by molecular docking and ADME/T analysis. Therefore, the results of this study may be useful for diagnosis and therapies of ccRCC patients who are also suffering from T2D.

Improving protein function prediction by learning and integrating representations of protein sequences and function labels.

As fewer than 1% of proteins have protein function information determined experimentally, computationally predicting the function of proteins is critical for obtaining functional information for most proteins and has been a major challenge in protein bioinformatics. Despite the significant progress made in protein function prediction by the community in the last decade, the general accuracy of protein function prediction is still not high, particularly for rare function terms associated with few proteins in the protein function annotation database such as the UniProt. We introduce TransFew, a new transformer model, to learn the representations of both protein sequences and function labels [Gene Ontology (GO) terms] to predict the function of proteins. TransFew leverages a large pre-trained protein language model (ESM2-t48) to learn function-relevant representations of proteins from raw protein sequences and uses a biological natural language model (BioBert) and a graph convolutional neural network-based autoencoder to generate semantic representations of GO terms from their textual definition and hierarchical relationships, which are combined together to predict protein function via the cross-attention. Integrating the protein sequence and label representations not only enhances overall function prediction accuracy, but delivers a robust performance of predicting rare function terms with limited annotations by facilitating annotation transfer between GO terms. https://github.com/BioinfoMachineLearning/TransFew.

Host and venom evolution in parasitoid wasps: does independently adapting to the same host shape the evolution of the venom gland transcriptome?

BackgroundVenoms have repeatedly evolved over 100 occasions throughout the animal tree of life, making them excellent systems for exploring convergent evolutionary novelty. Growing evidence supports that venom evolution is predominantly driven by prey or host-related selection pressures, and the expression patterns of venom glands reflect adaptive evolution. However, it remains elusive whether the evolution of expression patterns in venom glands is likewise a convergent evolution driven by their prey/host species.ResultsWe utilized parasitoid wasps that had independently adapted to Drosophila hosts as models to investigate the convergent evolution of venom gland transcriptomes in 19 hymenopteran species spanning ~ 200 million years of evolution. Comparative transcriptome analysis reveals that the global expression patterns among the venom glands of Drosophila parasitoid wasps do not achieve higher similarity compared to non-Drosophila parasitoid wasps. Further evolutionary analyses of expression patterns at the single gene, orthogroup, and Gene Ontology (GO) term levels indicate that some orthogroups/GO terms show correlation with the Drosophila parasitoid wasps. However, these groups rarely include genes highly expressed in venom glands or putative venom genes in the Drosophila parasitoid wasps.ConclusionsOur study suggests that convergent evolution may not play a predominant force shaping gene expression levels in the venom gland of the Drosophila parasitoid wasps, offering novel insights into the co-evolution between venom and prey/host.

Differential proteomic expression in non-functional pituitary neuroendocrine tumours and pituitary glands.

Pituitary neuroendocrine tumours (PitNETs) are common accounting for 10 to 25 % of all intracranial tumours. This project describes the feasibility of developing a novel membrane-based biomarker that could be used for fluorescent guided surgery. The aim was to catalogue the differential expression of membrane proteins between non-functional PitNETs and pituitary glands. Ten pituitary gland tissue specimens were obtained from the National Institute of Health (NIH) NeuroBio-Bank and twenty non-functional PitNETs were obtained from the Northwestern University Nervous System Tumour Bank. Mass spectrometry analysis using an Orbitrap Fusion Lumos Tribrid Mass Spectrometer linked to a Dionex Ultimate 3000 UPLC system was undertaken. Data Dependent Acquisition Mass Spectrometry and Data Independent Acquisition Mass Spectrometry was then completed. Pathway enrichment analysis was performed using clusterProfiler v4.6.0. Functional enrichment analysis was conducted using Gene Ontology terms and Reactome pathways. Differential expression analysis between the two groups revealed a total of 2110 significant differently expressed proteins (DEPs), with 1387 of these also having a Log2 fold change either greater than 1, or less than -1. Of the 2110 DEPs, 925 were upregulated in tumours compared to control, while 1185 were downregulated. We have demonstrated a proteomic comparison between non-functional PitNETs and normal pituitary glands. These results demonstrate differences consistent with contemporary literature but shows that NOTCH3 and PTPRJ are up-regulated in non-functional PitNETs compared to pituitary glands.

Integrative analysis of transcriptome-wide association study and mRNA expression profile identified risk genes for bipolar disorder

ObjectiveBipolar disorder (BD) is a debilitating neuropsychiatric disorder, which is associated with genetic variation through “vast but mixed” Genome-Wide Association Studies (GWAS). Transcriptome-Wide Association Study (TWAS) is more effective in explaining genetic factors that influence complex diseases and can help identifying risk genes more reliably. So, this study aims to identify potential BD risk genes in pedigrees with TWAS. MethodsWe conducted a TWAS analysis with expression quantitative trait loci (eQTL) analysis on extended BD pedigrees, and the BD genome-wide association study (GWAS) summary data acquired from the Psychiatric Genomics Consortium (PGC). Furthermore, the BD-associated genes identified by TWAS were validated by mRNA expression profiles from the Gene Expression Omnibus (GEO) Datasets (GSE23848 and GSE46416). Functional enrichment and annotation analysis were implemented by RStudio (version 4.2.0). ResultsTWAS identified 362 genes with P value < 0.05, and 18 genes remain significant after Bonferroni correction, such as SEMA3G (PTWAS=1.07 × 10-11), ALOX5AP (PTWAS=3.12 × 10-8), and PLEC (PTWAS=1.27 × 10-7). Further 6 overlapped genes were detected in integrative analysis, such as UQCRB (PTWAS=0.0020, PmRNA=0.0000), TMPRSS9 (PTWAS=0.0405, PmRNA=0.0032), and SNX10 (PTWAS=0.0104, PmRNA=0.0015). Using genes identified by TWAS, Gene Ontology (GO) enrichment analysis identified 40 significant GO terms, such as mitochondrial ATP synthesis coupled electron transport, mitochondrial respiratory, aerobic electron transport chain, oxidative phosphorylation, mitochondrial membrane proteins, and ubiquinone activity. The Kyoto Encyclopedia of Genes and Genomes (KEGG) Pathway enrichment analysis identified significant 15 pathways for BD, such as Oxidative phosphorylation, endocannabinoids signaling, neurodegeneration, and reactive oxide species. ConclusionsWe found a set of BD-associated genes and pathways, validating the important role of neurodevelopmental abnormalities, inflammatory responses, and mitochondrial dysfunction in the pathology of BD, offering novel information for comprehending the genetic basis of BD.

Transcriptomic Insights into the Physiological Aspects of the Saprotrophic Fungus Penicillium citrinum During the Spoilage of Tobacco Leaves

Penicillium citrinum is one of the most prevalent tobacco spoilage fungi. However, the mechanisms underlying fungal growth on tobacco leaves remain largely unknown. In this study, transcriptomic analyses were performed to reveal the genome-wide expression profiles of P. citrinum growing on tobacco leaves. First, a comparative analysis was conducted between two sets of transcriptomic data from P. citrinum growing on chemically defined media and tobacco leaves. Enrichment analyses showed that differentially regulated genes were mainly associated with carbohydrate degradation (e.g., cellulose, pectin, and xylan) and the catabolism of fatty acids and aromatic compounds. Comparative transcriptomic analyses between different time points indicated that the fungal transcriptome varied dynamically during the spoilage process, and the enriched terms were associated with small-molecule degradation and fungal development. Enrichment analyses indicated that more up-regulated genes appeared in all enriched Gene Ontology terms. Notably, more organelles significantly contributed to further fungal growth on tobacco leaves. In conclusion, P. citrinum activates a comprehensive transcriptome that changes dynamically when causing tobacco mildew.

Changes in gene expression due to aging in the hypothalamus of mice.

Aging generally affects food consumption and energy metabolism. Since the feeding center is located in the hypothalamus, it is a major target for understanding the mechanism of age-related changes in eating behavior and metabolism. To obtain insight into the age-related changes in gene expression in the hypothalamus, we investigated genes whose expression changes with age in the hypothalamus. A DNA microanalysis was performed using hypothalamus samples obtained from young (aged 24 weeks) and old male mice (aged 138 weeks). Gene Ontology (GO) analysis was performed using the identified differentially expressed genes. We observed that the expression of 377 probe sets was significantly altered with aging (177 were upregulated and 200 were downregulated in old mice). As a result of the GO analysis of these probe sets, 16 GO terms, including the neuropeptide signaling pathway, were obtained. Intriguingly, although the food intake in old mice was lower than that in young mice, we found that several neuropeptide genes, such as agouti-related neuropeptide (Agrp), neuropeptide Y (Npy), and pro-melanin-concentrating hormone (Pmch), all of which promote food intake, were upregulated in old mice. In conclusion, this suggests that the gene expression pattern in the hypothalamus is regulated to promote food intake.

ITRAQ-Based Serum Proteomic Analysis Reveals Multifactorial Cellular Function Impairment and Aggravated Systematic Inflammation in Drug-free Obsessive-Compulsive Disorders.

Obsessive-compulsive disorder (OCD) is a debilitating mental disorder with obvious difficulties in treatment. Its pathogenesis has not been fully elucidated. Further understanding of etiology and mechanism needs to be explored further. We employed the isobaric tag for relative and absolute quantitation (iTRAQ)-based proteomic analysis to compare serum proteome profile between OCD patients and healthy controls, in order to find out the possible mechanism of OCD in the downstream biological process. Eighty-one drug-free OCD patients and 78 healthy controls were enrolled. A total of 475 proteins were identified. Totally, 80 proteins with p < 0.05 were selected for gene set enrichment analysis (GSEA), and only those with a fold change ≥1.2 and q value <0.2 between groups were accepted as differentially expressed proteins (DEPs). We observed a significant enrichment of immuno-inflammation-related pathways, along with intriguing expression trends that immuno-inflammation-related proteins were upregulated in GSEA. After that, 2 up-regulated proteins and 13 down-regulated ones were accepted as DEP. According to the available literature, most of the DEPs have not been reported in OCD. These DEPs were enriched in 121 gene ontology (GO) terms, including hepatocyte growth factor receptor activity, angiogenin-PRI complex, and so on. DEPs were enriched in pathways including adherens junction in the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. Alterations in DEPs including STXBP5L, GRN, and ANG were validated in OCD animal models. Our study suggested that OCD patients manifested multifactorial impairment in neuronal or non-neuronal cellular function under the inflammatory background. Further research employing larger sample sizes, longitudinal design, stratified analysis, and multiomics methodology will be needed. Experiments in laboratories were essential in illuminating the mechanism.