SRX14313189 GSM5922114_r1 SRP361807 PRJNA811146 SRS12129336 GSM5922114 GSM5922114 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313190 GSM5922115_r1 SRP361807 PRJNA811146 SRS12129335 GSM5922115 GSM5922115 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313191 GSM5922116_r1 SRP361807 PRJNA811146 SRS12129337 GSM5922116 GSM5922116 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313192 GSM5922117_r1 SRP361807 PRJNA811146 SRS12129338 GSM5922117 GSM5922117 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313193 GSM5922118_r1 SRP361807 PRJNA811146 SRS12129340 GSM5922118 GSM5922118 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313194 GSM5922119_r1 SRP361807 PRJNA811146 SRS12129339 GSM5922119 GSM5922119 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313195 GSM5922120_r1 SRP361807 PRJNA811146 SRS12129342 GSM5922120 GSM5922120 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313196 GSM5922121_r1 SRP361807 PRJNA811146 SRS12129341 GSM5922121 GSM5922121 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313197 GSM5922122_r1 SRP361807 PRJNA811146 SRS12129344 GSM5922122 GSM5922122 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313198 GSM5922123_r1 SRP361807 PRJNA811146 SRS12129343 GSM5922123 GSM5922123 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313199 GSM5922172_r1 SRP361807 PRJNA811146 SRS12129345 GSM5922172 GSM5922172 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313200 GSM5922173_r1 SRP361807 PRJNA811146 SRS12129346 GSM5922173 GSM5922173 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313201 GSM5922174_r1 SRP361807 PRJNA811146 SRS12129347 GSM5922174 GSM5922174 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313202 GSM5922175_r1 SRP361807 PRJNA811146 SRS12129348 GSM5922175 GSM5922175 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313203 GSM5922176_r1 SRP361807 PRJNA811146 SRS12129349 GSM5922176 GSM5922176 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313204 GSM5922177_r1 SRP361807 PRJNA811146 SRS12129350 GSM5922177 GSM5922177 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313205 GSM5922178_r1 SRP361807 PRJNA811146 SRS12129351 GSM5922178 GSM5922178 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313206 GSM5922179_r1 SRP361807 PRJNA811146 SRS12129352 GSM5922179 GSM5922179 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313207 GSM5922180_r1 SRP361807 PRJNA811146 SRS12129353 GSM5922180 GSM5922180 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313208 GSM5922181_r1 SRP361807 PRJNA811146 SRS12129354 GSM5922181 GSM5922181 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313209 GSM5922182_r1 SRP361807 PRJNA811146 SRS12129355 GSM5922182 GSM5922182 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313210 GSM5922183_r1 SRP361807 PRJNA811146 SRS12129356 GSM5922183 GSM5922183 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313211 GSM5922184_r1 SRP361807 PRJNA811146 SRS12129357 GSM5922184 GSM5922184 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313212 GSM5922185_r1 SRP361807 PRJNA811146 SRS12129358 GSM5922185 GSM5922185 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313213 GSM5922186_r1 SRP361807 PRJNA811146 SRS12129360 GSM5922186 GSM5922186 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313214 GSM5922187_r1 SRP361807 PRJNA811146 SRS12129359 GSM5922187 GSM5922187 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313215 GSM5922284_r1 SRP361807 PRJNA811146 SRS12129361 GSM5922284 GSM5922284 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313216 GSM5922285_r1 SRP361807 PRJNA811146 SRS12129362 GSM5922285 GSM5922285 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313217 GSM5922286_r1 SRP361807 PRJNA811146 SRS12129363 GSM5922286 GSM5922286 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313218 GSM5922287_r1 SRP361807 PRJNA811146 SRS12129364 GSM5922287 GSM5922287 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313219 GSM5922288_r1 SRP361807 PRJNA811146 SRS12129366 GSM5922288 GSM5922288 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313220 GSM5922289_r1 SRP361807 PRJNA811146 SRS12129365 GSM5922289 GSM5922289 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313221 GSM5922290_r1 SRP361807 PRJNA811146 SRS12129368 GSM5922290 GSM5922290 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313222 GSM5922291_r1 SRP361807 PRJNA811146 SRS12129367 GSM5922291 GSM5922291 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313223 GSM5922292_r1 SRP361807 PRJNA811146 SRS12129369 GSM5922292 GSM5922292 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313224 GSM5922293_r1 SRP361807 PRJNA811146 SRS12129370 GSM5922293 GSM5922293 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313225 GSM5922294_r1 SRP361807 PRJNA811146 SRS12129372 GSM5922294 GSM5922294 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313226 GSM5922295_r1 SRP361807 PRJNA811146 SRS12129371 GSM5922295 GSM5922295 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313227 GSM5922296_r1 SRP361807 PRJNA811146 SRS12129373 GSM5922296 GSM5922296 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313228 GSM5922297_r1 SRP361807 PRJNA811146 SRS12129375 GSM5922297 GSM5922297 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313229 GSM5922298_r1 SRP361807 PRJNA811146 SRS12129374 GSM5922298 GSM5922298 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313230 GSM5922299_r1 SRP361807 PRJNA811146 SRS12129376 GSM5922299 GSM5922299 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313231 GSM5922060_r1 SRP361807 PRJNA811146 SRS12129377 GSM5922060 GSM5922060 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313232 GSM5922061_r1 SRP361807 PRJNA811146 SRS12129378 GSM5922061 GSM5922061 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313233 GSM5922062_r1 SRP361807 PRJNA811146 SRS12129380 GSM5922062 GSM5922062 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313234 GSM5922063_r1 SRP361807 PRJNA811146 SRS12129379 GSM5922063 GSM5922063 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313235 GSM5922064_r1 SRP361807 PRJNA811146 SRS12129381 GSM5922064 GSM5922064 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313236 GSM5922065_r1 SRP361807 PRJNA811146 SRS12129383 GSM5922065 GSM5922065 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313237 GSM5922066_r1 SRP361807 PRJNA811146 SRS12129382 GSM5922066 GSM5922066 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313238 GSM5922067_r1 SRP361807 PRJNA811146 SRS12129384 GSM5922067 GSM5922067 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313239 GSM5922068_r1 SRP361807 PRJNA811146 SRS12129386 GSM5922068 GSM5922068 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313240 GSM5922069_r1 SRP361807 PRJNA811146 SRS12129385 GSM5922069 GSM5922069 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313241 GSM5922070_r1 SRP361807 PRJNA811146 SRS12129387 GSM5922070 GSM5922070 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313242 GSM5922071_r1 SRP361807 PRJNA811146 SRS12129389 GSM5922071 GSM5922071 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313243 GSM5922072_r1 SRP361807 PRJNA811146 SRS12129388 GSM5922072 GSM5922072 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313244 GSM5922073_r1 SRP361807 PRJNA811146 SRS12129390 GSM5922073 GSM5922073 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313245 GSM5922074_r1 SRP361807 PRJNA811146 SRS12129392 GSM5922074 GSM5922074 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313246 GSM5922075_r1 SRP361807 PRJNA811146 SRS12129391 GSM5922075 GSM5922075 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313247 GSM5922076_r1 SRP361807 PRJNA811146 SRS12129393 GSM5922076 GSM5922076 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313248 GSM5922077_r1 SRP361807 PRJNA811146 SRS12129395 GSM5922077 GSM5922077 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313249 GSM5922078_r1 SRP361807 PRJNA811146 SRS12129394 GSM5922078 GSM5922078 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313250 GSM5922079_r1 SRP361807 PRJNA811146 SRS12129396 GSM5922079 GSM5922079 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313251 GSM5922080_r1 SRP361807 PRJNA811146 SRS12129397 GSM5922080 GSM5922080 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313252 GSM5922081_r1 SRP361807 PRJNA811146 SRS12129398 GSM5922081 GSM5922081 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313253 GSM5922082_r1 SRP361807 PRJNA811146 SRS12129399 GSM5922082 GSM5922082 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313254 GSM5922083_r1 SRP361807 PRJNA811146 SRS12129400 GSM5922083 GSM5922083 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313255 GSM5922084_r1 SRP361807 PRJNA811146 SRS12129401 GSM5922084 GSM5922084 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313256 GSM5922085_r1 SRP361807 PRJNA811146 SRS12129402 GSM5922085 GSM5922085 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313257 GSM5922086_r1 SRP361807 PRJNA811146 SRS12129403 GSM5922086 GSM5922086 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313258 GSM5922087_r1 SRP361807 PRJNA811146 SRS12129404 GSM5922087 GSM5922087 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313259 GSM5922088_r1 SRP361807 PRJNA811146 SRS12129405 GSM5922088 GSM5922088 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313260 GSM5922089_r1 SRP361807 PRJNA811146 SRS12129406 GSM5922089 GSM5922089 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313261 GSM5922090_r1 SRP361807 PRJNA811146 SRS12129407 GSM5922090 GSM5922090 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313262 GSM5922091_r1 SRP361807 PRJNA811146 SRS12129408 GSM5922091 GSM5922091 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313263 GSM5922124_r1 SRP361807 PRJNA811146 SRS12129409 GSM5922124 GSM5922124 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313264 GSM5922125_r1 SRP361807 PRJNA811146 SRS12129410 GSM5922125 GSM5922125 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313265 GSM5922126_r1 SRP361807 PRJNA811146 SRS12129411 GSM5922126 GSM5922126 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313266 GSM5922127_r1 SRP361807 PRJNA811146 SRS12129412 GSM5922127 GSM5922127 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313267 GSM5922128_r1 SRP361807 PRJNA811146 SRS12129413 GSM5922128 GSM5922128 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313268 GSM5922129_r1 SRP361807 PRJNA811146 SRS12129414 GSM5922129 GSM5922129 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313269 GSM5922130_r1 SRP361807 PRJNA811146 SRS12129415 GSM5922130 GSM5922130 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313270 GSM5922131_r1 SRP361807 PRJNA811146 SRS12129416 GSM5922131 GSM5922131 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313271 GSM5922132_r1 SRP361807 PRJNA811146 SRS12129417 GSM5922132 GSM5922132 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313272 GSM5922133_r1 SRP361807 PRJNA811146 SRS12129418 GSM5922133 GSM5922133 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313273 GSM5922134_r1 SRP361807 PRJNA811146 SRS12129420 GSM5922134 GSM5922134 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313274 GSM5922135_r1 SRP361807 PRJNA811146 SRS12129419 GSM5922135 GSM5922135 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313275 GSM5922136_r1 SRP361807 PRJNA811146 SRS12129421 GSM5922136 GSM5922136 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313276 GSM5922137_r1 SRP361807 PRJNA811146 SRS12129424 GSM5922137 GSM5922137 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313277 GSM5922138_r1 SRP361807 PRJNA811146 SRS12129422 GSM5922138 GSM5922138 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313278 GSM5922139_r1 SRP361807 PRJNA811146 SRS12129423 GSM5922139 GSM5922139 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313279 GSM5922188_r1 SRP361807 PRJNA811146 SRS12129425 GSM5922188 GSM5922188 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313280 GSM5922189_r1 SRP361807 PRJNA811146 SRS12129426 GSM5922189 GSM5922189 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313281 GSM5922190_r1 SRP361807 PRJNA811146 SRS12129427 GSM5922190 GSM5922190 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313282 GSM5922191_r1 SRP361807 PRJNA811146 SRS12129428 GSM5922191 GSM5922191 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313283 GSM5922192_r1 SRP361807 PRJNA811146 SRS12129429 GSM5922192 GSM5922192 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313284 GSM5922193_r1 SRP361807 PRJNA811146 SRS12129430 GSM5922193 GSM5922193 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313285 GSM5922194_r1 SRP361807 PRJNA811146 SRS12129431 GSM5922194 GSM5922194 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313286 GSM5922195_r1 SRP361807 PRJNA811146 SRS12129432 GSM5922195 GSM5922195 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313287 GSM5922196_r1 SRP361807 PRJNA811146 SRS12129433 GSM5922196 GSM5922196 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313288 GSM5922197_r1 SRP361807 PRJNA811146 SRS12129434 GSM5922197 GSM5922197 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313289 GSM5922198_r1 SRP361807 PRJNA811146 SRS12129435 GSM5922198 GSM5922198 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313290 GSM5922199_r1 SRP361807 PRJNA811146 SRS12129436 GSM5922199 GSM5922199 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313291 GSM5922200_r1 SRP361807 PRJNA811146 SRS12129437 GSM5922200 GSM5922200 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313292 GSM5922201_r1 SRP361807 PRJNA811146 SRS12129438 GSM5922201 GSM5922201 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313293 GSM5922202_r1 SRP361807 PRJNA811146 SRS12129439 GSM5922202 GSM5922202 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313294 GSM5922203_r1 SRP361807 PRJNA811146 SRS12129440 GSM5922203 GSM5922203 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313295 GSM5922236_r1 SRP361807 PRJNA811146 SRS12129441 GSM5922236 GSM5922236 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313296 GSM5922237_r1 SRP361807 PRJNA811146 SRS12129442 GSM5922237 GSM5922237 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313297 GSM5922238_r1 SRP361807 PRJNA811146 SRS12129443 GSM5922238 GSM5922238 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313298 GSM5922239_r1 SRP361807 PRJNA811146 SRS12129444 GSM5922239 GSM5922239 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313299 GSM5922240_r1 SRP361807 PRJNA811146 SRS12129445 GSM5922240 GSM5922240 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313300 GSM5922241_r1 SRP361807 PRJNA811146 SRS12129447 GSM5922241 GSM5922241 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313301 GSM5922242_r1 SRP361807 PRJNA811146 SRS12129446 GSM5922242 GSM5922242 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313302 GSM5922243_r1 SRP361807 PRJNA811146 SRS12129448 GSM5922243 GSM5922243 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313303 GSM5922244_r1 SRP361807 PRJNA811146 SRS12129450 GSM5922244 GSM5922244 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313304 GSM5922245_r1 SRP361807 PRJNA811146 SRS12129449 GSM5922245 GSM5922245 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313305 GSM5922246_r1 SRP361807 PRJNA811146 SRS12129451 GSM5922246 GSM5922246 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313306 GSM5922247_r1 SRP361807 PRJNA811146 SRS12129453 GSM5922247 GSM5922247 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313307 GSM5922248_r1 SRP361807 PRJNA811146 SRS12129452 GSM5922248 GSM5922248 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313308 GSM5922249_r1 SRP361807 PRJNA811146 SRS12129454 GSM5922249 GSM5922249 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313309 GSM5922250_r1 SRP361807 PRJNA811146 SRS12129455 GSM5922250 GSM5922250 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313310 GSM5922251_r1 SRP361807 PRJNA811146 SRS12129457 GSM5922251 GSM5922251 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313311 GSM5922300_r1 SRP361807 PRJNA811146 SRS12129456 GSM5922300 GSM5922300 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313312 GSM5922301_r1 SRP361807 PRJNA811146 SRS12129458 GSM5922301 GSM5922301 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313313 GSM5922302_r1 SRP361807 PRJNA811146 SRS12129459 GSM5922302 GSM5922302 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313314 GSM5922303_r1 SRP361807 PRJNA811146 SRS12129460 GSM5922303 GSM5922303 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313315 GSM5922304_r1 SRP361807 PRJNA811146 SRS12129462 GSM5922304 GSM5922304 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313316 GSM5922305_r1 SRP361807 PRJNA811146 SRS12129461 GSM5922305 GSM5922305 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313317 GSM5922306_r1 SRP361807 PRJNA811146 SRS12129463 GSM5922306 GSM5922306 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313318 GSM5922307_r1 SRP361807 PRJNA811146 SRS12129466 GSM5922307 GSM5922307 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313319 GSM5922308_r1 SRP361807 PRJNA811146 SRS12129464 GSM5922308 GSM5922308 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313320 GSM5922309_r1 SRP361807 PRJNA811146 SRS12129465 GSM5922309 GSM5922309 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313321 GSM5922310_r1 SRP361807 PRJNA811146 SRS12129468 GSM5922310 GSM5922310 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313322 GSM5922311_r1 SRP361807 PRJNA811146 SRS12129467 GSM5922311 GSM5922311 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313323 GSM5922312_r1 SRP361807 PRJNA811146 SRS12129469 GSM5922312 GSM5922312 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313324 GSM5922313_r1 SRP361807 PRJNA811146 SRS12129470 GSM5922313 GSM5922313 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313325 GSM5922314_r1 SRP361807 PRJNA811146 SRS12129472 GSM5922314 GSM5922314 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313326 GSM5922315_r1 SRP361807 PRJNA811146 SRS12129471 GSM5922315 GSM5922315 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313327 GSM5922028_r1 SRP361807 PRJNA811146 SRS12129475 GSM5922028 GSM5922028 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313328 GSM5922029_r1 SRP361807 PRJNA811146 SRS12129474 GSM5922029 GSM5922029 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313329 GSM5922030_r1 SRP361807 PRJNA811146 SRS12129473 GSM5922030 GSM5922030 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313330 GSM5922031_r1 SRP361807 PRJNA811146 SRS12129477 GSM5922031 GSM5922031 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313331 GSM5922032_r1 SRP361807 PRJNA811146 SRS12129476 GSM5922032 GSM5922032 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313332 GSM5922033_r1 SRP361807 PRJNA811146 SRS12129478 GSM5922033 GSM5922033 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313333 GSM5922034_r1 SRP361807 PRJNA811146 SRS12129479 GSM5922034 GSM5922034 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313334 GSM5922035_r1 SRP361807 PRJNA811146 SRS12129480 GSM5922035 GSM5922035 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313335 GSM5922036_r1 SRP361807 PRJNA811146 SRS12129481 GSM5922036 GSM5922036 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313336 GSM5922037_r1 SRP361807 PRJNA811146 SRS12129484 GSM5922037 GSM5922037 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313337 GSM5922038_r1 SRP361807 PRJNA811146 SRS12129482 GSM5922038 GSM5922038 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313338 GSM5922039_r1 SRP361807 PRJNA811146 SRS12129483 GSM5922039 GSM5922039 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313339 GSM5922040_r1 SRP361807 PRJNA811146 SRS12129485 GSM5922040 GSM5922040 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313340 GSM5922041_r1 SRP361807 PRJNA811146 SRS12129486 GSM5922041 GSM5922041 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313341 GSM5922042_r1 SRP361807 PRJNA811146 SRS12129487 GSM5922042 GSM5922042 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313342 GSM5922043_r1 SRP361807 PRJNA811146 SRS12129489 GSM5922043 GSM5922043 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313343 GSM5922332_r1 SRP361807 PRJNA811146 SRS12129488 GSM5922332 GSM5922332 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313344 GSM5922333_r1 SRP361807 PRJNA811146 SRS12129490 GSM5922333 GSM5922333 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313345 GSM5922334_r1 SRP361807 PRJNA811146 SRS12129491 GSM5922334 GSM5922334 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313346 GSM5922335_r1 SRP361807 PRJNA811146 SRS12129492 GSM5922335 GSM5922335 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313347 GSM5922336_r1 SRP361807 PRJNA811146 SRS12129495 GSM5922336 GSM5922336 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313348 GSM5922337_r1 SRP361807 PRJNA811146 SRS12129493 GSM5922337 GSM5922337 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313349 GSM5922140_r1 SRP361807 PRJNA811146 SRS12129494 GSM5922140 GSM5922140 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313350 GSM5922141_r1 SRP361807 PRJNA811146 SRS12129496 GSM5922141 GSM5922141 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313351 GSM5922142_r1 SRP361807 PRJNA811146 SRS12129498 GSM5922142 GSM5922142 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313352 GSM5922143_r1 SRP361807 PRJNA811146 SRS12129497 GSM5922143 GSM5922143 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313353 GSM5922144_r1 SRP361807 PRJNA811146 SRS12129499 GSM5922144 GSM5922144 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313354 GSM5922145_r1 SRP361807 PRJNA811146 SRS12129500 GSM5922145 GSM5922145 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313355 GSM5922146_r1 SRP361807 PRJNA811146 SRS12129501 GSM5922146 GSM5922146 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313356 GSM5922147_r1 SRP361807 PRJNA811146 SRS12129502 GSM5922147 GSM5922147 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313357 GSM5922148_r1 SRP361807 PRJNA811146 SRS12129503 GSM5922148 GSM5922148 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313358 GSM5922149_r1 SRP361807 PRJNA811146 SRS12129504 GSM5922149 GSM5922149 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313359 GSM5922150_r1 SRP361807 PRJNA811146 SRS12129505 GSM5922150 GSM5922150 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313360 GSM5922151_r1 SRP361807 PRJNA811146 SRS12129506 GSM5922151 GSM5922151 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313361 GSM5922152_r1 SRP361807 PRJNA811146 SRS12129508 GSM5922152 GSM5922152 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313362 GSM5922153_r1 SRP361807 PRJNA811146 SRS12129507 GSM5922153 GSM5922153 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313363 GSM5922154_r1 SRP361807 PRJNA811146 SRS12129509 GSM5922154 GSM5922154 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313364 GSM5922155_r1 SRP361807 PRJNA811146 SRS12129510 GSM5922155 GSM5922155 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313365 GSM5922204_r1 SRP361807 PRJNA811146 SRS12129511 GSM5922204 GSM5922204 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313366 GSM5922205_r1 SRP361807 PRJNA811146 SRS12129512 GSM5922205 GSM5922205 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313367 GSM5922206_r1 SRP361807 PRJNA811146 SRS12129513 GSM5922206 GSM5922206 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313368 GSM5922207_r1 SRP361807 PRJNA811146 SRS12129514 GSM5922207 GSM5922207 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313369 GSM5922208_r1 SRP361807 PRJNA811146 SRS12129515 GSM5922208 GSM5922208 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313370 GSM5922209_r1 SRP361807 PRJNA811146 SRS12129516 GSM5922209 GSM5922209 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313371 GSM5922210_r1 SRP361807 PRJNA811146 SRS12129517 GSM5922210 GSM5922210 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313372 GSM5922211_r1 SRP361807 PRJNA811146 SRS12129518 GSM5922211 GSM5922211 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313373 GSM5922212_r1 SRP361807 PRJNA811146 SRS12129520 GSM5922212 GSM5922212 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313374 GSM5922213_r1 SRP361807 PRJNA811146 SRS12129519 GSM5922213 GSM5922213 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313375 GSM5922214_r1 SRP361807 PRJNA811146 SRS12129521 GSM5922214 GSM5922214 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313376 GSM5922215_r1 SRP361807 PRJNA811146 SRS12129523 GSM5922215 GSM5922215 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313377 GSM5922216_r1 SRP361807 PRJNA811146 SRS12129522 GSM5922216 GSM5922216 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313378 GSM5922217_r1 SRP361807 PRJNA811146 SRS12129524 GSM5922217 GSM5922217 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313379 GSM5922218_r1 SRP361807 PRJNA811146 SRS12129526 GSM5922218 GSM5922218 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313380 GSM5922219_r1 SRP361807 PRJNA811146 SRS12129525 GSM5922219 GSM5922219 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313381 GSM5922220_r1 SRP361807 PRJNA811146 SRS12129527 GSM5922220 GSM5922220 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313382 GSM5922221_r1 SRP361807 PRJNA811146 SRS12129528 GSM5922221 GSM5922221 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313383 GSM5922222_r1 SRP361807 PRJNA811146 SRS12129529 GSM5922222 GSM5922222 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313384 GSM5922223_r1 SRP361807 PRJNA811146 SRS12129530 GSM5922223 GSM5922223 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313385 GSM5922224_r1 SRP361807 PRJNA811146 SRS12129531 GSM5922224 GSM5922224 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313386 GSM5922225_r1 SRP361807 PRJNA811146 SRS12129532 GSM5922225 GSM5922225 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313387 GSM5922226_r1 SRP361807 PRJNA811146 SRS12129533 GSM5922226 GSM5922226 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313388 GSM5922227_r1 SRP361807 PRJNA811146 SRS12129534 GSM5922227 GSM5922227 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313389 GSM5922228_r1 SRP361807 PRJNA811146 SRS12129535 GSM5922228 GSM5922228 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313390 GSM5922229_r1 SRP361807 PRJNA811146 SRS12129536 GSM5922229 GSM5922229 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313391 GSM5922230_r1 SRP361807 PRJNA811146 SRS12129538 GSM5922230 GSM5922230 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313392 GSM5922231_r1 SRP361807 PRJNA811146 SRS12129537 GSM5922231 GSM5922231 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313393 GSM5922232_r1 SRP361807 PRJNA811146 SRS12129539 GSM5922232 GSM5922232 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313394 GSM5922233_r1 SRP361807 PRJNA811146 SRS12129541 GSM5922233 GSM5922233 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313395 GSM5922234_r1 SRP361807 PRJNA811146 SRS12129540 GSM5922234 GSM5922234 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313396 GSM5922235_r1 SRP361807 PRJNA811146 SRS12129542 GSM5922235 GSM5922235 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313397 GSM5922252_r1 SRP361807 PRJNA811146 SRS12129544 GSM5922252 GSM5922252 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313398 GSM5922253_r1 SRP361807 PRJNA811146 SRS12129543 GSM5922253 GSM5922253 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313399 GSM5922254_r1 SRP361807 PRJNA811146 SRS12129545 GSM5922254 GSM5922254 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313400 GSM5922255_r1 SRP361807 PRJNA811146 SRS12129546 GSM5922255 GSM5922255 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313401 GSM5922256_r1 SRP361807 PRJNA811146 SRS12129547 GSM5922256 GSM5922256 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313402 GSM5922257_r1 SRP361807 PRJNA811146 SRS12129548 GSM5922257 GSM5922257 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313403 GSM5922258_r1 SRP361807 PRJNA811146 SRS12129549 GSM5922258 GSM5922258 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313404 GSM5922259_r1 SRP361807 PRJNA811146 SRS12129550 GSM5922259 GSM5922259 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313405 GSM5922260_r1 SRP361807 PRJNA811146 SRS12129551 GSM5922260 GSM5922260 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313406 GSM5922261_r1 SRP361807 PRJNA811146 SRS12129552 GSM5922261 GSM5922261 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313407 GSM5922262_r1 SRP361807 PRJNA811146 SRS12129553 GSM5922262 GSM5922262 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313408 GSM5922263_r1 SRP361807 PRJNA811146 SRS12129554 GSM5922263 GSM5922263 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313409 GSM5922264_r1 SRP361807 PRJNA811146 SRS12129555 GSM5922264 GSM5922264 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313410 GSM5922265_r1 SRP361807 PRJNA811146 SRS12129556 GSM5922265 GSM5922265 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313411 GSM5922266_r1 SRP361807 PRJNA811146 SRS12129557 GSM5922266 GSM5922266 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313412 GSM5922267_r1 SRP361807 PRJNA811146 SRS12129558 GSM5922267 GSM5922267 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313413 GSM5922316_r1 SRP361807 PRJNA811146 SRS12129559 GSM5922316 GSM5922316 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313414 GSM5922317_r1 SRP361807 PRJNA811146 SRS12129560 GSM5922317 GSM5922317 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313415 GSM5922318_r1 SRP361807 PRJNA811146 SRS12129561 GSM5922318 GSM5922318 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313416 GSM5922319_r1 SRP361807 PRJNA811146 SRS12129563 GSM5922319 GSM5922319 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313417 GSM5922320_r1 SRP361807 PRJNA811146 SRS12129562 GSM5922320 GSM5922320 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313418 GSM5922321_r1 SRP361807 PRJNA811146 SRS12129564 GSM5922321 GSM5922321 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313419 GSM5922322_r1 SRP361807 PRJNA811146 SRS12129565 GSM5922322 GSM5922322 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313420 GSM5922323_r1 SRP361807 PRJNA811146 SRS12129566 GSM5922323 GSM5922323 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313421 GSM5922324_r1 SRP361807 PRJNA811146 SRS12129567 GSM5922324 GSM5922324 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313422 GSM5922325_r1 SRP361807 PRJNA811146 SRS12129569 GSM5922325 GSM5922325 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313423 GSM5922326_r1 SRP361807 PRJNA811146 SRS12129568 GSM5922326 GSM5922326 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313424 GSM5922327_r1 SRP361807 PRJNA811146 SRS12129570 GSM5922327 GSM5922327 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313425 GSM5922328_r1 SRP361807 PRJNA811146 SRS12129572 GSM5922328 GSM5922328 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313426 GSM5922329_r1 SRP361807 PRJNA811146 SRS12129571 GSM5922329 GSM5922329 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313427 GSM5922330_r1 SRP361807 PRJNA811146 SRS12129573 GSM5922330 GSM5922330 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313428 GSM5922331_r1 SRP361807 PRJNA811146 SRS12129574 GSM5922331 GSM5922331 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313429 GSM5922044_r1 SRP361807 PRJNA811146 SRS12129575 GSM5922044 GSM5922044 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313430 GSM5922045_r1 SRP361807 PRJNA811146 SRS12129576 GSM5922045 GSM5922045 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313431 GSM5922046_r1 SRP361807 PRJNA811146 SRS12129577 GSM5922046 GSM5922046 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313432 GSM5922047_r1 SRP361807 PRJNA811146 SRS12129578 GSM5922047 GSM5922047 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313433 GSM5922048_r1 SRP361807 PRJNA811146 SRS12129579 GSM5922048 GSM5922048 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313434 GSM5922049_r1 SRP361807 PRJNA811146 SRS12129580 GSM5922049 GSM5922049 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313435 GSM5922050_r1 SRP361807 PRJNA811146 SRS12129581 GSM5922050 GSM5922050 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313436 GSM5922051_r1 SRP361807 PRJNA811146 SRS12129582 GSM5922051 GSM5922051 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313437 GSM5922052_r1 SRP361807 PRJNA811146 SRS12129583 GSM5922052 GSM5922052 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313438 GSM5922053_r1 SRP361807 PRJNA811146 SRS12129584 GSM5922053 GSM5922053 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313439 GSM5922054_r1 SRP361807 PRJNA811146 SRS12129585 GSM5922054 GSM5922054 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313440 GSM5922055_r1 SRP361807 PRJNA811146 SRS12129587 GSM5922055 GSM5922055 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313441 GSM5922056_r1 SRP361807 PRJNA811146 SRS12129586 GSM5922056 GSM5922056 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313442 GSM5922057_r1 SRP361807 PRJNA811146 SRS12129588 GSM5922057 GSM5922057 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313443 GSM5922058_r1 SRP361807 PRJNA811146 SRS12129589 GSM5922058 GSM5922058 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313444 GSM5922059_r1 SRP361807 PRJNA811146 SRS12129590 GSM5922059 GSM5922059 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313445 GSM5922092_r1 SRP361807 PRJNA811146 SRS12129591 GSM5922092 GSM5922092 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313446 GSM5922093_r1 SRP361807 PRJNA811146 SRS12129594 GSM5922093 GSM5922093 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313447 GSM5922094_r1 SRP361807 PRJNA811146 SRS12129593 GSM5922094 GSM5922094 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313448 GSM5922095_r1 SRP361807 PRJNA811146 SRS12129592 GSM5922095 GSM5922095 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313449 GSM5922096_r1 SRP361807 PRJNA811146 SRS12129595 GSM5922096 GSM5922096 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313450 GSM5922097_r1 SRP361807 PRJNA811146 SRS12129596 GSM5922097 GSM5922097 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313451 GSM5922098_r1 SRP361807 PRJNA811146 SRS12129597 GSM5922098 GSM5922098 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313452 GSM5922099_r1 SRP361807 PRJNA811146 SRS12129599 GSM5922099 GSM5922099 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313453 GSM5922100_r1 SRP361807 PRJNA811146 SRS12129598 GSM5922100 GSM5922100 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313454 GSM5922101_r1 SRP361807 PRJNA811146 SRS12129600 GSM5922101 GSM5922101 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313455 GSM5922102_r1 SRP361807 PRJNA811146 SRS12129601 GSM5922102 GSM5922102 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313456 GSM5922103_r1 SRP361807 PRJNA811146 SRS12129602 GSM5922103 GSM5922103 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313457 GSM5922104_r1 SRP361807 PRJNA811146 SRS12129603 GSM5922104 GSM5922104 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313458 GSM5922105_r1 SRP361807 PRJNA811146 SRS12129604 GSM5922105 GSM5922105 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313459 GSM5922106_r1 SRP361807 PRJNA811146 SRS12129605 GSM5922106 GSM5922106 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313460 GSM5922107_r1 SRP361807 PRJNA811146 SRS12129606 GSM5922107 GSM5922107 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313461 GSM5922156_r1 SRP361807 PRJNA811146 SRS12129607 GSM5922156 GSM5922156 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313462 GSM5922157_r1 SRP361807 PRJNA811146 SRS12129608 GSM5922157 GSM5922157 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313463 GSM5922158_r1 SRP361807 PRJNA811146 SRS12129609 GSM5922158 GSM5922158 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313464 GSM5922159_r1 SRP361807 PRJNA811146 SRS12129610 GSM5922159 GSM5922159 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313465 GSM5922160_r1 SRP361807 PRJNA811146 SRS12129611 GSM5922160 GSM5922160 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313466 GSM5922161_r1 SRP361807 PRJNA811146 SRS12129612 GSM5922161 GSM5922161 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313467 GSM5922162_r1 SRP361807 PRJNA811146 SRS12129613 GSM5922162 GSM5922162 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313468 GSM5922163_r1 SRP361807 PRJNA811146 SRS12129615 GSM5922163 GSM5922163 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313469 GSM5922164_r1 SRP361807 PRJNA811146 SRS12129614 GSM5922164 GSM5922164 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313470 GSM5922165_r1 SRP361807 PRJNA811146 SRS12129616 GSM5922165 GSM5922165 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313471 GSM5922166_r1 SRP361807 PRJNA811146 SRS12129618 GSM5922166 GSM5922166 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313472 GSM5922167_r1 SRP361807 PRJNA811146 SRS12129617 GSM5922167 GSM5922167 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313473 GSM5922168_r1 SRP361807 PRJNA811146 SRS12129619 GSM5922168 GSM5922168 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313474 GSM5922169_r1 SRP361807 PRJNA811146 SRS12129620 GSM5922169 GSM5922169 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313475 GSM5922170_r1 SRP361807 PRJNA811146 SRS12129621 GSM5922170 GSM5922170 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313476 GSM5922171_r1 SRP361807 PRJNA811146 SRS12129622 GSM5922171 GSM5922171 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313477 GSM5922268_r1 SRP361807 PRJNA811146 SRS12129623 GSM5922268 GSM5922268 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313478 GSM5922269_r1 SRP361807 PRJNA811146 SRS12129624 GSM5922269 GSM5922269 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313479 GSM5922270_r1 SRP361807 PRJNA811146 SRS12129625 GSM5922270 GSM5922270 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313480 GSM5922271_r1 SRP361807 PRJNA811146 SRS12129626 GSM5922271 GSM5922271 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313481 GSM5922272_r1 SRP361807 PRJNA811146 SRS12129627 GSM5922272 GSM5922272 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313482 GSM5922273_r1 SRP361807 PRJNA811146 SRS12129628 GSM5922273 GSM5922273 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313483 GSM5922274_r1 SRP361807 PRJNA811146 SRS12129629 GSM5922274 GSM5922274 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313484 GSM5922275_r1 SRP361807 PRJNA811146 SRS12129630 GSM5922275 GSM5922275 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313485 GSM5922276_r1 SRP361807 PRJNA811146 SRS12129631 GSM5922276 GSM5922276 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313486 GSM5922277_r1 SRP361807 PRJNA811146 SRS12129632 GSM5922277 GSM5922277 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313487 GSM5922278_r1 SRP361807 PRJNA811146 SRS12129633 GSM5922278 GSM5922278 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313488 GSM5922279_r1 SRP361807 PRJNA811146 SRS12129634 GSM5922279 GSM5922279 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313489 GSM5922280_r1 SRP361807 PRJNA811146 SRS12129635 GSM5922280 GSM5922280 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313490 GSM5922281_r1 SRP361807 PRJNA811146 SRS12129636 GSM5922281 GSM5922281 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313491 GSM5922282_r1 SRP361807 PRJNA811146 SRS12129637 GSM5922282 GSM5922282 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313492 GSM5922283_r1 SRP361807 PRJNA811146 SRS12129638 GSM5922283 GSM5922283 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313493 GSM5922108_r1 SRP361807 PRJNA811146 SRS12129639 GSM5922108 GSM5922108 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313494 GSM5922109_r1 SRP361807 PRJNA811146 SRS12129640 GSM5922109 GSM5922109 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313495 GSM5922110_r1 SRP361807 PRJNA811146 SRS12129641 GSM5922110 GSM5922110 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313496 GSM5922111_r1 SRP361807 PRJNA811146 SRS12129642 GSM5922111 GSM5922111 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313497 GSM5922112_r1 SRP361807 PRJNA811146 SRS12129643 GSM5922112 GSM5922112 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500 SRX14313498 GSM5922113_r1 SRP361807 PRJNA811146 SRS12129644 GSM5922113 GSM5922113 RNA-Seq TRANSCRIPTOMIC cDNA Flow citometry was used to sort tomato positive epidermal cells Full length single cell RNA-seq was performed using a modified version of the Smart-seq2 protocol. Briefly, individual cells are sorted into 96 well plates containing lysis buffer in presence of RNase inhibitor, dNTPs and OligodT. Reverse transcription of the polyadenylated RNA will be performed with SuperScriptII and Template Switching Oligos. The resulting cDNA will be amplified with 25 cycles of PCR and libraries will be prepared for sequencing with standard NexteraXT Illumina protocol. Libraries were sequenced on Illumina NextSeq 500 System (single-end 75bp reads). NextSeq 500