Optimal-Transport Analysis of Single-Cell Gene Expression Identifies Developmental Trajectories in Reprogramming.

Understanding the molecular programs that guide differentiation during development is a major challenge. Here, we introduce Waddington-OT, an approach for studying developmental time courses to infer ancestor-descendant fates and model the regulatory programs that underlie them. We apply the method to reconstruct the landscape of reprogramming from 315,000 single-cell RNA sequencing (scRNA-seq) profiles, collected at half-day intervals across 18 days. The results reveal a wider range of developmental programs than previously characterized. Cells gradually adopt either a terminal stromal state or a mesenchymal-to-epithelial transition state. The latter gives rise to populations related to pluripotent, extra-embryonic, and neural cells, with each harboring multiple finer subpopulations. The analysis predicts transcription factors and paracrine signals that affect fates and experiments validate that the TF Obox6 and the cytokine GDF9 enhance reprogramming efficiency. Our approach sheds light on the process and outcome of reprogramming and provides a framework applicable to diverse temporal processes in biology.

Type I interferon activates MHC class I-dressed CD11b+ conventional dendritic cells to promote protective anti-tumor CD8+ T cell immunity.

Tumor-infiltrating dendritic cells (DCs) assume varied functional states that impact anti-tumor immunity. To delineate the DC states associated with productive anti-tumor T cell immunity, we compared spontaneously regressing and progressing tumors. Tumor-reactive CD8+ T cell responses in Batf3-/- mice lacking type 1 DCs (DC1s) were lost in progressor tumors but preserved in regressor tumors. Transcriptional profiling of intra-tumoral DCs within regressor tumors revealed an activation state of CD11b+ conventional DCs (DC2s) characterized by expression of interferon (IFN)-stimulated genes (ISGs) (ISG+ DCs). ISG+ DC-activated CD8+ T cells ex vivo comparably to DC1. Unlike cross-presenting DC1, ISG+ DCs acquired and presented intact tumor-derived peptide-major histocompatibility complex class I (MHC class I) complexes. Constitutive type I IFN production by regressor tumors drove the ISG+ DC state, and activation of MHC class I-dressed ISG+ DCs by exogenous IFN-ß rescued anti-tumor immunity against progressor tumors in Batf3-/- mice. The ISG+ DC gene signature is detectable in human tumors. Engaging this functional DC state may present an approach for the treatment of human disease.

Multi-axis fields boost SABRE hyperpolarization.

The inherently low signal-to-noise ratio of NMR and MRI is now being addressed by hyperpolarization methods. For example, iridium-based catalysts that reversibly bind both parahydrogen and ligands in solution can hyperpolarize protons (SABRE) or heteronuclei (X-SABRE) on a wide variety of ligands, using a complex interplay of spin dynamics and chemical exchange processes, with common signal enhancements between 103 and 104. This does not approach obvious theoretical limits, and further enhancement would be valuable in many applications (such as imaging mM concentration species in vivo). Most SABRE/X-SABRE implementations require far lower fields (µT-mT) than standard magnetic resonance (>1T), and this gives an additional degree of freedom: the ability to fully modulate fields in three dimensions. However, this has been underexplored because the standard simplifying theoretical assumptions in magnetic resonance need to be revisited. Here, we take a different approach, an evolutionary strategy algorithm for numerical optimization, multi-axis computer-aided heteronuclear transfer enhancement for SABRE (MACHETE-SABRE). We find nonintuitive but highly efficient multiaxial pulse sequences which experimentally can produce a sevenfold improvement in polarization over continuous excitation. This approach optimizes polarization differently than traditional methods, thus gaining extra efficiency.

Optimal reaching subject to computational and physical constraints reveals structure of the sensorimotor control system.

Optimal feedback control provides an abstract framework describing the architecture of the sensorimotor system without prescribing implementation details such as what coordinate system to use, how feedback is incorporated, or how to accommodate changing task complexity. We investigate how such details are determined by computational and physical constraints by creating a model of the upper limb sensorimotor system in which all connection weights between neurons, feedback, and muscles are unknown. By optimizing these parameters with respect to an objective function, we find that the model exhibits a preference for an intrinsic (joint angle) coordinate representation of inputs and feedback and learns to calculate a weighted feedforward and feedback error. We further show that complex reaches around obstacles can be achieved by augmenting our model with a path-planner based on via points. The path-planner revealed "avoidance" neurons that encode directions to reach around obstacles and "placement" neurons that make fine-tuned adjustments to via point placement. Our results demonstrate the surprising capability of computationally constrained systems and highlight interesting characteristics of the sensorimotor system.

Hemodynamic regulation allows stable growth of microvascular networks.

How do vessels find optimal radii? Capillaries are known to adapt their radii to maintain the shear stress of blood flow at the vessel wall at a set point, yet models of adaptation purely based on average shear stress have not been able to produce complex loopy networks that resemble real microvascular systems. For narrow vessels where red blood cells travel in a single file, the shear stress on vessel endothelium peaks sharply when a red blood cell passes through. We show that stable shear-stress-based adaptation is possible if vessel shear stress set points are cued to the stress peaks. Model networks that respond to peak stresses alone can quantitatively reproduce the observed zebrafish trunk microcirculation, including its adaptive trajectory when hematocrit changes or parts of the network are amputated. Our work reveals the potential for mechanotransduction alone to generate stable hydraulically tuned microvascular networks.

Predicting resilience of migratory birds to environmental change.

The pace and scale of environmental change represent major challenges to many organisms. Animals that move long distances, such as migratory birds, are especially vulnerable to change since they need chains of intact habitat along their migratory routes. Estimating the resilience of such species to environmental changes assists in targeting conservation efforts. We developed a migration modeling framework to predict past (1960s), present (2010s), and future (2060s) optimal migration strategies across five shorebird species (Scolopacidae) within the East Asian-Australasian Flyway, which has seen major habitat deterioration and loss over the last century, and compared these predictions to empirical tracks from the present. Our model captured the migration strategies of the five species and identified the changes in migrations needed to respond to habitat deterioration and climate change. Notably, the larger species, with single or few major stopover sites, need to establish new migration routes and strategies, while smaller species can buffer habitat loss by redistributing their stopover areas to novel or less-used sites. Comparing model predictions with empirical tracks also indicates that larger species with the stronger need for adaptations continue to migrate closer to the optimal routes of the past, before habitat deterioration accelerated. Our study not only quantifies the vulnerability of species in the face of global change but also explicitly reveals the extent of adaptations required to sustain their migrations. This modeling framework provides a tool for conservation planning that can accommodate the future needs of migratory species.

Multiple and Optimal Screening Subset: a method selecting global characteristic congeners for robust foodomics analysis.

Metabolomics and foodomics shed light on the molecular processes within living organisms and the complex food composition by leveraging sophisticated analytical techniques to systematically analyze the vast array of molecular features. The traditional feature-picking method often results in arbitrary selections of the model, feature ranking, and cut-off, which may lead to suboptimal results. Thus, a Multiple and Optimal Screening Subset (MOSS) approach was developed in this study to achieve a balance between a minimal number of predictors and high predictive accuracy during statistical model setup. The MOSS approach compares five commonly used models in the context of food matrix analysis, specifically bourbons. These models include Student's t-test, receiver operating characteristic curve, partial least squares-discriminant analysis (PLS-DA), random forests, and support vector machines. The approach employs cross-validation to identify promising subset feature candidates that contribute to food characteristic classification. It then determines the optimal subset size by comparing it to the corresponding top-ranked features. Finally, it selects the optimal feature subset by traversing all possible feature candidate combinations. By utilizing MOSS approach to analyze 1406 mass spectral features from a collection of 122 bourbon samples, we were able to generate a subset of features for bourbon age prediction with 88% accuracy. Additionally, MOSS increased the area under the curve performance of sweetness prediction to 0.898 with only four predictors compared with the top-ranked four features at 0.681 based on the PLS-DA model. Overall, we demonstrated that MOSS provides an efficient and effective approach for selecting optimal features compared with other frequently utilized methods.

BayesKAT: bayesian optimal kernel-based test for genetic association studies reveals joint genetic effects in complex diseases.

Genome-wide Association Studies (GWAS) methods have identified individual single-nucleotide polymorphisms (SNPs) significantly associated with specific phenotypes. Nonetheless, many complex diseases are polygenic and are controlled by multiple genetic variants that are usually non-linearly dependent. These genetic variants are marginally less effective and remain undetected in GWAS analysis. Kernel-based tests (KBT), which evaluate the joint effect of a group of genetic variants, are therefore critical for complex disease analysis. However, choosing different kernel functions in KBT can significantly influence the type I error control and power, and selecting the optimal kernel remains a statistically challenging task. A few existing methods suffer from inflated type 1 errors, limited scalability, inferior power or issues of ambiguous conclusions. Here, we present a new Bayesian framework, BayesKAT (https://github.com/wangjr03/BayesKAT), which overcomes these kernel specification issues by selecting the optimal composite kernel adaptively from the data while testing genetic associations simultaneously. Furthermore, BayesKAT implements a scalable computational strategy to boost its applicability, especially for high-dimensional cases where other methods become less effective. Based on a series of performance comparisons using both simulated and real large-scale genetics data, BayesKAT outperforms the available methods in detecting complex group-level associations and controlling type I errors simultaneously. Applied on a variety of groups of functionally related genetic variants based on biological pathways, co-expression gene modules and protein complexes, BayesKAT deciphers the complex genetic basis and provides mechanistic insights into human diseases.

SPANN: annotating single-cell resolution spatial transcriptome data with scRNA-seq data.

MOTIVATION: The rapid development of spatial transcriptome technologies has enabled researchers to acquire single-cell-level spatial data at an affordable price. However, computational analysis tools, such as annotation tools, tailored for these data are still lacking. Recently, many computational frameworks have emerged to integrate single-cell RNA sequencing (scRNA-seq) and spatial transcriptomics datasets. While some frameworks can utilize well-annotated scRNA-seq data to annotate spatial expression patterns, they overlook critical aspects. First, existing tools do not explicitly consider cell type mapping when aligning the two modalities. Second, current frameworks lack the capability to detect novel cells, which remains a key interest for biologists. RESULTS: To address these problems, we propose an annotation method for spatial transcriptome data called SPANN. The main tasks of SPANN are to transfer cell-type labels from well-annotated scRNA-seq data to newly generated single-cell resolution spatial transcriptome data and discover novel cells from spatial data. The major innovations of SPANN come from two aspects: SPANN automatically detects novel cells from unseen cell types while maintaining high annotation accuracy over known cell types. SPANN finds a mapping between spatial transcriptome samples and RNA data prototypes and thus conducts cell-type-level alignment. Comprehensive experiments using datasets from various spatial platforms demonstrate SPANN's capabilities in annotating known cell types and discovering novel cell states within complex tissue contexts. AVAILABILITY: The source code of SPANN can be accessed at https://github.com/ddb-qiwang/SPANN-torch. CONTACT: dengmh@math.pku.edu.cn.

scDOT: enhancing single-cell RNA-Seq data annotation and uncovering novel cell types through multi-reference integration.

The proliferation of single-cell RNA-seq data has greatly enhanced our ability to comprehend the intricate nature of diverse tissues. However, accurately annotating cell types in such data, especially when handling multiple reference datasets and identifying novel cell types, remains a significant challenge. To address these issues, we introduce Single Cell annotation based on Distance metric learning and Optimal Transport (scDOT), an innovative cell-type annotation method adept at integrating multiple reference datasets and uncovering previously unseen cell types. scDOT introduces two key innovations. First, by incorporating distance metric learning and optimal transport, it presents a novel optimization framework. This framework effectively learns the predictive power of each reference dataset for new query data and simultaneously establishes a probabilistic mapping between cells in the query data and reference-defined cell types. Secondly, scDOT develops an interpretable scoring system based on the acquired probabilistic mapping, enabling the precise identification of previously unseen cell types within the data. To rigorously assess scDOT's capabilities, we systematically evaluate its performance using two diverse collections of benchmark datasets encompassing various tissues, sequencing technologies and diverse cell types. Our experimental results consistently affirm the superior performance of scDOT in cell-type annotation and the identification of previously unseen cell types. These advancements provide researchers with a potent tool for precise cell-type annotation, ultimately enriching our understanding of complex biological tissues.

Improving the use of plant genetic resources to sustain breeding programs' efficiency.

Genetic progress of crop plants is required to face human population growth and guarantee production stability in increasingly unstable environmental conditions. Breeding is accompanied by a loss in genetic diversity, which hinders sustainable genetic gain. Methodologies based on molecular marker information have been developed to manage diversity and proved effective in increasing long-term genetic gain. However, with realistic plant breeding population sizes, diversity depletion in closed programs appears ineluctable, calling for the introduction of relevant diversity donors. Although maintained with significant efforts, genetic resource collections remain underutilized, due to a large performance gap with elite germplasm. Bridging populations created by crossing genetic resources to elite lines prior to introduction into elite programs can manage this gap efficiently. To improve this strategy, we explored with simulations different genomic prediction and genetic diversity management options for a global program involving a bridging and an elite component. We analyzed the dynamics of quantitative trait loci fixation and followed the fate of allele donors after their introduction into the breeding program. Allocating 25% of total experimental resources to create a bridging component appears highly beneficial. We showed that potential diversity donors should be selected based on their phenotype rather than genomic predictions calibrated with the ongoing breeding program. We recommend incorporating improved donors into the elite program using a global calibration of the genomic prediction model and optimal cross selection maintaining a constant diversity. These approaches use efficiently genetic resources to sustain genetic gain and maintain neutral diversity, improving the flexibility to address future breeding objectives.

Internal feedback in the cortical perception-action loop enables fast and accurate behavior.

Animals move smoothly and reliably in unpredictable environments. Models of sensorimotor control, drawing on control theory, have assumed that sensory information from the environment leads to actions, which then act back on the environment, creating a single, unidirectional perception-action loop. However, the sensorimotor loop contains internal delays in sensory and motor pathways, which can lead to unstable control. We show here that these delays can be compensated by internal feedback signals that flow backward, from motor toward sensory areas. This internal feedback is ubiquitous in neural sensorimotor systems, and we show how internal feedback compensates internal delays. This is accomplished by filtering out self-generated and other predictable changes so that unpredicted, actionable information can be rapidly transmitted toward action by the fastest components, effectively compressing the sensory input to more efficiently use feedforward pathways: Tracts of fast, giant neurons necessarily convey less accurate signals than tracts with many smaller neurons, but they are crucial for fast and accurate behavior. We use a mathematically tractable control model to show that internal feedback has an indispensable role in achieving state estimation, localization of function (how different parts of the cortex control different parts of the body), and attention, all of which are crucial for effective sensorimotor control. This control model can explain anatomical, physiological, and behavioral observations, including motor signals in the visual cortex, heterogeneous kinetics of sensory receptors, and the presence of giant cells in the cortex of humans as well as internal feedback patterns and unexplained heterogeneity in neural systems.

Erroneous Compensation for Long-Latency Feedback Delays as Origin of Essential Tremor.

Essential tremor (ET), a movement disorder characterized by involuntary oscillations of the limbs during movement, remains to date not well understood. It has been recently suggested that the tremor originates from impaired delay compensation, affecting movement representation and online control. Here we tested this hypothesis directly with 24 ET patients (14 female; 10 male) and 28 neurologically intact (NI) human volunteers (17 female; 11 male) in an upper limb postural perturbation task. After maintaining their hand in a visual target, participants experienced perturbations of unpredictable direction and magnitude and were instructed to counter the perturbation and steer their hand back to the starting position. In comparison with NI volunteers, ET patients' early muscular responses (short and long-latency responses, 20-50 and 50-100 ms, respectively) were preserved or even slightly increased. However, they exhibited perturbation-dependent deficits when stopping and stabilizing their hand in the final target supporting the hypothesis that the tremor was generated by the feedback controller. We show in a computational model that errors in delay compensation accumulating over time produced the same small increase in initial feedback response followed by oscillations that scaled with the perturbation magnitude as observed in ET population. Our experimental results therefore validate the computational hypothesis that inaccurate delay compensation in long-latency pathways could be the origin of the tremor.

Lessons Learned from the ISCHEMIA Trial for the Management of Patients with Stable Ischemic Heart Disease.

The recent landmark International Study of Comparative Health Effectiveness With Medical and Invasive Approaches (ISCHEMIA) trial was undertaken to assess whether stable angina patients with moderate to severe baseline ischemia would benefit from an invasive approach with revascularization versus a conservative approach of intensive lifestyle intervention and pharmacologic secondary prevention. This trial addressed the hypothesis that treating ischemia with an invasive approach would reduce major adverse cardiac events more than a noninvasive pharmacologic and lifestyle approach. ISCHEMIA is discussed in detail, along with current implications for contemporary management of this very common cardiac disorder afflicting millions of patients worldwide.

Uncertainty directed factorial clinical trials.

The development and evaluation of novel treatment combinations is a key component of modern clinical research. The primary goals of factorial clinical trials of treatment combinations range from the estimation of intervention-specific effects, or the discovery of potential synergies, to the identification of combinations with the highest response probabilities. Most factorial studies use balanced or block randomization, with an equal number of patients assigned to each treatment combination, irrespective of the specific goals of the trial. Here, we introduce a class of Bayesian response-adaptive designs for factorial clinical trials with binary outcomes. The study design was developed using Bayesian decision-theoretic arguments and adapts the randomization probabilities to treatment combinations during the enrollment period based on the available data. Our approach enables the investigator to specify a utility function representative of the aims of the trial, and the Bayesian response-adaptive randomization algorithm aims to maximize this utility function. We considered several utility functions and factorial designs tailored to them. Then, we conducted a comparative simulation study to illustrate relevant differences of key operating characteristics across the resulting designs. We also investigated the asymptotic behavior of the proposed adaptive designs. We also used data summaries from three recent factorial trials in perioperative care, smoking cessation, and infectious disease prevention to define realistic simulation scenarios and illustrate advantages of the introduced trial designs compared to other study designs.

A fast and globally optimal solution for RNA-seq quantification.

Alignment-based RNA-seq quantification methods typically involve a time-consuming alignment process prior to estimating transcript abundances. In contrast, alignment-free RNA-seq quantification methods bypass this step, resulting in significant speed improvements. Existing alignment-free methods rely on the Expectation-Maximization (EM) algorithm for estimating transcript abundances. However, EM algorithms only guarantee locally optimal solutions, leaving room for further accuracy improvement by finding a globally optimal solution. In this study, we present TQSLE, the first alignment-free RNA-seq quantification method that provides a globally optimal solution for transcript abundances estimation. TQSLE adopts a two-step approach: first, it constructs a k-mer frequency matrix A for the reference transcriptome and a k-mer frequency vector b for the RNA-seq reads; then, it directly estimates transcript abundances by solving the linear equation ATAx = ATb. We evaluated the performance of TQSLE using simulated and real RNA-seq data sets and observed that, despite comparable speed to other alignment-free methods, TQSLE outperforms them in terms of accuracy. TQSLE is freely available at https://github.com/yhg926/TQSLE.

iEssLnc: quantitative estimation of lncRNA gene essentialities with meta-path-guided random walks on the lncRNA-protein interaction network.

Gene essentiality is defined as the extent to which a gene is required for the survival and reproductive success of a living system. It can vary between genetic backgrounds and environments. Essential protein coding genes have been well studied. However, the essentiality of non-coding regions is rarely reported. Most regions of human genome do not encode proteins. Determining essentialities of non-coding genes is demanded. We developed iEssLnc models, which can assign essentiality scores to lncRNA genes. As far as we know, this is the first direct quantitative estimation to the essentiality of lncRNA genes. By taking the advantage of graph neural network with meta-path-guided random walks on the lncRNA-protein interaction network, iEssLnc models can perform genome-wide screenings for essential lncRNA genes in a quantitative manner. We carried out validations and whole genome screening in the context of human cancer cell-lines and mouse genome. In comparisons to other methods, which are transferred from protein-coding genes, iEssLnc achieved better performances. Enrichment analysis indicated that iEssLnc essentiality scores clustered essential lncRNA genes with high ranks. With the screening results of iEssLnc models, we estimated the number of essential lncRNA genes in human and mouse. We performed functional analysis to find that essential lncRNA genes interact with microRNAs and cytoskeletal proteins significantly, which may be of interest in experimental life sciences. All datasets and codes of iEssLnc models have been deposited in GitHub (https://github.com/yyZhang14/iEssLnc).

Matching single cells across modalities with contrastive learning and optimal transport.

Understanding the interactions between the biomolecules that govern cellular behaviors remains an emergent question in biology. Recent advances in single-cell technologies have enabled the simultaneous quantification of multiple biomolecules in the same cell, opening new avenues for understanding cellular complexity and heterogeneity. Still, the resulting multimodal single-cell datasets present unique challenges arising from the high dimensionality and multiple sources of acquisition noise. Computational methods able to match cells across different modalities offer an appealing alternative towards this goal. In this work, we propose MatchCLOT, a novel method for modality matching inspired by recent promising developments in contrastive learning and optimal transport. MatchCLOT uses contrastive learning to learn a common representation between two modalities and applies entropic optimal transport as an approximate maximum weight bipartite matching algorithm. Our model obtains state-of-the-art performance on two curated benchmarking datasets and an independent test dataset, improving the top scoring method by 26.1% while preserving the underlying biological structure of the multimodal data. Importantly, MatchCLOT offers high gains in computational time and memory that, in contrast to existing methods, allows it to scale well with the number of cells. As single-cell datasets become increasingly large, MatchCLOT offers an accurate and efficient solution to the problem of modality matching.

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.

Gene Architectures that Minimize Cost of Gene Expression.

Gene expression burdens cells by consuming resources and energy. While numerous studies have investigated regulation of expression level, little is known about gene design elements that govern expression costs. Here, we ask how cells minimize production costs while maintaining a given protein expression level and whether there are gene architectures that optimize this process. We measured fitness of ∼14,000 E. coli strains, each expressing a reporter gene with a unique 5' architecture. By comparing cost-effective and ineffective architectures, we found that cost per protein molecule could be minimized by lowering transcription levels, regulating translation speeds, and utilizing amino acids that are cheap to synthesize and that are less hydrophobic. We then examined natural E. coli genes and found that highly expressed genes have evolved more forcefully to minimize costs associated with their expression. Our study thus elucidates gene design elements that improve the economy of protein expression in natural and heterologous systems.