Abstract PICK1 plays a crucial role in mammalian spermatogenesis. Here, we integrated single-molecule long-read and short-read sequencing to comprehensively examine PICK1 expression patterns in adult Baoshan pig (BS) testes. We identified the most important transcript ENSSSCT00000000120 of PICK1, obtaining its full-length coding sequence (CDS) spanning 1254 bp. Gene structure analysis located PICK1 on pig chromosome 5 with 14 exons. Protein structure analysis reflected that PICK1 consisted of 417 amino acids containing two conserved domains, PDZ and BAR_PICK1. Phylogenetic analysis underscored the evolutionary conservation and homology of PICK1 across different mammalian species. Evaluation of protein interaction network, KEGG, and GO pathways implied that interacted with 50 proteins, predominantly involved in glutamatergic synapses, amphetamine addiction, neuroactive ligand-receptor interactions, dopaminergic synapses, and synaptic vesicle recycling, and PICK1 exhibited significant correlation with DLG4 and TBC1D20. Functional annotation identified that PICK1 was involved in 9 GOs, including seven cellular components and two molecular functions. ceRNA network analysis suggested BS PICK1 was regulated by seven miRNA targets. Moreover, qPCR expression analysis across 15 tissues highlighted that PICK1 was highly expressed in the bulbourethral gland and testis. Subcellular localization analysis in ST (Swine Tesits) cells demonstrated that PICK1 significantly localized within the cytoplasm. Overall, our findings shed new light on PICK1’s role in BS reproduction, providing a foundation for further functional studies of PICK1. PICK spermatogenesis Here singlemolecule single molecule longread long read shortread short (BS testes ENSSSCT ENSSSCT0000000012 fulllength full length CDS (CDS 125 bp 1 exons 41 domains BARPICK1 BARPICK BAR BAR_PICK1 species KEGG proteins synapses addiction ligandreceptor ligand receptor interactions recycling DLG TBC1D20 TBCD TBC D GOs functions targets Moreover testis Swine Tesits cytoplasm Overall PICK1s PICKs s reproduction ENSSSCT000000001 12 4 BAR_PICK TBC1D2 ENSSSCT00000000 TBC1D ENSSSCT0000000 ENSSSCT000000 ENSSSCT00000 ENSSSCT0000 ENSSSCT000 ENSSSCT00 ENSSSCT0

Objective: The present study combined transcriptomic data and computational techniques based on gene expression signatures to identify new bioactive compounds or Food and Drug Administration-approved drugs for the treatment of bipolar disorder (BD). Methods: Five transcriptomic datasets containing 165 blood samples from individuals with BD were selected from the Gene Expression Omnibus (GEO). The number of participants varied from six to 60, with a mean age between 35 and 48 years and a gender difference between them. Most of these patients were receiving pharmacological treatment. Master regulator analysis (MRA) and gene set enrichment analysis (GSEA) were performed to identify genes that were significantly different between patients with BD and healthy controls and their associations with mood states in patients with BD. In addition, molecules that could reverse the transcriptomic profiles of BD-altered regulons were identified from the Library of Network-Based Cellular Signatures Consortium (LINCS) and the Broad Institute Connectivity Map Drug Repurposing Database (cMap) databases. Results: MRA identified 59 candidate master regulators (MRs) that modulate regulatory units enriched with BD-altered genes. In contrast, GSEA identified 134 enriched genes and 982 regulons whose activation state was determined. Both analyses revealed genes exclusively associated with mania, depression, or euthymia, and some genes were shared among these three mood states. We identified bioactive compounds and licensed drug candidates, including antihypertensives and antineoplastic agents, as promising candidates for the treatment of BD. However, experimental validation is essential to confirm these findings in further studies. Conclusion: Although our data are still preliminary, they provide some insights into the biological patterns of different mood states in patients with BD and their potential therapeutic targets. The strategy of transcriptomics plus bioinformatics offers a way to advance drug discovery and personalized medicine by using gene expression information. Objective Administrationapproved Administration approved . (BD) Methods 16 GEO. GEO (GEO) 60 3 4 them (MRA (GSEA addition BDaltered altered NetworkBased Network Based LINCS (LINCS cMap (cMap databases Results 5 MRs (MRs contrast 13 98 determined mania depression euthymia agents However studies Conclusion preliminary targets information (BD 1 (GEO 6 9

ABSTRACT Objective: To identify DNA methylation and gene expression profiles involved in obesity by implementing an integrated bioinformatics approach. Materials and methods: Gene expression (GSE94752, GSE55200, and GSE48964) and DNA methylation (GSE67024 and GSE111632) datasets were obtained from the GEO database. Differentially expressed genes (DEGs) and differentially methylated genes (DMGs) in subcutaneous adipose tissue of patients with obesity were identified using GEO2R. Methylation-regulated DEGs (MeDEGs) were identified by overlapping DEGs and DMGs. The protein–protein interaction (PPI) network was constructed with the STRING database and analyzed using Cytoscape. Functional modules and hub-bottleneck genes were identified by using MCODE and CytoHubba plugins. Functional enrichment analyses were performed based on Gene Ontology terms and KEGG pathways. To prioritize and identify candidate genes for obesity, MeDEGs were compared with obesity-related genes available at the DisGeNET database. Results: A total of 54 MeDEGs were identified after overlapping the lists of significant 274 DEGs and 11,556 DMGs. Of these, 25 were hypermethylated-low expression genes and 29 were hypomethylated-high expression genes. The PPI network showed three hub-bottleneck genes (PTGS2, TNFAIP3, and FBXL20) and one functional module. The 54 MeDEGs were mainly involved in the regulation of fibroblast growth factor production, the molecular function of arachidonic acid, and ubiquitin-protein transferase activity. Data collected from DisGeNET showed that 11 of the 54 MeDEGs were involved in obesity. Conclusion: This study identifies new MeDEGs involved in obesity and assessed their related pathways and functions. These results data may provide a deeper understanding of methylation-mediated regulatory mechanisms of obesity. Objective approach methods GSE94752, GSE94752 GSE (GSE94752 GSE55200 GSE48964 GSE67024 (GSE6702 GSE111632 (DEGs DMGs (DMGs GEO2R GEOR R Methylationregulated Methylation regulated (MeDEGs proteinprotein protein (PPI Cytoscape hubbottleneck hub bottleneck plugins obesityrelated Results 5 27 11556 556 11,55 these 2 hypermethylatedlow hypermethylated low hypomethylatedhigh hypomethylated high PTGS2, PTGS2 PTGS (PTGS2 TNFAIP3 TNFAIP FBXL20 FBXL module production acid ubiquitinprotein ubiquitin activity 1 Conclusion functions methylationmediated mediated GSE9475 (GSE9475 GSE5520 GSE4896 GSE6702 (GSE670 GSE11163 1155 55 11,5 (PTGS FBXL2 GSE947 (GSE947 GSE552 GSE489 GSE670 (GSE67 GSE1116 115 11, GSE94 (GSE94 GSE55 GSE48 GSE67 (GSE6 GSE111 GSE9 (GSE9 GSE5 GSE4 GSE6 (GSE GSE11 GSE1

Abstract This study aimed to explore the biomarkers associated with osteoarthritis (OA). Gene expression profile GSE16464 was downloaded from the Gene Expression Omnibus database. The differentially expressed genes (DEGs) in chondrocytes between OA patients and normal donors were analyzed, which were then subjected to Gene Ontology (GO) and pathway enrichment analyses, followed by microRNA (miRNA) and transcription factor (TF) prediction, and regulatory network construction. Finally, the copy number variants (CNVs) of target genes were searched. Total 79 up- and 147 down-regulated DEGs were identified. Nine GO terms were obtained and the biological process (BP) term related to neuron differentiation was enriched by 13 DEGs, such as GNAO1, POU3F4 andRPS27A. Pathway of regulation of actin cytoskeleton was enriched by six DEGs such as FGF18. Six miRNAs such as miRNA-524 and seven TFs, such as FOXO4 were detected. In the regulatory network, GNAO1, POU3F4 and RPS27A were key target genes and their CNVs were identified. Pathway of regulation of actin cytoskeleton and BP related to neuron differentiation may play important roles in OA progression. DEGs of FGF18, GNAO1 and POU3F4 as well as their regulatory factors such as FOXO4 and miRNA-524 may be potential biomarkers associated with OA.

Abstract Objective: Long Non-Coding RNAs (LncRNAs) act as an indispensable role in cancer development. The study aimed to investigate the role and mechanism of lncRNA Small Nucleolar RNA Host Gene 1 (SNHG1) in Bladder Cancer (BC) progression. Method: The expression, prognostic value, diagnostic value, and correlation of SNHG1, Enhancer of Zeste 2 polycomb repressive complex 2 subunit (EZH2), and Kruppel Like Factor 2 (KLF2) were analyzed through bioinformatics analysis. The expression was also validated in BC tissues and cell lines. Besides, their regulation and binding were tested via qPCR, Western blot, Dual-Luciferase Reporter Assay (DLRA), Argonaute RISC catalytic component 2-RNA Immunoprecipitation (AGO2-RIP), and Chromatin Immunoprecipitation (ChIP). A xenograft model in nude mice was also established. Results: SNHG1 was significantly overexpressed in BC tissues and cells. Importantly, SNHG1 was associated with poor survival, and ROC curves revealed high diagnostic values. Moreover, by CCK8, wound healing, transwell, and Western blot analysis, SNHG1 knockdown significantly inhibited the proliferation, migration, invasion, and epithelial-mesenchymal transition of BC cells. Additionally, in vivo experiments showed that silencing SNHG1 hindered tumorigenesis and tumor growth. Regarding mechanism, the results of AGO2-RIP, ChIP or DLRA showed that SNHG1 played different roles at diverse subcellular sites. In the cytoplasm, SNHG1 acted as a competing endogenous RNA for miR-137-3p to promote EZH2 expression. In the nucleus, SNHG1 could interact with EZH2 to inhibit KLF2 transcription. Conclusion: Our study elucidated that SNHG1 formed a regulatory network and played an oncogenic role in BC, which provided a novel therapeutic target for BC treatment.

BACKGROUND Healthcare-associated infections due to multidrug-resistant (MDR) bacteria such as Pseudomonas aeruginosa are significant public health issues worldwide. A system biology approach can help understand bacterial behaviour and provide novel ways to identify potential therapeutic targets and develop new drugs. Gene regulatory networks (GRN) are examples of in silico representation of interaction between regulatory genes and their targets. OBJECTIVES In this work, we update the MDR P. aeruginosa CCBH4851 GRN reconstruction and analyse and discuss its structural properties. METHODS We based this study on the gene orthology inference methodology using the reciprocal best hit method. The P. aeruginosa CCBH4851 genome and GRN, published in 2019, and the P. aeruginosa PAO1 GRN, published in 2020, were used for this update reconstruction process. FINDINGS Our result is a GRN with a greater number of regulatory genes, target genes, and interactions compared to the previous networks, and its structural properties are consistent with the complexity of biological networks and the biological features of P. aeruginosa. MAIN CONCLUSIONS Here, we present the largest and most complete version of P. aeruginosa GRN published to this date, to the best of our knowledge.

Abstract Background: Cytokinin signal transduction is mediated by a two-component system (TCS). Two-component systems are utilized in plant responses to hormones as well as to biotic and abiotic environmental stimuli. In plants, response regulatory genes (RRs) are one of the main members of the two-component system (TCS). Method: From the aspects of gene structure, evolution mode, expression type, regulatory network and gene function, the evolution process and role of RR genes in the evolution of the cotton genome were analyzed. Result: A total of 284 RR genes in four cotton species were identified. Including 1049 orthologous/paralogous gene pairs were identified, most of which were whole genome duplication (WGD). The RR genes promoter elements contain phytohormone responses and abiotic or biotic stress-related cis-elements. Expression analysis showed that RR genes family may be negatively regulate and involved in salt stress and drought stress in plants. Protein regulatory network analysis showed that RR family proteins are involved in regulating the DNA-binding transcription factor activity (COG5641) pathway and HP kinase pathways. VIGS analysis showed that the GhRR7 gene may be in the same regulatory pathway as GhAHP5 and GhPHYB, ultimately negatively regulating cotton drought stress by regulating POD, SOD, CAT, H2O2 and other reactive oxygen removal systems. Conclusion: This study is the first to gain insight into RR gene members in cotton. Our research lays the foundation for discovering the genes related to drought and salt tolerance and creating new cotton germplasm materials for drought and salt tolerance.

Resumo CKB3 é uma subunidade reguladora (beta) de CK2. Neste estudo, o mutante homozigoto ckb3 de Arabidopsis thaliana foi estudado para entender o papel da CKB3 na sinalização de ácido abscísico (ABA). Os resultados apresentados: CKB3 foi expresso em todos os órgãos e a maior expressão nas sementes, seguida pela raiz. Durante a germinação das sementes e o crescimento radicular, o mutante ckb3 mostrou sensibilidade reduzida ao ABA. O mutante ckb3 teve mais abertura estomática e aumento do acúmulo de prolina e perda de água nas folhas. Os níveis de expressão do número de genes na rede reguladora da ABA haviam mudado. Este estudo demonstra que CKB3 é um gene relacionado à sinalização ABA e pode desempenhar um papel positivo na sinalização ABA.

Abstract CKB3 is a regulatory (beta) subunit of CK2. In this study Arabidopsis thaliana homozygous T-DNA mutant ckb3 was studied to understand the role of CKB3 in abscisic acid (ABA) signaling. The results shown: CKB3 was expressed in all organs and the highest expression in the seeds, followed by the root. During seed germination and root growth the ckb3 mutant showed reduced sensitivity to ABA. The ckb3 mutant had more stomatal opening and increased proline accumulation and leaf water loss. The expression levels of number of genes in the ABA regulatory network had changed. This study demonstrates that CKB3 is an ABA signaling-related gene and may play a positive role in ABA signaling.

Abstract Breast cancer (BC) is a heterogeneous disease, and it is the leading cause of death among women. NORAD and HCG11 are highly similar lncRNAs that present binding sites for PUMILIO proteins. PUMILIO acts on hundreds of mRNA targets, contributing to the modulation of gene expression. We analyzed the expression levels of NORAD and HCG11 in the BC subtypes luminal A (LA) and basal-like (BL), and the regulatory networks associated with these lncRNAs. In the analysis of TCGA cohort (n=329) and Brazilian BC samples (n=44), NORAD was up-regulated in LA while HCG11 was up-regulated in BL subtype. An increased expression of NORAD is associated with reduced disease-free survival in basal-like patients (p = 0.002), which suggests that its prognostic value could be different in specific subtypes. The biological pathways observed for the HCG11 network are linked to the epithelial-to-mesenchymal transition; while NORAD associated pathways appear to be related to luminal epithelial cell transformation. NORAD and HCG11 regulons respectively present 36% and 21.5% of PUMILIO targets, which suggests that these lncRNAs act as a decoy for PUMILIO. These lncRNAs seem to work as players in the differentiation process that drives breast cells to acquire distinct phenotypes related to a specific BC subtype.

Abstract The agricultural yield of sugarcane (Saccharum spp.) is influenced by various abiotic stresses, including aluminum toxicity (Al3+). MicroRNAs (miRNAs) play a role in plant tolerance to such stresses by modulating the expression of several important target genes involved in plant growth. This study investigated the possible tolerance mechanisms of two sugarcane genotypes (CTC-2 and RB855453) under Al3+ stress through miRNA expression profiles and in silico analysis of target genes. The expression data obtained using RT-qPCR and co-expression network analysis identified two possible regulatory mechanisms in the tolerant genotype (CTC-2) under Al3+ stress. miR395 was involved in Al3+ detoxification, whereas miR160, miR6225-5p, and miR167 participated in the process of lateral root formation, conferring tolerance to the genotype. These findings might be useful for biotechnological strategies that aim for miRNA silencing or gene overexpression and provide subsidies for future genetic improvement programs aimed at the development of abiotic stress-tolerant sugarcane genotypes.

Abstract This study aimed to identify potential therapeutic targets in osteosarcoma (OS) through the network analysis of competing endogenous RNAs (ceRNAs). The differentially expressed miRNAs (DEMIs) and mRNAs (DEMs) were identified between OS cell lines and human mesenchymal stem cells (hMSCs) from the data deposited under GSE70415 using limma package. Functional analysis of DEMs was performed using DAVID and clusterProfiler to identify significantly enriched Gene Ontology biological processes and KEGG pathways, respectively. The DEMI-DEM interaction network was constructed using Cytoscape. LncRNA–miRNA interactions were predicted using starBase database. The ceRNA regulatory network was constructed by integrating mRNAs, miRNAs, and lncRNAs, and functional enrichment analysis was performed for the genes involved. The analysis revealed a total of 326 DEMs and 54 DEMIs between OS cells and hMSCs. We identified several novel therapeutic targets involved in the progression and metastasis of OS, such as CBX7, RAD9A, SNHG7 and miR-34a-5p. The miRNA, miR-543 (target gene: CBX7) was found to be associated with the pathway Mucin type O-glycan biosynthesis. Using the ceRNA network, we established the following regulatory interactions: NEAT1/miR-543/CBX7, SNHG7/miR-34a-5p/RAD9A, and XIST/miR-34a-5p/RAD9A. CBX7, RAD9A, lncRNA SNHG7, miR-543, and miR-34a-5p may be explored as novel therapeutic targets for treatment of OS.

Abstract Trichoderma reesei is the main filamentous fungus used in industry to produce cellulases. Here we investigated the role of CRZ1 and Ca2+signaling in the fungus T. reesei QM6a concerning holocellulases production. For this, we first searched for potential CRZ1 binding sites in promoter regions of key genes coding holocellulases, as well as transcriptional regulators and sugar and calcium transporters. Using a nearly constructed T. reeseiAcrz1 strain, we demonstrated that most of the genes expected to be regulated by CRZ1 were affected in the mutant strain induced with sugarcane bagasse (SCB) and cellulose. In particular, our data demonstrate that Ca2+ acts synergistically with CRZ1 to modulate gene expression, but also exerts CRZ1-independent regulatory role in gene expression in T. reesei, highlighting the role of the major regulator Ca2+ on the signaling for holocellulases transcriptional control in the most part of cellulases genes here investigated. This work presents new evidence on the regulatory role of CRZ1 and Ca2+ sensing in the regulation of cellulolytic enzymes in T. reesei, evidencing significant and previously unknown function of this Ca2+sensing system in the control key transcriptional regulators (XYR1 and CRE1) and on the expression of genes related to sugar and Ca2+ transport.

Abstract Anti-androgen therapies, including orchiectomy, are effective at promoting prostate cancer remission, but are followed by progression to the more aggressive castration-resistant prostate cancer (CRPC). Castration promotes gland and tumor shrinkage. However, prostate adaptation to androgen deprivation involves striking parallel events, all requiring changes in gene expression. We hypothesized that transcription factors (TF) and other transcription-related genes are needed to orchestrate those changes. In this work, downstream analysis using bioinformatic tools and published microarray data allowed us to identify sixty transcriptional regulators (including 10 TF) and to integrate their function in physiologically relevant networks. Functional associations revealed a connection between Arnt, Bhlhe41 and Dbp circadian rhythm genes with the Ar circuitry and a small gene network centered in Pex14, which might indicate a previously unanticipated metabolic shift. We have also identified human homologs and mapped the corresponding genes to human chromosome regions commonly affected in prostate cancer, with particular attention to the PTEN/HHEX/MXI1 cluster at 10q23-25 (frequently deleted in PCa) and to MAPK1 at 22q11.21 (delete in intermediate risk but not in high risk PCa). Twenty genes were found mutated or with copy number alterations in at least five percent of three cancer cohorts and six of them (PHOX2A, NFYC, EST2, EIF2S1, SSRP1 and PARP1) associated with impacted patient survival. These changes are specific to the adaptation to the hypoandrogen environment and seem important for the progression to CRPC when mutated.

Abstract Background: Hematoporphyrin derivative (HPD) has a sensibilization effect in lung adenocarcinoma. This study was conducted to identify the target genes of HPD in lung adenocarcinoma. Methods: RNA sequencing was performed using the lung adenocarcinoma cell line A549 after no treatment or treatment with X-ray or X-ray + HPD. The differentially expressed genes (DEGs) were screened using Mfuzz package by noise-robust soft clustering analysis. Enrichment analysis was carried out using “BioCloud” online tool. Protein–protein interaction (PPI) network and module analyses were performed using Cytoscape software. Using WebGestalt tool and integrated transcription factor platform (ITFP), microRNA target and transcription factor (TF) target pairs were separately predicted. An integrated regulatory network was visualized with Cytoscape software. Results: A total of 815 DEGs in the gene set G1 (continuously dysregulated genes along with changes in processing conditions [untreated—treated with X-ray—X-ray + treated with HPD]) and 464 DEGs in the gene set G2 (significantly dysregulated between X-ray + HPD-treated group and untreated/X-ray-treated group) were screened. The significant module identified from the PPI network for gene set G1 showed that ribosomal protein L3 (RPL3) gene could interact with heat shock protein 90 kDa alpha, class A member 1 (HSP90AA1). TFs AAA domain containing 2 (ATAD2) and protein inhibitor of activated STAT 1 (PIAS1) were separately predicted for the genes in gene set G1 and G2, respectively. In the integrated network for gene set G2, ubiquitin-specific peptidase 25 (USP25) was targeted by miR-200b, miR-200c, and miR-429. Conclusion: RPL3, HSP90AA1, ATAD2, and PIAS1 as well as USP25, which is targeted by miR-200b, miR-200c, and miR-429, may be the potential targets of HPD in lung adenocarcinoma.