RNA-Seq Validation of RNAi Identifies Additional Gene Connectivity in Tribolium castaneum (Coleoptera: Tenebrionidae).

RNA interference (RNAi) is a functional genomics tool to correlate genotype and phenotype by delivering targeted, gene-specific, and complementary dsRNA into a host via injection, feeding, or other means in order to reduce gene expression. In the red flour beetle, Tribolium castaneum, RNAi has been successful via injected dsRNA at all life stages. Traditionally, successful transcript knockdown has been quantified by qPCR on a gene-by-gene basis, where only expression of the target gene and normalization genes are evaluated. In this study, RNA-Seq was used to quantify transcript expression in larvae injected with dsRNA for aspartate 1-decarboxylase (ADC), which gives a reliable phenotype of an adult with a black cuticle instead of the wild-type red-brown. ANOVA of control, mock-injected, and ADC-dsRNA injected larvae indicated that target gene expression was significantly (P = 0.002) reduced 4-fold, and the black phenotype was achieved in all adults injected with ADC-dsRNA as larvae. In a pairwise analysis, significant (P < 0.05) differential expression of other genes in ADC-injected larvae suggested connections between gene pathways. One gene, dopamine receptor 2, was increased in expression 227-fold (P = 0.025), presumably connected to previous data that showed a reduction in expression of ADC results in increased levels of dopamine. To evaluate the hypothesis that increased dopamine levels can affect mobility, T. castaneum adults injected with ADC-dsRNA as larvae were significantly impaired in movement tests compared to controls, similar to black mutants in Drosophila melanogaster. The data demonstrate that RNA-Seq can reveal gene connectivity and provide more complete data validation and analysis compared to qPCR.

Nematode-resistance loci in upland cotton genomes are associated with structural differences.

Reniform and root-knot nematode are two of the most destructive pests of conventional upland cotton, Gossypium hirsutum L., and continue to be a major threat to cotton fiber production in semiarid regions of the Southern United States and Central America. Fortunately, naturally occurring tolerance to these nematodes has been identified in the Pima cotton species (Gossypium barbadense) and several upland cotton varieties (G. hirsutum), which has led to a robust breeding program that has successfully introgressed and stacked these independent resistant traits into several upland cotton lineages with superior agronomic traits, e.g. BAR 32-30 and BARBREN-713. This work identifies the genomic variations of these nematode-tolerant accessions by comparing their respective genomes to the susceptible, high-quality fiber-producing parental line of this lineage: Phytogen 355 (PSC355). We discover several large genomic differences within marker regions that harbor putative resistance genes as well as expression mechanisms shared by the two resistant lines, with respect to the susceptible PSC355 parental line. This work emphasizes the utility of whole-genome comparisons as a means of elucidating large and small nuclear differences by lineage and phenotype.

The Transcriptomic Response of the Boll Weevil, Anthonomus grandis grandis Boheman (Coleoptera: Curculionidae), following Exposure to the Organophosphate Insecticide Malathion.

Insecticide tolerance and resistance have evolved countless times in insect systems. Molecular drivers of resistance include mutations in the insecticide target site and/or gene duplication, and increased gene expression of detoxification enzymes. The boll weevil, Anthonomus grandis grandis Boheman (Coleoptera: Curculionidae), is a pest of commercial cotton and has developed resistance in the field to several insecticides; however, the current organophosphate insecticide, malathion, used by USA eradication programs remains effective despite its long-term use. Here, we present findings from an RNA-seq experiment documenting gene expression post-exposure to field-relevant concentrations of malathion, which was used to provide insight on the boll weevil's continued susceptibility to this insecticide. Additionally, we incorporated a large collection of boll weevil whole-genome resequencing data from nearly 200 individuals collected from three geographically distinct areas to determine SNP allele frequency of the malathion target site, as a proxy for directional selection in response to malathion exposure. No evidence was found in the gene expression data or SNP data consistent with a mechanism of enhanced tolerance or resistance adaptation to malathion in the boll weevil. Although this suggests continued effectiveness of malathion in the field, we identified important temporal and qualitative differences in gene expression between weevils exposed to two different concentrations of malathion. We also identified several tandem isoforms of the detoxifying esterase B1 and glutathione S-transferases, which are putatively associated with organophosphate resistance.

Insight into weevil biology from a reference quality genome of the boll weevil, Anthonomus grandis grandis Boheman (Coleoptera: Curculionidae).

The boll weevil, Anthonomus grandis grandis Boheman, is one of the most historically impactful insects due to its near destruction of the US cotton industry in the early 20th century. Contemporary efforts to manage this insect primarily use pheromone baited traps for detection and organophosphate insecticides for control, but this strategy is not sustainable due to financial and environmental costs. We present a high-quality boll weevil genome assembly, consisting of 306 scaffolds with approximately 24,000 annotated genes, as a first step in the identification of gene targets for novel pest control. Gene content and transposable element distribution are similar to those found in other Curculionidae genomes; however, this is the most contiguous and only assembly reported to date for a member in the species-rich genus Anthonomus. Transcriptome profiles across larval, pupal, and adult life stages led to identification of several genes and gene families that could present targets for novel control strategies.

A New qPCR Assay for the Rapid Diagnosis of Anthonomus grandis Subspecies.

Rapid and accurate identification of Anthonomus grandis subspecies is crucial for effective management and eradication. Current diagnostic methods have limitations in terms of time to diagnosis (up to seven days) and can yield ambiguous results. Here, we present the validation of a custom TaqMan SNP Genotyping Assay for the rapid and accurate identification of A. grandis grandis (boll weevil) and A. g. thurberiae (thurberia weevil) subspecies. To validate the assay, we conducted three main experiments: (1) a sensitivity test to determine the DNA concentration range at which the assay performs, (2) a non-target specificity test to ensure no amplification in non-target weevils (false positives), and (3) an accuracy test comparing the results of the new assay to previously established methods. These experiments were carried out in parallel at three independent facilities to confirm the robustness of the assay to variations in equipment and personnel. We used DNA samples from various sources, including field-collected specimens, museum specimens, and previously isolated DNA. The assay demonstrated high sensitivity (PCR success with ≥0.05 ng/µL DNA template), specificity (0.02 false positive rate), and accuracy (97.7%) in diagnosing boll weevil and thurberia weevil subspecies. The entire workflow, including DNA extraction, assay preparation, PCR run time, and data analysis, can be completed within a single workday (7-9 h) by a single technician. The deployment of this assay as a diagnostic tool could benefit boll weevil management and eradication programs by enabling same-day diagnosis of trap-captured or intercepted weevil specimens. Furthermore, it offers a more reliable method for identifying unknown specimens, contributing to the overall effectiveness of boll weevil research and control efforts.

Variants in the Mitochondrial Genome Sequence of Rhyzopertha dominica (Fabricius) (Coleoptera: Bostrycidae).

The lesser grain borer, Rhyzopertha dominica, is a coleopteran pest of stored grains and is mainly controlled by phosphine fumigation, but the increase in phosphine-resistant populations threatens efficacy. Some phosphine-resistant insects have reduced respiration, and thus studying the mitochondrial genome may provide additional information regarding resistance. Genomic DNA from an inbred laboratory strain of R. dominica was extracted and sequenced with both short (Illumina) and long (Pacific Biosciences) read technologies for whole genome sequence assembly and annotation. Short read sequences were assembled and annotated by open software to identify mitochondrial sequences, and the assembled sequence was manually annotated and verified by long read sequences. The mitochondrial genome sequence for R. dominica had a total length of 15,724 bp and encoded 22 trna genes, 2 rRNA genes, 13 protein coding genes (7 nad subunits, 3 cox, 2 atp, and 1 cytB), flanked by a long control region. We compared our predicted mitochondrial genome to that of another from a R. dominica strain from Jingziguan (China). While there was mostly agreement between the two assemblies, key differences will be further examined to determine if mutations in populations are related to insecticide control pressure, mainly that of phosphine. Differences in sequence data, assembly, and annotation also may result in different genome interpretations.

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

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.

Genome assembly of two nematode-resistant cotton lines (Gossypium hirsutum L.).

Upland cotton (Gossypium hirsutum L.) is susceptible to damage by the root-knot and the reniform nematodes, causing yield losses greater than 4% annually in the United States. In addition, these nematodes are synergistic with seeding disease and root rot pathogens that exacerbate diseases and subsequent yield losses. Production practices to minimize nematode damage include crop rotation and nematicides, but these techniques need to be repeated and are expensive. The use of resistant cultivars is deemed the most effective and economical approach for managing nematodes in cotton. Here, we describe the genomes of two nematode-resistant lines of cotton, BARBREN-713 and BAR 32-30. These genomes may expedite the development of DNA markers that can be used to efficiently introduce nematode resistance into commercially valuable Upland lines.

Population genomics and phylogeography of the boll weevil, Anthonomus grandis Boheman (Coleoptera: Curculionidae), in the United States, northern Mexico, and Argentina.

The boll weevil, Anthonomus grandis Boheman (Coleoptera: Curculionidae), is an important pest of commercial cotton across the Americas. In the United States, eradication of this species is complicated by re-infestations of areas where eradication has been previously successful and by the existence of morphologically similar variants that can confound identification efforts. To date, no study has applied a high-throughput sequencing approach to better understand the population genetic structure of the boll weevil. Furthermore, only a single study has investigated genetic relationships between populations in North and South America. We used double digest restriction site-associated DNA sequencing (ddRADseq) to resolve the population genomic structure of the boll weevil in the southern United States, northern Mexico, and Argentina. Additionally, we assembled the first complete mitochondrial genome for this species and generated a preliminary whole genome assembly, both of which were used to improve the identification of informative loci. Downstream analyses revealed two main lineages-one consisting of populations found geographically west of the Sierra Madre Occidental mountain range and the second consisting of populations found to the east-were revealed, and both were sub-structured. Population geographic structure was consistent with the isolation by distance model, indicating that geogrpahic distance is likely a primary mechanism driving divergence in this species. Boll weevil populations from Argentina were found to be more closely related to the eastern lineage, suggesting a recent colonization of South America by the eastern lineage, but additional sampling across Mexico, Central America and South America is needed to further clarify their origin. Finally, we uncovered an instance of population turnover or replacement, highlighting the temporal instability of population structure.

Transcriptome analysis of life stages of the house cricket, Acheta domesticus, to improve insect crop production.

To develop genetic resources for the improvement of insects as food, we sequenced transcripts from embryos, one-day hatchlings, three nymphal stages, and male and female adults of the house cricket, Acheta domesticus. A draft transcriptome was assembled from more than 138 million sequences combined from all life stages and sexes. The draft transcriptome assembly contained 45,866 contigs, and more than half were similar to sequences at NCBI (e value < e-3). The highest sequence identity was found in sequences from the termites Cryptotermes secundus and Zootermopsis nevadensis. Sequences with identity to Gregarina niphandrodes suggest that these crickets carry the parasite. Among all life stages, there were 5,042 genes with differential expression between life stages (significant at p < 0.05). An enrichment analysis of gene ontology terms from each life stage or sex highlighted genes that were important to biological processes in cricket development. We further characterized genes that may be important in future studies of genetically modified crickets for improved food production, including those involved in RNA interference, and those encoding prolixicin and hexamerins. The data represent an important first step in our efforts to provide genetically improved crickets for human consumption and livestock feed.

Gene expression in Tribolium castaneum life stages: Identifying a species-specific target for pest control applications.

The red flour beetle, Tribolium castaneum, is a major agricultural pest of post-harvest products and stored grain. Control of T. castaneum in stored products and grain is primarily by fumigants and sprays, but insecticide resistance is a major problem, and new control strategies are needed. T. castaneum is a genetic model for coleopterans, and the reference genome can be used for discovery of candidate gene targets for molecular-based control, such as RNA interference. Gene targets need to be pest specific, and ideally, they are expressed at low levels for successful control. Therefore, we sequenced the transcriptome of four major life stages of T. castaneum, sorted data into groups based on high or low expression levels, and compared relative gene expression among all life stages. We narrowed our candidate gene list to a cuticle protein gene (CPG) for further analysis. We found that the CPG sequence was unique to T. castaneum and expressed only in the larval stage. RNA interference targeting CPG in newly-emerged larvae caused a significant (p < 0.05) decrease in CPG expression (1,491-fold) compared to control larvae and 64% mortality over 18 d. RNA-Seq of survivors after 18 d identified changes in the expression of other genes as well, including 52 long noncoding RNAs. Expression of three additional cuticle protein genes were increased and two chitinase genes were decreased in response to injection of CPG dsRNA. The data demonstrate that RNA-Seq can identify genes important for insect survival and thus may be used to develop novel biologically-based insect control products.

A model species for agricultural pest genomics: the genome of the Colorado potato beetle, Leptinotarsa decemlineata (Coleoptera: Chrysomelidae).

The Colorado potato beetle is one of the most challenging agricultural pests to manage. It has shown a spectacular ability to adapt to a variety of solanaceaeous plants and variable climates during its global invasion, and, notably, to rapidly evolve insecticide resistance. To examine evidence of rapid evolutionary change, and to understand the genetic basis of herbivory and insecticide resistance, we tested for structural and functional genomic changes relative to other arthropod species using genome sequencing, transcriptomics, and community annotation. Two factors that might facilitate rapid evolutionary change include transposable elements, which comprise at least 17% of the genome and are rapidly evolving compared to other Coleoptera, and high levels of nucleotide diversity in rapidly growing pest populations. Adaptations to plant feeding are evident in gene expansions and differential expression of digestive enzymes in gut tissues, as well as expansions of gustatory receptors for bitter tasting. Surprisingly, the suite of genes involved in insecticide resistance is similar to other beetles. Finally, duplications in the RNAi pathway might explain why Leptinotarsa decemlineata has high sensitivity to dsRNA. The L. decemlineata genome provides opportunities to investigate a broad range of phenotypes and to develop sustainable methods to control this widely successful pest.

Gene Disruption Technologies Have the Potential to Transform Stored Product Insect Pest Control.

Stored product insects feed on grains and processed commodities manufactured from grain post-harvest, reducing the nutritional value and contaminating food. Currently, the main defense against stored product insect pests is the pesticide fumigant phosphine. Phosphine is highly toxic to all animals, but is the most effective and economical control method, and thus is used extensively worldwide. However, many insect populations have become resistant to phosphine, in some cases to very high levels. New, environmentally benign and more effective control strategies are needed for stored product pests. RNA interference (RNAi) may overcome pesticide resistance by targeting the expression of genes that contribute to resistance in insects. Most data on RNAi in stored product insects is from the coleopteran genetic model, Tribolium castaneum, since it has a strong RNAi response via injection of double stranded RNA (dsRNA) in any life stage. Additionally, Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) technology has been suggested as a potential resource for new pest control strategies. In this review we discuss background information on both gene disruption technologies and summarize the advances made in terms of molecular pest management in stored product insects, mainly T. castaneum, as well as complications and future needs.

Venom gland components of the ectoparasitoid wasp, Anisopteromalus calandrae.

The wasp Anisopteromalus calandrae is a small ectoparasitoid that attacks stored product pest beetle larvae that develop inside grain kernels, and is thus a potential insect control tool. The components of A. calandrae venom have not been studied, but venom from other organisms contains proteins with potential applications, such as pest management tools and treatments for human diseases. We dissected female A. calandrae and collected venom and associated glands. Using high throughput sequencing, a venom gland transcriptome was assembled that contained 45,432 contigs, 25,726 of which had BLASTx hits. The majority of hits were to Nasonia vitripennis, an ectoparasitoid from the same taxonomic family, as well as other bees, wasps, and ants. Gene ontology grouped sequences into eleven molecular functions, among which binding and catalytic activity had the most representatives. In this study, we highlighted the most abundant sequences, including those that are likely the functional components of the venom. Specifically, we focused on genes encoding proteins potentially involved in host developmental arrest, disrupting the host immune system, host paralysis, and transcripts that support these functions. Our report is the first to characterize components of the A. calandrae venom gland that may be useful as control tools for insect pests and other applications.