library_construction_protocol: RNA was extracted from ~1.5 million cells using the Zymo dual DNA/RNA extraction kit (D7001, Zymo, Irvine, CA, USA). All treatments and timepoints per individual were extracted in one batch resulting in 6 extraction batches of 12 samples corresponding to the 6 individuals. Median RIN scores were 9.35, 9.6, 9.3, 9.6, 9.65, and 9.55 across treatments for individuals 1-6, respectively. RNA-seq libraries were generated using 250 ng of RNA using the NEBNext® Poly(A) mRNA Magnetic Isolation Module kit (E7490L, Ipswich, MA, USA), NEBNext Ultra II RNA library prep with Sample Purification Beads kit (E7775K) and the NEBNext® Multiplex Oligos for Illumina® (96 Unique Dual Index Primer Pairs) kit (E6440S). Libraries were prepared in batches of 12 for each individual, following the manufacturer's protocol.