Deconvolution of blood microarray data identifies cellular activation patterns in systemic lupus erythematosus

Fourteen years of cellular deconvolution: methodology, applications, technical evaluation and outstanding challenges

Single-cell RNA sequencing (scRNA-Seq) is a recent technology that allows for the measurement of the expression of all genes in each individual cell contained in a sample. Information at the single-cell level has been shown to be extremely useful in many areas. However, performing single-cell experiments is expensive. Although cellular deconvolution cannot provide the same comprehensive information as single-cell experiments, it can extract cell-type information from bulk RNA data, and therefore it allows researchers to conduct studies at cell-type resolution from existing bulk datasets. For these reasons, a great effort has been made to develop such methods for cellular deconvolution. The large number of methods available, the requirement of coding skills, inadequate documentation, and lack of performance assessment all make it extremely difficult for life scientists to choose a suitable method for their experiment. This paper aims to fill this gap by providing a comprehensive review of 53 deconvolution methods regarding their methodology, applications, performance, and outstanding challenges. More importantly, the article presents a benchmarking of all these 53 methods using 283 cell types from 30 tissues of 63 individuals. We also provide an R package named DeconBenchmark that allows readers to execute and benchmark the reviewed methods (https://github.com/tinnlab/DeconBenchmark).

Regional and bilateral MRI and gene signatures in facioscapulohumeral dystrophy: implications for clinical trial design and mechanisms of disease progression

Identifying the aberrant expression of DUX4 in skeletal muscle as the cause of facioscapulohumeral dystrophy (FSHD) has led to rational therapeutic development and clinical trials. Several studies support the use of MRI characteristics and the expression of DUX4-regulated genes in muscle biopsies as biomarkers of FSHD disease activity and progression. We performed lower-extremity MRI and muscle biopsies in the mid-portion of the tibialis anterior (TA) muscles bilaterally in FSHD subjects and validated our prior reports of the strong association between MRI characteristics and expression of genes regulated by DUX4 and other gene categories associated with FSHD disease activity. We further show that measurements of normalized fat content in the entire TA muscle strongly predict molecular signatures in the mid-portion of the TA, indicating that regional biopsies can accurately measure progression in the whole muscle and providing a strong basis for inclusion of MRI and molecular biomarkers in clinical trial design. An unanticipated finding was the strong correlations of molecular signatures in the bilateral comparisons, including markers of B-cells and other immune cell populations, suggesting that a systemic immune cell infiltration of skeletal muscle might have a role in disease progression.

An Emerging Approach of Age-Related Hearing Loss Research: Application of Integrated Multi-Omics Analysis

As one of the most common otologic diseases in the elderly, age-related hearing loss (ARHL) usually characterized by hearing loss and cognitive disorders, which have a significant impact on the elderly's physical and mental health and quality of life. However, as a typical disease of aging, it is unclear why aging causes widespread hearing impairment in the elderly. As molecular biological experiments have been conducted for research recently, ARHL is gradually established at various levels with the application and development of integrated multi-omics analysis in the studies of ARHL. Here, the recent progress in the application of multi-omics analysis in the molecular mechanisms of ARHL development and therapeutic regimens, including the combined analysis of different omics, such as transcriptome, proteome, and metabolome, to screen for risk sites, risk genes, and differences in lipid metabolism, etc., is outlined and the integrated histological data further promote the profound understanding of the disease process as well as physiological mechanisms of ARHL. The advantages and disadvantages of multi-omics analysis in disease research are also discussed and the authors speculate on the future prospects and applications of this part-to-whole approach, which may provide more comprehensive guidance for ARHL and aging disease prevention and treatment.

CATD: a reproducible pipeline for selecting cell-type deconvolution methods across tissues

Motivation

Cell-type deconvolution methods aim to infer cell composition from bulk transcriptomic data. The proliferation of developed methods coupled with inconsistent results obtained in many cases, highlights the pressing need for guidance in the selection of appropriate methods. Additionally, the growing accessibility of single-cell RNA sequencing datasets, often accompanied by bulk expression from related samples enable the benchmark of existing methods.

Results

In this study, we conduct a comprehensive assessment of 31 methods, utilizing single-cell RNA-sequencing data from diverse human and mouse tissues. Employing various simulation scenarios, we reveal the efficacy of regression-based deconvolution methods, highlighting their sensitivity to reference choices. We investigate the impact of bulk-reference differences, incorporating variables such as sample, study and technology. We provide validation using a gold standard dataset from mononuclear cells and suggest a consensus prediction of proportions when ground truth is not available. We validated the consensus method on data from the stomach and studied its spillover effect. Importantly, we propose the use of the critical assessment of transcriptomic deconvolution (CATD) pipeline which encompasses functionalities for generating references and pseudo-bulks and running implemented deconvolution methods. CATD streamlines simultaneous deconvolution of numerous bulk samples, providing a practical solution for speeding up the evaluation of newly developed methods.

Availability and implementation

https://github.com/Papatheodorou-Group/CATD_snakemake.

Advances in cancer DNA methylation analysis with methPLIER: use of non-negative matrix factorization and knowledge-based constraints to enhance biological interpretability

DNA methylation is an epigenetic modification that results in dynamic changes during ontogenesis and cell differentiation. DNA methylation patterns regulate gene expression and have been widely researched. While tools for DNA methylation analysis have been developed, most of them have focused on intergroup comparative analysis within a dataset; therefore, it is difficult to conduct cross-dataset studies, such as rare disease studies or cross-institutional studies. This study describes a novel method for DNA methylation analysis, namely, methPLIER, which enables interdataset comparative analyses. methPLIER combines Pathway Level Information Extractor (PLIER), which is a non-negative matrix factorization (NMF) method, with regularization by a knowledge matrix and transfer learning. methPLIER can be used to perform intersample and interdataset comparative analysis based on latent feature matrices, which are obtained via matrix factorization of large-scale data, and factor-loading matrices, which are obtained through matrix factorization of the data to be analyzed. We used methPLIER to analyze a lung cancer dataset and confirmed that the data decomposition reflected sample characteristics for recurrence-free survival. Moreover, methPLIER can analyze data obtained via different preprocessing methods, thereby reducing distributional bias among datasets due to preprocessing. Furthermore, methPLIER can be employed for comparative analyses of methylation data obtained from different platforms, thereby reducing bias in data distribution due to platform differences. methPLIER is expected to facilitate cross-sectional DNA methylation data analysis and enhance DNA methylation data resources.

Modeling cancer progression: an integrated workflow extending data-driven kinetic models to bio-mechanical PDE models

Computational modeling of cancer can help unveil dynamics and interactions that are hard to replicate experimentally. Thanks to the advancement in cancer databases and data analysis technologies, these models have become more robust than ever. There are many mathematical models which investigate cancer through different approaches, from sub-cellular to tissue scale, and from treatment to diagnostic points of view. In this study, we lay out a step-by-step methodology for a data-driven mechanistic model of the tumor microenvironment. We discuss data acquisition strategies, data preparation, parameter estimation, and sensitivity analysis techniques. Furthermore, we propose a possible approach to extend mechanistic ordinary differential equation models to PDE models coupled with mechanical growth. The workflow discussed in this article can help understand the complex temporal and spatial interactions between cells and cytokines in the tumor microenvironment and their effect on tumor growth.

Knockdown of Hyaluronan synthase 2 suppresses liver fibrosis in mice via induction of transcriptomic changes similar to 4MU treatment

Hepatic fibrosis remains a significant clinical challenge due to ineffective treatments. 4-methylumbelliferone (4MU), a hyaluronic acid (HA) synthesis inhibitor, has proven safe in phase one clinical trials. In this study, we aimed to ameliorate liver fibrosis by inhibiting HA synthesis. We compared two groups of mice with CCl4-induced fibrosis, treated with 4-methylumbelliferone (4MU) and hyaluronan synthase 2 (HAS2) targeting siRNA (siHAS2). The administration of 4MU and siHAS2 significantly reduced collagen and HA deposition, as well as biochemical markers of hepatic damage induced by repeated CCl4 injections. The transcriptomic analysis revealed converging pathways associated with downstream HA signalling. 4MU- and siHAS2-treated fibrotic livers shared 405 upregulated and 628 downregulated genes. These genes were associated with xenobiotic and cholesterol metabolism, mitosis, endoplasmic reticulum stress, RNA processing, and myeloid cell migration. The functional annotation of differentially expressed genes (DEGs) in siHAS2-treated mice revealed attenuation of extracellular matrix-associated pathways. In comparison, in the 4MU-treated group, DEGs were related to lipid and bile metabolism pathways and cell cycle. These findings confirm that HAS2 is an important pharmacological target for suppressing hepatic fibrosis using siRNA.

Highly Accurate Estimation of Cell Type Abundance in Bulk Tissues Based on Single-Cell Reference and Domain Adaptive Matching

Accurately identifies the cellular composition of complex tissues, which is critical for understanding disease pathogenesis, early diagnosis, and prevention. However, current methods for deconvoluting bulk RNA sequencing (RNA-seq) typically rely on matched single-cell RNA sequencing (scRNA-seq) as a reference, which can be limiting due to differences in sequencing distribution and the potential for invalid information from single-cell references. Hence, a novel computational method named SCROAM is introduced to address these challenges. SCROAM transforms scRNA-seq and bulk RNA-seq into a shared feature space, effectively eliminating distributional differences in the latent space. Subsequently, cell-type-specific expression matrices are generated from the scRNA-seq data, facilitating the precise identification of cell types within bulk tissues. The performance of SCROAM is assessed through benchmarking against simulated and real datasets, demonstrating its accuracy and robustness. To further validate SCROAM's performance, single-cell and bulk RNA-seq experiments are conducted on mouse spinal cord tissue, with SCROAM applied to identify cell types in bulk tissue. Results indicate that SCROAM is a highly effective tool for identifying similar cell types. An integrated analysis of liver cancer and primary glioblastoma is then performed. Overall, this research offers a novel perspective for delivering precise insights into disease pathogenesis and potential therapeutic strategies.

A personalized network framework reveals predictive axis of anti-TNF response across diseases

Personalized treatment of complex diseases has been mostly predicated on biomarker identification of one drug-disease combination at a time. Here, we use a computational approach termed Disruption Networks to generate a data type, contextualized by cell-centered individual-level networks, that captures biology otherwise overlooked when performing standard statistics. This data type extends beyond the "feature level space", to the "relations space", by quantifying individual-level breaking or rewiring of cross-feature relations. Applying Disruption Networks to dissect high-dimensional blood data, we discover and validate that the RAC1-PAK1 axis is predictive of anti-TNF response in inflammatory bowel disease. Intermediate monocytes, which correlate with the inflammatory state, play a key role in the RAC1-PAK1 responses, supporting their modulation as a therapeutic target. This axis also predicts response in rheumatoid arthritis, validated in three public cohorts. Our findings support blood-based drug response diagnostics across immune-mediated diseases, implicating common mechanisms of non-response.

Challenges and opportunities to computationally deconvolve heterogeneous tissue with varying cell sizes using single-cell RNA-sequencing datasets

Deconvolution of cell mixtures in "bulk" transcriptomic samples from homogenate human tissue is important for understanding disease pathologies. However, several experimental and computational challenges impede transcriptomics-based deconvolution approaches using single-cell/nucleus RNA-seq reference atlases. Cells from the brain and blood have substantially different sizes, total mRNA, and transcriptional activities, and existing approaches may quantify total mRNA instead of cell type proportions. Further, standards are lacking for the use of cell reference atlases and integrative analyses of single-cell and spatial transcriptomics data. We discuss how to approach these key challenges with orthogonal "gold standard" datasets for evaluating deconvolution methods.

cellsig plug-in enhances CIBERSORTx signature selection for multidataset transcriptomes with sparse multilevel modelling

MOTIVATION

The precise characterization of cell-type transcriptomes is pivotal to understanding cellular lineages, deconvolution of bulk transcriptomes, and clinical applications. Single-cell RNA sequencing resources like the Human Cell Atlas have revolutionised cell-type profiling. However, challenges persist due to data heterogeneity and discrepancies across different studies. One limitation of prevailing tools such as CIBERSORTx is their inability to address hierarchical data structures and handle nonoverlapping gene sets across samples, relying on filtering or imputation.

RESULTS

Here, we present cellsig, a Bayesian sparse multilevel model designed to improve signature estimation by adjusting data for multilevel effects and modelling for gene-set sparsity. Our model is tailored to large-scale, heterogeneous pseudobulk and bulk RNA sequencing data collections with nonoverlapping gene sets. We tested the performances of cellsig on a novel curated Human Bulk Cell-type Catalogue, which harmonizes 1435 samples across 58 datasets. We show that cellsig significantly enhances cell-type marker gene ranking performance. This approach is valuable for cell-type signature selection, with implications for marker gene validation, single-cell annotation, and deconvolution benchmarks.

AVAILABILITY AND IMPLEMENTATION

Codes and the interactive app are available at https://github.com/stemangiola/cellsig; and the database is available at https://doi.org/10.5281/zenodo.7582421.

A modified immune cell infiltration score achieves ideal stratification for CD8+ T cell efficacy and immunotherapy benefit in hepatocellular carcinoma

Immunotherapy, which aims to enhance the function of T cells, has emerged as a novel therapeutic approach for hepatocellular carcinoma (HCC). Nevertheless, the clinical utility of using flow cytometry to assess immune cell infiltration (ICI) is hindered by its cumbersome procedures, prompting the need for more accessible methods. Here, we acquired gene expression profiles and survival data of HCC from TCGA and GSE10186 datasets. The patients were categorized into two clusters of ICI, and a set of 11 characteristic genes responsible for the differentiation performance of these ICI clusters were identified. Subsequently, we successfully developed a modified ICI score (mICIS) by utilizing the expression levels of these genes. The efficacy of our mICIS was confirmed via mass cytometry, flow cytometry, and immunohistochemistry. Our research indicated that the favorable overall survival (OS) rate could be attributed to the improved function of anti-tumor leukocytes rather than their infiltration. Furthermore, we observed that the low score group exhibited lower expression levels of T-cell exhaustion-associated genes, which was confirmed in both HCC tissues from patients and mice, which demonstrated that the benefits of the low scores were due to enhanced active/cytotoxic CD8+ T cells and reduced exhausted CD8+ T cells. Additionally, our mICIS stratified the benefits derived from immunotherapies. Lastly, we observed a misalignment between CD8+ T-cell infiltration and function in HCC. In summary, our mICIS demonstrated proficiency in assessing the OS rate of HCC and offering significant stratified data pertaining to distinct responses to immunotherapy.

Rat Deconvolution as Knowledge Miner for Immune Cell Trafficking from Toxicogenomics Databases

Toxicogenomics databases are useful for understanding biological responses in individuals because they include a diverse spectrum of biological responses. Although these databases contain no information regarding immune cells in the liver, which are important in the progression of liver injury, deconvolution that estimates cell-type proportions from bulk transcriptome could extend immune information. However, deconvolution has been mainly applied to humans and mice and less often to rats, which are the main target of toxicogenomics databases. Here, we developed a deconvolution method for rats to retrieve information regarding immune cells from toxicogenomics databases. The rat-specific deconvolution showed high correlations for several types of immune cells between spleen and blood, and between liver treated with toxicants compared with those based on human and mouse data. Additionally, we found 4 clusters of compounds in Open TG-GATEs database based on estimated immune cell trafficking, which are different from those based on transcriptome data itself. The contributions of this work are three-fold. First, we obtained the gene expression profiles of 6 rat immune cells necessary for deconvolution. Second, we clarified the importance of species differences on deconvolution. Third, we retrieved immune cell trafficking from toxicogenomics databases. Accumulated and comparable immune cell profiles of massive data of immune cell trafficking in rats could deepen our understanding of enable us to clarify the relationship between the order and the contribution rate of immune cells, chemokines and cytokines, and pathologies. Ultimately, these findings will lead to the evaluation of organ responses in Adverse Outcome Pathway.

Multi-omic identification of key transcriptional regulatory programs during endurance exercise training

Transcription factors (TFs) play a key role in regulating gene expression and responses to stimuli. We conducted an integrated analysis of chromatin accessibility, DNA methylation, and RNA expression across eight rat tissues following endurance exercise training (EET) to map epigenomic changes to transcriptional changes and determine key TFs involved. We uncovered tissue-specific changes and TF motif enrichment across all omic layers, differentially accessible regions (DARs), differentially methylated regions (DMRs), and differentially expressed genes (DEGs). We discovered distinct routes of EET-induced regulation through either epigenomic alterations providing better access for TFs to affect target genes, or via changes in TF expression or activity enabling target gene response. We identified TF motifs enriched among correlated epigenomic and transcriptomic alterations, DEGs correlated with exercise-related phenotypic changes, and EET-induced activity changes of TFs enriched for DEGs among their gene targets. This analysis elucidates the unique transcriptional regulatory mechanisms mediating diverse organ effects of EET.

The T Cell Immunoscore as a Reference for Biomarker Development Utilizing Real-World Data from Patients with Advanced Malignancies Treated with Immune Checkpoint Inhibitors

BACKGROUND

We aimed to determine the prognostic value of an immunoscore reflecting CD3+ and CD8+ T cell density estimated from real-world transcriptomic data of a patient cohort with advanced malignancies treated with immune checkpoint inhibitors (ICIs) in an effort to validate a reference for future machine learning-based biomarker development.

METHODS

Transcriptomic data was collected under the Total Cancer Care Protocol (NCT03977402) Avatar® project. The real-world immunoscore for each patient was calculated based on the estimated densities of tumor CD3+ and CD8+ T cells utilizing CIBERSORTx and the LM22 gene signature matrix. Then, the immunoscore association with overall survival (OS) was estimated using Cox regression and analyzed using Kaplan-Meier curves. The OS predictions were assessed using Harrell's concordance index (C-index). The Youden index was used to identify the optimal cut-off point. Statistical significance was assessed using the log-rank test.

RESULTS

Our study encompassed 522 patients with four cancer types. The median duration to death was 10.5 months for the 275 participants who encountered an event. For the entire cohort, the results demonstrated that transcriptomics-based immunoscore could significantly predict patients at risk of death (p-value < 0.001). Notably, patients with an intermediate-high immunoscore achieved better OS than those with a low immunoscore. In subgroup analysis, the prediction of OS was significant for melanoma and head and neck cancer patients but did not reach significance in the non-small cell lung cancer or renal cell carcinoma cohorts.

CONCLUSIONS

Calculating CD3+ and CD8+ T cell immunoscore using real-world transcriptomic data represents a promising signature for estimating OS with ICIs and can be used as a reference for future machine learning-based biomarker development.

Mapping the castor bean endosperm proteome revealed a metabolic interaction between plastid, mitochondria, and peroxisomes to optimize seedling growth

In this work, we studied castor-oil plant Ricinus communis as a classical system for endosperm reserve breakdown. The seeds of castor beans consist of a centrally located embryo with the two thin cotyledons surrounded by the endosperm. The endosperm functions as major storage tissue and is packed with nutritional reserves, such as oil, proteins, and starch. Upon germination, mobilization of the storage reserves requires inter-organellar interplay of plastids, mitochondria, and peroxisomes to optimize growth for the developing seedling. To understand their metabolic interactions, we performed a large-scale organellar proteomic study on castor bean endosperm. Organelles from endosperm of etiolated seedlings were isolated and subjected to liquid chromatography-tandem mass spectrometry (LC-MS/MS). Computer-assisted deconvolution algorithms were applied to reliably assign the identified proteins to their correct subcellular localization and to determine the abundance of the different organelles in the heterogeneous protein samples. The data obtained were used to build a comprehensive metabolic model for plastids, mitochondria, and peroxisomes during storage reserve mobilization in castor bean endosperm.

Identifying lupus Patient Subsets Through Immune Cell Deconvolution of Gene Expression Data in Two Atacicept Phase II Studies

OBJECTIVE

To use cell-based gene signatures to identify patients with systemic lupus erythematous (SLE) in the phase II/III APRIL-SLE and phase IIb ADDRESS II trials most likely to respond to atacicept.

METHODS

A published immune cell deconvolution algorithm based on Affymetrix gene array data was applied to whole blood gene expression from patients entering APRIL-SLE. Five distinct patient clusters were identified. Patient characteristics, biomarkers, and clinical response to atacicept were assessed per cluster. A modified immune cell deconvolution algorithm was developed based on RNA sequencing data and applied to ADDRESS II data to identify similar patient clusters and their responses.

RESULTS

Patients in APRIL-SLE (N = 105) were segregated into the following five clusters (P1-5) characterized by dominant cell subset signatures: high neutrophils, T helper cells and natural killer (NK) cells (P1), high plasma cells and activated NK cells (P2), high B cells and neutrophils (P3), high B cells and low neutrophils (P4), or high activated dendritic cells, activated NK cells, and neutrophils (P5). Placebo- and atacicept-treated patients in clusters P2,4,5 had markedly higher British Isles Lupus Assessment Group (BILAG) A/B flare rates than those in clusters P1,3, with a greater treatment effect of atacicept on lowering flares in clusters P2,4,5. In ADDRESS II, placebo-treated patients from P2,4,5 were less likely to be SLE Responder Index (SRI)-4, SRI-6, and BILAG-Based Combined Lupus Assessment responders than those in P1,3; the response proportions again suggested lower placebo effect and a greater treatment differential for atacicept in P2,4,5.

CONCLUSION

This exploratory analysis indicates larger differences between placebo- and atacicept-treated patients with SLE in a molecularly defined patient subset.

Performance of tumour microenvironment deconvolution methods in breast cancer using single-cell simulated bulk mixtures

Cells within the tumour microenvironment (TME) can impact tumour development and influence treatment response. Computational approaches have been developed to deconvolve the TME from bulk RNA-seq. Using scRNA-seq profiling from breast tumours we simulate thousands of bulk mixtures, representing tumour purities and cell lineages, to compare the performance of nine TME deconvolution methods (BayesPrism, Scaden, CIBERSORTx, MuSiC, DWLS, hspe, CPM, Bisque, and EPIC). Some methods are more robust in deconvolving mixtures with high tumour purity levels. Most methods tend to mis-predict normal epithelial for cancer epithelial as tumour purity increases, a finding that is validated in two independent datasets. The breast cancer molecular subtype influences this mis-prediction. BayesPrism and DWLS have the lowest combined numbers of false positives and false negatives, and have the best performance when deconvolving granular immune lineages. Our findings highlight the need for more single-cell characterisation of rarer cell types, and suggest that tumour cell compositions should be considered when deconvolving the TME.

TOX2 coordinates with TET2 to positively regulate central memory differentiation in human CAR T cells

Chimeric antigen receptor (CAR) T cell therapy is used in treating human hematological malignancies, but its efficacy is limited by T cell exhaustion (TEX). TEX arises at the expense of central memory T cells (TCM), which exhibit robust antitumor efficacy. Reduction of the TET2 gene led to increased TCM differentiation in a patient with leukemia who experienced a complete remission. We show that loss of TET2 led to increased chromatin accessibility at exhaustion regulators TOX and TOX2, plus increased expression of TOX2. Knockdown of TOX increased the percentage of TCM. However, unexpectedly, knockdown of TOX2 decreased TCM percentage and reduced proliferation. Consistently, a TCM gene signature was reduced in the TOX2 knockdown, and TOX2 bound to promoters of numerous TCM genes. Our results thus suggest a role for human TOX2, in contrast to exhaustion regulator TOX, as a potentiator of central memory differentiation of CAR T cells, with plausible utility in CAR T cell cancer therapy via modulated TOX2 expression.

Deconvolution of bulk gene expression profiles reveals the association between immune cell polarization and the prognosis of hepatocellular carcinoma patients

BACKGROUND

Many studies have utilized computational methods, including cell composition deconvolution (CCD), to correlate immune cell polarizations with the survival of cancer patients, including those with hepatocellular carcinoma (HCC). However, currently available cell deconvolution estimated (CDE) tools do not cover the wide range of immune cell changes that are known to influence tumor progression.

RESULTS

A new CCD tool, HCCImm, was designed to estimate the abundance of tumor cells and 16 immune cell types in the bulk gene expression profiles of HCC samples. HCCImm was validated using real datasets derived from human peripheral blood mononuclear cells (PBMCs) and HCC tissue samples, demonstrating that HCCImm outperforms other CCD tools. We used HCCImm to analyze the bulk RNA-seq datasets of The Cancer Genome Atlas (TCGA)-liver hepatocellular carcinoma (LIHC) samples. We found that the proportions of memory CD8+ T cells and Tregs were negatively associated with patient overall survival (OS). Furthermore, the proportion of naïve CD8+ T cells was positively associated with patient OS. In addition, the TCGA-LIHC samples with a high tumor mutational burden had a significantly high abundance of nonmacrophage leukocytes.

CONCLUSIONS

HCCImm was equipped with a new set of reference gene expression profiles that allowed for a more robust analysis of HCC patient expression data. The source code is provided at https://github.com/holiday01/HCCImm.

Identification of molecular mechanisms causing skin lesions of cutaneous leishmaniasis using weighted gene coexpression network analysis (WGCNA)

Leishmaniasis is an infectious disease, caused by a protozoan parasite. Its most common form is cutaneous leishmaniasis, which leaves scars on exposed body parts from bites by infected female phlebotomine sandflies. Approximately 50% of cases of cutaneous leishmaniasis fail to respond to standard treatments, creating slow-healing wounds which cause permanent scars on the skin. We performed a joint bioinformatics analysis to identify differentially expressed genes (DEGs) in healthy skin biopsies and Leishmania cutaneous wounds. DEGs and WGCNA modules were analyzed based on the Gene Ontology function, and the Cytoscape software. Among almost 16,600 genes that had significant expression changes on the skin surrounding Leishmania wounds, WGCNA determined that one of the modules, with 456 genes, has the strongest correlation with the size of the wounds. Functional enrichment analysis indicated that this module includes three gene groups with significant expression changes. These produce tissue-damaging cytokines or disrupt the production and activation of collagen, fibrin proteins, and the extracellular matrix, causing skin wounds or preventing them from healing. The hub genes of these groups are OAS1, SERPINH1, and FBLN1 respectively. This information can provide new ways to deal with unwanted and harmful effects of cutaneous leishmaniasis.

Single-cell multi-omic topic embedding reveals cell-type-specific and COVID-19 severity-related immune signatures

The advent of single-cell multi-omics sequencing technology makes it possible for re-searchers to leverage multiple modalities for individual cells and explore cell heterogeneity. However, the high dimensional, discrete, and sparse nature of the data make the downstream analysis particularly challenging. Most of the existing computational methods for single-cell data analysis are either limited to single modality or lack flexibility and interpretability. In this study, we propose an interpretable deep learning method called multi-omic embedded topic model (moETM) to effectively perform integrative analysis of high-dimensional single-cell multimodal data. moETM integrates multiple omics data via a product-of-experts in the encoder for efficient variational inference and then employs multiple linear decoders to learn the multi-omic signatures of the gene regulatory programs. Through comprehensive experiments on public single-cell transcriptome and chromatin accessibility data (i.e., scRNA+scATAC), as well as scRNA and proteomic data (i.e., CITE-seq), moETM demonstrates superior performance compared with six state-of-the-art single-cell data analysis methods on seven publicly available datasets. By applying moETM to the scRNA+scATAC data in human bone marrow mononuclear cells (BMMCs), we identified sequence motifs corresponding to the transcription factors that regulate immune gene signatures. Applying moETM analysis to CITE-seq data from the COVID-19 patients revealed not only known immune cell-type-specific signatures but also composite multi-omic biomarkers of critical conditions due to COVID-19, thus providing insights from both biological and clinical perspectives.

Revealing the immune cell subtype reconstitution profile in patients from the CLARITY study using deconvolution algorithms after cladribine tablets treatment

Immune Cell Deconvolution methods utilizing gene expression profiling to quantify immune cells in tissues and blood are an appealing alternative to flow cytometry. Our objective was to investigate the applicability of deconvolution approaches in clinical trial settings to better investigate the mode of action of drugs for autoimmune diseases. Popular deconvolution methods CIBERSORT and xCell were validated using gene expression from the publicly available GSE93777 dataset that has comprehensive matching flow cytometry. As shown in the online tool, ~ 50% of signatures show strong correlation (r > 0.5) with the remainder showing moderate correlation, or in a few cases, no correlation. Deconvolution methods were then applied to gene expression data from the phase III CLARITY study (NCT00213135) to evaluate the immune cell profile of relapsing multiple sclerosis patients treated with cladribine tablets. At 96 weeks after treatment, deconvolution scores showed the following changes vs placebo: naïve, mature, memory CD4⁺ and CD8⁺ T cells, non-class switched, and class switched memory B cells and plasmablasts were significantly reduced, naïve B cells and M2 macrophages were more abundant. Results confirm previously described changes in immune cell composition following cladribine tablets treatment and reveal immune homeostasis of pro- vs anti-inflammatory immune cell subtypes, potentially supporting long-term efficacy.

Predictive network analysis identifies JMJD6 and other potential key drivers in Alzheimer's disease

Despite decades of genetic studies on late-onset Alzheimer's disease, the underlying molecular mechanisms remain unclear. To better comprehend its complex etiology, we use an integrative approach to build robust predictive (causal) network models using two large human multi-omics datasets. We delineate bulk-tissue gene expression into single cell-type gene expression and integrate clinical and pathologic traits, single nucleotide variation, and deconvoluted gene expression for the construction of cell type-specific predictive network models. Here, we focus on neuron-specific network models and prioritize 19 predicted key drivers modulating Alzheimer's pathology, which we then validate by knockdown in human induced pluripotent stem cell-derived neurons. We find that neuronal knockdown of 10 of the 19 targets significantly modulates levels of amyloid-beta and/or phosphorylated tau peptides, most notably JMJD6. We also confirm our network structure by RNA sequencing in the neurons following knockdown of each of the 10 targets, which additionally predicts that they are upstream regulators of REST and VGF. Our work thus identifies robust neuronal key drivers of the Alzheimer's-associated network state which may represent therapeutic targets with relevance to both amyloid and tau pathology in Alzheimer's disease.

TSTEM-like CAR-T cells exhibit improved persistence and tumor control compared with conventional CAR-T cells in preclinical models

Patients who receive chimeric antigen receptor (CAR)-T cells that are enriched in memory T cells exhibit better disease control as a result of increased expansion and persistence of the CAR-T cells. Human memory T cells include stem-like CD8⁺ memory T cell progenitors that can become either functional stem-like T (T_STEM) cells or dysfunctional T progenitor exhausted (T_PEX) cells. To that end, we demonstrated that T_STEM cells were less abundant in infused CAR-T cell products in a phase 1 clinical trial testing Lewis Y-CAR-T cells (NCT03851146), and the infused CAR-T cells displayed poor persistence in patients. To address this issue, we developed a production protocol to generate T_STEM-like CAR-T cells enriched for expression of genes in cell replication pathways. Compared with conventional CAR-T cells, T_STEM-like CAR-T cells had enhanced proliferative capacity and increased cytokine secretion after CAR stimulation, including after chronic CAR stimulation in vitro. These responses were dependent on the presence of CD4⁺ T cells during T_STEM-like CAR-T cell production. Adoptive transfer of T_STEM-like CAR-T cells induced better control of established tumors and resistance to tumor rechallenge in preclinical models. These more favorable outcomes were associated with increased persistence of T_STEM-like CAR-T cells and an increased memory T cell pool. Last, T_STEM-like CAR-T cells and anti-programmed cell death protein 1 (PD-1) treatment eradicated established tumors, and this was associated with increased tumor-infiltrating CD8⁺CAR⁺ T cells producing interferon-γ. In conclusion, our CAR-T cell protocol generated T_STEM-like CAR-T cells with enhanced therapeutic efficacy, resulting in increased proliferative capacity and persistence in vivo.

Clonal associations of lymphocyte subsets and functional states revealed by single cell antigen receptor profiling of T and B cells in rheumatoid arthritis synovium

Rheumatoid arthritis (RA) is an autoimmune disease initiated by antigen-specific T cells and B cells, which promote synovial inflammation through a complex set of interactions with innate immune and stromal cells. To better understand the phenotypes and clonal relationships of synovial T and B cells, we performed single-cell RNA and repertoire sequencing on paired synovial tissue and peripheral blood samples from 12 donors with seropositive RA ranging from early to chronic disease. Paired transcriptomic-repertoire analyses highlighted 3 clonally distinct CD4 T cells populations that were enriched in RA synovium: T peripheral helper (Tph) and T follicular helper (Tfh) cells, CCL5+ T cells, and T regulatory cells (Tregs). Among these cells, Tph cells showed a unique transcriptomic signature of recent T cell receptor (TCR) activation, and clonally expanded Tph cells expressed an elevated transcriptomic effector signature compared to non-expanded Tph cells. CD8 T cells showed higher oligoclonality than CD4 T cells, and the largest CD8 T cell clones in synovium were highly enriched in GZMK+ cells. TCR analyses revealed CD8 T cells with likely viral-reactive TCRs distributed across transcriptomic clusters and definitively identified MAIT cells in synovium, which showed transcriptomic features of TCR activation. Among B cells, non-naive B cells including age-associated B cells (ABC), NR4A1+ activated B cells, and plasma cells, were enriched in synovium and had higher somatic hypermutation rates compared to blood B cells. Synovial B cells demonstrated substantial clonal expansion, with ABC, memory, and activated B cells clonally linked to synovial plasma cells. Together, these results reveal clonal relationships between functionally distinct lymphocyte populations that infiltrate RA synovium.

CeDAR: incorporating cell type hierarchy improves cell type-specific differential analyses in bulk omics data

Bulk high-throughput omics data contain signals from a mixture of cell types. Recent developments of deconvolution methods facilitate cell type-specific inferences from bulk data. Our real data exploration suggests that differential expression or methylation status is often correlated among cell types. Based on this observation, we develop a novel statistical method named CeDAR to incorporate the cell type hierarchy in cell type-specific differential analyses of bulk data. Extensive simulation and real data analyses demonstrate that this approach significantly improves the accuracy and power in detecting cell type-specific differential signals compared with existing methods, especially in low-abundance cell types.

The genetic and evolutionary basis of gene expression variation in East Africans

BACKGROUND

Mapping of quantitative trait loci (QTL) associated with molecular phenotypes is a powerful approach for identifying the genes and molecular mechanisms underlying human traits and diseases, though most studies have focused on individuals of European descent. While important progress has been made to study a greater diversity of human populations, many groups remain unstudied, particularly among indigenous populations within Africa. To better understand the genetics of gene regulation in East Africans, we perform expression and splicing QTL mapping in whole blood from a cohort of 162 diverse Africans from Ethiopia and Tanzania. We assess replication of these QTLs in cohorts of predominantly European ancestry and identify candidate genes under selection in human populations.

RESULTS

We find the gene regulatory architecture of African and non-African populations is broadly shared, though there is a considerable amount of variation at individual loci across populations. Comparing our analyses to an equivalently sized cohort of European Americans, we find that QTL mapping in Africans improves the detection of expression QTLs and fine-mapping of causal variation. Integrating our QTL scans with signatures of natural selection, we find several genes related to immunity and metabolism that are highly differentiated between Africans and non-Africans, as well as a gene associated with pigmentation.

CONCLUSION

Extending QTL mapping studies beyond European ancestry, particularly to diverse indigenous populations, is vital for a complete understanding of the genetic architecture of human traits and can reveal novel functional variation underlying human traits and disease.

Tissue-specific abundance of interferon-gamma drives regulatory T cells to restrain DC1-mediated priming of cytotoxic T cells against lung cancer

Local environmental factors influence CD8+ T cell priming in lymph nodes (LNs). Here, we sought to understand how factors unique to the tumor-draining mediastinal LN (mLN) impact CD8+ T cell responses toward lung cancer. Type 1 conventional dendritic cells (DC1s) showed a mLN-specific failure to induce robust cytotoxic T cells responses. Using regulatory T (Treg) cell depletion strategies, we found that Treg cells suppressed DC1s in a spatially coordinated manner within tissue-specific microniches within the mLN. Treg cell suppression required MHC II-dependent contact between DC1s and Treg cells. Elevated levels of IFN-γ drove differentiation Treg cells into Th1-like effector Treg cells in the mLN. In patients with cancer, Treg cell Th1 polarization, but not CD8+/Treg cell ratios, correlated with poor responses to checkpoint blockade immunotherapy. Thus, IFN-γ in the mLN skews Treg cells to be Th1-like effector Treg cells, driving their close interaction with DC1s and subsequent suppression of cytotoxic T cell responses.

Bioinformatics Tools for Bulk Gene Expression Deconvolution in Diabetic Retinopathy

Retinal neovascularization is one of the leading causes of vision loss and a hallmark of proliferative diabetic retinopathy (PDR). The immune system is observed to be involved in the pathogenesis of diabetic retinopathy (DR). The specific immune cell type that contributes to retinal neovascularization can be identified via a bioinformatics analysis of RNA sequencing (RNA-seq) data, known as deconvolution analysis. Previous study has identified the infiltration of macrophages in the retina of rats with hypoxia-induced retinal neovascularization and patients with PDR through a deconvolution algorithm, known as CIBERSORTx. Here, we describe the protocols of using CIBERSORTx to perform the deconvolution analysis and downstream analysis of RNA-seq data.

A large-scale screening and functional sorting of tumour microenvironment prognostic genes for breast cancer patients

PURPOSE

The aim of this study was to systematically establish a comprehensive tumour microenvironment (TME)-relevant prognostic gene and target miRNA network for breast cancer patients.

METHODS

Based on a large-scale screening of TME-relevant prognostic genes (760 genes) for breast cancer patients, the prognostic model was established. The primary TME prognostic genes were selected from the constructing database and verified in the testing database. The internal relationships between the potential TME prognostic genes and the prognosis of breast cancer patients were explored in depth. The associated miRNAs for the TME prognostic genes were generated, and the functions of each primary TME member were investigated in the breast cancer cell line.

RESULTS

Compared with sibling controls, breast cancer patients showed 55 differentially expressed TME prognostic genes, of which 31 were considered as protective genes, while the remaining 24 genes were considered as risk genes. According to the lambda values of the LASSO Cox analysis, the 15 potential TME prognostic genes were as follows: ENPEP, CCDC102B, FEZ1, NOS2, SCG2, RPLP2, RELB, RGS3, EMP1, PDLIM4, EPHA3, PCDH9, VIM, GFI1, and IRF1. Among these, there was a remarkable linear internal relationship for CCDC102B but non-linear relationships for others with breast cancer patient prognosis. Using the siRNA technique, we silenced the expression of each TME prognostic gene. Seven of the 15 TME prognostic genes (NOS2, SCG2, RGS3, EMP1, PDLIM4, PCDH9, and GFI1) were involved in enhancing cell proliferation, destroying cell apoptosis, promoting cell invasion, or migration in breast cancer. Six of them (CCDC102B, RPLP2, RELB, EPHA3, VIM, and IRF1) were favourable for maintaining cell invasion or migration. Only two of them (ENPEP and FEZ1) were favourable for the processes of cell proliferation and apoptosis.

CONCLUSIONS

This integrated study hypothesised an innovative TME-associated genetic functional network for breast cancer patients. The external relationships between these TME prognostic genes and the disease were measured. Meanwhile, the internal molecular mechanisms were also investigated.

Peripheral blood RNA biomarkers for cardiovascular disease from bench to bedside: a position paper from the EU-CardioRNA COST action CA17129

Despite significant advances in the diagnosis and treatment of cardiovascular diseases, recent calls have emphasized the unmet need to improve precision-based approaches in cardiovascular disease. Although some studies provide preliminary evidence of the diagnostic and prognostic potential of circulating coding and non-coding RNAs, the complex RNA biology and lack of standardization have hampered the translation of these markers into clinical practice. In this position paper of the CardioRNA COST action CA17129, we provide recommendations to standardize the RNA development process in order to catalyse efforts to investigate novel RNAs for clinical use. We list the unmet clinical needs in cardiovascular disease, such as the identification of high-risk patients with ischaemic heart disease or heart failure who require more intensive therapies. The advantages and pitfalls of the different sample types, including RNAs from plasma, extracellular vesicles, and whole blood, are discussed in the sample matrix, together with their respective analytical methods. The effect of patient demographics and highly prevalent comorbidities, such as metabolic disorders, on the expression of the candidate RNA is presented and should be reported in biomarker studies. We discuss the statistical and regulatory aspects to translate a candidate RNA from a research use only assay to an in-vitro diagnostic test for clinical use. Optimal planning of this development track is required, with input from the researcher, statistician, industry, and regulatory partners.

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.

Multiset multicover methods for discriminative marker selection

Markers are increasingly being used for several high-throughput data analysis and experimental design tasks. Examples include the use of markers for assigning cell types in scRNA-seq studies, for deconvolving bulk gene expression data, and for selecting marker proteins in single-cell spatial proteomics studies. Most marker selection methods focus on differential expression (DE) analysis. Although such methods work well for data with a few non-overlapping marker sets, they are not appropriate for large atlas-size datasets where several cell types and tissues are considered. To address this, we define the phenotype cover (PC) problem for marker selection and present algorithms that can improve the discriminative power of marker sets. Analysis of these sets on several marker-selection tasks suggests that these methods can lead to solutions that accurately distinguish different phenotypes in the data.

Knockdown of TTC9 inhibits the proliferation, migration and invasion, but induces the apoptosis of lung adenocarcinoma cells

Lung adenocarcinoma (LUAD) is one of the most commonly diagnosed subtypes of lung cancer, and one of the deadliest cancers. Tetratricopeptide repeat domain 9A (TTC9) is upregulated and has played an oncogenic role in some malignant tumors. However, the expression and role of TTC9 has not yet been elucidated in LUAD. Here, we investigated the expression profiles, biological functions and potential molecular mechanism of the TTC9 gene in LUAD. TTC9 expression was significantly overexpressed in LUAD tissues compared with that in normal lung tissues. TTC9 expression was closely correlated with gender, lymph node metastasis, and survival status in the TCGA-LUAD cohort. Subsequent cellular function assays demonstrated that knockdown of TTC9 promoted PC9 cell apoptosis and inhibited cell proliferation, migration and invasion, leading to cell cycle arrest in G2 phase. Moreover, inhibition of TTC9 suppressed the tumorigenicity of PC9 cells in nude mice. TTC9 might serve as oncogene in LUAD through cancer-related signaling pathways including p38 MAPK pathway. The expression of TTC9 gene might be modulated by DNA copy number variant and DNA methylation. TTC9 was significantly associated with tumor immune infiltration patterns. Accordingly, TTC9 may be a novel therapeutic target for the treatment of LUAD.

Tumor microenvironment: barrier or opportunity towards effective cancer therapy

Tumor microenvironment (TME) is a specialized ecosystem of host components, designed by tumor cells for successful development and metastasis of tumor. With the advent of 3D culture and advanced bioinformatic methodologies, it is now possible to study TME’s individual components and their interplay at higher resolution. Deeper understanding of the immune cell’s diversity, stromal constituents, repertoire profiling, neoantigen prediction of TMEs has provided the opportunity to explore the spatial and temporal regulation of immune therapeutic interventions. The variation of TME composition among patients plays an important role in determining responders and non-responders towards cancer immunotherapy. Therefore, there could be a possibility of reprogramming of TME components to overcome the widely prevailing issue of immunotherapeutic resistance. The focus of the present review is to understand the complexity of TME and comprehending future perspective of its components as potential therapeutic targets. The later part of the review describes the sophisticated 3D models emerging as valuable means to study TME components and an extensive account of advanced bioinformatic tools to profile TME components and predict neoantigens. Overall, this review provides a comprehensive account of the current knowledge available to target TME.

GEOMETRIC STRUCTURE GUIDED MODEL AND ALGORITHMS FOR COMPLETE DECONVOLUTION OF GENE EXPRESSION DATA

Complete deconvolution analysis for bulk RNA-seq data is important and helpful to distinguish whether the differences of disease-associated GEPs (gene expression profiles) in tissues of patients and normal controls are due to changes in cellular composition of tissue samples, or due to GEPs changes in specific cells. One of the major techniques to perform complete deconvolution is nonnegative matrix factorization (NMF), which also has a wide-range of applications in the machine learning community. However, the NMF is a well-known strongly ill-posed problem, so a direct application of NMF to RNA-seq data will suffer severe difficulties in the interpretability of solutions. In this paper, we develop an NMF-based mathematical model and corresponding computational algorithms to improve the solution identifiability of deconvoluting bulk RNA-seq data. In our approach, we combine the biological concept of marker genes with the solvability conditions of the NMF theories, and develop a geometric structures guided optimization model. In this strategy, the geometric structure of bulk tissue data is first explored by the spectral clustering technique. Then, the identified information of marker genes is integrated as solvability constraints, while the overall correlation graph is used as manifold regularization. Both synthetic and biological data are used to validate the proposed model and algorithms, from which solution interpretability and accuracy are significantly improved.

Transcriptional profiles predict treatment outcome in patients with tuberculosis and diabetes at diagnosis and at two weeks after initiation of anti-tuberculosis treatment

BACKGROUND

Globally, the tuberculosis (TB) treatment success rate is approximately 85%, with treatment failure, relapse and death occurring in a significant proportion of pulmonary TB patients. Treatment success is lower among people with diabetes mellitus (DM). Predicting treatment outcome early after diagnosis, especially in TB-DM patients, would allow early treatment adaptation for individuals and may improve global TB control.

METHODS

Samples were collected in a longitudinal cohort study of adult TB patients from South Africa (n  =  94) and Indonesia (n  =  81), who had concomitant DM (n  =  59), intermediate hyperglycaemia (n  =  79) or normal glycaemia/no DM (n  =  37). Treatment outcome was monitored, and patients were categorized as having a good (cured) or poor (failed, recurrence, died) outcome during treatment and 12 months follow-up. Whole blood transcriptional profiles before, during and at the end of TB treatment were characterized using unbiased RNA-Seq and targeted gene dcRT-MLPA.

FINDINGS

We report differences in whole blood transcriptome profiles, which were observed before initiation of treatment and throughout treatment, between patients with a good versus poor TB treatment outcome. An eight-gene and a 22-gene blood transcriptional signature distinguished patients with a good TB treatment outcome from patients with a poor TB treatment outcome at diagnosis (AUC = 0·815) or two weeks (AUC = 0·834) after initiation of TB treatment, respectively. High accuracy was obtained by cross-validating this signature in an external cohort (AUC = 0·749).

INTERPRETATION

These findings suggest that transcriptional profiles can be used as a prognostic biomarker for treatment failure and success, even in patients with concomitant DM.

FUNDING

The research leading to these results, as part of the TANDEM Consortium, received funding from the European Community's Seventh Framework Programme (FP7/2007-2013 Grant Agreement No. 305279) and the Netherlands Organization for Scientific Research (NWO-TOP Grant Agreement No. 91214038). The research leading to the results presented in the Indian validation cohort was supported by Research Council of Norway Global Health and Vaccination Research (GLOBVAC) projects: RCN 179342, 192534, and 248042, the University of Bergen (Norway).

More accurate estimation of cell composition in bulk expression through robust integration of single-cell information

Motivation

The rapid single-cell transcriptomic technology developments have led to an increasing interest in cellular heterogeneity within cell populations. Although cell-type proportions can be obtained directly from single-cell RNA sequencing (scRNA-seq), it is costly and not feasible in every study. Alternatively, with fewer experimental complications, cell-type compositions are characterized from bulk RNA-seq data. Many computational tools have been developed and reported in the literature. However, they fail to appropriately incorporate the covariance structures in both scRNA-seq and bulk RNA-seq datasets in use.

Results

We present a covariance-based single-cell decomposition (CSCD) method that estimates cell-type proportions in bulk data through building a reference expression profile based on a single-cell data, and learning gene-specific bulk expression transformations using a constrained linear inverse model. The approach is similar to Bisque, a cell-type decomposition method that was recently developed. Bisque is limited to a univariate model, thus unable to incorporate gene-gene correlations into the analysis. We introduce a more advanced model that successfully incorporates the covariance structures in both scRNA-seq and bulk RNA-seq datasets into the analysis, and fixes the collinearity issue by utilizing a linear shrinkage estimation of the corresponding covariance matrices. We applied CSCD to several publicly available datasets and measured the performance of CSCD, Bisque and six other common methods in the literature. Our results indicate that CSCD is more accurate and comprehensive than most of the existing methods.

Availability and implementation

The R package is available on https://github.com/empiricalbayes/CSCDRNA.

Tumor-induced erythroid precursor-differentiated myeloid cells mediate immunosuppression and curtail anti-PD-1/PD-L1 treatment efficacy

Despite the unprecedented success of immune checkpoint inhibitors (ICIs) as anti-cancer therapy, it remains a prevailing clinical need to identify additional mechanisms underlying ICI therapeutic efficacy and potential drug resistance. Here, using lineage tracking in cancer patients and tumor-bearing mice, we demonstrate that erythroid progenitor cells lose their developmental potential and switch to the myeloid lineage. Single-cell transcriptome analyses reveal that, notwithstanding quantitative differences in erythroid gene expression, erythroid differentiated myeloid cells (EDMCs) are transcriptionally indistinguishable from their myeloid-originated counterparts. EDMCs possess multifaceted machinery to curtail T cell-mediated anti-tumor responses. Consequently, EDMC content within tumor tissues is negatively associated with T cell inflammation for the majority of solid cancers; moreover, EDMC enrichment, in accordance with anemia manifestation, is predictive of poor prognosis in various cohorts of patients undergoing ICI therapy. Together, our findings reveal a feedforward mechanism by which tumors exploit anemia-triggered erythropoiesis for myeloid transdifferentiation and immunosuppression.

TumorDecon: A digital cytometry software

There are many experimental methods for characterizing immune profiles of tumors, such as flow and mass cytometry. However, these approaches are time and resource intensive. Thus, several "digital cytometry" methods have been developed to extract cell frequencies from RNA-seq data. Here, we introduce TumorDecon, named for its potential to deconvolve the distribution of cells from the gene expression levels of a bulk of cells, such as a tumor. The Python package provides an accessible way of applying these methods. It includes four deconvolution methods as well as several gene sets, signature matrices, and functions for generating custom signature matrices.

Macrophage inflammation resolution requires CPEB4-directed offsetting of mRNA degradation

Chronic inflammation is a major cause of disease. Inflammation resolution is in part directed by the differential stability of mRNAs encoding pro-inflammatory and anti-inflammatory factors. In particular, tristetraprolin (TTP)-directed mRNA deadenylation destabilizes AU-rich element (ARE)-containing mRNAs. However, this mechanism alone cannot explain the variety of mRNA expression kinetics that are required to uncouple degradation of pro-inflammatory mRNAs from the sustained expression of anti-inflammatory mRNAs. Here, we show that the RNA-binding protein CPEB4 acts in an opposing manner to TTP in macrophages: it helps to stabilize anti-inflammatory transcripts harboring cytoplasmic polyadenylation elements (CPEs) and AREs in their 3′-UTRs, and it is required for the resolution of the lipopolysaccharide (LPS)-triggered inflammatory response. Coordination of CPEB4 and TTP activities is sequentially regulated through MAPK signaling. Accordingly, CPEB4 depletion in macrophages impairs inflammation resolution in an LPS-induced sepsis model. We propose that the counterbalancing actions of CPEB4 and TTP, as well as the distribution of CPEs and AREs in their target mRNAs, define transcript-specific decay patterns required for inflammation resolution. Thus, these two opposing mechanisms provide a fine-tuning control of inflammatory transcript destabilization while maintaining the expression of the negative feedback loops required for efficient inflammation resolution; disruption of this balance can lead to disease.

Decosus: An R Framework for Universal Integration of Cell Proportion Estimation Methods

The assessment of the cellular heterogeneity and abundance in bulk tissue samples is essential for characterising cellular and organismal states. Computational approaches to estimate cellular abundance from bulk RNA-Seq datasets have variable performances, often requiring benchmarking matrices to select the best performing methods for individual studies. However, such benchmarking investigations are difficult to perform and assess in typical applications because of the absence of gold standard/ground-truth cellular measurements. Here we describe Decosus, an R package that integrates seven methods and signatures for deconvoluting cell types from gene expression profiles (GEP). Benchmark analysis on a range of datasets with ground-truth measurements revealed that our integrated estimates consistently exhibited stable performances across datasets than individual methods and signatures. We further applied Decosus to characterise the immune compartment of skin samples in different settings, confirming the well-established Th1 and Th2 polarisation in psoriasis and atopic dermatitis, respectively. Secondly, we revealed immune system-related UV-induced changes in sun-exposed skin. Furthermore, a significant motivation in the design of Decosus is flexibility and the ability for the user to include new gene signatures, algorithms, and integration methods at run time.

Comprehensive evaluation of deconvolution methods for human brain gene expression

Transcriptome deconvolution aims to estimate the cellular composition of an RNA sample from its gene expression data, which in turn can be used to correct for composition differences across samples. The human brain is unique in its transcriptomic diversity, and comprises a complex mixture of cell-types, including transcriptionally similar subtypes of neurons. Here, we carry out a comprehensive evaluation of deconvolution methods for human brain transcriptome data, and assess the tissue-specificity of our key observations by comparison with human pancreas and heart. We evaluate eight transcriptome deconvolution approaches and nine cell-type signatures, testing the accuracy of deconvolution using in silico mixtures of single-cell RNA-seq data, RNA mixtures, as well as nearly 2000 human brain samples. Our results identify the main factors that drive deconvolution accuracy for brain data, and highlight the importance of biological factors influencing cell-type signatures, such as brain region and in vitro cell culturing.

Transcriptome deconvolution aims to estimate cellular composition based on gene expression data. Here the authors evaluate deconvolution methods for human brain transcriptome and conclude that partial deconvolution algorithms work best, but that appropriate cell-type signatures are also important.

Whole-Tissue Deconvolution and scRNAseq Analysis Identify Altered Endometrial Cellular Compositions and Functionality Associated With Endometriosis

The uterine lining (endometrium) exhibits a pro-inflammatory phenotype in women with endometriosis, resulting in pain, infertility, and poor pregnancy outcomes. The full complement of cell types contributing to this phenotype has yet to be identified, as most studies have focused on bulk tissue or select cell populations. Herein, through integrating whole-tissue deconvolution and single-cell RNAseq, we comprehensively characterized immune and nonimmune cell types in the endometrium of women with or without disease and their dynamic changes across the menstrual cycle. We designed metrics to evaluate specificity of deconvolution signatures that resulted in single-cell identification of 13 novel signatures for immune cell subtypes in healthy endometrium. Guided by statistical metrics, we identified contributions of endometrial epithelial, endothelial, plasmacytoid dendritic cells, classical dendritic cells, monocytes, macrophages, and granulocytes to the endometrial pro-inflammatory phenotype, underscoring roles for nonimmune as well as immune cells to the dysfunctionality of this tissue.

A2Sign: Agnostic Algorithms for Signatures-a universal method for identifying molecular signatures from transcriptomic datasets prior to cell-type deconvolution

MOTIVATION

Molecular signatures are critical for inferring the proportions of cell types from bulk transcriptomics data. However, the identification of these signatures is based on a methodology that relies on prior biological knowledge of the cell types being studied. When working with less known biological material, a data-driven approach is required to uncover the underlying classes and generate ad hoc signatures from healthy or pathogenic tissue.

RESULTS

We present a new approach, A2Sign: Agnostic Algorithms for Signatures, based on a non-negative tensor factorization (NTF) strategy that allows us to identify cell-type-specific molecular signatures, greatly reduce collinearities and also account for inter-individual variability. We propose a global framework that can be applied to uncover molecular signatures for cell-type deconvolution in arbitrary tissues using bulk transcriptome data. We also present two new molecular signatures for deconvolution of up to 16 immune cell types using microarray or RNA-seq data.

AVAILABILITY AND IMPLEMENTATION

All steps of our analysis were implemented in annotated Python notebooks (https://github.com/paulfogel/A2SIGN). To perform NTF, we used the NMTF package, which can be downloaded using Python pip install.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Identification and Validation of Pathogenic Genes in Sepsis and Associated Diseases by Integrated Bioinformatics Approach

Sepsis is a clinical syndrome with high mortality and morbidity rates. In sepsis, the abrupt release of cytokines by the innate immune system may cause multiorgan failure, leading to septic shock and associated complications. In the presence of a number of systemic disorders, such as sepsis, infections, diabetes, and systemic lupus erythematosus (SLE), cardiorenal syndrome (CRS) type 5 is defined by concomitant cardiac and renal dysfunctions Thus, our study suggests that certain mRNAs and unexplored pathways may pave a way to unravel critical therapeutic targets in three debilitating and interrelated illnesses, namely, sepsis, SLE, and CRS. Sepsis, SLE, and CRS are closely interrelated complex diseases likely sharing an overlapping pathogenesis caused by erroneous gene network activities. We sought to identify the shared gene networks and the key genes for sepsis, SLE, and CRS by completing an integrative analysis. Initially, 868 DEGs were identified in 16 GSE datasets. Based on degree centrality, 27 hub genes were revealed. The gProfiler webtool was used to perform functional annotations and enriched molecular pathway analyses. Finally, core hub genes (EGR1, MMP9, and CD44) were validated using RT-PCR analysis. Our comprehensive multiplex network approach to hub gene discovery is effective, as evidenced by the findings. This work provides a novel research path for a new research direction in multi-omics biological data analysis.

Besca, a single-cell transcriptomics analysis toolkit to accelerate translational research

Single-cell RNA sequencing (scRNA-seq) revolutionized our understanding of disease biology. The promise it presents to also transform translational research requires highly standardized and robust software workflows. Here, we present the toolkit Besca, which streamlines scRNA-seq analyses and their use to deconvolute bulk RNA-seq data according to current best practices. Beyond a standard workflow covering quality control, filtering, and clustering, two complementary Besca modules, utilizing hierarchical cell signatures and supervised machine learning, automate cell annotation and provide harmonized nomenclatures. Subsequently, the gene expression profiles can be employed to estimate cell type proportions in bulk transcriptomics data. Using multiple, diverse scRNA-seq datasets, some stemming from highly heterogeneous tumor tissue, we show how Besca aids acceleration, interoperability, reusability and interpretability of scRNA-seq data analyses, meeting crucial demands in translational research and beyond.

How many markers are needed to robustly determine a cell's type?

Our understanding of cell types has advanced considerably with the publication of single-cell atlases. Marker genes play an essential role for experimental validation and computational analyses such as physiological characterization, annotation, and deconvolution. However, a framework for quantifying marker replicability and selecting replicable markers is currently lacking. Here, using high-quality data from the Brain Initiative Cell Census Network (BICCN), we systematically investigate marker replicability for 85 neuronal cell types. We show that, due to dataset-specific noise, we need to combine 5 datasets to obtain robust differentially expressed (DE) genes, particularly for rare populations and lowly expressed genes. We estimate that 10 to 200 meta-analytic markers provide optimal downstream performance and make available replicable marker lists for the 85 BICCN cell types. Replicable marker lists condense interpretable and generalizable information about cell types, opening avenues for downstream applications, including cell type annotation, selection of gene panels, and bulk data deconvolution.