Transcriptomic basis of sex loss in the pea aphid.

BACKGROUND: Transitions from sexual to asexual reproduction are common in eukaryotes, but the underlying mechanisms remain poorly known. The pea aphid-Acyrthosiphon pisum-exhibits reproductive polymorphism, with cyclical parthenogenetic and obligate parthenogenetic lineages, offering an opportunity to decipher the genetic basis of sex loss. Previous work on this species identified a single 840 kb region controlling reproductive polymorphism and carrying 32 genes. With the aim of identifying the gene(s) responsible for sex loss and the resulting consequences on the genetic programs controlling sexual or asexual embryogenesis, we compared the transcriptomic response to photoperiod shortening-the main sex-inducing cue-of a sexual and an obligate asexual lineage of the pea aphid, focusing on heads (where the photoperiodic cue is detected) and embryos (the final target of the cue). RESULTS: Our analyses revealed that four genes (one expressed in the head, and three in the embryos) of the region responded differently to photoperiod in the two lineages. We also found that the downstream genetic programs expressed during embryonic development of a future sexual female encompass ∼1600 genes, among which miRNAs, piRNAs and histone modification pathways are overrepresented. These genes mainly co-localize in two genomic regions enriched in transposable elements (TEs). CONCLUSIONS: Our results suggest that the causal polymorphism(s) in the 840 kb region somehow impair downstream epigenetic and post-transcriptional regulations in obligate asexual lineages, thereby sustaining asexual reproduction.

The locus of sexual selection: moving sexual selection studies into the post-genomics era.

Sexual selection drives fundamental evolutionary processes such as trait elaboration and speciation. Despite this importance, there are surprisingly few examples of genes unequivocally responsible for variation in sexually selected phenotypes. This lack of information inhibits our ability to predict phenotypic change due to universal behaviours, such as fighting over mates and mate choice. Here, we discuss reasons for this apparent gap and provide recommendations for how it can be overcome by adopting contemporary genomic methods, exploiting underutilized taxa that may be ideal for detecting the effects of sexual selection and adopting appropriate experimental paradigms. Identifying genes that determine variation in sexually selected traits has the potential to improve theoretical models and reveal whether the genetic changes underlying phenotypic novelty utilize common or unique molecular mechanisms. Such a genomic approach to sexual selection will help answer questions in the evolution of sexually selected phenotypes that were first asked by Darwin and can furthermore serve as a model for the application of genomics in all areas of evolutionary biology.

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 F(ST) 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.

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.

Inferring landscape effects on dispersal from genetic distances: how far can we go?

Functional connectivity affects demography and gene dynamics in fragmented populations. Besides species-specific dispersal ability, the connectivity between local populations is affected by the landscape elements encountered during dispersal. Documenting these effects is thus a central issue for the conservation and management of fragmented populations. In this study, we compare the power and accuracy of three methods (partial correlations, regressions and Approximate Bayesian Computations) that use genetic distances to infer the effect of landscape upon dispersal. We use stochastic individual-based simulations of fragmented populations surrounded by landscape elements that differ in their permeability to dispersal. The power and accuracy of all three methods are good when there is a strong contrast between the permeability of different landscape elements. The power and accuracy can be further improved by restricting analyses to adjacent pairs of populations. Landscape elements that strongly impede dispersal are the easiest to identify. However, power and accuracy decrease drastically when landscape complexity increases and the contrast between the permeability of landscape elements decreases. We provide guidelines for future studies and underline the needs to evaluate or develop approaches that are more powerful.

Predicting the deleterious effects of mutation load in fragmented populations.

Human-induced habitat fragmentation constitutes a major threat to biodiversity. Both genetic and demographic factors combine to drive small and isolated populations into extinction vortices. Nevertheless, the deleterious effects of inbreeding and drift load may depend on population structure, migration patterns, and mating systems and are difficult to predict in the absence of crossing experiments. We performed stochastic individual-based simulations aimed at predicting the effects of deleterious mutations on population fitness (offspring viability and median time to extinction) under a variety of settings (landscape configurations, migration models, and mating systems) on the basis of easy-to-collect demographic and genetic information. Pooling all simulations, a large part (70%) of variance in offspring viability was explained by a combination of genetic structure (F(ST)) and within-deme heterozygosity (H(S)). A similar part of variance in median time to extinction was explained by a combination of local population size (N) and heterozygosity (H(S)). In both cases the predictive power increased above 80% when information on mating systems was available. These results provide robust predictive models to evaluate the viability prospects of fragmented populations.

Habitat-quality effects on metapopulation dynamics in greater white-toothed shrews, Crocidura russula.

The effects of patch size and isolation on metapopulation dynamics have received wide empirical support and theoretical formalization. By contrast, the effects of patch quality seem largely underinvestigated, partly due to technical difficulties in properly assessing quality. Here we combine habitat-quality modeling with four years of demographic monitoring in a metapopulation of greater white-toothed shrews (Crocidura russula) to investigate the role of patch quality on metapopulation processes. Together, local patch quality and connectivity significantly enhanced local population sizes and occupancy rates (R2 = 14% and 19%, respectively). Accounting for the quality of patches connected to the focal one and acting as potential sources improved slightly the model explanatory power for local population sizes, pointing to significant source-sink dynamics. Local habitat quality, in interaction with connectivity, also increased colonization rate (R2 = 28%), suggesting the ability of immigrants to target high-quality patches. Overall, patterns were best explained when assuming a mean dispersal distance of 800 m, a realistic value for the species under study. Our results thus provide evidence that patch quality, in interaction with connectivity, may affect major demographic processes.

Isolation and characterization of microsatellite loci for the European tree frog (Hyla arborea).

We developed 11 new microsatellite markers for the European tree frog (Hyla arborea), and tested patterns of polymorphism in 54 adults (27 males and 27 females) from two ponds close to Lausanne (Western Switzerland). One marker was sex linked and two pairs displayed linkage disequilibrium. Comparisons of allele numbers with heterozygosity values support a stepwise-mutation model at neutral equilibrium, with mutation rates spanning nearly two orders of magnitude. These markers will prove useful for population genetic studies and fine-scale investigations requiring genetic assignment techniques.

Recombination is suppressed and variability reduced in a nascent Y chromosome.

Several hypotheses have been elaborated to account for the evolutionary decay commonly observed in full-fledged Y chromosomes. Enhanced drift, background selection and selective sweeps, which are expected to result from reduced recombination, may all share responsibilities in the initial decay of proto-Y chromosomes, but little empirical information has been gathered so far. Here we take advantage of three markers that amplify on both of the morphologically undifferentiated sex chromosomes of the European tree frog (Hyla arborea) to show that recombination is suppressed in males (the heterogametic sex) but not in females. Accordingly, genetic variability is reduced on the Y, but in a way that can be accounted for by merely the number of chromosome copies per breeding pair, without the need to invoke background selection or selective sweeps.

A sex-specific marker reveals male heterogamety in European tree frogs.

Most amphibians examined so far show undifferentiated sex chromosomes. The heterogametic sex's identity, usually revealed through indirect means, often varies among closely related species or even populations (as do sex-linkage groups), suggesting great evolutionary instability of the sex-determining genes. Here we take advantage of a sex-specific marker that amplifies in several related species of European tree frogs (Hyla arborea group) to disclose a homogeneous pattern of male heterogamety. Besides relevance for evolutionary studies of sex determination in amphibians, our results have potential for addressing practical issues in conservation biology because sex reversal by anthropogenic endocrine disruptors is considered one possible cause of amphibian decline.

Multilevel genetic analyses of two European supercolonies of the Argentine ant, Linepithema humile.

Some ants have an extraordinary unicolonial social organization, whereby individuals mix freely among physically separated nests. Recently, it was shown that the European population of Linepithema humile consisted of two enormous unicolonial supercolonies. Workers of the same supercolony are never aggressive to each other. In contrast, aggressiveness is invariably high between workers from different supercolonies. Here we investigated whether gene flow occurs between two supercolonies. We identified a contact zone in which we sampled 46 nests. For each nest, aggression tests were conducted against workers from reference nests from both supercolonies. Workers were always very aggressive towards workers of one of the supercolonies but not to workers of the other. Thus, all nests could be clearly assigned to one of the two supercolonies. For 22 of the 46 nests, we genotyped 15-16 workers at five microsatellite loci. A four-level hierarchical analysis of variance revealed very strong genetic differentiation between the two supercolonies (F(SUPERCOLONY-TOTAL) = 0.541) and low differentiation between sectors (i.e. group of nests connected together with trails) within supercolonies (F(SECTOR-SUPERCOLONY) = 0.064). The very high differentiation between the two supercolonies indicates a lack of ongoing gene flow, a conclusion further bolstered by the finding that the two supercolonies share no common alleles at two of the five microsatellite loci. A Bayesian clustering method also revealed the occurrence of two distinct clusters. These clusters exactly match the grouping obtained by aggression tests. None of the 332 genotyped individuals were admixed despite the fact that some nests of the two supercolonies were separated by less than 30 m. These results demonstrate that the two supercolonies have completely separate gene pools.