Transcriptome analysis of Gossypium hirsutum flower buds infested by cotton boll weevil (Anthonomus grandis) larvae

• Published in: BMC genomics Vol. 15; no. 1; p. 854
• Main Authors: Artico, Sinara, Ribeiro-Alves, Marcelo, Oliveira-Neto, Osmundo Brilhante, de Macedo, Leonardo Lima Pepino, Silveira, Sylvia, Grossi-de-Sa, Maria Fátima, Martinelli, Adriana Pinheiro, Alves-Ferreira, Marcio
• Format: Journal Article
• Language: English
• Published: England BioMed Central Ltd 04.10.2014; BioMed Central
• Subjects: Analysis; Animals; Anthonomus grandis; Cotton; Crops; DNA binding proteins; Enzymes; Flowers & plants; Flowers - genetics; Gene expression; Gene Expression Profiling - methods; Gene Expression Regulation, Plant; Genes; Genetic aspects; Genetic transcription; Genetics; Genomics; Gossypium - anatomy & histology; Gossypium - genetics; Gossypium hirsutum; Herbivory; Insects; Laboratories; Larva - pathogenicity; Mathematical models; Molecular biology; Phylogenetics; Phylogeny; Plant Proteins - genetics; Sequence Analysis, RNA - methods; Sequence Homology, Nucleic Acid; Signal transduction; Weevils - growth & development; Weevils - physiology
• ISSN: 1471-2164; 1471-2164
• DOI: 10.1186/1471-2164-15-854

Cotton is a major fibre crop grown worldwide that suffers extensive damage from chewing insects, including the cotton boll weevil larvae (Anthonomus grandis). Transcriptome analysis was performed to understand the molecular interactions between Gossypium hirsutum L. and cotton boll weevil larvae. The Illumina HiSeq 2000 platform was used to sequence the transcriptome of cotton flower buds infested with boll weevil larvae. The analysis generated a total of 327,489,418 sequence reads that were aligned to the G. hirsutum reference transcriptome. The total number of expressed genes was over 21,697 per sample with an average length of 1,063 bp. The DEGseq analysis identified 443 differentially expressed genes (DEG) in cotton flower buds infected with boll weevil larvae. Among them, 402 (90.7%) were up-regulated, 41 (9.3%) were down-regulated and 432 (97.5%) were identified as orthologues of A. thaliana genes using Blastx. Mapman analysis of DEG indicated that many genes were involved in the biotic stress response spanning a range of functions, from a gene encoding a receptor-like kinase to genes involved in triggering defensive responses such as MAPK, transcription factors (WRKY and ERF) and signalling by ethylene (ET) and jasmonic acid (JA) hormones. Furthermore, the spatial expression pattern of 32 of the genes responsive to boll weevil larvae feeding was determined by "in situ" qPCR analysis from RNA isolated from two flower structures, the stamen and the carpel, by laser microdissection (LMD). A large number of cotton transcripts were significantly altered upon infestation by larvae. Among the changes in gene expression, we highlighted the transcription of receptors/sensors that recognise chitin or insect oral secretions; the altered regulation of transcripts encoding enzymes related to kinase cascades, transcription factors, Ca2+ influxes, and reactive oxygen species; and the modulation of transcripts encoding enzymes from phytohormone signalling pathways. These data will aid in the selection of target genes to genetically engineer cotton to control the cotton boll weevil.