SRX3026717 GSM2714282 SRP113298 SRS2375808 GSM2714282 RNA-Seq TRANSCRIPTOMIC cDNA Worker and drone embryos were collected from the distinct type of frame laid by the same queen. We sampled two groups of worker and drone embryos laid by two different honeybee queens. Then a total of 36 embryos were collected for single-embryo RNA-seq. High-quality, full-length cDNA is generated directly from 1–1,000 cells by The SMART-Seq v4 Ultra Low Input RNA Kit for Sequencing (Cat. Nos. 634890). cDNA was amplified by LD PCR. RNA can be transcribed by T7 promoter from ds cDNA , RNA was treated with RQ1 DNase (promega) to remove DNA .The quantity of the purified DNA were determined by Qubit. For each sample, 200 ng RNA was used for RNA-seq library preparation. RNAs were iron fragmented at 95℃ followed by end repair and 5' adaptor ligation. Then reverse transcription was performed with RT primer harboring 3' adaptor sequence and randomized hexamer. The cDNAs were purified and amplified and PCR products corresponding to 200-500 bps were purified, quantified and stored at -80 ℃ until used for sequencing. NextSeq 500 gds 302714282 GSM2714282 GEO Accession GSM2714282 SRX3026718 GSM2714283 SRP113298 SRS2375809 GSM2714283 RNA-Seq TRANSCRIPTOMIC cDNA Worker and drone embryos were collected from the distinct type of frame laid by the same queen. We sampled two groups of worker and drone embryos laid by two different honeybee queens. Then a total of 36 embryos were collected for single-embryo RNA-seq. High-quality, full-length cDNA is generated directly from 1–1,000 cells by The SMART-Seq v4 Ultra Low Input RNA Kit for Sequencing (Cat. Nos. 634890). cDNA was amplified by LD PCR. RNA can be transcribed by T7 promoter from ds cDNA , RNA was treated with RQ1 DNase (promega) to remove DNA .The quantity of the purified DNA were determined by Qubit. For each sample, 200 ng RNA was used for RNA-seq library preparation. RNAs were iron fragmented at 95℃ followed by end repair and 5' adaptor ligation. Then reverse transcription was performed with RT primer harboring 3' adaptor sequence and randomized hexamer. The cDNAs were purified and amplified and PCR products corresponding to 200-500 bps were purified, quantified and stored at -80 ℃ until used for sequencing. NextSeq 500 gds 302714283 GSM2714283 GEO Accession GSM2714283 SRX3026719 GSM2714284 SRP113298 SRS2375810 GSM2714284 RNA-Seq TRANSCRIPTOMIC cDNA Worker and drone embryos were collected from the distinct type of frame laid by the same queen. We sampled two groups of worker and drone embryos laid by two different honeybee queens. Then a total of 36 embryos were collected for single-embryo RNA-seq. High-quality, full-length cDNA is generated directly from 1–1,000 cells by The SMART-Seq v4 Ultra Low Input RNA Kit for Sequencing (Cat. Nos. 634890). cDNA was amplified by LD PCR. RNA can be transcribed by T7 promoter from ds cDNA , RNA was treated with RQ1 DNase (promega) to remove DNA .The quantity of the purified DNA were determined by Qubit. For each sample, 200 ng RNA was used for RNA-seq library preparation. RNAs were iron fragmented at 95℃ followed by end repair and 5' adaptor ligation. Then reverse transcription was performed with RT primer harboring 3' adaptor sequence and randomized hexamer. The cDNAs were purified and amplified and PCR products corresponding to 200-500 bps were purified, quantified and stored at -80 ℃ until used for sequencing. NextSeq 500 gds 302714284 GSM2714284 GEO Accession GSM2714284 SRX3026720 GSM2714285 SRP113298 SRS2375811 GSM2714285 RNA-Seq TRANSCRIPTOMIC cDNA Worker and drone embryos were collected from the distinct type of frame laid by the same queen. We sampled two groups of worker and drone embryos laid by two different honeybee queens. Then a total of 36 embryos were collected for single-embryo RNA-seq. High-quality, full-length cDNA is generated directly from 1–1,000 cells by The SMART-Seq v4 Ultra Low Input RNA Kit for Sequencing (Cat. Nos. 634890). cDNA was amplified by LD PCR. RNA can be transcribed by T7 promoter from ds cDNA , RNA was treated with RQ1 DNase (promega) to remove DNA .The quantity of the purified DNA were determined by Qubit. For each sample, 200 ng RNA was used for RNA-seq library preparation. RNAs were iron fragmented at 95℃ followed by end repair and 5' adaptor ligation. Then reverse transcription was performed with RT primer harboring 3' adaptor sequence and randomized hexamer. The cDNAs were purified and amplified and PCR products corresponding to 200-500 bps were purified, quantified and stored at -80 ℃ until used for sequencing. NextSeq 500 gds 302714285 GSM2714285 GEO Accession GSM2714285 SRX3026721 GSM2714286 SRP113298 SRS2375812 GSM2714286 RNA-Seq TRANSCRIPTOMIC cDNA Worker and drone embryos were collected from the distinct type of frame laid by the same queen. We sampled two groups of worker and drone embryos laid by two different honeybee queens. Then a total of 36 embryos were collected for single-embryo RNA-seq. High-quality, full-length cDNA is generated directly from 1–1,000 cells by The SMART-Seq v4 Ultra Low Input RNA Kit for Sequencing (Cat. Nos. 634890). cDNA was amplified by LD PCR. RNA can be transcribed by T7 promoter from ds cDNA , RNA was treated with RQ1 DNase (promega) to remove DNA .The quantity of the purified DNA were determined by Qubit. For each sample, 200 ng RNA was used for RNA-seq library preparation. RNAs were iron fragmented at 95℃ followed by end repair and 5' adaptor ligation. Then reverse transcription was performed with RT primer harboring 3' adaptor sequence and randomized hexamer. The cDNAs were purified and amplified and PCR products corresponding to 200-500 bps were purified, quantified and stored at -80 ℃ until used for sequencing. NextSeq 500 gds 302714286 GSM2714286 GEO Accession GSM2714286 SRX3026722 GSM2714287 SRP113298 SRS2375813 GSM2714287 RNA-Seq TRANSCRIPTOMIC cDNA Worker and drone embryos were collected from the distinct type of frame laid by the same queen. We sampled two groups of worker and drone embryos laid by two different honeybee queens. Then a total of 36 embryos were collected for single-embryo RNA-seq. High-quality, full-length cDNA is generated directly from 1–1,000 cells by The SMART-Seq v4 Ultra Low Input RNA Kit for Sequencing (Cat. Nos. 634890). cDNA was amplified by LD PCR. RNA can be transcribed by T7 promoter from ds cDNA , RNA was treated with RQ1 DNase (promega) to remove DNA .The quantity of the purified DNA were determined by Qubit. For each sample, 200 ng RNA was used for RNA-seq library preparation. RNAs were iron fragmented at 95℃ followed by end repair and 5' adaptor ligation. Then reverse transcription was performed with RT primer harboring 3' adaptor sequence and randomized hexamer. The cDNAs were purified and amplified and PCR products corresponding to 200-500 bps were purified, quantified and stored at -80 ℃ until used for sequencing. NextSeq 500 gds 302714287 GSM2714287 GEO Accession GSM2714287 SRX3026723 GSM2714288 SRP113298 SRS2375814 GSM2714288 RNA-Seq TRANSCRIPTOMIC cDNA Worker and drone embryos were collected from the distinct type of frame laid by the same queen. We sampled two groups of worker and drone embryos laid by two different honeybee queens. Then a total of 36 embryos were collected for single-embryo RNA-seq. High-quality, full-length cDNA is generated directly from 1–1,000 cells by The SMART-Seq v4 Ultra Low Input RNA Kit for Sequencing (Cat. Nos. 634890). cDNA was amplified by LD PCR. RNA can be transcribed by T7 promoter from ds cDNA , RNA was treated with RQ1 DNase (promega) to remove DNA .The quantity of the purified DNA were determined by Qubit. For each sample, 200 ng RNA was used for RNA-seq library preparation. RNAs were iron fragmented at 95℃ followed by end repair and 5' adaptor ligation. Then reverse transcription was performed with RT primer harboring 3' adaptor sequence and randomized hexamer. The cDNAs were purified and amplified and PCR products corresponding to 200-500 bps were purified, quantified and stored at -80 ℃ until used for sequencing. NextSeq 500 gds 302714288 GSM2714288 GEO Accession GSM2714288 SRX3026724 GSM2714289 SRP113298 SRS2375815 GSM2714289 RNA-Seq TRANSCRIPTOMIC cDNA Worker and drone embryos were collected from the distinct type of frame laid by the same queen. We sampled two groups of worker and drone embryos laid by two different honeybee queens. Then a total of 36 embryos were collected for single-embryo RNA-seq. High-quality, full-length cDNA is generated directly from 1–1,000 cells by The SMART-Seq v4 Ultra Low Input RNA Kit for Sequencing (Cat. Nos. 634890). cDNA was amplified by LD PCR. RNA can be transcribed by T7 promoter from ds cDNA , RNA was treated with RQ1 DNase (promega) to remove DNA .The quantity of the purified DNA were determined by Qubit. For each sample, 200 ng RNA was used for RNA-seq library preparation. RNAs were iron fragmented at 95℃ followed by end repair and 5' adaptor ligation. Then reverse transcription was performed with RT primer harboring 3' adaptor sequence and randomized hexamer. The cDNAs were purified and amplified and PCR products corresponding to 200-500 bps were purified, quantified and stored at -80 ℃ until used for sequencing. NextSeq 500 gds 302714289 GSM2714289 GEO Accession GSM2714289 SRX3026725 GSM2714290 SRP113298 SRS2375816 GSM2714290 RNA-Seq TRANSCRIPTOMIC cDNA Worker and drone embryos were collected from the distinct type of frame laid by the same queen. We sampled two groups of worker and drone embryos laid by two different honeybee queens. Then a total of 36 embryos were collected for single-embryo RNA-seq. High-quality, full-length cDNA is generated directly from 1–1,000 cells by The SMART-Seq v4 Ultra Low Input RNA Kit for Sequencing (Cat. Nos. 634890). cDNA was amplified by LD PCR. RNA can be transcribed by T7 promoter from ds cDNA , RNA was treated with RQ1 DNase (promega) to remove DNA .The quantity of the purified DNA were determined by Qubit. For each sample, 200 ng RNA was used for RNA-seq library preparation. RNAs were iron fragmented at 95℃ followed by end repair and 5' adaptor ligation. Then reverse transcription was performed with RT primer harboring 3' adaptor sequence and randomized hexamer. The cDNAs were purified and amplified and PCR products corresponding to 200-500 bps were purified, quantified and stored at -80 ℃ until used for sequencing. NextSeq 500 gds 302714290 GSM2714290 GEO Accession GSM2714290 SRX3026726 GSM2714291 SRP113298 SRS2375817 GSM2714291 RNA-Seq TRANSCRIPTOMIC cDNA Worker and drone embryos were collected from the distinct type of frame laid by the same queen. We sampled two groups of worker and drone embryos laid by two different honeybee queens. Then a total of 36 embryos were collected for single-embryo RNA-seq. High-quality, full-length cDNA is generated directly from 1–1,000 cells by The SMART-Seq v4 Ultra Low Input RNA Kit for Sequencing (Cat. Nos. 634890). cDNA was amplified by LD PCR. RNA can be transcribed by T7 promoter from ds cDNA , RNA was treated with RQ1 DNase (promega) to remove DNA .The quantity of the purified DNA were determined by Qubit. For each sample, 200 ng RNA was used for RNA-seq library preparation. RNAs were iron fragmented at 95℃ followed by end repair and 5' adaptor ligation. Then reverse transcription was performed with RT primer harboring 3' adaptor sequence and randomized hexamer. The cDNAs were purified and amplified and PCR products corresponding to 200-500 bps were purified, quantified and stored at -80 ℃ until used for sequencing. NextSeq 500 gds 302714291 GSM2714291 GEO Accession GSM2714291 SRX3026727 GSM2714292 SRP113298 SRS2375818 GSM2714292 RNA-Seq TRANSCRIPTOMIC cDNA Worker and drone embryos were collected from the distinct type of frame laid by the same queen. We sampled two groups of worker and drone embryos laid by two different honeybee queens. Then a total of 36 embryos were collected for single-embryo RNA-seq. High-quality, full-length cDNA is generated directly from 1–1,000 cells by The SMART-Seq v4 Ultra Low Input RNA Kit for Sequencing (Cat. Nos. 634890). cDNA was amplified by LD PCR. RNA can be transcribed by T7 promoter from ds cDNA , RNA was treated with RQ1 DNase (promega) to remove DNA .The quantity of the purified DNA were determined by Qubit. For each sample, 200 ng RNA was used for RNA-seq library preparation. RNAs were iron fragmented at 95℃ followed by end repair and 5' adaptor ligation. Then reverse transcription was performed with RT primer harboring 3' adaptor sequence and randomized hexamer. The cDNAs were purified and amplified and PCR products corresponding to 200-500 bps were purified, quantified and stored at -80 ℃ until used for sequencing. NextSeq 500 gds 302714292 GSM2714292 GEO Accession GSM2714292 SRX3026728 GSM2714293 SRP113298 SRS2375819 GSM2714293 RNA-Seq TRANSCRIPTOMIC cDNA Worker and drone embryos were collected from the distinct type of frame laid by the same queen. We sampled two groups of worker and drone embryos laid by two different honeybee queens. Then a total of 36 embryos were collected for single-embryo RNA-seq. High-quality, full-length cDNA is generated directly from 1–1,000 cells by The SMART-Seq v4 Ultra Low Input RNA Kit for Sequencing (Cat. Nos. 634890). cDNA was amplified by LD PCR. RNA can be transcribed by T7 promoter from ds cDNA , RNA was treated with RQ1 DNase (promega) to remove DNA .The quantity of the purified DNA were determined by Qubit. For each sample, 200 ng RNA was used for RNA-seq library preparation. RNAs were iron fragmented at 95℃ followed by end repair and 5' adaptor ligation. Then reverse transcription was performed with RT primer harboring 3' adaptor sequence and randomized hexamer. The cDNAs were purified and amplified and PCR products corresponding to 200-500 bps were purified, quantified and stored at -80 ℃ until used for sequencing. NextSeq 500 gds 302714293 GSM2714293 GEO Accession GSM2714293 SRX3026729 GSM2714294 SRP113298 SRS2375820 GSM2714294 RNA-Seq TRANSCRIPTOMIC cDNA Worker and drone embryos were collected from the distinct type of frame laid by the same queen. We sampled two groups of worker and drone embryos laid by two different honeybee queens. Then a total of 36 embryos were collected for single-embryo RNA-seq. High-quality, full-length cDNA is generated directly from 1–1,000 cells by The SMART-Seq v4 Ultra Low Input RNA Kit for Sequencing (Cat. Nos. 634890). cDNA was amplified by LD PCR. RNA can be transcribed by T7 promoter from ds cDNA , RNA was treated with RQ1 DNase (promega) to remove DNA .The quantity of the purified DNA were determined by Qubit. For each sample, 200 ng RNA was used for RNA-seq library preparation. RNAs were iron fragmented at 95℃ followed by end repair and 5' adaptor ligation. Then reverse transcription was performed with RT primer harboring 3' adaptor sequence and randomized hexamer. The cDNAs were purified and amplified and PCR products corresponding to 200-500 bps were purified, quantified and stored at -80 ℃ until used for sequencing. NextSeq 500 gds 302714294 GSM2714294 GEO Accession GSM2714294 SRX3026730 GSM2714295 SRP113298 SRS2375821 GSM2714295 RNA-Seq TRANSCRIPTOMIC cDNA Worker and drone embryos were collected from the distinct type of frame laid by the same queen. We sampled two groups of worker and drone embryos laid by two different honeybee queens. Then a total of 36 embryos were collected for single-embryo RNA-seq. High-quality, full-length cDNA is generated directly from 1–1,000 cells by The SMART-Seq v4 Ultra Low Input RNA Kit for Sequencing (Cat. Nos. 634890). cDNA was amplified by LD PCR. RNA can be transcribed by T7 promoter from ds cDNA , RNA was treated with RQ1 DNase (promega) to remove DNA .The quantity of the purified DNA were determined by Qubit. For each sample, 200 ng RNA was used for RNA-seq library preparation. RNAs were iron fragmented at 95℃ followed by end repair and 5' adaptor ligation. Then reverse transcription was performed with RT primer harboring 3' adaptor sequence and randomized hexamer. The cDNAs were purified and amplified and PCR products corresponding to 200-500 bps were purified, quantified and stored at -80 ℃ until used for sequencing. NextSeq 500 gds 302714295 GSM2714295 GEO Accession GSM2714295 SRX3026731 GSM2714296 SRP113298 SRS2375822 GSM2714296 RNA-Seq TRANSCRIPTOMIC cDNA Worker and drone embryos were collected from the distinct type of frame laid by the same queen. We sampled two groups of worker and drone embryos laid by two different honeybee queens. Then a total of 36 embryos were collected for single-embryo RNA-seq. High-quality, full-length cDNA is generated directly from 1–1,000 cells by The SMART-Seq v4 Ultra Low Input RNA Kit for Sequencing (Cat. Nos. 634890). cDNA was amplified by LD PCR. RNA can be transcribed by T7 promoter from ds cDNA , RNA was treated with RQ1 DNase (promega) to remove DNA .The quantity of the purified DNA were determined by Qubit. For each sample, 200 ng RNA was used for RNA-seq library preparation. RNAs were iron fragmented at 95℃ followed by end repair and 5' adaptor ligation. Then reverse transcription was performed with RT primer harboring 3' adaptor sequence and randomized hexamer. The cDNAs were purified and amplified and PCR products corresponding to 200-500 bps were purified, quantified and stored at -80 ℃ until used for sequencing. NextSeq 500 gds 302714296 GSM2714296 GEO Accession GSM2714296 SRX3026732 GSM2714297 SRP113298 SRS2375824 GSM2714297 RNA-Seq TRANSCRIPTOMIC cDNA Worker and drone embryos were collected from the distinct type of frame laid by the same queen. We sampled two groups of worker and drone embryos laid by two different honeybee queens. Then a total of 36 embryos were collected for single-embryo RNA-seq. High-quality, full-length cDNA is generated directly from 1–1,000 cells by The SMART-Seq v4 Ultra Low Input RNA Kit for Sequencing (Cat. Nos. 634890). cDNA was amplified by LD PCR. RNA can be transcribed by T7 promoter from ds cDNA , RNA was treated with RQ1 DNase (promega) to remove DNA .The quantity of the purified DNA were determined by Qubit. For each sample, 200 ng RNA was used for RNA-seq library preparation. RNAs were iron fragmented at 95℃ followed by end repair and 5' adaptor ligation. Then reverse transcription was performed with RT primer harboring 3' adaptor sequence and randomized hexamer. The cDNAs were purified and amplified and PCR products corresponding to 200-500 bps were purified, quantified and stored at -80 ℃ until used for sequencing. NextSeq 500 gds 302714297 GSM2714297 GEO Accession GSM2714297 SRX3026733 GSM2714298 SRP113298 SRS2375823 GSM2714298 RNA-Seq TRANSCRIPTOMIC cDNA Worker and drone embryos were collected from the distinct type of frame laid by the same queen. We sampled two groups of worker and drone embryos laid by two different honeybee queens. Then a total of 36 embryos were collected for single-embryo RNA-seq. High-quality, full-length cDNA is generated directly from 1–1,000 cells by The SMART-Seq v4 Ultra Low Input RNA Kit for Sequencing (Cat. Nos. 634890). cDNA was amplified by LD PCR. RNA can be transcribed by T7 promoter from ds cDNA , RNA was treated with RQ1 DNase (promega) to remove DNA .The quantity of the purified DNA were determined by Qubit. For each sample, 200 ng RNA was used for RNA-seq library preparation. RNAs were iron fragmented at 95℃ followed by end repair and 5' adaptor ligation. Then reverse transcription was performed with RT primer harboring 3' adaptor sequence and randomized hexamer. The cDNAs were purified and amplified and PCR products corresponding to 200-500 bps were purified, quantified and stored at -80 ℃ until used for sequencing. NextSeq 500 gds 302714298 GSM2714298 GEO Accession GSM2714298 SRX3026734 GSM2714299 SRP113298 SRS2375826 GSM2714299 RNA-Seq TRANSCRIPTOMIC cDNA Worker and drone embryos were collected from the distinct type of frame laid by the same queen. We sampled two groups of worker and drone embryos laid by two different honeybee queens. Then a total of 36 embryos were collected for single-embryo RNA-seq. High-quality, full-length cDNA is generated directly from 1–1,000 cells by The SMART-Seq v4 Ultra Low Input RNA Kit for Sequencing (Cat. Nos. 634890). cDNA was amplified by LD PCR. RNA can be transcribed by T7 promoter from ds cDNA , RNA was treated with RQ1 DNase (promega) to remove DNA .The quantity of the purified DNA were determined by Qubit. For each sample, 200 ng RNA was used for RNA-seq library preparation. RNAs were iron fragmented at 95℃ followed by end repair and 5' adaptor ligation. Then reverse transcription was performed with RT primer harboring 3' adaptor sequence and randomized hexamer. The cDNAs were purified and amplified and PCR products corresponding to 200-500 bps were purified, quantified and stored at -80 ℃ until used for sequencing. NextSeq 500 gds 302714299 GSM2714299 GEO Accession GSM2714299 SRX3026735 GSM2714300 SRP113298 SRS2375825 GSM2714300 RNA-Seq TRANSCRIPTOMIC cDNA Worker and drone embryos were collected from the distinct type of frame laid by the same queen. We sampled two groups of worker and drone embryos laid by two different honeybee queens. Then a total of 36 embryos were collected for single-embryo RNA-seq. High-quality, full-length cDNA is generated directly from 1–1,000 cells by The SMART-Seq v4 Ultra Low Input RNA Kit for Sequencing (Cat. Nos. 634890). cDNA was amplified by LD PCR. RNA can be transcribed by T7 promoter from ds cDNA , RNA was treated with RQ1 DNase (promega) to remove DNA .The quantity of the purified DNA were determined by Qubit. For each sample, 200 ng RNA was used for RNA-seq library preparation. RNAs were iron fragmented at 95℃ followed by end repair and 5' adaptor ligation. Then reverse transcription was performed with RT primer harboring 3' adaptor sequence and randomized hexamer. The cDNAs were purified and amplified and PCR products corresponding to 200-500 bps were purified, quantified and stored at -80 ℃ until used for sequencing. NextSeq 500 gds 302714300 GSM2714300 GEO Accession GSM2714300 SRX3026736 GSM2714301 SRP113298 SRS2375827 GSM2714301 RNA-Seq TRANSCRIPTOMIC cDNA Worker and drone embryos were collected from the distinct type of frame laid by the same queen. We sampled two groups of worker and drone embryos laid by two different honeybee queens. Then a total of 36 embryos were collected for single-embryo RNA-seq. High-quality, full-length cDNA is generated directly from 1–1,000 cells by The SMART-Seq v4 Ultra Low Input RNA Kit for Sequencing (Cat. Nos. 634890). cDNA was amplified by LD PCR. RNA can be transcribed by T7 promoter from ds cDNA , RNA was treated with RQ1 DNase (promega) to remove DNA .The quantity of the purified DNA were determined by Qubit. For each sample, 200 ng RNA was used for RNA-seq library preparation. RNAs were iron fragmented at 95℃ followed by end repair and 5' adaptor ligation. Then reverse transcription was performed with RT primer harboring 3' adaptor sequence and randomized hexamer. The cDNAs were purified and amplified and PCR products corresponding to 200-500 bps were purified, quantified and stored at -80 ℃ until used for sequencing. NextSeq 500 gds 302714301 GSM2714301 GEO Accession GSM2714301 SRX3026737 GSM2714302 SRP113298 SRS2375828 GSM2714302 RNA-Seq TRANSCRIPTOMIC cDNA Worker and drone embryos were collected from the distinct type of frame laid by the same queen. We sampled two groups of worker and drone embryos laid by two different honeybee queens. Then a total of 36 embryos were collected for single-embryo RNA-seq. High-quality, full-length cDNA is generated directly from 1–1,000 cells by The SMART-Seq v4 Ultra Low Input RNA Kit for Sequencing (Cat. Nos. 634890). cDNA was amplified by LD PCR. RNA can be transcribed by T7 promoter from ds cDNA , RNA was treated with RQ1 DNase (promega) to remove DNA .The quantity of the purified DNA were determined by Qubit. For each sample, 200 ng RNA was used for RNA-seq library preparation. RNAs were iron fragmented at 95℃ followed by end repair and 5' adaptor ligation. Then reverse transcription was performed with RT primer harboring 3' adaptor sequence and randomized hexamer. The cDNAs were purified and amplified and PCR products corresponding to 200-500 bps were purified, quantified and stored at -80 ℃ until used for sequencing. NextSeq 500 gds 302714302 GSM2714302 GEO Accession GSM2714302 SRX3026738 GSM2714303 SRP113298 SRS2375829 GSM2714303 RNA-Seq TRANSCRIPTOMIC cDNA Worker and drone embryos were collected from the distinct type of frame laid by the same queen. We sampled two groups of worker and drone embryos laid by two different honeybee queens. Then a total of 36 embryos were collected for single-embryo RNA-seq. High-quality, full-length cDNA is generated directly from 1–1,000 cells by The SMART-Seq v4 Ultra Low Input RNA Kit for Sequencing (Cat. Nos. 634890). cDNA was amplified by LD PCR. RNA can be transcribed by T7 promoter from ds cDNA , RNA was treated with RQ1 DNase (promega) to remove DNA .The quantity of the purified DNA were determined by Qubit. For each sample, 200 ng RNA was used for RNA-seq library preparation. RNAs were iron fragmented at 95℃ followed by end repair and 5' adaptor ligation. Then reverse transcription was performed with RT primer harboring 3' adaptor sequence and randomized hexamer. The cDNAs were purified and amplified and PCR products corresponding to 200-500 bps were purified, quantified and stored at -80 ℃ until used for sequencing. NextSeq 500 gds 302714303 GSM2714303 GEO Accession GSM2714303 SRX3026739 GSM2714304 SRP113298 SRS2375830 GSM2714304 RNA-Seq TRANSCRIPTOMIC cDNA Worker and drone embryos were collected from the distinct type of frame laid by the same queen. We sampled two groups of worker and drone embryos laid by two different honeybee queens. Then a total of 36 embryos were collected for single-embryo RNA-seq. High-quality, full-length cDNA is generated directly from 1–1,000 cells by The SMART-Seq v4 Ultra Low Input RNA Kit for Sequencing (Cat. Nos. 634890). cDNA was amplified by LD PCR. RNA can be transcribed by T7 promoter from ds cDNA , RNA was treated with RQ1 DNase (promega) to remove DNA .The quantity of the purified DNA were determined by Qubit. For each sample, 200 ng RNA was used for RNA-seq library preparation. RNAs were iron fragmented at 95℃ followed by end repair and 5' adaptor ligation. Then reverse transcription was performed with RT primer harboring 3' adaptor sequence and randomized hexamer. The cDNAs were purified and amplified and PCR products corresponding to 200-500 bps were purified, quantified and stored at -80 ℃ until used for sequencing. NextSeq 500 gds 302714304 GSM2714304 GEO Accession GSM2714304 SRX3026740 GSM2714305 SRP113298 SRS2375831 GSM2714305 RNA-Seq TRANSCRIPTOMIC cDNA Worker and drone embryos were collected from the distinct type of frame laid by the same queen. We sampled two groups of worker and drone embryos laid by two different honeybee queens. Then a total of 36 embryos were collected for single-embryo RNA-seq. High-quality, full-length cDNA is generated directly from 1–1,000 cells by The SMART-Seq v4 Ultra Low Input RNA Kit for Sequencing (Cat. Nos. 634890). cDNA was amplified by LD PCR. RNA can be transcribed by T7 promoter from ds cDNA , RNA was treated with RQ1 DNase (promega) to remove DNA .The quantity of the purified DNA were determined by Qubit. For each sample, 200 ng RNA was used for RNA-seq library preparation. RNAs were iron fragmented at 95℃ followed by end repair and 5' adaptor ligation. Then reverse transcription was performed with RT primer harboring 3' adaptor sequence and randomized hexamer. The cDNAs were purified and amplified and PCR products corresponding to 200-500 bps were purified, quantified and stored at -80 ℃ until used for sequencing. NextSeq 500 gds 302714305 GSM2714305 GEO Accession GSM2714305 SRX3026741 GSM2714306 SRP113298 SRS2375832 GSM2714306 RNA-Seq TRANSCRIPTOMIC cDNA Worker and drone embryos were collected from the distinct type of frame laid by the same queen. We sampled two groups of worker and drone embryos laid by two different honeybee queens. Then a total of 36 embryos were collected for single-embryo RNA-seq. High-quality, full-length cDNA is generated directly from 1–1,000 cells by The SMART-Seq v4 Ultra Low Input RNA Kit for Sequencing (Cat. Nos. 634890). cDNA was amplified by LD PCR. RNA can be transcribed by T7 promoter from ds cDNA , RNA was treated with RQ1 DNase (promega) to remove DNA .The quantity of the purified DNA were determined by Qubit. For each sample, 200 ng RNA was used for RNA-seq library preparation. RNAs were iron fragmented at 95℃ followed by end repair and 5' adaptor ligation. Then reverse transcription was performed with RT primer harboring 3' adaptor sequence and randomized hexamer. The cDNAs were purified and amplified and PCR products corresponding to 200-500 bps were purified, quantified and stored at -80 ℃ until used for sequencing. NextSeq 500 gds 302714306 GSM2714306 GEO Accession GSM2714306 SRX3026742 GSM2714307 SRP113298 SRS2375833 GSM2714307 RNA-Seq TRANSCRIPTOMIC cDNA Worker and drone embryos were collected from the distinct type of frame laid by the same queen. We sampled two groups of worker and drone embryos laid by two different honeybee queens. Then a total of 36 embryos were collected for single-embryo RNA-seq. High-quality, full-length cDNA is generated directly from 1–1,000 cells by The SMART-Seq v4 Ultra Low Input RNA Kit for Sequencing (Cat. Nos. 634890). cDNA was amplified by LD PCR. RNA can be transcribed by T7 promoter from ds cDNA , RNA was treated with RQ1 DNase (promega) to remove DNA .The quantity of the purified DNA were determined by Qubit. For each sample, 200 ng RNA was used for RNA-seq library preparation. RNAs were iron fragmented at 95℃ followed by end repair and 5' adaptor ligation. Then reverse transcription was performed with RT primer harboring 3' adaptor sequence and randomized hexamer. The cDNAs were purified and amplified and PCR products corresponding to 200-500 bps were purified, quantified and stored at -80 ℃ until used for sequencing. NextSeq 500 gds 302714307 GSM2714307 GEO Accession GSM2714307 SRX3026743 GSM2714308 SRP113298 SRS2375834 GSM2714308 RNA-Seq TRANSCRIPTOMIC cDNA Worker and drone embryos were collected from the distinct type of frame laid by the same queen. We sampled two groups of worker and drone embryos laid by two different honeybee queens. Then a total of 36 embryos were collected for single-embryo RNA-seq. High-quality, full-length cDNA is generated directly from 1–1,000 cells by The SMART-Seq v4 Ultra Low Input RNA Kit for Sequencing (Cat. Nos. 634890). cDNA was amplified by LD PCR. RNA can be transcribed by T7 promoter from ds cDNA , RNA was treated with RQ1 DNase (promega) to remove DNA .The quantity of the purified DNA were determined by Qubit. For each sample, 200 ng RNA was used for RNA-seq library preparation. RNAs were iron fragmented at 95℃ followed by end repair and 5' adaptor ligation. Then reverse transcription was performed with RT primer harboring 3' adaptor sequence and randomized hexamer. The cDNAs were purified and amplified and PCR products corresponding to 200-500 bps were purified, quantified and stored at -80 ℃ until used for sequencing. NextSeq 500 gds 302714308 GSM2714308 GEO Accession GSM2714308 SRX3026744 GSM2714309 SRP113298 SRS2375835 GSM2714309 RNA-Seq TRANSCRIPTOMIC cDNA Worker and drone embryos were collected from the distinct type of frame laid by the same queen. We sampled two groups of worker and drone embryos laid by two different honeybee queens. Then a total of 36 embryos were collected for single-embryo RNA-seq. High-quality, full-length cDNA is generated directly from 1–1,000 cells by The SMART-Seq v4 Ultra Low Input RNA Kit for Sequencing (Cat. Nos. 634890). cDNA was amplified by LD PCR. RNA can be transcribed by T7 promoter from ds cDNA , RNA was treated with RQ1 DNase (promega) to remove DNA .The quantity of the purified DNA were determined by Qubit. For each sample, 200 ng RNA was used for RNA-seq library preparation. RNAs were iron fragmented at 95℃ followed by end repair and 5' adaptor ligation. Then reverse transcription was performed with RT primer harboring 3' adaptor sequence and randomized hexamer. The cDNAs were purified and amplified and PCR products corresponding to 200-500 bps were purified, quantified and stored at -80 ℃ until used for sequencing. NextSeq 500 gds 302714309 GSM2714309 GEO Accession GSM2714309 SRX3026745 GSM2714310 SRP113298 SRS2375837 GSM2714310 RNA-Seq TRANSCRIPTOMIC cDNA Worker and drone embryos were collected from the distinct type of frame laid by the same queen. We sampled two groups of worker and drone embryos laid by two different honeybee queens. Then a total of 36 embryos were collected for single-embryo RNA-seq. High-quality, full-length cDNA is generated directly from 1–1,000 cells by The SMART-Seq v4 Ultra Low Input RNA Kit for Sequencing (Cat. Nos. 634890). cDNA was amplified by LD PCR. RNA can be transcribed by T7 promoter from ds cDNA , RNA was treated with RQ1 DNase (promega) to remove DNA .The quantity of the purified DNA were determined by Qubit. For each sample, 200 ng RNA was used for RNA-seq library preparation. RNAs were iron fragmented at 95℃ followed by end repair and 5' adaptor ligation. Then reverse transcription was performed with RT primer harboring 3' adaptor sequence and randomized hexamer. The cDNAs were purified and amplified and PCR products corresponding to 200-500 bps were purified, quantified and stored at -80 ℃ until used for sequencing. NextSeq 500 gds 302714310 GSM2714310 GEO Accession GSM2714310 SRX3026746 GSM2714311 SRP113298 SRS2375836 GSM2714311 RNA-Seq TRANSCRIPTOMIC cDNA Worker and drone embryos were collected from the distinct type of frame laid by the same queen. We sampled two groups of worker and drone embryos laid by two different honeybee queens. Then a total of 36 embryos were collected for single-embryo RNA-seq. High-quality, full-length cDNA is generated directly from 1–1,000 cells by The SMART-Seq v4 Ultra Low Input RNA Kit for Sequencing (Cat. Nos. 634890). cDNA was amplified by LD PCR. RNA can be transcribed by T7 promoter from ds cDNA , RNA was treated with RQ1 DNase (promega) to remove DNA .The quantity of the purified DNA were determined by Qubit. For each sample, 200 ng RNA was used for RNA-seq library preparation. RNAs were iron fragmented at 95℃ followed by end repair and 5' adaptor ligation. Then reverse transcription was performed with RT primer harboring 3' adaptor sequence and randomized hexamer. The cDNAs were purified and amplified and PCR products corresponding to 200-500 bps were purified, quantified and stored at -80 ℃ until used for sequencing. NextSeq 500 gds 302714311 GSM2714311 GEO Accession GSM2714311 SRX3026747 GSM2714312 SRP113298 SRS2375838 GSM2714312 RNA-Seq TRANSCRIPTOMIC cDNA Worker and drone embryos were collected from the distinct type of frame laid by the same queen. We sampled two groups of worker and drone embryos laid by two different honeybee queens. Then a total of 36 embryos were collected for single-embryo RNA-seq. High-quality, full-length cDNA is generated directly from 1–1,000 cells by The SMART-Seq v4 Ultra Low Input RNA Kit for Sequencing (Cat. Nos. 634890). cDNA was amplified by LD PCR. RNA can be transcribed by T7 promoter from ds cDNA , RNA was treated with RQ1 DNase (promega) to remove DNA .The quantity of the purified DNA were determined by Qubit. For each sample, 200 ng RNA was used for RNA-seq library preparation. RNAs were iron fragmented at 95℃ followed by end repair and 5' adaptor ligation. Then reverse transcription was performed with RT primer harboring 3' adaptor sequence and randomized hexamer. The cDNAs were purified and amplified and PCR products corresponding to 200-500 bps were purified, quantified and stored at -80 ℃ until used for sequencing. NextSeq 500 gds 302714312 GSM2714312 GEO Accession GSM2714312 SRX3026748 GSM2714313 SRP113298 SRS2375839 GSM2714313 RNA-Seq TRANSCRIPTOMIC cDNA Worker and drone embryos were collected from the distinct type of frame laid by the same queen. We sampled two groups of worker and drone embryos laid by two different honeybee queens. Then a total of 36 embryos were collected for single-embryo RNA-seq. High-quality, full-length cDNA is generated directly from 1–1,000 cells by The SMART-Seq v4 Ultra Low Input RNA Kit for Sequencing (Cat. Nos. 634890). cDNA was amplified by LD PCR. RNA can be transcribed by T7 promoter from ds cDNA , RNA was treated with RQ1 DNase (promega) to remove DNA .The quantity of the purified DNA were determined by Qubit. For each sample, 200 ng RNA was used for RNA-seq library preparation. RNAs were iron fragmented at 95℃ followed by end repair and 5' adaptor ligation. Then reverse transcription was performed with RT primer harboring 3' adaptor sequence and randomized hexamer. The cDNAs were purified and amplified and PCR products corresponding to 200-500 bps were purified, quantified and stored at -80 ℃ until used for sequencing. NextSeq 500 gds 302714313 GSM2714313 GEO Accession GSM2714313 SRX3026749 GSM2714314 SRP113298 SRS2375840 GSM2714314 RNA-Seq TRANSCRIPTOMIC cDNA Worker and drone embryos were collected from the distinct type of frame laid by the same queen. We sampled two groups of worker and drone embryos laid by two different honeybee queens. Then a total of 36 embryos were collected for single-embryo RNA-seq. High-quality, full-length cDNA is generated directly from 1–1,000 cells by The SMART-Seq v4 Ultra Low Input RNA Kit for Sequencing (Cat. Nos. 634890). cDNA was amplified by LD PCR. RNA can be transcribed by T7 promoter from ds cDNA , RNA was treated with RQ1 DNase (promega) to remove DNA .The quantity of the purified DNA were determined by Qubit. For each sample, 200 ng RNA was used for RNA-seq library preparation. RNAs were iron fragmented at 95℃ followed by end repair and 5' adaptor ligation. Then reverse transcription was performed with RT primer harboring 3' adaptor sequence and randomized hexamer. The cDNAs were purified and amplified and PCR products corresponding to 200-500 bps were purified, quantified and stored at -80 ℃ until used for sequencing. NextSeq 500 gds 302714314 GSM2714314 GEO Accession GSM2714314 SRX3026750 GSM2714315 SRP113298 SRS2375841 GSM2714315 RNA-Seq TRANSCRIPTOMIC cDNA Worker and drone embryos were collected from the distinct type of frame laid by the same queen. We sampled two groups of worker and drone embryos laid by two different honeybee queens. Then a total of 36 embryos were collected for single-embryo RNA-seq. High-quality, full-length cDNA is generated directly from 1–1,000 cells by The SMART-Seq v4 Ultra Low Input RNA Kit for Sequencing (Cat. Nos. 634890). cDNA was amplified by LD PCR. RNA can be transcribed by T7 promoter from ds cDNA , RNA was treated with RQ1 DNase (promega) to remove DNA .The quantity of the purified DNA were determined by Qubit. For each sample, 200 ng RNA was used for RNA-seq library preparation. RNAs were iron fragmented at 95℃ followed by end repair and 5' adaptor ligation. Then reverse transcription was performed with RT primer harboring 3' adaptor sequence and randomized hexamer. The cDNAs were purified and amplified and PCR products corresponding to 200-500 bps were purified, quantified and stored at -80 ℃ until used for sequencing. NextSeq 500 gds 302714315 GSM2714315 GEO Accession GSM2714315 SRX3026751 GSM2714316 SRP113298 SRS2375842 GSM2714316 RNA-Seq TRANSCRIPTOMIC cDNA Worker and drone embryos were collected from the distinct type of frame laid by the same queen. We sampled two groups of worker and drone embryos laid by two different honeybee queens. Then a total of 36 embryos were collected for single-embryo RNA-seq. High-quality, full-length cDNA is generated directly from 1–1,000 cells by The SMART-Seq v4 Ultra Low Input RNA Kit for Sequencing (Cat. Nos. 634890). cDNA was amplified by LD PCR. RNA can be transcribed by T7 promoter from ds cDNA , RNA was treated with RQ1 DNase (promega) to remove DNA .The quantity of the purified DNA were determined by Qubit. For each sample, 200 ng RNA was used for RNA-seq library preparation. RNAs were iron fragmented at 95℃ followed by end repair and 5' adaptor ligation. Then reverse transcription was performed with RT primer harboring 3' adaptor sequence and randomized hexamer. The cDNAs were purified and amplified and PCR products corresponding to 200-500 bps were purified, quantified and stored at -80 ℃ until used for sequencing. NextSeq 500 gds 302714316 GSM2714316 GEO Accession GSM2714316 SRX3026752 GSM2714317 SRP113298 SRS2375843 GSM2714317 RNA-Seq TRANSCRIPTOMIC cDNA Worker and drone embryos were collected from the distinct type of frame laid by the same queen. We sampled two groups of worker and drone embryos laid by two different honeybee queens. Then a total of 36 embryos were collected for single-embryo RNA-seq. High-quality, full-length cDNA is generated directly from 1–1,000 cells by The SMART-Seq v4 Ultra Low Input RNA Kit for Sequencing (Cat. Nos. 634890). cDNA was amplified by LD PCR. RNA can be transcribed by T7 promoter from ds cDNA , RNA was treated with RQ1 DNase (promega) to remove DNA .The quantity of the purified DNA were determined by Qubit. For each sample, 200 ng RNA was used for RNA-seq library preparation. RNAs were iron fragmented at 95℃ followed by end repair and 5' adaptor ligation. Then reverse transcription was performed with RT primer harboring 3' adaptor sequence and randomized hexamer. The cDNAs were purified and amplified and PCR products corresponding to 200-500 bps were purified, quantified and stored at -80 ℃ until used for sequencing. NextSeq 500 gds 302714317 GSM2714317 GEO Accession GSM2714317