Gene Enrichment Analysis Research Articles

#1984 Extracellular vesicles drive loss of endothelial homeostasis in pediatric chronic kidney disease

Abstract Background and Aims Cardiovascular disease (CVD) is the main cause of mortality in chronic kidney disease (CKD). However, the pathogenesis of CVD in CKD remains incompletely understood. Extracellular vesicles (EVs) emerged as mediators of inter-organ cross-talk and endothelial EVs (EC-EVs) were previously associated with CVD. We hypothesized that CKD drives endothelial EV release in CKD and that these EC-EVs drive CVD in CKD. Method A cohort of 94 children (mean age 10.9 years) at different stages of CKD (patients with or without dialysis and after kidney transplantation (KTx)) and age-matched healthy donors was recruited, offering the unique opportunity to analyze cardiovascular effects of CKD and EV dynamics in absence of age-related comorbidities. Thorough characterization of plasma EVs by electron microscopy, nanoparticle tracking analysis (NTA), flow cytometry and miRNA sequencing was carried out. Functional properties of patient EVs and CKD-specific EV release mechanisms were tested in vitro. Results EC-EVs were elevated in hemodialysis (HD) patients (as measured by flow cytometry) while total EV counts showed higher concentrations in peritoneal dialysis patients (according to NTA). Longitudinal follow-ups revealed that EC-EVs decreased after receiving KTx. EV small RNA sequencing showed lower abundance of 30 microRNAs in EVs from CKD patients, most pronounced in HD patients compared to both healthy donors and KTx recipients. Downregulation of let-7d-5p, miR-19a-3p, miR-24-3p, miR-103a-3p and miR-142-3p was confirmed using RT-qPCR. Gene enrichment analyses predicted angiogenesis and smooth muscle cell (SMC) proliferation as targets of differentially regulated miRNAs in CKD. In vitro, angiogenesis of Human Umbilical Vein Endothelial Cells (HUVECs) was decreased upon treatment with CKD EVs as indicated by decreased vascular tube formation, endothelial migration and proliferation while SMC proliferation remained unchanged. To assess the putative role of uremic toxins for EV release in CKD, we focused on tryptophan metabolites and performed targeted plasma metabolomics that showed stage-dependent increases of uremic toxins of indole and kynurenine pathways in children with CKD, with highest levels in dialysis patients and almost normal levels after KTx. Uremic toxins indoxyl sulfate, formylkynurenine and xanthurenic acid correlated positively with EC-EVs in our cohort, but were not sufficient to trigger EV release from aortic endothelial cells alone. However, high shear stress together with indoxyl sulfate increased EV release from venous endothelial cells corresponding to higher endothelial EV concentrations in HD patients with arterio-venous fistulas (AVF) and recapitulated EV miRNA changes in vitro observed in CKD in vivo (lower abundance of let-7d-5p and miR-24-3p). Conclusion Endothelial EV release and decreased EV miRNAs associate with advanced CKD while increased endothelial shear stress and indoxyl sulfate form potential causative factors of altered EV phenotype. CKD EVs inhibit angiogenesis in vitro, potentially caused by dysregulated EV miRNAs that could form novel therapeutical targets to reduce CVD burden in CKD.

#2494 Pharmacological enhancement of circadian clock ameliorated fibrosis in adenine induced chronic kidney disease

Abstract Background and Aims Circadian clock regulates diverse biological processes. Recent studies demonstrating the important role of dysregulated molecular clock system in several metabolic diseases as well as in chronic kidney disease (CKD) suggest that modulation of clock machinery might be an attractive therapeutic target. Method We optimized a substance (compound X) that enhances the circadian molecular rhythm and validated its effect in HK-2 cells as well as in mouse model of adenine induced CKD. Results By using transcriptome analysis of whole kidneys, we identified that 19% of CKD associated genes showed rhythmic expressions in circadian manner. Gene enrichment analysis revealed that the rhythmic genes involved in extracellular matrix/collagen/myofibroblasts showed upregulation while those involved in transport/brush border/proximal tubule showed downregulation in CKD. Treatment with compound X in HK-2 cells resulted in marked attenuation of epithelial mesenchymal transition upon transforming growth factor-b (TGF- b) stimulation. Along with these in vitro data, we also observed the preserved kidney function and less interstitial fibrosis in compound X treated CKD mice. The protective effect was associated with suppression of inflammation, cell cycle arrest/cellular senescence as well as restoration of several transporters in proximal tubular cells. Conclusion This is the first study that showed that circadian rhythm might be a novel target for CKD progression by suppressing inflammation/fibrosis and enhancing tubular functions. Unlike other strategies targeting single pathway, targeting circadian rhythm seems to be more promising due to its multi-target effects. A better understanding of the complex role of circadian clock in kidney diseases may provide new perspectives for the prevention and treatment of CKD.

Cuproptosis-related Gene Signatures and Immunological Characterization in Sepsis-associated Acute Lung Injury.

Sepsis is a frequent cause of acute lung injury (ALI), characterized by immune dysregulation and a high mortality rate. The role of cuproptosis, a recently discovered cell death mechanism, in sepsis-associated ALI is still unclear. The study aimed to investigate the regulatory mechanisms and immune characteristics associated with cuproptosis in sepsisassociated ALI, with implications for novel diagnostic and therapeutic approaches. Data from the GEO database was utilized to conduct a comprehensive analysis of the cuproptosis-related genes (CRGs) in sepsis-associated ALI. Gene enrichment analysis, WGCNA, CIBERSORT algorithm, and consensus clustering were employed to investigate the associations between CRGs and immune cells. A predictive model for sepsis-associated ALI was developed based on key CRGs, and its diagnostic accuracy was assessed. Finally, qPCR was employed to validate alterations in the expression of CRGs in the sepsis-associated ALI cellular model. A total of 14 CRGs were identified in sepsis-associated ALI. Strong correlations between the CRGs and immune cells were observed, and two different CRG subtypes were identified. The expression of immune-related factors in both the CRG and gene clusters exhibited similarities, suggesting a connection between the subgroups and immune cells. The prediction model effectively forecasted the incidence of sepsis-associated ALI based on the expression of CRGs. Finally, qPCR analysis confirmed that the expressions of CRGs in the sepsis-associated ALI cell model closely matched those identified through bioinformatic analyses. The study comprehensively evaluated the complex relationship between cuproptosis and sepsis-associated ALI. CRGs were found to be significantly associated with the occurrence, immune characteristics, and biological processes of sepsis-associated ALI. These findings provide valuable new insights into the mechanisms underlying sepsis-associated ALI.

Identification of Autophagy-Related Genes Involved in Intervertebral Disc Degeneration by Microarray Data Analysis

Nucleus pulposus cells survive in a hypoxic, acidic, nutrient-poor, and hypotonic microenvironment. Consequently, they maintain low proliferation and undergo autophagy to protect themselves from cellular stress. Therefore, we aimed to identify autophagy-related biomarkers involved in intervertebral disc degeneration pathogenesis. Autophagy-related differentially expressed genes were derived from the intersection between the public GSE147383 microarray data set to identify differentially expressed genes and online databases to identify autophagy-related genes. Furthermore, we assessed their biological functions with gene annotation and enrichment analysis in the Metscape portal. Then, the STRING database and Cytoscape software allowed inferring a protein-protein interaction (PPI) network and identifying hub genes. In addition, to predict transcription factors that may regulate the hub genes, we used the GeneMANIA website. Finally, the competing endogenous RNA prediction tools and Cytoscape were also used to construct an mRNA-miRNA-lncRNA network. A total of 123 autophagy-related differentially expressed genes were identified, they were mainly involved in phosphoinositide 3-kinase-Akt signaling, autophagy animal, and apoptosis pathways. Nine were identified as hub genes (PTEN, MYC, CTNNB1, JUN, BECN1, ERBB2, FOXO3, ATM, and FN1) and 36 transcription factors were associated with them. Finally, an autophagy-associated competing endogenous RNA network was constructed based on the 9 hub genes. Nine hub genes were identified and a network of competing endogenous RNA associated with autophagy was established. They can be used as autophagy-related biomarkers of intervertebral disc degeneration and for further exploration.

Deciphering the phospho-signature induced by hepatitis B virus in primary human hepatocytes.

Phosphorylation is a major post-translation modification (PTM) of proteins which is finely tuned by the activity of several hundred kinases and phosphatases. It controls most if not all cellular pathways including anti-viral responses. Accordingly, viruses often induce important changes in the phosphorylation of host factors that can either promote or counteract viral replication. Among more than 500 kinases constituting the human kinome only few have been described as important for the hepatitis B virus (HBV) infectious cycle, and most of them intervene during early or late infectious steps by phosphorylating the viral Core (HBc) protein. In addition, little is known on the consequences of HBV infection on the activity of cellular kinases. The objective of this study was to investigate the global impact of HBV infection on the cellular phosphorylation landscape early after infection. For this, primary human hepatocytes (PHHs) were challenged or not with HBV, and a mass spectrometry (MS)-based quantitative phosphoproteomic analysis was conducted 2- and 7-days post-infection. The results indicated that while, as expected, HBV infection only minimally modified the cell proteome, significant changes were observed in the phosphorylation state of several host proteins at both time points. Gene enrichment and ontology analyses of up- and down-phosphorylated proteins revealed common and distinct signatures induced by infection. In particular, HBV infection resulted in up-phosphorylation of proteins involved in DNA damage signaling and repair, RNA metabolism, in particular splicing, and cytoplasmic cell-signaling. Down-phosphorylated proteins were mostly involved in cell signaling and communication. Validation studies carried out on selected up-phosphorylated proteins, revealed that HBV infection induced a DNA damage response characterized by the appearance of 53BP1 foci, the inactivation of which by siRNA increased cccDNA levels. In addition, among up-phosphorylated RNA binding proteins (RBPs), SRRM2, a major scaffold of nuclear speckles behaved as an antiviral factor. In accordance with these findings, kinase prediction analysis indicated that HBV infection upregulates the activity of major kinases involved in DNA repair. These results strongly suggest that HBV infection triggers an intrinsic anti-viral response involving DNA repair factors and RBPs that contribute to reduce HBV replication in cell culture models.

BR regulates wheat root salt tolerance by maintaining ROS homeostasis.

Trace amounts of epibrassinolide (EpiBL) could partially rescue wheat root length inhibition in salt-stressed situation by scavenging ROS, and ectopic expression of TaDWF4 or TaBAK1 enhances root salt tolerance in Arabidopsis by balancing ROS level. Salt stress often leads to ion toxicity and oxidative stress, causing cell structure damage and root development inhibition in plants. While prior research indicated the involvement of exogenous brassinosteroid (BR) in plant responses to salt stress, the precise cytological role and the function of BR in wheat root development under salt stress remain elusive. Our study demonstrates that 100mM NaCl solution inhibits wheat root development, but 5nM EpiBL partially rescues root length inhibition by decreasing H2O2 content, oxygen free radical (OFR) content, along with increasing the peroxidase (POD) and catalase (CAT) activities in salt-stressed roots. The qRT-PCR experiment also shows that expression of the ROS-scavenging genes (GPX2 and CAT2) increased in roots after applying BR, especially during salt stress situation. Transcriptional analysis reveals decreased expression of BR synthesis and root meristem development genes under salt stress in wheat roots. Differential expression gene (DEG) enrichment analysis highlights the significant impact of salt stress on various biological processes, particularly "hydrogen peroxide catabolic process" and "response to oxidative stress". Additionally, the BR biosynthesis pathway is enriched under salt stress conditions. Therefore, we investigated the involvement of wheat BR synthesis gene TaDWF4 and BR signaling gene TaBAK1 in salt stress responses in roots. Our results demonstrate that ectopic expression of TaDWF4 or TaBAK1 enhances salt tolerance in Arabidopsis by balancing ROS (Reactive oxygen species) levels in roots.

Epigenetics in etiopathology of hyperprolactinemia

Aim: Epigenetic alterations have been reported in patients with pituitary tumors and those on antipsychotic drugs, which are also responsible for hyperprolactinemia. This suggests a possible role of epigenetics in the etiopathology of hyperprolactinemia. Methods: The study recruited 83 hyperprolactinemia cases with prolactin > 100 ng/mL and 65 controls. Global DNA methylation status was studied by MethylFlash Methylated DNA Quantification Kit and genome-wide methylation analysis (GWMA) by Infinium Methylation EPIC BeadChip 850K array. Results: Hyperprolactinemia cases showed significant global DNA hypermethylation compared to controls. Around 66.67% of hypomethylated and 12.9% of hypermethylated cases were on antipsychotics. Gene enrichment analysis of 5-cytosine-phosphate-guanine-3 (CpG) site-associated genes demonstrated significantly enriched major histocompatibility complex (MHC)-related protein classes and cellular components. Conclusions: The study suggested the role of epigenetics in the etiopathology of hyperprolactinemia.

Studying Human Pathogenic Cryptococcus Gattii Lineages by Utilizing Simple Sequence Repeats to Create Diagnostic Markers and Analyzing Diversity.

In this study, we compared the occurrence, relative abundance (RA), and density (RD) of simple sequence repeats (SSRs) among the lineages of human pathogenicCryptococcus gattii using an in-silico approach to gain a deeper understanding of the structure and evolution of their genomes. C. gattii isolate MF34 showed the highest RA and RD of SSRs in both the genomic and transcriptomic sequences, followed by isolate WM276. In both the genomic (50%) and transcriptomic (65%) sequences, trinucleotide SSRs were the most common SSR class. A motif conservation study found that the isolates had stronger conservation (56.1%) of motifs, with isolate IND107 having the most (5.7%) unique motifs. We discovered the presence of SSRs in genes that are directly or indirectly associated with disease using gene enrichment analysis. Isolate-specific unique motifs identified in this study could be utilized as molecular probes for isolate identification. To improve genetic resources among C. gattii isolates, 6499 primers were developed. These genomic resources developed in this study could help with diversity analysis and the development of isolate-specific markers.

Human Conjunctival Transcriptome in Acanthamoeba Keratitis: An Exploratory Study.

The purpose of this study was to identify conjunctival transcriptome differences in patients with Acanthamoeba keratitis compared with keratitis with no known associated pathogen. The host conjunctival transcriptome of 9 patients with Acanthamoeba keratitis (AK) is compared with the host conjunctival transcriptome of 13 patients with pathogen-free keratitis. Culture and/or confocal confirmed Acanthamoeba in 8 of 9 participants with AK who underwent metagenomic RNA sequencing as the likely pathogen. Cultures were negative in all 13 cases where metagenomic RNA sequencing did not identify a pathogen. Transcriptome analysis identified 36 genes differently expressed between patients with AK and patients with presumed sterile, or pathogen-free, keratitis. Gene enrichment analysis revealed that some of these genes participate in several biologic pathways important for cellular signaling, ion transport and homeostasis, glucose transport, and mitochondrial metabolism. Notable relatively differentially expressed genes with potential relevance to Acanthamoeba infection included CPS1 , SLC35B4 , STEAP2 , ATP2B2 , NMNAT3 , and AKAP12 . This research suggests that the local transcriptome in Acanthamoeba keratitis may be sufficiently robust to be detected in the conjunctiva and that corneas infected with Acanthamoeba may be distinguished from the inflamed cornea where no pathogen was identified. Given the low sensitivity for corneal cultures, identification of differentially expressed genes may serve as a suggestive transcriptional signature allowing for a complementary diagnostic technique to identify this blinding parasite. Knowledge of differentially expressed genes may also direct investigation of disease pathophysiology and suggest novel pathways for therapeutic targets.

Revealing an adverse outcome pathway network for reproductive toxicity induced by atrazine, via oxidative stress

Adverse outcome pathway AOPs are conceptual frameworks that organize scientific knowledge about how stressors disrupt specific biological targets, pathways. AOP network consists of two or more AOPs that share common key events (KEs), including crucial events like molecular initiating events (MIEs) and adverse outcomes (AOs), which offer the opportunity to link toxicological pathways. Thus, to better understand the sequential series of KEs involved in the AOP 492 (https://aopwiki.org/aops/492), which is triggered by atrazine (ATZ), we first generated a Reproductive Toxicity via Oxidative Stress (RTOS) AOP network from individual AOPs published in the AOP-Wiki database, using this AOP as a seed. The KEs “Increased, Reactive oxygen species” and “Apoptosis” were considered the most common/highly connected KE within this network and an important point of divergence. Furthermore, “Increased, DNA damage and mutations” is a critical KE within the network, as it is highly connected and central, and represents a point of divergence. This suggests that these three KEs have a high predictive value and could, for example, serve as a basis for the development/selection of in vitro assays to assess reproductive toxicity. The in silico analyses revealed that the pivotal target proteins for ATZ-induced infertility via oxidative stress in humans are Tp53, Bcl2, Esr1, and Nos3, which interact indirectly with ATZ via intermediary factors such as Mapk3, Mapk1, and Cyp19a1. Further, the gene enrichment analyses indicate that these entities are involved in several biological processes and pathways directly associated with oxidative stress, DNA damage and apoptosis, further reinforcing the developed network.

Uncovering miRNA-mRNA Regulatory Networks Related to Olaparib Resistance and Resensitization of BRCA2MUT Ovarian Cancer PEO1-OR Cells with the ATR/CHK1 Pathway Inhibitors.

Resistance to olaparib is the major obstacle in targeted therapy for ovarian cancer (OC) with poly(ADP-ribose) polymerase inhibitors (PARPis), prompting studies on novel combination therapies to enhance olaparib efficacy. Despite identifying various mechanisms, understanding how OC cells acquire PARPi resistance remains incomplete. This study investigated microRNA (miRNA) expression in olaparib-sensitive (PEO1, PEO4) and previously established olaparib-resistant OC cell lines (PEO1-OR) using high-throughput RT-qPCR and bioinformatic analyses. The role of miRNAs was explored regarding acquired resistance and resensitization with the ATR/CHK1 pathway inhibitors. Differentially expressed miRNAs were used to construct miRNA-mRNA regulatory networks and perform functional enrichment analyses for target genes with miRNet 2.0. TCGA-OV dataset was analyzed to explore the prognostic value of selected miRNAs and target genes in clinical samples. We identified potential processes associated with olaparib resistance, including cell proliferation, migration, cell cycle, and growth factor signaling. Resensitized PEO1-OR cells were enriched in growth factor signaling via PDGF, EGFR, FGFR1, VEGFR2, and TGFβR, regulation of the cell cycle via the G2/M checkpoint, and caspase-mediated apoptosis. Antibody microarray analysis confirmed dysregulated growth factor expression. The addition of the ATR/CHK1 pathway inhibitors to olaparib downregulated FGF4, FGF6, NT-4, PLGF, and TGFβ1 exclusively in PEO1-OR cells. Survival and differential expression analyses for serous OC patients revealed prognostic miRNAs likely associated with olaparib resistance (miR-99b-5p, miR-424-3p, and miR-505-5p) and resensitization to olaparib (miR-324-5p and miR-424-3p). Essential miRNA-mRNA interactions were reconstructed based on prognostic miRNAs and target genes. In conclusion, our data highlight distinct miRNA profiles in olaparib-sensitive and olaparib-resistant cells, offering molecular insights into overcoming resistance with the ATR/CHK1 inhibitors in OC. Moreover, some miRNAs might serve as potential predictive signature molecules of resistance and therapeutic response.

Genome-wide association analysis explores the genetic loci of amino acid content in duck’s breast muscle

BackgroundAmino acids are the basic components of protein and an important index to evaluate meat quality. With the rapid development of genomics, candidate regions and genes affecting amino acid content in livestock and poultry have been gradually revealed. Hence, genome-wide association study (GWAS) can be used to screen candidate loci associated with amino acid content in duck meat.ResultIn the current study, the content of 16 amino acids was detected in 358 duck breast muscles. The proportion of Glu to the total amino acid content was relatively high, and the proportion was 0.14. However, the proportion of Met content was relatively low, at just 0.03. By comparative analysis, significant differences were found between males and females in 3 amino acids, including Ser, Met, and Phe. In addition, 12 SNPs were significantly correlated with Pro content by GWAS analysis, and these SNPs were annotated by 7 protein-coding genes; 8 significant SNPs were associated with Tyr content, and these SNPs were annotated by 6 protein-coding genes. At the same time, linkage disequilibrium (LD) analysis was performed on these regions with significant signals. The results showed that three SNPs in the 55–56 Mbp region of chromosome 3 were highly correlated with the leader SNP (chr3:55526954) that affected Pro content (r2 > 0.6). Similarly, LD analysis showed that there were three SNPs in the 21.2–21.6 Mbp region of chromosome 13, which were highly correlated with leader SNP (chr13:21421661) (r2 > 0.6). Moreover, Through functional enrichment analysis of all candidate genes. The results of GO enrichment analysis showed that several significant GO items were associated with amino acid transport function, including amino acid transmembrane transport and glutamine transport. The results further indicate that these candidate genes are closely associated with amino acid transport. Among them, key candidate genes include SLC38A1. For KEGG enrichment analysis, CACNA2D3 and CACNA1D genes were covered by significant pathways.ConclusionIn this study, GWAS analysis found a total of 28 significant SNPs affecting amino acid content. Through gene annotation, a total of 20 candidate genes were screened. In addition, Through LD analysis and enrichment analysis, we considered that SERAC1, CACNA2D3 and SLC38A1 genes are important candidate genes affecting amino acid content in duck breast muscle.

MRNA and miRNA expression profiles reveal the potential roles of RLRs signaling pathway and mitophagy in duck hepatitis A virus type 1 infection

Duck hepatitis A virus 1 (DHAV-1) is the primary cause of duck viral hepatitis, leading to sudden mortality in ducklings and significant economic losses in the duck industry. However, little is known about how DHAV-1 affects duckling liver at the molecular level. We conducted an analysis comparing the expression patterns of mRNAs and miRNAs in DHAV-1-infected duckling livers to understand the underlying mechanisms and dynamic changes. We identified 6,818 differentially expressed mRNAs (DEGs) and 144 differentially expressed microRNAs (DEMs) during DHAV-1 infection. Functional enrichment analysis of DEGs and miRNA target genes using gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) revealed their potential involvement in innate antiviral immunity, mitophagy, and pyroptosis. We constructed coexpression networks of mRNA-miRNA interactions and confirmed key DEMs (novel-mir333, novel-mir288, novel-mir197, and novel-mir71) using RT-qPCR. Further investigation demonstrated that DHAV-1 activates the RLRs signaling pathway, disrupts mitophagy, and induces pyroptosis. In conclusion, DHAV-1-induced antiviral immunity is closely linked to mitophagy, suggesting it could be a promising therapeutic target.

Combating anoikis resistance: bioactive compounds transforming prostate cancer therapy.

The study aims to discuss the challenges associated with treating prostate cancer (PCa), which is known for its complexity and drug resistance. It attempts to find differentially expressed genes (DEGs), such as those linked to anoikis resistance and circulating tumor cells, in PCa samples. This study involves analyzing the functional roles of these DEGs using gene enrichment analysis, and then screening of 102 bioactive compounds to identify a combination that can control the expression of the identified DEGs. In this study, 53 DEGs were identified from PCa samples including anoikis-resistant PCa cells and circulating tumor cells in PCa. Gene enrichment analysis with regards to functional enrichment of DEGs was performed. An inclusive screening process was carried out among 102 bioactive compounds to identify a combination capable of affecting and regulating the expression of selected DEGs. Eventually, gastrodin, nitidine chloride, chenodeoxycholic acid, and bilobalide were selected, as their combination demonstrated ability to modulate expression of 50 out of the 53 genes targeted. The subsequent analysis focused on investigating the biological pathways and processes influenced by this combination. The findings revealed a multifaceted and multidimensional approach to tumor regression. The combination of bioactive compounds exhibited effects on various genes including those related to production of inflammatory cytokines, cell proliferation, autophagy, apoptosis, angiogenesis, and metastasis. The current study has made a valuable contribution to the development of a combination of bioactive natural compounds that can significantly impede the development of treatment resistance in prostate tumor while countering the tumors' evasion of the immune system. The implications of this study are highly significant as it suggests the creation of an enhanced immunotherapeutic, natural therapeutic concoction with combinatorial potential.

Identification of microRNAs in black tiger shrimp (Penaeus monodon) under acute low-salinity stress

Salinity is a common abiotic stress in the culture of penaeid shrimp. Through post-transcriptional regulation of gene transcripts, microRNAs (miRNAs) play an important role in the adaptation to a stressful environment. However, the involvement of miRNAs in the salinity stress response of shrimp remains unclear. In the present study, the sequence and expression profile of miRNAs in the hepatopancreas of low-salinity-treated Penaeus monodon were obtained by the high-throughput sequencing technique. A total of 679 miRNAs were identified, including 167 miRNAs that were significantly differentially expressed after low-salinity exposure (p < 0.05). Remarkably, most of these miRNAs were downregulated, suggesting that a series of genes were activated to participate in stress response. In addition, 43 miRNAs differentially expressed at all treatment were selected as putative key modulators. Enrichment analysis of genes targeted by these miRNAs indicated that a network that consists of the nervous system, the immune system, and the endocrine system played a crucial role in maintaining the homeostasis of P. monodon under low-salinity stress. These findings may help contribute to a better understanding of the mechanism that regulates salinity tolerance in shrimp and provide valuable genetic information for subsequent studies.

Identification of circRNA–miRNA–mRNA regulatory network associated to the autism spectrum disorder in children through integrated bioinformatics analysis

BackgroundAutism spectrum disorder (ASD) is a complex neurological disability with multifactorial etiology. ASD is described by behavior, speech, language, and communication defects. CircRNA is a type of ceRNA that plays an important role in modulating microRNAs (miRNA) in several disorders. However, the potential role of the circRNA/miRNA/mRNA regulatory network in the pathogenesis of ASD is not fully understood. Therefore, this study aimed to create a circRNA/miRNA/mRNA network associated with ASD to cast light on the pathogenesis of ASD.MethodsCircRNA expression profile data were recruited from Gene Expression Omnibus datasets, and the differentially expressed circRNAs (DEcircRNAs) were identified. Then, miRNAs modulated by these circRNAs were predicted and overlapped with differentially expressed miRNAs. Next, the potentially involved genes were identified by overlapping predicted targets, and differentially expressed genes. The enrichment analysis was performed, and a PPI network was projected. Subsequently, ten key genes were selected from the network. Furthermore, a circRNA/miRNA/mRNA regulatory network was constructed, and probable molecules and drugs with potential anti-ASD effects were predicted.Results11 DEcircRNAs and 8 miRNAs regulated by 4 circRNAs were identified as being significantly involved. Subsequently, gene enrichment analysis of 71 overlapped mRNA regulated by these miRNAs showed that they are mostly associated with hippocampal synaptogenesis, neurogenesis, and axon guidance. Additionally, two high-score compounds, GSK3β inhibitor (SB216763) and dexamethasone, and three drugs (haloperidol, nystatin, paroxetine) were confirmed as potential therapeutic options for ASD.ConclusionThe results of this study may help gain deeper insight into the pathogenesis of the circRNA/miRNA/mRNA regulatory network in ASD, providing potential therapeutic management options.

Differential expression analysis identifies a prognostically significant extracellular matrix–enriched gene signature in hyaluronan-positive clear cell renal cell carcinoma

Hyaluronan (HA) accumulation in clear cell renal cell carcinoma (ccRCC) is associated with poor prognosis; however, its biology and role in tumorigenesis are unknown. RNA sequencing of 48 HA-positive and 48 HA-negative formalin-fixed paraffin-embedded (FFPE) samples was performed to identify differentially expressed genes (DEG). The DEGs were subjected to pathway and gene enrichment analyses. The Cancer Genome Atlas Kidney Renal Clear Cell Carcinoma (TCGA-KIRC) data and DEGs were used for the cluster analysis. In total, 129 DEGs were identified. HA-positive tumors exhibited enhanced expression of genes related to extracellular matrix (ECM) organization and ECM receptor interaction pathways. Gene set enrichment analysis showed that epithelial–mesenchymal transition-associated genes were highly enriched in the HA-positive phenotype. A protein–protein interaction network was constructed, and 17 hub genes were discovered. Heatmap analysis of TCGA-KIRC data identified two prognostic clusters corresponding to HA-positive and HA-negative phenotypes. These clusters were used to verify the expression levels and conduct survival analysis of the hub genes, 11 of which were linked to poor prognosis. These findings enhance our understanding of hyaluronan in ccRCC.

RRP8, associated with immune infiltration, is a prospective therapeutic target in hepatocellular carcinoma

BackgroundRibosomal RNA Processing 8 (RRP8) is a nucleolar Rossman fold-like methyltransferase that exhibits increased expression in many malignant tumours. However, the role of RRP8 in hepatocellular carcinoma (HCC) is still uncertain. We explored the relationships between RRP8 and prognosis and immune infiltration, as well as the putative pathological function and mechanism of RRP8 in HCC.MethodsAnalysis of RRP8 expression across cancers was performed by using multiple databases. Associations between RRP8 expression and clinicopathological factors were further examined. Gene enrichment analysis was used to identify various putative biological activities and regulatory networks of RRP8 in HCC. The relationship between RRP8 expression and immune infiltration was confirmed by single-sample gene set enrichment analysis (ssGSEA). Univariate and multivariate Cox regression analyses were conducted to assess the impact of clinical variables on patient outcomes. Furthermore, a nomogram was constructed to estimate survival probability based on multivariate Cox regression analysis. Functional validation of RRP8 in HCC was performed with two different systems: doxycycline-inducible shRNA knockdown and CRISPR-Cas9 knockout.ResultsRRP8 was markedly overexpressed in HCC clinical specimens compared to adjacent normal tissues. Further analysis demonstrated that RRP8 was directly connected to multiple clinical characteristics and strongly associated with various immune markers in HCC. Moreover, elevated RRP8 expression indicated an unfavourable prognosis. Our functional studies revealed that both knockdown and knockout of RRP8 dramatically attenuated liver cancer cells to proliferate and migrate. Knockout of RRP8 decreased the phosphorylation of MEK1/2 and β-catenin-(Y654) signalling pathway components; downregulated downstream signalling effectors, including Cyclin D1 and N-cadherin; and upregulated E-cadherin.ConclusionsRRP8 is strongly implicated in immune infiltration and could be a potential therapeutic target in HCC.

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

BackgroundIn 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.ResultsA 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.ConclusionThe 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.

SLAMF8 as a potential biomarker for rheumatoid arthritis identified by comparing peripheral blood mononuclear cells, fibroblast‐like synoviocytes, and synovial tissue using bioinformatics analysis

AbstractBackgroundRheumatoid arthritis (RA) is a chronic autoimmune inflammatory disease. Its pathogenesis is not fully understood, and early diagnosis is challenging owing to the lack of effective biomarkers. This study aimed to analyze different samples to identify potential biomarkers and therapeutic targets.MethodsMicroarray datasets of RA, osteoarthritis (OA), and healthy control (HC) were downloaded from the Gene Expression Omnibus database. R software was used to identify differentially expressed genes (DEGs), which were visualized using volcano and heat maps. Venn diagrams, principal component analysis, gene set enrichment analysis, gene ontology, and Kyoto Encyclopedia of genes and genomes were used to analyze the data. A protein–protein interaction network was constructed, and synovial tissues from patients with RA and OA were collected for verification using the collagen‐induced arthritis mouse model.ResultsMore DEGs were found in synovial tissues than in peripheral blood mononuclear cells or fibroblast‐like synoviocytes. Principal component analysis revealed significant differences between the RA and OA samples, highlighting the unique advantages of synovial tissue. Enrichment analysis revealed that metabolic and cytokine signaling pathways play crucial roles in the development of RA. Further analysis of the four synovial datasets identified 54 DEGs, of which signaling lymphocytic activation molecule family (SLAMF) 8 was identified as the key molecule. SLAMF8 levels were increased in the synovial tissue of patients with RA compared to those of patients with OA (0.38 ± 0.19 vs. 12.40 ± 1.66), and SLAMF8 levels were similarly elevated in collagen‐induced arthritis model mice compared with those in the healthy mice (1.13 ± 0.47 vs. 9.05 ± 2.52).ConclusionsThis study established the unique advantages of synovial tissue for RA research and identified metabolic and cytokine signaling pathways as important for RA development. Thus, SLAMF8 may be a potential therapeutic target for RA.