Co-option of immune effectors by the hormonal signalling system triggering metamorphosis in Drosophila melanogaster.

Insect metamorphosis is triggered by the production, secretion and degradation of 20-hydroxyecdysone (ecdysone). In addition to its role in developmental regulation, increasing evidence suggests that ecdysone is involved in innate immunity processes, such as phagocytosis and the induction of antimicrobial peptide (AMP) production. AMP regulation includes systemic responses as well as local responses at surface epithelia that contact with the external environment. At pupariation, Drosophila melanogaster increases dramatically the expression of three AMP genes, drosomycin (drs), drosomycin-like 2 (drsl2) and drosomycin-like 5 (drsl5). We show that the systemic action of drs at pupariation is dependent on ecdysone signalling in the fat body and operates via the ecdysone downstream target, Broad. In parallel, ecdysone also regulates local responses, specifically through the activation of drsl2 expression in the gut. Finally, we confirm the relevance of this ecdysone dependent AMP expression for the control of bacterial load by showing that flies lacking drs expression in the fat body have higher bacterial persistence over metamorphosis. In contrast, local responses may be redundant with the systemic effect of drs since reduction of ecdysone signalling or of drsl2 expression has no measurable negative effect on bacterial load control in the pupa. Together, our data emphasize the importance of the association between ecdysone signalling and immunity using in vivo studies and establish a new role for ecdysone at pupariation, which impacts developmental success by regulating the immune system in a stage-dependent manner. We speculate that this co-option of immune effectors by the hormonal system may constitute an anticipatory mechanism to control bacterial numbers in the pupa, at the core of metamorphosis evolution.

Genomic Analysis of European Drosophila melanogaster Populations Reveals Longitudinal Structure, Continent-Wide Selection, and Previously Unknown DNA Viruses.

Genetic variation is the fuel of evolution, with standing genetic variation especially important for short-term evolution and local adaptation. To date, studies of spatiotemporal patterns of genetic variation in natural populations have been challenging, as comprehensive sampling is logistically difficult, and sequencing of entire populations costly. Here, we address these issues using a collaborative approach, sequencing 48 pooled population samples from 32 locations, and perform the first continent-wide genomic analysis of genetic variation in European Drosophila melanogaster. Our analyses uncover longitudinal population structure, provide evidence for continent-wide selective sweeps, identify candidate genes for local climate adaptation, and document clines in chromosomal inversion and transposable element frequencies. We also characterize variation among populations in the composition of the fly microbiome, and identify five new DNA viruses in our samples.

Spider mites escape bacterial infection by avoiding contaminated food.

To fight infection, arthropods rely on the deployment of an innate immune response but also upon physical/chemical barriers and avoidance behaviours. However, most studies focus on immunity, with other defensive mechanisms being relatively overlooked. We have previously shown that the spider mite Tetranychus urticae does not mount an induced immune response towards systemic bacterial infections, entailing very high mortality rates. Therefore, we hypothesized that other defence mechanisms may be operating to minimize infection risk. Here, we test (a) if spider mites are also highly susceptible to other infection routes-spraying and feeding-and (b) if they display avoidance behaviours towards infected food. Individuals sprayed with or fed on Escherichia coli or Pseudomonas putida survived less than the control, pointing to a deficient capacity of the gut epithelium, and possibly of the cuticle, to contain bacteria. Additionally, we found that spider mites prefer uninfected food to food contaminated with bacteria, a choice that probably does not rely on olfactory cues. Our results suggest that spider mites may rely mostly on avoidance behaviours to minimize bacterial infection and highlight the multi-layered nature of immune strategies present in arthropods.

Drosophila Adaptation to Viral Infection through Defensive Symbiont Evolution.

Microbial symbionts can modulate host interactions with biotic and abiotic factors. Such interactions may affect the evolutionary trajectories of both host and symbiont. Wolbachia protects Drosophila melanogaster against several viral infections and the strength of the protection varies between variants of this endosymbiont. Since Wolbachia is maternally transmitted, its fitness depends on the fitness of its host. Therefore, Wolbachia populations may be under selection when Drosophila is subjected to viral infection. Here we show that in D. melanogaster populations selected for increased survival upon infection with Drosophila C virus there is a strong selection coefficient for specific Wolbachia variants, leading to their fixation. Flies carrying these selected Wolbachia variants have higher survival and fertility upon viral infection when compared to flies with the other variants. These findings demonstrate how the interaction of a host with pathogens shapes the genetic composition of symbiont populations. Furthermore, host adaptation can result from the evolution of its symbionts, with host and symbiont functioning as a single evolutionary unit.

Diverse Cis-Regulatory Mechanisms Contribute to Expression Evolution of Tandem Gene Duplicates.

Pairs of duplicated genes generally display a combination of conserved expression patterns inherited from their unduplicated ancestor and newly acquired domains. However, how the cis-regulatory architecture of duplicated loci evolves to produce these expression patterns is poorly understood. We have directly examined the gene-regulatory evolution of two tandem duplicates, the Drosophila Ly6 genes CG9336 and CG9338, which arose at the base of the drosophilids between 40 and 60 Ma. Comparing the expression patterns of the two paralogs in four Drosophila species with that of the unduplicated ortholog in the tephritid Ceratitis capitata, we show that they diverged from each other as well as from the unduplicated ortholog. Moreover, the expression divergence appears to have occurred close to the duplication event and also more recently in a lineage-specific manner. The comparison of the tissue-specific cis-regulatory modules (CRMs) controlling the paralog expression in the four Drosophila species indicates that diverse cis-regulatory mechanisms, including the novel tissue-specific enhancers, differential inactivation, and enhancer sharing, contributed to the expression evolution. Our analysis also reveals a surprisingly variable cis-regulatory architecture, in which the CRMs driving conserved expression domains change in number, location, and specificity. Altogether, this study provides a detailed historical account that uncovers a highly dynamic picture of how the paralog expression patterns and their underlying cis-regulatory landscape evolve. We argue that our findings will encourage studying cis-regulatory evolution at the whole-locus level to understand how interactions between enhancers and other regulatory levels shape the evolution of gene expression.

Tetranychus urticae mites do not mount an induced immune response against bacteria.

The genome of the spider mite Tetranychus urticae, a herbivore, is missing important elements of the canonical Drosophila immune pathways necessary to fight bacterial infections. However, it is not known whether spider mites can mount an immune response and survive bacterial infection. In other chelicerates, bacterial infection elicits a response mediated by immune effectors leading to the survival of infected organisms. In T. urticae, infection by either Escherichia coli or Bacillus megaterium did not elicit a response as assessed through genome-wide transcriptomic analysis. In line with this, spider mites died within days even upon injection with low doses of bacteria that are non-pathogenic to Drosophila Moreover, bacterial populations grew exponentially inside the infected spider mites. By contrast, Sancassania berlesei, a litter-dwelling mite, controlled bacterial proliferation and resisted infections with both Gram-negative and Gram-positive bacteria lethal to T. urticae This differential mortality between mite species was absent when mites were infected with heat-killed bacteria. Also, we found that spider mites harbour in their gut 1000-fold less bacteria than S. berlesei We show that T. urticae has lost the capacity to mount an induced immune response against bacteria, in contrast to other mites and chelicerates but similarly to the phloem feeding aphid Acyrthosiphon pisum Hence, our results reinforce the putative evolutionary link between ecological conditions regarding exposure to bacteria and the architecture of the immune response.

The genome of Tetranychus urticae reveals herbivorous pest adaptations.

The spider mite Tetranychus urticae is a cosmopolitan agricultural pest with an extensive host plant range and an extreme record of pesticide resistance. Here we present the completely sequenced and annotated spider mite genome, representing the first complete chelicerate genome. At 90 megabases T. urticae has the smallest sequenced arthropod genome. Compared with other arthropods, the spider mite genome shows unique changes in the hormonal environment and organization of the Hox complex, and also reveals evolutionary innovation of silk production. We find strong signatures of polyphagy and detoxification in gene families associated with feeding on different hosts and in new gene families acquired by lateral gene transfer. Deep transcriptome analysis of mites feeding on different plants shows how this pest responds to a changing host environment. The T. urticae genome thus offers new insights into arthropod evolution and plant-herbivore interactions, and provides unique opportunities for developing novel plant protection strategies.

Host adaptation to viruses relies on few genes with different cross-resistance properties.

Host adaptation to one parasite may affect its response to others. However, the genetics of these direct and correlated responses remains poorly studied. The overlap between these responses is instrumental for the understanding of host evolution in multiparasite environments. We determined the genetic and phenotypic changes underlying adaptation of Drosophila melanogaster to Drosophila C virus (DCV). Within 20 generations, flies selected with DCV showed increased survival after DCV infection, but also after cricket paralysis virus (CrPV) and flock house virus (FHV) infection. Whole-genome sequencing identified two regions of significant differentiation among treatments, from which candidate genes were functionally tested with RNAi. Three genes were validated--pastrel, a known DCV-response gene, and two other loci, Ubc-E2H and CG8492. Knockdown of Ubc-E2H and pastrel also led to increased sensitivity to CrPV, whereas knockdown of CG8492 increased susceptibility to FHV infection. Therefore, Drosophila adaptation to DCV relies on few major genes, each with different cross-resistance properties, conferring host resistance to several parasites.

Multispecies Analysis of Expression Pattern Diversification in the Recently Expanded Insect Ly6 Gene Family.

Gene families often consist of members with diverse expression domains reflecting their functions in a wide variety of tissues. However, how the expression of individual members, and thus their tissue-specific functions, diversified during the course of gene family expansion is not well understood. In this study, we approached this question through the analysis of the duplication history and transcriptional evolution of a rapidly expanding subfamily of insect Ly6 genes. We analyzed different insect genomes and identified seven Ly6 genes that have originated from a single ancestor through sequential duplication within the higher Diptera. We then determined how the original embryonic expression pattern of the founding gene diversified by characterizing its tissue-specific expression in the beetle Tribolium castaneum, the butterfly Bicyclus anynana, and the mosquito Anopheles stephensi and those of its duplicates in three higher dipteran species, representing various stages of the duplication history (Megaselia abdita, Ceratitis capitata, and Drosophila melanogaster). Our results revealed that frequent neofunctionalization episodes contributed to the increased expression breadth of this subfamily and that these events occurred after duplication and speciation events at comparable frequencies. In addition, at each duplication node, we consistently found asymmetric expression divergence. One paralog inherited most of the tissue-specificities of the founder gene, whereas the other paralog evolved drastically reduced expression domains. Our approach attests to the power of combining a well-established duplication history with a comprehensive coverage of representative species in acquiring unequivocal information about the dynamics of gene expression evolution in gene families.

Genetics of host-parasite interactions: towards a comprehensive dissection of Drosophila resistance to viral infection.

One of the major challenges in evolutionary biology is to unravel the genetic basis of adaptation. This issue has been gaining momentum in recent years with the accelerated development of novel genetic and genomic techniques and resources. In this issue of Molecular Ecology, Cogni et al. (2016) address the genetic basis of resistance to two viruses in Drosophila melanogaster using a panel of recombinant inbred lines with unprecedented resolution allowing detection of rare alleles and/or alleles of small effect. The study confirms the role of previously identified genes of major effect and adds novel regions with minor effect to the genetic basis of Drosophila resistance to the Drosophila C virus or the sigma virus. Additional analyses reveal the absence of cross-resistance and of epistasis between the various genomic regions. This detailed information on the genetic architecture of host resistance constitutes an important step towards the understanding of both the physiology of antiviral immunity and the evolution of host-parasite interactions.

Host adaptation is contingent upon the infection route taken by pathogens.

Evolution of pathogen virulence is affected by the route of infection. Also, alternate infection routes trigger different physiological responses on hosts, impinging on host adaptation and on its interaction with pathogens. Yet, how route of infection may shape adaptation to pathogens has not received much attention at the experimental level. We addressed this question through the experimental evolution of an outbred Drosophila melanogaster population infected by two different routes (oral and systemic) with Pseudomonas entomophila. The two selection regimes led to markedly different evolutionary trajectories. Adaptation to infection through one route did not protect from infection through the alternate route, indicating distinct genetic bases. Finally, relatively to the control population, evolved flies were not more resistant to bacteria other than Pseudomonas and showed higher susceptibility to viral infections. These specificities and trade-offs may contribute to the maintenance of genetic variation for resistance in natural populations. Our data shows that the infection route affects host adaptation and thus, must be considered in studies of host-pathogen interaction.

Steroid hormone signaling is essential to regulate innate immune cells and fight bacterial infection in Drosophila.

Coupling immunity and development is essential to ensure survival despite changing internal conditions in the organism. Drosophila metamorphosis represents a striking example of drastic and systemic physiological changes that need to be integrated with the innate immune system. However, nothing is known about the mechanisms that coordinate development and immune cell activity in the transition from larva to adult. Here, we reveal that regulation of macrophage-like cells (hemocytes) by the steroid hormone ecdysone is essential for an effective innate immune response over metamorphosis. Although it is generally accepted that steroid hormones impact immunity in mammals, their action on monocytes (e.g. macrophages and neutrophils) is still not well understood. Here in a simpler model system, we used an approach that allows in vivo, cell autonomous analysis of hormonal regulation of innate immune cells, by combining genetic manipulation with flow cytometry, high-resolution time-lapse imaging and tissue-specific transcriptomic analysis. We show that in response to ecdysone, hemocytes rapidly upregulate actin dynamics, motility and phagocytosis of apoptotic corpses, and acquire the ability to chemotax to damaged epithelia. Most importantly, individuals lacking ecdysone-activated hemocytes are defective in bacterial phagocytosis and are fatally susceptible to infection by bacteria ingested at larval stages, despite the normal systemic and local production of antimicrobial peptides. This decrease in survival is comparable to the one observed in pupae lacking immune cells altogether, indicating that ecdysone-regulation is essential for hemocyte immune functions and survival after infection. Microarray analysis of hemocytes revealed a large set of genes regulated at metamorphosis by EcR signaling, among which many are known to function in cell motility, cell shape or phagocytosis. This study demonstrates an important role for steroid hormone regulation of immunity in vivo in Drosophila, and paves the way for genetic dissection of the mechanisms at work behind steroid regulation of innate immune cells.

A discrete model of Drosophila eggshell patterning reveals cell-autonomous and juxtacrine effects.

The Drosophila eggshell constitutes a remarkable system for the study of epithelial patterning, both experimentally and through computational modeling. Dorsal eggshell appendages arise from specific regions in the anterior follicular epithelium that covers the oocyte: two groups of cells expressing broad (roof cells) bordered by rhomboid expressing cells (floor cells). Despite the large number of genes known to participate in defining these domains and the important modeling efforts put into this developmental system, key patterning events still lack a proper mechanistic understanding and/or genetic basis, and the literature appears to conflict on some crucial points. We tackle these issues with an original, discrete framework that considers single-cell models that are integrated to construct epithelial models. We first build a phenomenological model that reproduces wild type follicular epithelial patterns, confirming EGF and BMP signaling input as sufficient to establish the major features of this patterning system within the anterior domain. Importantly, this simple model predicts an instructive juxtacrine signal linking the roof and floor domains. To explore this prediction, we define a mechanistic model that integrates the combined effects of cellular genetic networks, cell communication and network adjustment through developmental events. Moreover, we focus on the anterior competence region, and postulate that early BMP signaling participates with early EGF signaling in its specification. This model accurately simulates wild type pattern formation and is able to reproduce, with unprecedented level of precision and completeness, various published gain-of-function and loss-of-function experiments, including perturbations of the BMP pathway previously seen as conflicting results. The result is a coherent model built upon rules that may be generalized to other epithelia and developmental systems.

Reversion of developmental mode in insects: evolution from long germband to short germband in the polyembrionic wasp Macrocentrus cingulum Brischke.

Germband size in insects has played a central role in our understanding of insect patterning mechanisms and their evolution. The polarity of evolutionary change in insect patterning has been viewed so far as the unidirectional shift from the ancestral short germband patterning of basal hemimetabolous insects to the long germband patterning observed in most modern Holometabola. However, some orders of holometabolic insects display both short and long germband development, though the absence of a clear phylogenetic context does not permit definite conclusions on the polarity of change. Derived hymenoptera, that is, bees and wasps, represent a classical textbook example of long germband development. Yet, in some wasps putative short germband development has been described correlating with lifestyle changes, namely with evolution of endoparasitism and polyembryony. To address the potential reversion from long to short germband, we focused on the family Braconidae, which displays ancestral long germband development, and examined the derived polyembryonic braconid Macrocentrus cingulum. Using SEM analysis of M. cingulum embryogenesis coupled with analyses of embryonic patterning markers, we show that this wasp evolved short germband embryogenesis secondarily, in a way that is reminiscent of embryogenesis in the beetle Tribolium castaneum. This work shows that the evolution of germband size in insects is a reversible process that may correlate with other life-history traits and suggests broader implications on the mechanisms and evolvability of insect development.

Wolbachia in the Malpighian tubules: evolutionary dead-end or adaptation?

Facultative endosymbionts, such as Wolbachia, perpetuate by vertical transmission mostly through colonization of the germline during embryogenesis. The remaining Wolbachia inside the embryo are internalized in progenitor cells of the somatic tissue. This perpetuation strategy triggers a cyclic bacterial bottleneck across host generations. However, throughout the host's life history (Drosophila, for example), some somatic tissues such as the Malpighian tubules (MTs) show large numbers of Wolbachia. It is assumed that Wolbachia present in the progenitor cells of the MTs are confined to this somatic tissue, implicitly considering MTs as an evolutionary dead-end for these bacteria. Nevertheless, the fact that bacteria can survive and proliferate inside MTs suggests a different fate as they may access the host's reproductive system and persist in the host population through vertical transmission. Indeed, based on the particular physiological and developmental characteristics of MT, as well as of Wolbachia, we argue the bacteria present in the MTs may constitute a secondary pool of vertically transmitted bacteria. Moreover, somatic pools of Wolbachia capable of reaching the gonads and insure vertical transmission may also provide an interesting element to the elucidation of horizontal transmission mechanisms. Finally, we also speculate that somatic pools of Wolbachia may play an important role in host fitness, namely during viral infections. In brief, we argue that the somatic pools of Wolbachia, with special emphasis on the MT subset, deserve experimental attention as putative players in the physiology and evolution of both bacteria and hosts.

Alternative splicing and gene duplication in the evolution of the FoxP gene subfamily.

The FoxP gene subfamily of transcription factors is defined by its characteristic 110 amino acid long DNA-binding forkhead domain and plays essential roles in vertebrate biology. Its four members, FoxP1-P4, have been extensively characterized functionally. FoxP1, FoxP2, and FoxP4 are involved in lung, heart, gut, and central nervous system (CNS) development. FoxP3 is necessary and sufficient for the specification of regulatory T cells (Tregs) of the adaptive immune system. In Drosophila melanogaster, in silico predictions identify one unique FoxP subfamily gene member (CG16899) with no described function. We characterized this gene and established that it generates by alternative splicing two isoforms that differ in the forkhead DNA-binding domain. In D. melanogaster, both isoforms are expressed in the embryonic CNS, but in hemocytes, only isoform A is expressed, hinting to a putative modulation through alternative splicing of FoxP1 function in immunity and/or other hemocyte-dependent processes. Furthermore, we show that in vertebrates, this novel alternative splicing pattern is conserved for FoxP1. In mice, this new FoxP1 isoform is expressed in brain, liver, heart, testes, thymus, and macrophages (equivalent in function to hemocytes). This alternative splicing pattern has arisen at the base of the Bilateria, probably through exon tandem duplication. Moreover, our phylogenetic analysis suggests that in vertebrates, FoxP1 is more related to the FoxP gene ancestral form and the other three paralogues, originated through serial duplications, which only retained one of the alternative exons. Also, the newly described isoform differs from the other in amino acids critical for DNA-binding specificity. The integrity of its fold is maintained, but the molecule has lost the direct hydrogen bonding to DNA bases leading to a putatively lower specificity and possibly affinity toward DNA. With the present comparative study, through the integration of experimental and in silico studies of the FoxP gene subfamily across the animal kingdom, we establish a new model for the FoxP gene in invertebrates and for the vertebrate FoxP1 paralogue. Furthermore, we present a scenario for the structural evolution of this gene class and reveal new previously unsuspected levels of regulation for FoxP1 in the vertebrate system.

Ploidy mosaicism and allele-specific gene expression differences in the allopolyploid Squalius alburnoides.

BACKGROUND: Squalius alburnoides is an Iberian cyprinid fish resulting from an interspecific hybridisation between Squalius pyrenaicus females (P genome) and males of an unknown Anaecypris hispanica-like species (A genome). S. alburnoides is an allopolyploid hybridogenetic complex, which makes it a likely candidate for ploidy mosaicism occurrence, and is also an interesting model to address questions about gene expression regulation and genomic interactions. Indeed, it was previously suggested that in S. alburnoides triploids (PAA composition) silencing of one of the three alleles (mainly of the P allele) occurs. However, not a whole haplome is inactivated but a more or less random inactivation of alleles varying between individuals and even between organs of the same fish was seen.In this work we intended to correlate expression differences between individuals and/or between organs to the occurrence of mosaicism, evaluating if mosaics could explain previous observations and its impact on the assessment of gene expression patterns. RESULTS: To achieve our goal, we developed flow cytometry and cell sorting protocols for this system generating more homogenous cellular and transcriptional samples. With this set-up we detected 10% ploidy mosaicism within the S. alburnoides complex, and determined the allelic expression profiles of ubiquitously expressed genes (rpl8; gapdh and ß-actin) in cells from liver and kidney of mosaic and non-mosaic individuals coming from different rivers over a wide geographic range. CONCLUSIONS: Ploidy mosaicism occurs sporadically within the S. alburnoides complex, but in a frequency significantly higher than reported for other organisms. Moreover, we could exclude the influence of this phenomenon on the detection of variable allelic expression profiles of ubiquitously expressed genes (rpl8; gapdh and ß-actin) in cells from liver and kidney of triploid individuals. Finally, we determined that the expression patterns previously detected only in a narrow geographic range is not a local restricted phenomenon but is pervasive in rivers where S. pyrenaicus is sympatric with S. alburnoides.We discuss mechanisms that could lead to the formation of mosaic S. alburnoides and hypothesise about a relaxation of the mechanisms that impose a tight control over mitosis and ploidy control in mixoploids.

Plant-microbe symbioses: new insights into common roots.

Arbuscular mycorrhiza (AM), a type of plant-fungal endosymbiosis, and nodulation, a bacterial-plant endosymbiosis, are the most ubiquitous symbioses on earth. Recent findings have established part of a shared genetic basis underlying these interactions. Here, we approach root endosymbioses through the lens of the homology and modularity concepts aiming at further clarifying the proximate and ultimate causes for the establishment of these biological systems. We review the genetics that underlie interspecific signaling and its concomitant shift in genetic programs for either partner. Also, through the comparative analysis of genetic modules shared by AM and nodulation symbioses, we identify fundamental nodes in these networks, suggesting the elemental steps that may have permitted symbiotic adaptation. Here, we show that this approach, allied to recent technical advances in the study of genetic systems architecture, can provide clear testable hypotheses for the advancement of our understanding on the evolution and development of symbiotic systems.

Endocrine regulation of immunity in insects.

Organisms have constant contact with potentially harmful agents that can compromise their fitness. However, most of the times these agents fail to cause serious disease by virtue of the rapid and efficient immune responses elicited in the host that can range from behavioural adaptations to immune system triggering. The immune system of insects does not comprise the adaptive arm, making it less complex than that of vertebrates, but key aspects of the activation and regulation of innate immunity are conserved across different phyla. This is the case for the hormonal regulation of immunity as a part of the broad organismal responses to external conditions under different internal states. In insects, depending on the physiological circumstances, distinct hormones either enhance or suppress the immune response integrating individual (and often collective) responses physiologically and behaviourally. In this review, we provide an overview of our current knowledge on the endocrine regulation of immunity in insects, its mechanisms and implications on metabolic adaptation and behaviour. We highlight the importance of this multilayered regulation of immunity in survival and reproduction (fitness) and its dependence on the hormonal integration with other mechanisms and life-history traits.

Wolbachia and host intrinsic reproductive barriers contribute additively to postmating isolation in spider mites.

Wolbachia are maternally-inherited bacteria that induce cytoplasmic incompatibility in many arthropod species. However, the ubiquity of this isolation mechanism for host speciation processes remains elusive, as only few studies have examined Wolbachia-induced incompatibilities when host populations are not genetically compatible. Here, we used three populations of two genetically differentiated colour forms of the haplodiploid spider mite Tetranychus urticae to dissect the interaction between Wolbachia-induced and host-associated incompatibilities, and their relative contribution to postmating isolation. We found that these two sources of incompatibility act through different mechanisms in an additive fashion. Host-associated incompatibility contributes 1.5 times more than Wolbachia-induced incompatibility in reducing hybrid production, the former through an overproduction of haploid sons at the expense of diploid daughters (ca. 75% decrease) and the latter by increasing the embryonic mortality of daughters (by ca. 49%). Furthermore, regardless of cross direction, we observed near-complete F1 hybrid sterility and complete F2 hybrid breakdown between populations of the two forms, but Wolbachia did not contribute to this outcome. We thus show mechanistic independence and an additive nature of host-intrinsic and Wolbachia-induced sources of isolation. Wolbachia may contribute to reproductive isolation in this system, thereby potentially affecting host differentiation and distribution in the field.