Sustained IFN signaling is associated with delayed development of SARS-CoV-2-specific immunity.

Plasma RNAemia, delayed antibody responses and inflammation predict COVID-19 outcomes, but the mechanisms underlying these immunovirological patterns are poorly understood. We profile 782 longitudinal plasma samples from 318 hospitalized patients with COVID-19. Integrated analysis using k-means reveals four patient clusters in a discovery cohort: mechanically ventilated critically-ill cases are subdivided into good prognosis and high-fatality clusters (reproduced in a validation cohort), while non-critical survivors segregate into high and low early antibody responders. Only the high-fatality cluster is enriched for transcriptomic signatures associated with COVID-19 severity, and each cluster has distinct RBD-specific antibody elicitation kinetics. Both critical and non-critical clusters with delayed antibody responses exhibit sustained IFN signatures, which negatively correlate with contemporaneous RBD-specific IgG levels and absolute SARS-CoV-2-specific B and CD4+ T cell frequencies. These data suggest that the "Interferon paradox" previously described in murine LCMV models is operative in COVID-19, with excessive IFN signaling delaying development of adaptive virus-specific immunity.

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.

Optimizing health and nutrition status of migrant construction workers consuming multiple micronutrient fortified rice in Singapore.

INTRODUCTION: A well-nourished workforce is instrumental in eradicating hunger, alleviating poverty, and spurring economic growth. A fifth of the total workforce in high-income countries are migrant workers. Despite the accessibility of nutritious foods in high-income countries, migrant workers often rely on nutrient-poor diets largely consisting of empty calories, which in turn leads to vitamin and mineral deficiency, also called hidden hunger, and resultant productivity loss. Here, we study the magnitude of hidden hunger in male migrant construction workers in Singapore and investigate the impact of consuming fortified rice for 6 consecutive months on the nutrition and health status of these workers. METHODS: 140 male migrant workers aged 20-51 years of either Bangladeshi or Indian ethnicity from a single dormitory in Singapore volunteered to participate in the study. In total, 133 blood samples were taken at the start of the study and were used to assess vitamin B12, hemoglobin, ferritin, folate, and zinc levels; a sub-sample underwent for homocysteine testing. Anthropometric measurements and vital signs, such as blood pressure, were recorded before and after the intervention. RESULTS: The results show that vitamin and mineral deficiency was present, especially folate (59% of workers deficient) and vitamin B12 (7% deficient, 31% marginally deficient). The consumption of fortified rice significantly improved the vitamin, iron and zinc level in the workers and significantly reduced the systolic blood pressure amongst the Bangladeshi migrant workers, specifically. CONCLUSION: Our study demonstrates that fortified rice may have a positive impact on male migrant construction worker health and nutrition status at the workplace.

Geometry- and Accuracy-Preserving Random Forest Proximities.

Random forests are considered one of the best out-of-the-box classification and regression algorithms due to their high level of predictive performance with relatively little tuning. Pairwise proximities can be computed from a trained random forest and measure the similarity between data points relative to the supervised task. Random forest proximities have been used in many applications including the identification of variable importance, data imputation, outlier detection, and data visualization. However, existing definitions of random forest proximities do not accurately reflect the data geometry learned by the random forest. In this paper, we introduce a novel definition of random forest proximities called Random Forest-Geometry- and Accuracy-Preserving proximities (RF-GAP). We prove that the proximity-weighted sum (regression) or majority vote (classification) using RF-GAP exactly matches the out-of-bag random forest prediction, thus capturing the data geometry learned by the random forest. We empirically show that this improved geometric representation outperforms traditional random forest proximities in tasks such as data imputation and provides outlier detection and visualization results consistent with the learned data geometry.

Geometry Regularized Autoencoders.

A fundamental task in data exploration is to extract low dimensional representations that capture intrinsic geometry in data, especially for faithfully visualizing data in two or three dimensions. Common approaches use kernel methods for manifold learning. However, these methods typically only provide an embedding of the input data and cannot extend naturally to new data points. Autoencoders have also become popular for representation learning. While they naturally compute feature extractors that are extendable to new data and invertible (i.e., reconstructing original features from latent representation), they often fail at representing the intrinsic data geometry compared to kernel-based manifold learning. We present a new method for integrating both approaches by incorporating a geometric regularization term in the bottleneck of the autoencoder. This regularization encourages the learned latent representation to follow the intrinsic data geometry, similar to manifold learning algorithms, while still enabling faithful extension to new data and preserving invertibility. We compare our approach to autoencoder models for manifold learning to provide qualitative and quantitative evidence of our advantages in preserving intrinsic structure, out of sample extension, and reconstruction. Our method is easily implemented for big-data applications, whereas other methods are limited in this regard.

Structural and developmental principles of neuropil assembly in C. elegans.

Neuropil is a fundamental form of tissue organization within the brain1, in which densely packed neurons synaptically interconnect into precise circuit architecture2,3. However, the structural and developmental principles that govern this nanoscale precision remain largely unknown4,5. Here we use an iterative data coarse-graining algorithm termed 'diffusion condensation'6 to identify nested circuit structures within the Caenorhabditis elegans neuropil, which is known as the nerve ring. We show that the nerve ring neuropil is largely organized into four strata that are composed of related behavioural circuits. The stratified architecture of the neuropil is a geometrical representation of the functional segregation of sensory information and motor outputs, with specific sensory organs and muscle quadrants mapping onto particular neuropil strata. We identify groups of neurons with unique morphologies that integrate information across strata and that create neural structures that cage the strata within the nerve ring. We use high resolution light-sheet microscopy7,8 coupled with lineage-tracing and cell-tracking algorithms9,10 to resolve the developmental sequence and reveal principles of cell position, migration and outgrowth that guide stratified neuropil organization. Our results uncover conserved structural design principles that underlie the architecture and function of the nerve ring neuropil, and reveal a temporal progression of outgrowth-based on pioneer neurons-that guides the hierarchical development of the layered neuropil. Our findings provide a systematic blueprint for using structural and developmental approaches to understand neuropil organization within the brain.

Imaging of PD-L1 in single cancer cells by SERS-based hyperspectral analysis.

We developed a hyperspectral imaging tool based on surface-enhanced Raman spectroscopy (SERS) probes to determine the expression level and visualize the distribution of PD-L1 in individual cells. Electron-microscopic analysis of PD-L1 antibody - gold nanorod conjugates demonstrated binding the cell surface and internalization into endosomal vesicles. Stimulation of cells with IFN-γ or metformin was used to confirm the ability of SERS probes to report treatment-induced changes. The multivariate curve resolution-alternating least squares (MCR-ALS) analysis of spectra provided a greater signal-noise ratio than single peak mapping. However, single peak mapping allowed a systematic subtraction of background and the removal of non-specific binding and endocytic SERS signals. The mean or maximum peak height in the cell or the mean peak height in the area of specific PD-L1 positive pixels was used to estimate the PD-L1 expression levels in single cells. The PD-L1 levels were significantly up-regulated by IFN-γ and inhibited by metformin in human lung cancer cells from the A549 cell line. In conclusion, the method of analyzing hyperspectral SERS imaging data together with systematic and comprehensive removal of non-specific signals allows SERS imaging to be a quantitative tool in the detection of the cancer biomarker, PD-L1.

Label-free discrimination and quantitative analysis of oxidative stress induced cytotoxicity and potential protection of antioxidants using Raman micro-spectroscopy and machine learning.

Diesel exhaust particles (DEPs) are major constituents of air pollution and associated with numerous oxidative stress-induced human diseases. In vitro toxicity studies are useful for developing a better understanding of species-specific in vivo conditions. Conventional in vitro assessments based on oxidative biomarkers are destructive and inefficient. In this study, Raman spectroscopy, as a non-invasive imaging tool, was used to capture the molecular fingerprints of overall cellular component responses (nucleic acid, lipids, proteins, carbohydrates) to DEP damage and antioxidant protection. We apply a novel data visualization algorithm called PHATE, which preserves both global and local structure, to display the progression of cell damage over DEP exposure time. Meanwhile, a mutual information (MI) estimator was used to identify the most informative Raman peaks associated with cytotoxicity. A health index was defined to quantitatively assess the protective effects of two antioxidants (resveratrol and mesobiliverdin IXα) against DEP induced cytotoxicity. In addition, a number of machine learning classifiers were applied to successfully discriminate different treatment groups with high accuracy. Correlations between Raman spectra and immunomodulatory cytokine and chemokine levels were evaluated. In conclusion, the combination of label-free, non-disruptive Raman micro-spectroscopy and machine learning analysis is demonstrated as a useful tool in quantitative analysis of oxidative stress induced cytotoxicity and for effectively assessing various antioxidant treatments, suggesting that this framework can serve as a high throughput platform for screening various potential antioxidants based on their effectiveness at battling the effects of air pollution on human health.

Single cell immune profiling of dengue virus patients reveals intact immune responses to Zika virus with enrichment of innate immune signatures.

The genus Flavivirus contains many mosquito-borne human pathogens of global epidemiological importance such as dengue virus, West Nile virus, and Zika virus, which has recently emerged at epidemic levels. Infections with these viruses result in divergent clinical outcomes ranging from asymptomatic to fatal. Myriad factors influence infection severity including exposure, immune status and pathogen/host genetics. Furthermore, pre-existing infection may skew immune pathways or divert immune resources. We profiled immune cells from dengue virus-infected individuals by multiparameter mass cytometry (CyTOF) to define functional status. Elevations in IFNß were noted in acute patients across the majority of cell types and were statistically elevated in 31 of 36 cell subsets. We quantified response to in vitro (re)infection with dengue or Zika viruses and detected a striking pattern of upregulation of responses to Zika infection by innate cell types which was not noted in response to dengue virus. Significance was discovered by statistical analysis as well as a neural network-based clustering approach which identified unusual cell subsets overlooked by conventional manual gating. Of public health importance, patient cells showed significant enrichment of innate cell responses to Zika virus indicating an intact and robust anti-Zika response despite the concurrent dengue infection.

Author Correction: Visualizing structure and transitions in high-dimensional biological data.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Visualizing structure and transitions in high-dimensional biological data.

The high-dimensional data created by high-throughput technologies require visualization tools that reveal data structure and patterns in an intuitive form. We present PHATE, a visualization method that captures both local and global nonlinear structure using an information-geometric distance between data points. We compare PHATE to other tools on a variety of artificial and biological datasets, and find that it consistently preserves a range of patterns in data, including continual progressions, branches and clusters, better than other tools. We define a manifold preservation metric, which we call denoised embedding manifold preservation (DEMaP), and show that PHATE produces lower-dimensional embeddings that are quantitatively better denoised as compared to existing visualization methods. An analysis of a newly generated single-cell RNA sequencing dataset on human germ-layer differentiation demonstrates how PHATE reveals unique biological insight into the main developmental branches, including identification of three previously undescribed subpopulations. We also show that PHATE is applicable to a wide variety of data types, including mass cytometry, single-cell RNA sequencing, Hi-C and gut microbiome data.

Exploring single-cell data with deep multitasking neural networks.

It is currently challenging to analyze single-cell data consisting of many cells and samples, and to address variations arising from batch effects and different sample preparations. For this purpose, we present SAUCIE, a deep neural network that combines parallelization and scalability offered by neural networks, with the deep representation of data that can be learned by them to perform many single-cell data analysis tasks. Our regularizations (penalties) render features learned in hidden layers of the neural network interpretable. On large, multi-patient datasets, SAUCIE's various hidden layers contain denoised and batch-corrected data, a low-dimensional visualization and unsupervised clustering, as well as other information that can be used to explore the data. We analyze a 180-sample dataset consisting of 11 million T cells from dengue patients in India, measured with mass cytometry. SAUCIE can batch correct and identify cluster-based signatures of acute dengue infection and create a patient manifold, stratifying immune response to dengue.

Dissecting alterations in human CD8+ T cells with aging by high-dimensional single cell mass cytometry.

We investigated the effect of aging on the multi-dimensional characteristics and heterogeneity of human peripheral CD8+ T cells defined by the expression of a set of molecules at the single cell level using the recently developed mass cytometry or Cytometry by Time-Of-Flight (CyTOF) and computational algorithms. CD8+ T cells of young and older adults had differential expression of molecules, especially those related to cell activation and migration, permitting the clustering of young and older adults through an unbiased approach. The changes in the expression of individual molecules were collectively reflected in the altered high-dimensional profiles of CD8+ T cells in older adults as visualized by the dimensionality reduction analysis tools principal component analysis (PCA) and t-distributed stochastic neighbor embedding (t-SNE). A combination of PhenoGraph clustering and t-SNE analysis revealed heterogeneous subsets of CD8+ T cells that altered with aging. Furthermore, intermolecular quantitative relationships in CD8+ T cells appeared to change with age as determined by the computational algorithm conditional-Density Resampled Estimate of Mutual Information (DREMI). The results of our study showed that heterogeneity, multidimensional characteristics, and intermolecular quantitative relationships in human CD8+ T cells altered with age, distinctively clustering young and older adults through an unbiased approach.

Coarse Graining of Data via Inhomogeneous Diffusion Condensation.

Big data often has emergent structure that exists at multiple levels of abstraction, which are useful for characterizing complex interactions and dynamics of the observations. Here, we consider multiple levels of abstraction via a multiresolution geometry of data points at different granularities. To construct this geometry we define a time-inhomogemeous diffusion process that effectively condenses data points together to uncover nested groupings at larger and larger granularities. This inhomogeneous process creates a deep cascade of intrinsic low pass filters on the data affinity graph that are applied in sequence to gradually eliminate local variability while adjusting the learned data geometry to increasingly coarser resolutions. We provide visualizations to exhibit our method as a "continuously-hierarchical" clustering with directions of eliminated variation highlighted at each step. The utility of our algorithm is demonstrated via neuronal data condensation, where the constructed multiresolution data geometry uncovers the organization, grouping, and connectivity between neurons.

Recovering Gene Interactions from Single-Cell Data Using Data Diffusion.

Single-cell RNA sequencing technologies suffer from many sources of technical noise, including under-sampling of mRNA molecules, often termed "dropout," which can severely obscure important gene-gene relationships. To address this, we developed MAGIC (Markov affinity-based graph imputation of cells), a method that shares information across similar cells, via data diffusion, to denoise the cell count matrix and fill in missing transcripts. We validate MAGIC on several biological systems and find it effective at recovering gene-gene relationships and additional structures. Applied to the epithilial to mesenchymal transition, MAGIC reveals a phenotypic continuum, with the majority of cells residing in intermediate states that display stem-like signatures, and infers known and previously uncharacterized regulatory interactions, demonstrating that our approach can successfully uncover regulatory relations without perturbations.

Ensemble Estimation of Information Divergence †.

Recent work has focused on the problem of nonparametric estimation of information divergence functionals between two continuous random variables. Many existing approaches require either restrictive assumptions about the density support set or difficult calculations at the support set boundary which must be known a priori. The mean squared error (MSE) convergence rate of a leave-one-out kernel density plug-in divergence functional estimator for general bounded density support sets is derived where knowledge of the support boundary, and therefore, the boundary correction is not required. The theory of optimally weighted ensemble estimation is generalized to derive a divergence estimator that achieves the parametric rate when the densities are sufficiently smooth. Guidelines for the tuning parameter selection and the asymptotic distribution of this estimator are provided. Based on the theory, an empirical estimator of Rényi-α divergence is proposed that greatly outperforms the standard kernel density plug-in estimator in terms of mean squared error, especially in high dimensions. The estimator is shown to be robust to the choice of tuning parameters. We show extensive simulation results that verify the theoretical results of our paper. Finally, we apply the proposed estimator to estimate the bounds on the Bayes error rate of a cell classification problem.

THE INTRINSIC VALUE OF HFO FEATURES AS A BIOMARKER OF EPILEPTIC ACTIVITY.

High frequency oscillations (HFOs) are a promising biomarker of epileptic brain tissue and activity. HFOs additionally serve as a prototypical example of challenges in the analysis of discrete events in high-temporal resolution, intracranial EEG data. Two primary challenges are 1) dimensionality reduction, and 2) assessing feasibility of classification. Dimensionality reduction assumes that the data lie on a manifold with dimension less than that of the features space. However, previous HFO analysis have assumed a linear manifold, global across time, space (i.e. recording electrode/channel), and individual patients. Instead, we assess both a) whether linear methods are appropriate and b) the consistency of the manifold across time, space, and patients. We also estimate bounds on the Bayes classification error to quantify the distinction between two classes of HFOs (those occurring during seizures and those occurring due to other processes). This analysis provides the foundation for future clinical use of HFO features and guides the analysis for other discrete events, such as individual action potentials or multi-unit activity.

Total building commissioning.

Healthcare facilities are generally commissioned according to standards laid down by specific stipulations or the consultant engineers. This commissioning process is normally only concerned with the physical functioning of the building services and is aimed at ensuring the installed services are operating according to the engineers' specifications. It is no longer acceptable just to conduct a thorough clean, commonly called a hospital clean, and consider that as adequate for patient occupation, writes Kevin Moon.