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

Complex evolution in Aphis gossypii group (Hemiptera: Aphididae), evidence of primary host shift and hybridization between sympatric species

Aphids provide a good model system to understand the ecological speciation concept, since the majority of the species are host-specific, and they spend their entire lifecycle on certain groups of host plants. Aphid species that apparently have wide host plant ranges have often turned out to be complexes of host-specialized biotypes. Here we investigated the various host-associated populations of the two recently diverged species, Aphis gossypii and A. rhamnicola, having multiple primary hosts, to understand the complex evolution with host-associated speciation. Using mitochondrial DNA marker and nine microsatellite loci, we reconstructed the haplotype network, and analyzed the genetic structure and relationships. Approximate Bayesian computation was also used to infer the ancestral primary host and host-associated divergence, which resulted in Rhamnus being the most ancestral host for A. gossypii and A. rhamnicola. As a result, Aphis gossypii and A. rhamnicola do not randomly use their primary and secondary host plants; rather, certain biotypes use only some secondary and specific primary hosts. Some biotypes are possibly in a diverging state through specialization to specific primary hosts. Our results also indicate that a new heteroecious race can commonly be derived from the heteroecious ancestor, showing strong evidence of ecological specialization through a primary host shift in both A. gossypii and A. rhamnicola. Interestingly, A. gossypii and A. rhamnicola shared COI haplotypes with each other, thus there is a possibility of introgression by hybridization between them by cross-sharing same primary hosts. Our results contribute to a new perspective in the study of aphid evolution by identifying complex evolutionary trends in the gossypii sensu lato complex.

Transcriptional responses of soybean aphids to sublethal insecticide exposure

Insecticides are a key tool in the management of many insect pests of agriculture, including soybean aphids. The selection imposed by insecticide use has often lead to the evolution of resistance by the target pest through enhanced detoxification mechanisms. We hypothesised that exposure of insecticide-susceptible aphids to sublethal doses of insecticides would result in the up-regulation of genes involved in detoxification of insecticides, revealing the genes upon which selection might act in the field. We used the soybean aphid biotype 1 reference genome, version 6.0 as a reference to analyze RNA-Seq data. We identified multiple genes with potential detoxification roles that were up-regulated 12 h after sublethal exposure to esfenvalerate or thiamethoxam. However, these genes were part of a dramatic burst of differential gene expression in which thousands of genes were up- or down-regulated, rather than a defined response to insecticides. Interestingly, the transcriptional burst observed at 12 h s declined dramatically by 24-hrs post-exposure, suggesting a general stress response that may become fine-tuned over time.

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.

Stress-induced changes in abundance differ among obligate and facultative endosymbionts of the soybean aphid

Bacterial endosymbionts can drive evolutionary novelty by conferring adaptive benefits under adverse environmental conditions. Among aphid species there is growing evidence that symbionts influence tolerance to various forms of stress. However, the extent to which stress inflicted on the aphid host has cascading effects on symbiont community dynamics remains poorly understood. Here we simultaneously quantified the effect of host-plant induced and xenobiotic stress on soybean aphid (Aphis glycines) fitness and relative abundance of its three bacterial symbionts. Exposure to soybean defensive stress (Rag1 gene) and a neurotoxic insecticide (thiamethoxam) substantially reduced aphid composite fitness (survival × reproduction) by 74 ± 10% and 92 ± 2%, respectively, which in turn induced distinctive changes in the endosymbiont microbiota. When challenged by host-plant defenses a 1.4-fold reduction in abundance of the obligate symbiont Buchnera was observed across four aphid clonal lines. Among facultative symbionts of Rag1-stressed aphids, Wolbachia abundance increased twofold and Arsenophonus decreased 1.5-fold. A similar pattern was observed under xenobiotic stress, with Buchnera and Arsenophonus titers decreasing (1.3-fold) and Wolbachia increasing (1.5-fold). Furthermore, variation in aphid virulence to Rag1 was positively correlated with changes in Arsenophonus titers, but not Wolbachia or Buchnera. A single Arsenophonus multi-locus genotype was found among aphid clonal lines, indicating strain diversity is not primarily responsible for correlated host-symbiont stress levels. Overall, our results demonstrate the nature of aphid symbioses can significantly affect the outcome of interactions under stress and suggests general changes in the microbiome can occur across multiple stress types.

Soybean aphid (Hemiptera: Aphididae) response to soybean plant defense: stress levels, tradeoffs, and cross-virulence

A variety of management methods to control the soybean aphid (Aphis glycines Matsumura) have been investigated since its invasion into North America in 2000, among them plant resistance has emerged as a viable option for reducing aphid damage to soybeans and preventing outbreaks. Plant resistance methods often use natural soybean plant defenses that impose stress on aphids by reducing fitness and altering behavior. Research efforts have heavily focused on identification and development of aphid resistant soybean varieties, leaving much unknown about soybean aphid response to stressful host plant defenses. In this study, we aimed to 1) evaluate lifetime fitness consequences and phenotypic variation in response to host plant-induced stress and 2) investigate whether trade-offs involving fitness costs and/or cross-virulence to multiple antibiotic soybean varieties exists. We compared aphid survival and reproduction during and after a short period of exposure to soybeans with the Rag2 resistance gene and measured aphid clonal variation in response to Rag2 soybeans. In addition, we measured the performance of Rag2 virulent and avirulent aphids on five soybean varieties with various forms of antibiotic resistance. Our results indicate that plant defenses impose high levels of stress and have long-term fitness consequences, even after aphids are removed from resistant plants. We identified one aphid clone that was able to colonize Rag2 among the seven clones tested, suggesting that virulent genotypes may be prevalent in natural populations. Finally, although we did not find evidence of cross-virulence to multiple antibiotic soybean varieties, our results suggest independent mechanisms of aphid virulence to Rag1 and Rag2 that may involve fitness costs.

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.

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.

Transcriptomic signatures of ash (Fraxinus spp.) phloem

BACKGROUND

Ash (Fraxinus spp.) is a dominant tree species throughout urban and forested landscapes of North America (NA). The rapid invasion of NA by emerald ash borer (Agrilus planipennis), a wood-boring beetle endemic to Eastern Asia, has resulted in the death of millions of ash trees and threatens billions more. Larvae feed primarily on phloem tissue, which girdles and kills the tree. While NA ash species including black (F. nigra), green (F. pennsylvannica) and white (F. americana) are highly susceptible, the Asian species Manchurian ash (F. mandshurica) is resistant to A. planipennis perhaps due to their co-evolutionary history. Little is known about the molecular genetics of ash. Hence, we undertook a functional genomics approach to identify the repertoire of genes expressed in ash phloem.

METHODOLOGY AND PRINCIPAL FINDINGS

Using 454 pyrosequencing we obtained 58,673 high quality ash sequences from pooled phloem samples of green, white, black, blue and Manchurian ash. Intriguingly, 45% of the deduced proteins were not significantly similar to any sequences in the GenBank non-redundant database. KEGG analysis of the ash sequences revealed a high occurrence of defense related genes. Expression analysis of early regulators potentially involved in plant defense (i.e. transcription factors, calcium dependent protein kinases and a lipoxygenase 3) revealed higher mRNA levels in resistant ash compared to susceptible ash species. Lastly, we predicted a total of 1,272 single nucleotide polymorphisms and 980 microsatellite loci, among which seven microsatellite loci showed polymorphism between different ash species.

CONCLUSIONS AND SIGNIFICANCE

The current transcriptomic data provide an invaluable resource for understanding the genetic make-up of ash phloem, the target tissue of A. planipennis. These data along with future functional studies could lead to the identification/characterization of defense genes involved in resistance of ash to A. planipennis, and in future ash breeding programs for marker development.

Ecology and management of the soybean aphid in North America

The soybean aphid, Aphis glycines Matsumura, has become the single most important arthropod pest of soybeans in North America. Native to Asia, this invasive species was first discovered in North America in July 2000 and has rapidly spread throughout the northcentral United States, much of southeastern Canada, and the northeastern United States. In response, important elements of the ecology of the soybean aphid in North America have been elucidated, with economic thresholds, sampling plans, and chemical control recommendations widely adopted. Aphid-resistant soybean varieties were available to growers in 2010. The preexisting community of aphid natural enemies has been highly effective in suppressing aphid populations in many situations, and classical biological control efforts have focused on the addition of parasitoids of Asian origin. The keys to sustainable management of this pest include understanding linkages between the soybean aphid and other introduced and native species in a landscape context along with continued development of aphid-resistant varieties.

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.

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.