A chromosome-level genome assembly of the soybean pod borer: insights into larval transcriptional response to transgenic soybean expressing the pesticidal Cry1Ac protein.

BACKGROUND: Genetically modified (GM) crop plants with transgenic expression of Bacillus thuringiensis (Bt) pesticidal proteins are used to manage feeding damage by pest insects. The durability of this technology is threatened by the selection for resistance in pest populations. The molecular mechanism(s) involved in insect physiological response or evolution of resistance to Bt is not fully understood. RESULTS: To investigate the response of a susceptible target insect to Bt, the soybean pod borer, Leguminivora glycinivorella (Lepidoptera: Tortricidae), was exposed to soybean, Glycine max, expressing Cry1Ac pesticidal protein or the non-transgenic parental cultivar. Assessment of larval changes in gene expression was facilitated by a third-generation sequenced and scaffolded chromosome-level assembly of the L. glycinivorella genome (657.4 Mb; 27 autosomes + Z chromosome), and subsequent structural annotation of 18,197 RefSeq gene models encoding 23,735 putative mRNA transcripts. Exposure of L. glycinivorella larvae to transgenic Cry1Ac G. max resulted in prediction of significant differential gene expression for 204 gene models (64 up- and 140 down-regulated) and differential splicing among isoforms for 10 genes compared to unexposed cohorts. Differentially expressed genes (DEGs) included putative peritrophic membrane constituents, orthologs of Bt receptor-encoding genes previously linked or associated with Bt resistance, and those involved in stress responses. Putative functional Gene Ontology (GO) annotations assigned to DEGs were significantly enriched for 36 categories at GO level 2, respectively. Most significantly enriched cellular component (CC), biological process (BP), and molecular function (MF) categories corresponded to vacuolar and microbody, transport and metabolic processes, and binding and reductase activities. The DEGs in enriched GO categories were biased for those that were down-regulated (≥ 0.783), with only MF categories GTPase and iron binding activities were bias for up-regulation genes. CONCLUSIONS: This study provides insights into pathways and processes involved larval response to Bt intoxication, which may inform future unbiased investigations into mechanisms of resistance that show no evidence of alteration in midgut receptors.

A draft Diabrotica virgifera virgifera genome: insights into control and host plant adaption by a major maize pest insect.

BACKGROUND: Adaptations by arthropod pests to host plant defenses of crops determine their impacts on agricultural production. The larval host range of western corn rootworm, Diabrotica virgifera virgifera (Coleoptera: Chrysomelidae), is restricted to maize and a few grasses. Resistance of D. v. virgifera to crop rotation practices and multiple insecticides contributes to its status as the most damaging pest of cultivated maize in North America and Europe. The extent to which adaptations by this pest contributes to host plant specialization remains unknown. RESULTS: A 2.42 Gb draft D. v. virgifera genome, Dvir_v2.0, was assembled from short shotgun reads and scaffolded using long-insert mate-pair, transcriptome and linked read data. K-mer analysis predicted a repeat content of ≥ 61.5%. Ortholog assignments for Dvir_2.0 RefSeq models predict a greater number of species-specific gene duplications, including expansions in ATP binding cassette transporter and chemosensory gene families, than in other Coleoptera. A majority of annotated D. v. virgifera cytochrome P450s belong to CYP4, 6, and 9 clades. A total of 5,404 transcripts were differentially-expressed between D. v. virgifera larvae fed maize roots compared to alternative host (Miscanthus), a marginal host (Panicum virgatum), a poor host (Sorghum bicolor) and starvation treatments; Among differentially-expressed transcripts, 1,908 were shared across treatments and the least number were between Miscanthus compared to maize. Differentially-expressed transcripts were enriched for putative spliceosome, proteosome, and intracellular transport functions. General stress pathway functions were unique and enriched among up-regulated transcripts in marginal host, poor host, and starvation responses compared to responses on primary (maize) and alternate hosts. CONCLUSIONS: Manual annotation of D. v. virgifera Dvir_2.0 RefSeq models predicted expansion of paralogs with gene families putatively involved in insecticide resistance and chemosensory perception. Our study also suggests that adaptations of D. v. virgifera larvae to feeding on an alternate host plant invoke fewer transcriptional changes compared to marginal or poor hosts. The shared up-regulation of stress response pathways between marginal host and poor host, and starvation treatments may reflect nutrient deprivation. This study provides insight into transcriptomic responses of larval feeding on different host plants and resources for genomic research on this economically significant pest of maize.

Up-regulation of apoptotic- and cell survival-related gene pathways following exposures of western corn rootworm to B. thuringiensis crystalline pesticidal proteins in transgenic maize roots.

BACKGROUND: Resistance of pest insect species to insecticides, including B. thuringiensis (Bt) pesticidal proteins expressed by transgenic plants, is a threat to global food security. Despite the western corn rootworm, Diabrotica virgifera virgifera, being a major pest of maize and having populations showing increasing levels of resistance to hybrids expressing Bt pesticidal proteins, the cell mechanisms leading to mortality are not fully understood. RESULTS: Twenty unique RNA-seq libraries from the Bt susceptible D. v. virgifera inbred line Ped12, representing all growth stages and a range of different adult and larval exposures, were assembled into a reference transcriptome. Ten-day exposures of Ped12 larvae to transgenic Bt Cry3Bb1 and Gpp34/Tpp35Ab1 maize roots showed significant differential expression of 1055 and 1374 transcripts, respectively, compared to cohorts on non-Bt maize. Among these, 696 were differentially expressed in both Cry3Bb1 and Gpp34/Tpp35Ab1 maize exposures. Differentially-expressed transcripts encoded protein domains putatively involved in detoxification, metabolism, binding, and transport, were, in part, shared among transcripts that changed significantly following exposures to the entomopathogens Heterorhabditis bacteriophora and Metarhizium anisopliae. Differentially expressed transcripts in common between Bt and entomopathogen treatments encode proteins in general stress response pathways, including putative Bt binding receptors from the ATP binding cassette transporter superfamily. Putative caspases, pro- and anti-apoptotic factors, as well as endoplasmic reticulum (ER) stress-response factors were identified among transcripts uniquely up-regulated following exposure to either Bt protein. CONCLUSIONS: Our study suggests that the up-regulation of genes involved in ER stress management and apoptotic progression may be important in determining cell fate following exposure of susceptible D. v. virgifera larvae to Bt maize roots. This study provides novel insights into insect response to Bt intoxication, and a possible framework for future investigations of resistance mechanisms.

Influence of host plant, geography and pheromone strain on genomic differentiation in sympatric populations of Ostrinia nubilalis.

Patterns of mating for the European corn borer (Ostrinia nubilalis) moth depend in part on variation in sex-pheromone blend. The ratio of (E)-11- and (Z)-11-tetradecenyl acetate (E11- and Z11-14:OAc) in the pheromone blend that females produce and males respond to differs between strains of O. nubilalis. Populations also vary in female oviposition preference for and larval performance on maize (C4) and nonmaize (C3) host plants. The relative contributions of sexual and ecological trait variation to the genetic structure of O. nubilalis remains unknown. Host-plant use (13 C/14 C ratios) and genetic differentiation were estimated among sympatric E and Z pheromone strain O. nubilalis males collected in sex-pheromone baited traps at 12 locations in Pennsylvania and New York between 2007 and 2010. Among genotypes at 65 single nucleotide polymorphism marker loci, variance at a position in the pheromone gland fatty acyl-reductase (pgfar) gene at the locus responsible for determining female pheromone ratio (Pher) explained 64% of the total genetic differentiation between males attracted to different pheromones (male response, Resp), providing evidence of sexual inter-selection at these unlinked loci. Principal coordinate, Bayesian clustering, and distance-based redundancy analysis (dbRDA) demonstrate that host plant history or geography does not significantly contribute to population variation or differentiation among males. In contrast, these analyses indicate that pheromone response and pgfar-defined strain contribute significantly to population genetic differentiation. This study suggests that behavioural divergence probably plays a larger role in driving genetic variation compared to host plant-defined ecological adaptation.

Cytochrome P450s Cyp4p1 and Cyp4p2 associated with the DDT tolerance in the Drosophila melanogaster strain 91-R.

Cytochrome P450s are part of a super-gene family that has undergone gene duplication, divergence, over-expression and, in some cases, loss of function. One such case is the 91-R and 91-C strains of common origin, in Drosophila melanogaster, whereby 91-R (DDT resistant strain) overexpresses Cyp4p1 and Cyp4p2 and both genes are lost in 91-C (DDT susceptible strain). In this study, we used a comparative approach to demonstrate that transcription of Cyp4p1 and Cyp4p2 were constitutively up-regulated in the Drosophila melanogaster strain 91-R as compared to another DDT susceptible strain Canton-S which does not have a loss of function of these genes. Furthermore, significantly increased expression of Cyp4p1 and Cyp4p2 was induced in 91-R in response to sublethal DDT exposure, however, such induction did not occur in the DDT treated Canton-S. Additionally, fixed nucleotide variation within putative transcription factor binding sites of Cyp4p1 and Cyp4p2 promoters were observed between 91-R and Canton-S, however, their impact on transcription remains to be determined. Two GAL4/UAS transgenic strains with integrated heat shock-inducible Cyp4p1- or Cyp4p2-RNAi constructs within wild-type genetic backgrounds were developed. Following heat shock induction of Cyp4p1 and Cyp4p2 knockdown, these transgenic lines showed increased DDT mortality as compared to their corresponding non-heat shock controls. These results provide a functional link of Cyp4p1 and Cyp4p2 in conferring tolerance to DDT exposure.

The invasive MED/Q Bemisia tabaci genome: a tale of gene loss and gene gain.

BACKGROUND: Sweetpotato whitefly, Bemisia tabaci MED/Q and MEAM1/B, are two economically important invasive species that cause considerable damages to agriculture crops through direct feeding and indirect vectoring of plant pathogens. Recently, a draft genome of B. tabaci MED/Q has been assembled. In this study, we focus on the genomic comparison between MED/Q and MEAM1/B, with a special interest in MED/Q's genomic signatures that may contribute to the highly invasive nature of this emerging insect pest. RESULTS: The genomes of both species share similarity in syntenic blocks, but have significant divergence in the gene coding sequence. Expansion of cytochrome P450 monooxygenases and UDP glycosyltransferases in MED/Q and MEAM1/B genome is functionally validated for mediating insecticide resistance in MED/Q using in vivo RNAi. The amino acid biosynthesis pathways in MED/Q genome are partitioned among the host and endosymbiont genomes in a manner distinct from other hemipterans. Evidence of horizontal gene transfer to the host genome may explain their obligate relationship. Putative loss-of-function in the immune deficiency-signaling pathway due to the gene loss is a shared ancestral trait among hemipteran insects. CONCLUSIONS: The expansion of detoxification genes families, such as P450s, may contribute to the development of insecticide resistance traits and a broad host range in MED/Q and MEAM1/B, and facilitate species' invasions into intensively managed cropping systems. Numerical and compositional changes in multiple gene families (gene loss and gene gain) in the MED/Q genome sets a foundation for future hypothesis testing that will advance our understanding of adaptation, viral transmission, symbiosis, and plant-insect-pathogen tritrophic interactions.

Variation in Mitochondria-Derived Transcript Levels Associated With DDT Resistance in the 91-R Strain of Drosophila melanogaster (Diptera: Drosophilidae).

The organochloride insecticide dichlorodiphenyltrichloroethane (DDT) and its metabolites can increase cellular levels of reactive oxygen species (ROS), cause mitochondrial dysfunction, and induce apoptosis. The highly DDT-resistant Drosophila melanogaster Meigen 1830 (Drosophila) strain, 91-R, and its susceptible control, 91-C, were used to investigate functional and structural changes among mitochondrial-derived pathways. Resequencing of mitochondrial genomes (mitogenomes) detected no structural differences between 91-R and 91-C, whereas RNA-seq suggested the differential expression of 221 mitochondrial-associated genes. Reverse transcriptase-quantitative PCR validation of 33 candidates confirmed that transcripts for six genes (Cyp12d1-p, Cyp12a4, cyt-c-d, COX5BL, COX7AL, CG17140) were significantly upregulated and two genes (Dif, Rel) were significantly downregulated in 91-R. Among the upregulated genes, four genes are duplicated within the reference genome (cyt-c-d, CG17140, COX5BL, and COX7AL). The predicted functions of the differentially expressed genes, or known functions of closely related genes, suggest that 91-R utilizes existing ROS regulation pathways of the mitochondria to combat increased ROS levels from exposure to DDT. This study represents, to our knowledge, the initial investigation of mitochondrial genome sequence variants and functional adaptations in responses to intense DDT selection and provides insights into potential adaptations of ROS management associated with DDT selection in Drosophila.

Binding affinity of five PBPs to Ostrinia sex pheromones.

BACKGROUND: Pheromone binding proteins (PBPs) of male Lepidoptera function in chemical communication, mate attraction and recognition. Directional selection was previously predicted between PBP3 orthologs of Ostrinia furnacalis and Ostrinia nubilalis were interpreted as being involved in sexual isolation. RESULTS: In vitro assays show that recombinant male OfurPBP3 bound O. furnacalis sex pheromones, Z-12-tetradecenyl acetate (Z12-14:OAc) and E-12-tetradecenyl acetate (E12-14:OAc), as well as to ECB pheromones Z11- and E11-14:OAc. Recombinant OfurPBP4 and OfurPBP5 bound E11- and Z11-14:OAc with greater affinity compared to Z12- and E12-14:OAc, and OfurPBP4 incapable of binding with E12-14:OAc. In silico molecular docking predicted OfurPBP3 residues Phe12, Ile52, Leu94, Ile113 within a hydrophobic ligand-binding pocket and may participate in E12- and Z12-14:OAc binding. Independent site-directed mutagenesis experiments demonstrated that Ser12, Asn52, Arg94, and Asn113 residues variants caused an approximately 1.7- to 4.6-fold reduction in OfurPBP3 affinity for Z12- and E12-14:OAc, and a 2.7- to 8.4-fold decrease in affinity towards E11- and Z11-14:OAc. CONCLUSIONS: Five PBPs of O. furnacalis play important functions in Ostrinia pheromones binding. These four amino acids may play a role in binding of sex pheromone, but this study does not address questions regarding specific response between males of O. furnacalis and O. nubilalis. Additional studies are required determine the role, if any, PBPs play in the evolution of sex pheromone communication.

Introgression between divergent corn borer species in a region of sympatry: Implications on the evolution and adaptation of pest arthropods.

The Asian corn borer, Ostrinia furnacalis, and European corn borer, O. nubilalis (Lepidoptera: Crambidae), cause damage to cultivated maize in spatially distinct geographies and have evolved divergent hydrocarbons as the basis of sexual communication. The Yili area of Xinjiang Uyghur Autonomous Region in China represents the only known region where O. furnacalis has invaded a native O. nubilalis range, and these two corn borer species have made secondary contact. Genetic differentiation was estimated between Ostrinia larvae collected from maize plants at 11 locations in Xinjiang and genotyped using high-throughput SNP and microsatellite markers. Maternal lineages were assessed by direct sequencing of mitochondrial cytochrome c oxidase subunit I and II haplotypes, and a high degree of genotypic diversity was demonstrated between lineages based on SNP genotypes. Furthermore, historical introgression was predicted among SNP genotypes only at sympatric locations in the Yili area, whereas in Xinjiang populations only O. furnacalis haplotypes were detected and no analogous introgressed genotypes were predicted. Our detection of putative hybrids and historical evidence of introgression defines Yili area as a hybrid zone between the species in normal ecological interactions and furthermore, might indicate that adaptive traits could spread even between seemingly divergent species through horizontal transmission. Results of this study indicate there may be a continuum in the degree of reproductive isolation between Ostrinia species and that the elegance of distinct and complete speciation based on modifications to the pheromone communication might need to be reconsidered.

A combination of sexual and ecological divergence contributes to rearrangement spread during initial stages of speciation.

Chromosomal rearrangements between sympatric species often contain multiple loci contributing to assortative mating, local adaptation and hybrid sterility. When and how these associations arise during the process of speciation remains a subject of debate. Here, we address the relative roles of local adaptation and assortative mating on the dynamics of rearrangement evolution by studying how a rearrangement covaries with sexual and ecological trait divergence within a species. Previously, a chromosomal rearrangement that suppresses recombination on the Z (sex) chromosome was identified in European corn borer moths (Ostrinia nubilalis). We further characterize this recombination suppressor and explore its association with variation in sex pheromone communication and seasonal ecological adaptation in pairs of populations that are divergent in one or both of these characteristics. Direct estimates of recombination suppression in pedigree mapping families indicated that more than 39% of the Z chromosome (encompassing up to ~10 megabases and ~300 genes) resides within a nonrecombining unit, including pheromone olfactory receptor genes and a major quantitative trait locus that contributes to ecotype differences (Pdd). Combining direct and indirect estimates of recombination suppression, we found that the rearrangement was occasionally present between sexually isolated strains (E vs. Z) and between divergent ecotypes (univoltine vs. bivoltine). However, it was only consistently present when populations differed in both sexual and ecological traits. Our results suggest that independent of the forces that drove the initial establishment of the rearrangement, a combination of sexual and ecological divergence is required for rearrangement spread during speciation.

Estimation of long terminal repeat element content in the Helicoverpa zea genome from high-throughput sequencing of bacterial artificial chromosome pools.

The lepidopteran pest insect Helicoverpa zea feeds on cultivated corn and cotton across the Americas where control remains challenging owing to the evolution of resistance to chemical and transgenic insecticidal toxins, yet genomic resources remain scarce for this species. A bacterial artificial chromosome (BAC) library having a mean genomic insert size of 145 ± 20 kbp was created from a laboratory strain of H. zea, which provides ∼12.9-fold coverage of a 362.8 ± 8.8 Mbp (0.37 ± 0.09 pg) flow cytometry estimated haploid genome size. Assembly of Illumina HiSeq 2000 reads generated from 14 pools that encompassed all BAC clones resulted in 165 485 genomic contigs (N50 = 3262 bp; 324.6 Mbp total). Long terminal repeat (LTR) protein coding regions annotated from 181 contigs included 30 Ty1/copia, 78 Ty3/gypsy, and 73 BEL/Pao elements, of which 60 (33.1%) encoded all five functional polyprotein (pol) domains. Approximately 14% of LTR elements are distributed non-randomly across pools of BAC clones.

Horizontal transfer of a non-autonomous Helitron among insect and viral genomes.

BACKGROUND: The movement of mobile elements among species by horizontal transposon transfer (HTT) influences the evolution of genomes through the modification of structure and function. Helitrons are a relatively new lineage of DNA-based (class II) transposable elements (TEs) that propagate by rolling-circle replication, and are capable of acquiring host DNA. The rapid spread of Helitrons among animal lineages by HTT is facilitated by shuttling in viral particles or by unknown mechanisms mediated by close organism associations (e.g. between hosts and parasites). RESULTS: A non-autonomous Helitron independently annotated as BmHel-2 from Bombyx mori and the MITE01 element from Ostrinia nubilalis was predicted in the genomes of 24 species in the insect Order Lepidoptera. Integrated Helitrons retained ≥ 65% sequence identity over a 250 bp consensus, and were predicted to retain secondary structures inclusive of a 3'-hairpin and a 5'-subterminal inverted repeat. Highly similar Hel-2 copies were predicted in the genomes of insects and associated viruses, which along with a previous documented case of real-time virus-insect cell line transposition suggests that this Helitron has likely propagated by HTT. CONCLUSIONS: These findings provide evidence that insect virus may mediate the HTT of Helitron-like TEs. This movement may facilitate the shuttling of DNA elements among insect genomes. Further sampling is required to determine the putative role of HTT in insect genome evolution.

Transcriptional analysis of susceptible and resistant European corn borer strains and their response to Cry1F protoxin.

BACKGROUND: Despite a number of recent reports of insect resistance to transgenic crops expressing insecticidal toxins from Bacillus thuringiensis (Bt), little is known about the mechanism of resistance to these toxins. The purpose of this study is to identify genes associated with the mechanism of Cry1F toxin resistance in European corn borer (Ostrinia nubilalis Hübner). For this, we compared the global transcriptomic response of laboratory selected resistant and susceptible O. nubilalis strain to Cry1F toxin. We further identified constitutive transcriptional differences between the two strains. RESULTS: An O. nubilalis midgut transcriptome of 36,125 transcripts was assembled de novo from 106 million Illumina HiSeq and Roche 454 reads and used as a reference for estimation of differential gene expression analysis. Evaluation of gene expression profiles of midgut tissues from the Cry1F susceptible and resistant strains after toxin exposure identified a suite of genes that responded to the toxin in the susceptible strain (n = 1,654), but almost 20-fold fewer in the resistant strain (n = 84). A total of 5,455 midgut transcripts showed significant constitutive expression differences between Cry1F susceptible and resistant strains. Transcripts coding for previously identified Cry toxin receptors, cadherin and alkaline phosphatase and proteases were also differentially expressed in the midgut of the susceptible and resistant strains. CONCLUSIONS: Our current study provides a valuable resource for further molecular characterization of Bt resistance and insect response to Cry1F toxin in O. nubilalis and other pest species.

The genetic structure of Asian corn borer, Ostrinia furnacalis, populations in China: haplotype variance in northern populations and potential impact on management of resistance to transgenic maize.

Asian corn borer, Ostrinia furnacalis (Guenée), is a severe pest that infests cultivated maize in the major production regions of China. Populations show genotype-by-environment variation in voltinism, such that populations with a single generation (univoltine) are fixed in Northern China where growing seasons are short. Low genetic differentiation was found among samples from 33 collection sites across China and one site from North Korea (n=1673) using variation at 6 nuclear microsatellite loci (ENA corrected global FST=0.020; P value<0.05). Analysis of molecular variance indicated that geographic region, number of generations or voltinism accounted for <0.38% of the total genetic variation at nuclear loci and was corroborated by clustering of co-ancestries among genotypes using the program STRUCTURE. In contrast, a mitochondrial haplotype network identified 4 distinct clusters, where 70.5% of samples from univoltine populations were within a single group. Univoltine populations were also placed into a unique cluster using Population Graph and Principal component analyses, which showed significant differentiation with multivoltine populations (φST=0.400; P value<0.01). This study suggests that gene flow among O. furnacalis in China may be high among regions, with the exception of northeastern localities. Haplotype variation may be due to random genetic drift resulting from partial reproductive isolation between univoltine and multivoltine O. furnacalis populations. Such reproductive isolation might impact the potential spread of alleles that confer resistance to transgenic maize in China.

Effect of maize lines on larval fitness costs of Cry1F resistance in the European corn borer (Lepidoptera: Crambidae).

Crops producing insecticidal toxins from the bacterium Bacillus thuringiensis (Bt) are widely planted and enable management of key insect pests while reducing the use of conventional insecticides. However, the evolution of Bt resistance could diminish these benefits. Fitness costs of Bt resistance occur in the absence of Bt toxin when individuals with resistance alleles show a reduction in fitness relative to susceptible individuals, and they can delay the evolution of resistance. Ecological factors including host-plant variety can affect the magnitude of fitness costs, and consequently, the degree to which fitness costs delay resistance. In this study, we measured fitness costs of resistance to Bt toxin Cry1F in the European corn borer Ostrinia nubilalis Hübner (Lepidoptera: Crambidae) using Cry1F-resistant and Cry1F-susceptible strains sharing a similar genetic background. Fitness costs were tested on three lines of maize, Zea mays L., by measuring larval survival and development in two greenhouse experiments with plants in either the vegetative or reproductive stage. Both experiments showed that maize line significantly affected larval survival and developmental rate. However, larval survival, mass, and developmental rate did not differ between the Cry1F-resistant and susceptible strains, indicating a lack of fitness costs of resistance to Cry1F for the larval fitness components measured in this experiment. Future experiments should test for fitness costs of Cry1F resistance affecting survival to adulthood and adult life-history parameters.

The Genome of Arsenophonus sp. and Its Potential Contribution in the Corn Planthopper, Peregrinus maidis.

The co-evolution between symbionts and their insect hosts has led to intricate functional interdependencies. Advances in DNA-sequencing technologies have not only reduced the cost of sequencing but, with the advent of highly accurate long-read methods, have also enabled facile genome assembly even using mixed genomic input, thereby allowing us to more easily assess the contribution of symbionts to their insect hosts. In this study, genomic data recently generated from Peregrinus maidis was used to assemble the genome of a bacterial symbiont, Pm Arsenophonus sp. This ~4.9-Mb assembly is one of the largest Arsenophonus genomes reported to date. The Benchmarking Universal Single-Copy Orthologs (BUSCO) result indicates that this Pm Arsenophonus assembly has a high degree of completeness, with 96% of the single-copy Enterobacterales orthologs found. The identity of the Pm Arsenophonus sp. was further confirmed by phylogenetic analysis. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis indicates a major contribution by Pm Arsenophonus sp. to the biosynthesis of B vitamins and essential amino acids in P. maidis, where threonine and lysine production is carried out solely by Pm Arsenophonus sp. This study not only provides deeper insights into the evolutionary relationships between symbionts and their insect hosts, but also adds to our understanding of insect biology, potentially guiding the development of novel pest control methods.

Evaluation of tolerance to Bacillus thuringiensis toxins among laboratory-reared western bean cutworm (Lepidoptera: Noctuidae).

The western bean cutworm, Striacosta albicosta (Smith) (Lepidoptera: Noctuidae), is a destructive insect pest of dry beans and corn within its native range of western Nebraska and eastern Colorado. However, since the initiation of an eastward range expansion of S. albicosta in the late 1990s, economic damage has been observed in the Midwest, and the species has now reached the Atlantic Coast and Quebec. Economic damage to corn occurs by larval feeding on ears, which is not controlled by commercial transgenic hybrids that express Bacillus thuringiensis (Bt) Cry1Ab, but partial control is observed by corn varieties that express Cry1 F toxins. Inadequate protocols for laboratory rearing of S. albicosta have hindered controlled efficacy experimentation in the laboratory and field. We report an S. albicosta rearing methodology used to maintain alaboratory colony for 12 continuous generations. Rearing procedures were adapted for Bt toxin diet-overlay assays, revealing that S. albicosta larvae exposed to Bt toxin for 14 d were insensitive to Cry1Ab concentrations up to 25,000 ng/cm2. In contrast, neonates exposed to Cry1 F toxin at > or = 250 ng/cm2, showed reduced developmental rates, with estimated effective concentrations of EC50 = 1,187.5 ng/cm2 and EC95 = 10,000.5 ng/cm2. The ability to mass produce this pest insect will enhance fundamental research, including evaluation of control tactics and toxin susceptibility.

A Chromosome-Scale Genome Assembly of a Helicoverpa zea Strain Resistant to Bacillus thuringiensis Cry1Ac Insecticidal Protein.

Helicoverpa zea (Lepidoptera: Noctuidae) is an insect pest of major cultivated crops in North and South America. The species has adapted to different host plants and developed resistance to several insecticidal agents, including Bacillus thuringiensis (Bt) insecticidal proteins in transgenic cotton and maize. Helicoverpa zea populations persist year-round in tropical and subtropical regions, but seasonal migrations into temperate zones increase the geographic range of associated crop damage. To better understand the genetic basis of these physiological and ecological characteristics, we generated a high-quality chromosome-level assembly for a single H. zea male from Bt-resistant strain, HzStark_Cry1AcR. Hi-C data were used to scaffold an initial 375.2 Mb contig assembly into 30 autosomes and the Z sex chromosome (scaffold N50 = 12.8 Mb and L50 = 14). The scaffolded assembly was error-corrected with a novel pipeline, polishCLR. The mitochondrial genome was assembled through an improved pipeline and annotated. Assessment of this genome assembly indicated 98.8% of the Lepidopteran Benchmark Universal Single-Copy Ortholog set were complete (98.5% as complete single copy). Repetitive elements comprised approximately 29.5% of the assembly with the plurality (11.2%) classified as retroelements. This chromosome-scale reference assembly for H. zea, ilHelZeax1.1, will facilitate future research to evaluate and enhance sustainable crop production practices.

Extended Sentinel Monitoring of Helicoverpa zea Resistance to Cry and Vip3Aa Toxins in Bt Sweet Corn: Assessing Changes in Phenotypic and Allele Frequencies of Resistance.

Transgenic corn and cotton that produce Cry and Vip3Aa toxins derived from Bacillus thuringiensis (Bt) are widely planted in the United States to control lepidopteran pests. The sustainability of these Bt crops is threatened because the corn earworm/bollworm, Helicoverpa zea (Boddie), is evolving a resistance to these toxins. Using Bt sweet corn as a sentinel plant to monitor the evolution of resistance, collaborators established 146 trials in twenty-five states and five Canadian provinces during 2020-2022. The study evaluated overall changes in the phenotypic frequency of resistance (the ratio of larval densities in Bt ears relative to densities in non-Bt ears) in H. zea populations and the range of resistance allele frequencies for Cry1Ab and Vip3Aa. The results revealed a widespread resistance to Cry1Ab, Cry2Ab2, and Cry1A.105 Cry toxins, with higher numbers of larvae surviving in Bt ears than in non-Bt ears at many trial locations. Depending on assumptions about the inheritance of resistance, allele frequencies for Cry1Ab ranged from 0.465 (dominant resistance) to 0.995 (recessive resistance). Although Vip3Aa provided high control efficacy against H. zea, the results show a notable increase in ear damage and a number of surviving older larvae, particularly at southern locations. Assuming recessive resistance, the estimated resistance allele frequencies for Vip3Aa ranged from 0.115 in the Gulf states to 0.032 at more northern locations. These findings indicate that better resistance management practices are urgently needed to sustain efficacy the of corn and cotton that produce Vip3Aa.

Distribution of genes and repetitive elements in the Diabrotica virgifera virgifera genome estimated using BAC sequencing.

Feeding damage caused by the western corn rootworm, Diabrotica virgifera virgifera, is destructive to corn plants in North America and Europe where control remains challenging due to evolution of resistance to chemical and transgenic toxins. A BAC library, DvvBAC1, containing 109,486 clones with 104 ± 34.5 kb inserts was created, which has an ~4.56X genome coverage based upon a 2.58 Gb (2.80 pg) flow cytometry-estimated haploid genome size. Paired end sequencing of 1037 BAC inserts produced 1.17 Mb of data (~0.05% genome coverage) and indicated ~9.4 and 16.0% of reads encode, respectively, endogenous genes and transposable elements (TEs). Sequencing genes within BAC full inserts demonstrated that TE densities are high within intergenic and intron regions and contribute to the increased gene size. Comparison of homologous genome regions cloned within different BAC clones indicated that TE movement may cause haplotype variation within the inbred strain. The data presented here indicate that the D. virgifera virgifera genome is large in size and contains a high proportion of repetitive sequence. These BAC sequencing methods that are applicable for characterization of genomes prior to sequencing may likely be valuable resources for genome annotation as well as scaffolding.