A modular framework for the development of targeted Covid-19 blood transcript profiling panels

An interactive web application for exploring systemic lupus erythematosus blood transcriptomic diversity

In the field of complex autoimmune diseases such as systemic lupus erythematosus (SLE), systems immunology approaches have proven invaluable in translational research settings. Large-scale datasets of transcriptome profiling have been collected and made available to the research community in public repositories, but remain poorly accessible and usable by mainstream researchers. Enabling tools and technologies facilitating investigators' interaction with large-scale datasets such as user-friendly web applications could promote data reuse and foster knowledge discovery. Microarray blood transcriptomic data from the LUPUCE cohort (publicly available on Gene Expression Omnibus, GSE49454), which comprised 157 samples from 62 adult SLE patients, were analyzed with the third-generation (BloodGen3) module repertoire framework, which comprises modules and module aggregates. These well-characterized samples corresponded to different levels of disease activity, different types of flares (including biopsy-proven lupus nephritis), different auto-antibody profiles and different levels of interferon signatures. A web application was deployed to present the aggregate-level, module-level and gene-level analysis results from LUPUCE dataset. Users can explore the similarities and heterogeneity of SLE samples, navigate through different levels of analysis, test hypotheses and generate custom fingerprint grids and heatmaps, which may be used in reports or manuscripts. This resource is available via this link: https://immunology-research.shinyapps.io/LUPUCE/. This web application can be employed as a stand-alone resource to explore changes in blood transcript profiles in SLE, and their relation to clinical and immunological parameters, to generate new research hypotheses.

A data browsing application for accessing gene and module-level blood transcriptome profiles of healthy pregnant women from high- and low-resource settings

Transcriptome profiling data, generated via RNA sequencing, are commonly deposited in public repositories. However, these data may not be easily accessible or usable by many researchers. To enhance data reuse, we present well-annotated, partially analyzed data via a user-friendly web application. This project involved transcriptome profiling of blood samples from 15 healthy pregnant women in a low-resource setting, taken at 6 consecutive time points beginning from the first trimester. Additional blood transcriptome profiles were retrieved from the National Center for Biotechnology Information (NCBI) Gene Expression Omnibus (GEO) public repository, representing a cohort of healthy pregnant women from a high-resource setting. We analyzed these datasets using the fixed BloodGen3 module repertoire. We deployed a web application, accessible at https://thejacksonlaboratory.shinyapps.io/BloodGen3_Pregnancy/which displays the module-level analysis results from both original and public pregnancy blood transcriptome datasets. Users can create custom fingerprint grid and heatmap representations via various navigation options, useful for reports and manuscript preparation. The web application serves as a standalone resource for exploring blood transcript abundance changes during pregnancy. Alternatively, users can integrate it with similar applications developed for earlier publications to analyze transcript abundance changes of a given BloodGen3 signature across a range of disease cohorts. Database URL: https://thejacksonlaboratory.shinyapps.io/BloodGen3_Pregnancy/.

Design of a targeted blood transcriptional panel for monitoring immunological changes accompanying pregnancy

Background

Immunomodulatory processes exert steering functions throughout pregnancy. Detecting diversions from this physiologic immune clock may help identify pregnant women at risk for pregnancy-associated complications. We present results from a data-driven selection process to develop a targeted panel of mRNAs that may prove effective in detecting pregnancies diverting from the norm.

Methods

Based on a de novo dataset from a resource-constrained setting and a dataset from a resource-rich area readily available in the public domain, whole blood gene expression profiles of uneventful pregnancies were captured at multiple time points during pregnancy. BloodGen3, a fixed blood transcriptional module repertoire, was employed to analyze and visualize gene expression patterns in the two datasets. Differentially expressed genes were identified by comparing their abundance to non-pregnant postpartum controls. The selection process for a targeted gene panel considered (i) transcript abundance in whole blood; (ii) degree of correlation with the BloodGen3 module; and (iii) pregnancy biology.

Results

We identified 176 transcripts that were complemented with eight housekeeping genes. Changes in transcript abundance were seen in the early stages of pregnancy and similar patterns were observed in both datasets. Functional gene annotation suggested significant changes in the lymphoid, prostaglandin and inflammation-associated compartments, when compared to the postpartum controls.

Conclusion

The gene panel presented here holds promise for the development of predictive, targeted, transcriptional profiling assays. Such assays might become useful for monitoring of pregnant women, specifically to detect potential adverse events early. Prospective validation of this targeted assay, in-depth investigation of functional annotations of differentially expressed genes, and assessment of common pregnancy-associated complications with the aim to identify these early in pregnancy to improve pregnancy outcomes are the next steps.

Harnessing large language models (LLMs) for candidate gene prioritization and selection

BACKGROUND

Feature selection is a critical step for translating advances afforded by systems-scale molecular profiling into actionable clinical insights. While data-driven methods are commonly utilized for selecting candidate genes, knowledge-driven methods must contend with the challenge of efficiently sifting through extensive volumes of biomedical information. This work aimed to assess the utility of large language models (LLMs) for knowledge-driven gene prioritization and selection.

METHODS

In this proof of concept, we focused on 11 blood transcriptional modules associated with an Erythroid cells signature. We evaluated four leading LLMs across multiple tasks. Next, we established a workflow leveraging LLMs. The steps consisted of: (1) Selecting one of the 11 modules; (2) Identifying functional convergences among constituent genes using the LLMs; (3) Scoring candidate genes across six criteria capturing the gene's biological and clinical relevance; (4) Prioritizing candidate genes and summarizing justifications; (5) Fact-checking justifications and identifying supporting references; (6) Selecting a top candidate gene based on validated scoring justifications; and (7) Factoring in transcriptome profiling data to finalize the selection of the top candidate gene.

RESULTS

Of the four LLMs evaluated, OpenAI's GPT-4 and Anthropic's Claude demonstrated the best performance and were chosen for the implementation of the candidate gene prioritization and selection workflow. This workflow was run in parallel for each of the 11 erythroid cell modules by participants in a data mining workshop. Module M9.2 served as an illustrative use case. The 30 candidate genes forming this module were assessed, and the top five scoring genes were identified as BCL2L1, ALAS2, SLC4A1, CA1, and FECH. Researchers carefully fact-checked the summarized scoring justifications, after which the LLMs were prompted to select a top candidate based on this information. GPT-4 initially chose BCL2L1, while Claude selected ALAS2. When transcriptional profiling data from three reference datasets were provided for additional context, GPT-4 revised its initial choice to ALAS2, whereas Claude reaffirmed its original selection for this module.

CONCLUSIONS

Taken together, our findings highlight the ability of LLMs to prioritize candidate genes with minimal human intervention. This suggests the potential of this technology to boost productivity, especially for tasks that require leveraging extensive biomedical knowledge.

Multi-objective optimization identifies a specific and interpretable COVID-19 host response signature

The identification of a COVID-19 host response signature in blood can increase the understanding of SARS-CoV-2 pathogenesis and improve diagnostic tools. Applying a multi-objective optimization framework to both massive public and new multi-omics data, we identified a COVID-19 signature regulated at both transcriptional and epigenetic levels. We validated the signature's robustness in multiple independent COVID-19 cohorts. Using public data from 8,630 subjects and 53 conditions, we demonstrated no cross-reactivity with other viral and bacterial infections, COVID-19 comorbidities, or confounders. In contrast, previously reported COVID-19 signatures were associated with significant cross-reactivity. The signature's interpretation, based on cell-type deconvolution and single-cell data analysis, revealed prominent yet complementary roles for plasmablasts and memory T cells. Although the signal from plasmablasts mediated COVID-19 detection, the signal from memory T cells controlled against cross-reactivity with other viral infections. This framework identified a robust, interpretable COVID-19 signature and is broadly applicable in other disease contexts. A record of this paper's transparent peer review process is included in the supplemental information.

An artificial neural network model based on autophagy-related genes in childhood systemic lupus erythematosus

Background

Childhood systemic lupus erythematosus (cSLE) is a multisystemic, life-threatening autoimmune disease. Compared to adults, SLE in childhood is more active, can cause multisystem involvement including renal, neurological and hematological, and can cause cumulative damage across systems more rapidly. Autophagy, one of the core functions of cells, is involved in almost every process of the immune response and has been shown to be associated with many autoimmune diseases, being a key factor in the interplay between innate and adaptive immunity. Autophagy influences the onset, progression and severity of SLE. This paper identifies new biomarkers for the diagnosis and treatment of childhood SLE based on an artificial neural network of autophagy-related genes.

Methods

We downloaded dataset GSE100163 from the Gene Expression Omnibus database and used Protein–protein Interaction Network (PPI) and Least Absolute Shrinkage and Selection Operator (LASSO) to screen the signature genes of autophagy-related genes in cSLE. A new artificial neural network model for cSLE diagnosis was constructed using the signature genes. The predictive efficiency of the model was also validated using the dataset GSE65391. Finally, "CIBERSORT" was used to calculate the infiltration of immune cells in cSLE and to analyze the relationship between the signature genes and the infiltration of immune cells.

Results

We identified 37 autophagy-related genes that differed in cSLE and normal samples, and finally obtained the seven most relevant signature genes for cSLE (DDIT3, GNB2L1, CTSD, HSPA8, ULK1, DNAJB1, CANX) by PPI and LASOO regression screening, and constructed an artificial neural network diagnostic model for cSLE. Using this model, we plotted the ROC curves for the training and validation group diagnoses with the area under the curve of 0.976 and 0.783, respectively. Finally, we performed immunoassays on cSLE samples, and the results showed that Plasma cells, Macrophages M0, Dendritic cells activated and Neutrophils were significantly infiltrated in cSLE.

Conclusion

We constructed an artificial neural network diagnostic model of seven autophagy-related genes that can be used for the diagnosis of cSLE. Meanwhile, the characteristic genes affect the immune infiltration of cSLE, which may provide new perspectives for the exploration of cSLE treatment and related mechanisms.

MX2: Identification and systematic mechanistic analysis of a novel immune-related biomarker for systemic lupus erythematosus

Background

Systemic lupus erythematosus (SLE) is an autoimmune disease that involves multiple organs. However, the current SLE-related biomarkers still lack sufficient sensitivity, specificity and predictive power for clinical application. Thus, it is significant to explore new immune-related biomarkers for SLE diagnosis and development.

Methods

We obtained seven SLE gene expression profile microarrays (GSE121239/11907/81622/65391/100163/45291/49454) from the GEO database. First, differentially expressed genes (DEGs) were screened using GEO2R, and SLE biomarkers were screened by performing WGCNA, Random Forest, SVM-REF, correlation with SLEDAI and differential gene analysis. Receiver operating characteristic curves (ROCs) and AUC values were used to determine the clinical value. The expression level of the biomarker was verified by RT‒qPCR. Subsequently, functional enrichment analysis was utilized to identify biomarker-associated pathways. ssGSEA, CIBERSORT, xCell and ImmuCellAI algorithms were applied to calculate the sample immune cell infiltration abundance. Single-cell data were analyzed for gene expression specificity in immune cells. Finally, the transcriptional regulatory network of the biomarker was constructed, and the corresponding therapeutic drugs were predicted.

Results

Multiple algorithms were screened together for a unique marker gene, MX2, and expression analysis of multiple datasets revealed that MX2 was highly expressed in SLE compared to the normal group (all P < 0.05), with the same trend validated by RT‒qPCR (P = 0.026). Functional enrichment analysis identified the main pathway of MX2 promotion in SLE as the NOD-like receptor signaling pathway (NES=2.492, P < 0.001, etc.). Immuno-infiltration analysis showed that MX2 was closely associated with neutrophils, and single-cell and transcriptomic data revealed that MX2 was specifically expressed in neutrophils. The NOD-like receptor signaling pathway was also remarkably correlated with neutrophils (r >0.3, P < 0.001, etc.). Most of the MX2-related interacting proteins were associated with SLE, and potential transcription factors of MX2 and its related genes were also significantly associated with the immune response.

Conclusion

Our study found that MX2 can serve as an immune-related biomarker for predicting the diagnosis and disease activity of SLE. It activates the NOD-like receptor signaling pathway and promotes neutrophil infiltration to aggravate SLE.

Transcriptomic profile investigations highlight a putative role for NUDT16 in sepsis

Sepsis is an aberrant systemic inflammatory response mediated by the acute activation of the innate immune system. Neutrophils are important contributors to the innate immune response that controls the infection, but harbour the risk of collateral tissue damage such as thrombosis and organ dysfunction. A better understanding of the modulations of cellular processes in neutrophils and other blood cells during sepsis is needed and can be initiated via transcriptomic profile investigations. To that point, the growing repertoire of publicly accessible transcriptomic datasets serves as a valuable resource for discovering and/or assessing the robustness of biomarkers. We employed systematic literature mining, reductionist approach to gene expression profile and empirical in vitro work to highlight the role of a Nudix hydrolase family member, NUDT16, in sepsis. The relevance and implication of the expression of NUDT16 under septic conditions and the putative functional roles of this enzyme are discussed.

Transcriptome and Literature Mining Highlight the Differential Expression of ERLIN1 in Immune Cells during Sepsis

Sepsis results from the dysregulation of the host immune system. This highly variable disease affects 19 million people globally, and accounts for 5 million deaths annually. In transcriptomic datasets curated from public repositories, we observed a consistent upregulation (3.26-5.29 fold) of ERLIN1-a gene coding for an ER membrane prohibitin and a regulator of inositol 1, 4, 5-trisphosphate receptors and sterol regulatory element-binding proteins-under septic conditions in healthy neutrophils, monocytes, and whole blood. In vitro expression of the ERLIN1 gene and proteins was measured by stimulating the whole blood of healthy volunteers to a combination of lipopolysaccharide and peptidoglycan. Septic stimulation induced a significant increase in ERLIN1 expression; however, ERLIN1 was differentially expressed among the immune blood cell subsets. ERLIN1 was uniquely increased in whole blood neutrophils, and confirmed in the differentiated HL60 cell line. The scarcity of ERLIN1 in sepsis literature indicates a knowledge gap between the functions of ERLIN1, calcium homeostasis, and cholesterol and fatty acid biosynthesis, and sepsis. In combination with experimental data, we bring forth the hypothesis that ERLIN1 is variably modulated among immune cells in response to cellular perturbations, and has implications for ER functions and/or ER membrane protein components during sepsis.

Development of a fixed module repertoire for the analysis and interpretation of blood transcriptome data

As the capacity for generating large-scale molecular profiling data continues to grow, the ability to extract meaningful biological knowledge from it remains a limitation. Here, we describe the development of a new fixed repertoire of transcriptional modules, BloodGen3, that is designed to serve as a stable reusable framework for the analysis and interpretation of blood transcriptome data. The construction of this repertoire is based on co-clustering patterns observed across sixteen immunological and physiological states encompassing 985 blood transcriptome profiles. Interpretation is supported by customized resources, including module-level analysis workflows, fingerprint grid plot visualizations, interactive web applications and an extensive annotation framework comprising functional profiling reports and reference transcriptional profiles. Taken together, this well-characterized and well-supported transcriptional module repertoire can be employed for the interpretation and benchmarking of blood transcriptome profiles within and across patient cohorts. Blood transcriptome fingerprints for the 16 reference cohorts can be accessed interactively via:  https://drinchai.shinyapps.io/BloodGen3Module/ .

Integrated transcriptional-phenotypic analysis captures systemic immunomodulation following antiangiogenic therapy in renal cell carcinoma patients

BACKGROUND

The combination of immune checkpoint blockade (ICB) with standard therapies is becoming a common approach for overcoming resistance to cancer immunotherapy in most human malignancies including metastatic renal cell carcinoma (mRCC). In this regard, insights into the immunomodulatory properties of antiangiogenic agents may help designing multidrug schedules based on specific immune synergisms.

METHODS

We used orthogonal transcriptomic and phenotyping platforms combined with functional analytic pipelines to elucidate the immunomodulatory effect of the antiangiogenic agent pazopanib in mRCC patients. Nine patients were studied longitudinally over a period of 6 months. We also analyzed transcriptional data from The Cancer Genome Atlas (TCGA) RCC cohort (N = 571) to assess the prognostic implications of our findings. The effect of pazopanib was assessed in vitro on NK cells and T cells. Additionally, myeloid-derived suppressor (MDSC)-like cells were generated from CD14+ monocytes transfected with mimics of miRNAs associated with MDSC function in the presence or absence of pazopanib.

RESULTS

Pazopanib administration caused a rapid and dramatic reshaping in terms of frequency and transcriptional activity of multiple blood immune cell subsets, with a downsizing of MDSC and regulatory T cells in favor of a strong enhancement in PD-1 expressing cytotoxic T and Natural Killer effectors. These changes were paired with an increase of the expression of transcripts reflecting activation of immune-effector functions. This immunomodulation was marked but transient, peaking at the third month of treatment. Moreover, the intratumoral expression level of a MDSC signature (MDSC INT) was strongly associated with poor prognosis in RCC patients. In vitro experiments indicate that the observed immunomodulation might be due to an inhibitory effect on MDSC-mediated suppression, rather than a direct effect on NK and T cells.

CONCLUSIONS

The marked but transient nature of this immunomodulation, peaking at the third month of treatment, provides the rationale for the use of antiangiogenics as a preconditioning strategy to improve the efficacy of ICB.

BloodGen3Module: Blood transcriptional module repertoire analysis and visualization using R

Motivation

We previously described the construction and characterization of fixed reusable blood transcriptional module repertoires. More recently we released a third iteration (‘BloodGen3’ module repertoire) that comprises 382 functionally annotated modules and encompasses 14 168 transcripts. Custom bioinformatic tools are needed to support downstream analysis, visualization and interpretation relying on such fixed module repertoires.

Results

We have developed and describe here an R package, BloodGen3Module. The functions of our package permit group comparison analyses to be performed at the module-level, and to display the results as annotated fingerprint grid plots. A parallel workflow for computing module repertoire changes for individual samples rather than groups of samples is also available; these results are displayed as fingerprint heatmaps. An illustrative case is used to demonstrate the steps involved in generating blood transcriptome repertoire fingerprints of septic patients. Taken together, this resource could facilitate the analysis and interpretation of changes in blood transcript abundance observed across a wide range of pathological and physiological states.

Availability and implementation

The BloodGen3Module package and documentation are freely available from Github: https://github.com/Drinchai/BloodGen3Module.

Supplementary information

Supplementary data are available at Bioinformatics online.