CuBlock: a cross-platform normalization method for gene-expression microarrays

Transcriptomic Analysis of Hub Genes Reveals Associated Inflammatory Pathways in Estrogen-Dependent Gynecological Diseases

Gynecological diseases are triggered by aberrant molecular pathways that alter gene expression, hormonal balance, and cellular signaling pathways, which may lead to long-term physiological consequences. This study was able to identify highly preserved modules and key hub genes that are mainly associated with gynecological diseases, represented by endometriosis (EM), ovarian cancer (OC), cervical cancer (CC), and endometrial cancer (EC), through the weighted gene co-expression network analysis (WGCNA) of microarray datasets sourced from the Gene Expression Omnibus (GEO) database. Five highly preserved modules were observed across the EM (GSE51981), OC (GSE63885), CC (GSE63514), and EC (GSE17025) datasets. The functional annotation and pathway enrichment analysis revealed that the highly preserved modules were heavily involved in several inflammatory pathways that are associated with transcription dysregulation, such as NF-kB signaling, JAK-STAT signaling, MAPK-ERK signaling, and mTOR signaling pathways. Furthermore, the results also include pathways that are relevant in gynecological disease prognosis through viral infections. Mutations in the ESR1 gene that encodes for ERα, which were shown to also affect signaling pathways involved in inflammation, further indicate its importance in gynecological disease prognosis. Potential drugs were screened through the Drug Repurposing Encyclopedia (DRE) based on the up-and downregulated hub genes, wherein a bacterial ribosomal subunit inhibitor and a benzodiazepine receptor agonist were the top candidates. Other drug candidates include a dihydrofolate reductase inhibitor, glucocorticoid receptor agonists, cholinergic receptor agonists, selective serotonin reuptake inhibitors, sterol demethylase inhibitors, a bacterial antifolate, and serotonin receptor antagonist drugs which have known anti-inflammatory effects, demonstrating that the gene network highlights specific inflammatory pathways as a therapeutic avenue in designing drug candidates for gynecological diseases.

Modular Hub Genes in DNA Microarray Suggest Potential Signaling Pathway Interconnectivity in Various Glioma Grades

Gliomas have displayed significant challenges in oncology due to their high degree of invasiveness, recurrence, and resistance to treatment strategies. In this work, the key hub genes mainly associated with different grades of glioma, which were represented by pilocytic astrocytoma (PA), oligodendroglioma (OG), anaplastic astrocytoma (AA), and glioblastoma multiforme (GBM), were identified through weighted gene co-expression network analysis (WGCNA) of microarray datasets retrieved from the Gene Expression Omnibus (GEO) database. Through this, four highly correlated modules were observed to be present across the PA (GSE50161), OG (GSE4290), AA (GSE43378), and GBM (GSE36245) datasets. The functional annotation and pathway enrichment analysis done through the Database for Annotation, Visualization, and Integrated Discovery (DAVID) showed that the modules and hub genes identified were mainly involved in signal transduction, transcription regulation, and protein binding, which collectively deregulate several signaling pathways, mainly PI3K/Akt and metabolic pathways. The involvement of several hub genes primarily linked to other signaling pathways, including the cAMP, MAPK/ERK, Wnt/β-catenin, and calcium signaling pathways, indicates potential interconnectivity and influence on the PI3K/Akt pathway and, subsequently, glioma severity. The Drug Repurposing Encyclopedia (DRE) was used to screen for potential drugs based on the up- and downregulated hub genes, wherein the synthetic progestin hormones norgestimate and ethisterone were the top drug candidates. This shows the potential neuroprotective effect of progesterone against glioma due to its influence on EGFR expression and other signaling pathways. Aside from these, several experimental and approved drug candidates were also identified, which include an adrenergic receptor antagonist, a PPAR-γ receptor agonist, a CDK inhibitor, a sodium channel blocker, a bradykinin receptor antagonist, and a dopamine receptor agonist, which further highlights the gene network as a potential therapeutic avenue for glioma.

A quantitative systems pharmacology model for certolizumab pegol treatment in moderate-to-severe psoriasis

Background

Psoriasis is a chronic immune-mediated inflammatory systemic disease with skin manifestations characterized by erythematous, scaly, itchy and/or painful plaques resulting from hyperproliferation of keratinocytes. Certolizumab pegol [CZP], a PEGylated antigen binding fragment of a humanized monoclonal antibody against TNF-alpha, is approved for the treatment of moderate-to-severe plaque psoriasis. Patients with psoriasis present clinical and molecular variability, affecting response to treatment. Herein, we utilized an in silico approach to model the effects of CZP in a virtual population (vPop) with moderate-to-severe psoriasis. Our proof-of-concept study aims to assess the performance of our model in generating a vPop and defining CZP response variability based on patient profiles.

Methods

We built a quantitative systems pharmacology (QSP) model of a clinical trial-like vPop with moderate-to-severe psoriasis treated with two dosing schemes of CZP (200 mg and 400 mg, both every two weeks for 16 weeks, starting with a loading dose of CZP 400 mg at weeks 0, 2, and 4). We applied different modelling approaches: (i) an algorithm to generate vPop according to reference population values and comorbidity frequencies in real-world populations; (ii) physiologically based pharmacokinetic (PBPK) models of CZP dosing schemes in each virtual patient; and (iii) systems biology-based models of the mechanism of action (MoA) of the drug.

Results

The combination of our different modelling approaches yielded a vPop distribution and a PBPK model that aligned with existing literature. Our systems biology and QSP models reproduced known biological and clinical activity, presenting outcomes correlating with clinical efficacy measures. We identified distinct clusters of virtual patients based on their psoriasis-related protein predicted activity when treated with CZP, which could help unravel differences in drug efficacy in diverse subpopulations. Moreover, our models revealed clusters of MoA solutions irrespective of the dosing regimen employed.

Conclusion

Our study provided patient specific QSP models that reproduced clinical and molecular efficacy features, supporting the use of computational methods as modelling strategy to explore drug response variability. This might shed light on the differences in drug efficacy in diverse subpopulations, especially useful in complex diseases such as psoriasis, through the generation of mechanistically based hypotheses.

Uniformly shaped harmonization combines human transcriptomic data from different platforms while retaining their biological properties and differential gene expression patterns

Introduction: Co-normalization of RNA profiles obtained using different experimental platforms and protocols opens avenue for comprehensive comparison of relevant features like differentially expressed genes associated with disease. Currently, most of bioinformatic tools enable normalization in a flexible format that depends on the individual datasets under analysis. Thus, the output data of such normalizations will be poorly compatible with each other. Recently we proposed a new approach to gene expression data normalization termed Shambhala which returns harmonized data in a uniform shape, where every expression profile is transformed into a pre-defined universal format. We previously showed that following shambhalization of human RNA profiles, overall tissue-specific clustering features are strongly retained while platform-specific clustering is dramatically reduced. Methods: Here, we tested Shambhala performance in retention of fold-change gene expression features and other functional characteristics of gene clusters such as pathway activation levels and predicted cancer drug activity scores. Results: Using 6,793 cancer and 11,135 normal tissue gene expression profiles from the literature and experimental datasets, we applied twelve performance criteria for different versions of Shambhala and other methods of transcriptomic harmonization with flexible output data format. Such criteria dealt with the biological type classifiers, hierarchical clustering, correlation/regression properties, stability of drug efficiency scores, and data quality for using machine learning classifiers. Discussion: Shambhala-2 harmonizer demonstrated the best results with the close to 1 correlation and linear regression coefficients for the comparison of training vs validation datasets and more than two times lesser instability for calculation of drug efficiency scores compared to other methods.

Transcriptomic Harmonization as the Way for Suppressing Cross-Platform Bias and Batch Effect

(1) Background: Emergence of methods interrogating gene expression at high throughput gave birth to quantitative transcriptomics, but also posed a question of inter-comparison of expression profiles obtained using different equipment and protocols and/or in different series of experiments. Addressing this issue is challenging, because all of the above variables can dramatically influence gene expression signals and, therefore, cause a plethora of peculiar features in the transcriptomic profiles. Millions of transcriptomic profiles were obtained and deposited in public databases of which the usefulness is however strongly limited due to the inter-comparison issues; (2) Methods: Dozens of methods and software packages that can be generally classified as either flexible or predefined format harmonizers have been proposed, but none has become to the date the gold standard for unification of this type of Big Data; (3) Results: However, recent developments evidence that platform/protocol/batch bias can be efficiently reduced not only for the comparisons of limited transcriptomic datasets. Instead, instruments were proposed for transforming gene expression profiles into the universal, uniformly shaped format that can support multiple inter-comparisons for reasonable calculation costs. This forms a basement for universal indexing of all or most of all types of RNA sequencing and microarray hybridization profiles; (4) Conclusions: In this paper, we attempted to overview the landscape of modern approaches and methods in transcriptomic harmonization and focused on the practical aspects of their application.

Shambhala-2: A Protocol for Uniformly Shaped Harmonization of Gene Expression Profiles of Various Formats

Uniformly shaped harmonization of gene expression profiles is central for the simultaneous comparison of multiple gene expression datasets. It is expected to operate with the gene expression data obtained using various experimental methods and equipment, and to return harmonized profiles in a uniform shape. Such uniformly shaped expression profiles from different initial datasets can be further compared directly. However, current harmonization techniques have strong limitations that prevent their broad use for bioinformatic applications. They can either operate with only up to two datasets/platforms or return data in a dynamic format that will be different for every comparison under analysis. This also does not allow for adding new data to the previously harmonized dataset(s), which complicates the analysis and increases calculation costs. We propose here a new method termed Shambhala-2 that can transform multi-platform expression data into a universal format that is identical for all harmonizations made using this technique. Shambhala-2 is based on sample-by-sample cubic conversion of the initial expression dataset into a preselected shape of the reference definitive dataset. Using 8390 samples of 12 healthy human tissue types and 4086 samples of colorectal, kidney, and lung cancer tissues, we verified Shambhala-2's capacity in restoring tissue-specific expression patterns for seven microarray and three RNA sequencing platforms. Shambhala-2 performed well for all tested combinations of RNAseq and microarray profiles, and retained gene-expression ranks, as evidenced by high correlations between different single- or aggregated gene expression metrics in pre- and post-Shambhalized samples, including preserving cancer-specific gene expression and pathway activation features. © 2022 Wiley Periodicals LLC. Basic Protocol: Shambhala-2 harmonizer Alternate Protocol 1: Linear Shambhala/Shambhala-1 Alternate Protocol 2: Alternative (flexible-format and uniformly shaped) normalization methods Support Protocol 1: Watermelon multisection (WM) Support Protocol 2: Calculation of cancer-to-normal log-fold-change (LFC) and pathway activation level (PAL).

Transcriptomic Portraits and Molecular Pathway Activation Features of Adult Spinal Intramedullary Astrocytomas

In this study, we report 31 spinal intramedullary astrocytoma (SIA) RNA sequencing (RNA-seq) profiles for 25 adult patients with documented clinical annotations. To our knowledge, this is the first clinically annotated RNA-seq dataset of spinal astrocytomas derived from the intradural intramedullary compartment. We compared these tumor profiles with the previous healthy central nervous system (CNS) RNA-seq data for spinal cord and brain and identified SIA-specific gene sets and molecular pathways. Our findings suggest a trend for SIA-upregulated pathways governing interactions with the immune cells and downregulated pathways for the neuronal functioning in the context of normal CNS activity. In two patient tumor biosamples, we identified diagnostic KIAA1549-BRAF fusion oncogenes, and we also found 16 new SIA-associated fusion transcripts. In addition, we bioinformatically simulated activities of targeted cancer drugs in SIA samples and predicted that several tyrosine kinase inhibitory drugs and thalidomide analogs could be potentially effective as second-line treatment agents to aid in the prevention of SIA recurrence and progression.