Microarray Gene Expression Datasets Research Articles

Deep learning assisted cancer disease prediction from gene expression data using WT-GAN

Several diverse fields including the healthcare system and drug development sectors have benefited immensely through the adoption of deep learning (DL), which is a subset of artificial intelligence (AI) and machine learning (ML). Cancer makes up a significant percentage of the illnesses that cause early human mortality across the globe, and this situation is likely to rise in the coming years, especially when non-communicable illnesses are not considered. As a result, cancer patients would greatly benefit from precise and timely diagnosis and prediction. Deep learning (DL) has become a common technique in healthcare due to the abundance of computational power. Gene expression datasets are frequently used in major DL-based applications for illness detection, notably in cancer therapy. The quantity of medical data, on the other hand, is often insufficient to fulfill deep learning requirements. Microarray gene expression datasets are used for training procedures despite their extreme dimensionality, limited volume of data samples, and sparsely available information. Data augmentation is commonly used to expand the training sample size for gene data. The Wasserstein Tabular Generative Adversarial Network (WT-GAN) model is used for the data augmentation process for generating synthetic data in this proposed work. The correlation-based feature selection technique selects the most relevant characteristics based on threshold values. Deep FNN and ML algorithms train and classify the gene expression samples. The augmented data give better classification results (> 97%) when using WT-GAN for cancer diagnosis.

Gene-expression profile analysis to disclose diagnostics and therapeutics biomarkers for thyroid carcinoma

The most frequent endocrine cancer of the head and neck is thyroid carcinoma (THCA). Although there is increasing evidence linking THCA to genetic alterations, the exact molecular mechanism behind this relationship is not yet completely known to the researchers. There is still much to learn about THCA's molecular roots and genetic biomarkers. Though drug therapies are the best choice after metastasis, unfortunately, the majority of the patients progressively develop resistance against the therapeutic drugs after receiving them for a few years. Therefore, multi-targeted different variants of therapeutic drugs may be essential for effective treatment against THCA. To understand molecular mechanisms of THCA development and progression and explore multi-targeted different variants of therapeutic drugs, we detected 80 common differentially expressed genes (cDEGs) between THCA and non-THCA samples from six microarray gene expression datasets using the statistical LIMMA approach. Through protein-protein interaction (PPI) network analysis, we identified the top-ranked eight differentially expressed genes (TIMP1, FN1, THBS1, RUNX2, SHANK2, TOP2A, LRP2, and ACTN1) as the THCA-causing key genes (KGs), where 6 KGs (TIMP1, TOP2A, FN1, ACTN1, RUNX2, THBS1) are upregulated and 2 KGs (LRP2, SHANK2) are downregulated. The expression pattern analysis of KGs with the independent TCGA database by Box plots also confirmed their upregulated and downregulated patterns. The expression analysis of KGs in different stages of THCA development indicated that these KGs might be utilized as early diagnostic and prognostic biomarkers. The pan-cancer analysis of KGs indicated a substantial correlation of KGs with multiple cancers, including THCA. Some transcription factors (TFs) and microRNAs were detected as the key transcriptional and post-transcriptional regulators of KGs using gene regulatory network (GRN) analysis. The enrichment analysis of the cDEGs revealed several key molecular functions, biological processes, cellular components, and pathways significantly associated with THCA. These findings highlight critical mechanisms influenced by the identified key genes (KGs), providing deeper insight into their roles in THCA development. Then we detected 6 repurposable drug molecules (Entrectinib, Imatinib, Ponatinib, Sorafenib, Retevmo, and Pazopanib) by molecular docking with KGs-mediated receptor proteins, ADME/T analysis, and cross-validation with the independent receptors. Therefore, these findings might be useful resources for wet lab researchers and clinicians to consider an effective treatment strategy against THCA.

Computational insights into irinotecan's interaction with UBE2I in ovarian and endometrial cancers

Endometrial and Ovarian cancers are two highly prevalent and fatal reproductive diseases with poor prognoses among women. Elevated estrogen levels in Ovarian Cancer (OC) stimulate the endometrium, causing Endometrial Cancer (EC). Although numerous studies have reported the crucial genes and pathways in this cancer, the pathogenesis of this disease remains unclear. In this study, used bioinformatics tools to analyse GSE63678, GSE115810, GSE36389, GSE26712, GSE36668, GSE27651, GSE6008, GSE69429, GSE69428, GSE18521, GSE185209, GSE54388 gene expression microarray datasets for both the cancers. We analyzed the differential gene expression, functional association, and structural studies. The analysis identified crucial differentially expressed genes (DEGs) in both cancers associated with DNA damage, DNA integrity, and cell-cycle checkpoint signaling pathways. CLDN7, UBE2I, WT1, JAM2, FOXL2, F11R, JAM3, ZFPM2, MEF2C, and PIAS1 are the top 10 hub genes commonly identified in both cancer types. Only CLDN7 and F11R are upregulated, whereas the remaining hub genes are downregulated in both cancers, suggesting a common framework for contributing to tumorigenesis. Molecular docking and dynamics were performed on the UBE2I protein with Irinotecan Hydrochloride, which could serve as the new approach for treating and managing both cancers. The study reveals the common molecular pathways, pointing out the role of cell cycle and DNA damage and integrity checkpoint signaling in the pathogenesis of both cancer types. This study explored the UBE2I gene as a potential biomarker in OC and EC. Further, this study concludes that the irinotecan hydrochloride drug has higher therapeutic effects on UBE2I protein through docking and dynamics studies.

Screening and validating genes associated with cuproptosis in systemic lupus erythematosus by expression profiling combined with machine learning.

Cell death has long been a focal point in life sciences research, and recently, scientists have discovered a novel form of cell death induced by copper, termed cuproptosis. This paper aimed to identify genes associated with cuproptosis in Systemic Lupus Erythematosus (SLE) through machine learning, combined with single-cell sequencing (scRNA-seq), to screen and validate related genes. The analytical results were then experimentally verified. Two published microarray gene expression datasets (GSE65391 and GSE61635) from SLE and control peripheral blood samples were downloaded from the GEO database. The GSE65391 dataset was used as the training group, while the GSE61635 dataset served as the validation group. Differentially expressed genes (DEGs) from GSE65391 identified 12 differential genes. Nine diagnostic genes, considered potential biomarkers, were selected using the least absolute shrinkage and selection operator (LASSO) and support vector machine recursive feature elimination (SVM-RFE) analysis. The Receiver Operating Characteristic (ROC) curves for both the training and validation groups were used to calculate the Area Under the Curve (AUC) to assess discriminatory properties. CIBERSORT was used to assess the relationship between these diagnostic genes and a reference set of infiltrating immune cells. Single-cell RNA sequencing data (GSE162577) from SLE patients were also obtained from the GEO database and analyzed. Experimental validation of the most important SLE biomarkers was performed. Twelve significantly different cuproptosis-related genes were identified in the GSE65391 training set. Immune cell analysis revealed 12 immune cell types and identified nine signature genes, including PDHB, GLS, DLAT, LIAS, MTF1, DLST, DLD, LIPT1, and FDX1. In the GSE61635 validation set, seven genes were weakly expressed, and two genes were strongly expressed in the treatment group. According to the ROC curves, PDHB and GLS demonstrated significant diagnostic value. Additionally, correlation analysis was conducted on the nine characteristic genes in relation to immune infiltration. The distribution of key genes in immune cells was determined using single-cell RNA sequencing data. Finally, the mRNA expression of the nine diagnostic genes was validated using qPCR.

Elucidating gastric cancer mechanisms and therapeutic potential of Adociaquinone A targeting EGFR: A genomic analysis and Computer Aided Drug Design (CADD) approach.

Gastric cancer predominantly adenocarcinoma, accounts for over 85% of gastric cancer diagnoses. Current therapeutic options are limited, necessitating the discovery of novel drug targets and effective treatments. The Affymetrix gene expression microarray dataset (GSE64951) was retrieved from NCBI-GEO data normalization and DEGs identification was done by using R-Bioconductor package. Gene Ontology (GO) analysis of DEGs was performed using DAVID. The protein-protein interaction network was constructed by STRING database plugin in Cytoscape. Subclusters/modules of important interacting genes in main network were extracted by using MCODE. The hub genes from in the network were identified by using Cytohubba. The miRNet tool built a hub gene/mRNA-miRNA network and Kaplan-Meier-Plotter conducted survival analysis. AutoDock Vina and GROMACS MD simulations were used for docking and stability analysis of marine compounds against the 5CNN protein. Total 734 DEGs (507 up-regulated and 228 down-regulated) were identified. Differentially expressed genes (DEGs) were enriched in processes like cell-cell adhesion and ATP binding. Eight hub genes (EGFR, HSPA90AA1, MAPK1, HSPA4, PPP2CA, CDKN2A, CDC20, and ATM) were selected for further analysis. A total of 23 miRNAs associated with hub genes were identified, with 12 of them targeting PPP2CA. EGFR displayed the highest expression and hazard rate in survival analyses. The kinase domain of EGFR (PDBID: 5CNN) was chosen as the drug target. Adociaquinone A from Petrosia alfiani, docked with 5CNN, showed the lowest binding energy with stable interactions across a 50 ns MD simulation, highlighting its potential as a lead molecule against EGFR. This study has identified crucial DEGs and hub genes in gastric cancer, proposing novel therapeutic targets. Specifically, Adociaquinone A demonstrates promising potential as a bioactive drug against EGFR in gastric cancer, warranting further investigation. The predicted miRNA against the hub gene/proteins can also be used as potential therapeutic targets.

Screening of differential gene expression patterns through survival analysis for diagnosis, prognosis and therapies of clear cell renal cell carcinoma.

Clear cell renal cell carcinoma (ccRCC) is the most prevalent subtype of kidney cancer. Although there is increasing evidence linking ccRCC to genetic alterations, the exact molecular mechanism behind this relationship is not yet completely known to the researchers. Though drug therapies are the best choice after the metastasis, unfortunately, the majority of the patients progressively develop resistance against the therapeutic drugs after receiving it for almost 2 years. In this case, multi-targeted different variants of therapeutic drugs are essential for effective treatment against ccRCC. To understand molecular mechanisms of ccRCC development and progression, and explore multi-targeted different variants of therapeutic drugs, it is essential to identify ccRCC-causing key genes (KGs). In order to obtain ccRCC-causing KGs, at first, we detected 133 common differentially expressed genes (cDEGs) between ccRCC and control samples based on nine (9) microarray gene-expression datasets with NCBI accession IDs GSE16441, GSE53757, GSE66270, GSE66272, GSE16449, GSE76351, GSE66271, GSE71963, and GSE36895. Then, we filtered these cDEGs through survival analysis with the independent TCGA and GTEx database and obtained 54 scDEGs having significant prognostic power. Next, we used protein-protein interaction (PPI) network analysis with the reduced set of 54 scDEGs to identify ccRCC-causing top-ranked eight KGs (PLG, ENO2, ALDOB, UMOD, ALDH6A1, SLC12A3, SLC12A1, SERPINA5). The pan-cancer analysis with KGs based on TCGA database showed the significant association with different subtypes of kidney cancers including ccRCC. The gene regulatory network (GRN) analysis revealed some crucial transcriptional and post-transcriptional regulators of KGs. The scDEGs-set enrichment analysis significantly identified some crucial ccRCC-causing molecular functions, biological processes, cellular components, and pathways that are linked to the KGs. The results of DNA methylation study indicated the hypomethylation and hyper-methylation of KGs, which may lead the development of ccRCC. The immune infiltrating cell types (CD8+ T and CD4+ T cell, B cell, neutrophil, dendritic cell and macrophage) analysis with KGs indicated their significant association in ccRCC, where KGs are positively correlated with CD4+ T cells, but negatively correlated with the majority of other immune cells, which is supported by the literature review also. Then we detected 10 repurposable drug molecules (Irinotecan, Imatinib, Telaglenastat, Olaparib, RG-4733, Sorafenib, Sitravatinib, Cabozantinib, Abemaciclib, and Dovitinib.) by molecular docking with KGs-mediated receptor proteins. Their ADME/T analysis and cross-validation with the independent receptors, also supported their potent against ccRCC. Therefore, these outputs might be useful inputs/resources to the wet-lab researchers and clinicians for considering an effective treatment strategy against ccRCC.

Identification of potential new T cell activation molecules: a Bioinformatic Approach.

T-cell activation is central for the initiation of T cell mediated adaptive immune response and is the result of the close communication between the Antigen Presenting Cell (APC) and the T lymphocyte. Although T-cell activation is currently well understood, and many intracellular pathways are well characterized, nevertheless new players are constantly identified, and this complements the known protein interactome. In this work we aimed to identify new proteins involved in T cell activation. We reviewed and analyzed results of microarray gene expression datasets reported in the public database GEO-NCBI. Using data from GSE136625, GSE50971, GSE13887, GSE11989 and GSE902 we performed different comparisons using R and other bioinformatic tools including GEO2R and we report here upregulated genes that have no previous reports in immune related functions and with potential participation upon T-cell activation. Our results indicate that RND3, SYT10, IgSF6 and PIN1 are potential new T-cell activation molecules.

Improving data interpretability with new differential sample variance gene set tests.

Gene set analysis methods have played a major role in generating biological interpretations from omics data such as gene expression datasets. However, most methods focus on detecting homogenous pattern changes in mean expression and methods detecting pattern changes in variance remain poorly explored. While a few studies attempted to use gene-level variance analysis, such approach remains under-utilized. When comparing two phenotypes, gene sets with distinct changes in subgroups under one phenotype are overlooked by available methods although they reflect meaningful biological differences between two phenotypes. Multivariate sample-level variance analysis methods are needed to detect such pattern changes. We use ranking schemes based on minimum spanning tree to generalize the Cramer-Von Mises and Anderson-Darling univariate statistics into multivariate gene set analysis methods to detect differential sample variance or mean. We characterize these methods in addition to two methods developed earlier using simulation results with different parameters. We apply the developed methods to microarray gene expression dataset of prednisolone-resistant and prednisolone-sensitive children diagnosed with B-lineage acute lymphoblastic leukemia and bulk RNA-sequencing gene expression dataset of benign hyperplastic polyps and potentially malignant sessile serrated adenoma/polyps. One or both of the two compared phenotypes in each of these datasets have distinct molecular subtypes that contribute to heterogeneous differences. Our results show that methods designed to detect differential sample variance are able to detect specific hallmark signaling pathways associated with the two compared phenotypes as documented in available literature. The results in this study demonstrate the usefulness of methods designed to detect differential sample variance in providing biological interpretations when biologically relevant but heterogeneous changes between two phenotypes are prevalent in specific signaling pathways. Software implementation of the developed methods is available with detailed documentation from Bioconductor package GSAR. The available methods are applicable to gene expression datasets in a normalized matrix form and could be used with other omics datasets in a normalized matrix form with available collection of feature sets.

A Bioinformatics Investigation of Hub Genes Involved in Treg Migration and Its Synergistic Effects, Using Immune Checkpoint Inhibitors for Immunotherapies.

This study aimed to identify hub genes involved in regulatory T cell (Treg) function and migration, offering insights into potential therapeutic targets for cancer immunotherapy. We performed a comprehensive bioinformatics analysis using three gene expression microarray datasets from the GEO database. Differentially expressed genes (DEGs) were identified to pathway enrichment analysis to explore their functional roles and potential pathways. A protein-protein interaction network was constructed to identify hub genes critical for Treg activity. We further evaluated the co-expression of these hub genes with immune checkpoint proteins (PD-1, PD-L1, CTLA4) and assessed their prognostic significance. Through this comprehensive analysis, we identified CCR8 as a key player in Treg migration and explored its potential synergistic effects with ICIs. Our findings suggest that CCR8-targeted therapies could enhance cancer immunotherapy outcomes, with breast invasive carcinoma (BRCA) emerging as a promising indication for combination therapy. This study highlights the potential of CCR8 as a biomarker and therapeutic target, contributing to the development of targeted cancer treatment strategies.

Gene selection based on recursive spider wasp optimizer guided by marine predators algorithm

Detecting tumors using gene analysis in microarray data is a critical area of research in artificial intelligence and bioinformatics. However, due to the large number of genes compared to observations, feature selection is a central process in microarray analysis. While various gene selection methods have been developed to select the most relevant genes, these methods’ efficiency and reliability can be improved. This paper proposes a new two-phase gene selection method that combines the ReliefF filter method with a novel version of the spider wasp optimizer (SWO) called RSWO-MPA. In the first phase, the ReliefF filter method is utilized to reduce the number of genes to a reasonable number. In the second phase, RSWO-MPA applies a recursive spider wasp optimizer guided by the marine predators algorithm (MPA) to select the most informative genes from the previously selected ones. The MPA is used in the initialization step of recursive SWO to narrow down the search space to the most relevant and accurate genes. The proposed RSWO-MPA has been implemented and validated through extensive experimentation using eight microarray gene expression datasets. The enhanced RSWO-MPA is compared with seven widely used and recently developed meta-heuristic algorithms, including Kepler optimization algorithm (KOA), marine predators algorithm (MPA), social ski-driver optimization (SSD), whale optimization algorithm (WOA), Harris hawks optimization (HHO), artificial bee colony (ABC) algorithm, and original SWO. The experimental results demonstrate that the developed method yields the highest accuracy, selects fewer features, and exhibits more stability than other compared algorithms and cutting-edge methods for all the datasets used. Specifically, it achieved an accuracy of 100.00%, 94.51%, 98.13%, 95.63%, 100.00%, 100.00%, 92.97%, and 100.00% for Yeoh, West, Chiaretti, Burcyznski, leukemia, ovarian cancer, central nervous system, and SRBCT datasets, respectively.

Retraction Note: Investigation of molecular biomarker candidates for diagnosis and prognosis of chronic periodontitis by bioinformatics analysis of pooled microarray gene expression datasets in Gene Expression Omnibus (GEO).

Retraction Note: Investigation of molecular biomarker candidates for diagnosis and prognosis of chronic periodontitis by bioinformatics analysis of pooled microarray gene expression datasets in Gene Expression Omnibus (GEO).

Abstract PO5-15-10: Plasticity marker FOXC1 expression accurately predicts efficacy of Adjuvant Tamoxifen + Chemotherapy in reducing all-cause mortality in ER+LN- Breast Cancer: Validation in the SCAN-B Prospective Study (NCT02306096)

Abstract INTRODUCTION: Despite the development and validation of multimarker gene panels like OncotypeDx®, Mammaprint® and EndoPredict®, an affordable and globally accessible biomarker approach that can predict increased risk of all-cause mortality in patients diagnosed with ER+LN- breast cancer, as well as the efficacy of adjuvant endocrine + chemotherapy in mitigating such elevated risk, continues to be an unmet medical need. Lack of such a pragmatic solution continues to be the cause of preventable recurrent/metastatic disease in patients diagnosed with ER+LN- breast cancer in resource-challenged settings around the world. We hypothesized that a predictive biomarker strategy that measures expression of the plasticity marker FOXC1, in combination with clinical parameters like tumor size (TS) and tumor grade (TG), may offer equivalent results to that of the above multimarker gene panels at a significantly lower economic cost without compromising accuracy of prediction results. METHODS: Pre-treatment tumor RNA data obtained from a training cohort (compendium of gene expression microarray datasets, n=2857) and a single large validation cohort (SCAN-B Prospective Multicenter Observational Study, n=3520) for patients diagnosed with ER+LN- breast cancer were analyzed for FOXC1 expression, tumor size (TS) and tumor grade (TG) and correlated with Recurrence-Free Survival (RFS) and Overall Survival (OS). Optimized biomarker cutoff values based on model area-under-curve were leave-one-out cross validated and Risk-of-Recurrence (ROR) prediction algorithm derived utilizing the (compendium) training dataset. The unmodified strategy was then validated in the independent (SCAN-B) validation dataset. RESULTS: A predetermined High ROR score calculated using FOXC1 expression, TS and TG (trained using the compendium dataset) predicted efficacy of adjuvant endocrine + chemotherapy over that of adjuvant endocrine therapy alone in ER+LN- patients in terms of statistically significant reduction in all-cause mortality at 8-years post-diagnosis in the SCAN-B validation dataset (4.66% vs 24.06%, n=326, OR=6.48, 95%CI [2.97-14.11], p< 0.0001). OncotypeDx® (5.7% vs 14.4%, n=762, OR=2.79, 95%CI [1.37-5.67], p=0.002), Mammaprint® (5.1% vs 18.4%, n=665, OR= 95%CI [2.04-8.43], p< 0.0001) and Endopredict® (4.8% vs 15.2%, n=923, OR=3.58 95%CI [1.78-7.21], p=0.0002) were also statistically significant predictors of the same. CONCLUSION: Pre-treatment tumor FOXC1 mRNA or protein expression (assessed using qRT-PCR or routine immunohistochemistry (IHC), respectively, when combined with TS and TG presents a unique and economical alternative solution to multimarker gene panel tests like OncotypeDx®, Mammaprint® or Endopredict®, for guiding therapy of patients diagnosed with ER+LN- breast cancer in resource challenged settings. Such an approach to identify elevated risk of recurrence in patients diagnosed with ER+LN- breast cancer and prevent the same by guiding adjuvant endocrine + chemotherapy decisions, could help to extend recurrence-free and overall survival. Such an approach merits testing in real world ER+LN- patient cohorts in resource-challenged settings to help support implementation of this FOXC1-driven predictive biomarker strategy in the clinic. Citation Format: Partha Ray, Tania Ray, Clive Taylor, Robert Hussa. Plasticity marker FOXC1 expression accurately predicts efficacy of Adjuvant Tamoxifen + Chemotherapy in reducing all-cause mortality in ER+LN- Breast Cancer: Validation in the SCAN-B Prospective Study (NCT02306096) [abstract]. In: Proceedings of the 2023 San Antonio Breast Cancer Symposium; 2023 Dec 5-9; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2024;84(9 Suppl):Abstract nr PO5-15-10.

A 50-gene biomarker identifies estrogen receptor-modulating chemicals in a microarray compendium

High throughput transcriptomics (HTTr) profiling has the potential to rapidly and comprehensively identify molecular targets of environmental chemicals that can be linked to adverse outcomes. We describe here the construction and characterization of a 50-gene expression biomarker designed to identify estrogen receptor (ER) active chemicals in HTTr datasets. Using microarray comparisons, the genes in the biomarker were identified as those that exhibited consistent directional changes when ER was activated (4 ER agonists; 4 ESR1 gene constitutively active mutants) and opposite directional changes when ER was suppressed (4 antagonist treatments; 4 ESR1 knockdown experiments). The biomarker was evaluated as a predictive tool using the Running Fisher algorithm by comparison to annotated gene expression microarray datasets including those evaluating the transcriptional effects of hormones and chemicals in MCF-7 cells. Depending on the reference dataset used, the biomarker had a predictive accuracy for activation of up to 96%. To demonstrate applicability for HTTr data analysis, the biomarker was used to identify ER activators in a set of 15 chemicals that are considered potential bisphenol A (BPA) alternatives examined at up to 10 concentrations in MCF-7 cells and analyzed by full-genome TempO-Seq. Using benchmark dose (BMD) modeling, the biomarker genes stratified the ER potency of BPA alternatives consistent with previous studies. These results demonstrate that the ER biomarker can be used to accurately identify ER activators in transcript profile data derived from MCF-7 cells.

LXRα Promotes Abdominal Aortic Aneurysm Formation Through UHRF1 Epigenetic Modification of miR-26b-3p.

Abdominal aortic aneurysm (AAA) is a severe aortic disease without effective pharmacological approaches. The nuclear hormone receptor LXRα (liver X receptor α), encoded by the NR1H3 gene, serves as a critical transcriptional mediator linked to several vascular pathologies, but its role in AAA remains elusive. Through integrated analyses of human and murine AAA gene expression microarray data sets, we identified NR1H3 as a candidate gene regulating AAA formation. To investigate the role of LXRα in AAA formation, we used global Nr1h3-knockout and vascular smooth muscle cell-specific Nr1h3-knockout mice in 2 AAA mouse models induced with angiotensin II (1000 ng·kg·min; 28 days) or calcium chloride (CaCl2; 0.5 mol/L; 42 days). Upregulated LXRα was observed in the aortas of patients with AAA and in angiotensin II- or CaCl2-treated mice. Global or vascular smooth muscle cell-specific Nr1h3 knockout inhibited AAA formation in 2 mouse models. Loss of LXRα function prevented extracellular matrix degeneration, inflammation, and vascular smooth muscle cell phenotypic switching. Uhrf1, an epigenetic master regulator, was identified as a direct target gene of LXRα by integrated analysis of transcriptome sequencing and chromatin immunoprecipitation sequencing. Susceptibility to AAA development was consistently enhanced by UHRF1 (ubiquitin-like containing PHD and RING finger domains 1) in both angiotensin II- and CaCl2-induced mouse models. We then determined the CpG methylation status and promoter accessibility of UHRF1-mediated genes using CUT&Tag (cleavage under targets and tagmentation), RRBS (reduced representation bisulfite sequencing), and ATAC-seq (assay for transposase-accessible chromatin with sequencing) in vascular smooth muscle cells, which revealed that the recruitment of UHRF1 to the promoter of miR-26b led to DNA hypermethylation accompanied by relatively closed chromatin states, and caused downregulation of miR-26b expression in AAA. Regarding clinical significance, we found that underexpression of miR-26b-3p correlated with high risk in patients with AAA. Maintaining miR-26b-3p expression prevented AAA progression and alleviated the overall pathological process. Our study reveals a pivotal role of the LXRα/UHRF1/miR-26b-3p axis in AAA and provides potential biomarkers and therapeutic targets for AAA.

Investigation of Common Pathways and Putative Biomarker Candidates of Colorectal Cancer and Insomnia by Using Integrative In-Silico Approaches.

Colorectal cancer (CRC) is one of the leading causes of cancer-related mortalities across the globe. Accumulating evidence shows that individuals having sleep disorders such as insomnia are at high risk of developing CRC, yet the association of sleep disorders with CRC risk is still unclear. Here, we investigated the potential molecular connections between CRC and insomnia using integrative in silico approaches. This study aims to explore the potential molecular connections between CRC and insomnia utilizing integrative in-silico methodologies. Gene expression microarray datasets for CRC and insomnia samples were retrieved from the NCBI-GEO database and analyzed using R. Functional enrichment analysis of common differentially expressed genes (DEGs) was performed by the g: Profiler tool. Cytoscape software was used to construct a protein-protein interaction network and hub gene identification. Expression profiles of hub genes in TCGA datasets were also determined, and predicted miRNAs targeting hub genes were analyzed by miRNA target prediction tools. Our results revealed a total of 113 shared DEGs between the CRC and insomnia datasets. Six genes (HSP8A, GAPDH, HSP90AA1, EEF1G, RPS6, and RPLP0), which were also differently expressed in TCGA datasets, were prioritized as hub genes and were found to be enriched in pathways related to protein synthesis. hsa-miR-324-3p, hsa-miR-769-3p, and hsa-miR-16-5p were identified as promising miRNA biomarkers for two diseases. Our in-silico analysis provides promising evidence of the molecular link between CRC and insomnia and highlights multiple potential molecular biomarkers and pathways. Validation of the results by wet lab work can be utilized for novel translational and precision medicine applications to alleviate the public health burden of CRC.

Abstract 6209: From data disparity to data harmony: A comprehensive pan-cancer omics data collection

Abstract In cancer research, the exponential growth of omics datasets offers a significant opportunity for scientific advancement. However, challenges such as the lack of uniform standards, in both clinical and omic data, hinder the effective utilization of these datasets, thus impeding our understanding of cancer biology and the development of innovative therapeutic approaches.Addressing these challenges, we have created a novel collection of pan-cancer omics datasets with extensive clinical data harmonization and consistent omic data normalization.Here, we focused on patient-derived gene expression microarray datasets from the Gene Expression Omnibus database. To navigate the complexities presented by the diverse clinical descriptions inherent in these datasets, we leveraged our proprietary ontology, machine learning models, and domain expert quality control processes to homogenize the clinical data elements.Datasets were then selected based on sample composition, molecular data compatibility, and clinical data availability, then passed through a uniform preprocessing and normalization pipeline to maximize data quality. Finally, gene names were aligned on a single annotation reference, and potential batch effects were adjusted before expression data were merged together.We obtained a total of 32,825 transcriptomic sample profiles from 470 datasets, covering 13,435 genes and 45 clinical data elements, across 30 cancer types. Healthy tissue was favored over adjacent tissue, to minimize the risk of introducing biases related to cancer patient background genomic profiles into downstream analyses. We compared our collection with The Cancer Genome Atlas (TCGA), the most commonly used RNA-seq transcriptomic dataset in cancer research. It covers 30 out of the 33 TCGA cancer types, with on average 4.2 times more samples per cancer type ([0.3; 45.5], median 3.4). Despite the two data collections being based on distinct technologies, we observed a Pearson correlation of 0.69 over the 11,753 genes in common, and a 100% overlap of the differentially expressed genes between genders. This consistency highlights cross-technology reliability and complementarity.We have built and continuously enriched a comprehensive dataset collection enabling the secondary analysis of high-quality omic data. This initial work - focused on microarray datasets - allows us to streamline design, exploration and validation of various omics data-driven studies in cancer research.Our ongoing efforts involve not only the continued integration of microarray datasets but also the integration of pan-cancer RNA-seq and single-cell data. This initiative is set to expand further, encompassing a broader range of omics datasets in the future. Citation Format: Lea Meunier, Guillaume Appe, Abdelkader Behdenna, Valentin Bernu, Helia Brull Corretger, Prashant Dhillon, Eleonore Fox, Julien Haziza, Charles Lescure, Camille Marijon, Clemence Petit, Solene Weill, Akpeli Nordor. From data disparity to data harmony: A comprehensive pan-cancer omics data collection [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 6209.

Differentially Expressed Genes in Atopic Dermatitis: Bioinformatics analysis of pooled microarray gene expression datasets in Gene Expression Omnibus

Introduction:
 Atopic dermatitis (AD) is a chronic and refractory inflammatory skin disease characterized by relapsing eczematous and pruritic skin lesions. Understanding the specific gene expression patterns associated with AD is crucial for advancing diagnosis and targeted treatment development. Using bioinformatics methods, candidate genes and biological pathways involved in AD pathogenesis were identified based on gene expression profiles in the Gene Expression Omnibus (GEO) database.
 Materials and Methods:
 A comprehensive analysis of four pooled transcriptomic datasets obtained from the Gene Expression Omnibus (GEO) database were conducted. Differential gene expression analysis was performed using the GEO2R. The differentially expressed genes (DEGs) between lesion skin of AD patients and normal skin of individuals were analyzed using the Gene Ontology (GO) term enrichment analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway, and protein-protein interaction (PPI) network to explore the functional roles of these genes.
 Results:
 Among the patient-level gene expression datasets, we identified 133 shared DEGs, consisting of 48 upregulated genes and 85 downregulated genes. GO analyses revealed these DEGs to be significantly enriched in biological processes including inflammatory responses, cytokine-mediated signaling pathway. These DEGs were also enriched in the KEGG pathway, including viral protein interaction with cytokine and cytokine receptor, C-type lectin receptor signaling pathway, cytokine-cytokine receptor interaction, JAK-STAT signaling pathway, and Adipocytokine signaling pathway.
 Conclusion:
 By comparing with other studies using the same method, we found that in addition to the already confirmed pathways such as inflammatory response, different studies have found changes in different hub genes and metabolic pathways, which prompts us to develop individualized treatments for AD.

Abstract B012: Pedigree and multi-omics analyses implicate nucleotide excision repair genes, XPA, ERCC6, and ERCC8, in susceptibility to breast cancer

Abstract Background: The majority of low- and moderate-risk [odds ratios (OR) ≤ 3] susceptibility genes for peri- and post-menopausal breast cancer (≥50 years of age), particularly those with minor allele frequency (MAF) of <30%, have not yet been identified and/or validated. We used a hypothesis-driven integrative approach to identify breast cancer susceptibility genes by analyzing pedigree and multi-omics (genomics, transcriptomics, mutation landscape) data. Methods: We conducted pedigree analysis in 31 families with DNA repair disorders, xeroderma pigmentosum (XP), trichothiodystrophy (TTD), Cockayne syndrome (CS), XP/TTD and XP/CS. These disorders can be caused by mutations in genes involved in the nucleotide excision repair (NER) pathway such as XPA-XPG, ERCC6 and ERCC8. We estimated the frequency of hereditary breast cancer, defined as the presence of at least two cases of premenopausal breast cancer (diagnosed <50 years of age) or one case of premenopausal breast cancer and one case of ovarian cancer in one side of the family, in our study. Furthermore, we conducted individual analysis of four genome-wide association studies (GWAS) datasets (totaling 2961 cases and 4335 controls), individual and meta-analysis of two gene expression microarray datasets (98 cases and 163 controls), and targeted analysis of The Cancer Genome Atlas (TCGA) (522 breast cancers). All analyses involved invasive ductal breast carcinoma in women diagnosed at ≥50 years of age. Results: In an XP family with causative mutations in XPA, there were four cases of premenopausal breast cancer (including heterozygote mutation carriers), one case of ovarian cancer, and one case of prostate cancer. Comprehensive sequencing failed to identify a BRCA1 or BRCA2 mutation in this family. Overall, 7 of 31 (22.6%) families with DNA repair disorders had a family history that was compatible with hereditary breast cancer. Hence, we selected a panel of NER genes as candidates for association studies with breast cancer. Our candidate-gene association analysis of GWAS datasets suggested an increased risk of breast cancer with ERCC6 (main effect:1.29 ≤ OR ≤ 2.91, 0.005 ≤ p ≤ 0.04, 11.8 ≤ MAF ≤ 29.7%) and implicated its interaction with ERCC8 (joint effect:3.03 ≤ OR ≤ 5.31, 0.01 ≤ pinteraction ≤ 0.03). Our transcriptomics analysis identified significant upregulation of ERCC6 (p=7.95 × 10–6) and ERCC8 (p=4.67 × 10–6) in breast cancer. Our analysis of TCGA breast tumor data showed similar mutation frequency of ERCC6 (1.8%) and ERCC8 (0.3%) as known breast cancer susceptibility genes such as BLM (1.9%) and LSP1 (0.3%). Conclusion: Our integrative analyses suggest that ERCC6 may be a previously unreported low- to moderate-risk breast cancer susceptibility gene, which may also interact with ERCC8. ERCC6 and ERCC8 code for the main components of the transcription-coupled (TC) repair sub-pathway of the NER. TC-NER protects against genotoxic agents, such as free radicals and reactive oxygen species. Our results provide a deeper biologic insight and leads for prevention of common breast cancer. Citation Format: Roxana Moslehi, Amiran Dzutsev. Pedigree and multi-omics analyses implicate nucleotide excision repair genes, XPA, ERCC6, and ERCC8, in susceptibility to breast cancer [abstract]. In: Proceedings of the AACR Special Conference in Cancer Research: DNA Damage Repair: From Basic Science to Future Clinical Application; 2024 Jan 9-11; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2024;84(1 Suppl):Abstract nr B012.

CFS‐MOES Ensemble Model on Metaheuristic Search‐Based Feature Selection

Cancer is one of the leading causes of death across the globe. There is a need for early diagnosis to improve the chance of successful treatment and reduce the mortality associated with cancer. Due to the availability of highly specialized cancer datasets, molecular classification of cancer by gene expression, machine learning, and deep learning, a part of artificial intelligence (AI) techniques is used in detecting the disease. The application of several classification and feature selection methods on microarray gene expression datasets helps learn models that are able to predict a given disease. However, the tremendous dimensionality of the microarray cancer dataset is the greatest challenge in interpreting the data. In this work, the optimal feature subsets are selected by combining the correlation‐based feature selection (CFS) technique with five distinct meta‐heuristic search methods: evolutionary search (ES), particle swarm optimization search (PSOS), genetic search (GS), harmony search (HS), and multiobject evolutionary search (MOES). Furthermore, a CFS‐MOES (correlation‐based feature selection—multiobject evolutionary search) ensemble model is proposed based on a majority voting mechanism to improve the classification performance. Six microarray cancer datasets are considered, and seven traditional classifiers are evaluated on those datasets. Three classifiers, namely, K‐nearest neighbour (KNN), multilayer perceptron (MLP), and random forest (RF), were chosen as the base classifiers based on their F‐measure score. The features chosen by our proposed CFS‐MOES method significantly improve the accuracy of the proposed model. Moreover, the proposed model has also been compared with the other ensemble models generated using CFS‐ES (correlation‐based feature selection —evolutionary search), CFS‐PSOS (correlation‐based feature selection—particle swarm optimization search), CFS‐GS (correlation‐based feature selection—genetic search), and CFS‐HS (correlation‐based feature selection—harmony search) feature selection methods, ensuring better classification accuracy with a reduced feature subset. This model is also evaluated using significant parameters such as precision, recall, F‐measure, accuracy, Matthews correlation coefficient (MCC), and mean absolute error (MAE). According to the experimental results, our proposed model has a remarkable accuracy of 98.83% for breast cancer and 98.79% for cervical cancer.

Microarray Gene Expression Dataset Feature Selection and Classification with Swarm Optimization to Diagnosis Diseases

Microarray Gene Expression Dataset Feature Selection and Classification with Swarm Optimization to Diagnosis Diseases