Search

The growing urbanization process is accompanied by the emergence of new habitats for wildlife, and cities are sometimes seen as refuges for pollinators such as wild bees compared to intensively cultivated rural habitats. However, the contrasting living conditions that combine high fragmentation, exposure to pollutants, and heat island effects, with low pesticide use and potentially high availability of resources, make it difficult to predict the overall effect of urban living on the health of wild bees. Moreover, if the responses of wild bee populations in terms of species richness and diversity have been the focus of many recent studies, individual responses to urbanization have been more rarely investigated. More specifically, data on the impacts on individual bee physiology and health are lacking. To help fill this gap, we collected red-tailed bumblebee ( Bombus lapidarius ) workers along a gradient of urbanization defined by the level of soil imperviousness, and estimated environmental (air quality) and ecological (pathogens' prevalence and loads; local competition) pressures they locally experienced. In parallel, we quantified the expression of selected immune marker genes. We measured how the immune system of bumblebees responds to urbanization gradient and which local parameters best explain the observed changes in immune gene expression. We evidenced three immune markers, tightly linked with cellular metabolism, whose expressions increase with the level of urbanization, independently of individual infection and pollution exposure. We suggest that induction of their expression reveals a shift in wild bee immunometabolism, supposedly in response to the stressful conditions experienced in areas with high built-up cover. The induction of these genes is likely at the root of any immune activation; they could thus be used as markers to estimate the levels of urban stress locally experienced by pollinators., Competing Interests: The authors declare no conflicts of interest., (© 2024 The Author(s). Ecology and Evolution published by John Wiley & Sons Ltd.)

Invasive vectors can induce dramatic changes in disease epidemiology. While viral emergence following geographical range expansion of a vector is well known, the influence a vector can have at the level of the host's pathobiome is less well understood. Taking advantage of the formerly heterogeneous spatial distribution of the ectoparasitic mite Varroa destructor that acts as potent virus vector among honeybees Apis mellifera , we investigated the impact of its recent global spread on the viral community of honeybees in a retrospective study of historical samples. We hypothesized that the vector has had an effect on the epidemiology of several bee viruses, potentially altering their transmissibility and/or virulence, and consequently their prevalence, abundance, or both. To test this, we quantified the prevalence and loads of 14 viruses from honeybee samples collected in mite-free and mite-infested populations in four independent geographical regions. The presence of the mite dramatically increased the prevalence and load of deformed wing virus, a cause of unsustainably high colony losses. In addition, several other viruses became more prevalent or were found at higher load in mite-infested areas, including viruses not known to be actively varroa-transmitted, but which may increase opportunistically in varroa-parasitized bees., Competing Interests: The authors declare no conflicts of interest., (© 2024 The Authors.)

The environmental factors affecting plant reproduction and effective dispersal, in particular biotic interactions, have a strong influence on plant expansion dynamics, but their demographic and genetic consequences remain an understudied body of theory. Here, we use a mathematical model in a one-dimensional space and on a single reproductive period to describe the joint effects of predispersal seed insect predators foraging strategy and plant reproduction strategy (masting) on the spatio-temporal dynamics of seed sources diversity in the colonisation front of expanding plant populations. We show that certain foraging strategies can result in a higher seed predation rate at the colonisation front compared to the core of the population, leading to an Allee effect. This effect promotes the contribution of seed sources from the core to the colonisation front, with long-distance dispersal further increasing this contribution. As a consequence, our study reveals a novel impact of the predispersal seed predation-induced Allee effect, which mitigates the erosion of diversity in expanding populations. We use rearrangement inequalities to show that masting has a buffering role: it mitigates this seed predation-induced Allee effect. This study shows that predispersal seed predation, plant reproductive strategies and seed dispersal patterns can be intermingled drivers of the diversity of seed sources in expanding plant populations, and opens new perspectives concerning the analysis of more complex models such as integro-difference or reaction-diffusion equations., (© 2023. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

Insects are under pressure from agricultural intensification. To protect pollinators, conservation measures such as the EU agri-environment schemes (AES) promote planting wildflowers along fields. However, this can potentially alter disease ecology by serving as transmission hubs or by diluting infections. We tested this by measuring plant-pollinator interactions and virus infections (DWV-A, DWV-B and ABPV) across pollinator communities in agricultural landscapes over a year. AES had a direct effect on DWV-B, reducing prevalence and load in honeybees, with a tentative general dilution effect on load in early summer. DWV-A prevalence was reduced both under AES and with increasing niche overlap between competent hosts, likely via a dilution effect. By contrast, AES had no impact on ABPV, its prevalence driven by the proportion of bumblebees in the community. Epidemiological differences were also reflected in the virus phylogenies, with DWV-B showing recent rapid expansion, while DWV-A and ABPV showed slower growth rates and geographical population structure. Phylogenies indicate that all three viruses freely circulate across their host populations. Our study illustrates how complex interactions between environmental, ecological and evolutionary factors may influence wildlife disease dynamics. Supporting pollinator nutrition can mitigate the transmission of important bee diseases, providing an unexpected boost to pollinator conservation. This article is part of the theme issue 'Infectious disease ecology and evolution in a changing world'.

Genetic variation for resistance and disease tolerance has been described in a range of species. In Drosophila melanogaster , genetic variation in mortality following systemic Drosophila C virus (DCV) infection is driven by large-effect polymorphisms in the restriction factor pastrel (pst) . However, it is unclear if pst contributes to disease tolerance. We investigated systemic DCV challenges spanning nine orders of magnitude, in males and females of 10 Drosophila Genetic Reference Panel lines carrying either a susceptible (S) or resistant (R) pst allele. We find among-line variation in fly survival, viral load and disease tolerance measured both as the ability to maintain survival (mortality tolerance) and reproduction (fecundity tolerance). We further uncover novel effects of pst on host vigour, as flies carrying the R allele exhibited higher survival and fecundity even in the absence of infection. Finally, we found significant genetic variation in the expression of the JAK-STAT ligand upd3 and the epigenetic regulator of JAK-STAT G9a. However, while G9a has been previously shown to mediate tolerance of DCV infection, we found no correlation between the expression of either upd3 or G9a on fly tolerance or resistance. Our work highlights the importance of both resistance and tolerance in viral defence.

Ecological restoration programs are established to reverse land degradation, mitigate biodiversity loss, and reinstate ecosystem services. Following recent agricultural intensification that led to a decrease in flower diversity and density in rural areas and subsequently to the decline of many insects, conservation measures targeted at pollinators have been established, including sown wildflower strips (WFS) along field margins. Historically successful in establishing a high density of generalist bees and increasing pollinator diversity, the impact of enhanced flower provision on wider ecological interactions and the structure of pollinator networks has been rarely investigated. Here, we tested the effects of increasing flower species richness and flower density in agricultural landscapes on bee-plant interaction networks. We measured plant species richness and flower density and surveyed honeybee and bumblebee visits on flowers across a range of field margins on 10 UK farms that applied different pollinator conservation measures. We found that both flower species richness and flower density significantly increased bee abundance, in early and late summer, respectively. At the network level, we found that higher flower species richness did not significantly alter bee species' generality indices, but significantly reduced network connectance and marginally reduced niche overlap across honeybees and bumblebee species, a proxy for insect competition. While higher connectance and niche overlap is believed to strengthen network robustness and often is the aim for the restoration of pollinator networks, we argue that carefully designed WFS may benefit bees by partitioning their foraging niche, limiting competition for resources and the potential for disease transmission via shared floral use. We also discuss the need to extend WFS and their positive effects into spring when wild bee populations are established., Competing Interests: None declared., (© 2022 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd.)

The OneHealth approach aims to further our understanding of the drivers of human, animal and environmental health, and, ultimately, to improve them by combining approaches and knowledge from medicine, biology and fields beyond. Wild and managed bees are essential pollinators of crops and wild flowers. Their health thus directly impacts on human and environmental health. At the same time, these bee species represent highly amenable and relevant model organisms for a OneHealth approach that aims to study fundamental epidemiological questions. In this review, we focus on how infectious diseases of wild and managed bees can be used as a OneHealth model system, informing fundamental questions on ecological immunology and disease transmission, while addressing how this knowledge can be used to tackle the issues facing pollinator health., (Copyright © 2020. Published by Elsevier Inc.)

The ratio of amino acids to carbohydrates (AA:C) that bumble bees consume has been reported to affect their survival. However, it is unknown how dietary AA:C ratio affects other bumble bee fitness traits (e.g., fecundity, condition) and possible trade-offs between them. Moreover, while individual AAs affect phenotype in many species, the effects of AA blend on bumble bee fitness and food intake are unclear. We test how the AA:C ratio that bumble bees (Bombus terrestris) consume affects their condition (abdomen lipid and dry mass), survival following food removal, and ovarian activation. We then compare ovarian activation and food intake in bees fed identical AA:C ratios, but where the blend of AAs in diets differ, i.e., diets contained the same 10 AAs in an equimolar ratio or in the same ratio as in bee collected pollen. We found that AA:C ratio did not significantly affect survival following food removal or ovarian activation; however, high AA intake increased body mass, which is positively correlated with multiple fitness traits in bumble bees. AA blend (i.e., equimolar versus pollen) did not significantly affect overall ovarian activation or consumption of each experimental diet. However, there was an interaction between AA mix and dietary AA:C ratio affecting survival during the feeding experiment, and signs that there may have been weak, interactive effects of AA mix and AA:C ratio on food consumption. These results suggest that the effect of total AA intake on bumble bee phenotype may depend on the blend of individual AAs in experimental diets. We suggest that research exploring how AA blend affects bumble bee performance and dietary intake is warranted, and highlight that comparing research on bee nutrition is complicated by even subtle variation in experimental diet composition., (© 2021. The Author(s).)

In the past centuries, viruses have benefited from globalization to spread across the globe, infecting new host species and populations. A growing number of viruses have been documented in the western honey bee, Apis mellifera . Several of these contribute significantly to honey bee colony losses. This review synthetizes the knowledge of the diversity and distribution of honey-bee-infecting viruses, including recent data from high-throughput sequencing (HTS). After presenting the diversity of viruses and their corresponding symptoms, we surveyed the scientific literature for the prevalence of these pathogens across the globe. The geographical distribution shows that the most prevalent viruses (deformed wing virus, sacbrood virus, black queen cell virus and acute paralysis complex) are also the most widely distributed. We discuss the ecological drivers that influence the distribution of these pathogens in worldwide honey bee populations. Besides the natural transmission routes and the resulting temporal dynamics, global trade contributes to their dissemination. As recent evidence shows that these viruses are often multihost pathogens, their spread is a risk for both the beekeeping industry and the pollination services provided by managed and wild pollinators.

Patterns and drivers of the spatio-temporal distribution of herbivores are key elements of their ecological and evolutionary impacts on plant populations. Herbivore spatial distributions may be influenced by increased (RCH: resource concentration hypothesis) or decreased (RDH: resource dilution hypothesis) resource densities, but the effect of temporal variations in resource densities on such distributions remains poorly documented. We used a survey of a masting tree species and its seed predators in Southeastern France to address the effect of a host's pulsed resource on the spatio-temporal distributions of highly specialized insect herbivores feeding on seeds. Variations in both resource and seed predator densities were assessed by estimating seed production and seed infestation rates in focus trees during 10 consecutive years. We found increasing seed infestation rates with decreasing host tree densities in years of low seed production, indicating a RDH pattern of seed predators. However, such pattern was not persistent in years of high seed production during which seed infestation rates did not depend on host tree densities. We showed that temporal variations in resource density can lead to transience of seed predator spatial distribution. This study highlights how predictions of plant-herbivore interactions in natural ecosystems may rely on temporal components underlying RCH and RDH hypotheses.

Plant-insect interactions are key model systems to assess how some species affect the distribution, the abundance, and the evolution of others. Tree reproductive structures represent a critical resource for many insect species, which can be likely drivers of demography, spatial distribution, and trait diversification of plants. In this review, we present the ecological implications of predispersal herbivory on tree reproductive structures by insects (PIHR) in forest ecosystems. Both insect's and tree's perspectives are addressed with an emphasis on how spatiotemporal variation and unpredictability in seed availability can shape such particular plant-animal interactions. Reproductive structure insects show strong trophic specialization and guild diversification. Insects evolved host selection and spatiotemporal dispersal strategies in response to variable and unpredictable abundance of reproductive structures in both space and time. If PIHR patterns have been well documented in numerous systems, evidences of the subsequent demographic and evolutionary impacts on tree populations are still constrained by time-scale challenges of experimenting on such long-lived organisms, and modeling approaches of tree dynamics rarely consider PIHR when including biotic interactions in their processes. We suggest that spatially explicit and mechanistic approaches of the interactions between individual tree fecundity and insect dynamics will clarify predictions of the demogenetic implications of PIHR in tree populations. In a global change context, further experimental and theoretical contributions to the likelihood of life-cycle disruptions between plants and their specialized herbivores, and to how these changes may generate novel dynamic patterns in each partner of the interaction are increasingly critical., (© 2017 Institute of Zoology, Chinese Academy of Sciences.)

Deformed wing virus (DWV) is an emerging infectious disease of the honey bee ( Apis mellifera ) that is considered a major cause of elevated losses of honey bee colonies. DWV comprises two widespread genotypes: the originally described genotype A, and genotype B. In adult honey bees, DWV-B has been shown to be more virulent than DWV-A. However, their comparative effects on earlier host developmental stages are unknown. Here, we experimentally inoculated honey bee pupae and tested for the relative impact of DWV-A versus DWV-B on mortality and wing deformities in eclosing adults. DWV-A and DWV-B caused similar, and only slightly elevated, pupal mortality (mean 18% greater mortality than control). Both genotypes caused similarly high wing deformities in eclosing adults (mean 60% greater wing deformities than control). Viral titer was high in all of the experimentally inoculated eclosing adults, and was independent of wing deformities, suggesting that the phenotype 'deformed wings' is not directly related to viral titer or viral genotype. These viral traits favor the emergence of both genotypes of DWV by not limiting the reproduction of its vector, the ectoparasitic Varroa destructor mite, in infected pupae, and thereby facilitating the spread of DWV in honey bees infested by the mite.

Bees are considered to be threatened globally, with severe overwinter losses of the most important commercial pollinator, the Western honeybee, a major concern in the Northern Hemisphere. Emerging infectious diseases have risen to prominence due to their temporal correlation with colony losses. Among these is Deformed wing virus (DWV), which has been frequently linked to colony mortality. We now provide evidence of a strong statistical association between overwintering colony decline in the field and the presence of DWV genotype-B (DWV-B), a genetic variant of DWV that has recently been shown to be more virulent than the original DWV genotype-A. We link the prevalence of DWV-B directly to a quantitative measure of overwinter decline (workforce mortality) of honeybee colonies in the field. We demonstrate that increased prevalence of virus infection in individual bees is associated with higher overwinter mortality. We also observed a substantial reduction of infected colonies in the spring, suggesting that virus-infected individuals had died during the winter. Our findings demonstrate that DWV-B, plus possible A/B recombinants exhibiting DWV-B at PCR primer binding sites, may be a major cause of elevated overwinter honeybee loss. Its potential emergence in naïve populations of bees may have far-reaching ecological and economic impacts.

Background: Organisms typically face infection by diverse pathogens, and hosts are thought to have developed specific responses to each type of pathogen they encounter. The advent of transcriptomics now makes it possible to test this hypothesis and compare host gene expression responses to multiple pathogens at a genome-wide scale. Here, we performed a meta-analysis of multiple published and new transcriptomes using a newly developed bioinformatics approach that filters genes based on their expression profile across datasets. Thereby, we identified common and unique molecular responses of a model host species, the honey bee (Apis mellifera), to its major pathogens and parasites: the Microsporidia Nosema apis and Nosema ceranae, RNA viruses, and the ectoparasitic mite Varroa destructor, which transmits viruses., Results: We identified a common suite of genes and conserved molecular pathways that respond to all investigated pathogens, a result that suggests a commonality in response mechanisms to diverse pathogens. We found that genes differentially expressed after infection exhibit a higher evolutionary rate than non-differentially expressed genes. Using our new bioinformatics approach, we unveiled additional pathogen-specific responses of honey bees; we found that apoptosis appeared to be an important response following microsporidian infection, while genes from the immune signalling pathways, Toll and Imd, were differentially expressed after Varroa/virus infection. Finally, we applied our bioinformatics approach and generated a gene co-expression network to identify highly connected (hub) genes that may represent important mediators and regulators of anti-pathogen responses., Conclusions: Our meta-analysis generated a comprehensive overview of the host metabolic and other biological processes that mediate interactions between insects and their pathogens. We identified key host genes and pathways that respond to phylogenetically diverse pathogens, representing an important source for future functional studies as well as offering new routes to identify or generate pathogen resilient honey bee stocks. The statistical and bioinformatics approaches that were developed for this study are broadly applicable to synthesize information across transcriptomic datasets. These approaches will likely have utility in addressing a variety of biological questions.

Regulation of gene expression in the brain plays an important role in behavioral plasticity and decision making in response to external stimuli. However, both can be severely affected by environmental factors, such as parasites and pathogens. In honey bees, the emergence and re-emergence of pathogens and potential for pathogen co-infection and interaction have been suggested as major components that significantly impaired social behavior and survival. To understand how the honey bee is affected and responds to interacting pathogens, we co-infected workers with two prevalent pathogens of different nature, the positive single strand RNA virus Black queen cell virus (BQCV), and the Microsporidia Nosema ceranae , and explored gene expression changes in brains upon single infections and co-infections. Our data provide an important resource for research on honey bee diseases, and more generally on insect host-pathogen and pathogen-pathogen interactions. Raw and processed data are publicly available in the NCBI/GEO database: (http://www.ncbi.nlm.nih.gov/geo/) under accession number GSE81664.

Emerging infectious diseases (EIDs) have contributed significantly to the current biodiversity crisis, leading to widespread epidemics and population loss. Owing to genetic variation in pathogen virulence, a complete understanding of species decline requires the accurate identification and characterization of EIDs. We explore this issue in the Western honeybee, where increasing mortality of populations in the Northern Hemisphere has caused major concern. Specifically, we investigate the importance of genetic identity of the main suspect in mortality, deformed wing virus (DWV), in driving honeybee loss. Using laboratory experiments and a systematic field survey, we demonstrate that an emerging DWV genotype (DWV-B) is more virulent than the established DWV genotype (DWV-A) and is widespread in the landscape. Furthermore, we show in a simple model that colonies infected with DWV-B collapse sooner than colonies infected with DWV-A. We also identify potential for rapid DWV evolution by revealing extensive genome-wide recombination in vivo The emergence of DWV-B in naive honeybee populations, including via recombination with DWV-A, could be of significant ecological and economic importance. Our findings emphasize that knowledge of pathogen genetic identity and diversity is critical to understanding drivers of species decline., (© 2016 The Authors.)

As pollinators, bees are cornerstones for terrestrial ecosystem stability and key components in agricultural productivity. All animals, including bees, are associated with a diverse community of microbes, commonly referred to as the microbiome. The bee microbiome is likely to be a crucial factor affecting host health. However, with the exception of a few pathogens, the impacts of most members of the bee microbiome on host health are poorly understood. Further, the evolutionary and ecological forces that shape and change the microbiome are unclear. Here, we discuss recent progress in our understanding of the bee microbiome, and we present challenges associated with its investigation. We conclude that global coordination of research efforts is needed to fully understand the complex and highly dynamic nature of the interplay between the bee microbiome, its host, and the environment. High-throughput sequencing technologies are ideal for exploring complex biological systems, including host-microbe interactions. To maximize their value and to improve assessment of the factors affecting bee health, sequence data should be archived, curated, and analyzed in ways that promote the synthesis of different studies. To this end, the BeeBiome consortium aims to develop an online database which would provide reference sequences, archive metadata, and host analytical resources. The goal would be to support applied and fundamental research on bees and their associated microbes and to provide a collaborative framework for sharing primary data from different research programs, thus furthering our understanding of the bee microbiome and its impact on pollinator health., (Copyright © 2016 Engel et al.)

Microbial pathogens are thought to have a profound impact on insect populations. Honey bees are suffering from elevated colony losses in the northern hemisphere possibly because of a variety of emergent microbial pathogens, with which pesticides may interact to exacerbate their impacts. To reveal such potential interactions, we administered at sublethal and field realistic doses one neonicotinoid pesticide (thiacloprid) and two common microbial pathogens, the invasive microsporidian Nosema ceranae and black queen cell virus (BQCV), individually to larval and adult honey bees in the laboratory. Through fully crossed experiments in which treatments were administered singly or in combination, we found an additive interaction between BQCV and thiacloprid on host larval survival likely because the pesticide significantly elevated viral loads. In adult bees, two synergistic interactions increased individual mortality: between N. ceranae and BQCV, and between N. ceranae and thiacloprid. The combination of two pathogens had a more profound effect on elevating adult mortality than N. ceranae plus thiacloprid. Common microbial pathogens appear to be major threats to honey bees, while sublethal doses of pesticide may enhance their deleterious effects on honey bee larvae and adults. It remains an open question as to whether these interactions can affect colony survival., (© 2014 Society for Applied Microbiology and John Wiley & Sons Ltd.)

There is increasing appreciation that hosts in natural populations are subject to infection by multiple parasite species. Yet the epidemiological and ecological processes determining the outcome of mixed infections are poorly understood. Here, we use two intracellular gut parasites (Microsporidia), one exotic and one co-evolved in the western honeybee (Apis mellifera), in an experiment in which either one or both parasites were administered either simultaneously or sequentially. We provide clear evidence of within-host competition; order of infection was an important determinant of the competitive outcome between parasites, with the first parasite significantly inhibiting the growth of the second, regardless of species. However, the strength of this 'priority effect' was highly asymmetric, with the exotic Nosema ceranae exhibiting stronger inhibition of Nosema apis than vice versa. Our results reveal an unusual asymmetry in parasite competition that is dependent on order of infection. When incorporated into a mathematical model of disease prevalence, we find asymmetric competition to be an important predictor of the patterns of parasite prevalence found in nature. Our findings demonstrate the wider significance of complex multi-host-multi-parasite interactions as drivers of host-pathogen community structure., (© 2014 The Author(s) Published by the Royal Society. All rights reserved.)

Two pathogens co-infecting a common host can either interact positively (facilitation), negatively (competition) or act independently. A correlative study has suggested that two pathogens of the honey bee, Nosema ceranae and Deformed wing virus (DWV), interact negatively within a host (Costa et al., 2011). To test this hypothesis, we sequentially co-infected honey bees with these pathogens in a reciprocally crossed experimental design. Prior establishment in the host ventriculus by N. ceranae inhibited DWV while prior infection by DWV did not impact N. ceranae, highlighting an asymmetry in the competitive interaction between these emerging pathogens., (Copyright © 2014 Elsevier Inc. All rights reserved.)

A faithful expression of the mitochondrial DNA is crucial for cell survival. Animal mitochondrial DNA (mtDNA) presents a highly compact gene organization. The typical 16.5 kbp animal mtDNA encodes 13 proteins, 2 rRNAs and 22 tRNAs. In the backyard pillbug Armadillidium vulgare, the rather small 13.9 kbp mtDNA encodes the same set of proteins and rRNAs as compared to animal kingdom mtDNA, but seems to harbor an incomplete set of tRNA genes. Here, we first confirm the expression of 13 tRNA genes in this mtDNA. Then we show the extensive repair of a truncated tRNA, the expression of tRNA involved in large gene overlaps and of tRNA genes partially or fully integrated within protein-coding genes in either direct or opposite orientation. Under selective pressure, overlaps between genes have been likely favored for strong genome size reduction. Our study underlines the existence of unknown biochemical mechanisms for the complete gene expression of A. vulgare mtDNA, and of co-evolutionary processes to keep overlapping genes functional in a compacted mitochondrial genome.

The aim of this study was to improve cage systems for maintaining adult honey bee (Apis mellifera L.) workers under in vitro laboratory conditions. To achieve this goal, we experimentally evaluated the impact of different cages, developed by scientists of the international research network COLOSS (Prevention of honey bee COlony LOSSes), on the physiology and survival of honey bees. We identified three cages that promoted good survival of honey bees. The bees from cages that exhibited greater survival had relatively lower titers of deformed wing virus, suggesting that deformed wing virus is a significant marker reflecting stress level and health status of the host. We also determined that a leak- and drip-proof feeder was an integral part of a cage system and a feeder modified from a 20-ml plastic syringe displayed the best result in providing steady food supply to bees. Finally, we also demonstrated that the addition of protein to the bees' diet could significantly increase the level ofvitellogenin gene expression and improve bees' survival. This international collaborative study represents a critical step toward improvement of cage designs and feeding regimes for honey bee laboratory experiments.

Mitochondrial DNA (mtDNA) is usually depicted as a circular molecule, however, there is increasing evidence that linearization of mtDNA evolved independently many times in organisms such as fungi, unicellular eukaryotes, and animals. Recent observations in various models with linear mtDNA revealed the presence of conserved inverted repeats (IR) at both ends that, when they become single-stranded, may be able to fold on themselves to create telomeric-hairpins involved in genome architecture conversions. The atypical mtDNA of terrestrial isopods (Crustacea: Oniscidea) composed of linear monomers and circular dimers is an interesting model to study genome architecture conversions. Here, we present the mtDNA control region sequences of two species of the genus Armadillidium: A. vulgare and A. pelagicum. All features of arthropods mtDNA control regions are present (origin of replication, poly-T stretch, GA and TA-rich blocks and one variable domain), plus a conserved IR. This IR can potentially fold into a hairpin structure and is present in two different orientations among the A. vulgare populations: either in one sense or in its reverse complement. This polymorphism, also observed in a single individual (heteroplasmy), might be a signature of genome architecture conversions from linear to circular monomeric mtDNA via successive opening and closing of the molecules.

Metazoan mitochondrial DNA (mtDNA) is generally composed of circular monomeric molecules. However, a few exceptions do exist and among them two terrestrial isopods Armadillidium vulgare and Porcellionides pruinosus have an atypical mtDNA composed of linear monomers associated with circular "head-to-head" dimers: a very unusual structure for animal mtDNA genome. To assess the distribution of this atypical mtDNA among isopods, we performed RFLP and Southern blot analyses on mtDNA of 16 terrestrial (Oniscidea family) and two aquatic isopod species: the marine Sphaeroma serratum (suborder Flabellifera, sister group of Oniscidea) and the freshwater Asellus aquaticus (Asellota, early derived taxon of isopod). The atypical mtDNA structure was observed in 15 terrestrial isopod species and A. aquaticus, suggesting a wide distribution of atypical mtDNA among isopods. However, a typical metazoan mtDNA structure was detected in the marine isopod S. serratum and the Oniscidea Ligia oceanica . Our results suggest two possible scenarios: an early origin of the atypical mtDNA in isopods followed by reversion to the typical ancestral mtDNA structure for several species, or a convergent appearance of the atypical mtDNA structure in two isopod suborders. We compare this distribution of the atypical mtDNA structure with the presence of a heteroplasmy also observed in the mtDNA of several terrestrial isopod species. We discuss if this transmitted heteroplasmy is vectored by the atypical mtDNA and its impact on the maintenance of the atypical mtDNA in isopods.

Wolbachia bacteria are obligate intracellular alpha-Proteobacteria of arthropods and nematodes. Although widespread among isopod crustaceans, they have seldom been found in non-isopod crustacean species. Here, we report Wolbachia infection in fourteen new crustacean species. Our results extend the range of Wolbachia infections in terrestrial isopods and amphipods (class Malacostraca). We report the occurrence of two different Wolbachia strains in two host species (a terrestrial isopod and an amphipod). Moreover, the discovery of Wolbachia in the goose barnacle Lepas anatifera (subclass Thecostraca) establishes Wolbachia infection in class Maxillopoda. The new bacterial strains are closely related to B-supergroup Wolbachia strains previously reported from crustacean hosts. Our results suggest that Wolbachia infection may be much more widespread in crustaceans than previously thought. The presence of related Wolbachia strains in highly divergent crustacean hosts suggests that Wolbachia endosymbionts can naturally adapt to a wide range of crustacean hosts. Given the ability of isopod Wolbachia strains to induce feminization of genetic males or cytoplasmic incompatibility, we speculate that manipulation of crustacean-borne Wolbachia bacteria might represent potential tools for controlling crustacean species of commercial interest and crustacean or insect disease vectors.

Due to essentially maternal inheritance and a bottleneck effect during early oogenesis, newly arising mitochondrial DNA (mtDNA) mutations segregate rapidly in metazoan female germlines. Consequently, heteroplasmy (i.e. the mixture of mtDNA genotypes within an organism) is generally resolved to homoplasmy within a few generations. Here, we report an exceptional transpecific heteroplasmy (predicting an alanine/valine alloacceptor tRNA change) that has been stably inherited in oniscid crustaceans for at least thirty million years. Our results suggest that this heteroplasmy is stably transmitted across generations because it occurs within mitochondria and therefore escapes the mtDNA bottleneck that usually erases heteroplasmy. Consistently, at least two oniscid species possess an atypical trimeric mitochondrial genome, which provides an adequate substrate for the emergence of a constitutive intra-mitochondrial heteroplasmy. Persistence of a mitochondrial polymorphism on such a deep evolutionary timescale suggests that balancing selection may be shaping mitochondrial sequence evolution in oniscid crustaceans.

Aphids exclusively feed on plant phloem sap that contains much sugar and some nonessential amino acids but is poor in lipids and proteins. Conventionally, it has been believed that aphids substantially have no intestinal digestion of proteins. However, we here report an unexpected finding that cysteine protease genes of the family cathepsin B are massively amplified in the lineage of aphids and that many of the protease genes exhibit gut-specific overexpression. By making use of expressed sequence tag data, sequenced cDNAs, and genomic trace sequences of the pea aphid Acyrthosiphon pisum, we identified a total of 28 cathepsin B-like gene copies in the genome of A. pisum. Phylogenetic analyses of all the cathepsin B genes in aphids revealed that genic expansion has continuously proceeded with basal, intermediary, and recent duplications. Estimation of molecular evolutionary rates indicated that major alterations of the rates often occurred after duplications. For example, a gene copy ("348") was shown to be slow evolving and close to genes of other insects like Drosophila melanogaster, whereas the other gene copies appeared to have evolved faster with higher ratios of nonsynonymous to synonymous substitutions. We identified a number of gene copies (16 in A. pisum) that contained a replacement at the site required for catalytic activity of the protease. Among these, 2 copies were pseudogenes, whereas the remaining copies were structurally intact and possibly acquired new functions. For example, a cluster of such gene copies ("1674") has been subjected to positive selection. Quantitative reverse transcriptase-polymerase chain reaction analyses revealed that the more conserved gene copy ("348") showed a constitutive expression, whereas 5 other forms ("84," "16," "16D," "1874," and "2744") were preferentially expressed in the gut of A. pisum. Putative biological roles of the diversified cathepsin B-like gene copies in aphids are discussed in relation to their nutritional physiology specialized for plant sap feeding lifestyle.

The crustacean isopod Armadillidium vulgare is characterized by an unusual approximately 42-kb-long mitochondrial genome consisting of two molecules co-occurring in mitochondria: a circular approximately 28-kb dimer formed by two approximately 14-kb monomers fused in opposite polarities and a linear approximately 14-kb monomer. Here we determined the nucleotide sequence of the fundamental monomeric unit of A. vulgare mitochondrial genome, to gain new insight into its structure and evolution. Our results suggest that the junction zone between monomers of the dimer structure is located in or near the control region. Direct sequencing indicated that the nucleotide sequences of the different monomer units are virtually identical. This suggests that gene conversion and/or replication processes play an important role in shaping nucleotide sequence variation in this mitochondrial genome. The only heteroplasmic site we identified predicts an alloacceptor tRNA change from tRNA(Ala) to tRNA(Val). Therefore, in A. vulgare, tRNA(Ala) and tRNA(Val) are found at the same locus in different monomers, ensuring that both tRNAs are present in mitochondria. The presence of this heteroplasmic site in all sequenced individuals suggests that the polymorphism is selectively maintained, probably because of the necessity of both tRNAs for maintaining proper mitochondrial functions. Thus, our results provide empirical evidence for the tRNA gene recruitment model of tRNA evolution. Moreover, interspecific comparisons showed that the A. vulgare mitochondrial gene order is highly derived compared to the putative ancestral arthropod type. By contrast, an overall high conservation of mitochondrial gene order is observed within crustacean isopods.