The discovery of upstream regulatory genes of a gene of interest still remains challenging. Here we applied a scalable computational method to unbiasedly predict candidate regulatory genes of critical transcription factors by searching the whole genome. We illustrated our approach with a case study on the master regulator FOXP3 of human primary regulatory T cells (Tregs). While target genes of FOXP3 have been identified, its upstream regulatory machinery still remains elusive. Our methodology selected five top-ranked candidates that were tested via proof-of-concept experiments. Following knockdown, three out of five candidates showed significant effects on the mRNA expression of FOXP3 across multiple donors. This provides insights into the regulatory mechanisms modulating FOXP3 transcriptional expression in Tregs. Overall, at the genome level this represents a high level of accuracy in predicting upstream regulatory genes of key genes of interest.

UNASSIGNED: Lymphangiogenesis (LYM) has an important role in tumor progression and is strongly associated with tumor metastasis. However, the clinical application of LYM has not progressed as expected. The potential value of LYM needs to be further developed in lung adenocarcinoma (LUAD) patients.
UNASSIGNED: The Sequencing data and clinical characteristics of LUAD patients were downloaded from The Cancer Genome Atlas and GEO databases. Multiple machine learning algorithms were used to screen feature genes and develop the LYM index. Immune cell infiltration, immune checkpoint expression, Tumor Immune Dysfunction and Exclusion (TIDE) algorithm and drug sensitivity analysis were used to explore the correlation of LYM index with immune profile and anti-tumor therapy.
UNASSIGNED: We screened four lymphangiogenic feature genes (PECAM1, TIMP1, CXCL5 and PDGFB) to construct LYM index based on multiple machine learning algorithms. We divided LUAD patients into the high LYM index group and the low LYM index group based on the median LYM index. LYM index is a risk factor for the prognosis of LUAD patients. In addition, there was a significant difference in immune profile between high LYM index and low LYM index groups. LUAD patients in the low LYM index group seemed to benefit more from immunotherapy based on the results of TIDE algorithm.
UNASSIGNED: Overall, we confirmed that the LYM index is a prognostic risk factor and a valuable predictor of immunotherapy response in LUAD patients, which provides new evidence for the potential application of LYM.

Histone deacetylases (HDACs) contribute significantly to the initiation, progression, and prognosis of colorectal adenocarcinoma (COAD). Additionally, HDACs regulate the tumor microenvironment, immune escape, and tumor stem cells, and are closely linked to COAD prognosis. We developed a prognostic model for COAD that incorporates HDACs to evaluate their specific roles. The COAD dataset containing clinical and mutation data was collected using the TCGA and GEO databases to obtain genes associated with HDAC. LASSO analysis and univariate and multivariate Cox regression analysis were used to determine the presence of prognostic genes. Multivariate Cox analysis was also used to determine risk scores for HDAC-related features. Furthermore, genomic alterations, immune infiltration, and drug response were compared between high- and low-risk groups. Cellular experiments validated the potential regulatory role of BRD3 on COAD proliferation, migration, and apoptosis. The median risk scores, calculated based on the characteristics, demonstrated a more significant prognostic improvement in patients in the low-risk group. Furthermore, HDAC-related features were identified as important independent prognostic factors for patients with COAD. Additionally, genomic mutation status, immune infiltration, and function, as well as response to immunotherapy and chemotherapy, were found to be associated with risk scores. Subgroup analyses indicate that anti-PD-1 therapy may be beneficial for patients in the low-risk group. Additionally, a decrease in risk score was associated with a decrease in immune infiltration. Finally, HCT116 and HT29 cells exhibited inhibition of BRD3 gene proliferation and migration, as well as promotion of apoptosis. In patients with COAD, HDAC-related characteristics may be useful in predicting survival and selecting treatment.

The Hippo pathway plays crucial roles in governing various biological processes during tumorigenesis and metastasis. Within this pathway, upstream signaling stimuli activate a core kinase cascade, involving MST1/2 and LATS1/2, that subsequently phosphorylates and inhibits the transcriptional co-activators YAP and its paralog TAZ. This inhibition modulates the transcriptional regulation of downstream target genes, impacting cell proliferation, migration, and death. Despite the acknowledged significance of protein kinases in the Hippo pathway, the regulatory influence of protein phosphatases remains largely unexplored. In this study, we conducted the first gain-of-functional screen for protein tyrosine phosphatases (PTPs) regulating the Hippo pathway. Utilizing a LATS kinase biosensor (LATS-BS), a YAP/TAZ activity reporter (STBS-Luc), and a comprehensive PTP library, we identified numerous novel PTPs that play regulatory roles in the Hippo pathway. Subsequent experiments validated PTPN12, a master regulator of oncogenic receptor tyrosine kinases (RTKs), as a previously unrecognized negative regulator of the Hippo pathway effectors, oncogenic YAP/TAZ, influencing breast cancer cell proliferation and migration. In summary, our findings offer valuable insights into the roles of PTPs in the Hippo signaling pathway, significantly contributing to our understanding of breast cancer biology and potential therapeutic strategies.

Deregulation of small non-coding RNAs (sncRNAs) has been associated with the onset of metastasis. We evaluated the expression of sncRNAs in patients with early-stage breast cancer, performing RNA sequencing in 60 patients for whom tumor and sentinel lymph node (SLN) samples were available, and conducting differential expression, gene ontology, enrichment and survival analyses. Sequencing annotation classified most of the sncRNAs into small nucleolar RNA (snoRNAs, 70%) and small nuclear RNA (snRNA, 13%). Our results showed no significant differences in sncRNA expression between tumor or SLNs obtained from the same patient. Differential expression analysis showed down-regulation (n = 21) sncRNAs and up-regulation (n = 2) sncRNAs in patients with locoregional metastasis. The expression of SNHG5, SNORD90, SCARNA2 and SNORD78 differentiated luminal A from luminal B tumors, whereas SNORD124 up-regulation was associated with luminal B HER2+ tumors. Discriminating analysis and receiver-operating curve analysis revealed a signature of six snoRNAs (SNORD93, SNORA16A, SNORD113-6, SNORA7A, SNORA57 and SNORA18A) that distinguished patients with locoregional metastasis and predicted patient outcome. Gene ontology and Reactome pathway analysis showed an enrichment of biological processes associated with translation initiation, protein targeting to specific cell locations, and positive regulation of Wnt and NOTCH signaling pathways, commonly involved in the promotion of metastases. Our results point to the potential of several sncRNAs as surrogate markers of lymph node metastases and patient outcome in early-stage breast cancer patients. Further preclinical and clinical studies are required to understand the biological significance of the most significant sncRNAs and to validate our results in a larger cohort of patients.

The aim of this study was to investigate gene expression alterations associated with overall survival (OS) in glioblastoma (GBM). Using the Nanostring nCounter platform, we identified four genes (COL1A2, IGFBP3, NGFR, and WIF1) that achieved statistical significance when comparing GBM with non-neoplastic brain tissue. The four genes were included in a multivariate Cox Proportional Hazard model, along with age, extent of resection, and O6-methylguanine-DNA methyltransferase (MGMT) promotor methylation, to create a unique glioblastoma prognostic index (GPI). The GPI score inversely correlated with survival: patient with a high GPI had a median OS of 7.5 months (18-month OS = 9.7%) whereas patients with a low GPI had a median OS of 20.1 months (18-month OS = 54.5%; log rank p-value = 0.004). The GPI score was then validated in 188 GBM patients from The Cancer Genome Atlas (TCGA) from a national data base; similarly, patients with a high GPI had a median OS of 10.5 months (18-month OS = 12.4%) versus 16.9 months (18-month OS = 41.5%) for low GPI (log rank p-value = 0.0003). We conclude that this novel mRNA-based prognostic index could be useful in classifying GBM patients into risk groups and refine prognosis estimates to better inform treatment decisions or stratification into clinical trials.

Multiple sclerosis (MS) is the most common chronic inflammatory disease of the central nervous system. Currently, the pathological mechanisms of MS are not fully understood, but research has suggested that iron metabolism disorder may be associated with the onset and clinical manifestations of MS.
The study utilized publicly available databases and bioinformatics techniques for gene expression data analysis, including differential expression analysis, weighted correlation network analysis, gene enrichment analysis, and construction of logistic regression models. Subsequently, Mendelian randomization was used to assess the causal relationship between different iron metabolism markers and MS.
This study identified IREB2, LAMP2, ISCU, ATP6V1G1, ATP13A2, and SKP1 as genes associated with multiple sclerosis (MS) and iron metabolism, establishing their multi-gene diagnostic value for MS with an AUC of 0.83. Additionally, Mendelian randomization analysis revealed a potential causal relationship between transferrin saturation and MS (p=2.22E-02; OR 95%CI=0.86 (0.75, 0.98)), as well as serum transferrin and MS (p=2.18E-04; OR 95%CI=1.22 (1.10, 1.36)).
This study comprehensively explored the relationship between iron metabolism and MS through integrated bioinformatics analysis and Mendelian randomization methods. The findings provide important insights for further research into the role of iron metabolism disorder in the pathogenesis of MS and offer crucial theoretical support for the treatment of MS.

Long noncoding RNAs (lncRNAs) are functional bridges connecting the genome with phenotypes by interacting with DNA, mRNA, and proteins. Using publically available acute heat stress (AHS)-related RNA-seq data, we discovered novel lncRNAs and tested their association with AHS along with ~ 8800 known lncRNAs and ~ 28,000 mRNA transcripts. Our pipeline discovered a total of 145 potentially novel-lncRNAs. One of them (Fishcomb_p-value = 0.06) along with another novel transcript (annotated as protein-coding; Fishcomb_p-value = 0.03) were identified as significantly associated with AHS. We found five known-lncRNAs and 134 mRNAs transcripts that were significantly associated with AHS. Four novel lncRNAs interact cis-regulated with 12 mRNA transcripts and are targeted by 11 miRNAs. Also six meta-lncRNAs associate with 134 meta-mRNAs through trans-acting co-expression, each targeted by 15 and 216 miRNAs, respectively. Three of the known-lncRNAs significantly co-expressed with almost 97 of the significant mRNAs (Pearson correlation p-value < 0.05). We report the mentioned three known-lncRNAs (ENSGALT00000099876, ENSGALT00000107573, and ENSGALT00000106323) as the most, significantly regulatory elements of AHS in chicken. It can be concluded that in order to alleviate the adverse effects of AHS on chicken, the manipulation of the three regulatory lncRNAs could lead to a more desirable result than the manipulation of the most significant mRNAs.

BACKGROUND: The nuclear distribution E homologue 1 (NDE1) is a crucial dynein binding partner. The NDE1 protein has the potential to disrupt the normal functioning of centrosomes, leading to a compromised ability to generate spindles and ensure precise separation of chromosomes during cell division. The potential consequences of this phenomenon include genomic instability, malignant transformation and the proliferation of neoplastic growths. However, studies examining the connection between NDE1 and cancer is still very rare.
METHODS: The expression level, prognostic impact, gene change, DNA methylation, protein interaction, mRNA m6A modification, ceRNA network, associated gene and function enrichment, and immune-related effects of NDE1 in pan-cancer were examined using a range of online analytic tools and the R software package. The CCK-8 test, transwell assay, scratch assay and colony formation assay were used to confirm the effects of NDE1 on the proliferation, invasion and metastasis of bladder cancer cells.
RESULTS: Numerous tumour types have elevated NDE1, which is linked to a bad prognosis. NDE1 is an excellent diagnostic tool for many different types of cancer. Numerous malignancies have been linked to genetic changes in NDE1. NDE1 was connected to TMB, MSI, several immunological checkpoint genes and immune cell infiltration. NDE1 is linked to a number of immunological subtypes. NDE1 could affect how well immunotherapy works to treat different types of cancer. NDE1 was mostly associated with cell cycle, chromosomal segregation, DNA replication and mitotic segregation, according to GO and KEGG analyses. NDE1 physically binds to PAFAH1B1 and DCTN1, respectively. The proliferation, invasion and metastasis of bladder cancer cells may be prevented by NDE1 knockdown. Furthermore, knockdown of NDE1 promoted the apoptosis of bladder cancer cells.
CONCLUSIONS: High expression of NDE1 is present in a variety of tumours, which is linked to a bad prognosis for cancer. Knockdown of NDE1 inhibited the proliferation, invasion and metastasis of bladder cancer cells, and promoted the apoptosis. For a number of malignancies, NDE1 may be a biomarker for immunotherapy and prognosis.

Gestational diabetes mellitus (GDM), a transient disease, may lead to short- or long-term adverse influences on maternal and fetal health. Therefore, its potential functions, mechanisms and related molecular biomarkers must be comprehended for the control, diagnosis and treatment of GDM.
The differentially expressed genes (DEGs) were identified using GSE49524 and GSE87295 associated with GDM from the Gene Expression Omnibus database, followed by function enrichment analysis, protein-protein interactions network construction, hub DEGs mining, diagnostic value evaluation and immune infiltration analysis. Finally, hub DEGs, the strongest related to immune infiltration, were screened as immune-related biomarkers.
A hundred and seven DEGs were identified between patients with GDM and healthy individuals. Six hub genes with high diagnostic values, including ALDH1A1, BMP4, EFNB2, MME, PLAUR and SLIT2, were identified. Among these, two immune-related genes (PLAUR and SLIT2) with the highest absolute correlation coefficient were considered immune-related biomarkers in GDM.
Our study provides a comprehensive analysis of GDM, which would provide a foundation for the development of diagnosis and treatment of GDM.