Comparative transcriptional analysis and identification of hub genes associated with wing differentiation of male in Aphis gossypii

Abstract Background Aphis gossypii Glover (Hemiptera: Aphididae), a worldwide polyphagous phloem-feeding agricultural pest, has three wing morphs (winged parthenogenetic female, gynopara, and male) in the life cycle. The exclusive males could fly from summer hosts to winter hosts, which are essential for gene exchanges of cotton aphid populations from different hosts or regions. However, the molecular mechanism of wing differentiation of male in A. gossypii remains unclear. Results Morphological observation of male A. gossypii showed that there is no distinct difference in the external morphologies of the 1st and 2nd instar nymphs. The obvious differentiation of wing buds started in the 3rd instar nymph and was visible via naked eyes in the 4th instar nymphal stage, then adult male emerged with full wings. According to morphological dynamic changes, the development of wings in males were divided into four stages: preliminary stage (the 1st instar to 2nd instar), prophase (the 3rd instar), metaphase (the 4th instar), anaphase (the 5th instar). Results of feeding behavior monitoring via EPG (electrical penetration graph) technology indicated that although the male cotton aphids had strong desire to feed (longer duration of C 55.24%, F 5.05% and Pd waves 2.56%), its feeding efficiency to summer host cotton was low (shorter E1 3.56% and E2 waves 2.63%). Dynamic transcriptome analysis of male aphid at 5 different developmental periods showed that in the 3rd instar nymph, the number of up-regulated DEGs was significant increased, and time-course gene transcriptional pattern analyses results also showed that numerous genes categorized in clusters 3, 5, and 8 had the highest expressed levels, which were consistent with morphological changes of wing buds. These results indicate that the 3rd instar nymph is the critical stage of wing bud differentiation in males. Furthermore, through pathway enrichment analysis of DEGs and WGCNA, it revealed that the neuroactive ligand-receptor interaction, Ras signaling pathway, dopaminergic synapse, circadian entrainment and the corresponding hub genes of PLK1, BUB1, SMC2, TUBG, ASPM, the kinesin family members (KIF23, KIF20, KIF18-19) and the novel subfamily of serine/threonine (Aurora kinase A and Aurora kinase B) probably played an important role in the critical stage of wing bud differentiation. Conclusion This study explored morphological changes and genes transcriptional dynamics males in cotton aphid, revealed the phenomenon of low feeding efficiency of winged males on summer host cotton, and identified key signaling pathways and potential hub genes potentially involved in wing bud differentiation of male in A. gossypii.