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/.

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.

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. 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. The BloodGen3Module package and documentation are freely available from Github: https://github.com/Drinchai/BloodGen3Module. Supplementary data are available at Bioinformatics online.

A potential role for the long-chain acyl-CoA synthetase family member 1 (ACSL1) in the immunobiology of sepsis was explored during a hands-on training workshop. Participants first assessed the robustness of the potential gap in biomedical knowledge identified via an initial screen of public transcriptome data and of the literature associated with ACSL1. Increase in ACSL1 transcript abundance during sepsis was confirmed in several independent datasets. Querying the ACSL1 literature also confirmed the absence of reports associating ACSL1 with sepsis. Inferences drawn from both the literature (via indirect associations) and public transcriptome data (via correlation) point to the likely participation of ACSL1 and ACSL4, another family member, in inflammasome activation in neutrophils during sepsis. Furthermore, available clinical data indicate that levels of ACSL1 and ACSL4 induction was significantly higher in fatal cases of sepsis. This denotes potential translational relevance and is consistent with involvement in pathways driving potentially deleterious systemic inflammation. Finally, while ACSL1 expression was induced in blood in vitro by a wide range of pathogen-derived factors as well as TNF, induction of ACSL4 appeared restricted to flagellated bacteria and pathogen-derived TLR5 agonists and IFNG. Taken together, this joint review of public literature and omics data records points to two members of the acyl-CoA synthetase family potentially playing a role in inflammasome activation in neutrophils. Translation relevance of these observations in the context of sepsis and other inflammatory conditions remain to be investigated.

Summary: According to publicly available transcriptome datasets, the abundance of Annexin A3 (ANXA3) is robustly increased during the course of sepsis; however, no studies have examined the biological significance or clinical relevance of ANXA3 in this pathology. Here we explored this interpretation gap and identified possible directions for future research. Based on reference transcriptome datasets, we found that ANXA3 expression is restricted to neutrophils, is upregulatedin vitroafter exposure to plasma obtained from septic patients, and is associated with adverse clinical outcomes. Secondly, an increase in ANXA3 transcript abundance was also observedin vivo, in the blood of septic patients in multiple independent studies. ANXA3 is known to mediate calcium‐dependent granules–phagosome fusion in support of microbicidal activity in neutrophils. More recent work has also shown that ANXA3 enhances proliferation and survival of tumour cells via a Caspase‐3‐dependent mechanism. And this same molecule is also known to play a critical role in regulation of apoptotic events in neutrophils. Thus, we posit that during sepsis ANXA3 might either play a beneficial role, by facilitating microbial clearance and resolution of the infection; or a detrimental role, by prolonging neutrophil survival, which is known to contribute to sepsis‐mediated organ damage.

Biomarkers to assess the risk of developing severe respiratory syncytial virus (RSV) infection are needed. We conducted a meta‐analysis of 490 unique profiles from six public RSV blood transcriptome datasets. A repertoire of 382 well‐characterized transcriptional modules was used to define dominant host responses to RSV infection. The consolidated RSV cohort was stratified according to four traits: “interferon response” (IFN), “neutrophil‐driven inflammation” (Infl), “cell cycle” (CC), and “erythrocytes” (Ery). We identified eight prevalent blood transcriptome phenotypes, of which three Ery+ phenotypes comprised higher proportions of patients requiring intensive care. This finding confirms on a larger scale data from one of our earlier reports describing an association between an erythrocyte signature and RSV disease severity. Further contextual interpretation made it possible to associate this signature with immunosuppressive states (late stage cancer, pharmacological immunosuppression), and with a population of fetal glycophorin A+ erythroid precursors. Furthermore, we posit that this erythrocyte cell signature may be linked to a population of immunosuppressive erythroid cells previously described in the literature, and that overabundance of this cell population in RSV patients may underlie progression to severe disease. These findings outline potential priority areas for biomarker development and investigations into the immune biology of RSV infection. The approach that we developed and employed here should also permit to delineate prevalent blood transcriptome phenotypes in other settings.