Differential expression of genes in the alate and apterous morphs of the brown citrus aphid, Toxoptera citricida
Identifieur interne : 000476 ( Pmc/Corpus ); précédent : 000475; suivant : 000477Differential expression of genes in the alate and apterous morphs of the brown citrus aphid, Toxoptera citricida
Auteurs : Feng Shang ; Bi-Yue Ding ; Ying Xiong ; Wei Dou ; Dong Wei ; Hong-Bo Jiang ; Dan-Dan Wei ; Jin-Jun WangSource :
- Scientific Reports [ 2045-2322 ] ; 2016.
Abstract
Winged and wingless morphs in insects represent a trade-off between dispersal ability and reproduction. We studied key genes associated with apterous and alate morphs in
Url:
DOI: 10.1038/srep32099
PubMed: 27577531
PubMed Central: 5006003
Links to Exploration step
PMC:5006003Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Differential expression of genes in the alate and apterous morphs of the brown citrus aphid, <italic>Toxoptera citricida</italic></title><author><name sortKey="Shang, Feng" sort="Shang, Feng" uniqKey="Shang F" first="Feng" last="Shang">Feng Shang</name><affiliation><nlm:aff id="a1"><institution>Key Laboratory of Entomology and Pest Control Engineering, College of Plant Protection, Southwest University</institution>, Chongqing 400716,<country>China</country></nlm:aff></affiliation></author><author><name sortKey="Ding, Bi Yue" sort="Ding, Bi Yue" uniqKey="Ding B" first="Bi-Yue" last="Ding">Bi-Yue Ding</name><affiliation><nlm:aff id="a1"><institution>Key Laboratory of Entomology and Pest Control Engineering, College of Plant Protection, Southwest University</institution>, Chongqing 400716,<country>China</country></nlm:aff></affiliation></author><author><name sortKey="Xiong, Ying" sort="Xiong, Ying" uniqKey="Xiong Y" first="Ying" last="Xiong">Ying Xiong</name><affiliation><nlm:aff id="a1"><institution>Key Laboratory of Entomology and Pest Control Engineering, College of Plant Protection, Southwest University</institution>, Chongqing 400716,<country>China</country></nlm:aff></affiliation></author><author><name sortKey="Dou, Wei" sort="Dou, Wei" uniqKey="Dou W" first="Wei" last="Dou">Wei Dou</name><affiliation><nlm:aff id="a1"><institution>Key Laboratory of Entomology and Pest Control Engineering, College of Plant Protection, Southwest University</institution>, Chongqing 400716,<country>China</country></nlm:aff></affiliation></author><author><name sortKey="Wei, Dong" sort="Wei, Dong" uniqKey="Wei D" first="Dong" last="Wei">Dong Wei</name><affiliation><nlm:aff id="a1"><institution>Key Laboratory of Entomology and Pest Control Engineering, College of Plant Protection, Southwest University</institution>, Chongqing 400716,<country>China</country></nlm:aff></affiliation></author><author><name sortKey="Jiang, Hong Bo" sort="Jiang, Hong Bo" uniqKey="Jiang H" first="Hong-Bo" last="Jiang">Hong-Bo Jiang</name><affiliation><nlm:aff id="a1"><institution>Key Laboratory of Entomology and Pest Control Engineering, College of Plant Protection, Southwest University</institution>, Chongqing 400716,<country>China</country></nlm:aff></affiliation></author><author><name sortKey="Wei, Dan Dan" sort="Wei, Dan Dan" uniqKey="Wei D" first="Dan-Dan" last="Wei">Dan-Dan Wei</name><affiliation><nlm:aff id="a1"><institution>Key Laboratory of Entomology and Pest Control Engineering, College of Plant Protection, Southwest University</institution>, Chongqing 400716,<country>China</country></nlm:aff></affiliation></author><author><name sortKey="Wang, Jin Jun" sort="Wang, Jin Jun" uniqKey="Wang J" first="Jin-Jun" last="Wang">Jin-Jun Wang</name><affiliation><nlm:aff id="a1"><institution>Key Laboratory of Entomology and Pest Control Engineering, College of Plant Protection, Southwest University</institution>, Chongqing 400716,<country>China</country></nlm:aff></affiliation></author></titleStmt><publicationStmt><idno type="wicri:source">PMC</idno><idno type="pmid">27577531</idno><idno type="pmc">5006003</idno><idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC5006003</idno><idno type="RBID">PMC:5006003</idno><idno type="doi">10.1038/srep32099</idno><date when="2016">2016</date><idno type="wicri:Area/Pmc/Corpus">000476</idno></publicationStmt><sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Differential expression of genes in the alate and apterous morphs of the brown citrus aphid, <italic>Toxoptera citricida</italic></title><author><name sortKey="Shang, Feng" sort="Shang, Feng" uniqKey="Shang F" first="Feng" last="Shang">Feng Shang</name><affiliation><nlm:aff id="a1"><institution>Key Laboratory of Entomology and Pest Control Engineering, College of Plant Protection, Southwest University</institution>, Chongqing 400716,<country>China</country></nlm:aff></affiliation></author><author><name sortKey="Ding, Bi Yue" sort="Ding, Bi Yue" uniqKey="Ding B" first="Bi-Yue" last="Ding">Bi-Yue Ding</name><affiliation><nlm:aff id="a1"><institution>Key Laboratory of Entomology and Pest Control Engineering, College of Plant Protection, Southwest University</institution>, Chongqing 400716,<country>China</country></nlm:aff></affiliation></author><author><name sortKey="Xiong, Ying" sort="Xiong, Ying" uniqKey="Xiong Y" first="Ying" last="Xiong">Ying Xiong</name><affiliation><nlm:aff id="a1"><institution>Key Laboratory of Entomology and Pest Control Engineering, College of Plant Protection, Southwest University</institution>, Chongqing 400716,<country>China</country></nlm:aff></affiliation></author><author><name sortKey="Dou, Wei" sort="Dou, Wei" uniqKey="Dou W" first="Wei" last="Dou">Wei Dou</name><affiliation><nlm:aff id="a1"><institution>Key Laboratory of Entomology and Pest Control Engineering, College of Plant Protection, Southwest University</institution>, Chongqing 400716,<country>China</country></nlm:aff></affiliation></author><author><name sortKey="Wei, Dong" sort="Wei, Dong" uniqKey="Wei D" first="Dong" last="Wei">Dong Wei</name><affiliation><nlm:aff id="a1"><institution>Key Laboratory of Entomology and Pest Control Engineering, College of Plant Protection, Southwest University</institution>, Chongqing 400716,<country>China</country></nlm:aff></affiliation></author><author><name sortKey="Jiang, Hong Bo" sort="Jiang, Hong Bo" uniqKey="Jiang H" first="Hong-Bo" last="Jiang">Hong-Bo Jiang</name><affiliation><nlm:aff id="a1"><institution>Key Laboratory of Entomology and Pest Control Engineering, College of Plant Protection, Southwest University</institution>, Chongqing 400716,<country>China</country></nlm:aff></affiliation></author><author><name sortKey="Wei, Dan Dan" sort="Wei, Dan Dan" uniqKey="Wei D" first="Dan-Dan" last="Wei">Dan-Dan Wei</name><affiliation><nlm:aff id="a1"><institution>Key Laboratory of Entomology and Pest Control Engineering, College of Plant Protection, Southwest University</institution>, Chongqing 400716,<country>China</country></nlm:aff></affiliation></author><author><name sortKey="Wang, Jin Jun" sort="Wang, Jin Jun" uniqKey="Wang J" first="Jin-Jun" last="Wang">Jin-Jun Wang</name><affiliation><nlm:aff id="a1"><institution>Key Laboratory of Entomology and Pest Control Engineering, College of Plant Protection, Southwest University</institution>, Chongqing 400716,<country>China</country></nlm:aff></affiliation></author></analytic><series><title level="j">Scientific Reports</title><idno type="eISSN">2045-2322</idno><imprint><date when="2016">2016</date></imprint></series></biblStruct></sourceDesc></fileDesc><profileDesc><textClass></textClass></profileDesc></teiHeader><front><div type="abstract" xml:lang="en"><p>Winged and wingless morphs in insects represent a trade-off between dispersal ability and reproduction. We studied key genes associated with apterous and alate morphs in <italic>Toxoptera citricida</italic> (Kirkaldy) using RNAseq, digital gene expression (DGE) profiling, and RNA interference. The <italic>de novo</italic> assembly of the transcriptome was obtained through Illumina short-read sequencing technology. A total of 44,199 unigenes were generated and 27,640 were annotated. The transcriptomic differences between alate and apterous adults indicated that 279 unigenes were highly expressed in alate adults, whereas 5,470 were expressed at low levels. Expression patterns of the top 10 highly expressed genes in alate adults agreed with wing bud development trends. Silencing of the lipid synthesis and degradation gene (3-ketoacyl-CoA thiolase, mitochondrial-like) and glycogen genes (Phosphoenolpyruvate carboxykinase [GTP]-like and Glycogen phosphorylase-like isoform 2) resulted in underdeveloped wings. This suggests that both lipid and glycogen metabolism provide energy for aphid wing development. The large number of sequences and expression data produced from the transcriptome and DGE sequencing, respectively, increases our understanding of wing development mechanisms.</p></div></front><back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Denno, R F" uniqKey="Denno R">R. F. Denno</name></author><author><name sortKey="Olmstead, K L" uniqKey="Olmstead K">K. L. Olmstead</name></author><author><name sortKey="Mccloud, E" uniqKey="Mccloud E">E. Mccloud</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zhao, Z" uniqKey="Zhao Z">Z. Zhao</name></author><author><name sortKey="Zera, A J" uniqKey="Zera A">A. J. Zera</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Tanaka, S" uniqKey="Tanaka S">S. Tanaka</name></author><author><name sortKey="Nishide, Y" uniqKey="Nishide Y">Y. Nishide</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Xu, X" uniqKey="Xu X">X. Xu</name></author><author><name sortKey="Liu, X" uniqKey="Liu X">X. Liu</name></author><author><name sortKey="Zhang, Q" uniqKey="Zhang Q">Q. Zhang</name></author><author><name sortKey="Wu, J" uniqKey="Wu J">J. Wu</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Yang, X" uniqKey="Yang X">X. Yang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Harrison, R G" uniqKey="Harrison R">R. G. Harrison</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Lees, A" uniqKey="Lees A">A. Lees</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Sutherland, O" uniqKey="Sutherland O">O. Sutherland</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Michaud, J" uniqKey="Michaud J">J. Michaud</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="De Barro, P" uniqKey="De Barro P">P. De Barro</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zhang, F" uniqKey="Zhang F">F. Zhang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Andersen, N M" uniqKey="Andersen N">N. M. Andersen</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Tanaka, S" uniqKey="Tanaka S">S. Tanaka</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Langellotto, G A" uniqKey="Langellotto G">G. A. Langellotto</name></author><author><name sortKey="Denno, R F" uniqKey="Denno R">R. F. Denno</name></author><author><name sortKey="Ott, J R" uniqKey="Ott J">J. R. Ott</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Braendle, C" uniqKey="Braendle C">C. Braendle</name></author><author><name sortKey="Davis, G K" uniqKey="Davis G">G. K. Davis</name></author><author><name sortKey="Brisson, J A" uniqKey="Brisson J">J. A. Brisson</name></author><author><name sortKey="Stern, D L" uniqKey="Stern D">D. L. Stern</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Braendle, C" uniqKey="Braendle C">C. Braendle</name></author><author><name sortKey="Friebe, I" uniqKey="Friebe I">I. Friebe</name></author><author><name sortKey="Caillaud, M C" uniqKey="Caillaud M">M. C. Caillaud</name></author><author><name sortKey="Stern, D L" uniqKey="Stern D">D. L. Stern</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Mackay, P A" uniqKey="Mackay P">P. A. Mackay</name></author><author><name sortKey="Wellington, W" uniqKey="Wellington W">W. Wellington</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Tsumuki, H" uniqKey="Tsumuki H">H. Tsumuki</name></author><author><name sortKey="Nagatsuka, H" uniqKey="Nagatsuka H">H. Nagatsuka</name></author><author><name sortKey="Kawada, K" uniqKey="Kawada K">K. Kawada</name></author><author><name sortKey="Kanehisa, K" uniqKey="Kanehisa K">K. Kanehisa</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hazell, S" uniqKey="Hazell S">S. Hazell</name></author><author><name sortKey="Gwynn, D" uniqKey="Gwynn D">D. Gwynn</name></author><author><name sortKey="Ceccarelli, S" uniqKey="Ceccarelli S">S. Ceccarelli</name></author><author><name sortKey="Fellowes, M" uniqKey="Fellowes M">M. Fellowes</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Fasulo, T R" uniqKey="Fasulo T">T. R. Fasulo</name></author><author><name sortKey="Halbert, S E" uniqKey="Halbert S">S. E. Halbert</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Halbert, S E" uniqKey="Halbert S">S. E. Halbert</name></author><author><name sortKey="Brown, L G" uniqKey="Brown L">L. G. Brown</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Consortium, I A G" uniqKey="Consortium I">I. A. G. Consortium</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Nicholson, S J" uniqKey="Nicholson S">S. J. Nicholson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hunter, W" uniqKey="Hunter W">W. Hunter</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Tang, L D" uniqKey="Tang L">L. D. Tang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zhu, J Y" uniqKey="Zhu J">J. Y. Zhu</name></author><author><name sortKey="Yang, P" uniqKey="Yang P">P. Yang</name></author><author><name sortKey="Zhang, Z" uniqKey="Zhang Z">Z. Zhang</name></author><author><name sortKey="Wu, G X" uniqKey="Wu G">G. X. Wu</name></author><author><name sortKey="Yang, B" uniqKey="Yang B">B. Yang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Li, Z Q" uniqKey="Li Z">Z.-Q. Li</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Kola, V S R" uniqKey="Kola V">V. S. R. Kola</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Yang, W J" uniqKey="Yang W">W. J. Yang</name></author><author><name sortKey="Xu, K K" uniqKey="Xu K">K. K. Xu</name></author><author><name sortKey="Cong, L" uniqKey="Cong L">L. Cong</name></author><author><name sortKey="Wang, J J" uniqKey="Wang J">J. J. Wang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="El Shesheny, I" uniqKey="El Shesheny I">I. El-Shesheny</name></author><author><name sortKey="Hajeri, S" uniqKey="Hajeri S">S. Hajeri</name></author><author><name sortKey="El Hawary, I" uniqKey="El Hawary I">I. El-Hawary</name></author><author><name sortKey="Gowda, S" uniqKey="Gowda S">S. Gowda</name></author><author><name sortKey="Killiny, N" uniqKey="Killiny N">N. Killiny</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Casta Eda, L E" uniqKey="Casta Eda L">L. E. Castañeda</name></author><author><name sortKey="Figueroa, C C" uniqKey="Figueroa C">C. C. Figueroa</name></author><author><name sortKey="Bacigalupe, L D" uniqKey="Bacigalupe L">L. D. Bacigalupe</name></author><author><name sortKey="Nespolo, R F" uniqKey="Nespolo R">R. F. Nespolo</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Kim, H" uniqKey="Kim H">H. Kim</name></author><author><name sortKey="Lee, S" uniqKey="Lee S">S. Lee</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Von Dohlen, C D" uniqKey="Von Dohlen C">C. D. von Dohlen</name></author><author><name sortKey="Rowe, C A" uniqKey="Rowe C">C. A. Rowe</name></author><author><name sortKey="Heie, O E" uniqKey="Heie O">O. E. Heie</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Rivers, D B" uniqKey="Rivers D">D. B. Rivers</name></author><author><name sortKey="Hink, W F" uniqKey="Hink W">W. F. Hink</name></author><author><name sortKey="Denlinger, D L" uniqKey="Denlinger D">D. L. Denlinger</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Harvey, J A" uniqKey="Harvey J">J. A. Harvey</name></author><author><name sortKey="Bezemer, T M" uniqKey="Bezemer T">T. M. Bezemer</name></author><author><name sortKey="Elzinga, J A" uniqKey="Elzinga J">J. A. Elzinga</name></author><author><name sortKey="Strand, M R" uniqKey="Strand M">M. R. Strand</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Brogiolo, W" uniqKey="Brogiolo W">W. Brogiolo</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Teleman, A" uniqKey="Teleman A">A. Teleman</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Xu, H J" uniqKey="Xu H">H. J. Xu</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Yang, P" uniqKey="Yang P">P. Yang</name></author><author><name sortKey="Chen, X M" uniqKey="Chen X">X.-M. Chen</name></author><author><name sortKey="Liu, W W" uniqKey="Liu W">W.-W. Liu</name></author><author><name sortKey="Feng, Y" uniqKey="Feng Y">Y. Feng</name></author><author><name sortKey="Sun, T" uniqKey="Sun T">T. Sun</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Brisson, J A" uniqKey="Brisson J">J. A. Brisson</name></author><author><name sortKey="Ishikawa, A" uniqKey="Ishikawa A">A. Ishikawa</name></author><author><name sortKey="Miura, T" uniqKey="Miura T">T. Miura</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wang, X" uniqKey="Wang X">X. Wang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wu, H J" uniqKey="Wu H">H.-J. Wu</name></author><author><name sortKey="Zhao, L Q" uniqKey="Zhao L">L.-Q. Zhao</name></author><author><name sortKey="Zeng, Y" uniqKey="Zeng Y">Y. Zeng</name></author><author><name sortKey="Zhu, D H" uniqKey="Zhu D">D.-H. Zhu</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Brisson, J A" uniqKey="Brisson J">J. A. Brisson</name></author><author><name sortKey="Davis, G K" uniqKey="Davis G">G. K. Davis</name></author><author><name sortKey="Stern, D L" uniqKey="Stern D">D. L. Stern</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Grabherr, M G" uniqKey="Grabherr M">M. G. Grabherr</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Deng, Y" uniqKey="Deng Y">Y. Deng</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Apweiler, R" uniqKey="Apweiler R">R. Apweiler</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Kanehisa, M" uniqKey="Kanehisa M">M. Kanehisa</name></author><author><name sortKey="Goto, S" uniqKey="Goto S">S. Goto</name></author><author><name sortKey="Kawashima, S" uniqKey="Kawashima S">S. Kawashima</name></author><author><name sortKey="Okuno, Y" uniqKey="Okuno Y">Y. Okuno</name></author><author><name sortKey="Hattori, M" uniqKey="Hattori M">M. Hattori</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ashburner, M" uniqKey="Ashburner M">M. Ashburner</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Tatusov, R L" uniqKey="Tatusov R">R. L. Tatusov</name></author><author><name sortKey="Galperin, M Y" uniqKey="Galperin M">M. Y. Galperin</name></author><author><name sortKey="Natale, D A" uniqKey="Natale D">D. A. Natale</name></author><author><name sortKey="Koonin, E V" uniqKey="Koonin E">E. V. Koonin</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Li, R" uniqKey="Li R">R. Li</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Mortazavi, A" uniqKey="Mortazavi A">A. Mortazavi</name></author><author><name sortKey="Williams, B A" uniqKey="Williams B">B. A. Williams</name></author><author><name sortKey="Mccue, K" uniqKey="Mccue K">K. McCue</name></author><author><name sortKey="Schaeffer, L" uniqKey="Schaeffer L">L. Schaeffer</name></author><author><name sortKey="Wold, B" uniqKey="Wold B">B. Wold</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Audic, S" uniqKey="Audic S">S. Audic</name></author><author><name sortKey="Claverie, J M" uniqKey="Claverie J">J.-M. Claverie</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Shang, F" uniqKey="Shang F">F. Shang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Livak, K J" uniqKey="Livak K">K. J. Livak</name></author><author><name sortKey="Schmittgen, T D" uniqKey="Schmittgen T">T. D. Schmittgen</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Luan, J B" uniqKey="Luan J">J. B. Luan</name></author><author><name sortKey="Ghanim, M" uniqKey="Ghanim M">M. Ghanim</name></author><author><name sortKey="Liu, S S" uniqKey="Liu S">S. S. Liu</name></author><author><name sortKey="Czosnek, H" uniqKey="Czosnek H">H. Czosnek</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Jia, S" uniqKey="Jia S">S. Jia</name></author><author><name sortKey="Wan, P J" uniqKey="Wan P">P.-J. Wan</name></author><author><name sortKey="Zhou, L T" uniqKey="Zhou L">L.-T. Zhou</name></author><author><name sortKey="Mu, L L" uniqKey="Mu L">L.-L. Mu</name></author><author><name sortKey="Li, G Q" uniqKey="Li G">G.-Q. Li</name></author></analytic></biblStruct></listBibl></div1></back></TEI><pmc article-type="research-article"><pmc-dir>properties open_access</pmc-dir>
<front><journal-meta><journal-id journal-id-type="nlm-ta">Sci Rep</journal-id><journal-id journal-id-type="iso-abbrev">Sci Rep</journal-id><journal-title-group><journal-title>Scientific Reports</journal-title></journal-title-group><issn pub-type="epub">2045-2322</issn><publisher><publisher-name>Nature Publishing Group</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="pmid">27577531</article-id><article-id pub-id-type="pmc">5006003</article-id><article-id pub-id-type="pii">srep32099</article-id><article-id pub-id-type="doi">10.1038/srep32099</article-id><article-categories><subj-group subj-group-type="heading"><subject>Article</subject></subj-group></article-categories><title-group><article-title>Differential expression of genes in the alate and apterous morphs of the brown citrus aphid, <italic>Toxoptera citricida</italic></article-title></title-group><contrib-group><contrib contrib-type="author"><name><surname>Shang</surname><given-names>Feng</given-names></name><xref ref-type="aff" rid="a1">1</xref></contrib><contrib contrib-type="author"><name><surname>Ding</surname><given-names>Bi-Yue</given-names></name><xref ref-type="aff" rid="a1">1</xref></contrib><contrib contrib-type="author"><name><surname>Xiong</surname><given-names>Ying</given-names></name><xref ref-type="aff" rid="a1">1</xref></contrib><contrib contrib-type="author"><name><surname>Dou</surname><given-names>Wei</given-names></name><xref ref-type="aff" rid="a1">1</xref></contrib><contrib contrib-type="author"><name><surname>Wei</surname><given-names>Dong</given-names></name><xref ref-type="aff" rid="a1">1</xref></contrib><contrib contrib-type="author"><name><surname>Jiang</surname><given-names>Hong-Bo</given-names></name><xref ref-type="aff" rid="a1">1</xref></contrib><contrib contrib-type="author"><name><surname>Wei</surname><given-names>Dan-Dan</given-names></name><xref ref-type="aff" rid="a1">1</xref></contrib><contrib contrib-type="author"><name><surname>Wang</surname><given-names>Jin-Jun</given-names></name><xref ref-type="corresp" rid="c1">a</xref><xref ref-type="aff" rid="a1">1</xref></contrib><aff id="a1"><label>1</label><institution>Key Laboratory of Entomology and Pest Control Engineering, College of Plant Protection, Southwest University</institution>, Chongqing 400716,<country>China</country></aff></contrib-group><author-notes><corresp id="c1"><label>a</label><email>wangjinjun@swu.edu.cn</email></corresp></author-notes><pub-date pub-type="epub"><day>31</day><month>08</month><year>2016</year></pub-date><pub-date pub-type="collection"><year>2016</year></pub-date><volume>6</volume><elocation-id>32099</elocation-id><history><date date-type="received"><day>03</day><month>03</month><year>2016</year></date><date date-type="accepted"><day>02</day><month>08</month><year>2016</year></date></history><permissions><copyright-statement>Copyright © 2016, The Author(s)</copyright-statement><copyright-year>2016</copyright-year><copyright-holder>The Author(s)</copyright-holder><license license-type="open-access" xlink:href="http://creativecommons.org/licenses/by/4.0/"><pmc-comment>author-paid</pmc-comment>
<license-p>This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit <ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/4.0/">http://creativecommons.org/licenses/by/4.0/</ext-link></license-p></license></permissions><abstract><p>Winged and wingless morphs in insects represent a trade-off between dispersal ability and reproduction. We studied key genes associated with apterous and alate morphs in <italic>Toxoptera citricida</italic> (Kirkaldy) using RNAseq, digital gene expression (DGE) profiling, and RNA interference. The <italic>de novo</italic> assembly of the transcriptome was obtained through Illumina short-read sequencing technology. A total of 44,199 unigenes were generated and 27,640 were annotated. The transcriptomic differences between alate and apterous adults indicated that 279 unigenes were highly expressed in alate adults, whereas 5,470 were expressed at low levels. Expression patterns of the top 10 highly expressed genes in alate adults agreed with wing bud development trends. Silencing of the lipid synthesis and degradation gene (3-ketoacyl-CoA thiolase, mitochondrial-like) and glycogen genes (Phosphoenolpyruvate carboxykinase [GTP]-like and Glycogen phosphorylase-like isoform 2) resulted in underdeveloped wings. This suggests that both lipid and glycogen metabolism provide energy for aphid wing development. The large number of sequences and expression data produced from the transcriptome and DGE sequencing, respectively, increases our understanding of wing development mechanisms.</p></abstract></article-meta></front><body><p>Many insect species have dispersing and non-dispersing morphs. These include short-winged and long-winged morphs in planthoppers<xref ref-type="bibr" rid="b1">1</xref>, crickets<xref ref-type="bibr" rid="b2">2</xref>, migratory locusts<xref ref-type="bibr" rid="b3">3</xref>, and wingless and winged morphs in aphids<xref ref-type="bibr" rid="b4">4</xref><xref ref-type="bibr" rid="b5">5</xref>. These phenotypes are associated with different dispersal capabilities. Typically, alate and apterous morphs develop in response to specific environmental conditions<xref ref-type="bibr" rid="b6">6</xref>. Wing dimorphism is clearly an adaptation to fluctuating environmental conditions, increasing the dispersal capability of the insects and allowing them to escape unfavorable conditions<xref ref-type="bibr" rid="b6">6</xref>. Many studies have focused on the abiotic and biotic factors influencing aphid wing dimorphism. For instance, increased population density triggers wing formation in most aphid species and in some species, a relatively small density increase is sufficient<xref ref-type="bibr" rid="b7">7</xref><xref ref-type="bibr" rid="b8">8</xref><xref ref-type="bibr" rid="b9">9</xref>. A decrease in plant quality can trigger wing induction in the bird cherry-oat aphid, <italic>Rhopalosiphum padi</italic><xref ref-type="bibr" rid="b10">10</xref>. Relatively higher temperatures favor wingless forms while lower temperatures facilitate wing induction. Studies on starvation combined with symbiosis indicate a significant role of symbionts in wing dimorphism of the English grain aphid, <italic>Sitobion avenae</italic><xref ref-type="bibr" rid="b11">11</xref>.</p><p>Wing dimorphism represents a life history trade-off between dispersal and reproduction<xref ref-type="bibr" rid="b1">1</xref><xref ref-type="bibr" rid="b12">12</xref><xref ref-type="bibr" rid="b13">13</xref><xref ref-type="bibr" rid="b14">14</xref>. In aphids, alate and apterous morphs differ in presence or absence of wings and also differ in body structure, fecundity, longevity and behavior. Flight morphs have heavier sclerotization of the head and thorax, more fully developed compound eyes and ocelli, longer antennae, more rhinaria, and sometimes larger siphunculi and cauda<xref ref-type="bibr" rid="b6">6</xref><xref ref-type="bibr" rid="b15">15</xref><xref ref-type="bibr" rid="b16">16</xref>. The morphological differences between apterous and alate morphs usually correlate with differences in life history. In general, the alate morph, compared to the apterous morph, has a longer nymphal development period, longer pre-reproductive adult period, lower offspring production, and greater longevity<xref ref-type="bibr" rid="b17">17</xref><xref ref-type="bibr" rid="b18">18</xref>. Alate morphs are also more resistant to starvation and they have acute sensory capability for detecting host plants<xref ref-type="bibr" rid="b19">19</xref>.</p><p>The brown citrus aphid, <italic>Toxoptera citricida</italic> (Kirkaldy) (Hemiptera: Aphididae), is an important citrus pest and the main vector of <italic>Citrus tristeza virus</italic> (CTV) worldwide. Like other aphids, <italic>T. citricida</italic> has wing dimorphism. The apterous morph is about 3.2 mm in body length, is shiny black, lacks wing remnants<xref ref-type="bibr" rid="b20">20</xref><xref ref-type="bibr" rid="b21">21</xref>, and possesses higher fecundity. The alate morph is smaller than the apterous morph, has 5-mm-long wings<xref ref-type="bibr" rid="b20">20</xref><xref ref-type="bibr" rid="b21">21</xref>, and strong flight muscles. This allows it to fly long distances with the wind and to spread CTV in citrus growing regions. These differences are likely associated with the gene expression of alate and apterous morphs. Unlike the pea aphid, <italic>Acythosiphon pisum</italic><xref ref-type="bibr" rid="b22">22</xref>, and Russian wheat aphid, <italic>Diuraphis noxia</italic><xref ref-type="bibr" rid="b23">23</xref>, the genome of <italic>T. citricida</italic> has not been sequenced or released. Hunter <italic>et al</italic>.<xref ref-type="bibr" rid="b24">24</xref> completed expressed sequence tag (EST) sequencing of alate <italic>T. citricida</italic>. About 5,180 cDNA clones were sequenced resulting in 4,263 ESTs<xref ref-type="bibr" rid="b24">24</xref>, but this was insufficient for gene identification and prediction of functions. Fortunately, the emergence of next-generation sequencing (high-throughput sequencing) technology has dramatically enhanced the efficiency and quantity of gene annotation. In insects, transcriptome sequencing combined with digital gene expression (DGE) analysis is a reliable and precise way to investigate transcriptomic characteristics, such as insecticide targets<xref ref-type="bibr" rid="b25">25</xref>, immune response<xref ref-type="bibr" rid="b26">26</xref>, and chemoreception<xref ref-type="bibr" rid="b27">27</xref>. RNA interference (RNAi) methods have been used to study the function of genes involved in insect molting, growth, and development<xref ref-type="bibr" rid="b28">28</xref><xref ref-type="bibr" rid="b29">29</xref><xref ref-type="bibr" rid="b30">30</xref>. Thus, we expected that the transcriptome and DGE analysis of alate and apterous adults combined with RNAi would greatly improve understanding of molecular level differences between alate and apterous morphs, especially for the genes associated with dispersion of the alate morph. Here we reported: (1) transcriptome analysis covering different developmental stages of <italic>T. citricida</italic>, (2) DGE profiling of alate and apterous adults, (3) validation of the fold change of the top ten highly expressed genes in alate adults by quantitative real-time PCR (qRT-PCR) and examination of gene expression levels in different developmental stages, wing morphs, and body parts, and (4) functional analysis, using RNAi, of the energy metabolism genes which may be involved in wing development. The results provide a resource for future functional studies on wing dimorphism and development of <italic>T. citricida</italic> and other aphids.</p><sec disp-level="1"><title>Results</title><sec disp-level="2"><title>Transcriptome sequencing and assembly</title><p>A library of different developmental stages of <italic>T. citricida</italic> was developed by Illumina sequencing in a single run, which generated 41,250,684 raw reads, containing about 8.33 Gb sequencing data. After initial adaptor trimming and quality filtering, 1,199,882 contigs with an N50 length of 101 bp were assembled from 37,182,535 clean reads. The sequencing quality of the clean reads was evaluated based on the base-calling quality scores of Illumina’s base-caller Bustard software. More than 96% of the clean reads had quality scores higher than the Q<sub>30</sub> level (an error probability of 0.1%) (see <xref ref-type="supplementary-material" rid="S1">Supplementary Fig. S1</xref>). The contigs were further assembled into 44,199 unigenes with a mean length of 749 bp by using paired-end joining and gap-filling methods (<xref ref-type="table" rid="t1">Table 1</xref>). Among the total unigenes, 9,468 unigenes (21.4%) were longer than 1,000 bp and 26,443 unigenes (58.8%) were among 200 bp to 500 bp, and no unigene was less than 200 bp (see <xref ref-type="supplementary-material" rid="S1">Supplementary Fig. S2</xref>).</p></sec><sec disp-level="2"><title>Annotation of predicted transcripts</title><p>Unigene sequences were annotated using BLASTX against the non-redundant (NR) NCBI protein database. The E-value distribution of the annotated unigenes showed that 53.5% of the sequences had high homology (<1.0 E-50), while 46.5% of the sequences ranged between 1.0 E-5 to 1.0 E-50 (<xref ref-type="fig" rid="f1">Fig. 1a</xref>). Further analysis of the similarity distribution indicated that 50.8% of the sequences had similarity higher than 80%, 21.9% of the sequences with a similarity between 60% to 80%, and only 7.4% of the sequences possessed a similarity lower than 40% (<xref ref-type="fig" rid="f1">Fig. 1b</xref>). According to the best hits in the NR database, about half of annotated unigenes (12,761 unigenes, 49.1%) exhibited homology with the pea aphid, <italic>A. pisum</italic>, and Russian wheat aphid, <italic>D. noxia</italic>, which reflected the close relationship among <italic>T. citricida</italic>, <italic>A. pisum</italic> and <italic>D. noxia,</italic> while 28.4% of them matched those in other species (<xref ref-type="fig" rid="f1">Fig. 1c</xref>).</p></sec><sec disp-level="2"><title>Unigene function annotation and metabolic pathway analysis</title><p>Assignments of Clusters of Orthologous Groups (COG) were used to predict and classify possible functions of the unigenes. Based on sequence homology, 11,052 unigenes were annotated and divided into 26 specific categories (see <xref ref-type="supplementary-material" rid="S1">Supplementary Fig. S3</xref>). The general function category, which contained 2,759 unigenes (25.0%), was the largest, followed by translation, ribosomal structure, and biogenesis (1,283, 11.6%), replication, recombination and repair (1,242, 11.2%), and transcription (1,196, 10.8%), respectively. Only 1 unigene (less than 1%) belonged to extracellular structures, which was the smallest group.</p><p>For Gene Ontology (GO) analysis, we categorized 12,189 unigenes into 59 function groups. Among these, metabolic process, cellular process, and catalytic activity were the 3 largest groups, containing 8,254, 7,361, and 7,164 unigenes, respectively. Only 1 unigene each was predicted to act in the component and function groups including protein tag, nutrient reservoir activity, morphogen activity, cell death, and carbon utilization (see <xref ref-type="supplementary-material" rid="S1">Supplementary Fig. S4</xref>).</p><p>We further mapped 10,522 unigenes to 202 pathways in the Kyoto Encyclopedia of Genes and Genomes (KEGG) annotation system. Metabolic pathways containing 2,180 unigenes (20.7%) were significantly larger than other pathways, such as protein processing in the endoplasmic reticulum (446, 4.2%), RNA transport (415, 3.9%), and spliceosome (403, 3.8%; see <xref ref-type="supplementary-material" rid="S1">Supplementary Table S1</xref>).</p></sec><sec disp-level="2"><title>Summary of digital gene expression sequencing</title><p>DGE analysis was performed to obtain a global view of the aphid transcriptome differences alate and apterous adults. To compare differentially expressed genes between the libraries (alate/apterous), the level of gene expression was determined by normalizing the number of unambiguous tags in each library to reads per kilobase mapped (RPKM). DGE gave high reproducibility among biological replicates, as indicated by the Pearson correlation coefficient (<italic>r</italic>) ≥ 0.98 (Alate1 versus Alate 2: 0.9925, Apterous 1 versus Apterous 2: 0.9881; see <xref ref-type="supplementary-material" rid="S1">Supplementary Fig. S5</xref>). In comparison with apterous adults, the expression levels of 5,749 unigenes was significantly different in alate adults (see <xref ref-type="supplementary-material" rid="S1">Supplementary Fig. S6</xref>) including 279 genes with higher expression (see <xref ref-type="supplementary-material" rid="S1">Supplementary Table S2</xref>) and 5,470 genes with lower expression (see <xref ref-type="supplementary-material" rid="S1">Supplementary Table S3</xref>). The fold changes (log<sub>2</sub>Ratio) was calculated based on the gene expression level of alate divided by that of apterous adults.</p><p>Using the GO classification, these differentially expressed genes were characterized into 3 groups: biological process, cellular component and molecular function (see <xref ref-type="supplementary-material" rid="S1">Supplementary Fig. S7</xref>). Among the biological process assignments, a high percentage of genes were concentrated in the cellular process category (1,024 genes, 19.2%), predominantly in the microtubule-based movement, spindle organization, and metabolic process category (969 genes, 18.2%). The cellular component group showed a significant percentage of genes assigned to cell part (1,074 genes, 25.3%) and organelle (737 genes, 17.3%) categories, where they were mainly concentrated in the microtubule dynein complex and mitochondrial respiratory chain.</p><p>KEGG pathway analysis of the highly expressed transcripts in alate adults was performed. Among the highly expressed transcripts, most genes could be annotated in KEGG orthology (KO) based on sequence homologies and mapped to 60 pathways, and the top 20 enriched metabolic pathways are listed in <xref ref-type="table" rid="t2">Table 2</xref>. Pathway-oriented analysis showed that most pathways were involved in energy resources and nutrient-sensing pathways, such as the citrate cycle (TCA cycle), fatty acid metabolism, fatty acid elongation, insulin signaling pathway, glyoxylate and dicarboxylate metabolism, and glycolysis/gluconeogenesis.</p></sec><sec disp-level="2"><title>qRT-PCR validation</title><p>To confirm the DGE results by Illumina sequencing, we selected the top 10 highly expressed genes in alate adults (Sodium-and chloride-dependent creatine transporter 1-like isoform 1 (<italic>SCDCT</italic>), 3-ketoacyl-CoA thiolase, mitochondrial-like (<italic>KAT</italic>), ATP-binding cassette sub-family G member 4-like isoform 1 (<italic>ABCG4</italic>), trifunctional enzyme subunit beta, mitochondrial-like (<italic>β-TE</italic>), short-chain specific acyl-CoA dehydrogenase, mitochondrial-like (<italic>SCAD</italic>), trifunctional enzyme subunit alpha, mitochondrial-like (<italic>α-TE</italic>), phosphoenolpyruvate carboxykinase [GTP]-like (<italic>PEPCK</italic>), hydroxysteroid dehydrogenase-like protein 2-like (<italic>HSD</italic>), E3 ubiquitin-protein ligase TRIM23-like (<italic>E3-TRIM23</italic>), glycogen phosphorylase-like isoform 2 (<italic>GP</italic>) for qRT-PCR confirmation. Their annotations and fold changes are summarized in <xref ref-type="table" rid="t3">Table 3</xref>. All 10 genes showed similar change trends as indicated by the DGE analysis (see <xref ref-type="supplementary-material" rid="S1">Supplementary Fig. S8</xref>), suggesting that the DGE approach can be a powerful tool for the discovery of candidate genes in a high-throughput manner. However, there were some fold change magnitude differences between qRT-PCR experiments and DGE analysis. For example, we found 37-fold in qRT-PCR versus 8-fold in DGE for <italic>SCDCT</italic> and 16-fold in qRT-PCR versus 5-fold in DGE for <italic>PEPCK</italic>.</p><p>The expression patterns of the top 10 highly expressed genes in alate adults were further determined in different developmental stages, wing morphs, and body parts, respectively.</p><p>All gene expression levels were sharply increased from fourth instar nymphs to adults and most showed no difference among the nymphal stages of alates. No significant differences existed among nymphal stages and adult in apterous aphids (<xref ref-type="fig" rid="f2">Fig. 2</xref>). The wing bud developed slowly in each nymphal instar until the fully formed wings unfolded after adult emergence. This means that the fourth instar nymph to adult transition was the key period of aphid wing development. The gene expression patterns agree with the trend in wing development. Thus, these genes might be involved in wing development.</p><p>Head, thorax, abdomen and wing were detached from fourth instar winged-nymphs, alate adults and apterous adults for expression pattern analyses. The results showed that all the genes were highly expressed in the thorax, and about 5-fold greater compared to other body parts of <italic>SCAD</italic> and <italic>HSD</italic> in alate adults. The changes were different for wing-related genes between fourth instar winged-nymphs and alate adults. <italic>ABCG4</italic>, <italic>PEPCK</italic>, <italic>GP</italic> were highly expressed in alate adults but <italic>KAT</italic>, <italic>SCAD</italic>, <italic>α-TE</italic>, <italic>β-TE</italic>, <italic>E3-TRIM23</italic>, and <italic>HSD</italic> exhibited higher expression levels in fourth instar winged-nymphs (<xref ref-type="fig" rid="f3">Fig. 3</xref>), suggesting the potential crucial roles of the genes in wing development.</p></sec><sec disp-level="2"><title>Functional analyses</title><p>RNAi was used to study the lipid synthesis/degradation gene (<italic>KAT</italic>) and glycogen related genes (<italic>PEPCK</italic> and <italic>GP</italic>), and their effects on wing development. Aphids were synchronized as new fourth instar winged-nymphs and fed with dsRNA. After 72 h, all aphids normally molt from fourth instar nymphs to adults, but RNAi treated aphids showed a variety of different phenotypes including those with under-developed, malformed, and normal wings. Because both malformed and under-developed wing aphids were unable to fly and the malformed wing aphid appeared at a very low rate (<5%), we categorized this type into the under-developed wing when counting the phenotype and collecting samples (<xref ref-type="fig" rid="f4">Fig. 4a</xref>). Silencing of these 3 genes resulted in under-developed wing aphids at rates of 58% (ds<italic>KAT</italic>) to 79% (ds<italic>PEPCK</italic>) compared to ds<italic>GFP</italic>, which was 100% normal (<xref ref-type="fig" rid="f4">Fig. 4b</xref>). When challenged at 72 h, the relative expression levels of these 3 genes in under-developed wing aphids were down-regulated by 50% (<italic>GP</italic>), 74% (<italic>PEPCK</italic>), and 95% (<italic>KAT</italic>). The relative expression levels of <italic>KAT</italic> (93%), <italic>PEPCK</italic> (75%) were decreased in normal wing aphids by feeding with dsRNA compared to ds<italic>GFP</italic> but the difference was not significant in the ds<italic>GP</italic> treatment (<xref ref-type="fig" rid="f4">Fig. 4c</xref>). These results indicate that the gene expression levels in under-developed wing aphids were all lower than normal when the aphids were treated with dsRNA.</p></sec></sec><sec disp-level="1"><title>Discussion</title><p>Flight capacity in insects is a key feature behind their ecological and evolutionary success. Flight benefits dispersal capacity balanced against potential metabolic, reproductive and survival costs<xref ref-type="bibr" rid="b14">14</xref><xref ref-type="bibr" rid="b31">31</xref>. In this study, we used the transcriptomic tools including RNA-seq and DGE, and RNAi technology, to analyze the molecular mechanisms of wing development of <italic>T. citricida</italic>.</p><p>We generated <italic>de novo</italic> assembly of the <italic>T. citricida</italic> transcriptome through short read sequencing technology (Illumina). In a single run, we obtained more than 41 million sequencing reads that were assembled into 44,199 unigenes by Trinity. Homology analysis of the unigenes demonstrated that only 49.1% showed the greatest similarity to <italic>A. pisum</italic> and <italic>D. noxia</italic>, the only 2 aphid species with a completely sequenced genome<xref ref-type="bibr" rid="b22">22</xref><xref ref-type="bibr" rid="b23">23</xref>. This rate is smaller than expected. <italic>A. pisum</italic> and <italic>D. noxia</italic> belong to Macrosiphini, but <italic>T. citricida</italic> belongs to Aphidini, and the host range of <italic>T. citricida</italic> is limited to citrus and a few close relatives<xref ref-type="bibr" rid="b32">32</xref><xref ref-type="bibr" rid="b33">33</xref>. The new transcriptome data reported here complements the scarce aphid sequence resources in the GenBank. We found that 12.3% of unigenes matched those in <italic>Microplitis demolitor</italic> and <italic>Nasonia vitripennis</italic>, which are endoparasitoids of Lepidopteran and Dipteran<xref ref-type="bibr" rid="b34">34</xref><xref ref-type="bibr" rid="b35">35</xref>, and we also found that 1.5% of the unigenes matched those in the fire ant <italic>Solenopsis invicta</italic>. Our laboratory will be focusing on this interesting phenomenon to explore the potential horizontal gene transfer among endoparasitoids, aphids, and ants.</p><p>The transcriptomic differences between alate and apterous adults were examined by DGE to explore the molecular mechanism associated with wing morphs, especially the genes influencing the development or behavior of alates. Through DGE analysis, we produced a total of 5,749 differentially expressed transcripts, of which 279 transcripts were highly expressed in alate adults. The highly expressed transcripts in alate adults belonged to 60 KEGG pathways, and most of the pathways were involved in energy resource and nutrient-sensing. The insulin signaling pathway was interesting pathway owing to its role in energy distribution, and this is an evolutionarily conserved nutrient-sensing pathway that modulates growth and development in metazoans<xref ref-type="bibr" rid="b36">36</xref><xref ref-type="bibr" rid="b37">37</xref>. In rice brown planthopper, <italic>Nilaparvata lugens</italic>, silencing of 2 insulin receptor genes (<italic>NlInR1</italic> and <italic>NlInR2</italic>) by RNAi was used to understand the regulation of wing dimorphism by the insulin/insulin-like signaling pathway<xref ref-type="bibr" rid="b38">38</xref>. This was the first evidence of a molecular basis for the regulation of wing polyphenism in a planthopper. The current results also implied that insulin/insulin-like signaling pathway maybe play an crucial role in alate aphid.</p><p>The qRT-PCR experiment confirmed the DGE results. All of the top 10 highly expressed genes in alate adults subjected to the qRT-PCR analysis showed changes similar to the DGE results but the fold differences varied. The discrepancies in the fold changes between qRT-PCR and DGE may have resulted from a sensitivity bias between the 2 methods or by virtue of different statistical methods and threshold values in qRT-PCR and DGE<xref ref-type="bibr" rid="b39">39</xref>. The analyses of gene expression patterns showed that almost all genes had low expression levels in nymphal stages but sharply increased in the transition from fourth instar winged nymphs (N4-WD) molting to alate adults. This result was consistent with findings that the transition period between the final instar winged-nymph to alate adults is the key stage of wing development<xref ref-type="bibr" rid="b40">40</xref>. The results indicated that these genes may be involved in wing development, and this was subsequently verified by dsRNA feeding mediated RNAi. Further analysis of gene expression levels among different body parts of N4-WD, alate, and apterous adults showed that all of the genes had their highest expression level in thorax except <italic>HSD</italic> suggesting that these genes may be involved in the formation or development of flight muscle<xref ref-type="bibr" rid="b41">41</xref><xref ref-type="bibr" rid="b42">42</xref>. We also found that gene expression differed in wings of N4-WD and alate adults, suggesting that they may play positive or negative roles in wing development<xref ref-type="bibr" rid="b30">30</xref>.</p><p>Of the 279 highly expressed unigenes in alate adults, the top 10 highly expressed genes were chosen to clarify their potential functions. <italic>KAT</italic>, <italic>SCAD</italic>, <italic>α-TE</italic>, and <italic>β-TE</italic> were “mitochondrial like” genes that may contribute to energy production and are thought to be involved in aspects of the citric acid cycle respiratory chain, ATP synthesis, and transport<xref ref-type="bibr" rid="b5">5</xref>. In <italic>Aphis gossypii</italic> and <italic>A. pisum</italic>, very long-chain specific, medium-chain, and short-chain specific acyl-CoA dehydrogenase, and <italic>KAT</italic> were highly expressed in alate adults. These genes, associated with lipid synthesis and degradation and genes involved in sugar metabolism are elevated in apterous adults, indicating that the aphids use lipid rather than sugar or glycogen as fuel during flight<xref ref-type="bibr" rid="b5">5</xref><xref ref-type="bibr" rid="b43">43</xref>. Of particular interest is that, in addition to the lipid synthesis/degradation related gene (<italic>KAT</italic>), <italic>PEPCK</italic>, and <italic>GP</italic> were also highly expressed in alate adults and these are linked with glycogenolysis. RNAi was further used to explore the function of lipid synthesis/degradation and glycogen related genes and their effects on aphid wing development. Silencing of <italic>KAT</italic>, <italic>PEPCK</italic>, and <italic>GP</italic> all resulted in under-developed wings and flight disruption of <italic>T. citricida</italic>. We found that the relative expression levels of <italic>KAT</italic> and <italic>PEPCK</italic> were significantly decreased in normal alates when challenged with dsRNA of target genes compared to ds<italic>GFP</italic>, and the relative expression levels of the target genes were also significantly different between alates with normal wings or under-developed wing for dsRNA treatments. One possible explanation is that the development period may vary among individuals from the fourth instar winged-nymph to the alate adult and the key period for wing development also may vary. Thus, although the genes were silenced, the phenotype remained unchanged. Another possibility is that there exists feedback regulation so when expression of <italic>KAT</italic> or <italic>PEPCK</italic> is decreased, the other gene might be able to keep the wing normal but when gene expression decreased below a certain level, this feedback will fail to work. The results suggested that both lipid and glycogen provide resources for wing development and dispersion in <italic>T. citricida</italic>, but the relative importance of these in dispersion is unclear. The role of lipid and glycogen during migration should be studied. Additional functions of these genes or their functions in apterous morphs require more study.</p><p>In conclusion, this study is the first report of genetic information on <italic>T. citricida</italic> from sequenced transcriptome and constructed DGE libraries. The data revealed a large number of genes, with both known and unknown functions, which greatly enrich the sequence information of aphids. In addition, we also identified genes that are candidates related to wing development in <italic>T. citricida</italic>. The RNAi results predict the possible function of lipid synthesis/degradation and glycogen related genes in wing development.</p></sec><sec disp-level="1"><title>Methods</title><sec disp-level="2"><title>Insect culture</title><p><italic>T. citricida</italic> used in this study were established in 2012 from a single collection of wild alate aphids from a screen house at Southwest University, Chongqing, China. Stock colonies were maintained on potted citrus seedlings (<italic>Citrus sinensis</italic>) in the laboratory at 25 ± 1 °C, 75–80% relative humidity, and a photoperiod of 14: 10 h (Light: Dark). Alate morphs were induced by high-density aphid rearing after transfer to fresh host plants<xref ref-type="bibr" rid="b9">9</xref><xref ref-type="bibr" rid="b24">24</xref>. All progeny were produced by parthenogenesis from the stock colony.</p></sec><sec disp-level="2"><title>RNA isolation, cDNA library preparation and Illumina sequencing</title><p>Thirty first instar nymphs (N1), second instar nymphs (N2), third instar wingless and winged nymphs (N3-WL and N3-WD), fourth instar wingless and winged nymphs (N4-WL and N4-WD), apterous adults (AP, with mature embryos) and alate adults (AL, with mature embryos) were all collected at about 14:00 of the day for total RNA isolation, respectively. Apterous and alate adults were collected separately within 48 h after the fourth instar nymph molt. Total RNA was isolated with a TRIzol kit (Invitrogen, Carlsbad, CA) according to manufacturer instructions. RNA was quantified by measuring absorbance at 260 nm using a NanoVue UV-Vis spectrophotometer (GE Healthcare Bio-Science, Uppsala, Sweden). The purity of all RNA samples was assessed at an absorbance ratio of OD<sub>260/280</sub> and OD<sub>260/230</sub>, and the integrity of RNA was checked on 1% agarose gel by electrophoresis. DNase I (Promega, Madison, WI) was used to remove genomic DNA from the samples.</p><p>The cDNA library construction and Illumina sequencing of the samples were performed at Beijing Biomarker Technologies, Beijing, China. Briefly, 12 μg of the RNA mixture (each development stages and wing morph were mixed at equal ratio) was used for cDNA library construction. Oligo (dT) magnetic beads were used to isolate poly (A) mRNA, which was then fragmented into small pieces by the addition of fragmentation buffer. These short fragments served as templates to synthesize first-strand cDNA using random hexamer-primers. Second-strand cDNA synthesis was performed using a buffer, dNTPs, RNaseH, and DNA polymerase I. Short fragments were purified using a QiaQuick PCR extraction kit (Qiagen) and then they were washed with ethidium bromide buffer for end preparation and single nucleotide adenine addition before being ligated to adaptors. Suitable fragments, as judged by agarose gel electrophoresis, were selected as templates for PCR amplification. The cDNA library was sequenced on Illumina Hiseq<sup>TM</sup> 2000 using paired-end technology in a single run. Through sequencing, a total of 41,250,684 pair-end reads and 8,332,638,168 bases were achieved and the raw reads were submitted to the NCBI Short Read Archive (SRA) database with accession number of SRR2123649.</p></sec><sec disp-level="2"><title><italic>De novo</italic> assembly and bioinformatics analysis</title><p>The raw reads were cleaned by removing adapter sequences, reads containing ambiguous bases, and low-quality reads. The sequencing quality of the clean reads was evaluated based on the base-calling quality scores of Illumina’s base-caller Bustard. The clean reads were assembled into contigs using the Trinity method (<ext-link ext-link-type="uri" xlink:href="https://github.com/trinityrnaseq/trinityrnaseq/wiki">https://github.com/trinityrnaseq/trinityrnaseq/wiki</ext-link>), which is efficient in reconstructing full-length transcripts across a broad range of expression levels and sequencing depths with an optimized k-mer length of 25 for <italic>de novo</italic> assembly<xref ref-type="bibr" rid="b44">44</xref>. Then the transcripts were clustered based on nucleotide sequence identity. If there was more than 1 transcript for a given gene, the longest transcript was used to calculate its expression level and coverage.</p><p>All assembled unigenes were determined by BLASTx against the NCBI non-redundant (NR) protein database<xref ref-type="bibr" rid="b45">45</xref>, the Swiss-Prot database<xref ref-type="bibr" rid="b46">46</xref>, the Kyoto Encyclopedia of Genes and Genomes (KEGG) database<xref ref-type="bibr" rid="b47">47</xref>, Gene Ontology (GO)<xref ref-type="bibr" rid="b48">48</xref> and Cluster of Orthologous Groups (COG)<xref ref-type="bibr" rid="b49">49</xref> with a cut-off <italic>E</italic>-value of 10<sup>−5</sup>. If the alignment results of different databases conflicted with each other, the priority order of NR, Swiss-Prot, KEGG, GO, and COG were followed. When a unigene did not align with any of the entries in these databases, ESTs were used to predict its coding regions and to determine its sequence direction.</p></sec><sec disp-level="2"><title>DGE library preparation and sequencing</title><p>A total of 100 apterous or alate adults (2 biological replicates for each morph and 50 aphids in each biological replicate) were collected separately within 48 h after fourth instar nymphs molting at the same time of the day according to the above description. RNA extraction was done using TRIzol kit in according to manufacturer instruction, as described above, and DNase I (Promega, Madison, WI) was also used to remove genomic DNA from the samples. Approximately 30 μg RNA from each sample was used to construct the DGE libraries. The mRNA was treated as described in cDNA library construction and enriched by PCR amplification. The library products were then ready for sequencing analysis via Illumina Hiseq<sup>TM</sup> 2000 (Beijing Genomics Institute, BGI, Shenzhen, China) using paired-end technology in a single run.</p></sec><sec disp-level="2"><title>Annotation of DGE tags and analysis of differentially expressed genes</title><p>All raw sequence reads were filtered using the Illumina pipeline before mapping reads to the reference transcriptome databases. Low-quality reads were omitted from data analysis. Low-quality reads were defined as those with the percentage of ambiguous bases >10% and reads in which >50% of the bases had a quality value ≤ 5. Clean reads were mapped to reference sequences (unigenes from the transcriptome data of different developmental stages was used as a reference) using SOAPaligner/soap2<xref ref-type="bibr" rid="b50">50</xref>. No more than 2 mismatches were allowed in the alignment. The level of gene expression was determined by normalizing the number of unambiguous tags in each sample to RPKM (reads per kilobase mapped)<xref ref-type="bibr" rid="b51">51</xref>. The relative gene expression in alate or apterous adult was equal to the mean of two biological samples of each morph. The fold change was calculated as fold change = log<sub>2</sub>Ratio. The differentially expressed transcripts between samples were identified using an Bayesian algorithm developed by Audic and Claverie<xref ref-type="bibr" rid="b52">52</xref>. The false discovery rate (FDR) method was used to determine the threshold of the <italic>P</italic>- value in multiple tests and analyses. The threshold FDR ≤ 0.001 and the absolute value of log<sub>2</sub>Ratio ≥ 1 were used to judge the significance of differences in gene expression. Then the genes that were expressed at different levels across the samples were further annotated by GO function and KEGG pathway analysis.</p></sec><sec disp-level="2"><title>Quantitative real-time PCR (qRT-PCR) analysis</title><p>To confirm the DGE results, the top 10 highly expressed genes in alate adults were chosen and the gene-specific primers designed for targets genes were listed in <xref ref-type="supplementary-material" rid="S1">Table S4</xref> (see <xref ref-type="supplementary-material" rid="S1">Supplementary Table S4</xref>). To determine the expression profiles of these genes in different development stages, body parts, and wing morphs, we collected the following for RNA isolation: 30 of each developmental stage (first instar nymph, second instar nymph, third instar nymph, fourth instar nymph, and adult) of apterous and alate aphids; 30 of each of the body parts (head, thorax, and abdomen) from fourth instar winged-nymphs, apterous and alate adults; and 30 wing buds from fourth instar winged-nymphs and alate adults. RNA was extracted as described for the DGE library preparation and sequencing with 4 biological replicates. The first strand cDNA was synthesized from 500 ng of DNA-free RNA using the PrimerScript<sup>®</sup> RT Reagent Kit Perfect Real time Kit (Takara, Dalian, China) according to the manufacturer instructions. Briefly, all reaction volumes were 10 μL and contained 500 ng of RNA, 200 pmol of Random hexamers, 2 μL of reverse transcription buffer, 0.5 μL of PrimerScript<sup>®</sup> RT Enzyme Mix I and RNase free H<sub>2</sub>O. The reverse transcription reaction was performed on a C1000<sup>TM</sup> Thermal Cycler (Bio-Rad, Hercules, CA) at 37 °C for 15 min followed by 85 °C for 5 s. The resulting cDNA were stored at −20 °C until further use.</p><p>Expression reactions of all genes were performed using an Mx3000P thermal cycler (Agilent Technologies, Inc., Wilmington, NC). Each 20 μL reaction mixture containing 10 μL of iQ<sup>TM</sup> SYBR Green Supermix (BIO-RAD, Hercules, CA), 7 μL nuclease-free water, 1 μL of each gene-specific primer (0.2 mM) and 1 μL of the first strand cDNA. The PCR parameters were as follows: 95 °C for 120 s, then 40 cycles at 95 °C for 30 s and 60 °C for 30 s, a final cycle of 60 °C for 30 s and 95 °C for 30 s. The elongation factor-1 alpha (<italic>EF1α</italic>) was used as an internal control<xref ref-type="bibr" rid="b53">53</xref> and the relative expression of genes were calculated using the 2<sup>−△△<italic>Ct</italic></sup>method<xref ref-type="bibr" rid="b54">54</xref>. The data was analyzed using one-way ANOVA followed by Tukey’s honestly significant difference test (<italic>P</italic> < 0.05) to determine the significant differences among various developmental stages and body parts (head, thorax, and abdomen). Student’s <italic>t</italic>-test was used to separate the means between the wing of N4-WD and alate adults.</p></sec><sec disp-level="2"><title>dsRNA synthesis and plant-mediated RNAi</title><p>RNAi was used to study the functions of lipid synthesis and the degradation gene (<italic>KAT</italic>) and glycogen related genes (<italic>PEPCK</italic> and <italic>GP</italic>) and the effect on wing development in <italic>T. citricida</italic>. The unique nucleotide regions of the genes were selected for specific dsRNA synthesis. The primers used to synthesize dsRNA are listed in <xref ref-type="supplementary-material" rid="S1">Table S4</xref>. All of the reagents and enzymes used for the dsRNA synthesis were from the TranscriptAid T7 High Yield Transcription Kit (Thermo Scientific, Wilmington, DE). The size of the dsRNA products was confirmed by electrophoresis on a 1% agarose gel and the final concentration of dsRNA was 1500 ng/μL.</p><p>We developed a method to silence gene activity by dsRNA feeding through a citrus stem based on the methods of RNAi used in <italic>Bemisia tabaci</italic><xref ref-type="bibr" rid="b55">55</xref> and <italic>Sogatella furcifera</italic><xref ref-type="bibr" rid="b56">56</xref> with slight modification. Briefly, an 8-cm-long citrus stem was detached from the citrus seeding and inserted into a 250-μL PCR tube containing 200 μL of dsRNA. Then the tube containing the dsRNA and stem was transferred into a 50 ml plastic tube. Thirty fourth instar winged-nymphs were collected at the same time of day as above description and released onto the stem and the aphids were reared under conditions described above. After feeding with dsRNA for 72 h, morphological changes were observed and photos were taken using a Leica M165C microscope (Leica Microsystems, Wetzlar, Germany). Total RNA was isolated from surviving aphids after feeding for 72 h and the levels of transcripts were measured using qRT-PCR as described above. Four biological replicates were performed for each sample and ds<italic>GFP</italic> was used as a control. One-way ANOVA followed by Tukey’s honestly significant difference test (<italic>P</italic> < 0.05) was used to determine the significant differences among treatments.</p></sec></sec><sec disp-level="1"><title>Additional Information</title><p><bold>How to cite this article</bold>: Shang, F. <italic>et al</italic>. Differential expression of genes in the alate and apterous morphs of the brown citrus aphid, <italic>Toxoptera citricida</italic>. <italic>Sci. Rep.</italic><bold>6</bold>, 32099; doi: 10.1038/srep32099 (2016).</p></sec><sec sec-type="supplementary-material" id="S1"><title>Supplementary Material</title><supplementary-material id="d33e23" content-type="local-data"><caption><title>Supplementary Information</title></caption><media xlink:href="srep32099-s1.pdf"></media></supplementary-material><supplementary-material id="d33e26" content-type="local-data"><caption><title>Supplementary dataset</title></caption><media xlink:href="srep32099-s2.xls"></media></supplementary-material></sec></body><back><ack><p>This research was supported in part by the Special Fund for Agro-scientific Research in the Public Interest (201203038), the earmarked fund for Modern Agro-industry (Citrus) Technology Research System of China (CARS-27), the Fundamental Research Funds for the Central Universities of China (XDJK2016D039), and the Innovation Fund for Graduate Students of Chongqing (CYB16073).</p></ack><ref-list><ref id="b1"><mixed-citation publication-type="journal"><name><surname>Denno</surname><given-names>R. F.</given-names></name>, <name><surname>Olmstead</surname><given-names>K. L.</given-names></name> & <name><surname>Mccloud</surname><given-names>E.</given-names></name><article-title>Reproductive cost of flight capability: a comparison of life history traits in wing dimorphic planthoppers</article-title>. <source>Ecol. Entomol.</source><volume>14</volume>, <fpage>31</fpage>–<lpage>44</lpage> (<year>1989</year>).</mixed-citation></ref><ref id="b2"><mixed-citation publication-type="journal"><name><surname>Zhao</surname><given-names>Z.</given-names></name> & <name><surname>Zera</surname><given-names>A. J.</given-names></name><article-title>Differential lipid biosynthesis underlies a tradeoff between reproduction and flight capability in a wing-polymorphic cricket</article-title>. <source>P. Natl. Acad. Sci. USA</source><volume>99</volume>, <fpage>16829</fpage>–<lpage>16834</lpage> (<year>2002</year>).</mixed-citation></ref><ref id="b3"><mixed-citation publication-type="journal"><name><surname>Tanaka</surname><given-names>S.</given-names></name> & <name><surname>Nishide</surname><given-names>Y.</given-names></name><article-title>First record of the occurrence and genetics of a short-winged morph in the migratory locust, <italic>Locusta migratoria</italic> (Orthoptera: Acrididae)</article-title>. <source>J. Orthoptera Res.</source><volume>21</volume>, <fpage>169</fpage>–<lpage>174</lpage> (<year>2012</year>).</mixed-citation></ref><ref id="b4"><mixed-citation publication-type="journal"><name><surname>Xu</surname><given-names>X.</given-names></name>, <name><surname>Liu</surname><given-names>X.</given-names></name>, <name><surname>Zhang</surname><given-names>Q.</given-names></name> & <name><surname>Wu</surname><given-names>J.</given-names></name><article-title>Morph-specific differences in biochemical composition related to flight capability in the wing-polyphenic <italic>Sitobion avenae</italic></article-title>. <source>Entomol. Exp. Appl.</source><volume>138</volume>, <fpage>128</fpage>–<lpage>136</lpage> (<year>2011</year>).</mixed-citation></ref><ref id="b5"><mixed-citation publication-type="journal"><name><surname>Yang</surname><given-names>X.</given-names></name><etal></etal>. <article-title>Gene expression profiling in winged and wingless cotton aphids, <italic>Aphis gossypii</italic> (Hemiptera: Aphididae)</article-title>. <source>Int. J. Biol. Sci.</source><volume>10</volume>, <fpage>257</fpage> (<year>2014</year>).<pub-id pub-id-type="pmid">24644424</pub-id></mixed-citation></ref><ref id="b6"><mixed-citation publication-type="journal"><name><surname>Harrison</surname><given-names>R. G.</given-names></name><article-title>Dispersal polymorphisms in insects</article-title>. <source>Annu. Rev. Eco. Evol. S.</source><volume>11</volume>, <fpage>95</fpage>–<lpage>118</lpage> (<year>1980</year>).</mixed-citation></ref><ref id="b7"><mixed-citation publication-type="journal"><name><surname>Lees</surname><given-names>A.</given-names></name><article-title>The production of the apterous and alate forms in the aphid <italic>Megoura viciae</italic> Buckton, with special reference to the role of crowding</article-title>. <source>J. Insect Physiol.</source><volume>13</volume>, <fpage>289</fpage>–<lpage>318</lpage> (<year>1967</year>).</mixed-citation></ref><ref id="b8"><mixed-citation publication-type="journal"><name><surname>Sutherland</surname><given-names>O.</given-names></name><article-title>The role of crowding in the production of winged forms by two strains of the pea aphid</article-title>, <source>Acyrthosiphon pisum. J. Insect Physiol.</source><volume>15</volume>, <fpage>1385</fpage>–<lpage>1410</lpage> (<year>1969</year>).</mixed-citation></ref><ref id="b9"><mixed-citation publication-type="journal"><name><surname>Michaud</surname><given-names>J.</given-names></name><article-title>Colony density and wing development in <italic>Toxoptera citricida</italic> (Homoptera: Aphididae)</article-title>. <source>Environ. Entomol.</source><volume>30</volume>, <fpage>1047</fpage>–<lpage>1051</lpage> (<year>2001</year>).</mixed-citation></ref><ref id="b10"><mixed-citation publication-type="journal"><name><surname>De Barro</surname><given-names>P.</given-names></name><article-title>The role of temperature, photoperiod, crowding and plant quality on the production of alate viviparous females of the bird cherry-oat aphid</article-title>, <source>Rhopalosiphum padi. Entomol. Exp. Appl.</source><volume>65</volume>, <fpage>205</fpage>–<lpage>214</lpage> (<year>1992</year>).</mixed-citation></ref><ref id="b11"><mixed-citation publication-type="journal"><name><surname>Zhang</surname><given-names>F.</given-names></name><etal></etal>. <article-title>Bacterial symbionts, <italic>Buchnera</italic>, and starvation on wing dimorphism in English grain aphid, <italic>Sitobion avenae</italic> (F.) (Homoptera: Aphididae)</article-title>. <source>Front. Physiol.</source><volume>6</volume>, <fpage>155</fpage> (<year>2015</year>).<pub-id pub-id-type="pmid">26042046</pub-id></mixed-citation></ref><ref id="b12"><mixed-citation publication-type="journal"><name><surname>Andersen</surname><given-names>N. M.</given-names></name><article-title>Seasonal polymorphism and developmental changes in organs of flight and reproduction in <italic>Bivoltine pondskaters</italic> (Hem. Gerridae)</article-title>. <source>Insect Syst. Evol.</source><volume>4</volume>, <fpage>1</fpage>–<lpage>20</lpage> (<year>1973</year>).</mixed-citation></ref><ref id="b13"><mixed-citation publication-type="journal"><name><surname>Tanaka</surname><given-names>S.</given-names></name><article-title>Wing polymorphism, egg production and adult longevity in Pteronemobius taprobanensis Walker (Orthoptera, Gryllidae)</article-title>. <source>Kontyu</source><volume>44</volume>, <fpage>327</fpage>–<lpage>333</lpage> (<year>1976</year>).</mixed-citation></ref><ref id="b14"><mixed-citation publication-type="journal"><name><surname>Langellotto</surname><given-names>G. A.</given-names></name>, <name><surname>Denno</surname><given-names>R. F.</given-names></name> & <name><surname>Ott</surname><given-names>J. R.</given-names></name><article-title>A trade-off between flight capability and reproduction in males of a wing-dimorphic insect</article-title>. <source>Ecology</source><volume>81</volume>, <fpage>865</fpage>–<lpage>875</lpage> (<year>2000</year>).</mixed-citation></ref><ref id="b15"><mixed-citation publication-type="journal"><name><surname>Braendle</surname><given-names>C.</given-names></name>, <name><surname>Davis</surname><given-names>G. K.</given-names></name>, <name><surname>Brisson</surname><given-names>J. A.</given-names></name> & <name><surname>Stern</surname><given-names>D. L.</given-names></name><article-title>Wing dimorphism in aphids</article-title>. <source>Heredity</source><volume>97</volume>, <fpage>192</fpage>–<lpage>199</lpage> (<year>2006</year>).<pub-id pub-id-type="pmid">16823401</pub-id></mixed-citation></ref><ref id="b16"><mixed-citation publication-type="journal"><name><surname>Braendle</surname><given-names>C.</given-names></name>, <name><surname>Friebe</surname><given-names>I.</given-names></name>, <name><surname>Caillaud</surname><given-names>M. C.</given-names></name> & <name><surname>Stern</surname><given-names>D. L.</given-names></name><article-title>Genetic variation for an aphid wing polyphenism is genetically linked to a naturally occurring wing polymorphism</article-title>. <source>P. Roy. Soc. B-Biol. Sci.</source><volume>272</volume>, <fpage>657</fpage>–<lpage>664</lpage> (<year>2005</year>).</mixed-citation></ref><ref id="b17"><mixed-citation publication-type="journal"><name><surname>Mackay</surname><given-names>P. A.</given-names></name> & <name><surname>Wellington</surname><given-names>W.</given-names></name><article-title>A comparison of the reproductive patterns of apterous and alate virginoparous <italic>Acyrthosiphon pisum</italic> (Homoptera: Aphididae)</article-title>. <source>Can. Entomol.</source><volume>107</volume>, <fpage>1161</fpage>–<lpage>1166</lpage> (<year>1975</year>).</mixed-citation></ref><ref id="b18"><mixed-citation publication-type="journal"><name><surname>Tsumuki</surname><given-names>H.</given-names></name>, <name><surname>Nagatsuka</surname><given-names>H.</given-names></name>, <name><surname>Kawada</surname><given-names>K.</given-names></name> & <name><surname>Kanehisa</surname><given-names>K.</given-names></name><article-title>Comparison of nutrient reservation in apterous and alate pea aphids, <italic>Acyrthosiphon pisum</italic> (Harris). 1. Developmental time and sugar content</article-title>. <source>Appl. Entomol. Zool.</source><volume>25</volume>, <fpage>215</fpage>–<lpage>221</lpage> (<year>1990</year>).</mixed-citation></ref><ref id="b19"><mixed-citation publication-type="journal"><name><surname>Hazell</surname><given-names>S.</given-names></name>, <name><surname>Gwynn</surname><given-names>D.</given-names></name>, <name><surname>Ceccarelli</surname><given-names>S.</given-names></name> & <name><surname>Fellowes</surname><given-names>M.</given-names></name><article-title>Competition and dispersal in the pea aphid: clonal variation and correlations across traits</article-title>. <source>Ecol. Entomol.</source><volume>30</volume>, <fpage>293</fpage>–<lpage>298</lpage> (<year>2005</year>).</mixed-citation></ref><ref id="b20"><mixed-citation publication-type="journal"><name><surname>Fasulo</surname><given-names>T. R.</given-names></name> & <name><surname>Halbert</surname><given-names>S. E.</given-names></name><article-title>Aphid pests of Florida citrus</article-title>. <source>University of Florida - Institute of Food and Agricultural Sciences. Extension Publication</source><fpage>ENY811</fpage> (<year>1993</year>).</mixed-citation></ref><ref id="b21"><mixed-citation publication-type="journal"><name><surname>Halbert</surname><given-names>S. E.</given-names></name> & <name><surname>Brown</surname><given-names>L. G.</given-names></name><article-title><italic>Toxoptera citricida</italic> (Kirkaldy), Brown citrus aphid – Identification, biology, and management strategies</article-title>. <source>Florida Department of Agriculture and Consumer Services - Division of Plant Industry. Entomology Circular.</source><volume>374</volume> (<year>1998</year>).</mixed-citation></ref><ref id="b22"><mixed-citation publication-type="journal"><name><surname>Consortium</surname><given-names>I. A. G.</given-names></name><article-title>Genome sequence of the pea aphid <italic>Acyrthosiphon pisum</italic></article-title>. <source>PLoS Biol.</source><volume>8</volume>, <fpage>e1000313</fpage> (<year>2010</year>).<pub-id pub-id-type="pmid">20186266</pub-id></mixed-citation></ref><ref id="b23"><mixed-citation publication-type="journal"><name><surname>Nicholson</surname><given-names>S. J.</given-names></name><etal></etal>. <article-title>The genome of <italic>Diuraphis noxia</italic>, a global aphid pest of small grains</article-title>. <source>BMC Genomics</source><volume>16</volume>, <fpage>1</fpage> (<year>2015</year>).<pub-id pub-id-type="pmid">25553907</pub-id></mixed-citation></ref><ref id="b24"><mixed-citation publication-type="journal"><name><surname>Hunter</surname><given-names>W.</given-names></name><etal></etal>. <article-title>Aphid biology: expressed genes from alate <italic>Toxoptera citricida</italic> (Kirkaldy), brown citrus aphid (Hemiptera: Aphididae)</article-title>. <source>J. Insect Sci.</source><volume>3</volume>, <fpage>23</fpage> (<year>2003</year>).<pub-id pub-id-type="pmid">15841239</pub-id></mixed-citation></ref><ref id="b25"><mixed-citation publication-type="journal"><name><surname>Tang</surname><given-names>L. D.</given-names></name><etal></etal>. <article-title><italic>De novo</italic> sequencing-based transcriptome and digital gene expression analysis reveals insecticide resistance-relevant genes in <italic>Propylaea japonica</italic> (Thunberg) (Coleoptea: Coccinellidae)</article-title>. <source>PLoS ONE</source><volume>9</volume>, <fpage>e100946</fpage> (<year>2014</year>).<pub-id pub-id-type="pmid">24959827</pub-id></mixed-citation></ref><ref id="b26"><mixed-citation publication-type="journal"><name><surname>Zhu</surname><given-names>J. Y.</given-names></name>, <name><surname>Yang</surname><given-names>P.</given-names></name>, <name><surname>Zhang</surname><given-names>Z.</given-names></name>, <name><surname>Wu</surname><given-names>G. X.</given-names></name> & <name><surname>Yang</surname><given-names>B.</given-names></name><article-title>Transcriptomic immune response of <italic>Tenebrio molitor</italic> pupae to parasitization by <italic>Scleroderma guani</italic></article-title>. <source>PLoS ONE</source><volume>8</volume>, <fpage>e54411</fpage> (<year>2013</year>).<pub-id pub-id-type="pmid">23342153</pub-id></mixed-citation></ref><ref id="b27"><mixed-citation publication-type="journal"><name><surname>Li</surname><given-names>Z.-Q.</given-names></name><etal></etal>. <article-title>First transcriptome and digital gene expression analysis in Neuroptera with an emphasis on chemoreception genes in <italic>Chrysopa pallens</italic> (Rambur)</article-title>. <source>PLoS ONE</source><volume>8</volume>, <fpage>e67151</fpage> (<year>2013</year>).<pub-id pub-id-type="pmid">23826220</pub-id></mixed-citation></ref><ref id="b28"><mixed-citation publication-type="journal"><name><surname>Kola</surname><given-names>V. S. R.</given-names></name><etal></etal>. <article-title>Silencing of <italic>CYP6</italic> and <italic>APN</italic> genes affects the growth and development of rice yellow stem borer, <italic>Scirpophaga incertulas</italic></article-title>. <source>Front. Physiol.</source><volume>7</volume>, <fpage>20</fpage> (<year>2016</year>).<pub-id pub-id-type="pmid">26903874</pub-id></mixed-citation></ref><ref id="b29"><mixed-citation publication-type="journal"><name><surname>Yang</surname><given-names>W. J.</given-names></name>, <name><surname>Xu</surname><given-names>K. K.</given-names></name>, <name><surname>Cong</surname><given-names>L.</given-names></name> & <name><surname>Wang</surname><given-names>J. J.</given-names></name><article-title>Identification, mRNA expression, and functional analysis of chitin synthase 1 gene and its two alternative splicing variants in oriental fruit fly</article-title>. <source>Bactrocera dorsalis. Int. J. Biol. Sci.</source><volume>9</volume>, <fpage>331</fpage>–<lpage>342</lpage> (<year>2013</year>).<pub-id pub-id-type="pmid">23569438</pub-id></mixed-citation></ref><ref id="b30"><mixed-citation publication-type="journal"><name><surname>El-Shesheny</surname><given-names>I.</given-names></name>, <name><surname>Hajeri</surname><given-names>S.</given-names></name>, <name><surname>El-Hawary</surname><given-names>I.</given-names></name>, <name><surname>Gowda</surname><given-names>S.</given-names></name> & <name><surname>Killiny</surname><given-names>N.</given-names></name><article-title>Silencing abnormal wing disc gene of the Asian citrus psyllid, <italic>Diaphorina citri</italic> disrupts adult wing development and increases nymph mortality</article-title>. <source>PLoS ONE</source><volume>8</volume>, <fpage>e65392</fpage> (<year>2013</year>).<pub-id pub-id-type="pmid">23734251</pub-id></mixed-citation></ref><ref id="b31"><mixed-citation publication-type="journal"><name><surname>Castañeda</surname><given-names>L. E.</given-names></name>, <name><surname>Figueroa</surname><given-names>C. C.</given-names></name>, <name><surname>Bacigalupe</surname><given-names>L. D.</given-names></name> & <name><surname>Nespolo</surname><given-names>R. F.</given-names></name><article-title>Effects of wing polyphenism, aphid genotype and host plant chemistry on energy metabolism of the grain aphid</article-title>. <source>Sitobion avenae. J. Insect Physiol.</source><volume>56</volume>, <fpage>1920</fpage>–<lpage>1924</lpage> (<year>2010</year>).<pub-id pub-id-type="pmid">20801126</pub-id></mixed-citation></ref><ref id="b32"><mixed-citation publication-type="journal"><name><surname>Kim</surname><given-names>H.</given-names></name> & <name><surname>Lee</surname><given-names>S.</given-names></name><article-title>A molecular phylogeny of the tribe Aphidini (Insecta: Hemiptera: Aphididae) based on the mitochondrial tRNA/COII, 12S/16S and the nuclear EF1α genes</article-title>. <source>Syst. Entomol.</source><volume>33</volume>, <fpage>711</fpage>–<lpage>721</lpage> (<year>2008</year>).</mixed-citation></ref><ref id="b33"><mixed-citation publication-type="journal"><name><surname>von Dohlen</surname><given-names>C. D.</given-names></name>, <name><surname>Rowe</surname><given-names>C. A.</given-names></name> & <name><surname>Heie</surname><given-names>O. E.</given-names></name><article-title>A test of morphological hypotheses for tribal and subtribal relationships of Aphidinae (Insecta: Hemiptera: Aphididae) using DNA sequences</article-title>. <source>Mol. Phylogenet. Evol.</source><volume>38</volume>, <fpage>316</fpage>–<lpage>329</lpage> (<year>2006</year>).<pub-id pub-id-type="pmid">16368250</pub-id></mixed-citation></ref><ref id="b34"><mixed-citation publication-type="journal"><name><surname>Rivers</surname><given-names>D. B.</given-names></name>, <name><surname>Hink</surname><given-names>W. F.</given-names></name> & <name><surname>Denlinger</surname><given-names>D. L.</given-names></name><article-title>Toxicity of the venom from <italic>Nasonia vitripennis</italic> (Hymenoptera: Pteromalidae) toward fly hosts, nontarget insects, different developmental stages, and cultured insect cells</article-title>. <source>Toxicon</source><volume>31</volume>, <fpage>755</fpage>–<lpage>765</lpage> (<year>1993</year>).<pub-id pub-id-type="pmid">8342173</pub-id></mixed-citation></ref><ref id="b35"><mixed-citation publication-type="journal"><name><surname>Harvey</surname><given-names>J. A.</given-names></name>, <name><surname>Bezemer</surname><given-names>T. M.</given-names></name>, <name><surname>Elzinga</surname><given-names>J. A.</given-names></name> & <name><surname>Strand</surname><given-names>M. R.</given-names></name><article-title>Development of the solitary endoparasitoid <italic>Microplitis demolitor</italic>: host quality does not increase with host age and size</article-title>. <source>Ecol. Entomol.</source><volume>29</volume>, <fpage>35</fpage>–<lpage>43</lpage> (<year>2004</year>).</mixed-citation></ref><ref id="b36"><mixed-citation publication-type="journal"><name><surname>Brogiolo</surname><given-names>W.</given-names></name><etal></etal>. <article-title>An evolutionarily conserved function of the <italic>Drosophila</italic> insulin receptor and insulin-like peptides in growth control</article-title>. <source>Curr. Biol.</source><volume>11</volume>, <fpage>213</fpage>–<lpage>221</lpage> (<year>2001</year>).<pub-id pub-id-type="pmid">11250149</pub-id></mixed-citation></ref><ref id="b37"><mixed-citation publication-type="journal"><name><surname>Teleman</surname><given-names>A.</given-names></name><article-title>Molecular mechanisms of metabolic regulation by insulin in <italic>Drosophila</italic></article-title>. <source>Biochem. J.</source><volume>425</volume>, <fpage>13</fpage>–<lpage>26</lpage> (<year>2010</year>).<pub-id pub-id-type="pmid">20001959</pub-id></mixed-citation></ref><ref id="b38"><mixed-citation publication-type="journal"><name><surname>Xu</surname><given-names>H. J.</given-names></name><etal></etal>. <article-title>Two insulin receptors determine alternative wing morphs in planthoppers</article-title>. <source>Nature</source><volume>519</volume>, <fpage>464</fpage>–<lpage>467</lpage> (<year>2015</year>).<pub-id pub-id-type="pmid">25799997</pub-id></mixed-citation></ref><ref id="b39"><mixed-citation publication-type="journal"><name><surname>Yang</surname><given-names>P.</given-names></name>, <name><surname>Chen</surname><given-names>X.-M.</given-names></name>, <name><surname>Liu</surname><given-names>W.-W.</given-names></name>, <name><surname>Feng</surname><given-names>Y.</given-names></name> & <name><surname>Sun</surname><given-names>T.</given-names></name><article-title>Transcriptome analysis of sexually dimorphic Chinese white wax scale insects reveals key differences in developmental programs and transcription factor expression</article-title>. <source>Sci. Rep.</source><volume>5</volume>, <fpage>8141</fpage> (<year>2015</year>).<pub-id pub-id-type="pmid">25634031</pub-id></mixed-citation></ref><ref id="b40"><mixed-citation publication-type="journal"><name><surname>Brisson</surname><given-names>J. A.</given-names></name>, <name><surname>Ishikawa</surname><given-names>A.</given-names></name> & <name><surname>Miura</surname><given-names>T.</given-names></name><article-title>Wing development genes of the pea aphid and differential gene expression between winged and unwinged morphs</article-title>. <source>Insect Mol. Biol.</source><volume>19</volume>, <fpage>63</fpage>–<lpage>73</lpage> (<year>2010</year>).<pub-id pub-id-type="pmid">20482640</pub-id></mixed-citation></ref><ref id="b41"><mixed-citation publication-type="journal"><name><surname>Wang</surname><given-names>X.</given-names></name><etal></etal>. <article-title>The locust genome provides insight into swarm formation and long-distance flight</article-title>. <source>Nat. Commun.</source><volume>5</volume>, <fpage>2957</fpage> (<year>2014</year>).<pub-id pub-id-type="pmid">24423660</pub-id></mixed-citation></ref><ref id="b42"><mixed-citation publication-type="journal"><name><surname>Wu</surname><given-names>H.-J.</given-names></name>, <name><surname>Zhao</surname><given-names>L.-Q.</given-names></name>, <name><surname>Zeng</surname><given-names>Y.</given-names></name> & <name><surname>Zhu</surname><given-names>D.-H.</given-names></name><article-title>Effect of food stress on flight muscle and reproduction development in a wing dimorphic cricket, <italic>Velarifictorus ornatus</italic></article-title>. <source>Acta Ecologica Sinica</source><volume>34</volume>, <fpage>2963</fpage>–<lpage>2970</lpage> (<year>2014</year>).</mixed-citation></ref><ref id="b43"><mixed-citation publication-type="journal"><name><surname>Brisson</surname><given-names>J. A.</given-names></name>, <name><surname>Davis</surname><given-names>G. K.</given-names></name> & <name><surname>Stern</surname><given-names>D. L.</given-names></name><article-title>Common genome‐wide patterns of transcript accumulation underlying the wing polyphenism and polymorphism in the pea aphid (<italic>Acyrthosiphon pisum</italic>)</article-title>. <source>Evol. Dev.</source><volume>9</volume>, <fpage>338</fpage>–<lpage>346</lpage> (<year>2007</year>).<pub-id pub-id-type="pmid">17651358</pub-id></mixed-citation></ref><ref id="b44"><mixed-citation publication-type="journal"><name><surname>Grabherr</surname><given-names>M. G.</given-names></name><etal></etal>. <article-title>Full-length transcriptome assembly from RNA-Seq data without a reference genome</article-title>. <source>Nat. Biotechnol.</source><volume>29</volume>, <fpage>644</fpage>–<lpage>652</lpage> (<year>2011</year>).<pub-id pub-id-type="pmid">21572440</pub-id></mixed-citation></ref><ref id="b45"><mixed-citation publication-type="journal"><name><surname>Deng</surname><given-names>Y.</given-names></name><etal></etal>. <article-title>Integrated nr Database in protein annotation system and its localization</article-title>. <source>Comput. Eng.</source><volume>32</volume>, <fpage>71</fpage>–<lpage>76</lpage> (<year>2006</year>).</mixed-citation></ref><ref id="b46"><mixed-citation publication-type="journal"><name><surname>Apweiler</surname><given-names>R.</given-names></name><etal></etal>. <article-title>UniProt: the universal protein knowledgebase</article-title>. <source>Nucleic Acids Res.</source><volume>32</volume>, <fpage>D115</fpage>–<lpage>D119</lpage> (<year>2004</year>).<pub-id pub-id-type="pmid">14681372</pub-id></mixed-citation></ref><ref id="b47"><mixed-citation publication-type="journal"><name><surname>Kanehisa</surname><given-names>M.</given-names></name>, <name><surname>Goto</surname><given-names>S.</given-names></name>, <name><surname>Kawashima</surname><given-names>S.</given-names></name>, <name><surname>Okuno</surname><given-names>Y.</given-names></name> & <name><surname>Hattori</surname><given-names>M.</given-names></name><article-title>The KEGG resource for deciphering the genome</article-title>. <source>Nucleic Acids Res.</source><volume>32</volume>, <fpage>D277</fpage>–<lpage>D280</lpage> (<year>2004</year>).<pub-id pub-id-type="pmid">14681412</pub-id></mixed-citation></ref><ref id="b48"><mixed-citation publication-type="journal"><name><surname>Ashburner</surname><given-names>M.</given-names></name><etal></etal>. <article-title>Gene Ontology: tool for the unification of biology</article-title>. <source>Nat. Genet.</source><volume>25</volume>, <fpage>25</fpage>–<lpage>29</lpage> (<year>2000</year>).<pub-id pub-id-type="pmid">10802651</pub-id></mixed-citation></ref><ref id="b49"><mixed-citation publication-type="journal"><name><surname>Tatusov</surname><given-names>R. L.</given-names></name>, <name><surname>Galperin</surname><given-names>M. Y.</given-names></name>, <name><surname>Natale</surname><given-names>D. A.</given-names></name> & <name><surname>Koonin</surname><given-names>E. V.</given-names></name><article-title>The COG database: a tool for genome-scale analysis of protein functions and evolution</article-title>. <source>Nucleic Acids Res.</source><volume>28</volume>, <fpage>33</fpage>–<lpage>36</lpage> (<year>2000</year>).<pub-id pub-id-type="pmid">10592175</pub-id></mixed-citation></ref><ref id="b50"><mixed-citation publication-type="journal"><name><surname>Li</surname><given-names>R.</given-names></name><etal></etal>. <article-title>SOAP2: an improved ultrafast tool for short read alignment</article-title>. <source>Bioinf.</source><volume>25</volume>, <fpage>1966</fpage>–<lpage>1967</lpage> (<year>2009</year>).</mixed-citation></ref><ref id="b51"><mixed-citation publication-type="journal"><name><surname>Mortazavi</surname><given-names>A.</given-names></name>, <name><surname>Williams</surname><given-names>B. A.</given-names></name>, <name><surname>McCue</surname><given-names>K.</given-names></name>, <name><surname>Schaeffer</surname><given-names>L.</given-names></name> & <name><surname>Wold</surname><given-names>B.</given-names></name><article-title>Mapping and quantifying mammalian transcriptomes by RNA-Seq</article-title>. <source>Nat. Methods</source><volume>5</volume>, <fpage>621</fpage>–<lpage>628</lpage> (<year>2008</year>).<pub-id pub-id-type="pmid">18516045</pub-id></mixed-citation></ref><ref id="b52"><mixed-citation publication-type="journal"><name><surname>Audic</surname><given-names>S.</given-names></name> & <name><surname>Claverie</surname><given-names>J.-M.</given-names></name><article-title>The significance of digital gene expression profiles</article-title>. <source>Genome Res.</source><volume>7</volume>, <fpage>986</fpage>–<lpage>995</lpage> (<year>1997</year>).<pub-id pub-id-type="pmid">9331369</pub-id></mixed-citation></ref><ref id="b53"><mixed-citation publication-type="journal"><name><surname>Shang</surname><given-names>F.</given-names></name><etal></etal>. <article-title>Reference gene validation for quantitative PCR under various biotic and abiotic stress conditions in <italic>Toxoptera citricida</italic> (Hemiptera, Aphididae)</article-title>. <source>J. Econ. Entomol.</source><volume>108</volume>, <fpage>2040</fpage>–<lpage>2047</lpage> (<year>2015</year>).<pub-id pub-id-type="pmid">26470351</pub-id></mixed-citation></ref><ref id="b54"><mixed-citation publication-type="journal"><name><surname>Livak</surname><given-names>K. J.</given-names></name> & <name><surname>Schmittgen</surname><given-names>T. D.</given-names></name><article-title>Analysis of relative gene expression data using real-time quantitative PCR and the 2<sup>−ΔΔ<italic>C</italic>T</sup> method</article-title>. <source>Methods</source><volume>25</volume>, <fpage>402</fpage>–<lpage>408</lpage> (<year>2001</year>).<pub-id pub-id-type="pmid">11846609</pub-id></mixed-citation></ref><ref id="b55"><mixed-citation publication-type="journal"><name><surname>Luan</surname><given-names>J. B.</given-names></name>, <name><surname>Ghanim</surname><given-names>M.</given-names></name>, <name><surname>Liu</surname><given-names>S. S.</given-names></name> & <name><surname>Czosnek</surname><given-names>H.</given-names></name><article-title>Silencing the ecdysone synthesis and signaling pathway genes disrupts nymphal development in the whitefly</article-title>. <source>Insect Biochem. Mol. Biol.</source><volume>43</volume>, <fpage>740</fpage>–<lpage>746</lpage> (<year>2013</year>).<pub-id pub-id-type="pmid">23748027</pub-id></mixed-citation></ref><ref id="b56"><mixed-citation publication-type="journal"><name><surname>Jia</surname><given-names>S.</given-names></name>, <name><surname>Wan</surname><given-names>P.-J.</given-names></name>, <name><surname>Zhou</surname><given-names>L.-T.</given-names></name>, <name><surname>Mu</surname><given-names>L.-L.</given-names></name> & <name><surname>Li</surname><given-names>G.-Q.</given-names></name><article-title>Molecular cloning and RNA interference-mediated functional characterization of a Halloween gene spook in the white-backed planthopper <italic>Sogatella furcifera</italic></article-title>. <source>BMC Mol. Biol.</source><volume>14</volume>, <fpage>19</fpage> (<year>2013</year>).<pub-id pub-id-type="pmid">24007644</pub-id></mixed-citation></ref></ref-list><fn-group><fn><p><bold>Author Contributions</bold> F.S., D.D.W. and J.J.W. conceived and designed the experiments. F.S., B.Y.D. and Y.X. performed the experiments. F.S., D.W. and H.B.J. analyzed the data. F.S., W.D. and J.J.W. wrote the paper, and all authors reviewed the manuscript.</p></fn></fn-group></back><floats-group><fig id="f1"><label>Figure 1</label><caption><title>Homology analysis of unigenes of <italic>Toxoptera citricida</italic>.</title><p>(<bold>a</bold>) E-value distribution of BLAST hits for each unigene with a cut-off of 1.0E-5; (<bold>b</bold>) Similarity distribution of the top BLAST hits for each sequence; (<bold>c</bold>) Species distribution. The species distribution is shown as the percentage of the total homologous sequences in the NCBI Nr protein database with an E-value < 10-5.</p></caption><graphic xlink:href="srep32099-f1"></graphic></fig><fig id="f2"><label>Figure 2</label><caption><title>Relative expression levels of top 10 highly expressed genes in alates in different developmental stages of <italic>Toxoptera citricida</italic>.</title><p>N1, first instar nymphs; N2, second instar nymphs; N3, third instar nymphs; N4, fourth instar nymphs. Bars with a slash line indicate apterous aphids and bars with scatter indicate alate aphids in each stage, apterous and alate aphids cannot be distinguished by examining morphology of first and second instar nymphs. Data are means ± SE of 4 biological replicates and the elongation factor-1 alpha (<italic>EF1α</italic>) was used as a reference gene. Lowercase letter (<bold>a–c</bold>) above each bar indicates a significant difference among different developmental stages and wing morphs using one-way ANOVA followed by Tukey’s honestly significant difference test (<italic>P</italic> < 0.05). The relative expression was calculated using the 2<sup>−△△<italic>Ct</italic></sup>method based on the value of the first instar nymph expression, which was ascribed an arbitrary value of 1.</p></caption><graphic xlink:href="srep32099-f2"></graphic></fig><fig id="f3"><label>Figure 3</label><caption><title>Relative expression levels of the top 10 highly expressed genes in alate adults in different body parts of <italic>Toxoptera citricida</italic>.</title><p>The relative expression was calculated using the 2<sup>−△△<italic>Ct</italic></sup>method based on the value of the head of fourth instar winged-nymphs expression, which was assigned an arbitrary value of 1. Uppercase letter A, B, and C above each bar indicates a significant difference among different developmental stages and wing morphs within same body parts using one-way ANOVA followed by Tukey’s honestly significant difference test (<italic>P</italic> < 0.05). Lowercase letter a, b, and c above each bar indicates a significant difference among different body parts within same developmental stages and wing morphs using one-way ANOVA followed by Tukey’s honestly significant difference test (<italic>P</italic> < 0.05). A significant difference of gene expression level in wing/wing bud between fourth instar winged-nymphs and alate adults is indicated with asterisks (*<italic>P</italic> < 0.05; **<italic>P</italic> < 0.01; ***<italic>P</italic> < 0.001, Student’s <italic>t</italic> test).</p></caption><graphic xlink:href="srep32099-f3"></graphic></fig><fig id="f4"><label>Figure 4</label><caption><title>dsRNA feeding mediated RNAi and its effects on wing development of <italic>Toxoptera citricida</italic>.</title><p>(<bold>a</bold>) Representative phenotypes of <italic>Toxoptera citricida</italic> after feeding with dsRNA of <italic>KAT</italic>, <italic>PEPCK</italic>, and <italic>GP</italic>. (<bold>b</bold>) The rate of under-developed wing aphid after dsRNA feeding at 72 h. (<bold>c</bold>) Relative expression level of <italic>KAT</italic>, <italic>PEPCK</italic>, and <italic>GP</italic> after RNAi. The relative expression was calculated using the 2<sup>−△△<italic>Ct</italic></sup>method based on the value of ds<italic>GFP</italic> (control) at 72 h that was ascribed an arbitrary value of 1. Data are means ± SE of 4 replicates and the elongation factor-1 alpha (<italic>EF1α</italic>) was used as a reference gene. Lowercase letters above each bar indicate significant differences among different samples using one-way ANOVA followed by Tukey’s honestly significant difference test (<italic>P</italic> < 0.05).</p></caption><graphic xlink:href="srep32099-f4"></graphic></fig><table-wrap position="float" id="t1"><label>Table 1</label><caption><title>Sequence summary of the Illumina sequencing from <italic>Toxoptera citricida</italic> transcriptome.</title></caption><table frame="hsides" rules="groups" border="1"><colgroup><col align="left"></col><col align="center"></col></colgroup><thead valign="bottom"><tr><th colspan="2" align="left" valign="top" charoff="50">Sequencing</th></tr></thead><tbody valign="top"><tr><td align="left" valign="top" charoff="50">Total number of reads</td><td align="center" valign="top" charoff="50">41,250,684</td></tr><tr><td align="left" valign="top" charoff="50">Total number of clean reads</td><td align="center" valign="top" charoff="50">37,182,535</td></tr><tr><td align="left" valign="top" charoff="50">Total clean nucleotides (bp)</td><td align="center" valign="top" charoff="50">7,510,872,070</td></tr><tr><td align="left" valign="top" charoff="50">Q<sub>30</sub> percentage (%)</td><td align="center" valign="top" charoff="50">96.17</td></tr><tr><td align="left" valign="top" charoff="50">GC percentage (%)</td><td align="center" valign="top" charoff="50">41.97</td></tr><tr><td align="left" valign="top" charoff="50">Number of contigs</td><td align="center" valign="top" charoff="50">1,199,882</td></tr><tr><td align="left" valign="top" charoff="50">Length of contigs (bp)</td><td align="center" valign="top" charoff="50">85,765,139</td></tr><tr><td align="left" valign="top" charoff="50">N50 length of contigs (bp)</td><td align="center" valign="top" charoff="50">101</td></tr><tr><td align="left" valign="top" charoff="50">Number of unigenes</td><td align="center" valign="top" charoff="50">44,199</td></tr><tr><td align="left" valign="top" charoff="50">Length of unigenes (bp)</td><td align="center" valign="top" charoff="50">33,087,880</td></tr><tr><td align="left" valign="top" charoff="50">Mean length of unigenes (bp)</td><td align="center" valign="top" charoff="50">749</td></tr><tr><td align="left" valign="top" charoff="50">Unigenes annotations against NR</td><td align="center" valign="top" charoff="50">25,682</td></tr><tr><td align="left" valign="top" charoff="50">Unigenes annotations against Swiss-Prot</td><td align="center" valign="top" charoff="50">18,070</td></tr><tr><td align="left" valign="top" charoff="50">Unigenes annotations against KEGG</td><td align="center" valign="top" charoff="50">10,522</td></tr><tr><td align="left" valign="top" charoff="50">Unigenes annotations against GO</td><td align="center" valign="top" charoff="50">12,189</td></tr><tr><td align="left" valign="top" charoff="50">Unigenes annotations against COG</td><td align="center" valign="top" charoff="50">11,052</td></tr></tbody></table></table-wrap><table-wrap position="float" id="t2"><label>Table 2</label><caption><title>The top twenty enriched metabolic pathways among highly expressed transcripts (HETs) in alate adults.</title></caption><table frame="hsides" rules="groups" border="1"><colgroup><col align="left"></col><col align="left"></col><col align="center"></col><col align="center"></col><col align="center"></col></colgroup><thead valign="bottom"><tr><th align="left" valign="top" charoff="50">Pathway ID</th><th align="left" valign="top" charoff="50">Description</th><th align="center" valign="top" charoff="50">HETs</th><th align="center" valign="top" charoff="50">P-value</th><th align="center" valign="top" charoff="50">Q-value</th></tr></thead><tbody valign="top"><tr><td align="left" valign="top" charoff="50">ko01100</td><td align="left" valign="top" charoff="50">Metabolic pathways</td><td align="center" valign="top" charoff="50">33</td><td align="center" valign="top" charoff="50">0.01704447</td><td align="center" valign="top" charoff="50">8.59E-02</td></tr><tr><td align="left" valign="top" charoff="50">ko00071</td><td align="left" valign="top" charoff="50">Fatty acid metabolism</td><td align="center" valign="top" charoff="50">10</td><td align="center" valign="top" charoff="50">1.62E-08</td><td align="center" valign="top" charoff="50">1.96E-06</td></tr><tr><td align="left" valign="top" charoff="50">ko05164</td><td align="left" valign="top" charoff="50">Influenza A</td><td align="center" valign="top" charoff="50">9</td><td align="center" valign="top" charoff="50">0.000529084</td><td align="center" valign="top" charoff="50">1.12E-02</td></tr><tr><td align="left" valign="top" charoff="50">ko05130</td><td align="left" valign="top" charoff="50">Pathogenic Escherichia coli infection</td><td align="center" valign="top" charoff="50">8</td><td align="center" valign="top" charoff="50">2.69E-05</td><td align="center" valign="top" charoff="50">1.61E-03</td></tr><tr><td align="left" valign="top" charoff="50">ko05134</td><td align="left" valign="top" charoff="50">Legionellosis</td><td align="center" valign="top" charoff="50">7</td><td align="center" valign="top" charoff="50">3.98E-05</td><td align="center" valign="top" charoff="50">1.61E-03</td></tr><tr><td align="left" valign="top" charoff="50">ko00280</td><td align="left" valign="top" charoff="50">Valine, leucine and isoleucine degradation</td><td align="center" valign="top" charoff="50">7</td><td align="center" valign="top" charoff="50">0.000159331</td><td align="center" valign="top" charoff="50">4.82E-03</td></tr><tr><td align="left" valign="top" charoff="50">ko04145</td><td align="left" valign="top" charoff="50">Phagosome</td><td align="center" valign="top" charoff="50">7</td><td align="center" valign="top" charoff="50">0.006746751</td><td align="center" valign="top" charoff="50">4.30E-02</td></tr><tr><td align="left" valign="top" charoff="50">ko05162</td><td align="left" valign="top" charoff="50">Measles</td><td align="center" valign="top" charoff="50">6</td><td align="center" valign="top" charoff="50">0.000596582</td><td align="center" valign="top" charoff="50">1.12E-02</td></tr><tr><td align="left" valign="top" charoff="50">ko00020</td><td align="left" valign="top" charoff="50">Citrate cycle (TCA cycle)</td><td align="center" valign="top" charoff="50">6</td><td align="center" valign="top" charoff="50">0.000738194</td><td align="center" valign="top" charoff="50">1.12E-02</td></tr><tr><td align="left" valign="top" charoff="50">ko00564</td><td align="left" valign="top" charoff="50">Glycerophospholipid metabolism</td><td align="center" valign="top" charoff="50">6</td><td align="center" valign="top" charoff="50">0.001874586</td><td align="center" valign="top" charoff="50">2.27E-02</td></tr><tr><td align="left" valign="top" charoff="50">ko04144</td><td align="left" valign="top" charoff="50">Endocytosis</td><td align="center" valign="top" charoff="50">6</td><td align="center" valign="top" charoff="50">0.09049112</td><td align="center" valign="top" charoff="50">3.22E-01</td></tr><tr><td align="left" valign="top" charoff="50">ko04612</td><td align="left" valign="top" charoff="50">Antigen processing and presentation</td><td align="center" valign="top" charoff="50">5</td><td align="center" valign="top" charoff="50">0.002141193</td><td align="center" valign="top" charoff="50">2.36E-02</td></tr><tr><td align="left" valign="top" charoff="50">ko05020</td><td align="left" valign="top" charoff="50">Prion diseases</td><td align="center" valign="top" charoff="50">5</td><td align="center" valign="top" charoff="50">0.003264736</td><td align="center" valign="top" charoff="50">3.29E-02</td></tr><tr><td align="left" valign="top" charoff="50">ko04540</td><td align="left" valign="top" charoff="50">Gap junction</td><td align="center" valign="top" charoff="50">5</td><td align="center" valign="top" charoff="50">0.004379756</td><td align="center" valign="top" charoff="50">3.79E-02</td></tr><tr><td align="left" valign="top" charoff="50">ko05145</td><td align="left" valign="top" charoff="50">Toxoplasmosis</td><td align="center" valign="top" charoff="50">5</td><td align="center" valign="top" charoff="50">0.005449289</td><td align="center" valign="top" charoff="50">4.22E-02</td></tr><tr><td align="left" valign="top" charoff="50">ko03320</td><td align="left" valign="top" charoff="50">PPAR signaling pathway</td><td align="center" valign="top" charoff="50">5</td><td align="center" valign="top" charoff="50">0.006048017</td><td align="center" valign="top" charoff="50">4.22E-02</td></tr><tr><td align="left" valign="top" charoff="50">ko04530</td><td align="left" valign="top" charoff="50">Tight junction</td><td align="center" valign="top" charoff="50">5</td><td align="center" valign="top" charoff="50">0.1223737</td><td align="center" valign="top" charoff="50">3.54E-01</td></tr><tr><td align="left" valign="top" charoff="50">ko05169</td><td align="left" valign="top" charoff="50">Epstein-Barr virus infection</td><td align="center" valign="top" charoff="50">5</td><td align="center" valign="top" charoff="50">0.2186322</td><td align="center" valign="top" charoff="50">4.99E-01</td></tr><tr><td align="left" valign="top" charoff="50">ko04141</td><td align="left" valign="top" charoff="50">Protein processing in endoplasmic reticulum</td><td align="center" valign="top" charoff="50">5</td><td align="center" valign="top" charoff="50">0.2950153</td><td align="center" valign="top" charoff="50">5.58E-01</td></tr><tr><td align="left" valign="top" charoff="50">ko04910</td><td align="left" valign="top" charoff="50">Insulin signaling pathway</td><td align="center" valign="top" charoff="50">4</td><td align="center" valign="top" charoff="50">0.2388871</td><td align="center" valign="top" charoff="50">5.26E-01</td></tr></tbody></table></table-wrap><table-wrap position="float" id="t3"><label>Table 3</label><caption><title>The top ten highly expressed transcripts in alate adults of <italic>T. citricida</italic>.</title></caption><table frame="hsides" rules="groups" border="1"><colgroup><col align="left"></col><col align="left"></col><col align="center"></col><col align="center"></col></colgroup><thead valign="bottom"><tr><th align="left" valign="top" charoff="50">Gene ID</th><th align="left" valign="top" charoff="50">Annotation</th><th align="center" valign="top" charoff="50">Abbreviation</th><th align="center" valign="top" charoff="50">Change fold</th></tr></thead><tbody valign="top"><tr><td align="left" valign="top" charoff="50">c12058.graph_c0</td><td align="left" valign="top" charoff="50">Sodium- and chloride-dependent creatine transporter 1-like isoform 1</td><td align="center" valign="top" charoff="50"><italic>SCDCT</italic></td><td align="center" valign="top" charoff="50">8.43</td></tr><tr><td align="left" valign="top" charoff="50">c11063.graph_c0</td><td align="left" valign="top" charoff="50">3-ketoacyl-CoA thiolase, mitochondrial-like</td><td align="center" valign="top" charoff="50"><italic>KAT</italic></td><td align="center" valign="top" charoff="50">8.07</td></tr><tr><td align="left" valign="top" charoff="50">c9682.graph_c0</td><td align="left" valign="top" charoff="50">ATP-binding cassette sub-family G member 4-like isoform 1</td><td align="center" valign="top" charoff="50"><italic>ABCG4</italic></td><td align="center" valign="top" charoff="50">7.74</td></tr><tr><td align="left" valign="top" charoff="50">c14944.graph_c0</td><td align="left" valign="top" charoff="50">Trifunctional enzyme subunit beta, mitochondrial-like</td><td align="center" valign="top" charoff="50"><italic>β-TE</italic></td><td align="center" valign="top" charoff="50">5.76</td></tr><tr><td align="left" valign="top" charoff="50">c15715.graph_c0</td><td align="left" valign="top" charoff="50">Short-chain specific acyl-CoA dehydrogenase, mitochondrial-like</td><td align="center" valign="top" charoff="50"><italic>SCAD</italic></td><td align="center" valign="top" charoff="50">5.54</td></tr><tr><td align="left" valign="top" charoff="50">c15078.graph_c0</td><td align="left" valign="top" charoff="50">Trifunctional enzyme subunit alpha, mitochondrial-like</td><td align="center" valign="top" charoff="50"><italic>α-TE</italic></td><td align="center" valign="top" charoff="50">5.31</td></tr><tr><td align="left" valign="top" charoff="50">c15087.graph_c0</td><td align="left" valign="top" charoff="50">Phosphoenolpyruvate carboxykinase [GTP]-like</td><td align="center" valign="top" charoff="50"><italic>PEPCK</italic></td><td align="center" valign="top" charoff="50">5.17</td></tr><tr><td align="left" valign="top" charoff="50">c15313.graph_c0</td><td align="left" valign="top" charoff="50">Hydroxysteroid dehydrogenase-like protein 2-like</td><td align="center" valign="top" charoff="50"><italic>HSD</italic></td><td align="center" valign="top" charoff="50">4.72</td></tr><tr><td align="left" valign="top" charoff="50">c13625.graph_c0</td><td align="left" valign="top" charoff="50">E3 ubiquitin-protein ligase TRIM23-like</td><td align="center" valign="top" charoff="50"><italic>E3-TRIM23</italic></td><td align="center" valign="top" charoff="50">4.31</td></tr><tr><td align="left" valign="top" charoff="50">c11587.graph_c1</td><td align="left" valign="top" charoff="50">Glycogen phosphorylase-like isoform 2</td><td align="center" valign="top" charoff="50"><italic>GP</italic></td><td align="center" valign="top" charoff="50">4.29</td></tr></tbody></table></table-wrap></floats-group></pmc></record>Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Wicri/Bois/explor/OrangerV1/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 000476 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000476 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien
|wiki= Wicri/Bois
|area= OrangerV1
|flux= Pmc
|étape= Corpus
|type= RBID
|clé= PMC:5006003
|texte= Differential expression of genes in the alate and apterous morphs of the brown citrus aphid, Toxoptera citricida
}}
Pour générer des pages wiki
HfdIndexSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/RBID.i -Sk "pubmed:27577531" \
| HfdSelect -Kh $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd \
| NlmPubMed2Wicri -a OrangerV1
| This area was generated with Dilib version V0.6.25. |  |