Search

Abstract Background The high invasiveness of phytophagous insects is related to their adaptability to various environments, that can be influenced by their associated microbial community. Microbial symbionts are known to play a key role in the biology, ecology, and evolution of phytophagous insects, but their abundance and diversity are suggested to be influenced by environmental stressors. In this work, using 16 S rRNA metabarcoding we aim to verify (1) if laboratory rearing affects microbial symbiont communities of Zeugodacus cucurbitae females, a cosmopolitan pest of cucurbitaceous crops (2) if temperature, diet quality, and antibiotic treatments affect microbial symbiont communities of both laboratory and wild populations, and (3) if changes in microbial symbiont communities due to temperature, diet and antibiotic affect longevity and fecundity of Z. cucurbitae. Results The results showed that microbial diversity, particularly the β-diversity was significantly affected by insect origin, temperature, diet quality, and antibiotic treatment. The alteration of gut microbial symbionts, specifically Enterobacteriaceae, was associated with low fecundity and longevity of Z. cucurbitae females feeding on optimal diet only. Fecundity reduction in antibiotic treated females was more pronounced when flies were fed on a poor diet without protein. Conclusions our study proves the relationship between gut microbiome and host fitness under thermal and diet fluctuation highlighting the importance of microbial community in the adaptation of Z. cucurbitae to environmental stress. Clinical trial number Not applicable.

Abstract The Mediterranean fruit fly (Ceratitis capitata) is a globally invasive species and an economically significant pest of fruit crops. Understanding the evolutionary history and local climatic adaptation of this species is crucial for developing effective pest management strategies. We conducted a comprehensive investigation using whole genome sequencing to explore (i) the invasion history of C. capitata with an emphasis on historical admixture and (ii) local climatic adaptation across African, European, Central, and South American populations of C. capitata. Our results suggest a stepwise colonization of C. capitata in Europe and Latin America in which Mediterranean and Central American populations share an ancestral lineage. Conversely, South American invasion history is more complex, and our results partly suggest an old secondary invasion into South America from Europe or a colonization of South America directly from Africa, followed by admixture with an European lineage. Throughout its invasive range, C. capitata is challenged with diverse climatic regimes. A genome wide association study identified a relationship between allele frequency changes and specific bioclimatic variables. Notably, we observed a significant allele frequency shift related to adaptation to cold stress (BIO6), highlighting the species’ ability to rapidly adapt to seasonal variations in colder climates.

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.

Intentionally preserved biological material in natural history collections represents a vast repository of biodiversity. Advances in laboratory and sequencing technologies have made these specimens increasingly accessible for genomic analyses, offering a window into the genetic past of species and often permitting access to information that can no longer be sampled in the wild. Due to their age, preparation and storage conditions, DNA retrieved from museum and herbarium specimens is often poor in yield, heavily fragmented and biochemically modified. This not only poses methodological challenges in recovering nucleotide sequences, but also makes such investigations susceptible to environmental and laboratory contamination. In this paper, we review the practical challenges associated with making the recovery of DNA sequence data from museum collections more routine. We first review key operational principles and issues to address, to guide the decision-making process and dialogue between researchers and curators about when and how to sample museum specimens for genomic analyses. We then outline the range of steps that can be taken to reduce the likelihood of contamination including laboratory set-ups, workflows and working practices. We finish by presenting a series of case studies, each focusing on protocol practicalities for the application of different mainstream methodologies to museum specimens including: (i) shotgun sequencing of insect mitogenomes, (ii) whole genome sequencing of insects, (iii) genome skimming to recover plant plastid genomes from herbarium specimens, (iv) target capture of multi-locus nuclear sequences from herbarium specimens, (v) RAD-sequencing of bird specimens and (vi) shotgun sequencing of ancient bovid bone samples.

Gut microbial communities are critical in determining the evolutive success of fruit fly phytophagous pests (Diptera, Tephritidae), facilitating their adaptation to suboptimal environmental conditions and to plant allelochemical defences. An important source of variation for the microbial diversity of fruit flies is represented by the crop on which larvae are feeding. However, a "crop effect" is not always the main driver of microbial patterns, and it is often observed in combination with other and less obvious processes. In this work, we aim at verifying if environmental stress and, by extension, changing environmental conditions, can promote microbial diversity in Zeugodacus cucurbitae (Coquillett), a cosmopolitan pest of cucurbit crops. With this objective, 16S rRNA metabarcoding was used to test differences in the microbial profiles of wild fly populations in a large experimental setup in Eastern Central Tanzania. The analysis of 2,973 unique ASV, which were assigned to 22 bacterial phyla, 221 families and 590 putative genera, show that microbial α diversity (as estimated by Abundance Coverage Estimator, Faith's Phylogenetic Diversity, Shannon-Weiner and the Inverse Simpson indexes) as well as β microbial diversity (as estimated by Compositional Data analysis of ASVs and of aggregated genera) significantly change as the species gets closer to its altitudinal limits, in farms where pesticides and agrochemicals are used. Most importantly, the multivariate dispersion of microbial patterns is significantly higher in these stressful environmental conditions thus indicating that Anna Karenina effects contribute to the microbial diversity of Z. cucurbitae. The crop effect was comparably weaker and detected as non-consistent changes across the experimental sites. We speculate that the impressive adaptive potential of polyphagous fruit flies is, at least in part, related to the Anna Karenina principle, which promotes stochastic changes in the microbial diversity of fly populations exposed to suboptimal environmental conditions.

Abstract An increasing number of invasive fruit fly pests are colonizing new grounds. With this study, we aimed to uncover the invasion pathways of the oriental fruit fly, Bactrocera dorsalis into the islands of the Indian Ocean. By using genome‐wide SNP data and a multipronged approach consisting of PCA, ancestry analysis, phylogenetic inference, and kinship networks, we were able to resolve two independent invasion pathways. A western invasion pathway involved the stepping‐stone migration of B. dorsalis from the east African coast into the Comoros, along Mayotte and into Madagascar with a decreasing genetic diversity. The Mascarene islands (Reunion and Mauritius), on the contrary, were colonized directly from Asia and formed a distinct cluster. The low nucleotide diversity suggests that only a few genotypes invaded the Mascarenes. The presence of many long runs of homozygosity (ROH) in the introduced populations is indicative of population bottlenecks, with evidence of a more severe bottleneck for populations along the western migration pathway than on the Mascarene islands. More strict phytosanitary regulations are recommended in order to prevent the further spread of B. dorsalis.

Abstract Background Gastropod snails remain strongly understudied, despite their important role in transmitting parasitic diseases. Knowledge of their distribution and population dynamics increases our understanding of the processes driving disease transmission. We report the first study to use high-throughput sequencing (HTS) to elucidate the population genetic structure of the hermaphroditic snail Bulinus truncatus (Gastropoda, Heterobranchia) on a regional (17–150 km) and inter-regional (1000–5400 km) scale. This snail species acts as an intermediate host of Schistosoma haematobium and Schistosoma bovis, which cause human and animal schistosomiasis respectively. Methods Bulinus truncatus snails were collected in Senegal, Cameroon, Egypt and France and identified through DNA barcoding. A single-end genotyping-by-sequencing (GBS) library, comprising 87 snail specimens from the respective countries, was built and sequenced on an Illumina HiSeq 2000 platform. Reads were mapped against S. bovis and S. haematobium reference genomes to identify schistosome infections, and single nucleotide polymorphisms (SNPs) were scored using the Stacks pipeline. These SNPs were used to estimate genetic diversity, assess population structure and construct phylogenetic trees of B. truncatus. Results A total of 10,750 SNPs were scored and used in downstream analyses. The phylogenetic analysis identified five clades, each consisting of snails from a single country but with two distinct clades within Senegal. Genetic diversity was low in all populations, reflecting high selfing rates, but varied between locations due to habitat variability. Significant genetic differentiation and isolation by distance patterns were observed at both spatial scales, indicating that gene flow is not strong enough to counteract the effects of population bottlenecks, high selfing rates and genetic drift. Remarkably, the population genetic differentiation on a regional scale (i.e. within Senegal) was as large as that between populations on an inter-regional scale. The blind GBS technique was able to pick up parasite DNA in snail tissue, demonstrating the potential of HTS techniques to further elucidate the role of snail species in parasite transmission. Conclusions HTS techniques offer a valuable toolbox to further investigate the population genetic patterns of intermediate schistosome host snails and the role of snail species in parasite transmission. Graphical Abstract

Abstract While contemporary changes in feeding preferences have been documented in phytophagous insects, the mechanisms behind these processes remain to be fully clarified. In this context, the insect gut microbiome plays a central role in adaptation to novel host plants. The cucurbit frugivorous fruit fly Zeugodacus cucurbitae (Diptera, Tephritidae) has occasionally been reported on “unconventional” host plants from different families, including Solanaceae. In this study, we focus on wild parental (F0) adults and semiwild first filial (F1) larvae of Z. cucurbitae from multiple sites in La Réunion and explore how the gut microbiome composition changes when this fly is feeding on a noncucurbit host (Solanum melongena). Our analyses show nonobvious gut microbiome responses following the F0–F1 host shift and the importance of not just diet but also local effects, which heavily affected the diversity and composition of microbiomes. We identified the main bacterial genera responsible for differences between treatments. These data further stress the importance of a careful approach when drawing general conclusions based on laboratory populations or inadequately replicated field samples.

Abstract Background Clupeid fisheries in Lake Tanganyika (East Africa) provide food for millions of people in one of the world’s poorest regions. Due to climate change and overfishing, the clupeid stocks of Lake Tanganyika are declining. We investigate the population structure of the Lake Tanganyika sprat Stolothrissa tanganicae, using for the first time a genomic approach on this species. This is an important step towards knowing if the species should be managed separately or as a single stock. Population structure is important for fisheries management, yet understudied for many African freshwater species. We hypothesize that distinct stocks of S. tanganicae could be present due to the large size of the lake (isolation by distance), limnological variation (adaptive evolution), or past separation of the lake (historical subdivision). On the other hand, high mobility of the species and lack of obvious migration barriers might have resulted in a homogenous population. Results We performed a population genetic study on wild-caught S. tanganicae through a combination of mitochondrial genotyping (96 individuals) and RAD sequencing (83 individuals). Samples were collected at five locations along a north-south axis of Lake Tanganyika. The mtDNA data had low global FST and, visualised in a haplotype network, did not show phylogeographic structure. RAD sequencing yielded a panel of 3504 SNPs, with low genetic differentiation (FST = 0.0054; 95% CI: 0.0046–0.0066). PCoA, fineRADstructure and global FST suggest a near-panmictic population. Two distinct groups are apparent in these analyses (FST = 0.1338 95% CI: 0.1239,0.1445), which do not correspond to sampling locations. Autocorrelation analysis showed a slight increase in genetic difference with increasing distance. No outlier loci were detected in the RADseq data. Conclusion Our results show at most very weak geographical structuring of the stock and do not provide evidence for genetic adaptation to historical or environmental differences over a north-south axis. Based on these results, we advise to manage the stock as one population, integrating one management strategy over the four riparian countries. These results are a first comprehensive study on the population structure of these important fisheries target species, and can guide fisheries management.

The family Tephritidae includes some of the most notorious insect pests of agricultural and horticultural crops in tropical and sub-tropical regions. Despite the interest in the study of their gut microbiome, our present knowledge is largely based on the analysis of laboratory strains. In this study, we present a first comparative analysis of the gut microbiome profiles of field populations of ten African and Mediterranean tephritid pests. For each species, third instar larvae were sampled from different locations and host fruits and compared using 16S rRNA amplicon sequencing and a multi-factorial sampling design. We observed considerable variation in gut microbiome diversity and composition both between and within fruit fly species. A “core” microbiome, shared across all targeted species, could only be identified at most at family level (Enterobacteriaceae). At genus level only a few bacterial genera (Klebsiella, Enterobacter, and Bacillus) were present in most, but not all, samples, with high variability in their relative abundance. Higher relative abundances were found for seven bacterial genera in five of the fruit fly species considered. These were Erwinia in Bactrocera oleae, Lactococcus in B. zonata, Providencia in Ceratitis flexuosa, Klebsiella, and Rahnella in C. podocarpi and Acetobacter and Serratia in C. rosa. With the possible exception of C. capitata and B. dorsalis (the two most polyphagous species considered) we could not detect obvious relationships between fruit fly dietary breadth and microbiome diversity or abundance patterns. Similarly, our results did not suggest straightforward differences between the microbiome profiles of species belonging to Ceratitis and the closely related Bactrocera/Zeugodacus. These results provide a first comparative analysis of the gut microbiomes of field populations of multiple economically relevant tephritids and provide base line information for future studies that will further investigate the possible functional role of the observed associations.

Abstract The PREDICTS project—Projecting Responses of Ecological Diversity In Changing Terrestrial Systems (https://www.nhm.ac.uk/our-science/our-work/biodiversity/predicts.html)—has collated from published studies a large, reasonably representative database of comparable samples of biodiversity from multiple sites that differ in the nature or intensity of human impacts relating to land use. We have used this evidence base to develop global and regional statistical models of how local biodiversity responds to these measures. We describe and make freely available this 2016 release of the database, containing more than 3.2 million records sampled at over 26,000 locations and representing over 47,000 species. We outline how the database can help in answering a range of questions in ecology and conservation biology. To our knowledge, this is the largest and most geographically and taxonomically representative database of spatial comparisons of biodiversity that has been collated to date; it will be useful to researchers and international efforts wishing to model and understand the global status of biodiversity.

High-throughput sequencing (HTS) techniques are of great value for the investigation of microbial communities, and have been extensively used to study the gut microbiome. While most studies focus on the human gut, many others have investigated insects. However, because of the rapid spread of HTS techniques, a lot of variation exists in the protocols for sample preparation. In the present study, we investigated the impact of two widely adopted sample-processing procedures preceding library preparation, i.e., preservation of insect tissue in 70% ethanol (EtOH) and sample dissection. We used the fruit fly Ceratitis capitata (Diptera: Tephritidae) as a model organism and set up two experiments, one comparing the effects of sample manipulation and preservation across life stages and the other across fruit samples from different sources. The results of this study showed no major effects of dissection on the outcome of HTS. However, EtOH preservation did have effects on the recovered gut microbiome, the main effect being a significant reduction of the dominant genus, Providencia, in EtOH-preserved samples. Less abundant bacterial groups were also affected resulting in altered microbial profiles obtained from samples preserved in 70% EtOH. These results have important implications for the planning of future studies and when comparing studies that used different sample preparation protocols.

Here, we present novel high quality genome assemblies for five invasive tephritid species of agricultural concern: Ceratitis capitata, C. quilicii, C. rosa, Zeugodacus cucurbitae and Bactrocera zonata (read depths between 65 and 78x). Three assemblies ( C. capitata, C. quilicii and Z. cucurbitae) were scaffolded with chromosome conformation data and annotated using RNAseq reads. For some species this is the first reference genome available ( B. zonata, C. quilicii and C. rosa), for others we have published improved annotated genomes ( C. capitata and Z. cucurbitae). Together, the new references provide an important resource to advance research on genetic techniques for population control, develop rapid species identification methods, and explore eco-evolutionary studies.

Identification by DNA barcoding is more likely to be erroneous when it is based on a large distance between the query (the barcode sequence of the specimen to identify) and its best match in a reference barcode library. The number of such false positive identifications can be decreased by setting a distance threshold above which identification has to be rejected. To this end, we proposed recently to use an ad hoc distance threshold producing identifications with an estimated relative error probability that can be fixed by the user (e.g. 5%). Here we introduce two R functions that automate the calculation of ad hoc distance thresholds for reference libraries of DNA barcodes. The scripts of both functions, a user manual and an example file are available on the JEMU website (http://jemu.myspecies.info/computer-programs) as well as on the comprehensive R archive network (CRAN, http://cran.r-project.org).

This paper describes five new Ceratitis species from the eastern and southern parts of the Afrotropical Region: C. (Pterandrus) quilicii De Meyer, Mwatawala & Virgilio sp. nov.; C. (Ceratalaspis) pallidula De Meyer, Mwatawala & Virgilio sp. nov.; C. (Ceratalaspis) taitaensis De Meyer & Copeland sp. nov.; C. (Ceratalaspis) sawahilensis De Meyer & Virgilio sp. nov.; and C. (Ceratalaspis) flavipennata De Meyer & Virgilio sp. nov. Their relationships with closely allied species within their respective subgenera are discussed where appropriate, and diagnostic characters are given. DNA barcodes are provided for all new species. In addition, the hitherto unknown male of C. (Pardalaspis) serrata De Meyer, 1996 is described, based on material collected in the Democratic Republic of Congo. Recognition of these new species and sexes is the result of an integrative approach using morphological characters and DNA data.

The identification of Afrotropical hoverflies is very difficult because of limited recent taxonomic revisions and the lack of comprehensive identification keys. In order to assist in their identification, and to improve the taxonomy of this group, we constructed a reference dataset of 513 COI barcodes of 90 of the more common nominal species from Ghana, Togo, Benin and Nigeria (W Africa) and added ten publically available COI barcodes from nine nominal Afrotropical species to this (total: 523 COI barcodes; 98 nominal species; 26 genera). The identification accuracy of this dataset was evaluated with three methods (K2P distance-based, Neighbor-Joining (NJ) / Maximum Likelihood (ML) analysis, and using SpeciesIdentifier). Results of the three methods were highly congruent and showed a high identification success. Nine species pairs showed a low (< 0.03) mean interspecific K2P distance that resulted in several incorrect identifications. A high (> 0.03) maximum intraspecific K2P distance was observed in eight species and barcodes of these species not always formed single clusters in the NJ / ML analayses which may indicate the occurrence of cryptic species. Optimal K2P thresholds to differentiate intra- from interspecific K2P divergence were highly different among the three subfamilies (Eristalinae: 0.037, Syrphinae: 0.06, Microdontinae: 0.007-0.02), and among the different general suggesting that optimal thresholds are better defined at the genus level. In addition to providing an alternative identification tool, our study indicates that DNA barcoding improves the taxonomy of Afrotropical hoverflies by selecting (groups of) taxa that deserve further taxonomic study, and by attributing the unknown sex to species for which only one of the sexes is known.

We propose a general working strategy to deal with incomplete reference libraries in the DNA barcoding identification of species. Considering that (1) queries with a large genetic distance with their best DNA barcode match are more likely to be misidentified and (2) imposing a distance threshold profitably reduces identification errors, we modelled relationships between identification performances and distance thresholds in four DNA barcode libraries of Diptera (n = 4270), Lepidoptera (n = 7577), Hymenoptera (n = 2067) and Tephritidae (n = 602 DNA barcodes). In all cases, more restrictive distance thresholds produced a gradual increase in the proportion of true negatives, a gradual decrease of false positives and more abrupt variations in the proportions of true positives and false negatives. More restrictive distance thresholds improved precision, yet negatively affected accuracy due to the higher proportions of queries discarded (viz. having a distance query-best match above the threshold). Using a simple linear regression we calculated an ad hoc distance threshold for the tephritid library producing an estimated relative identification error