Abstract: 【Aim】 This study aims to understand the transcriptome features of Myzus persicae of different wing types (winged and apterous), to enrich its transcriptome data, and to reveal the molecular mechanisms of wing dimorphism of this insect. 【Methods】 The differences in transcriptome expression between nymphs, apterous adults, and winged adults of M. persicae were detected by high-throughput sequencing technology Illumina HiSeq^TM2000. The expression patterns of wing type related genes were analyzed by real-time quantitative PCR (qPCR). 【Results】 The percentage of nucleotides with the quality score over 20 was not less than 94% (GenBank accession numbers: SRR6112438, SRR6112436 and SRR6112437). All reads were assembled into 128 065 unigenes with an average length of 607.08 bp. The resulting sequences were annotated to databases NR, NT, KOG, GO, KEGG, etc., and a total of 128 065 unigenes were annotated. Gene expression analysis showed that the genes of flight proteins, peptide hormones, receptor proteins, etc., which are related to wing formation during growth, were differentially expressed between winged adults and apterous adults or nymphs. qPCR analysis showed that there were significant differences in the expression levels of wing type differentiation related genes including Bursicon (bur- α and bur-β), flightin (Fln), Wntless (Wnt), Distal-less (Dll), Decapentaplegic (Dpp), juvenile hormone epoxide hydrolase gene (JHEH) and ecdysone receptor gene (EcR) between winged adults and apterous adults. 【Conclusion】 The expression levels of bur-α, bur-β , Fln, Wnt, Dll, Dpp, JHEH and EcR change during wing formation in M. persicae. Our results show the general expression characteristics of wing formation genes and provide a basis for further research of wing formation in M. persicae.

Key words: Myzus persicae; transcriptome, wing development, differentially expressed gene, high-throughput sequencing