A survey of some tensor analysis techniques for biological systems

The structure is the message: Preserving experimental context through tensor decomposition

Recent biological studies have been revolutionized in scale and granularity by multiplex and high-throughput assays. Profiling cell responses across several experimental parameters, such as perturbations, time, and genetic contexts, leads to richer and more generalizable findings. However, these multidimensional datasets necessitate a reevaluation of the conventional methods for their representation and analysis. Traditionally, experimental parameters are merged to flatten the data into a two-dimensional matrix, sacrificing crucial experiment context reflected by the structure. As Marshall McLuhan famously stated, "the medium is the message." In this work, we propose that the experiment structure is the medium in which subsequent analysis is performed, and the optimal choice of data representation must reflect the experiment structure. We review how tensor-structured analyses and decompositions can preserve this information. We contend that tensor methods are poised to become integral to the biomedical data sciences toolkit.

Integrative, high-resolution analysis of single cells across experimental conditions with PARAFAC2

Effective tools for exploration and analysis are needed to extract insights from large-scale single-cell measurement data. However, current techniques for handling single-cell studies performed across experimental conditions (e.g., samples, perturbations, or patients) require restrictive assumptions, lack flexibility, or do not adequately deconvolute condition-to-condition variation from cell-to-cell variation. Here, we report that the tensor decomposition method PARAFAC2 (Pf2) enables the dimensionality reduction of single-cell data across conditions. We demonstrate these benefits across two distinct contexts of single-cell RNA-sequencing (scRNA-seq) experiments of peripheral immune cells: pharmacologic drug perturbations and systemic lupus erythematosus (SLE) patient samples. By isolating relevant gene modules across cells and conditions, Pf2 enables straightforward associations of gene variation patterns across specific patients or perturbations while connecting each coordinated change to certain cells without pre-defining cell types. The theoretical grounding of Pf2 suggests a unified framework for many modeling tasks associated with single-cell data. Thus, Pf2 provides an intuitive universal dimensionality reduction approach for multi-sample single-cell studies across diverse biological contexts.

Censored Least Squares for Imputing Missing Values in PARAFAC Tensor Factorization

Tensor factorization is a dimensionality reduction method applied to multidimensional arrays. These methods are useful for identifying patterns within a variety of biomedical datasets due to their ability to preserve the organizational structure of experiments and therefore aid in generating meaningful insights. However, missing data in the datasets being analyzed can impose challenges. Tensor factorization can be performed with some level of missing data and reconstruct a complete tensor. However, while tensor methods may impute these missing values, the choice of fitting algorithm may influence the fidelity of these imputations. Previous approaches, based on alternating least squares with prefilled values or direct optimization, suffer from introduced bias or slow computational performance. In this study, we propose that censored least squares can better handle missing values with data structured in tensor form. We ran censored least squares on four different biological datasets and compared its performance against alternating least squares with prefilled values and direct optimization. We used the error of imputation and the ability to infer masked values to benchmark their missing data performance. Censored least squares appeared best suited for the analysis of high-dimensional biological data by accuracy and convergence metrics across several studies.

Frataxin Deficit Leads to Reduced Dynamics of Growth Cones in Dorsal Root Ganglia Neurons of Friedreich’s Ataxia YG8sR Model: A Multilinear Algebra Approach

Computational techniques for analyzing biological images offer a great potential to enhance our knowledge of the biological processes underlying disorders of the nervous system. Friedreich’s Ataxia (FRDA) is a rare progressive neurodegenerative inherited disorder caused by the low expression of frataxin, which is a small mitochondrial protein. In FRDA cells, the lack of frataxin promotes primarily mitochondrial dysfunction, an alteration of calcium (Ca²⁺) homeostasis and the destabilization of the actin cytoskeleton in the neurites and growth cones of sensory neurons. In this paper, a computational multilinear algebra approach was used to analyze the dynamics of the growth cone and its function in control and FRDA neurons. Computational approach, which includes principal component analysis and a multilinear algebra method, is used to quantify the dynamics of the growth cone (GC) morphology of sensory neurons from the dorsal root ganglia (DRG) of the YG8sR humanized murine model for FRDA. It was confirmed that the dynamics and patterns of turning were aberrant in the FRDA growth cones. In addition, our data suggest that other cellular processes dependent on functional GCs such as axonal regeneration might also be affected. Semiautomated computational approaches are presented to quantify differences in GC behaviors in neurodegenerative disease. In summary, the deficiency of frataxin has an adverse effect on the formation and, most importantly, the growth cones’ function in adult DRG neurons. As a result, frataxin deficient DRG neurons might lose the intrinsic capability to grow and regenerate axons properly due to the dysfunctional GCs they build.

Identification of genes associated with altered gene expression and m6A profiles during hypoxia using tensor decomposition based unsupervised feature extraction

Although hypoxia is a critical factor that can drive the progression of various diseases, the mechanism underlying hypoxia itself remains unclear. Recently, m6A has been proposed as an important factor driving hypoxia. Despite successful analyses, potential genes were not selected with statistical significance but were selected based solely on fold changes. Because the number of genes is large while the number of samples is small, it was impossible to select genes using conventional feature selection methods with statistical significance. In this study, we applied the recently proposed principal component analysis (PCA), tensor decomposition (TD), and kernel tensor decomposition (KTD)-based unsupervised feature extraction (FE) to a hypoxia data set. We found that PCA, TD, and KTD-based unsupervised FE could successfully identify a limited number of genes associated with altered gene expression and m6A profiles, as well as the enrichment of hypoxia-related biological terms, with improved statistical significance.