A multi-genome analysis approach enables tracking of the invasion of a single Russian wheat aphid (Diuraphis noxia) clone throughout the New World

A layover in Europe: Reconstructing the invasion route of asexual lineages of a New Zealand snail to North America

Non-native invasive species are threatening ecosystems and biodiversity worldwide. High genetic variation is thought to be a critical factor for invasion success. Accordingly, the global invasion of a few clonal lineages of the gastropod Potamopyrgus antipodarum is thus both puzzling and has the potential to help illuminate why some invasions succeed while others fail. Here, we used SNP markers and a geographically broad sampling scheme (N = 1617) including native New Zealand populations and invasive North American and European populations to provide the first widescale population genetic assessment of the relationships between and among native and invasive P. antipodarum. We used a combination of traditional and Bayesian molecular analyses to demonstrate that New Zealand populations harbour very high diversity relative to the invasive populations and are the source of the two main European genetic lineages. One of these two European lineages was in turn the source of at least one of the two main North American genetic clusters of invasive P. antipodarum, located in Lake Ontario. The other widespread North American group had a more complex origin that included the other European lineage and two New Zealand clusters. Altogether, our analyses suggest that just a small handful of clonal lineages of P. antipodarum were responsible for invasion across continents. Our findings provide critical information for prevention of additional invasions and control of existing invasive populations and are of broader relevance towards understanding the establishment and evolution of asexual populations and the forces driving biological invasion.

Range-wide population genomics of the Mexican fruit fly: Toward development of pathway analysis tools

Recurrently invading pests provide unique challenges for pest management, but also present opportunities to utilize genomics to understand invasion dynamics and inform regulatory management through pathway analysis. In the southern United States, the Mexican fruit fly Anastrepha ludens is such a pest, and its incursions into Texas and California represent major threats to the agricultural systems of those regions. We developed a draft genome assembly for A. ludens, conducted range-wide population genomics using restriction site-associated DNA sequencing, and then developed and demonstrated a panel of highly differentiated diagnostic SNPs for source determination of intercepted flies in this system. Using 2,081 genomewide SNPs, we identified four populations across the range of A. ludens, corresponding to western Mexico, eastern Mexico/Texas, Guatemala/Belize/Honduras, and Costa Rica/Panama, with some intergradation present between clusters, particularly in Central America. From this population genomics framework, we developed a diagnostic panel of 28 highly differentiated SNPs that were able to recreate the genomewide population structure in this species. We demonstrated this panel on a set of test specimens, including specimens intercepted as part of regular trapping surveillance in Texas and California, and we were able to predict populations of origin for these specimens. This methodology presents a highly applied use of genomic techniques and can be implemented in any group of recurrently invading pests.

The potential distribution of the Russian wheat aphid (Diuraphis noxia): an updated distribution model including irrigation improves model fit for predicting potential spread

Diuraphis noxia (Kurdjumov), Russian wheat aphid, is one of the world's most invasive and economically important agricultural pests of wheat and barley. In May 2016, it was found for the first time in Australia, with further sampling confirming it was widespread throughout south-eastern regions. Russian wheat aphid is not yet present in New Zealand. The impacts of this pest if it establishes in New Zealand, could result in serious control problems in wheat- and barley-growing regions. To evaluate whether D. noxia could establish populations in New Zealand we used the climate modelling software CLIMEX to locate where potential viable populations might occur. We re-parameterised the existing CLIMEX model by Hughes and Maywald (1990) by improving the model fit using currently known distribution records of D. noxia, and we also considered the role of irrigation into the potential spread of this invasive insect. The updated model now fits the current known distribution better than the previous Hughes and Maywald CLIMEX model, particularly in temperate and Mediterranean areas in Australia and Europe; and in more semi-arid areas in north-western China and Middle Eastern countries. Our model also highlights new climatically suitable areas for the establishment of D. noxia, not previously reported, including parts of France, the UK and New Zealand. Our results suggest that, when suitable host plants are present, Russian wheat aphid could establish in these regions. The new CLIMEX projections in the present study are useful tools to inform risk assessments and target surveillance and monitoring efforts for identifying susceptible areas to invasion by Russian wheat aphid.

Resistance to insect pests in wheat-rye and Aegilops speltoides Tausch translocation and substitution lines

Various insect pests attack wheat (Triticum aestivum L.) that can cause significant grain yield losses to the crop. Farmers usually depend on pesticides, however, smallholder farmers often have limited and ill-timed access to control methods, including insecticides. Host plant resistance is an alternative to protect grain yield and reduce costs to farmers. Three of the most serious pests of wheat are Diuraphis noxia (Kurdjumov), Mayetiola destructor (Say), and Cephus pygmeaus L. These pests occur in most of the wheat growing areas. However, they are of high importance in North Africa and West Asia. The aim of this study was to evaluate a set of wheat-alien translocations for resistance against D. noxia, M. destructor and C. pygmeaus. Genotypes of this germplasm set have already been reported to carry resistance against certain wheat aphid species. Genotypes 1RS_am.1AL and MA1S.1RL_e(1B), displayed high levels of resistance against D. noxia and C. pygmeaus, respectively. While three genotypes showed resistance reaction against M. destructor: 1R_e(1D), 7A.7S-L5, and 7A.7S-Gb5. Except for the resistant genotype to C. pygmeaus, the other four genotypes were previously reported to carry resistance against Sitobion avenae Fabricius, Rhopalosiphum padi L. and Schizaphis graminum (Rondani). These five wheat-alien translocations are currently being used in the bread-wheat breeding programs at CIMMYT and ICARDA to transfer the multiple pest resistance in elite germplasm.

Population structure of a global agricultural invasive pest, Bactrocera dorsalis (Diptera: Tephritidae)

Bactrocera dorsalis, the Oriental fruit fly, is one of the world's most destructive agricultural insect pests and a major impediment to international fresh commodity trade. The genetic structuring of the species across its entire geographic range has never been undertaken, because under a former taxonomy B. dorsalis was divided into four distinct taxonomic entities, each with their own, largely non-overlapping, distributions. Based on the extensive sampling of six a priori groups from 63 locations, genetic and geometric morphometric datasets were generated to detect macrogeographic population structure, and to determine prior and current invasion pathways of this species. Weak population structure and high genetic diversity were detected among Asian populations. Invasive populations in Africa and Hawaii are inferred to be the result of separate, single invasions from South Asia, while South Asia is also the likely source of other Asian populations. The current northward invasion of B. dorsalis into Central China is the result of multiple, repeated dispersal events, most likely related to fruit trade. Results are discussed in the context of global quarantine, trade, and management of this pest. The recent expansion of the fly into temperate China, with very few associated genetic changes, clearly demonstrates the threat posed by this pest to ecologically similar areas in Europe and North America.

Genetic origin and dispersal of the invasive soybean aphid inferred from population genetic analysis and approximate Bayesian computation

Biological invasion is considered among the most important global environmental problems. Knowledge of the source and dispersal routes of invasion could facilitate the eradication and control of invasive species. Soybean aphid, Aphis glycines, is among the most destructive soybean pests. For effective management of this pest, we conducted genetic analyses and approximate Bayesian computation analysis to determine the origins and dispersal of the aphid species, as well as the source of its invasion in the USA, using 8 microsatellite loci and the mitochondrial cytochrome c oxidase subunit I gene. We were able to identify a significant isolation by distance pattern and 3 genetic lineages in the microsatellite data but not in the mtDNA dataset. The genetic structure showed that the USA population had the closest relationship with those from Korea and Japan, indicating that the 2 latter populations might be the sources of the invasion to the USA. Both population genetic analyses and approximate Bayesian computation showed that the northeastern populations in China were the possible sources of the further spread of A. glycines to Indonesia. The dispersal history of this aphid can provide useful information for pest management strategies and can further help predict areas at risk of invasion.

Biases of STRUCTURE software when exploring introduction routes of invasive species

Population genetic methods are widely used to retrace the introduction routes of invasive species. The unsupervised Bayesian clustering algorithm implemented in STRUCTURE is amongst the most frequently used of these methods, but its ability to provide reliable information about introduction routes has never been assessed. We simulated microsatellite datasets to evaluate the extent to which the results provided by STRUCTURE were misleading for the inference of introduction routes. We focused on an invasion scenario involving one native and two independently introduced populations, because it is the sole scenario that can be rejected when obtaining a particular clustering with a STRUCTURE analysis at K = 2 (two clusters). Results were classified as "misleading" or "non-misleading". We investigated the influence of effective size, bottleneck severity and number of loci on the type and frequency of misleading results. We showed that misleading STRUCTURE results were obtained for 10% of all simulated datasets. Our results highlighted two categories of misleading output. The first occurs when the native population has a low level of diversity. In this case, the two introduced populations may be very similar, despite their independent introduction histories. The second category results from convergence issues in STRUCTURE for K = 2, with strong bottleneck severity and/or large numbers of loci resulting in high levels of differentiation between the three populations. Overall, the risk of being misled by STRUCTURE in the context of introduction routes inferences is moderate, but it is important to remain cautious when low genetic diversity or genuine multimodality between runs are involved.

Low genetic diversity but strong population structure reflects multiple introductions of western flower thrips (Thysanoptera: Thripidae) into China followed by human-mediated spread

Historical invasion scenarios based on observational records are usually incomplete and biased, but these can be supplemented by population genetic data. The western flower thrips (WFT), Frankliniella occidentalis, invaded China in the last 13 years and has rapidly become one of the most serious pests in the country. To assess whether this invasion involved a single event or multiple events, we examined patterns of genetic diversity and population structure of WFT across 12 Chinese populations and a native US population based on mitochondrial DNA and/or 18 microsatellite loci. The average allelic richness and haplotype diversity in Chinese populations were significantly lower than in a population from its native range. The distribution of mitochondrial haplotypes suggested multiple independent invasions of WFT into China, including two invasions into the Beijing region. Based on microsatellite data, two distinct clusters were identified, with both of them splitting further into two clusters; in the Beijing region, the microsatellite data also provided evidence for two introductions. Both the absence of isolation by distance and the fact that distant populations were similar genetically suggest patterns of WFT movement linked to human activities. Our study therefore suggests multiple introductions of WFT into China and human-assisted spread.

The genome of Diuraphis noxia, a global aphid pest of small grains

BACKGROUND

The Russian wheat aphid, Diuraphis noxia Kurdjumov, is one of the most important pests of small grains throughout the temperate regions of the world. This phytotoxic aphid causes severe systemic damage symptoms in wheat, barley, and other small grains as a direct result of the salivary proteins it injects into the plant while feeding.

RESULTS

We sequenced and de novo assembled the genome of D. noxia Biotype 2, the strain most virulent to resistance genes in wheat. The assembled genomic scaffolds span 393 MB, equivalent to 93% of its 421 MB genome, and contains 19,097 genes. D. noxia has the most AT-rich insect genome sequenced to date (70.9%), with a bimodal CpG(O/E) distribution and a complete set of methylation related genes. The D. noxia genome displays a widespread, extensive reduction in the number of genes per ortholog group, including defensive, detoxification, chemosensory, and sugar transporter groups in comparison to the Acyrthosiphon pisum genome, including a 65% reduction in chemoreceptor genes. Thirty of 34 known D. noxia salivary genes were found in this assembly. These genes exhibited less homology with those salivary genes commonly expressed in insect saliva, such as glucose dehydrogenase and trehalase, yet greater conservation among genes that are expressed in D. noxia saliva but not detected in the saliva of other insects. Genes involved in insecticide activity and endosymbiont-derived genes were also found, as well as genes involved in virus transmission, although D. noxia is not a viral vector.

CONCLUSIONS

This genome is the second sequenced aphid genome, and the first of a phytotoxic insect. D. noxia's reduced gene content of may reflect the influence of phytotoxic feeding in shaping the D. noxia genome, and in turn in broadening its host range. The presence of methylation-related genes, including cytosine methylation, is consistent with other parthenogenetic and polyphenic insects. The D. noxia genome will provide an important contrast to the A. pisum genome and advance functional and comparative genomics of insects and other organisms.

Change in Biotypic Diversity of Russian Wheat Aphid (Hemiptera: Aphididae) Populations in the United States

A key component of Russian wheat aphid, Diuraphis noxia (Kurdjumov), management has been through planting resistant wheat cultivars. A new biotype, RWA2, appeared in 2003 which caused widespread damage to wheat cultivars containing the Dn4 gene. Biotypic diversity in Russian wheat aphid populations has not been addressed since 2005 when RWA2 dominated the biotype complex. Our objectives were to determine the biotypic diversity in the Central Great Plains and Colorado Plateau at regional (2010, 2011, 2013) and local (2012) levels and detect the presence of new Russian wheat aphid biotypes. Regional and within-field aphid collections were screened against Russian wheat aphid-resistant wheat genotypes containing genes Dn3, Dn4, Dn6, Dn7, Dn9, CI2401; and resistant barley STARS 9301B. In 2010, all aphid collections from Texas were avirulent to the Dn4 resistance gene in wheat. Regional results revealed Dn4 avirulent RWA6 was widespread (55-84%) in populations infesting wheat in both regions. Biotypes RWA1, 2, and 3/7 were equally represented with percentages<20% each while RWA8 was rarely detected. Combining percentages of RWA1, 6, and 8 across regions to estimate avirulence to Dn4 gene revealed high percentages for both 2011 (64-80%) and 2013 (69-90%). In contrast, the biotype structure at the local level differed where biotype percentages varied up to ≥2-fold between fields. No new biotypes were detected; therefore, Dn7, CI2401, and STARS9301B remained resistant to all known Russian wheat aphid biotypes. This study documents a shift to Dn4 avirulent biotypes and serves as a valuable baseline for biotypic diversity in Russian wheat aphid populations prior to the deployment of new Russian wheat aphid-resistant wheat cultivars.

What we still don't know about invasion genetics

Publication of The Genetics of Colonizing Species in 1965 launched the field of invasion genetics and highlighted the value of biological invasions as natural ecological and evolutionary experiments. Here, we review the past 50 years of invasion genetics to assess what we have learned and what we still don't know, focusing on the genetic changes associated with invasive lineages and the evolutionary processes driving these changes. We also suggest potential studies to address still-unanswered questions. We now know, for example, that rapid adaptation of invaders is common and generally not limited by genetic variation. On the other hand, and contrary to prevailing opinion 50 years ago, the balance of evidence indicates that population bottlenecks and genetic drift typically have negative effects on invasion success, despite their potential to increase additive genetic variation and the frequency of peak shifts. Numerous unknowns remain, such as the sources of genetic variation, the role of so-called expansion load and the relative importance of propagule pressure vs. genetic diversity for successful establishment. While many such unknowns can be resolved by genomic studies, other questions may require manipulative experiments in model organisms. Such studies complement classical reciprocal transplant and field-based selection experiments, which are needed to link trait variation with components of fitness and population growth rates. We conclude by discussing the potential for studies of invasion genetics to reveal the limits to evolution and to stimulate the development of practical strategies to either minimize or maximize evolutionary responses to environmental change.

Biological invasions, climate change and genomics

The rate of biological invasions is expected to increase as the effects of climate change on biological communities become widespread. Climate change enhances habitat disturbance which facilitates the establishment of invasive species, which in turn provides opportunities for hybridization and introgression. These effects influence local biodiversity that can be tracked through genetic and genomic approaches. Metabarcoding and metagenomic approaches provide a way of monitoring some types of communities under climate change for the appearance of invasives. Introgression and hybridization can be followed by the analysis of entire genomes so that rapidly changing areas of the genome are identified and instances of genetic pollution monitored. Genomic markers enable accurate tracking of invasive species' geographic origin well beyond what was previously possible. New genomic tools are promoting fresh insights into classic questions about invading organisms under climate change, such as the role of genetic variation, local adaptation and climate pre-adaptation in successful invasions. These tools are providing managers with often more effective means to identify potential threats, improve surveillance and assess impacts on communities. We provide a framework for the application of genomic techniques within a management context and also indicate some important limitations in what can be achieved.