Reconstructing the global migration history of Phytophthora infestans towards Colombia.

The evolution of new variants of plant pathogens is one of the biggest challenges to controlling and managing plant diseases. Of the forces driving these evolutionary processes, global migration events are particularly important for widely distributed diseases such as the potato late blight, caused by the oomycete Phytophthora infestans (Mont.) de Bary. However, little is known about its migration routes outside North America and Europe. This work used genotypic data from population studies to elucidate the migration history originating the Colombian P. infestans population. For this purpose, a dataset of 1706 P. infestans genotypes was recollected, representing North and South America, Europe, and Asia. Descriptive analysis and historical records from North America and Europe were used to propose three global migration hypotheses, differing on the origin of the disease (Mexico or Peru) and the hypothesis that it came back to South America from Europe. These scenarios were tested using Approximate Bayesian Computation. According to this analysis, the most probable scenario (posterior probability = 0.631) was the one proposing a Peruvian origin for P. infestans, an initial migration towards Colombia and Mexico, and a later event from Mexico to the United States, and then to Europe and Asia, with no return to northern South America. In Colombia, the scenario considering a single migration from Peru and posterior migrations within Colombia was the most probable with a posterior probability of 0.640. The obtained results support the hypothesis of a Peruvian origin for P. infestans followed by rare colonization events worldwide.

Quantitative trait loci involved in the reproductive success of a parasitoid wasp.

Dissecting the genetic basis of intraspecific variations in life history traits is essential to understand their evolution, notably for potential biocontrol agents. Such variations are observed in the endoparasitoid Cotesia typhae (Hymenoptera: Braconidae), specialized on the pest Sesamia nonagrioides (Lepidoptera: Noctuidae). Previously, we identified two strains of C. typhae that differed significantly for life history traits on an allopatric host population. To investigate the genetic basis underlying these phenotypic differences, we used a quantitative trait locus (QTL) approach based on restriction site-associated DNA markers. The characteristic of C. typhae reproduction allowed us generating sisters sharing almost the same genetic content, named clonal sibship. Crosses between individuals from the two strains were performed to generate F2 and F8 recombinant CSS. The genotypes of 181 clonal sibships were determined as well as the phenotypes of the corresponding 4,000 females. Informative markers were then used to build a high-quality genetic map. These 465 markers spanned a total length of 1,300 cM and were organized in 10 linkage groups which corresponded to the number of C. typhae chromosomes. Three QTLs were detected for parasitism success and two for offspring number, while none were identified for sex ratio. The QTLs explained, respectively, 27.7% and 24.5% of the phenotypic variation observed. The gene content of the genomic intervals was investigated based on the genome of C. congregata and revealed 67 interesting candidates, as potentially involved in the studied traits, including components of the venom and of the symbiotic virus (bracovirus) shown to be necessary for parasitism success in related wasps.

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.

Critical limb ischaemia and the response to bone marrow-derived cell therapy according to tcPO2 measurement.

BACKGROUND: Cell therapy is an emerging potential biotherapy for critical limb ischaemia (CLI) patients who are not eligible for revascularization. However, the findings on this technique's efficacy are inconsistent. Trials investigating this topic focused on the more severe CLI patients who were often beyond any therapy. Therefore, identifying those who may truly benefit from cell transplantation is now warranted. To this end, we studied the prognostic value of tcPO2 for major amputation after 1 year in patients treated with bone marrow-derived cells. PATIENTS AND METHODS: CLI patients ineligible for revascularization were included in a cell-therapy pilot study. On inclusion, patients underwent tcPO2 measurement in supine and sitting positions. For a tcPO2 < 10 mmHg in the supine position, the vascular reserve was defined by tcPO2 > 30 mmHg in the sitting position. Patients were administered intramuscular injections of mononuclear cells derived from aspirated bone marrow. RESULTS: In total, 25 patients (a lower limbs) were included for analysis. At inclusion, 11 lower limbs had tcPO2 at rest > 10 mmHg, and 16 lower limbs had a tcPO2 < 10 mmHg. The success probability for cell therapy was 0.79 (95 % CI 0.38-0.94) and 0.44 (95 % CI 0.18-0.67), respectively (p = 0.1). Of the 16 limbs with tcPO2 < 10 mmHg, the success rate was considerably higher in patients demonstrating a tcPO2 increase in a sitting position of over 30 mmHg (6/8, success probability 0.71, 95 % CI 0.26-0.92) compared to those without (2/8, success probability 0.15, 95 % CI 0.01-0.48, p = 0.03). CONCLUSIONS: For patients with chronic CLI for whom cellular therapy is a therapeutic option, a tcPO2 < 10 mmHg at rest, without vascular reserve (i. e. < 30 mmHg when sitting), is a prognostic indicator for poor outcome.
.

High-NaCl perception in Drosophila melanogaster.

Salt is a fundamental nutrient that is required for many physiological processes, including electrolyte homeostasis and neuronal activity. In mammals and Drosophila, the detection of NaCl induces two different behaviors: low-salt concentrations provide an attractive stimulus, whereas high-salt concentrations are avoided. We identified the gene called serrano (sano) as being expressed in the sensory organs of Drosophila larvae. A transgenic reporter line showed that sano was coexpressed with Gr66a in a subset of gustatory neurons in the terminal organ of third-instar larvae. The disruption of sano gene expression in gustatory neurons led to the specific loss of high-salt concentration avoidance in larvae, whereas the detection of other attractive or aversive substances was unaffected. Moreover, using a cellular marker sensitive to calcium levels, Sano function was shown to be required for neuronal activity in response to high-salt concentrations. In these neurons, the loss of the DEG/ENaC channel PPK19 function also eliminated the cellular response to high-salt concentrations. Our study revealed that PPK19 and Sano are required in the neurons of the larval gustatory organs for the detection of high-salt concentrations.

Expression of a desaturase gene, desat1, in neural and nonneural tissues separately affects perception and emission of sex pheromones in Drosophila.

Animals often use sex pheromones for mate choice and reproduction. As for other signals, the genetic control of the emission and perception of sex pheromones must be tightly coadapted, and yet we still have no worked-out example of how these two aspects interact. Most models suggest that emission and perception rely on separate genetic control. We have identified a Drosophila melanogaster gene, desat1, that is involved in both the emission and the perception of sex pheromones. To explore the mechanism whereby these two aspects of communication interact, we investigated the relationship between the molecular structure, tissue-specific expression, and pheromonal phenotypes of desat1. We characterized the five desat1 transcripts-all of which yielded the same desaturase protein-and constructed transgenes with the different desat1 putative regulatory regions. Each region was used to target reporter transgenes with either (i) the fluorescent GFP marker to reveal desat1 tissue expression, or (ii) the desat1 RNAi sequence to determine the effects of genetic down-regulation on pheromonal phenotypes. We found that desat1 is expressed in a variety of neural and nonneural tissues, most of which are involved in reproductive functions. Our results suggest that distinct desat1 putative regulatory regions independently drive the expression in nonneural and in neural cells, such that the emission and perception of sex pheromones are precisely coordinated in this species.

Development of a viral biopesticide for the control of the Guatemala potato tuber moth Tecia solanivora.

The Guatemala potato tuber moth Tecia solanivora (Povolny) (Lep. Gelechiidae) is an invasive species from Mesoamerica that has considerably extended its distribution area in recent decades. While this species is considered to be a major potato pest in Venezuela, Colombia, and Ecuador, currently no specific control methods are available for farmers. To address this issue we developed a biopesticide formulation to be used in integrated pest management of T. solanivora, following three steps. First, search for entomopathogenic viruses were carried out through extensive bioprospections in 12 countries worldwide. As a result, new Phthorimaea operculella granulovirus (PhopGV) isolates were found in T. solanivora and five other gelechid species. Second, twenty PhopGV isolates, including both previously known and newly found isolates, were genetically and/or biologically characterized in order to choose the best candidate for a biopesticide formulation. Sequence data were obtained for the ecdysteroid UDP-glucosyltransferase (egt) gene, a single copy gene known to play a role in pathogenicity. Three different sizes (1086, 1305 and 1353 bp) of egt were found among the virus isolates analyzed. Unexpectedly, no obvious correlation between egt size and pathogenicity was found. Bioassays on T. solanivora neonates showed a maximum of a 14-fold difference in pathogenicity among the eight PhopGV isolates tested. The most pathogenic PhopGV isolate, JLZ9f, had a medium lethal concentration (LC(50)) of 10 viral occlusion bodies per square mm of consumed tuber skin. Third, we tested biopesticide dust formulations by mixing a dry carrier (calcium carbonate) with different adjuvants (magnesium chloride or an optical brightener or soya lecithin) and different specific amounts of JLZ9f. During laboratory experiments, satisfactory control of the pest (>98% larva mortality compared to untreated control) was achieved with a formulation containing 10 macerated JLZ9f-dead T. solanivora larvae per kg of calcium carbonate mixed with 50 mL/kg of soya lecithin. The final product provides an interesting alternative to chemical pesticides for Andean farmers affected by this potato pest.

A minisatellite-based MLVA for deciphering the global epidemiology of the bacterial cassava pathogen Xanthomonas phaseoli pv. manihotis.

Cassava Bacterial Blight (CBB) is a destructive disease widely distributed in the different areas where this crop is grown. Populations studies have been performed at local and national scales revealing a geographical genetic structure with temporal variations. A global epidemiology analysis of its causal agent Xanthomonas phaseoli pv. manihotis (Xpm) is needed to better understand the expansion of the disease for improving the monitoring of CBB. We targeted new tandem repeat (TR) loci with large repeat units, i.e. minisatellites, that we multiplexed in a scheme of Multi-Locus Variable number of TR Analysis (MLVA-8). This genotyping scheme separated 31 multilocus haplotypes in three clusters of single-locus variants and a singleton within a worldwide collection of 93 Xpm strains isolated over a period of fifty years. The major MLVA-8 cluster 1 grouped strains originating from all countries, except the unique Chinese strain. On the contrary, all the Xpm strains genotyped using the previously developed MLVA-14 microsatellite scheme were separated as unique haplotypes. We further propose an MLVA-12 scheme which takes advantage of combining TR loci with different mutation rates: the eight minisatellites and four faster evolving microsatellite markers, for global epidemiological surveillance. This MLVA-12 scheme identified 78 haplotypes and separated most of the strains in groups of double-locus variants (DLV) supporting some phylogenetic relationships. DLV groups were subdivided into closely related clusters of strains most often sharing the same geographical origin and isolated over a short period, supporting epidemiological relationships. The main MLVA-12 DLV group#1 was composed by strains from South America and all the African strains. The MLVA-12 scheme combining both minisatellite and microsatellite loci with different discriminatory power is expected to increase the accuracy of the phylogenetic signal and to minimize the homoplasy effects. Further investigation of the global epidemiology of Xpm will be helpful for a better control of CBB worldwide.

Evolutionary mechanisms driving the evolution of a large polydnavirus gene family coding for protein tyrosine phosphatases.

BACKGROUND: Gene duplications have been proposed to be the main mechanism involved in genome evolution and in acquisition of new functions. Polydnaviruses (PDVs), symbiotic viruses associated with parasitoid wasps, are ideal model systems to study mechanisms of gene duplications given that PDV genomes consist of virulence genes organized into multigene families. In these systems the viral genome is integrated in a wasp chromosome as a provirus and virus particles containing circular double-stranded DNA are injected into the parasitoids' hosts and are essential for parasitism success. The viral virulence factors, organized in gene families, are required collectively to induce host immune suppression and developmental arrest. The gene family which encodes protein tyrosine phosphatases (PTPs) has undergone spectacular expansion in several PDV genomes with up to 42 genes. RESULTS: Here, we present strong indications that PTP gene family expansion occurred via classical mechanisms: by duplication of large segments of the chromosomally integrated form of the virus sequences (segmental duplication), by tandem duplications within this form and by dispersed duplications. We also propose a novel duplication mechanism specific to PDVs that involves viral circle reintegration into the wasp genome. The PTP copies produced were shown to undergo conservative evolution along with episodes of adaptive evolution. In particular recently produced copies have undergone positive selection in sites most likely involved in defining substrate selectivity. CONCLUSION: The results provide evidence about the dynamic nature of polydnavirus proviral genomes. Classical and PDV-specific duplication mechanisms have been involved in the production of new gene copies. Selection pressures associated with antagonistic interactions with parasitized hosts have shaped these genes used to manipulate lepidopteran physiology with evidence for positive selection involved in adaptation to host targets.

Intraspecific specialization of the generalist parasitoid Cotesia sesamiae revealed by polyDNAvirus polymorphism and associated with different Wolbachia infection.

As a result of an intense host-parasite evolutionary arms race, parasitic wasps frequently display high levels of specialization on very few host species. For instance, in braconid wasps very few generalist species have been described. However, within this family, Cotesia sesamiae is a generalist species that is widespread in sub-Saharan Africa and develops on several lepidopteran hosts. In this study, we tested the hypothesis that C. sesamiae may be a cryptic specialist when examined at the intraspecific level. We sequenced exon 2 of CrV1, a gene of the symbiotic polyDNAvirus that is integrated into the wasp genome and is associated with host immune suppression. We found that CrV1 genotype was more closely associated with the host in which the parasitoid developed than any abiotic environmental factor tested. We also tested a correlation between CrV1 genotype and an infection with Wolbachia bacteria, which are known for their ability to induce reproductive isolation. The Wolbachia bacteria infection polymorphism was also found as a major factor explaining the genetic structure of CrV1, and, in addition, the best model explaining CrV1 genetic structure involved an interaction between Wolbachia infection and host species. We suggest that Wolbachia could act as an agent capable of maintaining advantageous alleles for host specialization in different populations of C. sesamiae. This mechanism could be applicable to other insect models because of the high prevalence of Wolbachia in insects.

Parasitism of lepidopterous stem borers in cultivated and natural habitats.

Plant infestation, stem borer density, parasitism, and parasitoid abundance were assessed during two years in two host plants, Zea mays (L.) (Cyperales: Poaceae) and Sorghum bicolor (L.) (Cyperales: Poaceae), in cultivated habitats. The four major host plants (Cyperus spp., Panicum spp., Pennisetum spp., and Sorghum spp.) found in natural habitats were also assessed, and both the cultivated and natural habitat species occurred in four agroecological zones in Kenya. Across habitats, plant infestation (23.2%), stem borer density (2.2 per plant), and larval parasitism (15.0%) were highest in maize in cultivated habitats. Pupal parasitism was not higher than 4.7% in both habitats, and did not vary with locality during each season or with host plant between each season. Cotesia sesamiae (Cameron) and C. flavipes Cameron (Hymenoptera: Braconidae) were the key parasitoids in cultivated habitats (both species accounted for 76.4% of parasitized stem borers in cereal crops), but not in natural habitats (the two Cotesia species accounted for 14.5% of parasitized stem borers in wild host plants). No single parasitoid species exerted high parasitism rates on stem borer populations in wild host plants. Low stem borer densities across seasons in natural habitats indicate that cereal stem borer pests do not necessarily survive the non-cropping season feeding actively in wild host plants. Although natural habitats provided refuges for some parasitoid species, stem borer parasitism was generally low in wild host plants. Overall, because parasitoids contribute little in reducing cereal stem borer pest populations in cultivated habitats, there is need to further enhance their effectiveness in the field to regulate these pests.

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.

Collisional mechanism of ligand release by Bombyxmori JHBP, a member of the TULIP / Takeout family of lipid transporters.

Juvenile hormones (JHs) regulate important processes in insects, such as postembryonic development and reproduction. In the hemolymph of Lepidoptera, these lipophilic sesquiterpenic hormones are transported from their site of synthesis to target tissues by high affinity carriers, the juvenile hormone binding proteins (JHBPs). Lepidopteran JHBPs belong to a recently uncovered, yet very ancient family of proteins sharing a common lipid fold (TULIP domain) and involved in shuttling various lipid ligands. One important, but poorly understood aspect of JHs action, is the mechanism of hormone transfer to or through the plasma membranes of target cells. Since many membrane-active peptides and proteins, such as the pore-forming bacterial toxins, are activated by low pH or interaction with phospholipid membranes, we have examined the effect of these factors on JH binding by JHBPs. The affinity of Bombyx mori and Manduca sexta JHBPs for JH III was determined by the DCC assay, equilibrium dialysis, and isothermal titration calorimetry, and found to be greatly reduced at low pH, in agreement with previous observations. Loss of binding was accompanied by changes in fluorescence and near-UV CD spectra, indicating significant changes in protein structure in the environment of aromatic residues. The apparent dissociation rate constant (koff) of the JHBP-JH III complex was greater at acidic pH, suggesting that low pH favors ligand release by opening of the binding pocket. The affinity of recombinant B. mori JHBP (rBmJHBP) was also decreased in the presence of anionic phospholipid vesicles. Measurements of steady-state fluorescence anisotropy with the lipophilic probe TMA-DPH demonstrated that rBmJHBP specifically interacts with anionic membranes. These results suggest the existence of a collisional mechanism for ligand release that may be important for delivery of JHs to the target cells, and could be relevant to the function of related members of this emerging family of lipid-transport proteins.

Maintenance of adaptive differentiation by Wolbachia induced bidirectional cytoplasmic incompatibility: the importance of sib-mating and genetic systems.

BACKGROUND: Bacteria of the genus Wolbachia are reproductive parasites widespread among arthropods. The most common effect arising from the presence of Wolbachia in a population is Cytoplasmic Incompatibility (CI), whereby postmating reproductive isolation occurs in crosses between an infected male and an uninfected female, or when a male is infected with a different strain of Wolbachia to that of the female (bidirectional CI). Previous theoretical models have demonstrated that bidirectional CI can contribute to the genetic divergence of populations in haploid and diploid organisms. However, haplodiploid organisms were not considered in these models even though they include Nasonia parasitoid wasps - the best example of the implication of Wolbachia in ongoing speciation. Moreover, previous work did not investigate inbreeding mating systems, which are frequently observed in arthropod species. RESULTS: We developed a stochastic two-island model which simulated three genetic scenarios, diploidy, haploidy, and haplodiploidy, with two CI phenotypes being considered for the latter: (1) male development of female progeny; and (2) mortality of fertilized eggs. We also investigated the effect of varying the proportion of sib mating. In the model each allopatric population was initially fixed for a single allele at a nuclear locus under positive selection and infected with one strain of Wolbachia. Each simulation presupposed that the two populations were fixed for a different allele and a different strain of Wolbachia. The degree of genetic differentiation observed in the locus under selection due to bidirectional CI was much lower for the two haplodiploid phenotypes than for either diploids or haploids. Furthermore, we demonstrated that sib-mating may compensate for the lower efficiency of bidirectional CI in haplodiploids by maintaining genetic divergence. CONCLUSION: Our model suggests that maintenance of genetic differentiation facilitated by Wolbachia is more likely to occur in diploids and haploids than in haplodiploids. However, increasing the level of sib-mating may compensate for the weak effect of bidirectional CI in haplodiploids. Our work therefore gives a potential explanation for why the haplodiploid Nasonia species, which are infected with bidirectionally incompatible Wolbachia strains and undergo sib-mating, have differentiated genetically and maintained this differentiation without premating isolation.

Specific interaction between Sam68 and neuronal mRNAs: implication for the activity-dependent biosynthesis of elongation factor eEF1A.

In cultured hippocampal neurons and in adult brain, the splicing regulatory protein Sam68 is partially relocated to the somatodendritic domain and associates with dendritic polysomes. Transfer to the dendrites is activity-dependent. We have investigated the repertoire of neuronal mRNAs to which Sam68 binds in vivo. By using coimmunoprecipitation and microarray screening techniques, Sam68 was found to associate with a number of plasticity-related mRNA species, including Eef1a1, an activity-responsive mRNA coding for translation elongation factor eEF1A. In cortical neuronal cultures, translation of the Eef1a1 mRNA was strongly induced by neuronal depolarisation and correlated with enhanced association of Sam68 with polysomal mRNAs. The possible function of Sam68 in Eef1a1 mRNA utilization was studied by expressing a dominant-negative, cytoplasmic Sam68 mutant (GFP-Sam68DeltaC) in cultured hippocampal neurons. The level of eEF1A was lower in neurons expressing GFP-Sam68DeltaC than in control neurons, supporting the proposal that endogenous Sam68 may contribute to the translational efficiency of the Eef1a1 mRNA. These findings are discussed in the light of the complex, potentially crucial regulation of eEF1A biosynthesis during long-term synaptic change.

Viral cystatin evolution and three-dimensional structure modelling: a case of directional selection acting on a viral protein involved in a host-parasitoid interaction.

BACKGROUND: In pathogens, certain genes encoding proteins that directly interact with host defences coevolve with their host and are subject to positive selection. In the lepidopteran host-wasp parasitoid system, one of the most original strategies developed by the wasps to defeat host defences is the injection of a symbiotic polydnavirus at the same time as the wasp eggs. The virus is essential for wasp parasitism success since viral gene expression alters the immune system and development of the host. As a wasp mutualist symbiont, the virus is expected to exhibit a reduction in genome complexity and evolve under wasp phyletic constraints. However, as a lepidopteran host pathogenic symbiont, the virus is likely undergoing strong selective pressures for the acquisition of new functions by gene acquisition or duplication. To understand the constraints imposed by this particular system on virus evolution, we studied a polydnavirus gene family encoding cyteine protease inhibitors of the cystatin superfamily. RESULTS: We show that cystatins are the first bracovirus genes proven to be subject to strong positive selection within a host-parasitoid system. A generated three-dimensional model of Cotesia congregata bracovirus cystatin 1 provides a powerful framework to position positively selected residues and reveal that they are concentrated in the vicinity of actives sites which interact with cysteine proteases directly. In addition, phylogenetic analyses reveal two different cystatin forms which evolved under different selective constraints and are characterized by independent adaptive duplication events. CONCLUSION: Positive selection acts to maintain cystatin gene duplications and induces directional divergence presumably to ensure the presence of efficient and adapted cystatin forms. Directional selection has acted on key cystatin active sites, suggesting that cystatins coevolve with their host target. We can strongly suggest that cystatins constitute major virulence factors, as was already proposed in previous functional studies.

The Quetzal Coalescence template library: A C++ programmers resource for integrating distributional, demographic and coalescent models.

Genetic samples can be used to understand and predict the behaviour of species living in a fragmented and temporally changing environment. In this regard, models of coalescence conditioned to an environment through an explicit modelling of population growth and migration have been developed in recent years, and simulators implementing these models have been developed, enabling biologists to estimate parameters of interest with Approximate Bayesian Computation techniques. However, model choice remains limited, and developing new coalescence simulators is extremely time consuming because code re-use is limited. We present Quetzal, a C++ library composed of re-usable components, which is sufficiently general to efficiently implement a wide range of spatially explicit coalescence-based environmental models of population genetics and to embed the simulation in an Approximate Bayesian Computation framework. Quetzal is not a simulation program, but a toolbox for programming simulators aimed at the community of scientific coders and research software engineers in molecular ecology and phylogeography. This new code resource is open-source and available at https://becheler.github.io/pages/quetzal.html along with other documentation resources.

Evolution of a polydnavirus gene in relation to parasitoid-host species immune resistance.

CrV1, a polydisperse DNA virus (polydnavirus or PDV) gene contributes to the suppression of host immunity in Cotesia genus parasitoids. Its molecular evolution was analyzed in relation to levels of resistance in the sympatric host species. Natural selection for nonsynonymous substitutions (positive Darwinian selection) was observed at specific amino acid sites among CrV1 variants; particularly, between parasitoid strains immune suppressive and nonimmune suppressive to the main resistant stem borer host, Busseola fusca. In Cotesia sesamiae, geographic distribution of CrV1 alleles in Kenya was correlated to the relative abundance of B. fusca. These results suggest that PDV genes evolve through natural selection and are genetically linked to factors of suppression of local host resistance. We discuss the forces driving the evolution of CrV1 and its use as a marker to understand parasitoid adaptation to host resistance in biological control.

An isometric virus of the potato tuber moth Tecia solanivora (Povolny) (Lepidoptera: Gelechiidae) has a tri-segmented RNA genome.

A small isometric virus has been isolated from larvae of the Guatemala potato tuber moth, Tecia solanivora (Povolny), collected in Ecuador. It was designated the Anchilibi virus (AnchV). The non-enveloped viral particles have an estimated diameter of 32+/-2 nm. Three major proteins were found in virions, with estimated sizes of 102.0+/-2.1, 95.8+/-2.0 and 92.4+/-1.5 kDa for AnchV as determined by polyacrylamide gel electrophoresis. After denaturing agarose gel electrophoresis, the genome of AnchV appeared to be a tri-segmented single-stranded RNA with fragment sizes of 4.1+/-0.2, 2.8+/-0.2 and 1.65+/-0.2 kb. In addition to a high virulence towards its original host, AnchV also caused high mortality in larvae of two other potato tuber moth species, Phthorimaea operculella (Zeller) and Symmetrischema (tangolias) plaesiosema (Turner). Electron microscopy confirmed that AnchV replication occurs in the cell cytoplasm, mainly in vesicles. Several important characteristics exhibited by this virus differ from those reported for known families of insect viruses. Thus, AnchV might be member of a new taxonomic group.

Genome ecosystem and transposable elements species.

Transposable elements are known to be "selfish DNA" sequences able to spread and be maintained in all genomes analyzed so far. Their evolution depends on the interaction they have with the other components of the genome, including genes and other transposable elements. These relationships are complex and have often been compared to those of species living and competing in an ecosystem. The aim of this current work is a proposition to fill the conceptual gap existing between genome biology and ecology, assuming that genomic components, such as transposable elements families, can be compared to species interacting in an ecosystem. Using this framework, some of the main models defined in the population genetics of transposable elements can then been reformulated, and some new kinds of realistic relationships, such as symbiosis between different genomic components, can then be modelled and explored.