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

snpAIMeR: R package for evaluating ancestry informative marker contributions in non-model population diagnostics

MOTIVATION

Single nucleotide polymorphism (SNP) markers are increasingly popular for population genomics and inferring ancestry for individuals of unknown origin. Because large SNP datasets are impractical for rapid and routine analysis, diagnostics rely on panels of highly informative markers. Strategies exist for selecting these markers, however, resources for efficiently evaluating their performance are limited for non-model systems.

RESULTS

snpAIMeR is a user-friendly R package that evaluates the efficacy of genomic markers for the cluster assignment of unknown individuals. It is intended to help minimize panel size and genotyping effort by determining the informativeness of candidate diagnostic markers. Provided genotype data from individuals of known origin, it uses leave-one-out cross-validation to determine population assignment rates for individual markers and marker combinations.

AVAILABILITY AND IMPLEMENTATION

snpAIMeR is available on CRAN (https://CRAN.R-project.org/package=snpAIMeR).

Genomes of the cosmopolitan fruit pest Bactrocera dorsalis (Diptera: Tephritidae) reveal its global invasion history and thermal adaptation

INTRODUCTION

The oriental fruit fly Bactrocera dorsalis is one of the most destructive agricultural pests worldwide, with highly debated species delimitation, origin, and global spread routes.

OBJECTIVES

Our study intended to (i) resolve the taxonomic uncertainties between B. dorsalis and B. carambolae, (ii) reveal the population structure and global invasion routes of B. dorsalis across Asia, Africa, and Oceania, and (iii) identify genomic regions that are responsible for the thermal adaptation of B. dorsalis.

METHODS

Based on a high-quality chromosome-level reference genome assembly, we explored the population relationship using a genome-scale single nucleotide polymorphism dataset generated from the resequencing data of 487 B. dorsalis genomes and 25 B. carambolae genomes. Genome-wide association studies and silencing using RNA interference were used to identify and verify the candidate genes associated with extreme thermal stress.

RESULTS

We showed that B. dorsalis originates from the Southern India region with three independent invasion and spread routes worldwide: (i) from Northern India to Northern Southeast Asia, then to Southern Southeast Asia; (ii) from Northern India to Northern Southeast Asian, then to China and Hawaii; and (iii) from Southern India toward the African mainland, then to Madagascar, which is mainly facilitated by human activities including trade and immigration. Twenty-seven genes were identified by a genome-wide association study to be associated with 11 temperature bioclimatic variables. The Cyp6a9 gene may enhance the thermal adaptation of B. dorsalis and thus boost its invasion, which tended to be upregulated at a hardening temperature of 38 °C. Functional verification using RNA interference silencing against Cyp6a9, led to the specific decrease in Cyp6a9 expression, reducing the survival rate of dsRNA-feeding larvae exposed to extreme thermal stress of 45 °C after heat hardening treatments in B. dorsalis.

CONCLUSION

This study provides insights into the evolutionary history and genetic basis of temperature adaptation in B. dorsalis.

From chip to SNP: Rapid development and evaluation of a targeted capture genotyping-by-sequencing approach to support research and management of a plaguing rodent

The management of invasive species has been greatly enhanced by population genetic analyses of multilocus single-nucleotide polymorphism (SNP) datasets that provide critical information regarding pest population structure, invasion pathways, and reproductive biology. For many applications there is a need for protocols that offer rapid, robust and efficient genotyping on the order of hundreds to thousands of SNPs, that can be tailored to specific study populations and that are scalable for long-term monitoring schemes. Despite its status as a model laboratory species, there are few existing resources for studying wild populations of house mice (Mus musculus spp.) that strike this balance between data density and laboratory efficiency. Here we evaluate the utility of a custom targeted capture genotyping-by-sequencing approach to support research on plaguing house mouse populations in Australia. This approach utilizes 3,651 hybridization capture probes targeting genome-wide SNPs identified from a sample of mice collected in grain-producing regions of southeastern Australia genotyped using a commercially available microarray platform. To assess performance of the custom panel, we genotyped wild caught mice (N = 320) from two adjoining farms and demonstrate the ability to correctly assign individuals to source populations with high confidence (mean >95%), as well as robust kinship inference within sites. We discuss these results in the context of proposed applications for future genetic monitoring of house mice in Australia.

Divergent east-west lineages in an Australian fruit fly, (Bactrocera jarvisi), associated with the Carpentaria Basin divide

Bactrocera jarvisi is an endemic Australian fruit fly species (Diptera: Tephritidae). It occurs commonly across tropical and subtropical coastal Australia, from far-northern Western Australia, across the 'Top End' of the Northern Territory, and then down the Queensland east coast. Across this range, its distribution crosses several well documented biogeographic barriers. In order to better understand factors leading to the divergence of Australian fruit fly lineages, we carried out a population genetic study of B. jarvisi from across its range using genome-wide SNP analysis, utilising adult specimens gained from trapping and fruit rearing. Populations from the Northern Territory (NT) and Western Australia were genetically similar to each other, but divergent from the genetically uniform east-coast (= Queensland, QLD) population. Phylogenetic analysis demonstrated that the NT population derived from the QLD population. We infer a role for the Carpentaria Basin as a biogeographic barrier restricting east-west gene flow. The QLD populations were largely panmictic and recognised east-coast biogeographic barriers play no part in north-south population structuring. While the NT and QLD populations were genetically distinct, there was evidence for the historically recent translocation of flies from each region to the other. Flies reared from different host fruits collected in the same location showed no genetic divergence. While a role for the Carpentaria Basin as a barrier to gene flow for Australian fruit flies agrees with existing work on the related B. tryoni, the reason(s) for population panmixia for B. jarvisi (and B. tryoni) over the entire Queensland east coast, a linear north-south distance of >2000km, remains unknown.

Nuclear genetic diversity and structure of Anastrepha ludens wild populations evidenced by microsatellite markers

The Mexican fruit fly, Anastrepha ludens, is an important pest that causes widespread damage to a number of fruit crops in Mexico. The sterile insect technique (SIT) is commonly used for its control. However, the existence of natural barriers can give rise to a population structure in neutral loci and possibly behavioral or adaptive traits that interfere with SIT. For this reason, it is important to understand the genetic diversity and structure of A. ludens populations and to better understand the evolutionary ecology and population processes in view of possible expansions and possible host shifts due to climate change. We genotyped nine nuclear DNA (nDNA) microsatellite loci among fruit fly populations collected from five biogeographic areas within Mexico, namely, the Mexican Plateau, the Northeastern Coastal Plain, the Pacific Coast, the Gulf Coast of Mexico, and the Soconusco, and a laboratory strain. The nuclear genetic diversity was moderate (from He = 0.34 to He = 0.39) within the wild mexfly population. We found that populations were clustered in three genetic groups (K = 3). The diversity and the genetic structure of A. ludens are determined by environmental and geological conditions, as well as local conditions like anthropogenic perturbation, which would produce population expansion and the existence of possible predators that would affect the population density. Gene flow showed recent migration among populations. The laboratory strain showed fewer diversity than the wild samples. Large values of current and ancestral population size suggest high resistance to climatic changes, probably due to biological attributes, such as its polyphagous, multivoltine, and high dispersal characteristics. In particular, ecosystem fragmentation and perturbation as well as the existence of new plant hosts would probably increase the abundance of flies.

Invasive potential of tropical fruit flies in temperate regions under climate change

Tropical fruit flies are considered among the most economically important invasive species detected in temperate areas of the United States and the European Union. Detections often trigger quarantine and eradication programs that are conducted without a holistic understanding of the threat posed. Weather-driven physiologically-based demographic models are used to estimate the geographic range, relative abundance, and threat posed by four tropical tephritid fruit flies (Mediterranean fruit fly, melon fly, oriental fruit fly, and Mexican fruit fly) in North and Central America, and the European-Mediterranean region under extant and climate change weather (RCP8.5 and A1B scenarios). Most temperate areas under tropical fruit fly propagule pressure have not been suitable for establishment, but suitability is predicted to increase in some areas with climate change. To meet this ongoing challenge, investments are needed to collect sound biological data to develop mechanistic models to predict the geographic range and relative abundance of these and other invasive species, and to put eradication policies on a scientific basis.

The Sterile Insect Technique: Success and Perspectives in the Neotropics

The sterile insect technique (SIT), an environmentally friendly means of control, is currently used against plant, animal, and human pests under the area-wide integrated pest management. It consists in the mass production, sterilization, and release of insects in an affected area where sterile males mate with wild females leading to no reproduction. Here, we review SIT in the Neotropics and focus on particular recent successful cases of eradication of the Mediterranean fruit fly, Ceratitis capitata (Wiedemann), as well as effective programs used against the Mexican fruit fly Anastrepha ludens (Loew), the New World screwworm fly Cochliomyia hominivorax (Coquerel)), and the Cactus moth Cactoblastis cactorum (Berg). We examine when SIT does not work and innovations that have made SIT more efficient and also highlight complimentary techniques that can be used in conjunction. We address potential candidate species that could be controlled through SIT, for example Philornis downsi Dodge & Aitken. Finally, we consider the impact of climate change in the context of the use of the SIT against these pests. Given the recent dramatic decline in insect biodiversity, investing in environmentally friendly means of pest control should be a priority. We conclude that SIT should be promoted in the region, and leadership and political will is needed for continued success of SIT in the Neotropics.

Disjunction between canola distribution and the genetic structure of its recently described pest, the canola flower midge (Contarinia brassicola)

Population genomics is a useful tool to support integrated pest management as it can elucidate population dynamics, demography, and histories of invasion. Here, we use a restriction site-associated DNA sequencing approach combined with whole-genome amplification (WGA) to assess genomic population structure of a newly described pest of canola, the diminutive canola flower midge, Contarinia brassicola. Clustering analyses recovered little geographic structure across the main canola production region but differentiated several geographically disparate populations at edges of the agricultural zone. Given a lack of alternative hypotheses for this pattern, we suggest these data support alternative hosts for this species and thus our canola-centric view of this midge as a pest has limited our understanding of its biology. These results speak to the need for increased surveying efforts across multiple habitats and other potential hosts within Brassicaceae to improve both our ecological and evolutionary knowledge of this species and contribute to effective management strategies. We additionally found that use of WGA prior to library preparation was an effective method for increasing DNA quantity of these small insects prior to restriction site-associated DNA sequencing and had no discernible impact on genotyping consistency for population genetic analysis; WGA is therefore likely to be tractable for other similar studies that seek to randomly sample markers across the genome in small organisms.

Genetic and Ecological Relationships of Anastrepha ludens (Diptera: Tephritidae) Populations in Southern Mexico

Knowledge of the influence of evolutionary factors that promote either the differentiation or cohesion of pest insect populations is critical for the improvement of control strategies. Here, we explore the extent to which genetic differentiation occurs between populations of the Mexican fruit fly, Anastrepha ludens, in association with four plant hosts (Citrus sinensis, C. paradisi, Mangifera indica and Casimiroa edulis) in the Soconusco region of Chiapas (Mexico). Using variants from six enzymatic loci, we obtained measures of genetic diversity for three sample arrangements: (1) by sex per locality, (2) by locality and (3) by host. The extent of genetic differentiation in populations was assessed using the Analyses of Molecular Variance (AMOVA) method for each array of samples, and moderate to high levels of genetic variation were observed between the sexes, as well as among localities and host plants. A Bayesian approach was then used to assess any population structure underlying the genetic data we obtained, but this analysis showed no significant structuring due to locality or host plant. We also considered whether the observed genotypic frequencies in male and females matched those expected under a hypothesis of random mating. Here we found significant deviations from expected genotypic frequencies, suggesting that sexual selection is acting on these populations. Overall, our results indicate that sexual selection, along with the presence of some heterogeneity in environments provided by both geographical factors and availability of host plants, has influenced the evolution of pest populations in this region of Mexico. Implications for area-wide pest management strategies are discussed.

From Nucleotides to Satellite Imagery: Approaches to Identify and Manage the Invasive Pathogen Xylella fastidiosa and Its Insect Vectors in Europe

Biological invasions represent some of the most severe threats to local communities and ecosystems. Among invasive species, the vector-borne pathogen Xylella fastidiosa is responsible for a wide variety of plant diseases and has profound environmental, social and economic impacts. Once restricted to the Americas, it has recently invaded Europe, where multiple dramatic outbreaks have highlighted critical challenges for its management. Here, we review the most recent advances on the identification, distribution and management of X. fastidiosa and its insect vectors in Europe through genetic and spatial ecology methodologies. We underline the most important theoretical and technological gaps that remain to be bridged. Challenges and future research directions are discussed in the light of improving our understanding of this invasive species, its vectors and host–pathogen interactions. We highlight the need of including different, complimentary outlooks in integrated frameworks to substantially improve our knowledge on invasive processes and optimize resources allocation. We provide an overview of genetic, spatial ecology and integrated approaches that will aid successful and sustainable management of one of the most dangerous threats to European agriculture and ecosystems.