Host plant-related genomic differentiation in the European cherry fruit fly, Rhagoletis cerasi

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.

Genomic Traces of the Fruit Fly Anastrepha obliqua Associated with Its Polyphagous Nature

Anastrepha obliqua (Macquart) (Diptera: Tephritidae) is an important pest in the neotropical region. It is considered a polyphagous insect, meaning it infests plants of different taxonomic families and readily colonizes new host plants. The change to new hosts can lead to diversification and the formation of host races. Previous studies investigating the effect of host plants on population structure and selection in Anastrepha obliqua have focused on the use of data from the mitochondrial DNA sequence and microsatellite markers of nuclear DNA, and there are no analyses at the genomic level. To better understand this issue, we used a pooled restriction site-associated DNA sequencing (pooled RAD-seq) approach to assess genomic differentiation and population structure across sympatric populations of Anastrepha obliqua that infest three host plants-Spondias purpurea (red mombin), Mangifera indica (mango) of the family Anacardiaceae and Averrhoa carambola (carambola) of the family Oxalidaceae-in sympatric populations of the species Anastrepha obliqua of Inter-Andean Valley of the Cauca River in southwestern Colombia. Our results show genomic differentiation of populations from carambola compared to mango and red mombin populations, but the genetic structure was mainly established by geography rather than by the host plant. On the other hand, we identified 54 SNPs in 23 sequences significantly associated with the use of the host plant. Of these 23 sequences, we identified 17 candidate genes and nine protein families, of which four protein families are involved in the nutrition of these flies. Future studies should investigate the adaptive processes undergone by phytophagous insects in the Neotropics, using fruit flies as a model and state-of-the-art molecular tools.

Genome analyses of four Wolbachia strains and associated mitochondria of Rhagoletis cerasi expose cumulative modularity of cytoplasmic incompatibility factors and cytoplasmic hitchhiking across host populations

Background

The endosymbiont Wolbachia can manipulate arthropod reproduction and invade host populations by inducing cytoplasmic incompatibility (CI). Some host species are coinfected with multiple Wolbachia strains which may have sequentially invaded host populations by expressing different types of modular CI factor (cif) genes. The tephritid fruit fly Rhagoletis cerasi is a model for CI and Wolbachia population dynamics. It is associated with at least four Wolbachia strains in various combinations, with demonstrated (wCer2, wCer4), predicted (wCer1) or unknown (wCer5) CI phenotypes.

Results

We sequenced and assembled the draft genomes of the Wolbachia strains wCer1, wCer4 and wCer5, and compared these with the previously sequenced genome of wCer2 which currently invades R. cerasi populations. We found complete cif gene pairs in all strains: four pairs in wCer2 (three Type I; one Type V), two pairs in wCer1 (both Type I) and wCer4 (one Type I; one Type V), and one pair in wCer5 (Type IV). Wolbachia genome variant analyses across geographically and genetically distant host populations revealed the largest diversity of single nucleotide polymorphisms (SNPs) in wCer5, followed by wCer1 and then wCer2, indicative of their different lengths of host associations. Furthermore, mitogenome analyses of the Wolbachia genome-sequenced individuals in combination with SNP data from six European countries revealed polymorphic mitogenome sites that displayed reduced diversity in individuals infected with wCer2 compared to those without.

Conclusions

Coinfections with Wolbachia are common in arthropods and affect options for Wolbachia-based management strategies of pest and vector species already infected by Wolbachia. Our analyses of Wolbachia genomes of a host naturally coinfected by several strains unravelled signatures of the evolutionary dynamics in both Wolbachia and host mitochondrial genomes as a consequence of repeated invasions. Invasion of already infected populations by new Wolbachia strains requires new sets of functionally different cif genes and thereby may select for a cumulative modularity of cif gene diversity in invading strains. Furthermore, we demonstrated at the mitogenomic scale that repeated CI-driven Wolbachia invasions of hosts result in reduced mitochondrial diversity and hitchhiking effects. Already resident Wolbachia strains may experience similar cytoplasmic hitchhiking effects caused by the invading Wolbachia strain.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07906-6.

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.

Wolbachia-Mitochondrial DNA Associations in Transitional Populations of Rhagoletis cerasi

The endosymbiont Wolbachia can manipulate arthropod host reproduction by inducing cytoplasmic incompatibility (CI), which results in embryonic mortality when infected males mate with uninfected females. A CI-driven invasion of Wolbachia can result in a selective sweep of associated mitochondrial haplotype. The co-inheritance of Wolbachia and host mitochondrial DNA can therefore provide significant information on the dynamics of an ongoing Wolbachia invasion. Therefore, transition zones (i.e., regions where a Wolbachia strain is currently spreading from infected to uninfected populations) represent an ideal area to investigate the relationship between Wolbachia and host mitochondrial haplotype. Here, we studied Wolbachia-mitochondrial haplotype associations in the European cherry fruit fly, Rhagoletis cerasi, in two transition zones in the Czech Republic and Hungary, where the CI-inducing strain wCer2 is currently spreading. The wCer2-infection status of 881 individuals was compared with the two known R. cerasi mitochondrial haplotypes, HT1 and HT2. In accordance with previous studies, wCer2-uninfected individuals were associated with HT1, and wCer2-infected individuals were mainly associated with HT2. We found misassociations only within the transition zones, where HT2 flies were wCer2-uninfected, suggesting the occurrence of imperfect maternal transmission. We did not find any HT1 flies that were wCer2-infected, suggesting that Wolbachia was not acquired horizontally. Our study provides new insights into the dynamics of the early phase of a Wolbachia invasion.

Species complex diversification by host plant use in an herbivorous insect: The source of Puerto Rican cactus mealybug pest and implications for biological control

Cryptic taxa have often been observed in the form of host‐associated species that diverged as the result of adaptation to alternate host plants. Untangling cryptic diversity in species complexes that encompass invasive species is a mandatory task for pest management. Moreover, investigating the evolutionary history of a species complex may help to understand the drivers of their diversification. The mealybug Hypogeococcus pungens was believed to be a polyphagous species from South America and has been reported as a pest devastating native cacti in Puerto Rico, also threatening cactus diversity in the Caribbean and North America. There is neither certainty about the identity of the pest nor the source population from South America. Recent studies pointed to substantial genetic differentiation among local populations, suggesting that H. pungens is a species complex. In this study, we used a combination of genome‐wide SNPs and mtDNA variation to investigate species diversity within H. pungens sensu lato to establish host plant ranges of each one of the putative members of the complex, to evaluate whether the pattern of host plant association drove diversification in the species complex, and to determine the source population of the Puerto Rican cactus pest. Our results suggested that H. pungens comprises at least five different species, each one strongly associated with specific host plants. We also established that the Puerto Rican cactus pest derives from southeastern Brazilian mealybugs. This is an important achievement because it will help to design reliable strategies for biological control using natural enemies of the pest from its native range.

Host plant-related genomic differentiation in the European cherry fruit fly, Rhagoletis cerasi

Elucidating the mechanisms and conditions facilitating the formation of biodiversity are central topics in evolutionary biology. A growing number of studies imply that divergent ecological selection may often play a critical role in speciation by counteracting the homogenising effects of gene flow. Several examples involve phytophagous insects, where divergent selection pressures associated with host plant shifts may generate reproductive isolation, promoting speciation. Here, we use ddRADseq to assess the population structure and to test for host‐related genomic differentiation in the European cherry fruit fly, Rhagoletis cerasi (L., 1758) (Diptera: Tephritidae). This tephritid is distributed throughout Europe and western Asia, and has adapted to two different genera of host plants, Prunus spp. (cherries) and Lonicera spp. (honeysuckle). Our data imply that geographic distance and geomorphic barriers serve as the primary factors shaping genetic population structure across the species range. Locally, however, flies genetically cluster according to host plant, with consistent allele frequency differences displayed by a subset of loci between Prunus and Lonicera flies across four sites surveyed in Germany and Norway. These 17 loci display significantly higher F_ST values between host plants than others. They also showed high levels of linkage disequilibrium within and between Prunus and Lonicera flies, supporting host‐related selection and reduced gene flow. Our findings support the existence of sympatric host races in R. cerasi embedded within broader patterns of geographic variation in the fly, similar to the related apple maggot, Rhagoletis pomonella, in North America.