Does Host Plant Drive Variation in Microbial Gut Communities in a Recently Shifted Pest?

Biotic interactions can modulate the responses of organisms to environmental stresses, including diet changes. Gut microbes have substantial effects on diverse ecological and evolutionary traits of their hosts, and microbial communities can be highly dynamic within and between individuals in space and time. Modulations of the gut microbiome composition and their potential role in the success of a species to maintain itself in a new environment have been poorly studied to date. Here we examine this question in a large wood-boring beetle Cacosceles newmannii (Cerambycidae), that was recently found thriving on a newly colonized host plant. Using 16S metabarcoding, we assessed the gut bacterial community composition of larvae collected in an infested field and in "common garden" conditions, fed under laboratory-controlled conditions on four either suspected or known hosts (sugarcane, tea tree, wattle, and eucalyptus). We analysed microbiome variation (i.e. diversity and differentiation), measured fitness-related larval growth, and studied host plant lignin and cellulose contents, since their degradation is especially challenging for wood-boring insects. We show that sugarcane seems to be a much more favourable host for larval growth. Bacterial diversity level was the highest in field-collected larvae, whereas lab-reared larvae fed on sugarcane showed a relatively low level of diversity but very specific bacterial variants. Bacterial communities were mainly dominated by Proteobacteria, but were significantly different between sugarcane-fed lab-reared larvae and any other hosts or field-collected larvae. We identified changes in the gut microbiome associated with different hosts over a short time frame, which support the hypothesis of a role of the microbiome in host switches.

A young age of subspecific divergence in the desert locust inferred by ABC random forest.

Dating population divergence within species from molecular data and relating such dating to climatic and biogeographic changes is not trivial. Yet it can help formulating evolutionary hypotheses regarding local adaptation and future responses to changing environments. Key issues include statistical selection of a demographic and historical scenario among a set of possible scenarios, and estimation of the parameter(s) of interest under the chosen scenario. Such inferences greatly benefit from (a) independent information on evolutionary rate and pattern at genetic markers; and (b) new statistical approaches, such as approximate Bayesian computation-random forest (ABC-RF), which provides reliable inference at a low computational cost and the possibility to measure prediction quality at the exact position of the observed data set. Here, we show full potential of the ABC-RF approach including prior knowledge on microsatellite genetic markers to decipher the evolutionary history of the African arid-adapted pest locust, Schistocerca gregaria, with support for a southern colonization of Africa, from a low number of founders of northern origin, dating back 2.6 Ky (90% CI: 0.9-6.6 Ky). We verify that this divergence time estimate accurately reflected true divergence time values by computing accuracy at a local posterior scale from simulated pseudo-observed data sets. The inferred divergence history is better explained by the peculiar biology of S. gregaria, which involves a density-dependent swarming phase with some exceptional spectacular migrations, rather than a continuous colonization resulting from the continental expansion of open vegetation habitats during more ancient Quaternary glacial climatic episodes.

On the relative role of climate change and management in the current desert locust outbreak in East Africa.

While large-scale monitoring, early detection and control can greatly reduce desert locust invasions, global change is most likely to affect conditions that promote the transition from solitary to gregarious populations. Although climate change scenarios point to an increase in aridity and further desertification in vast areas of Africa, some regions that have been at the origin of past outbreaks are likely to see a reversed trend (i.e., increase in frequency and intensity of rains), potentially favoring the formation of swarms. This makes reinforcing early detection and keeping a sustained monitoring effort in place even more important under climate change.

Climate-driven geographic distribution of the desert locust during recession periods: Subspecies' niche differentiation and relative risks under scenarios of climate change.

The desert locust is an agricultural pest that is able to switch from a harmless solitarious stage, during recession periods, to swarms of gregarious individuals that disperse long distances and affect areas from western Africa to India during outbreak periods. Large outbreaks have been recorded through centuries, and the Food and Agriculture Organization keeps a long-term, large-scale monitoring survey database in the area. However, there is also a much less known subspecies that occupies a limited area in Southern Africa. We used large-scale climatic and occurrence data of the solitarious phase of each subspecies during recession periods to understand whether both subspecies climatic niches differ from each other, what is the current potential geographical distribution of each subspecies, and how climate change is likely to shift their potential distribution with respect to current conditions. We evaluated whether subspecies are significantly specialized along available climate gradients by using null models of background climatic differences within and between southern and northern ranges and applying niche similarity and niche equivalency tests. The results point to climatic niche conservatism between the two clades. We complemented this analysis with species distribution modeling to characterize current solitarious distributions and forecast potential recession range shifts under two extreme climate change scenarios at the 2050 and 2090 time horizon. Projections suggest that, at a global scale, the northern clade could contract its solitarious recession range, while the southern clade is likely to expand its recession range. However, local expansions were also predicted in the northern clade, in particular in southern and northern margins of the current geographical distribution. In conclusion, monitoring and management practices should remain in place in northern Africa, while in Southern Africa the potential for the subspecies to pose a threat in the future should be investigated more closely.

Spatial heterogeneity in landscape structure influences dispersal and genetic structure: empirical evidence from a grasshopper in an agricultural landscape.

Dispersal may be strongly influenced by landscape and habitat characteristics that could either enhance or restrict movements of organisms. Therefore, spatial heterogeneity in landscape structure could influence gene flow and the spatial structure of populations. In the past decades, agricultural intensification has led to the reduction in grassland surfaces, their fragmentation and intensification. As these changes are not homogeneously distributed in landscapes, they have resulted in spatial heterogeneity with generally less intensified hedged farmland areas remaining alongside streams and rivers. In this study, we assessed spatial pattern of abundance and population genetic structure of a flightless grasshopper species, Pezotettix giornae, based on the surveys of 363 grasslands in a 430-km² agricultural landscape of western France. Data were analysed using geostatistics and landscape genetics based on microsatellites markers and computer simulations. Results suggested that small-scale intense dispersal allows this species to survive in intensive agricultural landscapes. A complex spatial genetic structure related to landscape and habitat characteristics was also detected. Two P. giornae genetic clusters bisected by a linear hedged farmland were inferred from clustering analyses. This linear hedged farmland was characterized by high hedgerow and grassland density as well as higher grassland temporal stability that were suspected to slow down dispersal. Computer simulations demonstrated that a linear-shaped landscape feature limiting dispersal could be detected as a barrier to gene flow and generate the observed genetic pattern. This study illustrates the relevance of using computer simulations to test hypotheses in landscape genetics studies.

Demographic processes shaping genetic variation of the solitarious phase of the desert locust.

Between plagues, the solitarious desert locust (Schistocerca gregaria) is generally thought to exist as small populations, which are particularly prone to extinction events in arid regions of Africa and Asia. Given the high genetic structuring observed in one geographical area (the Eritrean coast) by former authors, a metapopulation dynamics model involving repeated extinction and colonization events was favoured. In this study, we assessed the validity of a demographic scenario involving temporary populations of the solitarious phase of the desert locust by analysing large-scale population genetic data. We scored 24 microsatellites in 23 solitarious population samples collected over most of the species range during remission. We found very little genetic structuring and little evidence of declining genetic diversity. A Bayesian clustering method distinguished four genetically differentiated units. Three groups were largely consistent with three population samples which had undergone recent bottleneck events. Nevertheless, the last genetically homogeneous unit included all individuals from the remaining 18 population samples and did not show evidence of demographic disequilibrium. An approximate Bayesian computation treatment indicated a large population size for this main genetic group, moderately reduced between plague and remission but still containing tens of thousands of individuals. Our results diverge from the hypothesis of a classical metapopulation dynamics model. They instead support the scenario in which large populations persist in the solitarious phase of the desert locust.

Evaluation of 96-well high-throughput DNA extraction methods for 16S rRNA gene metabarcoding.

Gaining meaningful insights into bacterial communities associated with animal hosts requires the provision of high-quality nucleic acids. Although many studies have compared DNA extraction methods for samples with low bacterial biomass (e.g. water) or specific PCR inhibitors (e.g. plants), DNA extraction bias in samples without inherent technical constraint (e.g. animal samples) has received little attention. Furthermore, there is an urgent need to identify a DNA extraction methods in a high-throughput format that decreases the cost and time for processing large numbers of samples. We here evaluated five DNA extraction protocols, using silica membrane-based spin columns and a 96-well microplate format and based on either mechanical or enzymatic lysis or a combination of both, using three bacterial mock communities and Illumina sequencing of the V4 region of the 16SrRNA gene. Our results showed that none of the DNA extraction methods fully eliminated bias associated with unequal lysis efficiencies. However, we identified a DNA extraction method with a lower bias for each mock community standard. Of these methods, those including an enzymatic lysis showed biases specific to some bacteria. Altogether, these results again demonstrate the importance of DNA extraction standardization to be able to compare the microbiome results of different samples. In this attempt, we advise for the use of the 96-well DNeasy Blood and Tissue kit (Qiagen) with a zirconia bead-beating procedure, which optimizes altogether the cost, handling time and bacteria-specific effects associated with enzymatic lysis.

Domestication provides durum wheat with protection from locust herbivory.

Lower plant resistance to herbivores following domestication has been suggested as the main cause for higher feeding damage in crops than in wild progenitors. While herbivore compensatory feeding has also been proposed as a possible mechanism for raised damage in crops with low nutritional quality, predictions regarding the effects of plant domestication on nutritional quality for herbivores remain unclear. In particular, data on primary metabolites, even major macronutrients, measured in the organs consumed by herbivores, are scarce. In this study, we used a collection of 10 accessions of wild ancestors and 10 accessions of modern progenies of Triticum turgidum to examine whether feeding damage and selectivity by nymphs of Locusta migratoria primarily depended on five leaf traits related to structural resistance or nutrient profiles. Our results unexpectedly showed that locusts favored wild ancestors over domesticated accessions and that leaf toughness and nitrogen and soluble protein contents increased with the domestication process. Furthermore, the quantitative relationship between soluble protein and digestible carbohydrates was found to poorly meet the specific requirements of the herbivore, in all wheat accessions, both wild and modern. The increase in leaf structural resistance to herbivores in domesticated tetraploid wheat accessions suggested that resource allocation trade-offs between growth and herbivory resistance may have been disrupted by domestication in the vegetative organs of this species. Since domestication did not result in a loss of nutritional quality in the leaves of the tetraploid wheat, our results rather provides evidence for a role of the content of plants in nonnutritive nitrogenous secondary compounds, possibly deterrent or toxic, at least for grasshopper herbivores.

Evaluation of identification methods for cryptic Bactrocera dorsalis (Diptera: Tephritidae) specimens: combining morphological and molecular techniques.

The potential for population genomics to elucidate invasion pathways of a species is limited by taxonomic identification issues. The Oriental fruit fly pest, Bactrocera dorsalis (Hendel) belongs to a complex in which several sympatric species are attracted to the same lure used in trapping and are morphologically cryptic and/or reported to hybridize. In this study, we evaluated the taxonomic ambiguity between B. dorsalis and 2 major cryptic species, based on morphological expertise and 289 target specimens sampled across the whole distribution range. Specimens were then subjected to DNA sequence analyses of the COI mitochondrial barcode and the EIF3L nuclear marker to evaluate the potential for molecular identification, in particular for specimens for which morphological identification was inconclusive. To this aim, we produced reference datasets with DNA sequences from target specimens whose morphological identification was unambiguous, which we complemented with 56 new DNA sequences from closest relatives and 76 published and curated DNA sequences of different species in the complex. After the necessary morphological observation, about 3.5% of the target dataset and 47.6% of the specimens from Southeast Asian islands displayed ambiguous character states shared with B. carambolae and/or B. occipitalis. Critical interpretation of DNA sequence data solved morphological ambiguities only when combining both mitochondrial and nuclear markers. COI discriminated B. dorsalis from 5 species; EIF3L and ITS from another species. We recommend this procedure to ensure correct identification of B. dorsalis specimens in population genetics studies and surveillance programs.

Mitochondrial genomes reveal the global phylogeography and dispersal routes of the migratory locust.

The migratory locust, Locusta migratoria, is the most widely distributed grasshopper species in the world. However, its global genetic structure and phylogeographic relationships have not been investigated. In this study, we explored the worldwide genetic structure and phylogeography of the locust populations based on the sequence information of 65 complete mitochondrial genomes and three mitochondrial genes of 263 individuals from 53 sampling sites. Although this locust can migrate over long distances, our results revealed high genetic differentiation among the geographic populations. The populations can be divided into two different lineages: the Northern lineage, which includes individuals from the temperate regions of the Eurasian continent, and the Southern lineage, which includes individuals from Africa, southern Europe, the Arabian region, India, southern China, South-east Asia and Australia. An analysis of population genetic diversity indicated that the locust species originated from Africa. Ancestral populations likely separated into Northern and Southern lineages 895 000 years ago by vicariance events associated with Pleistocene glaciations. These two lineages evolved in allopatry and occupied their current distributions in the world via distinct southern and northern dispersal routes. Genetic differences, caused by the long-term independent diversification of the two lineages, along with other factors, such as geographic barriers and temperature limitations, may play important roles in maintaining the present phylogeographic patterns. Our phylogeographic evidence challenged the long-held view of multiple subspecies in the locust species and tentatively divided it into two subspecies, L. m. migratoria and L. m. migratorioides.

Long microsatellites and unusually high levels of genetic diversity in the Orthoptera.

Much remains to be learned about the mutational processes governing the evolution of microsatellite repeat regions and the associated levels of genetic diversity observed at microsatellite markers across populations or species. An extensive survey of microsatellite variation in 210 insect species from six major orders revealed that within Orthopterans, which are characterized by giant genomes, levels of genetic diversity were ~20% higher and microsatellite repeat arrays were longer than in any other group. Because of the mutation dependence on repeat length, this result suggests a higher microsatellite loci mutation rate in the Orthoptera. We deem it plausible that differences among insect orders, either in mismatch repair systems or in abundance of transposable element-derived microsatellites, can shape the size distribution of both genomes and microsatellite repeat regions. Our findings emphasise that observed levels of genetic diversity can greatly vary across species (orders at least) because of molecular differences in the mechanisms that determine microsatellite size, and are therefore critical to conservation and population genetics studies, where microsatellite repeat variability is primarily interpreted in terms of population demography and history.

Isolation and characterization of twelve polymorphic microsatellite Loci for the cocoa mirid bug Sahlbergella singularis.

Mirids are the primary pests affecting cocoa production in Africa, but no genetic studies have been conducted on these insects. Here we report the isolation and characterization of 12 polymorphic microsatellite loci for Sahlbergella singularis. A microsatellite-enriched genomic DNA library was developed and screened to identify marker loci. Twelve polymorphic loci were identified by screening 28 individuals collected from one presumed population in cocoa plantations in Southern Cameroon. The number of alleles ranged from 5 to 25, whereas the observed and the expected heterozygosities ranged from 0.179 to 0.786 and from 0.671 to 0.946, respectively. Tests showed significant deviations from HW equilibrium for four loci, but no linkage disequilibrium was detected at any of the loci. No cross-species amplification was observed in two other mirid pests in Africa.

Microsatellite markers for the Chameleon grasshopper (Kosciuscola tristis) (Orthoptera: Acrididae), an Australian Alpine Specialist.

A set of polymorphic loci was characterised using an enrichment library for the Australian alpine specialist, the chameleon grasshopper (Kosciuscola tristis), an atypical grasshopper known for its remarkable temperature-controlled colour change. The number of alleles per locus ranged from three to 20 and observed heterozygosity from 0.16 to 0.76. These are the first microsatellite markers for a non-endangered Australian alpine animal and will inform questions of gene flow across the sky islands of this unique and threatened region.

Faeces' odours attract gregarious locust hoppers.

Collective motion is one of the most impressive common features of gregarious locusts: once formed, bands and swarms get moving for long distances. It was shown that visual perception of neighbours plays a key role in maintaining marching behaviour at a local scale. But at a larger scale, mechanisms underlying band cohesion are less understood. It was shown in several field studies that individuals separated from the band were able to get back to the group, even after being separated since a night. In this context, faeces' odours could be a possible indicator of the recent passage of a group. In this study, we tested if nymphs are attracted by faeces' odours and if this effect is modulated by the age of the faeces. To this end, we conducted individual olfactometric behavioural assays of 3rd instar hoppers of desert locust, Schistocerca gregaria, exposed to odours of 1 h-old and 24 h-old faeces. We also used Gas Chromatography-Mass Spectrometry (GC-MS) to identify odours' volatile organic compounds from faeces. The results of behavioural assays indicated a strong attractive effect of faeces, with no preference for one of the two faecal age classes. Nymphs spent significantly more time in the side of the olfactometer where the faeces' odours came from, and 72.7% of tested individuals chose this side first. We filtered and annotated 11 volatile organic compounds present in both fresh and old faeces in GC-MS analyses, including guaiacol and phenol, which are known to cause an aggregative effect on desert locusts. As the attractive effect lasted over 24 h, band's faeces could still have an attractive effect when individuals are separated from the band since one day. In this situation, latecomers individuals would be able to get back to the group by following the traces of their predecessors.

Assessment and validation of a suite of reverse transcription-quantitative PCR reference genes for analyses of density-dependent behavioural plasticity in the Australian plague locust.

BACKGROUND: The Australian plague locust, Chortoicetes terminifera, is among the most promising species to unravel the suites of genes underling the density-dependent shift from shy and cryptic solitarious behaviour to the highly active and aggregating gregarious behaviour that is characteristic of locusts. This is because it lacks many of the major phenotypic changes in colour and morphology that accompany phase change in other locust species. Reverse transcription-quantitative polymerase chain reaction (RT-qPCR) is the most sensitive method available for determining changes in gene expression. However, to accurately monitor the expression of target genes, it is essential to select an appropriate normalization strategy to control for non-specific variation between samples. Here we identify eight potential reference genes and examine their expression stability at different rearing density treatments in neural tissue of the Australian plague locust. RESULTS: Taking advantage of the new orthologous DNA sequences available in locusts, we developed primers for genes encoding 18SrRNA, ribosomal protein L32 (RpL32), armadillo (Arm), actin 5C (Actin), succinate dehydrogenase (SDHa), glyceraldehyde-3P-dehydrogenase (GAPDH), elongation factor 1 alpha (EF1a) and annexin IX (AnnIX). The relative transcription levels of these eight genes were then analyzed in three treatment groups differing in rearing density (isolated, short- and long-term crowded), each made up of five pools of four neural tissue samples from 5th instar nymphs. SDHa and GAPDH, which are both involved in metabolic pathways, were identified as the least stable in expression levels, challenging their usefulness in normalization. Based on calculations performed with the geNorm and NormFinder programs, the best combination of two genes for normalization of gene expression data following crowding in the Australian plague locust was EF1a and Arm. We applied their use to studying a target gene that encodes a Ca2+ binding glycoprotein, SPARC, which was previously found to be up-regulated in brains of gregarious desert locusts, Schistocerca gregaria. Interestingly, expression of this gene did not vary with rearing density in the same way in brains of the two locust species. Unlike S. gregaria, there was no effect of any crowding treatment in the Australian plague locust. CONCLUSION: Arm and EF1a is the most stably expressed combination of two reference genes of the eight examined for reliable normalization of RT-qPCR assays studying density-dependent behavioural change in the Australian plague locust. Such normalization allowed us to show that C. terminifera crowding did not change the neuronal expression of the SPARC gene, a gregarious phase-specific gene identified in brains of the desert locust, S. gregaria. Such comparative results on density-dependent gene regulation provide insights into the evolution of gregarious behaviour and mass migration of locusts. The eight identified genes we evaluated are also candidates as normalization genes for use in experiments involving other Oedipodinae species, but the rank order of gene stability must necessarily be determined on a case-by-case basis.

Challenges to assessing connectivity between massive populations of the Australian plague locust.

Linking demographic and genetic dispersal measures is of fundamental importance for movement ecology and evolution. However, such integration can be difficult, particularly for highly fecund species that are often the target of management decisions guided by an understanding of population movement. Here, we present an example of how the influence of large population sizes can preclude genetic approaches from assessing demographic population structuring, even at a continental scale. The Australian plague locust, Chortoicetes terminifera, is a significant pest, with populations on the eastern and western sides of Australia having been monitored and managed independently to date. We used microsatellites to assess genetic variation in 12 C. terminifera population samples separated by up to 3000 km. Traditional summary statistics indicated high levels of genetic diversity and a surprising lack of population structure across the entire range. An approximate Bayesian computation treatment indicated that levels of genetic diversity in C. terminifera corresponded to effective population sizes conservatively composed of tens of thousands to several million individuals. We used these estimates and computer simulations to estimate the minimum rate of dispersal, m, that could account for the observed range-wide genetic homogeneity. The rate of dispersal between both sides of the Australian continent could be several orders of magnitude lower than that typically considered as required for the demographic connectivity of populations.

Additive genetic variance for traits least related to fitness increases with environmental stress in the desert locust, Schistocerca gregaria.

Under environmental stress, previously hidden additive genetic variation can be unmasked and exposed to selection. The amount of hidden variation is expected to be higher for life history traits, which strongly correlate to individual fitness, than for morphological traits, in which fitness effects are more ambiguous. However, no consensual pattern has been recovered yet, and this idea is still debated in the literature. Here, we hypothesize that the classical categorization of traits (i.e., life history and morphology) may fail to capture their proximity to fitness. In the desert locust, Schistocerca gregaria, a model organism for the study of insect polyphenism, we quantified changes in additive genetic variation elicited by lifetime thermal stress for ten traits, in which evolutionary significance is known. Irrespective of their category, traits under strong stabilizing selection showed genetic invariance with environmental stress, while traits more loosely associated with fitness showed a marked increase in additive genetic variation in the stressful environment. Furthermore, traits involved in adaptive phenotypic plasticity (growth compensation) showed either no change in additive genetic variance or a change of moderate magnitude across thermal environments. We interpret this mitigated response of plastic traits in the context of integrated evolution to adjust the entire phenotype in heterogeneous environments (i.e., adaptiveness of initial plasticity, compromise of phenotypic compensation with stress, and shared developmental pathway). Altogether, our results indicate, in agreement with theoretical expectations, that environmental stress can increase available additive genetic variance in some desert locust traits, but those closely linked to fitness are largely unaffected. Our study also highlights the importance of assessing the proximity to fitness of a trait on a case-by-case basis and in an ecologically relevant context, as well as considering the processes of canalization and plasticity, involved in the control of phenotypic variation.

Integrative taxonomy confirms that Gregarina garnhami and G. acridiorum (Apicomplexa, Gregarinidae), parasites of Schistocerca gregaria and Locusta migratoria (Insecta, Orthoptera), are distinct species.

Orthoptera are infected by about 60 species of gregarines assigned to the genus Gregarina Dufour, 1828. Among these species, Gregarina garnhami Canning, 1956 from Schistocerca gregaria (Forsskål, 1775) was considered by Lipa et al. in 1996 to be synonymous with Gregarina acridiorum (Léger 1893), a parasite of several orthopteran species including Locusta migratoria (Linné, 1758). Here, a morphological study and molecular analyses of the SSU rDNA marker demonstrate that specimens of S. gregaria and specimens of L. migratoria are infected by two distinct Gregarina species, G. garnhami and G. acridiorum, respectively. Validation of the species confirms that molecular analyses provide useful taxonomical information. Phenotypic plasticity was clearly observed in the case of G. garnhami: the morphology of its trophozoites, gamonts and syzygies varied according to the geographical location of S. gregaria and the subspecies infected.

TITLE: La taxonomie intégrative confirme que Gregarina garnhami et G. acridiorum (Apicomplexa, Gregarinidae), parasites de Schistocerca gregaria et Locusta migratoria (Insecta, Orthoptera), sont des espèces distinctes. ABSTRACT: Les orthoptères sont parasités par environ soixante espèces de grégarines affiliées au genre Gregarina Dufour, 1828. Parmi ces espèces Gregarina garnhami Canning, 1956 décrite chez Schistocerca gregaria (Forskål, 1775), a été mise en synonymie par Lipa et al. en 1996 avec Gregarina acridiorum (Léger 1893), parasite de plusieurs espèces d'orthoptères dont Locusta migratoria (Linné, 1758). Ici, une étude morphologique et des analyses moléculaires du marqueur SSU rDNA démontrent que les spécimens de S. gregaria et ceux de L. migratoria sont infectés par 2 espèces distinctes de grégarines, Gregarina garnhami et Gregarina acridiorum, respectivement. La validation de ces espèces confirme l'importance des informations fournies par les analyses moléculaires dans les études taxonomiques. Une plasticité phénotypique a été clairement observée dans le cas de G. garnhami : la morphologie de ses trophozoïtes, gamontes et syzygies varie selon la localisation géographique et la sous-espèce de S. gregaria infectée.

Outbreaks, gene flow and effective population size in the migratory locust, Locusta migratoria: a regional-scale comparative survey.

The potential effect of population outbreaks on within and between genetic variation of populations in pest species has rarely been assessed. In this study, we compare patterns of genetic variation in different sets of historically frequently outbreaking and rarely outbreaking populations of an agricultural pest of major importance, the migratory locust, Locusta migratoria. We analyse genetic variation within and between 24 populations at 14 microsatellites in Western Europe, where only ancient and low-intensity outbreaks have been reported (non-outbreaking populations), and in Madagascar and Northern China, where frequent and intense outbreak events have been recorded over the last century (outbreaking populations). Our comparative survey shows that (i) the long-term effective population size is similar in outbreaking and non-outbreaking populations, as evidenced by similar estimates of genetic diversity, and (ii) gene flow is substantially larger among outbreaking populations than among non-outbreaking populations, as evidenced by a fourfold to 30-fold difference in FST values. We discuss the implications for population dynamics and the consequences for management strategies of the observed patterns of genetic variation in L. migratoria populations with contrasting historical outbreak frequency and extent.

Deciphering host-parasitoid interactions and parasitism rates of crop pests using DNA metabarcoding.

An accurate estimation of parasitism rates and diversity of parasitoids of crop insect pests is a prerequisite for exploring processes leading to efficient natural biocontrol. Traditional methods such as rearing have been often limited by taxonomic identification, insect mortality and intensive work, but the advent of high-throughput sequencing (HTS) techniques, such as DNA metabarcoding, is increasingly seen as a reliable and powerful alternative approach. Little has been done to explore the benefits of such an approach for estimating parasitism rates and parasitoid diversity in an agricultural context. In this study, we compared the composition of parasitoid species and parasitism rates between rearing and DNA metabarcoding of host eggs and larvae of the millet head miner, Heliocheilus albipunctella De Joannis (Lepidoptera, Noctuidae), collected from millet fields in Senegal. We first assessed the detection threshold for the main ten endoparasitoids, by sequencing PCR products obtained from artificial dilution gradients of the parasitoid DNAs in the host moth. We then assessed the potential of DNA metabarcoding for diagnosing parasitism rates in samples collected from the field. Under controlled conditions, our results showed that relatively small quantities of parasitoid DNA (0.07 ng) were successfully detected within an eight-fold larger quantity of host DNA. Parasitoid diversity and parasitism rate estimates were always higher for DNA metabarcoding than for host rearing. Furthermore, metabarcoding detected multi-parasitism, cryptic parasitoid species and differences in parasitism rates between two different sampling sites. Metabarcoding shows promise for gaining a clearer understanding of the importance and complexity of host-parasitoid interactions in agro-ecosystems, with a view to improving pest biocontrol strategies.