1 of 90

The Physics of Transcriptomes

MCB137L/237L

Spring 2025

2 of 90

What Does the mRNA Distribution Tell Us About How Transcription Happens?

Zenklusen et al. (2008)

3 of 90

Testing the Null Hypothesis�Deviations from Poisson Reveal Molecular Mechanism

Zenklusen et al. (2008)

4 of 90

And, Yet, Our Ignorance is Vast

5 of 90

Regulatory Ignorance Throughout Our�Model Organisms

6 of 90

The Cell as a Bag of RNA

7 of 90

Technologies Driving the DNA Sequencing Revolution

8 of 90

From Gene-Wide to�Genome-Wide Studies

9 of 90

A Deluge of Sequencing Data

10 of 90

A Feeling for the Scale of Our Sequencing Data

11 of 90

Estimating Book Lengths

12 of 90

Our Collective Memory:�The Library of Congress

13 of 90

14 of 90

15 of 90

Number of Letters and Their Meaning

16 of 90

The Sequence Read Archive Versus Shakespeare

17 of 90

18 of 90

19 of 90

Fidelity in biological polymerization: Key question, are we surprised?

20 of 90

The Insufficiency of Equilibrium Molecular Recognition

21 of 90

A Toy Model of Translation

22 of 90

The Kinetic Proofreading Idea: Energy to Fuel Error Correction

23 of 90

It’s Not Just About Information Amount, It’s Also About Information Density

24 of 90

Predictive Understanding of Cellular Decision Making Through the Theory-Experiment Dialogue

25 of 90

Precision Measurements to Fuel the Theory-Experiment Dialogue:�Measuring Protein Expression

26 of 90

Precision Measurements to Fuel the Theory-Experiment Dialogue:�Measuring mRNA Expression

27 of 90

Perrin’s Take on Precision Measurements and Reproducibility

28 of 90

The Meaning of Precision Measurements

  • Example of LIGO

29 of 90

Demanding Quantitative Agreement Between Measurements:�The Example of Mass Spectrometry

30 of 90

Demanding Quantitative Agreement Between Measurements:�Flow Cytometry Vs. Microscopy

31 of 90

Demanding Quantitative Agreement:�smFISH vs. Enzymatic Assays vs. Microscopy

32 of 90

Querying the Transcriptome at the Single Cell Level

33 of 90

Querying the Transcriptome at the Single Cell Level

34 of 90

The Single-Cell Sequencing Revolution

  • 44,494 cells and 10,000 genes measures
  • How do you reduce this 10,000 dimensional data to 2 dimensions?
  • How do you identify cell types?

The Tabula muris project

35 of 90

Assume That We Have a Constitutive Promoter

36 of 90

The mRNA Distribution in Space and Time

37 of 90

The Poisson Distribution Is Fully Determined by One Parameter

38 of 90

A Physical Model of the Single-Cell Sequencing Process

39 of 90

A Dishonest Coin Flip Decides Whether Each mRNA Will Be Sequenced

40 of 90

The Statistics of Coin Flips

41 of 90

The Statistics of Coin Flips

42 of 90

The Order of Coin Flips Doesn’t Matter

43 of 90

The Statistics of Coin Flips

44 of 90

The Binomial Distribution�One of the Great Probability Distributions

45 of 90

Add Savage Rosenfeld and other explanations of the Binomial distribution

46 of 90

What Happens With the mRNA Molecules That Were Not Captured?

47 of 90

A Measure of Our Precision: The Debate over Zero Inflation

48 of 90

A Challenge to Quantitative Single-Cell RNA Sequencing:�Zero Inflation and Dropout Probability

49 of 90

SANITY: Assigning Error Bars to scRNA-Seq Data

50 of 90

Querying the Transcriptome

51 of 90

Querying the Transcriptome at the Single Cell Level

52 of 90

The Single-Cell Sequencing Revolution

  • 44,494 cells and 10,000 genes measures
  • How do you reduce this 10,000 dimensional data to 2 dimensions?
  • How do you identify cell types?

The Tabula muris project

53 of 90

Cellular Decisions Are Often Driven by a Handful of Genes

54 of 90

Our Toy Model: 2D Synthetic Transcriptome

55 of 90

Identifying Cell Types in a 2D Synthetic Transcriptome

56 of 90

Finding the Right Coordinate System to Describe our Data

57 of 90

Finding the Right Coordinate System to Describe our Data

58 of 90

Key Idea: Finding the “Right” Coordinates

59 of 90

Key Idea: Finding the “Right” Coordinates

60 of 90

A Toy Model From Mechanics of the Key Idea: Finding the “Right” Coordinates For Two Coupled Oscillators

(Berman et al.)

61 of 90

Solving the Coupled Oscillators in a Bad Coordinate System

(Berman et al.)

62 of 90

Plotting The Two Coordinates Together Reveals Structure

63 of 90

Finding the “Right” Coordinates

(Berman et al.)

64 of 90

The “Right” Coordinates Reveal the Natural Variables of the System

65 of 90

Several Ways of Looking at the Problem: One from Mechanics, One as Data

(Berman et al.)

66 of 90

The Covariance Matrix of Our Rotated Data

67 of 90

Eigenvectors and Eigenvalues

68 of 90

The “Right” Coordinates Reveal the Natural Variables of the System

69 of 90

The Eigenvectors of the Covariance Matrix Define the Normal Modes (or Principal Components)

70 of 90

Your Turn: A Synthetic Transcriptome Made of Two Constitutive Promoters

71 of 90

Your Turn: Creating a Synthetic Transcriptome

72 of 90

Your Turn: Creating a Synthetic Transcriptome

73 of 90

Your Turn: Finding the Right Coordinate System to Describe our Data

74 of 90

Projecting Data Using the Dot Product

75 of 90

Finding the Right Coordinate System to Describe our Data

76 of 90

The Eigenvalues Report on the Spread of the Data Along Each Dimension

77 of 90

The Error in a Reduced Description of our System

78 of 90

The 3D Synthetic Transcriptome

79 of 90

Finding the Natural Coordinates of Our Synthetic Transcriptome

80 of 90

Dimensionality Reduction�Most of the Relevant Information Lives on a Plane

81 of 90

The Error in a Reduced Description of our System

82 of 90

Homework�Adding Downstream Genes

83 of 90

Homework�Adding Downstream Genes

84 of 90

Homework�Adding Downstream Genes

85 of 90

Homework�Adding Noisy, Uncorrelated Genes

86 of 90

A More Common Definition of PCA

87 of 90

Quantifying C. elegans movement and shape

88 of 90

The Eigenworm!

89 of 90

A Simple Synthetic Transcriptome

  • Each gene is driven by a constitutive promoter with low/high expression levels of 15/30 mRNA per cell

90 of 90

Finding Cell Types in the Transcriptome Using k-means Clustering