Multiscale PHATE identifies multimodal signatures of COVID-19

Single-cell spatial transcriptomics in cardiovascular development, disease, and medicine

Cardiovascular diseases (CVDs) impose a significant burden worldwide. Despite the elucidation of the etiology and underlying molecular mechanisms of CVDs by numerous studies and recent discovery of effective drugs, their morbidity, disability, and mortality are still high. Therefore, precise risk stratification and effective targeted therapies for CVDs are warranted. Recent improvements in single-cell RNA sequencing and spatial transcriptomics have improved our understanding of the mechanisms and cells involved in cardiovascular phylogeny and CVDs. Single-cell RNA sequencing can facilitate the study of the human heart at remarkably high resolution and cellular and molecular heterogeneity. However, this technique does not provide spatial information, which is essential for understanding homeostasis and disease. Spatial transcriptomics can elucidate intracellular interactions, transcription factor distribution, cell spatial localization, and molecular profiles of mRNA and identify cell populations causing the disease and their underlying mechanisms, including cell crosstalk. Herein, we introduce the main methods of RNA-seq and spatial transcriptomics analysis and highlight the latest advances in cardiovascular research. We conclude that single-cell RNA sequencing interprets disease progression in multiple dimensions, levels, perspectives, and dynamics by combining spatial and temporal characterization of the clinical phenome with multidisciplinary techniques such as spatial transcriptomics. This aligns with the dynamic evolution of CVDs (e.g., "angina-myocardial infarction-heart failure" in coronary artery disease). The study of pathways for disease onset and mechanisms (e.g., age, sex, comorbidities) in different patient subgroups should improve disease diagnosis and risk stratification. This can facilitate precise individualized treatment of CVDs.

Multi-omic analysis of Huntington's disease reveals a compensatory astrocyte state

The mechanisms underlying the selective regional vulnerability to neurodegeneration in Huntington's disease (HD) have not been fully defined. To explore the role of astrocytes in this phenomenon, we used single-nucleus and bulk RNAseq, lipidomics, HTT gene CAG repeat-length measurements, and multiplexed immunofluorescence on HD and control post-mortem brains. We identified genes that correlated with CAG repeat length, which were enriched in astrocyte genes, and lipidomic signatures that implicated poly-unsaturated fatty acids in sensitizing neurons to cell death. Because astrocytes play essential roles in lipid metabolism, we explored the heterogeneity of astrocytic states in both protoplasmic and fibrous-like (CD44+) astrocytes. Significantly, one protoplasmic astrocyte state showed high levels of metallothioneins and was correlated with the selective vulnerability of distinct striatal neuronal populations. When modeled in vitro, this state improved the viability of HD-patient-derived spiny projection neurons. Our findings uncover key roles of astrocytic states in protecting against neurodegeneration in HD.

Building and analyzing metacells in single-cell genomics data

The advent of high-throughput single-cell genomics technologies has fundamentally transformed biological sciences. Currently, millions of cells from complex biological tissues can be phenotypically profiled across multiple modalities. The scaling of computational methods to analyze and visualize such data is a constant challenge, and tools need to be regularly updated, if not redesigned, to cope with ever-growing numbers of cells. Over the last few years, metacells have been introduced to reduce the size and complexity of single-cell genomics data while preserving biologically relevant information and improving interpretability. Here, we review recent studies that capitalize on the concept of metacells-and the many variants in nomenclature that have been used. We further outline how and when metacells should (or should not) be used to analyze single-cell genomics data and what should be considered when analyzing such data at the metacell level. To facilitate the exploration of metacells, we provide a comprehensive tutorial on the construction and analysis of metacells from single-cell RNA-seq data ( https://github.com/GfellerLab/MetacellAnalysisTutorial ) as well as a fully integrated pipeline to rapidly build, visualize and evaluate metacells with different methods ( https://github.com/GfellerLab/MetacellAnalysisToolkit ).

Single Cell Profiling Distinguishes Leukemia-Selective Chemotypes

A central challenge in chemical biology is to distinguish molecular families in which small structural changes trigger large changes in cell biology. Such families might be ideal scaffolds for developing cell-selective chemical effectors - for example, molecules that activate DNA damage responses in malignant cells while sparing healthy cells. Across closely related structural variants, subtle structural changes have the potential to result in contrasting bioactivity patterns across different cell types. Here, we tested a 600-compound Diversity Set of screening molecules from the Boston University Center for Molecular Discovery (BU-CMD) in a novel phospho-flow assay that tracked fundamental cell biological processes, including DNA damage response, apoptosis, M-phase cell cycle, and protein synthesis in MV411 leukemia cells. Among the chemotypes screened, synthetic congeners of the rocaglate family were especially bioactive. In follow-up studies, 37 rocaglates were selected and deeply characterized using 12 million additional cellular measurements across MV411 leukemia cells and healthy peripheral blood mononuclear cells. Of the selected rocaglates, 92% displayed significant bioactivity in human cells, and 65% selectively induced DNA damage responses in leukemia and not healthy human blood cells. Furthermore, the signaling and cell-type selectivity were connected to structural features of rocaglate subfamilies. In particular, three rocaglates from the rocaglate pyrimidinone (RP) structural subclass were the only molecules that activated exceptional DNA damage responses in leukemia cells without activating a detectable DNA damage response in healthy cells. These results indicate that the RP subset should be extensively characterized for anticancer therapeutic potential as it relates to the DNA damage response. This single cell profiling approach advances a chemical biology platform to dissect how systematic variations in chemical structure can profoundly and differentially impact basic functions of healthy and diseased cells.

Archetype analysis and the PHATE algorithm as methods to describe and visualize pregnant women's levels of physical activity knowledge

The knowledge of physical activity (PA) recommended for pregnant women and practical application of it has positive impact on the outcome. Nevertheless, it is estimated that in high-income countries over 40% of pregnant women are insufficiently physically active. One of the reasons is insufficient knowledge pregnant women have about allowed effort during pregnancy and both recommended and not recommended physical activities. Description of knowledge about physical activity the women have and distinguishing patterns of their knowledge is becoming an increasingly important issue. A common approach to handle survey data that reflect knowledge involves clustering methods or Principal Component Analysis (PCA). Nevertheless, new procedures of data analysis are still being sought. Using survey data collected by the Institute of Mother and Child Archetypal analysis has been applied to detect levels of knowledge reflected by answers given in a questionnaire and to derive patterns of knowledge contained in the data. Next, PHATE (Potential of Heat-diffusion for Affinity-based Trajectory Embedding) algorithm has been used to visualize the results and to get a deeper insight into the data structure. The results were compared with picture derived from PCA. Three archetypes representing three patterns of knowledge have been distinguished and described. The presentation of complex data in a low dimension was obtained with help of PHATE. The formations revealed by PHATE have been successfully described in terms of knowledge levels reflected by the survey. Finally, comparison of PHATE with PCA has been shown. Archetype analysis combined with PHATE provides novel opportunities in examining nonlinear structure of survey data and allows for visualization that captures complex relations in the data. PHATE has made it possible to distinguish sets of objects that have common features but were captured neither by Archetypal analysis nor PCA. Moreover, for our data, PHATE provides an image of data structure which is more detailed than interpretation of PCA.

Starfysh integrates spatial transcriptomic and histologic data to reveal heterogeneous tumor-immune hubs

Spatially resolved gene expression profiling provides insight into tissue organization and cell-cell crosstalk; however, sequencing-based spatial transcriptomics (ST) lacks single-cell resolution. Current ST analysis methods require single-cell RNA sequencing data as a reference for rigorous interpretation of cell states, mostly do not use associated histology images and are not capable of inferring shared neighborhoods across multiple tissues. Here we present Starfysh, a computational toolbox using a deep generative model that incorporates archetypal analysis and any known cell type markers to characterize known or new tissue-specific cell states without a single-cell reference. Starfysh improves the characterization of spatial dynamics in complex tissues using histology images and enables the comparison of niches as spatial hubs across tissues. Integrative analysis of primary estrogen receptor (ER)-positive breast cancer, triple-negative breast cancer (TNBC) and metaplastic breast cancer (MBC) tissues led to the identification of spatial hubs with patient- and disease-specific cell type compositions and revealed metabolic reprogramming shaping immunosuppressive hubs in aggressive MBC.

Utilizing multimodal AI to improve genetic analyses of cardiovascular traits

Electronic health records, biobanks, and wearable biosensors contain multiple high-dimensional clinical data (HDCD) modalities (e.g., ECG, Photoplethysmography (PPG), and MRI) for each individual. Access to multimodal HDCD provides a unique opportunity for genetic studies of complex traits because different modalities relevant to a single physiological system (e.g., circulatory system) encode complementary and overlapping information. We propose a novel multimodal deep learning method, M-REGLE, for discovering genetic associations from a joint representation of multiple complementary HDCD modalities. We showcase the effectiveness of this model by applying it to several cardiovascular modalities. M-REGLE jointly learns a lower representation (i.e., latent factors) of multimodal HDCD using a convolutional variational autoencoder, performs genome wide association studies (GWAS) on each latent factor, then combines the results to study the genetics of the underlying system. To validate the advantages of M-REGLE and multimodal learning, we apply it to common cardiovascular modalities (PPG and ECG), and compare its results to unimodal learning methods in which representations are learned from each data modality separately, but the downstream genetic analyses are performed on the combined unimodal representations. M-REGLE identifies 19.3% more loci on the 12-lead ECG dataset, 13.0% more loci on the ECG lead I + PPG dataset, and its genetic risk score significantly outperforms the unimodal risk score at predicting cardiac phenotypes, such as atrial fibrillation (Afib), in multiple biobanks.

Toward the novel AI tasks in infection biology

Machine learning and artificial intelligence (AI) are becoming more common in infection biology laboratories around the world. Yet, as they gain traction in research, novel frontiers arise. Novel artificial intelligence algorithms are capable of addressing advanced tasks like image generation and question answering. However, similar algorithms can prove useful in addressing advanced questions in infection biology like prediction of host-pathogen interactions or inferring virus protein conformations. Addressing such tasks requires large annotated data sets, which are often scarce in biomedical research. In this review, I bring together several successful examples where such tasks were addressed. I underline the importance of formulating novel AI tasks in infection biology accompanied by freely available benchmark data sets to address these tasks. Furthermore, I discuss the current state of the field and potential future trends. I argue that one such trend involves AI tools becoming more versatile.

Detection of PatIent-Level distances from single cell genomics and pathomics data with Optimal Transport (PILOT)

Although clinical applications represent the next challenge in single-cell genomics and digital pathology, we still lack computational methods to analyze single-cell or pathomics data to find sample-level trajectories or clusters associated with diseases. This remains challenging as single-cell/pathomics data are multi-scale, i.e., a sample is represented by clusters of cells/structures, and samples cannot be easily compared with each other. Here we propose PatIent Level analysis with Optimal Transport (PILOT). PILOT uses optimal transport to compute the Wasserstein distance between two individual single-cell samples. This allows us to perform unsupervised analysis at the sample level and uncover trajectories or cellular clusters associated with disease progression. We evaluate PILOT and competing approaches in single-cell genomics or pathomics studies involving various human diseases with up to 600 samples/patients and millions of cells or tissue structures. Our results demonstrate that PILOT detects disease-associated samples from large and complex single-cell or pathomics data. Moreover, PILOT provides a statistical approach to find changes in cell populations, gene expression, and tissue structures related to the trajectories or clusters supporting interpretation of predictions.

Gaining Biological Insights through Supervised Data Visualization

Dimensionality reduction-based data visualization is pivotal in comprehending complex biological data. The most common methods, such as PHATE, t-SNE, and UMAP, are unsupervised and therefore reflect the dominant structure in the data, which may be independent of expert-provided labels. Here we introduce a supervised data visualization method called RF-PHATE, which integrates expert knowledge for further exploration of the data. RF-PHATE leverages random forests to capture intricate featurelabel relationships. Extracting information from the forest, RF-PHATE generates low-dimensional visualizations that highlight relevant data relationships while disregarding extraneous features. This approach scales to large datasets and applies to classification and regression. We illustrate RF-PHATE's prowess through three case studies. In a multiple sclerosis study using longitudinal clinical and imaging data, RF-PHATE unveils a sub-group of patients with non-benign relapsingremitting Multiple Sclerosis, demonstrating its aptitude for time-series data. In the context of Raman spectral data, RF-PHATE effectively showcases the impact of antioxidants on diesel exhaust-exposed lung cells, highlighting its proficiency in noisy environments. Furthermore, RF-PHATE aligns established geometric structures with COVID-19 patient outcomes, enriching interpretability in a hierarchical manner. RF-PHATE bridges expert insights and visualizations, promising knowledge generation. Its adaptability, scalability, and noise tolerance underscore its potential for widespread adoption.

Benchmarking differential abundance methods for finding condition-specific prototypical cells in multi-sample single-cell datasets

BACKGROUND

To analyze the large volume of data generated by single-cell technologies and to identify cellular correlates of particular clinical or experimental outcomes, differential abundance analyses are often applied. These algorithms identify subgroups of cells whose abundances change significantly in response to disease progression, or to an experimental perturbation. Despite the effectiveness of differential abundance analyses in identifying critical cell-states, there is currently no systematic benchmarking study to compare their applicability, usefulness, and accuracy in practice across single-cell modalities.

RESULTS

Here, we perform a comprehensive benchmarking study to objectively evaluate and compare the benefits and potential downsides of current state-of-the-art differential abundance testing methods. We benchmarked six single-cell testing methods on several practical tasks, using both synthetic and real single-cell datasets. The tasks evaluated include effectiveness in identifying true differentially abundant subpopulations, accuracy in the adequate handling of batch effects, runtime efficiency, and hyperparameter usability and robustness. Based on various evaluation results, this paper gives dataset-specific suggestions for the practical use of differential abundance testing approaches.

CONCLUSIONS

Based on our benchmarking study, we provide a set of recommendations for the optimal usage of single-cell DA testing methods in practice, particularly with respect to factors such as the presence of technical noise (for example batch effects), dataset size, and hyperparameter sensitivity.

Deep generative model deciphers derailed trajectories in acute myeloid leukemia

Single-cell genomics has the potential to map cell states and their dynamics in an unbiased way in response to perturbations like disease. However, elucidating the cell-state transitions from healthy to disease requires analyzing data from perturbed samples jointly with unperturbed reference samples. Existing methods for integrating and jointly visualizing single-cell datasets from distinct contexts tend to remove key biological differences or do not correctly harmonize shared mechanisms. We present Decipher, a model that combines variational autoencoders with deep exponential families to reconstruct derailed trajectories ( https://github.com/azizilab/decipher ). Decipher jointly represents normal and perturbed single-cell RNA-seq datasets, revealing shared and disrupted dynamics. It further introduces a novel approach to visualize data, without the need for methods such as UMAP or TSNE. We demonstrate Decipher on data from acute myeloid leukemia patient bone marrow specimens, showing that it successfully characterizes the divergence from normal hematopoiesis and identifies transcriptional programs that become disrupted in each patient when they acquire NPM1 driver mutations.

Multi-OMIC analysis of Huntington disease reveals a neuroprotective astrocyte state

Huntington disease (HD) is an incurable neurodegenerative disease characterized by neuronal loss and astrogliosis. One hallmark of HD is the selective neuronal vulnerability of striatal medium spiny neurons. To date, the underlying mechanisms of this selective vulnerability have not been fully defined. Here, we employed a multi-omic approach including single nucleus RNAseq (snRNAseq), bulk RNAseq, lipidomics, HTT gene CAG repeat length measurements, and multiplexed immunofluorescence on post-mortem brain tissue from multiple brain regions of HD and control donors. We defined a signature of genes that is driven by CAG repeat length and found it enriched in astrocytic and microglial genes. Moreover, weighted gene correlation network analysis showed loss of connectivity of astrocytic and microglial modules in HD and identified modules that correlated with CAG-repeat length which further implicated inflammatory pathways and metabolism. We performed lipidomic analysis of HD and control brains and identified several lipid species that correlate with HD grade, including ceramides and very long chain fatty acids. Integration of lipidomics and bulk transcriptomics identified a consensus gene signature that correlates with HD grade and HD lipidomic abnormalities and implicated the unfolded protein response pathway. Because astrocytes are critical for brain lipid metabolism and play important roles in regulating inflammation, we analyzed our snRNAseq dataset with an emphasis on astrocyte pathology. We found two main astrocyte types that spanned multiple brain regions; these types correspond to protoplasmic astrocytes, and fibrous-like - CD44-positive, astrocytes. HD pathology was differentially associated with these cell types in a region-specific manner. One protoplasmic astrocyte cluster showed high expression of metallothionein genes, the depletion of this cluster positively correlated with the depletion of vulnerable medium spiny neurons in the caudate nucleus. We confirmed that metallothioneins were increased in cingulate HD astrocytes but were unchanged or even decreased in caudate astrocytes. We combined existing genome-wide association studies (GWAS) with a GWA study conducted on HD patients from the original Venezuelan cohort and identified a single-nucleotide polymorphism in the metallothionein gene locus associated with delayed age of onset. Functional studies found that metallothionein overexpressing astrocytes are better able to buffer glutamate and were neuroprotective of patient-derived directly reprogrammed HD MSNs as well as against rotenone-induced neuronal death in vitro. Finally, we found that metallothionein-overexpressing astrocytes increased the phagocytic activity of microglia in vitro and increased the expression of genes involved in fatty acid binding. Together, we identified an astrocytic phenotype that is regionally-enriched in less vulnerable brain regions that can be leveraged to protect neurons in HD.

A unique cytotoxic CD4+ T cell-signature defines critical COVID-19

Objectives

SARS-CoV-2 infection causes a spectrum of clinical disease presentation, ranging from asymptomatic to fatal. While neutralising antibody (NAb) responses correlate with protection against symptomatic and severe infection, the contribution of the T-cell response to disease resolution or progression is still unclear. As newly emerging variants of concern have the capacity to partially escape NAb responses, defining the contribution of individual T-cell subsets to disease outcome is imperative to inform the development of next-generation COVID-19 vaccines.

Methods

Immunophenotyping of T-cell responses in unvaccinated individuals was performed, representing the full spectrum of COVID-19 clinical presentation. Computational and manual analyses were used to identify T-cell populations associated with distinct disease states.

Results

Critical SARS-CoV-2 infection was characterised by an increase in activated and cytotoxic CD4+ lymphocytes (CTL). These CD4+ CTLs were largely absent in asymptomatic to severe disease states. In contrast, non-critical COVID-19 was associated with high frequencies of naïve T cells and lack of activation marker expression.

Conclusion

Highly activated and cytotoxic CD4+ T-cell responses may contribute to cell-mediated host tissue damage and progression of COVID-19. Induction of these potentially detrimental T-cell responses should be considered when developing and implementing effective COVID-19 control strategies.

Multiscale geometric and topological analyses for characterizing and predicting immune responses from single cell data

Single cell genomics has revolutionized our ability to map immune heterogeneity and responses. With the influx of large-scale data sets from diverse modalities, the resolution achieved has supported the long-held notion that immune cells are naturally organized into hierarchical relationships, characterized at multiple levels. Such a multigranular structure corresponds to key geometric and topological features. Given that differences between an effective and ineffective immunological response may not be found at one level, there is vested interest in characterizing and predicting outcomes from such features. In this review, we highlight single cell methods and principles for learning geometric and topological properties of data at multiple scales, discussing their contributions to immunology. Ultimately, multiscale approaches go beyond classical clustering, revealing a more comprehensive picture of cellular heterogeneity.

HSV-2 triggers upregulation of MALAT1 in CD4+ T cells and promotes HIV latency reversal

Herpes simplex virus type 2 (HSV-2) coinfection is associated with increased HIV-1 viral loads and expanded tissue reservoirs, but the mechanisms are not well defined. HSV-2 recurrences result in an influx of activated CD4+ T cells to sites of viral replication and an increase in activated CD4+ T cells in peripheral blood. We hypothesized that HSV-2 induces changes in these cells that facilitate HIV-1 reactivation and replication and tested this hypothesis in human CD4+ T cells and 2D10 cells, a model of HIV-1 latency. HSV-2 promoted latency reversal in HSV-2-infected and bystander 2D10 cells. Bulk and single-cell RNA-Seq studies of activated primary human CD4+ T cells identified decreased expression of HIV-1 restriction factors and increased expression of transcripts including MALAT1 that could drive HIV replication in both the HSV-2-infected and bystander cells. Transfection of 2D10 cells with VP16, an HSV-2 protein that regulates transcription, significantly upregulated MALAT1 expression, decreased trimethylation of lysine 27 on histone H3 protein, and triggered HIV latency reversal. Knockout of MALAT1 from 2D10 cells abrogated the response to VP16 and reduced the response to HSV-2 infection. These results demonstrate that HSV-2 contributes to HIV-1 reactivation through diverse mechanisms, including upregulation of MALAT1 to release epigenetic silencing.

linc-mipep and linc-wrb encode micropeptides that regulate chromatin accessibility in vertebrate-specific neural cells

Thousands of long intergenic non-coding RNAs (lincRNAs) are transcribed throughout the vertebrate genome. A subset of lincRNAs enriched in developing brains have recently been found to contain cryptic open-reading frames and are speculated to encode micropeptides. However, systematic identification and functional assessment of these transcripts have been hindered by technical challenges caused by their small size. Here, we show that two putative lincRNAs (linc-mipep, also called lnc-rps25, and linc-wrb) encode micropeptides with homology to the vertebrate-specific chromatin architectural protein, Hmgn1, and demonstrate that they are required for development of vertebrate-specific brain cell types. Specifically, we show that NMDA receptor-mediated pathways are dysregulated in zebrafish lacking these micropeptides and that their loss preferentially alters the gene regulatory networks that establish cerebellar cells and oligodendrocytes - evolutionarily newer cell types that develop postnatally in humans. These findings reveal a key missing link in the evolution of vertebrate brain cell development and illustrate a genetic basis for how some neural cell types are more susceptible to chromatin disruptions, with implications for neurodevelopmental disorders and disease.

Interferon lambda as a potential treatment for COVID-19

INTRODUCTION

Pegylated interferon lambda substantially reduced the risk of COVID-19-related hospitalizations or emergency room visits in a recent phase 3, multi-center, randomized, double-blind, placebo-controlled study of high-risk, non-hospitalized adult patients with SARS-CoV-2 infection compared to treatment with placebo.

AREAS COVERED

Interferons are a family of signaling molecules produced as part of the innate immune response to viral infections. The administration of exogenous interferon may limit disease progression in patients with COVID-19.

EXPERT OPINION

Interferons have been used to treat viral infections, including hepatitis B and hepatitis C, and malignancies such as non-Hodgkin's lymphoma, as well as the autoimmune condition multiple sclerosis. This manuscript examines what is known about the role of interferon lambda in the treatment of COVID-19, including potential limitations, and explores how this approach may be used in the future.

Mother-to-infant microbiota transmission and infant microbiota development across multiple body sites

Early-life microbiota seeding and subsequent development is crucial to future health. Cesarean-section (CS) birth, as opposed to vaginal delivery, affects early mother-to-infant transmission of microbes. Here, we assess mother-to-infant microbiota seeding and early-life microbiota development across six maternal and four infant niches over the first 30 days of life in 120 mother-infant pairs. Across all infants, we estimate that on average 58.5% of the infant microbiota composition can be attributed to any of the maternal source communities. All maternal source communities seed multiple infant niches. We identify shared and niche-specific host/environmental factors shaping the infant microbiota. In CS-born infants, we report reduced seeding of infant fecal microbiota by maternal fecal microbes, whereas colonization with breastmilk microbiota is increased when compared with vaginally born infants. Therefore, our data suggest auxiliary routes of mother-to-infant microbial seeding, which may compensate for one another, ensuring that essential microbes/microbial functions are transferred irrespective of disrupted transmission routes.

Artificial intelligence in clinical multiparameter flow cytometry and mass cytometry-key tools and progress

There are many research studies and emerging tools using artificial intelligence (AI) and machine learning to augment flow and mass cytometry workflows. Emerging AI tools can quickly identify common cell populations with continuous improvement of accuracy, uncover patterns in high-dimensional cytometric data that are undetectable by human analysis, facilitate the discovery of cell subpopulations, perform semi-automated immune cell profiling, and demonstrate potential to automate aspects of clinical multiparameter flow cytometric (MFC) diagnostic workflow. Utilizing AI in the analysis of cytometry samples can reduce subjective variability and assist in breakthroughs in understanding diseases. Here we review the diverse types of AI that are being applied to clinical cytometry data and how AI is driving advances in data analysis to improve diagnostic sensitivity and accuracy. We review supervised and unsupervised clustering algorithms for cell population identification, various dimensionality reduction techniques, and their utilities in visualization and machine learning pipelines, and supervised learning approaches for classifying entire cytometry samples.Understanding the AI landscape will enable pathologists to better utilize open source and commercially available tools, plan exploratory research projects to characterize diseases, and work with machine learning and data scientists to implement clinical data analysis pipelines.

Multi-view manifold learning of human brain-state trajectories

The complexity of the human brain gives the illusion that brain activity is intrinsically high-dimensional. Nonlinear dimensionality-reduction methods such as uniform manifold approximation and t-distributed stochastic neighbor embedding have been used for high-throughput biomedical data. However, they have not been used extensively for brain activity data such as those from functional magnetic resonance imaging (fMRI), primarily due to their inability to maintain dynamic structure. Here we introduce a nonlinear manifold learning method for time-series data-including those from fMRI-called temporal potential of heat-diffusion for affinity-based transition embedding (T-PHATE). In addition to recovering a low-dimensional intrinsic manifold geometry from time-series data, T-PHATE exploits the data's autocorrelative structure to faithfully denoise and unveil dynamic trajectories. We empirically validate T-PHATE on three fMRI datasets, showing that it greatly improves data visualization, classification, and segmentation of the data relative to several other state-of-the-art dimensionality-reduction benchmarks. These improvements suggest many potential applications of T-PHATE to other high-dimensional datasets of temporally diffuse processes.

The race to understand immunopathology in COVID-19: Perspectives on the impact of quantitative approaches to understand within-host interactions

The COVID-19 pandemic has revealed the need for the increased integration of modelling and data analysis to public health, experimental, and clinical studies. Throughout the first two years of the pandemic, there has been a concerted effort to improve our understanding of the within-host immune response to the SARS-CoV-2 virus to provide better predictions of COVID-19 severity, treatment and vaccine development questions, and insights into viral evolution and the impacts of variants on immunopathology. Here we provide perspectives on what has been accomplished using quantitative methods, including predictive modelling, population genetics, machine learning, and dimensionality reduction techniques, in the first 26 months of the COVID-19 pandemic approaches, and where we go from here to improve our responses to this and future pandemics.

Benchmarking differential abundance methods for finding condition-specific prototypical cells in multi-sample single-cell datasets

Modern single-cell data analysis relies on statistical testing (e.g. differential expression testing) to identify genes or proteins that are up-or down-regulated in relation to cell-types or clinical outcomes. However, existing algorithms for such statistical testing are often limited by technical noise and cellular heterogeneity, which lead to false-positive results. To constrain the analysis to a compact and phenotype-related cell population, differential abundance (DA) testing methods were employed to identify subgroups of cells whose abundance changed significantly in response to disease progression, or experimental perturbation. Despite the effectiveness of DA testing algorithms of identifying critical cell-states, there are no systematic benchmarking or comparative studies to compare their usages in practice. Herein, we performed the first comprehensive benchmarking study to objectively evaluate and compare the benefits and potential downsides of current state-of-the-art DA testing methods. We benchmarked six DA testing methods on several practical tasks, using both synthetic and real single-cell datasets. The task evaluated include, recognizing true DA subpopulations, appropriate handing of batch effects, runtime efficiency, and hyperparameter usability and robustness. Based on various evaluation results, this paper gives dataset-specific suggestions for the usage of DA testing methods.

Innate and Adaptive Immunity during SARS-CoV-2 Infection: Biomolecular Cellular Markers and Mechanisms

The coronavirus 2019 (COVID-19) pandemic was caused by a positive sense single-stranded RNA (ssRNA) severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). However, other human coronaviruses (hCoVs) exist. Historical pandemics include smallpox and influenza, with efficacious therapeutics utilized to reduce overall disease burden through effectively targeting a competent host immune system response. The immune system is composed of primary/secondary lymphoid structures with initially eight types of immune cell types, and many other subtypes, traversing cell membranes utilizing cell signaling cascades that contribute towards clearance of pathogenic proteins. Other proteins discussed include cluster of differentiation (CD) markers, major histocompatibility complexes (MHC), pleiotropic interleukins (IL), and chemokines (CXC). The historical concepts of host immunity are the innate and adaptive immune systems. The adaptive immune system is represented by T cells, B cells, and antibodies. The innate immune system is represented by macrophages, neutrophils, dendritic cells, and the complement system. Other viruses can affect and regulate cell cycle progression for example, in cancers that include human papillomavirus (HPV: cervical carcinoma), Epstein-Barr virus (EBV: lymphoma), Hepatitis B and C (HB/HC: hepatocellular carcinoma) and human T cell Leukemia Virus-1 (T cell leukemia). Bacterial infections also increase the risk of developing cancer (e.g., Helicobacter pylori). Viral and bacterial factors can cause both morbidity and mortality alongside being transmitted within clinical and community settings through affecting a host immune response. Therefore, it is appropriate to contextualize advances in single cell sequencing in conjunction with other laboratory techniques allowing insights into immune cell characterization. These developments offer improved clarity and understanding that overlap with autoimmune conditions that could be affected by innate B cells (B1+ or marginal zone cells) or adaptive T cell responses to SARS-CoV-2 infection and other pathologies. Thus, this review starts with an introduction into host respiratory infection before examining invaluable cellular messenger proteins and then individual immune cell markers.

NK- and T-cell granzyme B and K expression correlates with age, CMV infection and influenza vaccine-induced antibody titres in older adults

Granzymes are a family of serine-proteases that act as critical mediators in the cytolytic and immunomodulatory activities of immune cells such as CD8⁺ T-cells and natural killer (NK) cells. Previous work indicates that both granzyme B (GZB) and K (GZK) are increased with age in CD8⁺ T-cells, and in the case of GZB, contribute to dysfunctional immune processes observed in older adults. Here, we sought to determine how GZB and GZK expression in NK-cells, and CD4⁺, CD8⁺, and gamma-delta T-cells, quantified in terms of positive cell frequency and mean fluorescence intensity (MFI), differed with age, age-related health-traits and the antibody response to high-dose influenza vaccine. We found that the frequency and MFI of GZB-expressing NK-cells, and CD8⁺ and Vδ1+ T-cells, and GZK-expressing CD8⁺ T-cells was significantly higher in older (66-97 years old; n = 75) vs. younger (24-37 years old; n = 10) adults by up to 5-fold. There were no significant associations of GZB/GZK expression with sex, frailty or plasma levels of TNF or IL-6 in older adults, but those who were seropositive for cytomegalovirus (CMV) exhibited significantly higher frequencies of GZB+ NK-cells, and CD4⁺, CD8⁺ and Vδ1+ T-cells, and GZK+ CD8⁺ T-cells (Cohen's d = .5-1.5). Pre-vaccination frequencies of GZB+ NK-cells were positively correlated with vaccine antibody responses against A/H3N2 (d = .17), while the frequencies of GZK+ NK and CD8⁺ T-cells were inversely associated with A/H1N1 (d = -0.18 to -0.20). Interestingly, GZK+ NK-cell frequency was inversely correlated with pre-vaccination A/H1N1 antibody titres, as well as those measured over the previous 4 years, further supporting a role for this subset in influencing vaccine antibody-responses. These findings further our understanding of how granzyme expression in different lymphoid cell-types may change with age, while suggesting that they influence vaccine responsiveness in older adults.

Dissecting CD8+ T cell pathology of severe SARS-CoV-2 infection by single-cell immunoprofiling

Introduction

SARS-CoV-2 infection results in varying disease severity, ranging from asymptomatic infection to severe illness. A detailed understanding of the immune response to SARS-CoV-2 is critical to unravel the causative factors underlying differences in disease severity and to develop optimal vaccines against new SARS-CoV-2 variants.

Methods

We combined single-cell RNA and T cell receptor sequencing with CITE-seq antibodies to characterize the CD8⁺ T cell response to SARS-CoV-2 infection at high resolution and compared responses between mild and severe COVID-19.

Results

We observed increased CD8⁺ T cell exhaustion in severe SARS-CoV-2 infection and identified a population of NK-like, terminally differentiated CD8⁺ effector T cells characterized by expression of FCGR3A (encoding CD16). Further characterization of NK-like CD8⁺ T cells revealed heterogeneity among CD16⁺ NK-like CD8⁺ T cells and profound differences in cytotoxicity, exhaustion, and NK-like differentiation between mild and severe disease conditions.

Discussion

We propose a model in which differences in the surrounding inflammatory milieu lead to crucial differences in NK-like differentiation of CD8⁺ effector T cells, ultimately resulting in the appearance of NK-like CD8⁺ T cell populations of different functionality and pathogenicity. Our in-depth characterization of the CD8⁺ T cell-mediated response to SARS-CoV-2 infection provides a basis for further investigation of the importance of NK-like CD8⁺ T cells in COVID-19 severity.

Domestic dog lineages reveal genetic drivers of behavioral diversification

Selective breeding of domestic dogs has generated diverse breeds often optimized for performing specialized tasks. Despite the heritability of breed-typical behavioral traits, identification of causal loci has proven challenging due to the complexity of canine population structure. We overcome longstanding difficulties in identifying genetic drivers of canine behavior by developing a framework for understanding relationships between breeds and the behaviors that define them, utilizing genetic data for over 4,000 domestic, semi-feral, and wild canids and behavioral survey data for over 46,000 dogs. We identify ten major canine genetic lineages and their behavioral correlates and show that breed diversification is predominantly driven by non-coding regulatory variation. We determine that lineage-associated genes converge in neurodevelopmental co-expression networks, identifying a sheepdog-associated enrichment for interrelated axon guidance functions. This work presents a scaffold for canine diversification that positions the domestic dog as an unparalleled system for revealing the genetic origins of behavioral diversity.

Immune cells and their inflammatory mediators modify β cells and cause checkpoint inhibitor-induced diabetes

Checkpoint inhibitors (CPIs) targeting programmed death 1 (PD-1)/programmed death ligand 1 (PD-L1) and cytotoxic T lymphocyte antigen 4 (CTLA-4) have revolutionized cancer treatment but can trigger autoimmune complications, including CPI-induced diabetes mellitus (CPI-DM), which occurs preferentially with PD-1 blockade. We found evidence of pancreatic inflammation in patients with CPI-DM with shrinkage of pancreases, increased pancreatic enzymes, and in a case from a patient who died with CPI-DM, peri-islet lymphocytic infiltration. In the NOD mouse model, anti-PD-L1 but not anti-CTLA-4 induced diabetes rapidly. RNA sequencing revealed that cytolytic IFN-γ+CD8+ T cells infiltrated islets with anti-PD-L1. Changes in β cells were predominantly driven by IFN-γ and TNF-α and included induction of a potentially novel β cell population with transcriptional changes suggesting dedifferentiation. IFN-γ increased checkpoint ligand expression and activated apoptosis pathways in human β cells in vitro. Treatment with anti-IFN-γ and anti-TNF-α prevented CPI-DM in anti-PD-L1-treated NOD mice. CPIs targeting the PD-1/PD-L1 pathway resulted in transcriptional changes in β cells and immune infiltrates that may lead to the development of diabetes. Inhibition of inflammatory cytokines can prevent CPI-DM, suggesting a strategy for clinical application to prevent this complication.

Analyzing network diversity of cell–cell interactions in COVID-19 using single-cell transcriptomics

Cell–cell interactions (CCI) play significant roles in manipulating biological functions of cells. Analyzing the differences in CCI between healthy and diseased conditions of a biological system yields greater insight than analyzing either conditions alone. There has been a recent and rapid growth of methods to infer CCI from single-cell RNA-sequencing (scRNA-seq), revealing complex CCI networks at a previously inaccessible scale. However, the majority of current CCI analyses from scRNA-seq data focus on direct comparisons between individual CCI networks of individual samples from patients, rather than “group-level” comparisons between sample groups of patients comprising different conditions. To illustrate new biological features among different disease statuses, we investigated the diversity of key network features on groups of CCI networks, as defined by different disease statuses. We considered three levels of network features: node level, as defined by cell type; node-to-node level; and network level. By applying these analysis to a large-scale single-cell RNA-sequencing dataset of coronavirus disease 2019 (COVID-19), we observe biologically meaningful patterns aligned with the progression and subsequent convalescence of COVID-19.

Mapping Phenotypic Plasticity upon the Cancer Cell State Landscape Using Manifold Learning

ABSTRACT

Phenotypic plasticity describes the ability of cancer cells to undergo dynamic, nongenetic cell state changes that amplify cancer heterogeneity to promote metastasis and therapy evasion. Thus, cancer cells occupy a continuous spectrum of phenotypic states connected by trajectories defining dynamic transitions upon a cancer cell state landscape. With technologies proliferating to systematically record molecular mechanisms at single-cell resolution, we illuminate manifold learning techniques as emerging computational tools to effectively model cell state dynamics in a way that mimics our understanding of the cell state landscape. We anticipate that "state-gating" therapies targeting phenotypic plasticity will limit cancer heterogeneity, metastasis, and therapy resistance.

SIGNIFICANCE

Nongenetic mechanisms underlying phenotypic plasticity have emerged as significant drivers of tumor heterogeneity, metastasis, and therapy resistance. Herein, we discuss new experimental and computational techniques to define phenotypic plasticity as a scaffold to guide accelerated progress in uncovering new vulnerabilities for therapeutic exploitation.

Heterogeneity in IgG-CD16 signaling in infectious disease outcomes

In this review, we discuss how IgG antibodies can modulate inflammatory signaling during viral infections with a focus on CD16a-mediated functions. We describe the structural heterogeneity of IgG antibody ligands, including subclass and glycosylation that impact binding by and downstream activity of CD16a, as well as the heterogeneity of CD16a itself, including allele and expression density. While inflammation is a mechanism required for immune homeostasis and resolution of acute infections, we focus here on two infectious diseases that are driven by pathogenic inflammatory responses during infection. Specifically, we review and discuss the evolving body of literature showing that afucosylated IgG immune complex signaling through CD16a contributes to the overwhelming inflammatory response that is central to the pathogenesis of severe forms of dengue disease and coronavirus disease 2019 (COVID-19).

A diminished immune response underlies age-related SARS-CoV-2 pathologies

Morbidity and mortality in response to SARS-CoV-2 infection are significantly elevated in people of advanced age. To understand the underlying biology of this phenotype, we utilize the golden hamster model to compare how the innate and adaptive immune responses to SARS-CoV-2 infection differed between younger and older animals. We find that while both hamster cohorts showed similar virus kinetics in the lungs, the host response in older animals was dampened, with diminished tissue repair in the respiratory tract post-infection. Characterization of the adaptive immune response also revealed age-related differences, including fewer germinal center B cells in older hamsters, resulting in reduced potency of neutralizing antibodies. Moreover, older animals demonstrate elevated suppressor T cells and neutrophils in the respiratory tract, correlating with an increase in TGF-β and IL-17 induction. Together, these data support that diminished immunity is one of the underlying causes of age-related morbidity.