Generative Distribution Embeddings: Lifting autoencoders to the space of distributions for multiscale representation learning

TL;DR

The paper introduces generative distribution embeddings (GDE), a framework that extends autoencoders to operate on distributions, enabling multiscale representation learning and applications in computational biology.

Contribution

GDE lifts autoencoders to the space of distributions, coupling generative models with encoders to learn distributional representations satisfying distributional invariance.

Findings

GDEs learn predictive sufficient statistics in Wasserstein space.

Latent GDE distances approximate the Wasserstein-2 distance.

GDEs outperform existing methods on synthetic benchmarks.

Abstract

Many real-world problems require reasoning across multiple scales, demanding models which operate not on single data points, but on entire distributions. We introduce generative distribution embeddings (GDE), a framework that lifts autoencoders to the space of distributions. In GDEs, an encoder acts on sets of samples, and the decoder is replaced by a generator which aims to match the input distribution. This framework enables learning representations of distributions by coupling conditional generative models with encoder networks which satisfy a criterion we call distributional invariance. We show that GDEs learn predictive sufficient statistics embedded in the Wasserstein space, such that latent GDE distances approximately recover the $W_2$ distance, and latent interpolation approximately recovers optimal transport trajectories for Gaussian and Gaussian mixture distributions. We systematically benchmark GDEs against existing approaches on synthetic datasets, demonstrating consistently stronger performance. We then apply GDEs to six key problems in computational biology: learning donor-level representations from single-nuclei RNA sequencing data (6M cells), capturing clonal dynamics in lineage-traced RNA sequencing data (150K cells), predicting perturbation effects on transcriptomes (1M cells), predicting perturbation effects on cellular phenotypes (20M single-cell images), designing synthetic yeast promoters (34M sequences), and spatiotemporal modeling of viral protein sequences (1M sequences).