Organoids derived from two independent donors were microinjected with oocyst or control media. Total RNA was extracted from collected organoids at 24 and 72 hr post-injection using RNeasy Micro Kit (Qiagen) according to manufacturer’s protocol including DNaseI treatment. Quality and quantity of isolated RNA was checked and measured with Bioanalyzer 2100 RNA Nano 6000 chips (Agilent). Libraries were prepared as described previously47 (link) and sequenced on an Illumina NextSeq500 by using 75-bp paired-end sequencing. Paired-end reads from Illumina sequencing were aligned to the human transcriptome genome and C. parvum transcriptome genome (Iowa strain) by BWA48 (link). DeSeq (v1.18.0) was used for read normalization and differential expression analysis. Gene set enrichment analysis (GSEA) was performed using gene lists for type I interferon response and regulation against normalized RNA-seq reads of injected SI and lung organoids using GSEA software v3.0 beta2.