Genetic footprints of adaptive divergence in the bracovirus of Cotesia sesamiae identified by targeted resequencing.

The African parasitoid wasp Cotesia sesamiae is a generalist species structured in locally adapted populations showing differences in host range. The recent discovery of Cotesia typhae, a specialist, sister species to C. sesamiae, provides a good framework to study the genetic determinants of parasitoid host range. To investigate the genomic bases of divergence between these populations and species, we used a targeted sequencing approach on 24 samples. We targeted the bracovirus genomic region encoding virulence genes involved in the interaction with the lepidopteran hosts of the wasps. High sequencing coverage was obtained for all samples, allowing the study of genetic variation between wasp populations and species. By combining population genetic estimations, such as nucleotide diversity (π), relative differentiation (FST ) and absolute divergence (dxy ), with branch-site dN/dS measures, we identified six of 98 bracovirus genes showing significant divergence and evidence of positive selection. These genes, belonging to different gene families, are potentially involved in host adaptation and in the specialization process. Fine-scale analyses of genetic variation also revealed mutations and large deletions in certain genes inducing pseudogenization and loss of function. The image emerging from these results is that adaptation mediated by bracovirus genes happens through selection of particularly adaptive alleles and loss of nonadaptive genes. These results highlight the central role of the bracovirus in the molecular interactions between the wasps and their hosts and in the evolutionary processes of specialization.

The discovery of Halictivirus resolves the Sinaivirus phylogeny.

By providing pollination services, bees are among the most important insects, both in ecological and economical terms. Combined next-generation and classical sequencing approaches were applied to discover and study new insect viruses potentially harmful to bees. A bioinformatics virus discovery pipeline was used on individual Illumina transcriptomes of 13 wild bees from three species from the genus Halictus and 30 ants from six species of the genera Messor and Aphaenogaster. This allowed the discovery and description of three sequences of a new virus termed Halictus scabiosae Adlikon virus (HsAV). Phylogenetic analyses of ORF1, RNA-dependent RNA-polymerase (RdRp) and capsid genes showed that HsAV is closely related to (+)ssRNA viruses of the unassigned Sinaivirus genus but distant enough to belong to a different new genus we called Halictivirus. In addition, our study of ant transcriptomes revealed the first four sinaivirus sequences from ants (Messor barbarus, M. capitatus and M. concolor). Maximum likelihood phylogenetic analyses were performed on a 594 nt fragment of the ORF1/RdRp region from 84 sinaivirus sequences, including 31 new Lake Sinai viruses (LSVs) from honey bees collected in five countries across the globe and the four ant viral sequences. The phylogeny revealed four main clades potentially representing different viral species infecting honey bees. Moreover, the ant viruses belonged to the LSV4 clade, suggesting a possible cross-species transmission between bees and ants. Lastly, wide honey bee screening showed that all four LSV clades have worldwide distributions with no obvious geographical segregation.

Molecular Evolution of Freshwater Snails with Contrasting Mating Systems.

Because mating systems affect population genetics and ecology, they are expected to impact the molecular evolution of species. Self-fertilizing species experience reduced effective population size, recombination rates, and heterozygosity, which in turn should decrease the efficacy of natural selection, both adaptive and purifying, and the strength of meiotic drive processes such as GC-biased gene conversion. The empirical evidence is only partly congruent with these predictions, depending on the analyzed species, some, but not all, of the expected effects have been observed. One possible reason is that self-fertilization is an evolutionary dead-end, so that most current selfers recently evolved self-fertilization, and their genome has not yet been strongly impacted by selfing. Here, we investigate the molecular evolution of two groups of freshwater snails in which mating systems have likely been stable for several millions of years. Analyzing coding sequence polymorphism, divergence, and expression levels, we report a strongly reduced genetic diversity, decreased efficacy of purifying selection, slower rate of adaptive evolution, and weakened codon usage bias/GC-biased gene conversion in the selfer Galba compared with the outcrosser Physa, in full agreement with theoretical expectations. Our results demonstrate that self-fertilization, when effective in the long run, is a major driver of population genomic and molecular evolutionary processes. Despite the genomic effects of selfing, Galba truncatula seems to escape the demographic consequences of the genetic load. We suggest that the particular ecology of the species may buffer the negative consequences of selfing, shedding new light on the dead-end hypothesis.

Reference-free population genomics from next-generation transcriptome data and the vertebrate-invertebrate gap.

In animals, the population genomic literature is dominated by two taxa, namely mammals and drosophilids, in which fully sequenced, well-annotated genomes have been available for years. Data from other metazoan phyla are scarce, probably because the vast majority of living species still lack a closely related reference genome. Here we achieve de novo, reference-free population genomic analysis from wild samples in five non-model animal species, based on next-generation sequencing transcriptome data. We introduce a pipe-line for cDNA assembly, read mapping, SNP/genotype calling, and data cleaning, with specific focus on the issue of hidden paralogy detection. In two species for which a reference genome is available, similar results were obtained whether the reference was used or not, demonstrating the robustness of our de novo inferences. The population genomic profile of a hare, a turtle, an oyster, a tunicate, and a termite were found to be intermediate between those of human and Drosophila, indicating that the discordant genomic diversity patterns that have been reported between these two species do not reflect a generalized vertebrate versus invertebrate gap. The genomic average diversity was generally higher in invertebrates than in vertebrates (with the notable exception of termite), in agreement with the notion that population size tends to be larger in the former than in the latter. The non-synonymous to synonymous ratio, however, did not differ significantly between vertebrates and invertebrates, even though it was negatively correlated with genetic diversity within each of the two groups. This study opens promising perspective regarding genome-wide population analyses of non-model organisms and the influence of population size on non-synonymous versus synonymous diversity.

Three infectious viral species lying in wait in the banana genome.

Plant pararetroviruses integrate serendipitously into their host genomes. The banana genome harbors integrated copies of banana streak virus (BSV) named endogenous BSV (eBSV) that are able to release infectious pararetrovirus. In this investigation, we characterized integrants of three BSV species-Goldfinger (eBSGFV), Imove (eBSImV), and Obino l'Ewai (eBSOLV)-in the seedy Musa balbisiana Pisang klutuk wulung (PKW) by studying their molecular structure, genomic organization, genomic landscape, and infectious capacity. All eBSVs exhibit extensive viral genome duplications and rearrangements. eBSV segregation analysis on an F1 population of PKW combined with fluorescent in situ hybridization analysis showed that eBSImV, eBSOLV, and eBSGFV are each present at a single locus. eBSOLV and eBSGFV contain two distinct alleles, whereas eBSImV has two structurally identical alleles. Genotyping of both eBSV and viral particles expressed in the progeny demonstrated that only one allele for each species is infectious. The infectious allele of eBSImV could not be identified since the two alleles are identical. Finally, we demonstrate that eBSGFV and eBSOLV are located on chromosome 1 and eBSImV is located on chromosome 2 of the reference Musa genome published recently. The structure and evolution of eBSVs suggest sequential integration into the plant genome, and haplotype divergence analysis confirms that the three loci display differential evolution. Based on our data, we propose a model for BSV integration and eBSV evolution in the Musa balbisiana genome. The mutual benefits of this unique host-pathogen association are also discussed.

Evidence for widespread positive and purifying selection across the European rabbit (Oryctolagus cuniculus) genome.

The nearly neutral theory of molecular evolution predicts that the efficacy of both positive and purifying selection is a function of the long-term effective population size (N(e)) of a species. Under this theory, the efficacy of natural selection should increase with N(e). Here, we tested this simple prediction by surveying ~1.5 to 1.8 Mb of protein coding sequence in the two subspecies of the European rabbit (Oryctolagus cuniculus algirus and O. c. cuniculus), a mammal species characterized by high levels of nucleotide diversity and N(e) estimates for each subspecies on the order of 1 × 10(6). When the segregation of slightly deleterious mutations and demographic effects were taken into account, we inferred that >60% of amino acid substitutions on the autosomes were driven to fixation by positive selection. Moreover, we inferred that a small fraction of new amino acid mutations (<4%) are effectively neutral (defined as 0 < N(e)s < 1) and that this fraction was negatively correlated with a gene's expression level. Consistent with models of recurrent adaptive evolution, we detected a negative correlation between levels of synonymous site polymorphism and the rate of protein evolution, although the correlation was weak and nonsignificant. No systematic X chromosome-autosome difference was found in the efficacy of selection. For example, the proportion of adaptive substitutions was significantly higher on the X chromosome compared with the autosomes in O. c. algirus but not in O. c. cuniculus. Our findings support widespread positive and purifying selection in rabbits and add to a growing list of examples suggesting that differences in N(e) among taxa play a substantial role in determining rates and patterns of protein evolution.

Building bridges: formin1 of Arabidopsis forms a connection between the cell wall and the actin cytoskeleton.

Actin microfilament (MF) organization and remodelling is critical to cell function. The formin family of actin binding proteins are involved in nucleating MFs in Arabidopsis thaliana. They all contain formin homology domains in the intracellular, C-terminal half of the protein that interacts with MFs. Formins in class I are usually targeted to the plasma membrane and this is true of Formin1 (AtFH1) of A. thaliana. In this study, we have investigated the extracellular domain of AtFH1 and we demonstrate that AtFH1 forms a bridge from the actin cytoskeleton, across the plasma membrane and is anchored within the cell wall. AtFH1 has a large, extracellular domain that is maintained by purifying selection and that contains four conserved regions, one of which is responsible for immobilising the protein. Protein anchoring within the cell wall is reduced in constructs that express truncations of the extracellular domain and in experiments in protoplasts without primary cell walls. The 18 amino acid proline-rich extracellular domain that is responsible for AtFH1 anchoring has homology with cell-wall extensins. We also have shown that anchoring of AtFH1 in the cell wall promotes actin bundling within the cell and that overexpression of AtFH1 has an inhibitory effect on organelle actin-dependant dynamics. Thus, the AtFH1 bridge provides stable anchor points for the actin cytoskeleton and is probably a crucial component of the signalling response and actin-remodelling mechanisms.

Genomic and evolutionary features of the SPI-1 type III secretion system that is present in Xanthomonas albilineans but is not essential for xylem colonization and symptom development of sugarcane leaf scald.

Xanthomonas albilineans is the causal agent of sugarcane leaf scald. Interestingly, this bacterium, which is not known to be insect or animal associated, possesses a type III secretion system (T3SS) belonging to the injectisome family Salmonella pathogenicity island 1 (SPI-1). The T3SS SPI-1 of X. albilineans shares only low similarity with other available T3SS SPI-1 sequences. Screening of a collection of 128 plant-pathogenic bacteria revealed that this T3SS SPI-1 is present in only two species of Xanthomonas: X. albilineans and X. axonopodis pv. phaseoli. Inoculation of sugarcane with knockout mutants showed that this system is not required by X. albilineans to spread within xylem vessels and to cause disease symptoms. This result was confirmed by the absence of this T3SS SPI-1 in an X. albilineans strain isolated from diseased sugarcane. To investigate the importance of the T3SS SPI-1 during the life cycle of X. albilineans, we analyzed T3SS SPI-1 sequences from 11 strains spanning the genetic diversity of this species. No nonsense mutations or frameshifting indels were observed in any of these strains, suggesting that the T3SS SPI-1 system is maintained within the species X. albilineans. Evolutionary features of T3SS SPI-1 based on phylogenetic, recombination, and selection analyses are discussed in the context of the possible functional importance of T3SS SPI-1 in the ecology of X. albilineans.

Evolution of endogenous sequences of banana streak virus: what can we learn from banana (Musa sp.) evolution?

Endogenous plant pararetroviruses (EPRVs) are viral sequences of the family Caulimoviridae integrated into the nuclear genome of numerous plant species. The ability of some endogenous sequences of Banana streak viruses (eBSVs) in the genome of banana (Musa sp.) to induce infections just like the virus itself was recently demonstrated (P. Gayral et al., J. Virol. 83:6697-6710, 2008). Although eBSVs probably arose from accidental events, infectious eBSVs constitute an extreme case of parasitism, as well as a newly described strategy for vertical virus transmission in plants. We investigated the early evolutionary stages of infectious eBSV for two distinct BSV species-GF (BSGFV) and Imové (BSImV)-through the study of their distribution, insertion polymorphism, and structure evolution among selected banana genotypes representative of the diversity of 60 wild Musa species and genotypes. To do so, the historical frame of host evolution was analyzed by inferring banana phylogeny from two chloroplast regions-matK and trnL-trnF-as well as from the nuclear genome, using 19 microsatellite loci. We demonstrated that both BSV species integrated recently in banana evolution, circa 640,000 years ago. The two infectious eBSVs were subjected to different selective pressures and showed distinct levels of rearrangement within their final structure. In addition, the molecular phylogenies of integrated and nonintegrated BSVs enabled us to establish the phylogenetic origins of eBSGFV and eBSImV.

Chromosomal scale assembly of parasitic wasp genome reveals symbiotic virus colonization.

Endogenous viruses form an important proportion of eukaryote genomes and a source of novel functions. How large DNA viruses integrated into a genome evolve when they confer a benefit to their host, however, remains unknown. Bracoviruses are essential for the parasitism success of parasitoid wasps, into whose genomes they integrated ~103 million years ago. Here we show, from the assembly of a parasitoid wasp genome at a chromosomal scale, that bracovirus genes colonized all ten chromosomes of Cotesia congregata. Most form clusters of genes involved in particle production or parasitism success. Genomic comparison with another wasp, Microplitis demolitor, revealed that these clusters were already established ~53 mya and thus belong to remarkably stable genomic structures, the architectures of which are evolutionary constrained. Transcriptomic analyses highlight temporal synchronization of viral gene expression without resulting in immune gene induction, suggesting that no conflicts remain between ancient symbiotic partners when benefits to them converge.

A four-partner plant­virus interaction: enemies can also come from within.

Plant viruses are disseminated by either vertical (vegetative multiplication or sexual reproduction) or horizontal (vector-mediated) propagation. Plant pararetroviruses—members of the Caulimoviridae family—have developed an alternative strategy for vertical propagation via integration within the host plant genome, although integration is not required for viral replication. Integrated endogenous pararetrovirus (EPRV) sequences have undergone extensive viral genome rearrangements and contain more than one copy of the viral genome. Furthermore, EPRV can become infectious upon spontaneous escape of active virus following stresses such as wounding, tissue culture, or interspecific crosses. Such infectious EPRV are of great importance, not only in terms of their ability to precipitate epidemic outbreaks but also because of their effect on breeding of numerous plant genomes in temperate and tropical crops. This is especially true for banana, a crop susceptible to banana streak viruses, the causative agents of banana streak disease. Thus, the classical three-component banana–Banana streak virus (BSV)–mealybug pathosystem can be expanded to include endogenous BSV as an alternative source of active virions. The BSV-banana pathosystem is one of only three pathosystems known to date to harbor this remarkable feature, and the present review focuses exclusively on it to illustrate this four-partner interaction.

Phylogeny of Banana Streak Virus reveals recent and repetitive endogenization in the genome of its banana host (Musa sp.).

Banana streak virus (BSV) is a plant dsDNA pararetrovirus (family Caulimoviridae, genus badnavirus). Although integration is not an essential step in the BSV replication cycle, the nuclear genome of banana (Musa sp.) contains BSV endogenous pararetrovirus sequences (BSV EPRVs). Some BSV EPRVs are infectious by reconstituting a functional viral genome. Recent studies revealed a large molecular diversity of episomal BSV viruses (i.e., nonintegrated) while others focused on BSV EPRV sequences only. In this study, the evolutionary history of badnavirus integration in banana was inferred from phylogenetic relationships between BSV and BSV EPRVs. The relative evolution rates and selective pressures (d(N)/d(S) ratio) were also compared between endogenous and episomal viral sequences. At least 27 recent independent integration events occurred after the divergence of three banana species, indicating that viral integration is a recent and frequent phenomenon. Relaxation of selective pressure on badnaviral sequences that experienced neutral evolution after integration in the plant genome was recorded. Additionally, a significant decrease (35%) in the EPRV evolution rate was observed compared to BSV, reflecting the difference in the evolution rate between episomal dsDNA viruses and plant genome. The comparison of our results with the evolution rate of the Musa genome and other reverse-transcribing viruses suggests that EPRVs play an active role in episomal BSV diversity and evolution.

A single Banana streak virus integration event in the banana genome as the origin of infectious endogenous pararetrovirus.

Sequencing of plant nuclear genomes reveals the widespread presence of integrated viral sequences known as endogenous pararetroviruses (EPRVs). Banana is one of the three plant species known to harbor infectious EPRVs. Musa balbisiana carries integrated copies of Banana streak virus (BSV), which are infectious by releasing virions in interspecific hybrids. Here, we analyze the organization of the EPRV of BSV Goldfinger (BSGfV) present in the wild diploid M. balbisiana cv. Pisang Klutuk Wulung (PKW) revealed by the study of Musa bacterial artificial chromosome resources and interspecific genetic cross. cv. PKW contains two similar EPRVs of BSGfV. Genotyping of these integrants and studies of their segregation pattern show an allelic insertion. Despite the fact that integrated BSGfV has undergone extensive rearrangement, both EPRVs contain the full-length viral genome. The high degree of sequence conservation between the integrated and episomal form of the virus indicates a recent integration event; however, only one allele is infectious. Analysis of BSGfV EPRV segregation among an F1 population from an interspecific genetic cross revealed that these EPRV sequences correspond to two alleles originating from a single integration event. We describe here for the first time the full genomic and genetic organization of the two EPRVs of BSGfV present in cv. PKW in response to the challenge facing both scientists and breeders to identify and generate genetic resources free from BSV. We discuss the consequences of this unique host-pathogen interaction in terms of genetic and genomic plant defenses versus strategies of infectious BSGfV EPRVs.

Viruses in the Invasive Hornet Vespa velutina.

The Asian yellow-legged hornet Vespa velutina nigrithorax, a major predator of honeybees, is spreading in Europe in part due to a lack of efficient control methods. In this study, as a first step to identify biological control agents, we characterized viral RNA sequences present in asymptomatic or symptomatic hornets. Among 19 detected viruses, the honey bee virus Deformed wing virus-B was predominant in all the samples, particularly in muscles from the symptomatic hornet, suggesting a putative cause of the deformed wing symptom. Interestingly, two new viruses closely related to Acyrthosiphon pisumvirus and Himetobi Pvirus and viruses typically associated with honey bees, Acute bee paralysis virus and Black queen cell virus, were detected in the brain and muscles, and may correspond to the circulation and possible replication forms of these viruses in the hornet. Aphid lethal paralysis virus, Bee Macula-like virus, and Moku virus, which are known to infect honey bees, were also identified in the gut virus metagenome of hornets. Therefore, our study underlined the urgent need to study the host range of these newly discovered viruses in hornets to determine whether they represent a new threat for honey bees or a hope for the biocontrol of V. velutina.

Gall Wasp Transcriptomes Unravel Potential Effectors Involved in Molecular Dialogues With Oak and Rose.

To gain insight into wasp factors that might be involved in the initial induction of galls on woody plants, we performed high throughput (454) transcriptome analysis of ovaries and venom glands of two cynipid gall wasps, Biorhiza pallida and Diplolepis rosae, inducing galls on oak and rose, respectively. De novo assembled and annotated contigs were compared to sequences from phylogenetically related parasitoid wasps. The relative expression levels of contigs were estimated to identify the most expressed gene sequences in each tissue. We identify for the first time a set of maternally expressed gall wasp proteins potentially involved in the interaction with the plant. Some genes highly expressed in venom glands and ovaries may act to suppress early plant defense signaling. We also identify gall wasp cellulases that could be involved in observed local lysis of plant tissue following oviposition, and which may have been acquired from bacteria by horizontal gene transfer. We find no evidence of virus-related gene expression, in contrast to many non-cynipid parasitoid wasps. By exploring gall wasp effectors, this study is a first step toward understanding the molecular mechanisms underlying cynipid gall induction in woody plants, and the recent sequencing of oak and rose genomes will enable study of plant responses to these factors.

Occurrence and Molecular Phylogeny of Honey Bee Viruses in Vespids.

Since the discovery that honey bee viruses play a role in colony decline, researchers have made major breakthroughs in understanding viral pathology and infection processes in honey bees. Work on virus transmission patterns and virus vectors, such as the mite Varroa destructor, has prompted intense efforts to manage honey bee health. However, little is known about the occurrence of honey bee viruses in bee predators, such as vespids. In this study, we characterized the occurrence of 11 honey bee viruses in five vespid species and one wasp from four provinces in China and two vespid species from four locations in France. The results showed that all the species from China carried certain honey bee viruses, notably Apis mellifera filamentous virus (AmFV), Deformed wing virus (DWV), and Israeli acute paralysis virus (IAPV); furthermore, in some vespid colonies, more than three different viruses were identified. In France, DWV was the most common virus; Sacbrood virus (SBV) and Black queen cell virus (BQCV) were observed in one and two samples, respectively. Phylogenetic analyses of IAPV and BQCV sequences indicated that most of the IAPV sequences belonged to a single group, while the BQCV sequences belonged to several groups. Additionally, our study is the first to detect Lake Sinai virus (LSV) in a hornet from China. Our findings can guide further research into the origin and transmission of honey bee viruses in Vespidae, a taxon of ecological, and potentially epidemiological, relevance.

Discovery of Culex pipiens associated tunisia virus: a new ssRNA(+) virus representing a new insect associated virus family.

In the global context of arboviral emergence, deep sequencing unlocks the discovery of new mosquito-borne viruses. Mosquitoes of the species Culex pipiens, C. torrentium, and C. hortensis were sampled from 22 locations worldwide for transcriptomic analyses. A virus discovery pipeline was used to analyze the dataset of 0.7 billion reads comprising 22 individual transcriptomes. Two closely related 6.8 kb viral genomes were identified in C. pipiens and named as Culex pipiens associated tunisia virus (CpATV) strains Ayed and Jedaida. The CpATV genome contained four ORFs. ORF1 possessed helicase and RNA-dependent RNA polymerase (RdRp) domains related to new viral sequences recently found mainly in dipterans. ORF2 and 4 contained a capsid protein domain showing strong homology with Virgaviridae plant viruses. ORF3 displayed similarities with eukaryotic Rhoptry domain and a merozoite surface protein (MSP7) domain only found in mosquito-transmitted Plasmodium, suggesting possible interactions between CpATV and vertebrate cells. Estimation of a strong purifying selection exerted on each ORFs and the presence of a polymorphism maintained in the coding region of ORF3 suggested that both CpATV sequences are genuine functional viruses. CpATV is part of an entirely new and highly diversified group of viruses recently found in insects, and that bears the genomic hallmarks of a new viral family.

Comparative Analysis of Transcriptomes from Secondary Reproductives of Three Reticulitermes Termite Species.

Termites are eusocial insects related to cockroaches that feed on lignocellulose. These insects are key species in ecosystems since they recycle a large amount of nutrients but also are pests, exerting major economic impacts. Knowledge on the molecular pathways underlying reproduction, caste differentiation or lignocellulose digestion would largely benefit from additional transcriptomic data. This study focused on transcriptomes of secondary reproductive females (nymphoid neotenics). Thirteen transcriptomes were used: 10 of Reticulitermes flavipes and R. grassei sequenced from a previous study, and two transcriptomes of R. lucifugus sequenced for the present study. After transcriptome assembly and read mapping, we examined interspecific variations of genes expressed by termites or gut microorganisms. A total of 18,323 orthologous gene clusters were detected. Functional annotation and taxonomic assignment were performed on a total of 41,287 predicted contigs in the three termite species. Between the termite species studied, functional categories of genes were comparable. Gene ontology (GO) terms analysis allowed the discovery of 9 cellulases and a total of 79 contigs potentially involved in 11 enzymatic activities used in wood metabolism. Altogether, results of this study illustrate the strong potential for the use of comparative interspecific transcriptomes, representing a complete resource for future studies including differentially expressed genes between castes or SNP analysis for population genetics.

Adaptive selection on bracovirus genomes drives the specialization of Cotesia parasitoid wasps.

The geographic mosaic of coevolution predicts parasite virulence should be locally adapted to the host community. Cotesia parasitoid wasps adapt to local lepidopteran species possibly through their symbiotic bracovirus. The virus, essential for the parasitism success, is at the heart of the complex coevolutionary relationship linking the wasps and their hosts. The large segmented genome contained in the virus particles encodes virulence genes involved in host immune and developmental suppression. Coevolutionary arms race should result in the positive selection of particular beneficial alleles. To understand the global role of bracoviruses in the local adaptation or specialization of parasitoid wasps to their hosts, we studied the molecular evolution of four bracoviruses associated with wasps of the genus Cotesia, including C congregata, C vestalis and new data and annotation on two ecologically differentiated populations of C sesamie, Kitale and Mombasa. Paired orthologs analyses revealed more genes under positive selection when comparing the two C sesamiae bracoviruses belonging to the same species, and more genes under strong evolutionary constraint between species. Furthermore branch-site evolutionary models showed that 17 genes, out of the 54 currently available shared by the four bracoviruses, harboured sites under positive selection including: the histone H4-like, a C-type lectin, two ep1-like, ep2, a viral ankyrin, CrV1, a ben-domain, a Serine-rich, and eight unknown genes. Lastly the phylogenetic analyses of the histone, ep2 and CrV1 genes in different African C sesamiae populations showed that each gene described differently the individual relationships. In particular we found recombination had happened between the ep2 and CrV1 genes, which are localized 37.5 kb apart on the wasp chromosomes. Involved in multidirectional coevolutionary interactions, C sesamiae wasps rely on different bracovirus mediated molecular pathways to overcome local host resistance.