Multi-Omics Profiling Reveals Phenotypic and Functional Heterogeneity of Neutrophils in COVID-19

Accumulating evidence has revealed unexpected phenotypic heterogeneity and diverse functions of neutrophils in several diseases. Coronavirus disease (COVID-19) can alter the leukocyte phenotype based on disease severity, including neutrophil activation in severe cases. However, the plasticity of neutrophil phenotypes and their relative impact on COVID-19 pathogenesis has not been well addressed. This study aimed to identify and validate the heterogeneity of neutrophils in COVID-19 and evaluate the functions of each subpopulation. We analyzed public single-cell RNA-seq, bulk RNA-seq, and proteome data from healthy donors and patients with COVID-19 to investigate neutrophil subpopulations and their response to disease pathogenesis. We identified eight neutrophil subtypes: pro-neutrophil, pre-neutrophil, immature neutrophil, and five mature neutrophil subpopulations. The subtypes exhibited distinct features, including diverse activation signatures and multiple enriched pathways. The pro-neutrophil subtype was associated with severe and fatal disease, while the pre-neutrophil subtype was particularly abundant in mild/moderate disease. One of the mature neutrophil subtypes showed consistently large fractions in patients with different disease severity. Bulk RNA-seq dataset analyses using a cellular deconvolution approach validated the relative abundances of neutrophil subtypes and the expansion of pro-neutrophils in severe COVID-19 patients. Cell-cell communication analysis revealed representative ligand-receptor interactions among the identified neutrophil subtypes. Further investigation into transcription factors and differential protein abundance revealed the regulatory network differences between healthy donors and patients with severe COVID-19. Overall, we demonstrated the complex interactions among heterogeneous neutrophil subtypes and other blood cell types during COVID-19 disease. Our work has great value in terms of both clinical and public health as it furthers our understanding of the phenotypic and functional heterogeneity of neutrophils and other cell populations in multiple diseases.

Cellular abundance-based prognostic model associated with deregulated gene expression of leukemic stem cells in acute myeloid leukemia

Background: Previous studies have reported that genes highly expressed in leukemic stem cells (LSC) may dictate the survival probability of patients and expression-based cellular deconvolution may be informative in forecasting prognosis. However, whether the prognosis of acute myeloid leukemia (AML) can be predicted using gene expression and deconvoluted cellular abundances is debatable. Methods: Nine different cell-type abundances of a training set composed of the AML samples of 422 patients, were used to build a model for predicting prognosis by least absolute shrinkage and selection operator Cox regression. This model was validated in two different validation sets, TCGA-LAML and Beat AML (n = 179 and 451, respectively). Results: We introduce a new prognosis predicting model for AML called the LSC activity (LSCA) score, which incorporates the abundance of 5 cell types, granulocyte-monocyte progenitors, common myeloid progenitors, CD45RA + cells, megakaryocyte-erythrocyte progenitors, and multipotent progenitors. Overall survival probabilities between the high and low LSCA score groups were significantly different in TCGA-LAML and Beat AML cohorts (log-rank p-value = 3.3 × 10 - 4 and 4.3 × 10 - 3 , respectively). Also, multivariate Cox regression analysis on these two validation sets shows that LSCA score is independent prognostic factor when considering age, sex, and cytogenetic risk (hazard ratio, HR = 2.17; 95% CI 1.40-3.34; p < 0.001 and HR = 1.20; 95% CI 1.02-1.43; p < 0.03, respectively). The performance of the LSCA score was comparable to other prognostic models, LSC17, APS, and CTC scores, as indicated by the area under the curve. Gene set variation analysis with six LSC-related functional gene sets indicated that high and low LSCA scores are associated with upregulated and downregulated genes in LSCs. Conclusion: We have developed a new prognosis prediction scoring system for AML patients, the LSCA score, which uses deconvoluted cell-type abundance only.

Calculating detection limits and uncertainty of reference-based deconvolution of whole-blood DNA methylation data

DNA methylation (DNAm)-based cell mixture deconvolution (CMD) has become a quintessential part of epigenome-wide association studies where DNAm is profiled in heterogeneous tissue types. Despite being introduced over a decade ago, detection limits, which represent the smallest fraction of a cell type in a mixed biospecimen that can be reliably detected, have yet to be determined in the context of DNAm-based CMD. Moreover, there has been little attention given to approaches for quantifying the uncertainty associated with DNAm-based CMD. Here, analytical frameworks for determining both cell-specific limits of detection and quantification of uncertainty associated with DNAm-based CMD are described. This work may contribute to improved rigor, reproducibility and replicability of epigenome-wide association studies involving CMD.

Characterization of Continuous Transcriptional Heterogeneity in High-Risk Blastemal-Type Wilms' Tumors Using Unsupervised Machine Learning

Wilms' tumors are pediatric malignancies that are thought to arise from faulty kidney development. They contain a wide range of poorly differentiated cell states resembling various distorted developmental stages of the fetal kidney, and as a result, differ between patients in a continuous manner that is not well understood. Here, we used three computational approaches to characterize this continuous heterogeneity in high-risk blastemal-type Wilms' tumors. Using Pareto task inference, we show that the tumors form a triangle-shaped continuum in latent space that is bounded by three tumor archetypes with "stromal", "blastemal", and "epithelial" characteristics, which resemble the un-induced mesenchyme, the cap mesenchyme, and early epithelial structures of the fetal kidney. By fitting a generative probabilistic "grade of membership" model, we show that each tumor can be represented as a unique mixture of three hidden "topics" with blastemal, stromal, and epithelial characteristics. Likewise, cellular deconvolution allows us to represent each tumor in the continuum as a unique combination of fetal kidney-like cell states. These results highlight the relationship between Wilms' tumors and kidney development, and we anticipate that they will pave the way for more quantitative strategies for tumor stratification and classification.

Smoking induces shifts in cellular composition and transcriptome within the bronchial mucus barrier

BACKGROUND AND OBJECTIVE

Smoking disturbs the bronchial-mucus-barrier. This study assesses the cellular composition and gene expression shifts of the bronchial-mucus-barrier with smoking to understand the mechanism of mucosal damage by cigarette smoke exposure. We explore whether single-cell-RNA-sequencing (scRNA-seq) based cellular deconvolution (CD) can predict cell-type composition in RNA-seq data.

METHODS

RNA-seq data of bronchial biopsies from three cohorts were analysed using CD. The cohorts included 56 participants with chronic obstructive pulmonary disease [COPD] (38 smokers; 18 ex-smokers), 77 participants without COPD (40 never-smokers; 37 smokers) and 16 participants who stopped smoking for 1 year (11 COPD and 5 non-COPD-smokers). Differential gene expression was used to investigate gene expression shifts. The CD-derived goblet cell ratios were validated by correlating with staining-derived goblet cell ratios from the COPD cohort. Statistics were done in the R software (false discovery rate p-value < 0.05).

RESULTS

Both CD methods indicate a shift in bronchial-mucus-barrier cell composition towards goblet cells in COPD and non-COPD-smokers compared to ex- and never-smokers. It shows that the effect was reversible within a year of smoking cessation. A reduction of ciliated and basal cells was observed with current smoking, which resolved following smoking cessation. The expression of mucin and sodium channel (ENaC) genes, but not chloride channel genes, were altered in COPD and current smokers compared to never smokers or ex-smokers. The goblet cell-derived staining scores correlate with CD-derived goblet cell ratios.

CONCLUSION

Smoking alters bronchial-mucus-barrier cell composition, transcriptome and increases mucus production. This effect is partly reversible within a year of smoking cessation. CD methodology can predict goblet-cell percentages from RNA-seq.

Spatially resolved deconvolution of the fibrotic niche in lung fibrosis

A defining pathological feature of human lung fibrosis is localized tissue heterogeneity, which challenges the interpretation of transcriptomic studies that typically lose spatial information. Here we investigate spatial gene expression in diagnostic tissue using digital profiling technology. We identify distinct, region-specific gene expression signatures as well as shared gene signatures. By integration with single-cell data, we spatially map the cellular composition within and distant from the fibrotic niche, demonstrating discrete changes in homeostatic and pathologic cell populations even in morphologically preserved lung, while through ligand-receptor analysis, we investigate cellular cross-talk within the fibrotic niche. We confirm findings through bioinformatic, tissue, and in vitro analyses, identifying that loss of NFKB inhibitor zeta in alveolar epithelial cells dysregulates the TGFβ/IL-6 signaling axis, which may impair homeostatic responses to environmental stress. Thus, spatially resolved deconvolution advances understanding of cell composition and microenvironment in human lung fibrogenesis.

Unraveling Targetable Systemic and Cell-Type-Specific Molecular Phenotypes of Alzheimer's and Parkinson's Brains With Digital Cytometry

Alzheimer's disease (AD) and Parkinson's disease (PD) are the two most common neurodegenerative disorders worldwide, with age being their major risk factor. The increasing worldwide life expectancy, together with the scarcity of available treatment choices, makes it thus pressing to find the molecular basis of AD and PD so that the causing mechanisms can be targeted. To study these mechanisms, gene expression profiles have been compared between diseased and control brain tissues. However, this approach is limited by mRNA expression profiles derived for brain tissues highly reflecting their degeneration in cellular composition but not necessarily disease-related molecular states. We therefore propose to account for cell type composition when comparing transcriptomes of healthy and diseased brain samples, so that the loss of neurons can be decoupled from pathology-associated molecular effects. This approach allowed us to identify genes and pathways putatively altered systemically and in a cell-type-dependent manner in AD and PD brains. Moreover, using chemical perturbagen data, we computationally identified candidate small molecules for specifically targeting the profiled AD/PD-associated molecular alterations. Our approach therefore not only brings new insights into the disease-specific and common molecular etiologies of AD and PD but also, in these realms, foster the discovery of more specific targets for functional and therapeutic exploration.