Optimizing a rapid LAMP assay for discrimination of Drosophila suzukii (Diptera: Drosophilidae) from common drosophilids captured in monitoring traps from the Midwest, United States

Spotted-wing Drosophila, Drosophila suzukii, is an economically important pest of small fruits worldwide. Currently, the timing of management strategies relies on detection of adult flies captured in baited monitoring traps; however, identifying D. suzukii in trap catch based on morphology can be challenging for growers. DNA-based diagnostic methods such as loop-mediated isothermal amplification (LAMP) have the potential to improve D. suzukii detection. This study evaluated a LAMP assay as a diagnostic tool to discriminate between D. suzukii and closely related drosophilid species found commonly in monitoring traps in the Midwestern United States. Targeting the mitochondrial cytochrome oxidase I (COI) gene, we found the LAMP assay accurately detected D. suzukii with as little as 0.1 ng/μl of DNA at 63 °C for 50 min. Under these optimal incubation conditions, D. suzukii could be discriminated from D. affinis and D. simulans consistently, when specimens collected from liquid monitoring traps were tested independently. Compared to other DNA-based diagnostic tools for D. suzukii, LAMP offers unique benefits: DNA extraction is not required, testing occurs at one temperature in less than 1 h, and positive results are visible as a colorimetric change from pink to yellow. The LAMP assay for D. suzukii can reduce reliance on morphological identification, enhance the adoption of monitoring tools, and improve accuracy of detection. Further optimization can be conducted to evaluate the accuracy and sensitivity of results when a mixture of DNA from both D. suzukii and congener flies are tested in a single LAMP reaction.

Genetic mutations linked to field-evolved Cry1Fa-resistance in the European corn borer, Ostrinia nubilalis

Transgenic corn, Zea mays (L.), expressing insecticidal toxins such as Cry1Fa, from Bacillus thuringiensis (Bt corn) targeting Ostrinia nubilalis (Hübner) (Lepidoptera: Crambidae) resulted in over 20 years of management success. The first case of practical field-evolved resistance by O. nubilalis to a Bt corn toxin, Cry1Fa, was discovered in Nova Scotia, Canada, in 2018. Laboratory-derived Cry1Fa-resistance by O. nubilalis was linked to a genome region encoding the ATP Binding Cassette subfamily C2 (ABCC2) gene; however, the involvement of ABCC2 and specific mutations in the gene leading to resistance remain unknown. Using a classical candidate gene approach, we report on O. nubilalis ABCC2 gene mutations linked to laboratory-derived and field-evolved Cry1Fa-resistance. Using these mutations, a DNA-based genotyping assay was developed to test for the presence of the Cry1Fa-resistance alleles in O. nubilalis strains collected in Canada. Screening data provide strong evidence that field-evolved Cry1Fa-resistance in O. nubilalis maps to the ABCC2 gene and demonstrates the utility of this assay for detecting the Cry1Fa resistance allele in O. nubilalis. This study is the first to describe mutations linked to Bt resistance in O. nubilalis and provides a DNA-based detection method that can be used for monitoring.

Presence-Absence Sampling Plans for Stink Bugs (Hemiptera: Pentatomidae) in the Midwest Region of the United States

Stink bugs represent an increasing risk to soybean production in the Midwest region of the United States. The current sampling protocol for stink bugs in this region is tailored for population density estimation and thus is more relevant to research purposes. A practical decision-making framework with more efficient sampling effort for management of herbivorous stink bugs is needed. Therefore, a binomial sequential sampling plan was developed for herbivorous stink bugs in the Midwest region. A total of 146 soybean fields were sampled across 11 states using sweep nets in 2016, 2017, and 2018. The binomial sequential sampling plans were developed using combinations of five tally thresholds at two proportion infested action thresholds to identify those that provided the best sampling outcomes. Final assessment of the operating characteristic curves for each plan indicated that a tally threshold of 3 stink bugs per 25 sweeps, and proportion infested action thresholds of 0.75 and 0.95 corresponding to the action thresholds of 5 and 10 stink bugs per 25 sweeps, provided the optimal balance between highest probability of correct decisions (≥ 99%) and lowest probability of incorrect decisions (≤ 1%). In addition, the average sample size for both plans (18 and 12 sets of 25 sweeps, respectively) was lower than that for the other proposed plans. The binomial sequential sampling plan can reduce the number of sample units required to achieve a management decision, which is important because it can potentially reduce risk/cost of management for stink bugs in soybean in this region.

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.

Hybridization and introgression between Helicoverpa armigera and H. zea: an adaptational bridge

BACKGROUND

Invasion of organisms into new ecosystems is increasingly common, due to the global trade in commodities. One of the most complex post-invasion scenarios occurs when an invasive species is related to a native pest, and even more so when they can hybridize and produce fertile progeny. The global pest Helicoverpa armigera was first detected in Brazil in 2013 and generated a wave of speculations about the possibility of hybridization with the native sister taxon Helicoverpa zea. In the present study, we used genome-wide single nucleotide polymorphisms from field-collected individuals to estimate hybridization between H. armigera and H. zea in different Brazilian agricultural landscapes.

RESULTS

The frequency of hybridization varied from 15 to 30% depending on the statistical analyses. These methods showed more congruence in estimating that hybrids contained approximately 10% mixed ancestry (i.e. introgression) from either species. Hybridization also varied considerably depending on the geographic locations where the sample was collected, forming a 'mosaic' hybrid zone where introgression may be facilitated by environmental and landscape variables. Both landscape composition and bioclimatic variables indicated that maize and soybean cropland are the main factors responsible for high levels of introgression in agricultural landscapes. The impact of multiple H. armigera incursions is reflected in the structured and inbred pattern of genetic diversity.

CONCLUSIONS

Our data showed that the landscape composition and bioclimatic variables influence the introgression rate between H. armigera and H. zea in agricultural areas. Continuous monitoring of the hybridization process in the field is necessary, since agricultural expansion, climatic fluctuations, changing composition of crop species and varieties, and dynamic planting seasons are some factors in South America that could cause a sudden alteration in the introgression rate between Helicoverpa species. Introgression between invasive and native pests can dramatically impact the evolution of host ranges and resistance management.

Spatial Patterns and Sequential Sampling Plans for Estimating Densities of Stink Bugs (Hemiptera: Pentatomidae) in Soybean in the North Central Region of the United States

Stink bugs are an emerging threat to soybean (Fabales: Fabaceae) in the North Central Region of the United States. Consequently, region-specific scouting recommendations for stink bugs are needed. The aim of this study was to characterize the spatial pattern and to develop sampling plans to estimate stink bug population density in soybean fields. In 2016 and 2017, 125 fields distributed across nine states were sampled using sweep nets. Regression analyses were used to determine the effects of stink bug species [Chinavia hilaris (Say) (Hemiptera: Pentatomidae) and Euschistus spp. (Hemiptera: Pentatomidae)], life stages (nymphs and adults), and field locations (edge and interior) on spatial pattern as represented by variance-mean relationships. Results showed that stink bugs were aggregated. Sequential sampling plans were developed for each combination of species, life stage, and location and for all the data combined. Results for required sample size showed that an average of 40-42 sample units (sets of 25 sweeps) would be necessary to achieve a precision of 0.25 for stink bug densities commonly encountered across the region. However, based on the observed geographic gradient of stink bug densities, more practical sample sizes (5-10 sample units) may be sufficient in states in the southeastern part of the region, whereas impractical sample sizes (>100 sample units) may be required in the northwestern part of the region. Our findings provide research-based sampling recommendations for estimating densities of these emerging pests in soybean.

Community Composition, Abundance, and Phenology of Stink Bugs (Hemiptera: Pentatomidae) in Soybean in the North Central Region of the United States

Stink bugs (Hemiptera: Pentatomidae) are an increasing threat to soybean (Fabales: Fabaceae) production in the North Central Region of the United States, which accounts for 80% of the country's total soybean production. Characterization of the stink bug community is essential for development of management programs for these pests. However, the composition of the stink bug community in the region is not well defined. This study aimed to address this gap with a 2-yr, 9-state survey. Specifically, we characterized the relative abundance, richness, and diversity of taxa in this community, and assessed phenological differences in abundance of herbivorous and predatory stink bugs. Overall, the stink bug community was dominated by Euschistus spp. (Hemiptera: Pentatomidae) and Chinavia hilaris (Say) (Hemiptera: Pentatomidae). Euschistus variolarius (Palisot de Beauvois) (Hemiptera: Pentatomidae), C. hilaris and Halyomorpha halys (Stål) (Hemiptera: Pentatomidae) were more abundant in the northwestern, southeastern and eastern parts, respectively, of the North Central Region of the United States. Economically significant infestations of herbivorous species occurred in fields in southern parts of the region. Species richness differed across states, while diversity was the same across the region. Herbivorous and predatory species were more abundant during later soybean growth stages. Our results represent the first regional characterization of the stink bug community in soybean fields and will be fundamental for the development of state- and region-specific management programs for these pests in the North Central Region of the United States.

Weekly Survivorship Curves of Soybean Aphid Biotypes 1 and 4 on Insecticidal Seed-Treated Soybean

Thiamethoxam, an insecticide used in soybean seed treatments, effectively suppresses soybean aphids (Aphis glycines) Matsumura (Hemiptera: Aphididae) for a short time after planting. However, exactly when and how quickly soybean aphid populations could increase is unknown. Likewise, we lack data on virulent soybean aphid biotypes (that can overcome soybean resistance) when fed on seed-treated soybean. Determining the survival of soybean aphids over time on insecticidal seed-treated soybean is critical for improving soybean aphid management and may provide insights to manage aphid virulence to aphid resistant-soybean. In greenhouse and field experiments, aphid-susceptible soybean plants (with and without an insecticidal seed treatment) were infested at 7, 14, 21, 28, 35, and 42 days after planting (DAP). We compared aphid survival among biotypes 1 (avirulent) and 4 (virulent) and insecticide treatment 72 h after infestation. We also measured thiamethoxam concentrations in plant tissue using liquid chromatography-tandem mass spectrometry. As expected, soybean aphid survival was significantly lower on seed-treated soybean up to 35 DAP for both biotypes, which correlates with the decrease of thiamethoxam in the plant over time. Moreover, we found no significant difference between avirulent and virulent biotype survivorship on insecticidal seed-treated soybean plants, although we did find significantly greater survival for the virulent biotype compared with the avirulent biotype on untreated soybean in the field. In conclusion, our study further characterized the relative short duration of seed treatment effectiveness on soybean aphid and showed that survivorship of virulent aphids on seed-treated soybean is similar to avirulent aphids.

Inward rectifier potassium (Kir) channels in the soybean aphid Aphis glycines: Functional characterization, pharmacology, and toxicology

Inward rectifier K+ (Kir) channels contribute to a variety of physiological processes in insects and are emerging targets for insecticide development. Previous studies on insect Kir channels have primarily focused on dipteran species (e.g., mosquitoes, fruit flies). Here we identify and functionally characterize Kir channel subunits in a hemipteran insect, the soybean aphid Aphis glycines, which is an economically important insect pest and vector of soybeans. From the transcriptome and genome of Ap. glycines we identified two cDNAs, ApKir1 and ApKir2, encoding Kir subunits that were orthologs of insect Kir1 and Kir2, respectively. Notably, a gene encoding a Kir3 subunit was absent from the transcriptome and genome of Ap. glycines, similar to the pea aphid Acyrthosiphon pisum. Heterologous expression of ApKir1 and ApKir2 in Xenopus laevis oocytes enhanced K+-currents in the plasma membrane; these currents were inhibited by barium and the small molecule VU041. Compared to ApKir2, ApKir1 mediated currents that were larger in magnitude, more sensitive to barium, and less inhibited by small molecule VU041. Moreover, ApKir1 exhibited stronger inward rectification compared to ApKir2. Topical application of VU041 in adult aphids resulted in dose-dependent mortality within 24 h that was more efficacious than flonicamid, an established insecticide also known to inhibit Kir channels. We conclude that despite the apparent loss of Kir3 genes in aphid evolution, Kir channels are important to aphid survival and represent a promising target for the development of new insecticides.

Influence of historical land use and modern agricultural expansion on the spatial and ecological divergence of sugarcane borer, Diatraea saccharalis (Lepidoptera: Crambidae) in Brazil

Human-mediated changes in landscapes can facilitate niche expansion and accelerate the adaptation of insect species. The interaction between the evolutionary history of the sugarcane borer, Diatraea saccharalis Fabricius, and historical and modern agricultural activity in Brazil shaped its spatial genetic structure, facilitating ecological divergence and incipient host shifting. Based on microsatellite data, STRUCTURE analyses identified two (K = 2) and three (K = 3) significant genetic clusters that corresponded to: (a) a strong signal of spatial genetic structure and, (b) a cryptic signal of host differentiation. We inferred that K = 2 reflects the footprint of agricultural activity, such as expansion of crop production (sugarcane and maize), unintentional dispersion of pests, and management practices. In contrast, K = 3 indicated incipient host differentiation between larvae collected from sugarcane or maize. Our estimates of population size changes indicated that a historical bottleneck was associated with a reduction of sugarcane production ≈200 years ago. However, a more recent population expansion was detected (>1950s), associated with agricultural expansion of large crop production into previously unfarmed land. Partial Mantel tests supported our hypothesis of incipient host adaptation, and identified isolation-by-environment (e.g., host plant) in São Paulo and Minas Gerais states, where sugarcane has been traditionally produced in Brazil. The impact of agricultural production on D. saccharalis may continue, as the current population structure may hinder the efficacy of refuge plants in delaying insect resistance evolution to Bt toxin.

Pan-American Similarities in Genetic Structures of Helicoverpa armigera and Helicoverpa zea (Lepidoptera: Noctuidae) With Implications for Hybridization

The genus Helicoverpa (Lepidoptera: Noctuidae) includes phytophagous and polyphagous agricultural insect pests. In the Americas, a native pest, Helicoverpa zea (Boddie), and an invasive pest, Helicoverpa armigera (Hübner), are causing severe damage in vegetable and agronomic crops. The population structure of both species in South America is poorly understood, and the phylogenetic relatedness of H. armigera and H. zea suggests natural interspecific gene flow between these species. Using microsatellite loci, we investigated: 1) the genetic diversity and gene flow of H. armigera specimens from Brazil; 2) the genetic diversity and gene flow between H. zea specimens from Brazil and the United States; and 3) the possibility of interspecific gene flow and the frequency of putative hybrids in Brazil. We detected high intraspecific gene flow among populations collected in the same country. However, there is a geographic limit to gene flow among H. zea individuals from South and North America. Pairwise Fst and private alleles showed that H. armigera is more similar to H. zea from Brazil than H. zea from the United States. A comparative STRUCTURE analysis suggests natural hybridization between H. armigera and H. zea in Brazil. High gene flow and natural hybridization are key traits to population adaptation in new and disturbed environments, which can influence the management of these pests in the American continent.

Transcriptomic dynamics in soybean near-isogenic lines differing in alleles for an aphid resistance gene, following infestation by soybean aphid biotype 2

BACKGROUND

Genetic resistance of soybean [Glycine max (L.) Merr] against Aphis glycines provides effective management of this invasive pest, though the underlying molecular mechanisms are largely unknown. This study aimed to investigate genome-wide changes in gene expressions of soybean near-isogenic lines (NILs) either with the Rag5 allele for resistance or the rag5 allele for susceptibility to the aphid following infestation with soybean aphid biotype 2.

RESULTS

The resistant (R)-NIL responded more rapidly to aphid infestation than the susceptible (S)-NIL, with differential expressions of 2496 genes during first 12 h of infestation (hai), compared to the aphid-free control. Although the majority of the differentially expressed genes (DEGs) in the R-NIL also responded to aphid infestation in S-NIL, overall the response time was longer and/or the magnitude of change was smaller in the S-NIL. In addition, 915 DEGs in R-NIL continued to be regulated at all time points (0, 6, 12, and 48 hai), while only 20 DEGs did so in S-NIL. Enriched gene ontology of the 2496 DEGs involved in plant defense responses including primary metabolite catalysis, oxidative stress reduction, and phytohormone-related signaling. By comparing R- vs. S-NIL, a total of 556 DEGs were identified. Of the 13 genes annotated in a 120-kb window of the Rag5 locus, two genes (Glyma.13 g190200 and Glyma.13 g190600) were differentially expressed (upregulated in S- or R-NIL), and another gene (Glyma.13 g190500) was induced up to 4-fold in the R-NIL at 6 and 12 h following aphid infestation.

CONCLUSIONS

This study strengthens our understanding of the defense dynamics in compatible and incompatible interactions of soybean and soybean aphid biotype 2. Several DEGs (e.g., Glyma.13 g190200, Glyma.13 g190500, and Glyma.13 g190600) near the Rag5 locus are strong candidate genes for further investigations.

Restriction site associated DNA (RAD) for de novo sequencing and marker discovery in sugarcane borer, Diatraea saccharalis Fab. (Lepidoptera: Crambidae)

We present the development of a genomic library using RADseq (restriction site associated DNA sequencing) protocol for marker discovery that can be applied on evolutionary studies of the sugarcane borer Diatraea saccharalis, an important South American insect pest. A RADtag protocol combined with Illumina paired-end sequencing allowed de novo discovery of 12 811 SNPs and a high-quality assembly of 122.8M paired-end reads from six individuals, representing 40 Gb of sequencing data. Approximately 1.7 Mb of the sugarcane borer genome distributed over 5289 minicontigs were obtained upon assembly of second reads from first reads RADtag loci where at least one SNP was discovered and genotyped. Minicontig lengths ranged from 200 to 611 bp and were used for functional annotation and microsatellite discovery. These markers will be used in future studies to understand gene flow and adaptation to host plants and control tactics.

Whole Genome Sequence of the Soybean Aphid Endosymbiont Buchnera aphidicola and Genetic Differentiation among Biotype-Specific Strains

Endosymbiosis with microorganisms is common in insects, with more than 10% of species requiring the metabolic capabilities of intracellular bacteria for their nutrient acquisition. Aphids harbor an obligate mutualism with the vertically transferred endosymbiont, Buchnera aphidicola, which produces key nutrients lacking in the aphid's phloem-based diet that are necessary for normal development and reproduction. It is thought that, in some groups of insects, bacterial symbionts may play key roles in biotype evolution against host-plant resistance. The genome of Buchnera has been sequenced in several aphid strains but little genomic data is currently available for the soybean aphid (Aphis glycines), one of the most important pests of soybean in North America. In this study, DNA sequencing was used to assemble and annotate the genome sequence of the Buchnera A. glycines strain and to reconstruct phylogenetic relationships among different strains. In addition, we identified several fixed Buchnera SNPs between Aphis glycines biotypes that were avirulent or virulent to a soybean aphid resistance gene (Rag1). The results of this study describe the genetic and evolutionary relationships of the Buchnera A. glycines strain, and begin to define the roles of an aphid symbiont in host-plant resistance.

Shifts in Buchnera aphidicola density in soybean aphids (Aphis glycines) feeding on virus-infected soybean

Vertically transmitted bacterial symbionts are common in arthropods. Aphids undergo an obligate symbiosis with Buchnera aphidicola, which provides essential amino acids to its host and contributes directly to nymph growth and reproduction. We previously found that newly adult Aphis glycines feeding on soybean infected with the beetle-transmitted Bean pod mottle virus (BPMV) had significantly reduced fecundity. We hypothesized that the reduced fecundity was attributable to detrimental impacts of the virus on the aphid microbiome, namely Buchnera. To test this, mRNA sequencing and quantitative real-time PCR were used to assay Buchnera transcript abundance and titre in A. glycines feeding on Soybean mosaic virus-infected, BPMV-infected, and healthy soybean for up to 14 days. Our results indicated that Buchnera density was lower and ultimately suppressed in aphids feeding on virus-infected soybean. While the decreased Buchnera titre may be associated with reduced aphid fecundity, additional mechanisms are probably involved. The present report begins to describe how interactions among insects, plants, and plant pathogens influence endosymbiont population dynamics.

Molecular interactions and immune responses between Maize fine streak virus and the leafhopper vector Graminella nigrifrons through differential expression and RNA interference

Graminella nigrifrons is the only known vector for Maize fine streak virus (MFSV). In this study, we used real-time quantitative PCR to compare the expression profiles of transcripts that putatively function in the insect immune response: four peptidoglycan recognition proteins (PGRP-SB1, -SD, -LC and LB), Toll, spaetzle, defensin, Dicer-2 (Dcr-2), Argonaut-2 (Ago-2) and Arsenic resistance protein 2 (Ars-2). Except for PGRP-LB and defensin, transcripts involved in humoral pathways were significantly suppressed in G. nigrifrons fed on MFSV-infected maize. The abundance of three RNA interference (RNAi) pathway transcripts (Dcr-2, Ago-2, Ars-2) was significantly lower in nontransmitting relative to transmitting G. nigrifrons. Injection with double-stranded RNA (dsRNA) encoding segments of the PGRP-LC and Dcr-2 transcripts effectively reduced transcript levels by 90 and 75% over 14 and 22 days, respectively. MFSV acquisition and transmission were not significantly affected by injection of either dsRNA. Knock-down of PGRP-LC resulted in significant mortality (greater than 90%) at 27 days postinjection, and resulted in more abnormal moults relative to those injected with Dcr-2 or control dsRNA. The use of RNAi to silence G. nigrifrons transcripts will facilitate the study of gene function and pathogen transmission, and may provide approaches for developing novel targets of RNAi-based pest control.

Identification and mapping of quantitative trait loci (QTL) conferring resistance to Fusarium graminearum from soybean PI 567301B

KEY MESSAGE

A major novel QTL was identified in a recombinant inbred line population derived from a cross of 'Wyandot' × PI 567301B for Fusarium graminearum, a seed and seedling pathogen of soybean. Fusarium graminearum is now recognized as a primary pathogen of soybean, causing root, seed rot and seedling damping-off in North America. In a preliminary screen, 'Wyandot' and PI 567301B were identified with medium and high levels of partial resistance to F. graminearum, respectively. The objective of this study was to characterise resistance towards F. graminearum using 184 recombinant inbred lines (RILs) derived from a cross of 'Wyandot' × PI 567301B. The parents and the RILs of the mapping population were evaluated for resistance towards F. graminearum using the rolled towel assay in a randomized incomplete block design. A genetic map was constructed from 2545 SNP markers and 2 SSR markers by composite interval mapping. One major and one minor QTL were identified on chromosomes 8 and 6, respectively, which explained 38.5 and 8.1 % of the phenotypic variance. The major QTL on chromosome 8 was mapped to a 300 kb size genomic region of the Williams 82 sequence. Annotation of this region indicates that there are 39 genes including the Rhg4 locus for soybean cyst nematode (SCN) resistance. Based on previous screens, PI 567301B is susceptible to SCN. Fine mapping of this locus will assist in cloning these candidate genes as well as identifying DNA markers flanking the QTL that can be used in marker-assisted breeding to develop cultivars with high levels of resistance to F. graminearum.

Evidence for trade-offs in detoxification and chemosensation gene signatures in Plutella xylostella

BACKGROUND

Detoxification genes have been associated with insecticide adaptation in the diamondback moth, Plutella xylostella. The link between chemosensation genes and adaptation, however, remains unexplored. To gain a better understanding of the involvement of these genes in insecticide adaptation, the authors exposed lines of P. xylostella to either high uniform (HU) or low heterogeneous (LH) concentrations of permethrin, expecting primarily physiological or behavioral selection respectively. Initially, 454 pyrosequencing was applied, followed by an examination of expression profiles of candidate genes that responded to selection [cytochrome P450 (CYP), glutathione S-transferase (GST), carboxylesterase (CarE), chemosensory protein (CSP) and odorant-binding protein (OBP)] by quantitative PCR in the larvae. Toxicity and behavioral assays were also conducted to document the effects of the two forms of exposure.

RESULTS

Pyrosequencing of the P. xylostella transcriptome from adult heads and third instars produced 198,753 reads with 52,752,486 bases. Quantitative PCR revealed overexpression of CYP4M14, CYP305B1 and CSP8 in HU larvae. OBP13, however, was highest in LH. Larvae from LH and HU lines had up to five- and 752-fold resistance levels respectively, which could be due to overexpression of P450s. However, the behavioral responses of all lines to a series of permethrin concentrations did not vary significantly in any of the generations examined, in spite of the observed upregulation of CSP8 and OBP13.

CONCLUSION

Expression patterns from the target genes provide insights into behavioral and physiological responses to permethrin and suggest a new avenue of research on the role of chemosensation genes in insect adaptation to toxins.

Genetic insights into Graminella nigrifrons Competence for maize fine streak virus infection and transmission

BACKGROUND

Most plant-infecting rhabdoviruses are transmitted by one or a few closely related insect species. Additionally, intraspecific differences in transmission efficacy often exist among races/biotypes within vector species and among strains within a virus species. The black-faced leafhopper, Graminella nigrifrons, is the only known vector of the persistent propagative rhabdovirus Maize fine streak virus (MFSV). Only a small percentage of leafhoppers are capable of transmitting the virus, although the mechanisms underlying vector competence are not well understood.

METHODOLOGY

RNA-Seq was carried out to explore transcript expression changes and sequence variation in G. nigrifrons and MFSV that may be associated with the ability of the vector to acquire and transmit the virus. RT-qPCR assays were used to validate differential transcript accumulation.

RESULTS/SIGNIFICANCE

Feeding on MFSV-infected maize elicited a considerable transcriptional response in G. nigrifrons, with increased expression of cytoskeleton organization and immunity transcripts in infected leafhoppers. Differences between leafhoppers capable of transmitting MFSV, relative to non-transmitting but infected leafhoppers were more limited, which may reflect difficulties discerning between the two groups and/or the likelihood that the transmitter phenotype results from one or a few genetic differences. The ability of infected leafhoppers to transmit MFSV did not appear associated with virus transcript accumulation in the infected leafhoppers or sequence polymorphisms in the viral genome. However, the non-structural MFSV 3 gene was expressed at unexpectedly high levels in infected leafhoppers, suggesting it plays an active role in the infection of the insect host. The results of this study begin to define the functional roles of specific G. nigrifrons and MFSV genes in the viral transmission process.

The crypt-dwelling primary bacterial symbiont of the polyphagous pentatomid pest Halyomorpha halys (Hemiptera: Pentatomidae)

A recent invader of North America, the brown marmorated stink bug (Halyomorpha halys Stål) is a polyphagous pentatomid that harbors a gammaproteobacterial mutualist in the crypts of specialized midgut gastric caeca (region V4). Histological analyses revealed a single rod-shaped morphology abundant in distal V4 midgut caecal crypts. A strong fluorescence signal was detected when thin sections of these tissues were hybridized with a fluorescently-labeled, Enterobacteriaceae-specific oligonucleotide probe. A single operational taxonomic unit (OTU) assigned to the Pantoea genus represented >99% of 3,454 16S rDNA amplicons obtained from midgut V4 tissues and egg samples. Detection of H. halys primary symbiont in DNA extracted from eggs suggested vertical maternal inheritance as the mode of intergenerational transmission. Consistent detection of the bacterial symbiont in geographically distinct H. halys populations strongly supports an intimate association between these two organisms. An inferred phylogeny of gammaproteobacterial symbionts of pentatomids placed the Pantoea-assigned OTU from H. halys within a clade distinct from primary bacterial symbionts of related stink bugs, Nezara viridula (L.) and Eurydema rugosa Motschulsky. Given these data, Candidatus "Pantoea carbekii" is proposed as the name of the primary bacterial symbiont of H. halys.

Reduction in fecundity and shifts in cellular processes by a native virus on an invasive insect

Pathogens and their vectors have coevolutionary histories that are intricately intertwined with their ecologies, environments, and genetic interactions. The soybean aphid, Aphis glycines, is native to East Asia but has quickly become one of the most important aphid pests in soybean-growing regions of North America. In this study, we used bioassays to examine the effects of feeding on soybean infected with a virus it vectors (Soybean mosaic virus [SMV]) and a virus it does not vector (Bean pod mottle virus [BPMV]) have on A. glycines survival and fecundity. The genetic underpinnings of the observed changes in fitness phenotype were explored using RNA-Seq. Aphids fed on SMV-infected soybean had transcriptome and fitness profiles that were similar to that of aphids fed on healthy control plants. Strikingly, a significant reduction in fecundity was seen in aphids fed on BPMV-infected soybean, concurrent with a large and persistent downregulation of A. glycines transcripts involved in regular cellular activities. Although molecular signatures suggested a small regulatory RNA pathway defense response was repressed in aphids feeding on infected plants, BPMV did not appear to be replicating in the vector. These results suggest that incompatibilities with BPMV or the effects of BPMV infection on soybean caused A. glycines to allot available energy resources to survival rather than reproduction and other core cellular processes. Ultimately, the detrimental impacts to A. glycines may reflect the short tritrophic evolutionary histories between the insect, plant, and virus.

Virus-independent and common transcriptome responses of leafhopper vectors feeding on maize infected with semi-persistently and persistent propagatively transmitted viruses

BACKGROUND

Insects are the most important epidemiological factors for plant virus disease spread, with >75% of viruses being dependent on insects for transmission to new hosts. The black-faced leafhopper (Graminella nigrifrons Forbes) transmits two viruses that use different strategies for transmission: Maize chlorotic dwarf virus (MCDV) which is semi-persistently transmitted and Maize fine streak virus (MFSV) which is persistently and propagatively transmitted. To date, little is known regarding the molecular and cellular mechanisms in insects that regulate the process and efficiency of transmission, or how these mechanisms differ based on virus transmission strategy.

RESULTS

RNA-Seq was used to examine transcript changes in leafhoppers after feeding on MCDV-infected, MFSV-infected and healthy maize for 4 h and 7 d. After sequencing cDNA libraries constructed from whole individuals using Illumina next generation sequencing, the Rnnotator pipeline in Galaxy was used to reassemble the G. nigrifrons transcriptome. Using differential expression analyses, we identified significant changes in transcript abundance in G. nigrifrons. In particular, transcripts implicated in the innate immune response and energy production were more highly expressed in insects fed on virus-infected maize. Leafhoppers fed on MFSV-infected maize also showed an induction of transcripts involved in hemocoel and cell-membrane linked immune responses within four hours of feeding. Patterns of transcript expression were validated for a subset of transcripts by quantitative real-time reverse transcription polymerase chain reaction using RNA samples collected from insects fed on healthy or virus-infected maize for between a 4 h and seven week period.

CONCLUSIONS

We expected, and found, changes in transcript expression in G. nigrifrons feeding of maize infected with a virus (MFSV) that also infects the leafhopper, including induction of immune responses in the hemocoel and at the cell membrane. The significant induction of the innate immune system in G. nigrifrons fed on a foregut-borne virus (MCDV) that does not infect leafhoppers was less expected. The changes in transcript accumulation that occur independent of the mode of pathogen transmission could be key for identifying insect factors that disrupt vector-mediated plant virus transmission.

Genetic mapping of three quantitative trait loci for soybean aphid resistance in PI 567324

Host-plant resistance is an effective method for controlling soybean aphid (Aphis glycines Matsumura), the most damaging insect pest of soybean (Glycine max (L.) Merr.) in North America. Recently, resistant soybean lines have been discovered and at least four aphid resistance genes (Rag1, Rag2, Rag3 and rag4) have been mapped on different soybean chromosomes. However, the evolution of new soybean aphid biotypes capable of defeating host-plant resistance conferred by most single genes demonstrates the need for finding germplasm with multigenic resistance to the aphid. This study was conducted to map quantitative trait loci (QTL) for aphid resistance in PI 567324. We identified two major QTL (QTL_13_1 and QTL_13_2) for aphid resistance on soybean chromosome 13 using 184 recombinant inbred lines from a 'Wyandot' × PI 567324 cross. QTL_13_1 was located close to the previously reported Rag2 gene locus, and QTL_13_2 was close to the rag4 locus. A minor QTL (QTL_6_1) was also detected on chromosome 6, where no gene for soybean aphid resistance has been reported so far. These results indicate that PI 567324 possesses oligogenic resistance to the soybean aphid. The molecular markers closely linked to the QTL reported here will be useful for development of cultivars with oligogenic resistance that are expected to provide broader and more durable resistance against soybean aphids compared with cultivars with monogenic resistance.

Population genetic structure and genetic diversity of soybean aphid collections from the USA, South Korea, and Japan

Following its recent invasion of North America, the soybean aphid (Aphis glycines Matsumura) has become the number one insect pest of soybean (Glycine max L. Merr.) in the north central states of the USA. A few studies have been conducted on the population genetic structure and genetic diversity of the soybean aphid and the source of its invasion in North America. Molecular markers, such as simple sequence repeats (SSRs) are very useful in the evaluation of population structure and genetic diversity. We used 18 SSR markers to assess the genetic diversity of soybean aphid collections from the USA, South Korea, and Japan. The aphids were collected from two sites in the USA (Indiana and South Dakota), two sites in South Korea (Yeonggwang district and Cheonan city), and one site in Japan (Utsunomiya). The SSR markers were highly effective in differentiating among aphid collections from different countries. The level of differentiation within each population and among populations from the same country was limited, even in the case of the USA where the two collection sites were more than 1200 km apart.

Permanent genetic resources added to Molecular Ecology Resources Database 1 August 2012 - 30 September 2012

This article documents the addition of 83 microsatellite marker loci and 96 pairs of single-nucleotide polymorphism (SNP) sequencing primers to the Molecular Ecology Resources Database. Loci were developed for the following species: Bembidion lampros, Inimicus japonicus, Lymnaea stagnalis, Panopea abbreviata, Pentadesma butyracea, Sycoscapter hirticola and Thanatephorus cucumeris (anamorph: Rhizoctonia solani). These loci were cross-tested on the following species: Pentadesma grandifolia and Pentadesma reyndersii. This article also documents the addition of 96 sequencing primer pairs and 88 allele-specific primers or probes for Plutella xylostella.

Genetic mapping of the powdery mildew resistance gene in soybean PI 567301B

Powdery mildew (PMD) of soybean [Glycine max (L.) Merr.] is caused by the fungus Microsphaera diffusa. Severe infection of PMD on susceptible varieties often causes premature defoliation and chlorosis of the leaves, which can result in considerable yield losses under favorable environmental conditions for disease development in the field. A total of 334 F(7)-derived recombinant inbred lines (RILs) from a cross of a PMD susceptible soybean cultivar Wyandot and PMD-resistant PI 567301B were used for genetic mapping of PMD resistance in PI 567301B and for development of molecular markers tightly linked to the gene. The result of the PMD screening for each line in the field was in agreement with that in the greenhouse test. The genetic map containing the PMD resistance gene was constructed in a 3.3 cM interval flanked by two simple sequence repeat (SSR) markers on chromosome 16. The PMD resistance gene was mapped at the same location with SSR marker BARCSOYSSR_16_1291, indicating that there was no recombination between the 334 RILs and this marker. In addition, a single nucleotide polymorphism (SNP) marker developed by high-resolution melting curve analysis and a cleaved amplified polymorphic sequence (CAPS) marker with Rsa1 recognition site were used for the genetic mapping. These two markers were also mapped to the same genomic location with the PMD resistance gene. We validated three tightly linked markers to the PMD resistance gene using 38 BC(6)F(2) lines and corresponding BC(6)F(2:3) families. The three marker genotypes of the backcross lines predicted the observed PMD phenotypes of the lines with complete accuracy. We have mapped a putatively novel single dominant PMD resistance gene in PI 567301B and developed three new molecular markers closely linked to the gene. Molecular markers developed from this study may be used for high-throughput marker-assisted breeding for PMD resistance with the gene from PI 567301B.

Maintaining genetic diversity and population panmixia through dispersal and not gene flow in a holocyclic heteroecious aphid species

Heteroecious holocyclic aphids exhibit both sexual and asexual reproduction and alternate among primary and secondary hosts. Most of these aphids can feed on several related hosts, and invasions to new habitats may limit the number of suitable hosts. For example, the aphid specialist Aphis glycines survives only on the primary host buckthorn (Rhamnus spp.) and the secondary host soybean (Glycine max) in North America where it is invasive. Owing to this specialization and sparse primary host distribution, host colonization events could be localized and involve founder effects, impacting genetic diversity, population structure and adaptation. We characterized changes in the genetic diversity and structure across time among A. glycines populations. Populations were sampled from secondary hosts twice in the same geographical location: once after secondary colonization (early season), and again immediately before primary host colonization (late season). We tested for evidence of founder effects and genetic isolation in early season populations, and whether or not late-season dispersal restored genetic diversity and reduced fragmentation. A total of 24 single-nucleotide polymorphisms and 6 microsatellites were used for population genetic statistics. We found significantly lower levels of genotypic diversity and more genetic isolation among early season collections, indicating secondary host colonization occurred locally and involved founder effects. Pairwise F(ST) decreased from 0.046 to 0.017 in early and late collections, respectively, and while genetic relatedness significantly decreased with geographical distance in early season collections, no spatial structure was observed in late-season collections. Thus, late-season dispersal counteracts the secondary host colonization through homogenization and increases genetic diversity before primary host colonization.

Genetic mapping revealed two loci for soybean aphid resistance in PI 567301B

The soybean aphid (Aphis glycines Matsumura) is the most damaging insect pest of soybean [Glycine max (L.) Merr.] in North America. New soybean aphid biotypes have been evolving quickly and at least three confirmed biotypes have been reported in USA. These biotypes are capable of defeating most known aphid resistant soybean genes indicating the need for identification of new genes. Plant Introduction (PI) 567301B was earlier identified to have antixenosis resistance against biotype 1 and 2 of the soybean aphid. Two hundred and three F(7:9) recombinant inbred lines (RILs) developed from a cross of soybean aphid susceptible cultivar Wyandot and resistant PI 567301B were used for mapping aphid resistance genes using the quantitative trait loci (QTL) mapping approach. A subset of 94 RILs and 516 polymorphic SNP makers were used to construct a genome-wide molecular linkage map. Two candidate QTL regions for aphid resistance were identified on this linkage map. Fine mapping of the QTL regions was conducted with SSR markers using all 203 RILs. A major gene on chromosome 13 was mapped near the previously identified Rag2 gene. However, an earlier study revealed that the detached leaves of PI 567301B had no resistance against the soybean aphids while the detached leaves of PI 243540 (source of Rag2) maintained aphid resistance. These results and the earlier finding that PI 243540 showed antibiosis resistance and PI 567301B showed antixenosis type resistance, indicating that the aphid resistances in the two PIs are not controlled by the same gene. Thus, we have mapped a new gene near the Rag2 locus for soybean aphid resistance that should be useful in breeding for new aphid-resistant soybean cultivars. Molecular markers closely linked to this gene are available for marker-assisted breeding. Also, the minor locus found on chromosome 8 represents the first reported soybean aphid-resistant locus on this chromosome.

Cloning and characterization of mariner-like elements in the soybean aphid, Aphis glycines Matsumura

Soybean aphid, Aphis glycines Matsumura (Hemiptera: Aphididae), is currently the most important insect pest of soybean (Glycine max (L.) Merr.) in the United States and causes significant economic damage worldwide, but little is known about the aphid at the molecular level. Mariner-like transposable elements (MLEs) are ubiquitous within the genomes of arthropods and various other invertebrates. In this study, we report the cloning of MLEs from the soybean aphid genome using degenerate PCR primers designed to amplify conserved regions of mariner transposases. Two of the ten sequenced clones (designated as Agmar1 and Agmar2) contained partial but continuous open reading frames, which shared high levels of homology at the protein level with other mariner transposases from insects and other taxa. Phylogenetic analysis revealed Agmar1 to group within the irritans subfamily of MLEs and Agmar2 within the mellifera subfamily. Southern blot analysis and quantitative PCR analysis indicated a low copy number for Agmar1-like elements within the soybean aphid genome. These results suggest the presence of at least two different putative mariner-like transposases encoded by the soybean aphid genome. Both Agmar1 and Agmar2 could play influential roles in the architecture of the soybean aphid genome. Transposable elements are also thought to potentially mediate resistance in insects through changes in gene amplification and mutations in coding sequences. Finally, Agmar1 and Agmar2 may represent useful genetic tools and provide insights on A. glycines adaptation.

Permanent Genetic Resources added to Molecular Ecology Resources Database 1 June 2011-31 July 2011

This article documents the addition of 112 microsatellite marker loci and 24 pairs of single nucleotide polymorphism (SNP) sequencing primers to the Molecular Ecology Resources Database. Loci were developed for the following species: Agelaius phoeniceus, Austrolittorina cincta, Circus cyaneus, Circus macrourus, Circus pygargus, Cryptocoryne × purpurea Ridl. nothovar. purpurea, Mya arenaria, Patagioenas squamosa, Prochilodus mariae, Scylla serrata and Scytalopus speluncae. These loci were cross-tested on the following species: Cryptocoryne × purpurea nothovar. purpurea, Cryptocoryne affinis, Cryptocoryne ciliata, Cryptocoryne cordata var. cordata, Cryptocoryne elliptica, Cryptocoryne griffithii, Cryptocoryne minima, Cryptocoryne nurii and Cryptocoryne schulzei. This article also documents the addition of 24 sequencing primer pairs and 24 allele-specific primers or probes for Aphis glycines.

Development of soybean aphid genomic SSR markers using next generation sequencing

Simple sequence repeats (SSRs) or microsatellites are very useful molecular markers, owing to their locus-specific codominant and multiallelic nature, high abundance in the genome, and high rates of transferability across species. The soybean aphid (Aphis glycines Matsumura) has become the most damaging insect pest of soybean (Glycine max (L.) Merr.) in North America, since it was first found in the Midwest of the United States in 2000. Biotypes of the soybean aphid capable of colonizing newly developed aphid-resistant soybean cultivars have been recently discovered. Genetic resources, including molecular markers, to study soybean aphids are severely lacking. Recently developed next generation sequencing platforms offer opportunities for high-throughput and inexpensive genome sequencing and rapid marker development. The objectives of this study were (i) to develop and characterize genomic SSR markers from soybean aphid genomic sequences generated by next generation sequencing technology and (ii) to evaluate the utility of the SSRs for genetic diversity or relationship analyses. In total 128 SSR primer pairs were designed from sequences generated by Illumina GAII from a reduced representation library of A. glycines. Nearly 94% (120) of the primer pairs amplified SSR alleles of expected size and 24 SSR loci were polymorphic among three aphid samples from three populations. The polymorphic SSRs were successfully used to differentiate among 24 soybean aphids from Ohio and South Dakota. Sequencing of PCR products of two SSR markers from four aphid samples revealed that the allelic polymorphism was due to variation in the SSR repeats among the aphids. These markers should be particularly useful for genetic differentiation among aphids collected from soybean fields at different localities and regions. These SSR markers provide the soybean aphid research community with the first set of PCR-based codominant markers developed from the genomic sequences of A. glycines.

Genetic diversity and differentiation among laboratory and field populations of the soybean aphid, Aphis glycines

The soybean aphid, Aphis glycines Matsumura, is a recent invasive pest of soybean in North America. Currently, much research is focused on developing and characterizing soybean cultivars expressing host-plant resistance. During the initial phases of host-plant resistance screening, many of these studies use soybean aphid laboratory populations. Previous studies in other systems have documented substantial differences among laboratory and field populations. Whether or not this pattern exists in A. glycines is unknown, but it is extremely important when estimating the level of selection and virulence to host-plant resistant soybeans. In this study, we used seven microsatellite markers to estimate and compare genetic diversity and differentiation among five laboratory and 12 field populations. Our results indicate that soybean aphid laboratory populations are severely lacking in genotypic diversity and show extreme genetic differentiation among each other and to field populations. Continued use of laboratory populations for initial soybean aphid resistance screening could lead to erroneous estimations of the potential success for host-plant resistance.

Combining next-generation sequencing strategies for rapid molecular resource development from an invasive aphid species, Aphis glycines

Background

Aphids are one of the most important insect taxa in terms of ecology, evolutionary biology, genetics and genomics, and interactions with endosymbionts. Additionally, many aphids are serious pest species of agricultural and horticultural plants. Recent genetic and genomic research has expanded molecular resources for many aphid species, including the whole genome sequencing of the pea aphid, Acrythosiphon pisum. However, the invasive soybean aphid, Aphis glycines, lacks in any significant molecular resources.

Methodology/Principal Findings

Two next-generation sequencing technologies (Roche-454 and Illumina GA-II) were used in a combined approach to develop both transcriptomic and genomic resources, including expressed genes and molecular markers. Over 278 million bp were sequenced among the two methods, resulting in 19,293 transcripts and 56,688 genomic sequences. From this data set, 635 SNPs and 1,382 microsatellite markers were identified. For each sequencing method, different soybean aphid biotypes were used which revealed potential biotype specific markers. In addition, we uncovered 39,822 bp of sequence that were related to the obligatory endosymbiont, Buchnera aphidicola, as well as sequences that suggest the presence of Hamiltonella defensa, a facultative endosymbiont.

Conclusions and Significance

Molecular resources for an invasive, non-model aphid species were generated. Additionally, the power of next-generation sequencing to uncover endosymbionts was demonstrated. The resources presented here will complement ongoing molecular studies within the Aphididae, including the pea aphid whole genome, lead to better understanding of aphid adaptation and evolution, and help provide novel targets for soybean aphid control.

Detached leaf and whole plant assays for soybean aphid resistance: differential responses among resistance sources and biotypes

The soybean aphid, Aphis glycines Matsumura, is a pest of cultivated soybean, Glycine max (L.) Merr., in North America. Recent developments in host plant resistance studies have identified at least four soybean aphid resistance genes (Rag1-4) and two soybean aphid biotypes (biotype 1 and 2), defined by differential survivability on resistant soybean. Detached soybean leaves were tested as a more rapid and practical assay to assess host plant resistance and virulence. Two susceptible lines ('Wyandot' and 'Williams 82') and two resistant lines (PI 243540 and PI 567301B) were examined. Various life history characteristics were compared among aphids on whole plants and detached leaves. Results indicated that resistance to soybean aphid was lost using detached leaves of PI 567301B but retained with PI 243540. To test for aphid virulence, net fecundities were compared among biotype 1 and biotype 2 after rearing on detached leaves of the resistant 'Jackson' (to which biotype 2 is virulent). A significant difference was detected in net fecundities among biotypes on detached leaves of Jackson and used to predict growth rates and virulence from 30 field-collected individuals of unknown virulence. No field individuals matched biotype 2 predictions, but four individuals had higher net fecundities than biotype 2 predictions (13%) and could be considered moderately virulent. The results indicated that the retention of soybean aphid resistance in detached leaves is dependent on PI and resistant source, but if resistance is retained, detached leaves could be used to determine soybean aphid virulence.

Population genetic structure of Aphis glycines

The soybean aphid (Aphis glycines Matsumura) is an invasive pest of cultivated soybean (Glycine max L.) in North America. After the initial invasion in 2000, the aphid has quickly spread across most of the United States and Canada, suggesting large-scale dispersal and rapid adaptation to new environments. Using microsatellite markers from closely related species, we compared the genetic diversity and the amount of genetic differentiation within and among 2 South Korean and 10 North American populations. Overall allelic polymorphism was low, never exceeding four alleles per locus. However, differences in genetic diversity were seen among South Korean and North American populations in terms of heterozygote excesses and genotypic richness. Within North America, two populations (Michigan and Ontario), had lower genetic diversities and exhibited high genetic differentiation compared with the remaining eight populations. The earlier collection time of Michigan and Ontario samples explained the genetic differences better than geographic subdivisions. These data indicate a pattern of small colonizing populations on soybeans, followed by rapid clonal amplification and subsequent large-scale dispersal across North America.

Cross-species amplification and polymorphism of microsatellite loci in the soybean aphid, Aphis glycines

We tested the utility of 18 previously characterized Aphis spp. microsatellite loci for polymorphism and differentiation among populations of the soybean aphid, Aphis glycines. Loci were chosen from a closely related species (Aphis gossypii) and a more distantly related species (Aphis fabae). We found nine loci to be polymorphic among Korean and North American populations. Overall expected heterozygosity was moderate (average: 0.47; range: 0-1), although populations substantially differed in deviations from Hardy-Weinberg equilibrium. These loci will be valuable in characterizing population differentiation, migration and adaptation in an economically important pest of soybeans.

The genetic structure of hawthorn-infesting Rhagoletis pomonella populations in Mexico: implications for sympatric host race formation

The genetic origins of species may not all trace to the same time and place as the proximate cause(s) for population divergence. Moreover, inherent gene-flow barriers separating populations may not all have evolved under the same geographical circumstances. These considerations have lead to a greater appreciation of the plurality of speciation: that one geographical mode for divergence may not always be sufficient to describe a speciation event. The apple maggot fly, Rhagoletis pomonella, a model system for sympatric speciation via host-plant shifting, has been a surprising contributor to the concept of speciation mode plurality. Previous studies have suggested that past introgression of inversion polymorphism from a hawthorn-fly population in the trans-Mexican volcanic belt (EVTM) introduced diapause life-history variation into a more northern fly population that subsequently contributed to sympatric host race formation and speciation in the United States (US). Here, we report results from a microsatellite survey implying (i) that volcanic activity in the eastern EVTM may have been responsible for the initial geographical isolation of the Mexican and northern hawthorn-fly populations c. 1.57 mya; and (ii) that flies in the Sierra Madre Oriental Mountains (SMO) likely served as a conduit for past gene flow from the EVTM into the US. Indeed, the microsatellite data suggest that the current US population may represent a range expansion from the northern SMO. We discuss the implications of these findings for sympatric race formation in Rhagoletis and speciation theory.

Hawthorn-infesting populations of Rhagoletis pomonella in Mexico and speciation mode plurality

Categorizing speciation into dichotomous allopatric versus nonallopatric modes may not always adequately describe the geographic context of divergence for taxa. If some of the genetic changes generating inherent barriers to gene flow between populations evolved in geographic isolation, whereas others arose in sympatry, then the mode of divergence would be mixed. The apple maggot fly, Rhagoletis pomonella, has contributed to this emerging concept of a mixed speciation mode "plurality." Genetic studies have implied that a source of diapause life-history variation associated with inversions and contributing to sympatric host race formation and speciation for R. pomonella in the United States may have introgressed from the Eje Volcanico Trans Mexicano (EVTM; a.k.a. the Altiplano) in the past. A critical unresolved issue concerning the introgression hypothesis is how past gene flow occurred given the current 1200-km disjunction in the ranges of hawthorn-infesting flies in the EVTM region of Mexico and the southern extreme of the U.S. population in Texas. Here, we report the discovery of a hawthorn-infesting population of R. pomonella in the Sierra Madre Oriental Mountains (SMO) of Mexico. Sequence data from 15 nuclear loci and mitochondrial DNA imply that the SMO flies are related to, but still different from, U.S. and EVTM flies. The host affiliations, diapause characteristics, and phylogeography of the SMO population are consistent with it having served as a conduit for gene flow between Mexico and the United States. We also present evidence suggesting greater permeability of collinear versus rearranged regions of the genome to introgression, in accord with recent models of chromosomal speciation. We discuss the implications of the results in the context of speciation mode plurality. We do not argue for abandoning the terms sympatry or allopatry, but caution that categorizing divergence into either/or geographic modes may not describe the genetic origins of all species. For R. pomonella in the United States, the proximate selection pressures triggering race formation and speciation stem from sympatric host shifts. However, some of the phenological variation contributing to host-related ecological adaptation and reproductive isolation in sympatry at the present time appears to have an older history, having originated and become packaged into inversion polymorphism in allopatry.

Rangewide population genetic structure of the African malaria vector Anopheles funestus

Anopheles funestus is a primary vector of malaria in Africa south of the Sahara. We assessed its rangewide population genetic structure based on samples from 11 countries, using 10 physically mapped microsatellite loci, two per autosome arm and the X (N = 548), and 834 bp of the mitochondrial ND5 gene (N = 470). On the basis of microsatellite allele frequencies, we found three subdivisions: eastern (coastal Tanzania, Malawi, Mozambique and Madagascar), western (Burkina Faso, Mali, Nigeria and western Kenya), and central (Gabon, coastal Angola). A. funestus from the southwest of Uganda had affinities to all three subdivisions. Mitochondrial DNA (mtDNA) corroborated this structure, although mtDNA gene trees showed less resolution. The eastern subdivision had significantly lower diversity, similar to the pattern found in the codistributed malaria vector Anopheles gambiae. This suggests that both species have responded to common geographic and/or climatic constraints. The western division showed signatures of population expansion encompassing Kenya west of the Rift Valley through Burkina Faso and Mali. This pattern also bears similarity to A. gambiae, and may reflect a common response to expanding human populations following the development of agriculture. Due to the presumed recent population expansion, the correlation between genetic and geographic distance was weak. Mitochondrial DNA revealed further cryptic subdivision in A. funestus, not detected in the nuclear genome. Mozambique and Madagascar samples contained two mtDNA lineages, designated clade I and clade II, that were separated by two fixed differences and an average of 2% divergence, which implies that they have evolved independently for approximately 1 million years. Clade I was found in all 11 locations, whereas clade II was sampled only on Madagascar and Mozambique. We suggest that the latter clade may represent mtDNA capture by A. funestus, resulting from historical gene flow either among previously isolated and divergent populations or with a related species.

Divergence with gene flow in Anopheles funestus from the Sudan Savanna of Burkina Faso, West Africa

Anopheles funestus is a major vector of malaria across Africa. Understanding its complex and nonequilibrium population genetic structure is an important challenge that must be overcome before vector populations can be successfully perturbed for malaria control. Here we examine the role of chromosomal inversions in structuring genetic variation and facilitating divergence in Burkina Faso, West Africa, where two incipient species (chromosomal forms) of A. funestus, defined principally by rearrangements of chromosome 3R, have been hypothesized. Sampling across an approximately 300-km east-west transect largely contained within the Sudan-Savanna ecoclimatic zone, we analyzed chromosomal inversions, 16 microsatellite loci distributed genomewide, and 834 bp of the mtDNA ND5 gene. Both molecular markers revealed high genetic diversity, nearly all of which was accounted for by within-population differences among individuals, owing to recent population expansion. Across the study area there was no correlation between genetic and geographic distance. Significant genetic differentiation found between chromosomal forms on the basis of microsatellites was not genomewide but could be explained by chromosome 3R alone on the basis of loci inside and near inversions. These data are not compatible with complete reproductive isolation but are consistent with differential introgression and sympatric divergence between the chromosomal forms, facilitated by chromosome 3R inversions.

Genetic structure of Anopheles gambiae populations on islands in northwestern Lake Victoria, Uganda

Background

Alternative means of malaria control are urgently needed. Evaluating the effectiveness of measures that involve genetic manipulation of vector populations will be facilitated by identifying small, genetically isolated vector populations. The study was designed to use variation in microsatellite markers to look at genetic structure across four Lake Victoria islands and two surrounding mainland populations and for evidence of any restriction to free gene flow.

Methods

Four Islands (from 20–50 km apart) and two surrounding mainland populations (96 km apart) were studied. Samples of indoor resting adult mosquitoes, collected over two consecutive years, were genotyped at microsatellite loci distributed broadly throughout the genome and analysed for genetic structure, effective migration (Nem) and effective population size (Ne).

Results

Ne estimates showed island populations to consist of smaller demes compared to the mainland ones. Most populations were significantly differentiated geographically, and from one year to the other. Average geographic pair-wise FST ranged from 0.014–0.105 and several pairs of populations had Ne m < 3. The loci showed broad heterogeneity at capturing or estimating population differences.

Conclusion

These island populations are significantly genetically differentiated. Differences reoccurred over the study period, between the two mainland populations and between each other. This appears to be the product of their separation by water, dynamics of small populations and local adaptation. With further characterisation these islands could become possible sites for applying measures evaluating effectiveness of control by genetic manipulation.

Molecular differentiation between chromosomally defined incipient species of Anopheles funestus

Anopheles funestus Giles is one of the most important vectors of malaria in sub-Saharan Africa. The population structure of this mosquito in Burkina Faso, West Africa based on chromosomal inversion data led to the description of two chromosomal forms, Kiribina and Folonzo. Because both forms co-occur in the same locales yet differ significantly, both in the frequency of inverted arrangements on chromosome arms 3R and 2R and in vectorial capacity, they were hypothesized to be emerging species with at least partial barriers to gene flow. This hypothesis would be strengthened by molecular evidence of differentiation between Kiribina and Folonzo at loci outside chromosomal inversions. We surveyed molecular variation in sympatric populations of the two forms using sequences from the mitochondrial ND5 gene and genotypes at sixteen microsatellite loci distributed across the genome. Both classes of marker revealed slight but significant differentiation between the two forms (mtDNA F(ST) = 0.023, P < 0.001; microsatellite F(ST) = 0.004, P < 0.001; R(st) = 0.009, P = 0.002). Locus-by-locus analysis of the microsatellite data showed that significant differentiation was not genome-wide, but could be attributed to five loci on chromosome 3R (F(ST) = 0.010, P < 0.001; R(st) = 0.016, P = 0.002). Importantly, three of these loci are outside of, and in linkage equilibrium with, chromosomal inversions, suggesting that differentiation between chromosomal forms extends beyond the inversions themselves. The slight overall degree of differentiation indicated by both marker classes is likely an underestimate because of recent population expansion inferred for both Folonzo and Kiribina. The molecular evidence from this study is consistent with the hypothesis of incipient speciation between Kiribina and Folonzo.