Lucie Mieuzet

Host range expansion of an introduced insect pest through multiple colonizations of specialized clones

Asexuality confers demographic advantages to invasive taxa, but generally limits adaptive potential for colonizing of new habitats. Therefore, pre‐existing adaptations and habitat tolerance are essential in the success of asexual invaders. We investigated these key factors of invasiveness by assessing reproductive modes and host‐plant adaptations in the pea aphid, Acyrthosiphon pisum, a pest recently introduced into Chile. The pea aphid encompasses lineages differing in their reproductive mode, ranging from obligatory cyclical parthenogenesis to fully asexual reproduction. This species also shows variation in host use, with distinct biotypes specialized on different species of legumes as well as more polyphagous populations. In central Chile, microsatellite genotyping of pea aphids sampled on five crops and wild legumes revealed three main clonal genotypes, which showed striking associations with particular host plants rather than sampling locations. Phenotypic analyses confirmed their strong host specialization and demonstrated parthenogenesis as their sole reproductive mode. The genetic relatedness of these clonal genotypes with corresponding host‐specialized populations from the Old World indicated that each clone descended from a particular Eurasian biotype, which involved at least three successful introduction events followed by spread on different crops. This study illustrates that multiple introductions of highly specialized clones, rather than local evolution in resource use and/or selection of generalist genotypes, can explain the demographic success of a strictly asexual invader.

Permanent Genetic Resources added to Molecular Ecology Resources Database 1 April 2010-31 May 2010

This article documents the addition of 396 microsatellite marker loci to the Molecular Ecology Resources Database. Loci were developed for the following species: Anthocidaris crassispina, Aphis glycines, Argyrosomus regius, Astrocaryum sciophilum, Dasypus novemcinctus, Delomys sublineatus, Dermatemys mawii, Fundulus heteroclitus, Homalaspis plana, Jumellea rossii, Khaya senegalensis, Mugil cephalus, Neoceratitis cyanescens, Phalacrocorax aristotelis, Phytophthora infestans, Piper cordulatum, Pterocarpus indicus, Rana dalmatina, Rosa pulverulenta, Saxifraga oppositifolia, Scomber colias, Semecarpus kathalekanensis, Stichopus monotuberculatus, Striga hermonthica, Tarentola boettgeri and Thermophis baileyi. These loci were cross‐tested on the following species: Aphis gossypii, Sooretamys angouya, Euryoryzomys russatus, Fundulus notatus, Fundulus olivaceus, Fundulus catenatus, Fundulus majalis, Jumellea fragrans, Jumellea triquetra Jumellea recta, Jumellea stenophylla, Liza richardsonii, Piper marginatum, Piper aequale, Piper darienensis, Piper dilatatum, Rana temporaria, Rana iberica, Rana pyrenaica, Semecarpus anacardium, Semecarpus auriculata, Semecarpus travancorica, Spondias acuminata, Holigarna grahamii, Holigarna beddomii, Mangifera indica, Anacardium occidentale, Tarentola delalandii, Tarentola caboverdianus and Thermophis zhaoermii.

Genome scans reveal candidate regions involved in the adaptation to host plant in the pea aphid complex

A major goal in evolutionary biology is to uncover the genetic basis of adaptation. Divergent selection exerted on ecological traits may result in adaptive population differentiation and reproductive isolation and affect differentially the level of genetic divergence along the genome. Genome‐wide scan of large sets of individuals from multiple populations is a powerful approach to identify loci or genomic regions under ecologically divergent selection. Here, we focused on the pea aphid, a species complex of divergent host races, to explore the organization of the genomic divergence associated with host plant adaptation and ecological speciation. We analysed 390 microsatellite markers located at variable distances from predicted genes in replicate samples of sympatric populations of the pea aphid collected on alfalfa, red clover and pea, which correspond to three common host‐adapted races reported in this species complex. Using a method that accounts for the hierarchical structure of our data set, we found a set of 11 outlier loci that show higher genetic differentiation between host races than expected under the null hypothesis of neutral evolution. Two of the outliers are close to olfactory receptor genes and three other nearby genes encoding salivary proteins. The remaining outliers are located in regions with genes of unknown functions, or which functions are unlikely to be involved in interactions with the host plant. This study reveals genetic signatures of divergent selection across the genome and provides an inventory of candidate genes responsible for plant specialization in the pea aphid, thereby setting the stage for future functional studies.

Bacterial communities associated with host-adapted populations of pea aphids revealed by deep sequencing of 16S ribosomal DNA.

Associations between microbes and animals are ubiquitous and hosts may benefit from harbouring microbial communities through improved resource exploitation or resistance to environmental stress. The pea aphid, Acyrthosiphon pisum, is the host of heritable bacterial symbionts, including the obligate endosymbiont Buchnera aphidicola and several facultative symbionts. While obligate symbionts supply aphids with key nutrients, facultative symbionts influence their hosts in many ways such as protection against natural enemies, heat tolerance, color change and reproduction alteration. The pea aphid also encompasses multiple plant-specialized biotypes, each adapted to one or a few legume species. Facultative symbiont communities differ strongly between biotypes, although bacterial involvement in plant specialization is uncertain. Here, we analyse the diversity of bacterial communities associated with nine biotypes of the pea aphid complex using amplicon pyrosequencing of 16S rRNA genes. Combined clustering and phylogenetic analyses of 16S sequences allowed identifying 21 bacterial OTUs (Operational Taxonomic Unit). More than 98% of the sequencing reads were assigned to known pea aphid symbionts. The presence of Wolbachia was confirmed in A. pisum while Erwinia and Pantoea, two gut associates, were detected in multiple samples. The diversity of bacterial communities harboured by pea aphid biotypes was very low, ranging from 3 to 11 OTUs across samples. Bacterial communities differed more between than within biotypes but this difference did not correlate with the genetic divergence between biotypes. Altogether, these results confirm that the aphid microbiota is dominated by a few heritable symbionts and that plant specialization is an important structuring factor of bacterial communities associated with the pea aphid complex. However, since we examined the microbiota of aphid samples kept a few generations in controlled conditions, it may be that bacterial diversity was underestimated due to the possible loss of environmental or transient taxa.

Masculinization of the X Chromosome in the Pea Aphid

Evolutionary theory predicts that sexually antagonistic mutations accumulate differentially on the X chromosome and autosomes in species with an XY sex-determination system, with effects (masculinization or feminization of the X) depending on the dominance of mutations. Organisms with alternative modes of inheritance of sex chromosomes offer interesting opportunities for studying sexual conflicts and their resolution, because expectations for the preferred genomic location of sexually antagonistic alleles may differ from standard systems. Aphids display an XX/X0 system and combine an unusual inheritance of the X chromosome with the alternation of sexual and asexual reproduction. In this study, we first investigated theoretically the accumulation of sexually antagonistic mutations on the aphid X chromosome. Our results show that i) the X is always more favourable to the spread of male-beneficial alleles than autosomes, and should thus be enriched in sexually antagonistic alleles beneficial for males, ii) sexually antagonistic mutations beneficial for asexual females accumulate preferentially on autosomes, iii) in contrast to predictions for standard systems, these qualitative results are not affected by the dominance of mutations. Under the assumption that sex-biased gene expression evolves to solve conflicts raised by the spread of sexually antagonistic alleles, one expects that male-biased genes should be enriched on the X while asexual female-biased genes should be enriched on autosomes. Using gene expression data (RNA-Seq) in males, sexual females and asexual females of the pea aphid, we confirm these theoretical predictions. Although other mechanisms than the resolution of sexual antagonism may lead to sex-biased gene expression, we argue that they could hardly explain the observed difference between X and autosomes. On top of reporting a strong masculinization of the aphid X chromosome, our study highlights the relevance of organisms displaying an alternative mode of sex chromosome inheritance to understanding the forces shaping chromosome evolution.

Accelerated Evolution of Sex Chromosomes in Aphids, an X0 System

Sex chromosomes play a role in many important biological processes, including sex determination, genomic conflicts, imprinting, and speciation. In particular, they exhibit several unusual properties such as inheritance pattern, hemizygosity, and reduced recombination, which influence their response to evolutionary factors (e.g., drift, selection, and demography). Here, we examine the evolutionary forces driving X chromosome evolution in aphids, an XO system where females are homozygous (XX) and males are hemizygous (X0) at sex chromosomes. We show by simulations that the unusual mode of transmission of the X chromosome in aphids, coupled with cyclical parthenogenesis, results in similar effective population sizes and predicted levels of genetic diversity for X chromosomes and autosomes under neutral evolution. These results contrast with expectations from standard XX/XY or XX/X0 systems (where the effective population size of the X is three-fourths that of autosomes) and have deep consequences for aphid X chromosome evolution. We then localized 52 microsatellite markers on the X and 351 on autosomes. We genotyped 167 individuals with 356 of these loci and found similar levels of allelic richness on the X and on the autosomes, as predicted by our simulations. In contrast, we detected higher dN and dN/dS ratio for X-linked genes compared with autosomal genes, a pattern compatible with either positive or relaxed selection. Given that both types of chromosomes have similar effective population sizes and that the single copy of the X chromosome of male aphids exposes its recessive genes to selection, some degree of positive selection seems to best explain the higher rates of evolution of X-linked genes. Overall, this study highlights the particular relevance of aphids to study the evolutionary factors driving sex chromosomes and genome evolution.

Genetic Control of Contagious Asexuality in the Pea Aphid

Although evolutionary transitions from sexual to asexual reproduction are frequent in eukaryotes, the genetic bases of such shifts toward asexuality remain largely unknown. We addressed this issue in an aphid species where both sexual and obligate asexual lineages coexist in natural populations. These sexual and asexual lineages may occasionally interbreed because some asexual lineages maintain a residual production of males potentially able to mate with the females produced by sexual lineages. Hence, this species is an ideal model to study the genetic basis of the loss of sexual reproduction with quantitative genetic and population genomic approaches. Our analysis of the co-segregation of ∼300 molecular markers and reproductive phenotype in experimental crosses pinpointed an X-linked region controlling obligate asexuality, this state of character being recessive. A population genetic analysis (>400-marker genome scan) on wild sexual and asexual genotypes from geographically distant populations under divergent selection for reproductive strategies detected a strong signature of divergent selection in the genomic region identified by the experimental crosses. These population genetic data confirm the implication of the candidate region in the control of reproductive mode in wild populations originating from 700 km apart. Patterns of genetic differentiation along chromosomes suggest bidirectional gene flow between populations with distinct reproductive modes, supporting contagious asexuality as a prevailing route to permanent parthenogenesis in pea aphids. This genetic system provides new insights into the mechanisms of coexistence of sexual and asexual aphid lineages.

Genomic regions repeatedly involved in divergence among plant‐specialized pea aphid biotypes

Understanding the genetic bases of biological diversification is a long‐standing goal in evolutionary biology. Here, we investigate whether replicated cases of adaptive divergence involve the same genomic regions in the pea aphid, Acyrthosiphon pisum, a large complex of genetically differentiated biotypes, each specialized on different species of legumes. A previous study identified genomic regions putatively involved in host‐plant adaptation and/or reproductive isolation by performing a hierarchical genome scan in three biotypes. This led to the identification of 11 FST outliers among 390 polymorphic microsatellite markers. In this study, the outlier status of these 11 loci was assessed in eight biotypes specialized on other host plants. Four of the 11 previously identified outliers showed greater genetic differentiation among these additional biotypes than expected under the null hypothesis of neutral evolution (α < 0.01). Whether these hotspots of genomic divergence result from adaptive events, intrinsic barriers or reduced recombination is discussed.

Disentangling the Causes for Faster-X Evolution in Aphids

Abstract The faster evolution of X chromosomes has been documented in several species, and results from the increased efficiency of selection on recessive alleles in hemizygous males and/or from increased drift due to the smaller effective population size of X chromosomes. Aphids are excellent models for evaluating the importance of selection in faster-X evolution because their peculiar life cycle and unusual inheritance of sex chromosomes should generally lead to equivalent effective population sizes for X and autosomes. Because we lack a high-density genetic map for the pea aphid, whose complete genome has been sequenced, we first assigned its entire genome to the X or autosomes based on ratios of sequencing depth in males (X0) to females (XX). Then, we computed nonsynonymous to synonymous substitutions ratios (dN/dS) for the pea aphid gene set and found faster evolution of X-linked genes. Our analyses of substitution rates, together with polymorphism and expression data, showed that relaxed selection is likely to be the greatest contributor to faster-X because a large fraction of X-linked genes are expressed at low rates and thus escape selection. Yet, a minor role for positive selection is also suggested by the difference between substitution rates for X and autosomes for male-biased genes (but not for asexual female-biased genes) and by lower Tajima’s D for X-linked compared with autosomal genes with highly male-biased expression patterns. This study highlights the relevance of organisms displaying alternative chromosomal inheritance to the understanding of forces shaping genome evolution.

Identifying genomic hotspots of differentiation and candidate genes involved in the adaptive divergence of pea aphid host races

Identifying the genomic bases of adaptation to novel environments is a long‐term objective in evolutionary biology. Because genetic differentiation is expected to increase between locally adapted populations at the genes targeted by selection, scanning the genome for elevated levels of differentiation is a first step towards deciphering the genomic architecture underlying adaptive divergence. The pea aphid Acyrthosiphon pisum is a model of choice to address this question, as it forms a large complex of plant‐specialized races and cryptic species, resulting from recent adaptive radiation. Here, we characterized genomewide polymorphisms in three pea aphid races specialized on alfalfa, clover and pea crops, respectively, which we sequenced in pools (poolseq). Using a model‐based approach that explicitly accounts for selection, we identified 392 genomic hotspots of differentiation spanning 47.3 Mb and 2,484 genes (respectively, 9.12% of the genome size and 8.10% of its genes). Most of these highly differentiated regions were located on the autosomes, and overall differentiation was weaker on the X chromosome. Within these hotspots, high levels of absolute divergence between races suggest that these regions experienced less gene flow than the rest of the genome, most likely by contributing to reproductive isolation. Moreover, population‐specific analyses showed evidence of selection in every host race, depending on the hotspot considered. These hotspots were significantly enriched for candidate gene categories that control host–plant selection and use. These genes encode 48 salivary proteins, 14 gustatory receptors, 10 odorant receptors, five P450 cytochromes and one chemosensory protein, which represent promising candidates for the genetic basis of host–plant specialization and ecological isolation in the pea aphid complex. Altogether, our findings open new research directions towards functional studies, for validating the role of these genes on adaptive phenotypes.