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.

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.

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.

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.

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.

The oriental armyworm genome yields insights into the long-distance migration of noctuid moths.

The oriental armyworm, Mythimna separata, is known for its long-distance seasonal migration and environment-dependent phase polymorphisms. Here, we present a chromosome-level genome reference and integrate multi-omics, functional genetics, and behavioral assays to explore the genetic bases of the hallmark traits of M. separata migration. Gene family comparisons show expansion of gustatory receptor genes in this cereal crop pest. Functional investigation of magnetoreception-related genes and associated flight behaviors suggest that M. separata may use the geomagnetic field to guide orientation in its nocturnal flight. Comparative transcriptome characterizes a suite of genes that may confer the observed plasticity between phases, including genes involved in protein processing, hormone regulation, and dopamine metabolism. We further report molecular signatures that underlie the dynamic regulation of a migratory syndrome coordinating reproduction and flight. Our study yields insights into environment-dependent developmental plasticity in moths and advances our understanding of long-distance migration in nocturnal insect pests.

Association of voltage-gated sodium channel mutations with field-evolved pyrethroid resistant phenotypes in soybean aphid and genetic markers for their detection.

The frequent use of insecticides to manage soybean aphids, Aphis glycines (Hemiptera: Aphididae), in the United States has contributed to field-evolved resistance. Pyrethroid-resistant aphids have nonsynonymous mutations in the voltage-gated sodium channel (vgsc). We identified a leucine to phenylalanine mutation at position 1014 (L1014F) and a methionine to isoleucine mutation (M918I) of the A. glycines vgsc, both suspected of conferring knockdown resistance (kdr) to lambda-cyhalothrin. We developed molecular markers to identify these mutations in insecticide-resistant aphids. We determined that A. glycines which survived exposure to a diagnostic concentration of lambda-cyhalothrin and bifenthrin via glass-vial bioassays had these mutations, and showed significant changes in the resistance allele frequency between samples collected before and after field application of lambda-cyhalothrin. Thus, a strong association was revealed between aphids with L1014F and M918I vgsc mutations and survival following exposure to pyrethroids. Specifically, the highest survival was observed for aphids with the kdr (L1014F) and heterozygote super-kdr (L1014F + M918I) genotypes following laboratory bioassays and in-field application of lambda-cyhalothrin. These genetic markers could be used as a diagnostic tool for detecting insecticide-resistant A. glycines and monitoring the geographic distribution of pyrethroid resistance. We discuss how generating these types of data could improve our efforts to mitigate the effects of pyrethroid resistance on crop production.

Genome size evolution in the beetle genus Diabrotica.

Diabrocite corn rootworms are one of the most economically significant pests of maize in the United States and Europe and an emerging model for insect-plant interactions. Genome sizes of several species in the genus Diabrotica were estimated using flow cytometry along with that of Acalymma vittatum as an outgroup. Genome sizes ranged between 1.56 and 1.64 gigabase pairs and between 2.26 and 2.59 Gb, respectively, for the Diabrotica subgroups fucata and virgifera; the Acalymma vittatum genome size was around 1.65 Gb. This result indicated that a substantial increase in genome size occurred in the ancestor of the virgifera group. Further analysis of the fucata group and the virgifera group genome sequencing reads indicated that the genome size difference between the Diabrotica subgroups could be attributed to a higher content of transposable elements, mostly miniature inverted-transposable elements and gypsy-like long terminal repeat retroelements.

Consequences of coupled barriers to gene flow for the build-up of genomic differentiation.

Theory predicts that when different barriers to gene flow become coincident, their joint effects enhance reproductive isolation and genomic divergence beyond their individual effects, but empirical tests of this "coupling" hypothesis are rare. Here, we analyze patterns of gene exchange among populations of European corn borer moths that vary in the number of acting barriers, allowing for comparisons of genomic variation when barrier traits or loci are in coincident or independent states. We find that divergence is mainly restricted to barrier loci when populations differ by a single barrier, whereas the coincidence of temporal and behavioral barriers is associated with divergence of two chromosomes harboring barrier loci. Furthermore, differentiation at temporal barrier loci increases in the presence of behavioral divergence and differentiation at behavioral barrier loci increases in the presence of temporal divergence. Our results demonstrate how the joint action of coincident barrier effects leads to levels of genomic differentiation that far exceed those of single barriers acting alone, consistent with theory arguing that coupling allows indirect selection to combine with direct selection and thereby lead to a stronger overall barrier to gene flow. Thus, the state of barriers-independent or coupled-strongly influences the accumulation of genomic differentiation.

Sequences Encoding a Novel Toursvirus Identified from Southern and Northern Corn Rootworms (Coleoptera: Chrysomelidae).

Sequences derived from a novel toursvirus were identified from pooled genomic short read data from U.S. populations of southern corn rootworm (SCR, Diabrotica undecimpunctata howardi Barber) and northern corn rootworm (NCR, Diabrotica barberi Smith & Lawrence). Most viral sequences were identified from the SCR genomic dataset. As proteins encoded by toursvirus sequences from SCR and NCR were almost identical, the contig sets from SCR and NCR were combined to generate 26 contigs. A total of 108,176 bp were assembled from these contigs, with 120 putative toursviral ORFs identified indicating that most of the viral genome had been recovered. These ORFs included all 40 genes that are common to members of the Ascoviridae. Two genes typically present in Ascoviridae (ATP binding cassette transport system permeases and Baculovirus repeated open reading frame), were not detected. There was evidence for transposon insertion in viral sequences at different sites in the two host species. Phylogenetic analyses based on a concatenated set of 45 translated protein sequences clustered toursviruses into a distinct clade. Based on the combined evidence, we propose taxonomic separation of toursviruses from Ascoviridae.

Evidence of enhanced reproductive performance and lack-of-fitness costs among soybean aphids, Aphis glycines, with varying levels of pyrethroid resistance.

BACKGROUND: Foliar application of insecticides is the main strategy to manage soybean aphid, Aphis glycines (Hemiptera: Aphididae), in the northcentral United States. Subpopulations of A. glycines have multiple nonsynonymous mutations in the voltage-gated sodium channel (vgsc) genes that are associated with pyrethroid resistance. We explored if fitness costs are associated with phenotypes conferred by vgsc mutations using life table analyses. We predicted that there would be significant differences between pyrethroid susceptibility and field-collected, parthenogenetic isofemale clones with differing, nonsynonymous mutations in vgsc genes. RESULTS: Estimated resistance ratios for the pyrethroid-resistant clones ranged from 3.1 to 37.58 and 5.6 to 53.91 for lambda-cyhalothrin and bifenthrin, respectively. Although life table analyses revealed some biological and demographic parameters to be significantly different among the clonal lines, there was no association between levels of pyrethroid resistance and a decline in fitness. By contrast, one of the most resistant clonal lines (SBA-MN1-2017) had a significantly higher finite rate of increase, intrinsic rate of increase and greater overall fitness compared to the susceptible control and other pyrethroid-resistant clonal lines. CONCLUSIONS: Our life history analysis suggests that there are no negative pleotropic effects associated with the pyrethroid resistance in the clonal A. glycines lines used in this study. We discuss the potential impact of these results on efficacies of insecticide resistance management (IRM) and integrated pest management (IPM) plans directed at delaying the spread of pyrethroid-resistant A. glycines.

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.

The USDA-ARS Ag100Pest Initiative: High-Quality Genome Assemblies for Agricultural Pest Arthropod Research.

The phylum Arthropoda includes species crucial for ecosystem stability, soil health, crop production, and others that present obstacles to crop and animal agriculture. The United States Department of Agriculture's Agricultural Research Service initiated the Ag100Pest Initiative to generate reference genome assemblies of arthropods that are (or may become) pests to agricultural production and global food security. We describe the project goals, process, status, and future. The first three years of the project were focused on species selection, specimen collection, and the construction of lab and bioinformatics pipelines for the efficient production of assemblies at scale. Contig-level assemblies of 47 species are presented, all of which were generated from single specimens. Lessons learned and optimizations leading to the current pipeline are discussed. The project name implies a target of 100 species, but the efficiencies gained during the project have supported an expansion of the original goal and a total of 158 species are currently in the pipeline. We anticipate that the processes described in the paper will help other arthropod research groups or other consortia considering genome assembly at scale.

bric à brac controls sex pheromone choice by male European corn borer moths.

The sex pheromone system of ~160,000 moth species acts as a powerful form of assortative mating whereby females attract conspecific males with a species-specific blend of volatile compounds. Understanding how female pheromone production and male preference coevolve to produce this diversity requires knowledge of the genes underlying change in both traits. In the European corn borer moth, pheromone blend variation is controlled by two alleles of an autosomal fatty-acyl reductase gene expressed in the female pheromone gland (pgFAR). Here we show that asymmetric male preference is controlled by cis-acting variation in a sex-linked transcription factor expressed in the developing male antenna, bric à brac (bab). A genome-wide association study of preference using pheromone-trapped males implicates variation in the 293 kb bab intron 1, rather than the coding sequence. Linkage disequilibrium between bab intron 1 and pgFAR further validates bab as the preference locus, and demonstrates that the two genes interact to contribute to assortative mating. Thus, lack of physical linkage is not a constraint for coevolutionary divergence of female pheromone production and male behavioral response genes, in contrast to what is often predicted by evolutionary theory.

Nudivirus Sequences Identified from the Southern and Western Corn Rootworms (Coleoptera: Chrysomelidae).

Analysis of pooled genomic short read sequence data revealed the presence of nudivirus-derived sequences from U.S. populations of both southern corn rootworm (SCR, Diabrotica undecimpunctata howardi Barber) and western corn rootworm (WCR, Diabrotica virgifera virgifera LeConte). A near complete nudivirus genome sequence was assembled from sequence data for an SCR population with relatively high viral titers. A total of 147,179 bp was assembled from five contigs that collectively encode 109 putative open reading frames (ORFs) including 20 nudivirus core genes. In contrast, genome sequence recovery was incomplete for a second nudivirus from WCR, although sequences derived from this virus were present in three geographically dispersed populations. Only 48,989 bp were assembled with 48 putative ORFs including 13 core genes, representing about 20% of a typical nudivirus genome. Phylogenetic analysis indicated that both corn rootworm nudiviruses grouped with the third known nudivirus of beetles, Oryctes rhinoceros nudivirus in the genus Alphanudivirus. On the basis of phylogenetic and additional analyses, we propose further taxonomic separation of nudiviruses within Alphanudivirus and Betanudivirus into two subfamilies and five genera. Identification of nudivirus-derived sequences from two species of corn rootworm highlights the diversity of viruses associated with these agricultural insect pests.

Prediction of mitochondrial genome-wide variation through sequencing of mitochondrion-enriched extracts.

Although mitochondrial DNA (mtDNA) haplotype variation is often applied for estimating population dynamics and phylogenetic relationships, economical and generalized methods for entire mtDNA genome enrichment prior to high-throughput sequencing are not readily available. This study demonstrates the utility of differential centrifugation to enrich for mitochondrion within cell extracts prior to DNA extraction, short-read sequencing, and assembly using exemplars from eight maternal lineages of the insect species, Ostrinia nubilalis. Compared to controls, enriched extracts showed a significant mean increase of 48.2- and 86.1-fold in mtDNA based on quantitative PCR, and proportion of subsequent short sequence reads that aligned to the O. nubilalis reference mitochondrial genome, respectively. Compared to the reference genome, our de novo assembled O. nubilalis mitochondrial genomes contained 82 intraspecific substitution and insertion/deletion mutations, and provided evidence for correction of mis-annotated 28 C-terminal residues within the NADH dehydrogenase subunit 4. Comparison to a more recent O. nubilalis mtDNA assembly from unenriched short-read data analogously showed 77 variant sites. Twenty-eight variant positions, and a triplet ATT codon (Ile) insertion within ATP synthase subunit 8, were unique within our assemblies. This study provides a generalizable pipeline for whole mitochondrial genome sequence acquisition adaptable to applications across a range of taxa.

Post-transcriptional modulation of cytochrome P450s, Cyp6g1 and Cyp6g2, by miR-310s cluster is associated with DDT-resistant Drosophila melanogaster strain 91-R.

The role of miRNAs in mediating insecticide resistance remains largely unknown, even for the model species Drosophila melanogaster. Building on prior research, this study used microinjection of synthetic miR-310s mimics into DDT-resistant 91-R flies and observed both a significant transcriptional repression of computationally-predicted endogenous target P450 detoxification genes, Cyp6g1 and Cyp6g2, and also a concomitant increase in DDT susceptibility. Additionally, co-transfection of D. melanogaster S2 cells with dual luciferase reporter constructs validated predictions that miR-310s bind to target binding sites in the 3' untranslated regions (3'-UTR) of both Cyp6g1 and Cyp6g2 in vitro. Findings in the current study provide empirical evidence for a link between reduced miRNA expression and an insecticidal resistance phenotype through reduced targeted post-transcriptional suppression of transcripts encoding proteins involved in xenobiotic detoxification. These insights are important for understanding the breadth of adaptive molecular changes that have contributed to the evolution of DDT resistance in D. melanogaster.

Geographic Distribution of Bacillus thuringiensis Cry1F Toxin Resistance in Western Bean Cutworm (Lepidoptera: Noctuidae) Populations in the United States.

The western bean cutworm (WBC), Striacosta albicosta (Lepidoptera: Noctuidae), can be a severe pest of transgenic corn in the western Plains and Great Lakes regions of North America, including on hybrids expressing the Bacillus thuringiensis (Bt) Cry1F toxin. The level and geographic distribution of Cry1F resistance are not completely known. Neonate S. albicosta from 10 locations between Nebraska and New York state were subjected to dose-response trypsin-activated native Cry1F toxin overlay bioassays. In 2017, the mean estimated lethal concentration causing 50% larval mortality (LC50) ranged from 15.1 to 18.4 µg Cry1F cm-2, and were not significantly different among locations. In 2018, LC50 estimates at Scottsbluff, NE (22.0 µg Cry1F cm-2) and Watertown, NY (21.7 µg Cry1F cm-2) were significantly higher when compared to locations in Michigan (15.8 µg Cry1F cm-2). Significantly lower 14-day larval weight among survivors was correlated with higher Cry1F dose. Results from this study indicate that S. albicosta survivorship on purified Bt Cry1F toxin shows a relatively even distribution across the native and range expansion areas where seasonal field infestations typically occur.

Endogenous viral elements integrated into the genome of the soybean aphid, Aphis glycines.

Sequence analysis of the genomic DNA isolated from four biotypes of the soybean aphid, Aphis glycines (AG), revealed that in addition to the commonly observed retrovirus-related retrotransposons, viral sequences derived from multiple RNA and DNA viruses have integrated into the genome. Notably, sequences of more than 60 nudiviral genes were identified from de novo assembled DNA contigs, and mapped to assembled genomic scaffolds of AG, indicating that an ancient nudivirus, named Aphis glycines endogenous nudivirus (AgENV), had integrated into the AG genome. Furthermore, sequences derived from a similar endogenous nudivirus, Melanaphis sacchari endogenous nudivirus (MsENV), were identified from the genomic scaffolds of the sugarcane aphid, Melanaphis sacchari. Analysis of transcriptome and small RNA sequence data derived from AG did not provide evidence for transcription of the integrated AgENV genes. Hence, the genes of AgENV may be present as pseudogenes. Phylogenetic analysis based on nudivirus core genes indicated that these aphid ENVs belong to the genus Alphanudivirus.

Genome scan detection of selective sweeps among biotypes of the soybean aphid, Aphis glycines, with differing virulence to resistance to A. glycines (Rag) traits in soybean, Glycine max.

Multiple biotypes of soybean aphid, Aphis glycines, occur in North America adapted for survival (virulence) on soybean, Glycine max, with one or more different resistance to A. glycines (Rag) traits. The degree of genome-wide variance between biotypes and the basis of virulence remains unknown, but the latter is hypothesized to involve secreted effector proteins. Between 167,249 and 217,750 single nucleotide polymorphisms (SNPs) were predicted from whole genome re-sequencing of A. glycines avirulent biotype 1 (B1) and virulent B2, B3 and B4 colony-derived iso-female lines when compared to the draft B1 genome assembly, Ag_bt1_v6.0. Differences in nucleotide diversity indices (π) estimated within 1000 bp sliding windows demonstrated that 226 of 353 (64.0%) regions most differentiated between B1 and ≥ 2 virulent biotypes, representing < 0.1% of the 308 Mb assembled genome size, are located on 15 unordered scaffolds. Furthermore, these 226 intervals were coincident and show a significant association with 326 of 508 SNPs with significant locus-by-locus FST estimates between biotype populations (r = 0.6271; F1,70 = 45.36, P < 0.001) and genes showing evidence of directions selection (πN/πS > 2.0; r = 0.6233; F1,70 = 50.20, P < 0.001). A putative secreted effector glycoprotein is encoded in proximity to genome intervals of low estimated π (putative selective sweep) within avirulent B1 compared to all three virulent biotypes. Additionally, SNPs are clustered in or in proximity to genes putatively involved in intracellular protein cargo transport and the regulation of secretion. Results of this study indicate that factors on a small number of scaffolds of the A. glycines genome may contribute to variance in virulence towards Rag traits in G. max.