Next Article in Journal
Andean Plants Essential Oils: A Scented Alternative to Synthetic Insecticides for the Control of Blowflies
Next Article in Special Issue
Genus-Wide Characterization of Nuclear Mitochondrial DNAs in Bumblebee (Hymenoptera: Apidae) Genomes
Previous Article in Journal / Special Issue
De Novo Genome Assembly of Chinese Plateau Honeybee Unravels Intraspecies Genetic Diversity in the Eastern Honeybee, Apis cerana
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Identification of Boll Weevil, Anthonomus grandis grandis (Coleoptera: Curculionidae), Genes Involved in Pheromone Production and Pheromone Biosynthesis

Insect and Cotton Disease Research Unit, USDA Agricultural Research Service, College Station, TX 77845, USA
*
Author to whom correspondence should be addressed.
Insects 2021, 12(10), 893; https://doi.org/10.3390/insects12100893
Submission received: 1 September 2021 / Revised: 20 September 2021 / Accepted: 24 September 2021 / Published: 1 October 2021
(This article belongs to the Special Issue Insect Genomics)

Abstract

:

Simple Summary

The boll weevil is a destructive pest of commercial cotton throughout the Americas. An eradication program in the United States has removed the boll weevil from most of its range. However, weevil populations in South Texas remain a threat to eradicated areas. Pheromone traps are used to monitor boll weevil activity, and when a weevil is captured, eradication programs rely on malathion for control. However, the effectiveness of pheromone traps in detecting incipient boll weevil populations is reduced during certain times of the year. Additionally, human safety and environmental concerns, as well as the potential development of malathion-resistant populations, have prompted program managers to seek alternative control methods. We sequenced and compared pheromone-producing and non-pheromone-producing weevils to identify genes involved in pheromone production, which, in turn, could be an environmentally friendly way to target gene-level pest control that is specific to the boll weevil. Our results revealed genes involved in pheromone production, as well as insect development and immunity, which may be targeted for boll weevil suppression.

Abstract

Eradication programs for the boll weevil, Anthonomus grandis grandis Boheman (Coleoptera: Curculionidae), rely almost exclusively on pheromone traps to indicate the need for insecticide applications. However, the effectiveness of traps in detecting weevil populations is reduced during certain times of the year, particularly when cotton is actively fruiting. Consequently, this could result in fields becoming heavily infested with weevils. It is widely speculated that the lack of weevil captures in traps during this period is largely due to the overwhelming amount of pheromone released by weevils in the field, which outcompete the pheromone released from traps. Thus, this work sought to identify genes involved in pheromone production so that new control methods that target these genes can be explored. We conducted an RNA-seq experiment that revealed 2479 differentially expressed genes between pheromone-producing and non-pheromone-producing boll weevils. Of those genes, 1234 were up-regulated, and 1515 were down-regulated, and most had gene annotations associated with pheromone production, development, or immunity. This work advances our understanding of boll weevil pheromone production and brings us one step closer to developing gene-level control strategies for this cotton pest.

1. Introduction

The boll weevil, Anthonomus grandis grandis Boheman (Coleoptera: Curculionidae), is a pest of commercial cotton (Gossypium hirsutum) throughout the Americas [1]. Efforts to eradicate the boll weevil in the United States were initiated in the late 1970s, and since that time, the boll weevil has been eliminated from most of its U.S. range [2,3,4]. However, boll weevil populations remaining in the Lower Rio Grande Valley (LRGV) production area of Texas continue to pose a threat to other cotton-growing areas in the state and neighboring states (https://www.txbollweevil.org.html accessed on 19 September 2021). In fact, recent infestations in eradicated areas in Texas have resulted in multi-year infestations, which have cost millions of dollars to remediate (P. Burson, personal communication).
Eradication programs rely exclusively on pheromone traps to detect and monitor boll weevil populations [5]. Pheromone traps are typically placed next to cotton fields and monitored weekly during the active phase of eradication. Each trap is baited with a pheromone lure containing a nominal dose of 10 mg of grandlure (synthesized boll weevil pheromone), which is composed of four pheromone components: two terpene alcohols (components I and II; (+)-cis-2-isopropenyl-1-methyl cyclobutaneethanol and cis-3,3-dimethyl-Δ 1,β-cyclohexaneethanol, respectively) and two terpene aldehydes (components III and IV; cis-3,3-dimethyl- Δ 1,α-cyclohexaneacetaldehyde and trans-3,3-dimethyl- Δ 1,α-cyclohexaneacetaldehyde) [6]. The pheromone is produced by male weevils and is attractive to both sexes and subsequently considered an aggregation rather than a sex pheromone [7]. During post-eradication, programs may place traps in close proximity to cotton fields or in a grid-like manner irrespective of the location of cotton fields. Traps are monitored every two or three weeks, and use lures dosed with 25 mg of grandlure and 30 mg of eugenol. Regardless, both active and post-eradication programs rely on traps to indicate the presence of weevils and the need for remedial actions.
It is generally recognized that weevil captures in traps vary seasonally and decline rapidly once cotton begins to produce fruiting structures [8,9,10]. In fact, some cotton fields may be heavily infested with boll weevils, yet adjacent traps may not capture any weevils [11]. Such scenarios could result in fields serving as nurseries for boll weevils, producing thousands to millions of boll weevils throughout the growing season if left untreated. Several reasons for the reduced performance of pheromone traps during the fruiting stage of cotton development have been proposed [8,9,10,12]. The most widely accepted reason is that weevils producing pheromone in the field outcompete the pheromone released from traps [13,14,15,16,17]. Nevertheless, once a weevil is detected in a trap, a malathion application is triggered for the respective field or growing area. Malathion treatments may continue weekly until weevils are no longer detected. Consequently, a field may receive as many as 25 applications in a single growing season (P. Burson, personal communication). Thus, new control methods are needed that are cleaner and safer for humans and the environment surrounding cotton fields and, ideally, specific to the boll weevil.
One way to avoid trap competition in the field is to develop a gene-based control strategy that suppresses the production of pheromone in the male weevil. This could increase the attractiveness of the pheromone lure in traps allowing for more reliable and early detection of weevils; thus, leading to more timely insecticide applications. Suppression of pheromone-producing weevils could also reduce and/or delay early-season colonization of cotton fields by overwintered weevils. Such delays or reductions could potentially prevent fields from becoming heavily infested [18], which in turn could reduce malathion applications per season. Here we present a boll weevil transcriptome created using next-generation RNA sequencing with a BUSCO score of 99.27%. We used the transcriptome to conduct an RNA-seq experiment to identify genes involved in pheromone production and pheromone biosynthesis of male boll weevils.

2. Materials and Methods

2.1. Boll Weevil Collections and Experimental Design

Boll weevils were initially collected as larvae in infested squares from a field northwest of Edinburg, TX (26.5°N, 98.2°W). The major life stages (early instar larvae, late instar larvae, pupae, adult male, and adult female) were placed in RNAlater (Invitrogen, Carlsbad, CA, USA), and stored at −20 °C until total RNA was extracted for the transcriptome assembly. A group of pupae from the same field was kept on vermiculite in a 29.4 ± 1 °C degree incubator to obtain adults for the pheromone experiment. Newly eclosed adults (<24 h old) were sexed, and only the males were kept. Forty males were divided into two groups. One group was held without food for 24 h. The other group was fed squares (6–9 mm in diameter; bracts intact) daily for 6–8 days and provided water on a section of cotton dental wick saturated with deionized water to promote the production of all four pheromone components [19,20]. Weevils in both groups were tested for pheromone production using gas chromatography-mass spectrometry and procedures described by Suh and Spurgeon (2016) [21]. Ten weevils in each group (i.e., pheromone-producing or non-pheromone-producing) were placed in 2-mL snap cap tubes with RNAlater and stored at −20 °C until nucleic acid was extracted.

2.2. RNA Extraction and Sequencing

RNA was extracted from weevils in each life stage using the RNeasy extraction kit (Qiagen, Hilden, Germany). Samples were checked for quantity and quality using a Tapestation 4200 (Agilent, Santa Clara, CA, USA). The total RNA from 12 samples at different life stages (two second instar larvae; two third instar larvae; two pupae; two adult females; and two adult males) and 20 weevils from the pheromone experiment (10 pheromone-producing males and 10 non-pheromone-producing males) were submitted to Texas A&M AgriLife Genomics and Bioinformatics Service (TxGen) for purification, library preparation, and sequencing on an Illumina NovaSeq (Illumina, San Diego, CA, USA), paired-end with a read length of 2 × 150 base pairs.

2.3. Transcriptome Assembly and Annotation

A total of approximately 30 million reads per sample were downloaded from TxGen and trimmed of adaptors and low-quality sequences using Trimmomatic 0.38 [22]. Approximately 97% of all life stage reads passed quality control measures resulting in a total of 352,189,343 sequences. Raw sequences are available at NCBI SRA (PRJNA734329). The trimmed reads were used to assemble a de novo transcriptome using the TRINITY 2.12.0 assembler [23]. The default parameters were used, including the options for pair-wise reads, non-strand specific assembly, minimum codon length of 200 bp, and the construction of super contigs. CD-HIT was used to cluster the sequences and reduce the number of isomers [24]. Completeness of the transcriptome was assessed using Benchmarking Universal Single-Copy Orthologs (BUSCO) [25,26,27]; selected from OrthoDB using the insecta database [28]. The transcriptome was annotated using CloudBLAST and InterProScan via OmicsBox (Valencia, Spain). Parameters were limited to the insecta database and E-value ≤ e × 10−5.

2.4. RNA-Seq Analysis

A total of approximately 30 million reads per sample (586,574,829 total reads) for the pheromone and non-pheromone samples were used in the RNA-seq analysis, which was accomplished using edgeR as part of the OmicsBox software [29]. Because no reference genome was available for pair-wise differential expression analysis, the reads were mapped to the de novo transcriptome using RNA-seq by Expectation-Maximization (RSEM). The workflow also included quantifying the mapped reads at each locus and statistically assessing fold change (False Discovery Rate (FDR) < 0.05) differences between pheromone-producing and non-pheromone-producing weevils [30]. Individual genes with relevant Gene Ontology (GO) terms were further examined manually. The complete list of differentially expressed genes (DEGs) can be found in Supplemental Table S1.

3. Results and Discussion

3.1. Transcriptome

The final compiled transcriptome had 195,128 total transcripts with an average length of 922.89 bp, Table 1. We assessed the completeness of the transcriptome by comparing our assembled transcripts to conserved orthologs found in other insects using BUSCO [25,26,27]. The transcriptome had a BUSCO score of 99.27% with only 0.59% and 0.15% fragmented or missing, respectively, Table 1.
The BLAST annotation revealed most sequences aligned to two weevil species, Dendroctonus ponderosae Hopkins and Sitophilus oryzae Linnaeus, followed by the tenebrionid, Tribolium castaneum Herbst. Nearly all sequences in the dataset aligned to other coleopterans, Figure 1.
Gene Ontology (GO) was performed for all three major categories: molecular function (MF), biological process (BP), and cellular component (CC). Standard gene groups were represented in all three categories, Figure 2. The largest categories for MF were protein binding (24%) and nucleic acid binding (13%), Figure 2a. Top ontologies for BP were cellular macromolecule biosynthetic process (17%) and organonitrogen compound biosynthetic process (10%), Figure 2b. CC GO terms were the highest for cytoplasm (30%) and an integral component of membrane (17%), Figure 2c.

3.2. RNA-Seq Analysis

RNA-seq analysis identified a total of 2479 DEGs between the males that produced pheromone compared to those that did not produce pheromone. Of those 2479 genes, 1234 genes were up-regulated, and 1515 genes were down-regulated, Figure 3 and Supplemental Table S1.
Evidence from other coleopterans shows pheromone is made by utilizing the first section of the juvenile hormone pathway, called the mevalonate pathway, which uses byproducts to produce the monoterpenoid pheromone component [31]. The remaining steps in the juvenile hormone pathway produce juvenile hormones, which are critical in the regulation of molting, development, and induction of pheromone production. Thus, pheromone production and insect development are correlated phenotypes relying, at least partially, on the same biochemical pathway. Additionally, we established the pheromone-producing and non-pheromone-producing phenotypes by (1) withholding food from newly eclosed weevils and (2) allowing weevils that were given food 6–8 days to stimulate the production of pheromone [20]. Thus, differences in the age and nutrient availability between treatment groups also likely contributed to the differences observed between pheromone-producing and non-pheromone-producing boll weevils. With that in mind, we first discuss genes that have a direct function in pheromone production and pheromone biosynthesis that were significantly DE between treatment groups, Table 2 and Figure 4. We then discuss other interesting genes found within the top 50 significantly up- and down-regulated DEG genes (based on FDR value) and make predictions on their function based on GO term descriptions and the literature, Table 2 and Table 3.
The mevalonate pathway in insects is a metabolic pathway initiated by reductive polymerization of acetyl-CoA and leads to a diversity of isoprenoid compounds [32]. Information on pheromone production and the mevalonate pathway in insects has been studied mostly in pine beetles and Drosophila; however, more studies working to elucidate the steps in the pathway are becoming available [31,33]. In this study, we detected eight genes directly involved in the mevalonate pathway that were up-regulated in pheromone-producing males, Table 2 and Figure 4. For example, hydroxymethylglutaryl-CoA synthase 1 (HMG-S), 3-hydroxy-3-methylglutaryl-coenzyme A reductase (HMG-R), and phosphomevalonate kinase were all highly up-regulated and found in the first steps of the mevalonate pathway converting acetyl-CoA to mevalonate 5-PP, the final step before isopentenyl synthesis, Figure 4. Additionally, our study identified ATP-citrate synthase, which converts citrate to acetyl-CoA, and is needed to initiate the mevalonate pathway, Table 2 and Figure 4.
Our study found up-regulation in genes in the isopentenyl, dolicol, and farnesyl portions of the pathway as well, Table 2 and Figure 4. We also found a gene involved in the processes that transfers sugar to glycoproteins during the dolichol pathway, dolichyl pyrophosphate. Similar molecules have also been shown to be part of the production of isoprenoid compounds in blowfly larvae [34]. Several up-regulated genes had links to the accumulation of lipids or lipid metabolism, fatty acid synthase, and glycerol-3-phosphate phosphatase, and most likely function to increase the size of the fat body and hormone production. The increase in a particular type of fat body is a feature that correlates with reproductive maturity in the boll weevil [20]. Furthermore, Tittiger and Blomquist (2016) reported that pine bark beetle pheromone also consists of components derived from fatty acids, which could also be the case in the boll weevil [35]. There were seven juvenile hormones (JH) up-regulated in pheromone-producing weevils. The final steps of the mevalonate pathway result in the production of JHs. While JHs have many functions in insects, they have been shown to regulate isoprenoid pheromone production in the midgut of I. pini males and are speculated to regulate the mevalonate pathway [36].
We also identified two genes that may function to stimulate hormone production in pheromone-producing male boll weevils, shuttle craft-like and regucalcin-like, Table 2. The gene shuttle craft has mostly been studied in its role in the embryonic central nervous system development in D. melanogaster [37], for which it is required. It also has been shown to have an expression in the nervous system and reproductive system in adult male and female flies. It has been speculated that shuttle craft regulates genes that are either involved in axon growth or the secretion of molecules needed for chemoattraction [38]. Interestingly, it was up-regulated in pheromone-producing male boll weevils and may function in laying down the neurons needed to produce pheromones. Likewise, the protein regucalcin-like has been studied very little in insects but has been characterized in mammals as a calcium-binding protein that regulates many cell functions needed for hormonal stimulation [39]. Thus, both genes seem to be involved in hormone stimulation that may be needed to activate the pheromone production pathway in the boll weevil.
Aside from genes directly involved in pheromone production, we also found genes up-regulated in associated processes, such as digestion and immunity. The most highly up-regulated gene in the dataset was a serine protease, transmembrane protease serine 9-like, Table 3. In fact, there were six serine protease genes up-regulated in the 50 most highly up-regulated gene list, Supplemental Table S1. Serine proteases serve many cellular functions, including development, apoptosis, and digestion [40,41]. The most highly expressed serine peptidase in this study had a high sequence identity to a gut serine peptidase in T. castaneum that was highly induced in response to feeding [42]. Thus, it is likely that transmembrane protease serine 9-like is also involved in digestion in the boll weevil. We found six other DE serine peptidase sequences, outside of the top 50 gene list, that had A. grandis as a top BLAST hit and were previously identified as midgut serine proteases [43]. In that study, two of these genes, Agser2p and Agser5p, were induced in response to feeding. The up-regulation of digestive peptidase genes in our dataset is fitting because pheromone-producing males were fed fresh squares daily whereas, the non-pheromone-producing males were not given food. The remaining up-regulated serine peptidase genes could be involved in digestion as well but may also have roles in development and apoptosis, both of which may be necessary for boll weevil sexual maturity and pheromone production.
Two immune genes were differentially expressed protein spaetzle 3 and DNA/RNA non-specific nuclease 1. In Drosophila, spaetzle is involved in embryonic development but also has a role in the innate immune response. Spaetzle proteins bind to Toll receptors and activate antimicrobial peptide gene expression [44]. The second gene, DNA/RNA non-specific nuclease 1, was shown to be in the midgut of Bombyx mori and interferes with the RNA interference response. Because of its propensity to degrade dsRNA, it has been suggested the enzyme might be involved in the innate immune response against invading nucleic acids, such as RNA viruses [45]. The functionality and increased expression in innate immune genes in pheromone-producing weevils is not well understood, but some studies suggest that male pheromone quality is a direct assessment of overall male health [46]. In fact, a study using Tenebrio molitor, showed males were more attractive to females if they were immunocompetent [47]. Thus, having a good immune system may affect the overall quality of pheromones.
Many of the genes that were significantly down-regulated in pheromone-producing weevils had functions consistent with the cuticle and some were directly involved in cuticle tanning and sclerotization, Table 3. This is the process where the cuticle darkens (melanin production), and the exoskeleton becomes hard [48,49]. This process is achieved through the tyrosine metabolic pathway and has been characterized mainly in the model beetle, T. castaneum. The gene laccase 2 was found to play a major role in cuticle tanning for all life stages except the egg stage [50]. Likewise, this gene was decreased nearly 19-fold in this study, Table 3. Dopamine also plays a critical role in the conjugation of β-alanine to produce Ν-β -alanyldopamine, which directly leads to pigmentation precursors [51]. We found two dopamine-related genes in our dataset, dopamine D2-like receptor and dopamine N-acetyltransferase, that may play a part in weevil cuticle tanning. Additionally, we found four cuticle genes, two chitin synthase genes, and many other cuticle genes outside of the top 50 DEGs, all down-regulated in pheromone-producing weevils, Table 3 and Supplemental Table S1. This pathway is most likely highlighted in our data because the pheromone-producing weevils were 6 to 8 days older compared to the non-pheromone-producing weevils. Thus, pheromone-producing weevils were likely finished with the tanning process, while non-pheromone-producing males, which were less than 24 h old, still had soft cuticles that were not fully sclerotized and had not become dark.
The gene transcription factor ken, was also significantly down-regulated in pheromone-producing boll weevils, Table 3. Interestingly, the ortholog in D. melanogaster, a transcription factor named ken and barbie, has a distinct role in the development of male and female external genitalia and is expressed in the highest amounts during embryo and pupal stages of development. In fact, mutants of this gene gave a range of deformities, including no external genitalia [52]. This is another example of an interesting gene to target with gene-based pest control strategies because reduced expression at key periods of insect development could lead to males with completely absent external reproductive organs, thus reducing successful mating in the field.

4. Conclusions

We provide a boll weevil de novo transcriptome with a BUSCO score of 99.27%. We used this assembly to perform an RNA-seq analysis comparing gene expression between pheromone-producing and non-pheromone-producing boll weevils. We ultimately identified significant DEGs that likely have key functions in pheromone production and biosynthesis. In subsequent experiments, successful knockdown of these DEGs should result in boll weevils with no or limited pheromone production, which could result in increased detection of weevils with traps. Because eradication programs in the U.S. rely almost exclusively on pheromone traps to indicate the need for insecticide applications, improved detection of boll weevils with traps could lead to more timely application of insecticides. Furthermore, disruption of boll weevil pheromone production may reduce early-season colonization of weevils that rely on pheromone to find cotton. Finally, the other DEGs identified in this study involved in immunity, cuticle formation, and digestion may also serve as targets for gene-level pest control. Thus, this RNA-seq analysis provides the foundation for the production of gene-level control strategies that may target pheromone production and other life functions in the boll weevil. Control strategies that incorporate these methods, such as RNA interference or CRISPR, would not only be specific to the boll weevil but would also be safe and environmentally friendly.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/insects12100893/s1, Table S1: Complete list of 192 differentially expressed genes between pheromone-producing and non-producing boll weevils.

Author Contributions

The author contributions for this manuscript are as follows: conceptualization, L.C.P. and C.P.-C.S.; methodology, L.C.P. and J.L.P.; validation, J.L.P. and C.P.-C.S.; formal analysis, L.C.P.; data curation, L.C.P.; writing—original draft preparation, L.C.P.; writing—review and editing, L.C.P., J.L.P. and C.P.-C.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Cotton Incorporated, Project # 20-216.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The submission of raw sequences to NCBI SRA has been started under project number PRJNA734329.

Acknowledgments

We thank Mercedes Perez for mobile lab space and Adriane Fifield for help editing the manuscript. This research was funded through Cotton Inc. (Project # 20-216). The mention of trade names or commercial products does not imply recommendation or endorsement by the USDA. USDA is an equal opportunity provider and employer.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Paula, D.P.; Claudino, D.; Timbó, R.V.; Miranda, J.E.; Bemquerer, M.; Ribeiro, A.; Sujii, E.; Fontes, E.; Pires, C. Reproductive dormancy in boll-weevil from populations of the Midwest of Brazil. J. Econ. Ѐntomol. 2013, 106, 86–96. [Google Scholar] [CrossRef]
  2. Smith, J.W. Boll weevil eradication: Area-wide pest management. Ann. Ѐntomol. Soc. Am. 1998, 91, 239–247. [Google Scholar] [CrossRef]
  3. Carter, F.L.; Nelson, T.C.; Jordan, A.G.; Smith, J.R. US cotton declares war on the boll weevil. In Boll Weevil Eradication in the United States through 1999; The Cotton Foundation: Memphis, TN, USA, 1999; pp. 25–54. [Google Scholar]
  4. El-Lissy, O.A.; Grefenstette, W.J. Progress of boll weevil eradication in the U.S., 2005. In Proceedings of the Beltwide Cotton Conferences, San Antonio, TX, USA, 3–6 January 2006; pp. 1266–1276. [Google Scholar]
  5. Texas Boll Weevil Eradication Foundation. Trapping. Available online: https://www.txbollweevil.org/trapping.html (accessed on 7 May 2021).
  6. Tumlinson, J.H.; Hardee, D.D.; Gueldner, R.C.; Thompson, A.C.; Hedin, P.A.; Minyard, J.P. Sex pheromones produced by male boll weevil: Isolation, identification, and synthesis. Science 1969, 166, 1010–1012. [Google Scholar] [CrossRef]
  7. Hardee, D.D.; McKibben, G.H.; Gueldner, R.C.; Mitchell, E.B.; Tumlinson, J.H.; Cross, W.H. Boll weevils in nature respond to grandlure, a synthetic pheromone. J. Econ. Ѐntomol. 1972, 65, 97–100. [Google Scholar]
  8. Rogers, C.; Oakes, S.; Rummel, D. Evaluation of infield pheromone traps for boll weevil suppression in the Texas Rolling Plains. Res. Monogr -Tex. Agric. Exp. Station 1976, 8, 45–52. [Google Scholar]
  9. Rummel, D.R.; Bottrell, D.G. Seasonally related decline in response of boll weevils to pheromone traps during mid-season. Environ. Ѐntomol. 1976, 5, 783–787. [Google Scholar] [CrossRef] [Green Version]
  10. Wolfenbarger, D.A.; Graham, H.M.; Parker, R.D.; Davis, J.W. Boll weevil: Seasonal patterns of response to traps baited with grandlure in the Lower Rio Grande Valley. Environ. Ѐntomol. 1976, 5, 403–408. [Google Scholar] [CrossRef]
  11. Suh, C.P.-C.; Westbrook, J.K. Failure of pheromone traps in detecting incipient populations of boll weevils (Coleoptera: Curculiondae): Investigation of two potential contributing factors. J. Ѐntomol. Sci. 2014, 49, 211–214. [Google Scholar] [CrossRef]
  12. Rummel, D.; Bottrell, D. Relationship of overwintered boll weevil response to pheromone traps and natural entry into cotton. Res. Monogr.-Tex. Agric. Exp. Stn. 1976, 8, 26–31. [Google Scholar]
  13. Hardee, D.D.; Cross, W.H.; Huddleston, P.M.; Davich, T.B. Survey and control of the boll weevil in West Texas with traps baited with males. J. Econ. Ѐntomol. 1970, 63, 1041–1048. [Google Scholar] [CrossRef]
  14. Boyd, F., Jr.; Brazzel, J.; Helms, W.; Moritz, R.; Edwards, R. Spring destruction of overwintered boll weevils in West Texas with wing traps. J. Econ. Entomol. 1973, 66, 507–510. [Google Scholar] [CrossRef] [Green Version]
  15. Scott, W.P.; Lloyd, E.P.; Bryson, J.O.; Davich, T.B. Trap plots for suppression of low-density overwintered populations of boll weevils. J. Econ. Ѐntomol. 1974, 67, 281–283. [Google Scholar] [CrossRef]
  16. Cherry, E.T. Monitoring boll weevil movement with pheromone traps. Tenn. Farm Home Sci. Prog. Rep. 1974, 90, 27–29. [Google Scholar]
  17. Segers, J.; Urban, T.; George, D.; Benedict, J.; Walmsley, M.; Pieters, E. Seasonal numbers, sex diapause states of boll weevils captured in pheromone traps in the Lower Gulf Coast of Texas. Southwest. Entomol. 1988, 12, 311–316. [Google Scholar]
  18. White, J.R.; Rummel, D.R. Emergence profile of overwintered boll weevils and entry into cotton. Environ. Ѐntomol. 1978, 7, 7–14. [Google Scholar] [CrossRef] [Green Version]
  19. Spurgeon, D.W. Age dependence of pheromone production by the boll weevil (Coleoptera: Curculionidae). Environ. Ѐntomol. 2003, 32, 31–38. [Google Scholar] [CrossRef]
  20. Spurgeon, D.W.; Suh, C.P.-C. Pheromone production by the boll weevil (Coleoptera: Curculionidae) fed cotton squares and bolls. J. Ѐntomol. Sci. 2009, 44, 209–221. [Google Scholar] [CrossRef] [Green Version]
  21. Suh, C.-C.; Spurgeon, D. Continued pheromone release by boll weevils (Coleoptera: Curculionidae) following host removal. J. Ѐntomol. Sci. 2016, 51, 332–335. [Google Scholar] [CrossRef]
  22. Bolger, A.M.; Lohse, M.; Usadel, B.J.B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [Green Version]
  23. Grabherr, M.G.; Haas, B.J.; Yassour, M.; Levin, J.; Thompson, D.A.; Amit, I.; Adiconis, X.; Fan, L.; Raychowdhury, R.; Zeng, Q.; et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 2011, 29, 644–652. [Google Scholar] [CrossRef] [Green Version]
  24. Li, W.; Godzik, A. Cd-hit: A fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics 2006, 22, 1658–1659. [Google Scholar] [CrossRef] [Green Version]
  25. Simão, F.A.; Waterhouse, R.; Ioannidis, P.; Kriventseva, E.V.; Zdobnov, E.M. BUSCO: Assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics 2015, 31, 3210–3212. [Google Scholar] [CrossRef] [Green Version]
  26. Waterhouse, R.; Seppey, M.; Simão, F.A.; Manni, M.; Ioannidis, P.; Klioutchnikov, G.; Kriventseva, E.V.; Zdobnov, E.M. BUSCO: Applications from quality assessments to gene prediction and phylogenomics. Mol. Biol. Evol. 2018, 35, 543–548. [Google Scholar] [CrossRef] [Green Version]
  27. Seppey, M.; Manni, M.; Zdobnov, E.M. BUSCO: Assessing genome assembly and annotation completeness. Bioinformatics 2019, 1962, 227–245. [Google Scholar]
  28. Kriventseva, E.V.; Kuznetsov, D.; Tegenfeldt, F.; Manni, M.; Dias, R.; Simão, F.A.; Zdobnov, E.M. OrthoDB v10: Sampling the diversity of animal, plant, fungal, protist, bacterial and viral genomes for evolutionary and functional annotations of orthologs. Nucleic Acids Res. 2019, 47, D807–D811. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Robinson, M.D.; McCarthy, D.; Smyth, G. edgeR: A bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 2010, 26, 139–140. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Li, B.; Dewey, C.N. RSEM: Accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinform. 2011, 12, 323. [Google Scholar] [CrossRef] [Green Version]
  31. Keeling, C.I.; Blomquist, G.J.; Tittiger, C. Coordinated gene expression for pheromone biosynthesis in the pine engraver beetle, Ips pini (Coleoptera: Scolytidae). Naturwissenschaften 2004, 91, 324–328. [Google Scholar] [CrossRef]
  32. Bellés, X.; Martín, D.; Piulachs, M.-D. The mevalonate pathway and the synthesis of juvenile hormone in insects. Annu. Rev. Ѐntomol. 2005, 50, 181–199. [Google Scholar] [CrossRef] [Green Version]
  33. González-Caballero, N.; Rodríguez-Vega, A.; Dias-Lopes, G.; Valenzuela, J.G.; Ribeiro, J.M.; Carvalho, P.C.; Valente, R.H.; Brazil, R.P.; Cuervo, P. Expression of the mevalonate pathway enzymes in the Lutzomyia longipalpis (Diptera: Psychodidae) sex pheromone gland demonstrated by an integrated proteomic approach. J. Proteom. 2014, 96, 117–132. [Google Scholar] [CrossRef] [Green Version]
  34. Beedle, A.; Walton, M.; Goodwin, T. Isoprenoid biosynthesis in aseptic larvae of Calliphora erythrocephala. Insect Biochem. 1975, 5, 465–472. [Google Scholar] [CrossRef]
  35. Tittiger, C.; Blomquist, G. Pheromone production in pine bark beetles. Adv. Insect Physiol. 2016, 50, 235–263. [Google Scholar]
  36. Tillman, J.A.; Holbrook, G.L.; Dallara, P.L.; Schal, C.; Wood, D.L.; Blomquist, G.J.; Seybold, S.J. Endocrine regulation of de novo aggregation pheromone biosynthesis in the pine engraver, Ips pini (Say) (Coleoptera: Scolytidae). Insect Biochem. Mol. Biol. 1998, 28, 705–715. [Google Scholar] [CrossRef]
  37. Stroumbakis, N.D.; Li, Z.; Tolias, P.P. A homolog of human transcription factor NF-X1 encoded by the Drosophila shuttle craft gene is required in the embryonic central nervous system. Mol. Cell. Biol. 1996, 16, 192–201. [Google Scholar] [CrossRef] [Green Version]
  38. Goodman, C.S. The likeness of being: Phylogenetically conserved molecular mechanisms of growth cone guidance. Cell 1994, 78, 353–356. [Google Scholar] [CrossRef]
  39. Yamaguchi, M. The transcriptional regulation of regucalcin gene expression. Mol. Cell. Biochem. 2011, 346, 147–171. [Google Scholar] [CrossRef]
  40. Barrett, A.; Rawlings, N. Families and Clans of Serine Peptidases. Arch. Biochem. Biophys. 1995, 318, 247–250. [Google Scholar] [CrossRef]
  41. Rao, M.B.; Tanksale, A.M.; Ghatge, M.S.; Deshpande, V.V. Molecular and biotechnological aspects of microbial proteases. Microbiol. Mol. Biol. Rev. 1998, 62, 597–635. [Google Scholar] [CrossRef] [Green Version]
  42. Perkin, L.; Elpidina, E.; Oppert, B. RNA interference and dietary inhibitors induce a similar compensation response in Tribolium castaneum larvae. Insect Mol. Biol. 2017, 26, 35–45. [Google Scholar] [CrossRef]
  43. Oliveira-Neto, O.B.; Batista, J.A.; Rigden, D.J.; Fragoso, R.R.; Silva, R.O.; Gomes, E.A.; Franco, O.L.; Dias, S.C.; Cordeiro, C.M.; Monnerat, R.G.; et al. A diverse family of serine proteinase genes expressed in cotton boll weevil (Anthonomus grandis): Implications for the design of pest-resistant transgenic cotton plants. Insect Biochem. Mol. Biol. 2004, 34, 903–918. [Google Scholar] [CrossRef] [PubMed]
  44. Chowdhury, M.; Li, C.-F.; He, Z.; Lu, Y.; Liu, X.-S.; Wang, Y.-F.; Ip, Y.T.; Strand, M.R.; Yu, X.-Q. Toll family members bind multiple Spätzle proteins and activate antimicrobial peptide gene expression in Drosophila. J. Biol. Chem. 2019, 294, 10172–10181. [Google Scholar] [CrossRef] [PubMed]
  45. Liu, J.; Swevers, L.; Iatrou, K.; Huvenne, H.; Smagghe, G. Bombyx mori DNA/RNA non-specific nuclease: Expression of isoforms in insect culture cells, subcellular localization and functional assays. J. Insect Physiol. 2012, 58, 1166–1176. [Google Scholar] [CrossRef] [PubMed]
  46. Johansson, B.G.; Jones, T.M. The role of chemical communication in mate choice. Biol. Rev. 2007, 82, 265–289. [Google Scholar] [CrossRef] [PubMed]
  47. Rantala, M.J.; Jokinen, I.; Kortet, R.; Vainikka, A.; Suhonen, J. Do pheromones reveal male immunocompetence? Proc. Biol. Sci. 2002, 269, 1681–1685. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Sugumaran, M. Comparative biochemistry of eumelanogenesis and the protective roles of phenoloxidase and melanin in insects. Pigment. Cell Res. 2002, 15, 2–9. [Google Scholar] [CrossRef]
  49. Stoehr, A. Costly melanin ornaments: The importance of taxon? Func. Ecol. 2006, 20, 276–281. [Google Scholar] [CrossRef]
  50. Arakane, Y.; Muthukrishnan, S.; Beeman, R.W.; Kanost, M.R.; Kramer, K.J. Laccase 2 is the phenoloxidase gene required for beetle cuticle tanning. Proc. Natl. Acad. Sci. USA 2005, 102, 11337–11342. [Google Scholar] [CrossRef] [Green Version]
  51. Arakane, Y.; Lomakin, J.; Beeman, R.W.; Muthukrishnan, S.; Gehrke, S.H.; Kanost, M.R.; Kramer, K.J. Molecular and functional analyses of amino acid decarboxylases involved in cuticle tanning in Tribolium castaneum. J. Biol. Chem. 2009, 284, 16584–16594. [Google Scholar] [CrossRef] [Green Version]
  52. Lukacsovich, T.; Yuge, K.; Awano, W.; Asztalos, Z.; Kondo, S.; Juni, N.; Yamamoto, D. The ken and barbie gene encoding a putative transcription factor with a BTB domain and three zinc finger motifs functions in terminalia development of Drosophila. Arch. Insect Biochem. Physiol. 2003, 54, 77–94. [Google Scholar] [CrossRef] [PubMed]
Figure 1. BLAST top species hits for the annotated Anthonomus grandis grandis de novo transcriptome. Top hit species are listed on the y-axis, and the number of hits for each species is on the x-axis.
Figure 1. BLAST top species hits for the annotated Anthonomus grandis grandis de novo transcriptome. Top hit species are listed on the y-axis, and the number of hits for each species is on the x-axis.
Insects 12 00893 g001
Figure 2. Gene ontology categories for the annotated Anthonomus grandis grandis de novo transcriptome indicating (a) molecular function (MF), (b) biological process (BP), and (c) cellular component (CC).
Figure 2. Gene ontology categories for the annotated Anthonomus grandis grandis de novo transcriptome indicating (a) molecular function (MF), (b) biological process (BP), and (c) cellular component (CC).
Insects 12 00893 g002
Figure 3. Volcano Plot showing the 2479 significant differentially expressed genes between Anthonomus grandis grandis males that produced pheromone and those that did not produce pheromone. Of these genes, 1234 were up-regulated (red; FDR < 0.05, logFC > 0.01) and 1515 genes were down-regulated (green; FDR < 0.05, logFC < 0.01). Grey dots represent genes with FDR < 0.01, and black dots represent the majority of genes with FDR > 0.05. The blue vertical lines show the −1 and 1 log fold change threshold. The y-axis shows the negative log10 FDR value, and the x-axis shows the log2 fold change in gene expression.
Figure 3. Volcano Plot showing the 2479 significant differentially expressed genes between Anthonomus grandis grandis males that produced pheromone and those that did not produce pheromone. Of these genes, 1234 were up-regulated (red; FDR < 0.05, logFC > 0.01) and 1515 genes were down-regulated (green; FDR < 0.05, logFC < 0.01). Grey dots represent genes with FDR < 0.01, and black dots represent the majority of genes with FDR > 0.05. The blue vertical lines show the −1 and 1 log fold change threshold. The y-axis shows the negative log10 FDR value, and the x-axis shows the log2 fold change in gene expression.
Insects 12 00893 g003
Figure 4. Mevalonate pathway in insects. Genes that were up-regulated in pheromone-producing Anthonomus grandis grandis are highlighted in red and accompanied by an up-facing arrow. This figure was modified from Bellès et al., 2005 [32].
Figure 4. Mevalonate pathway in insects. Genes that were up-regulated in pheromone-producing Anthonomus grandis grandis are highlighted in red and accompanied by an up-facing arrow. This figure was modified from Bellès et al., 2005 [32].
Insects 12 00893 g004
Table 1. Anthonomus grandis grandis de novo transcriptome statistics and BUSCO scores.
Table 1. Anthonomus grandis grandis de novo transcriptome statistics and BUSCO scores.
Raw sequences352,189,343
Total transcripts195,128
Average sequence length922.89 bp
N502040 bp
BUSCO
Complete99.27%
Complete single copy97.81%
Complete duplicated1.46%
Fragmented0.59%
Missing0.15%
Underline shows that the information below describes the BUSCO score.
Table 2. List of genes consistent with pheromone production and biosynthesis that were up-regulated in pheromone-producing Anthonomus grandis grandis compared to non-pheromone-producing boll weevils.
Table 2. List of genes consistent with pheromone production and biosynthesis that were up-regulated in pheromone-producing Anthonomus grandis grandis compared to non-pheromone-producing boll weevils.
Sequence NameSequence DescriptionBLAST Top HitFold ChangeUP/DOWNp-ValueFDR
Mevalonate pathway or associated with pheromone production and biosynthesis
TRINITY_DN11108hydroxymethylglutaryl-CoA synthase 1 (HMG-S)D. ponderosae23.081.45 × 10−372.67 × 10−34
TRINITY_DN13668dolichyl pyrophosphate Man9GlcNAc2 α-1,3-glucosyltransferaseD. ponderosae2.014.44 × 10−051.82 × 10−03
TRINITY_DN5050farnesyl pyrophosphate synthase-likeI. pini11.771.30 × 10−044.59 × 10−03
TRINITY_DN5422phosphomevalonate kinaseD. ponderosae3.059.80 × 10−121.51 × 10−09
TRINITY_DN8104farnesyl pyrophosphate synthaseI. pini34.051.10 × 10−055.33 × 10−04
TRINITY_DN8499isopentenyl-diphosphate Delta-isomerase 1A. glabripennis37.641.57 × 10−091.69 × 10−07
TRINITY_DN8733-hydroxy-3-methylglutaryl-coenzyme A reductase-like (HMG-R)I. paraconfusus16.579.18 × 10−299.18 × 10−26
TRINITY_DN4352ATP-citrate synthaseD. ponderosae6.661.89 × 10−253.07 × 10−32
Juvenile hormone
TRINITY_DN11992juvenile hormone inducible proteinD. ponderosae7.757.40 × 10−131.39 × 10−10
TRINITY_DN4315juvenile hormone esterase-likeD. ponderosae2.985.14 × 10−117.24 × 10−09
TRINITY_DN874juvenile hormone esterase-likeD. ponderosae2.711.99 × 10−081.79 × 10−06
TRINITY_DN4184juvenile hormone binding proteinR. ferrugineus5.642.69 × 10−082.34 × 10−06
TRINITY_DN9052juvenile hormone inducible proteinD. ponderosae4.537.72 × 10−074.94 × 10−05
TRINITY_DN115965putative juvenile hormone inducible proteinD. ponderosae3.261.10 × 10−055.33 × 10−04
TRINITY_DN11903juvenile hormone acid O-methyltransferaseD. ponderosae2.895.87 × 10−052.31 × 10−03
Fatty acid metabolism
TRINITY_DN58fatty acid synthaseD. ponderosae4.922.20 × 10−261.74 × 10−23
TRINITY_DN5986glycerol-3-phosphate phosphatase-likeS. oryzae7.503.68 × 10−262.89 × 10−23
Neurological and hormonal regulation
TRINITY_DN4411protein shuttle craft likeD. ponderosae5.383.60 × 10−293.64 × 10−26
TRINITY_DN72610regucalcin-likeD. ponderosae5.281.45 × 10−251.11 × 10−22
Table 3. List of top up- and down-regulated genes in pheromone-producing Anthonomus grandis grandis compared to non-pheromone-producing boll weevils. Rows highlighted in light grey represent genes that were DE but not part of the top 50 DEGs.
Table 3. List of top up- and down-regulated genes in pheromone-producing Anthonomus grandis grandis compared to non-pheromone-producing boll weevils. Rows highlighted in light grey represent genes that were DE but not part of the top 50 DEGs.
Sequence NameSequence DescriptionBlast Top HitFold ChangeUP/DOWNp-ValueFDR
Peptidase activity
TRINITY_DN62502transmembrane protease serine 9-likeZ. cucurbitae27.181.17 × 10−537.28 × 10−50
TRINITY_DN35065cathepsin L1-likeR. ferrugineus14.423.42 × 10−511.73 × 10−47
TRINITY_DN2409carboxypeptidase B-likeD. ponderosae12.749.29 × 10412.28 × 1037
TRINITY_DN5051venom serine carboxypeptidase-likeS. oryzae26.912.87 × 10−313.46 × 10−28
TRINITY_DN13263trypsin-like serine prtoeaseA. grandis15.384.09 × 10−263.18 × 10−23
TRINITY_DN4893venom serine protease-likeD. ponderosae10.492.70 × 10−221.45 × 10−19
TRINITY_DN34699trypsin alpha 3-like (Agser2p)—induced by eatingA. grandis3.836.99 × 10−151.76 × 10−12
TRINITY_DN35978brachyurin-like (Agser9p)A. grandis3.369.32 × 10−121.46 × 10−09
TRINITY_DN4825brachyurin-like (Agser5p)—induced by eatingA. grandis2.265.60 × 10−062.93 × 10−04
TRINITY_DN49721brachyurin-like (Agser9p)A. grandis3.569.43 × 10−131.74 × 10−10
TRINITY_DN63273trypsin alpha-3-like (Agser2p)—induced by feedingA. grandis3.572.63 × 10−135.21 × 10−11
TRINITY_DN69718serine protease (Agser12p)A. grandis2.272.86 × 10−051.24 × 10−03
Immune response
TRINITY_DN4819protein spaetzle 3D. ponderosae25.783.30 × 10−451.11 × 10−41
TRINITY_DN3567DNA/RNA non-specific nuclease 1A. grandis13.543.13 × 10−396.34 × 10−36
Sexual maturity
TRINITY_DN10824transcription factor Ken 7.861.35 × 10−444.39 × 10−41
Cuticle tanning and sclerotization
TRINITY_DN16013pupal cuticle protein 36-likeD. ponderosae89.536.46 × 10−892.61 × 10−84
TRINITY_DN11117cuticular protein 141D. ponderosae13.541.68 × 10−682.72 × 10−64
TRINITY_DN3947chitin deacetylase 4D. ponderosae20.481.14 × 10−631.16 × 10−59
TRINITY_DN4706chitin synthase 1A. grandis5.845.48 × 10−358.37 × 10−32
TRINITY_DN1999cuticle protein 19.8-likeS. oryzae15.151.49 × 10−321.89 × 10−29
TRINITY_DN4303cuticular protein analogous to peritrophinsD. ponderosae5.841.53 × 10−321.91 × 10−29
TRINITY_DN2312laccase 2A. eugenii18.719.87 × 10−832.66 × 10−78
TRINITY_DN22343dopamine D2-like receptorS. oryzae4.854.95 × 10−094.90 × 10−07
TRINITY_DN6738dopamine N-acetyltransferaseD. ponderosae2.034.40 × 10−041.30 × 10−02
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Perkin, L.C.; Perez, J.L.; Suh, C.P.-C. The Identification of Boll Weevil, Anthonomus grandis grandis (Coleoptera: Curculionidae), Genes Involved in Pheromone Production and Pheromone Biosynthesis. Insects 2021, 12, 893. https://doi.org/10.3390/insects12100893

AMA Style

Perkin LC, Perez JL, Suh CP-C. The Identification of Boll Weevil, Anthonomus grandis grandis (Coleoptera: Curculionidae), Genes Involved in Pheromone Production and Pheromone Biosynthesis. Insects. 2021; 12(10):893. https://doi.org/10.3390/insects12100893

Chicago/Turabian Style

Perkin, Lindsey C., Jose L. Perez, and Charles P.-C. Suh. 2021. "The Identification of Boll Weevil, Anthonomus grandis grandis (Coleoptera: Curculionidae), Genes Involved in Pheromone Production and Pheromone Biosynthesis" Insects 12, no. 10: 893. https://doi.org/10.3390/insects12100893

APA Style

Perkin, L. C., Perez, J. L., & Suh, C. P. -C. (2021). The Identification of Boll Weevil, Anthonomus grandis grandis (Coleoptera: Curculionidae), Genes Involved in Pheromone Production and Pheromone Biosynthesis. Insects, 12(10), 893. https://doi.org/10.3390/insects12100893

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop