Restriction site-associated DNA sequencing, genotyping error estimation and de novo assembly optimization for population genetic inference.

Restriction site-associated DNA sequencing (RADseq) provides researchers with the ability to record genetic polymorphism across thousands of loci for nonmodel organisms, potentially revolutionizing the field of molecular ecology. However, as with other genotyping methods, RADseq is prone to a number of sources of error that may have consequential effects for population genetic inferences, and these have received only limited attention in terms of the estimation and reporting of genotyping error rates. Here we use individual sample replicates, under the expectation of identical genotypes, to quantify genotyping error in the absence of a reference genome. We then use sample replicates to (i) optimize de novo assembly parameters within the program Stacks, by minimizing error and maximizing the retrieval of informative loci; and (ii) quantify error rates for loci, alleles and single-nucleotide polymorphisms. As an empirical example, we use a double-digest RAD data set of a nonmodel plant species, Berberis alpina, collected from high-altitude mountains in Mexico.

Changes in reproductive roles are associated with changes in gene expression in fire ant queens.

Abstract In species with social hierarchies, the death of dominant individuals typically upheaves the social hierarchy and provides an opportunity for subordinate individuals to become reproductives. Such a phenomenon occurs in the monogyne form of the fire ant, Solenopsis invicta, where colonies typically contain a single wingless reproductive queen, thousands of workers and hundreds of winged nonreproductive virgin queens. Upon the death of the mother queen, many virgin queens shed their wings and initiate reproductive development instead of departing on a mating flight. Workers progressively execute almost all of them over the following weeks. To identify the molecular changes that occur in virgin queens as they perceive the loss of their mother queen and begin to compete for reproductive dominance, we collected virgin queens before the loss of their mother queen, 6 h after orphaning and 24 h after orphaning. Their RNA was extracted and hybridized against microarrays to examine the expression levels of approximately 10 000 genes. We identified 297 genes that were consistently differentially expressed after orphaning. These include genes that are putatively involved in the signalling and onset of reproductive development, as well as genes underlying major physiological changes in the young queens.

Mediterranean populations of the lesser white-toothed shrew (Crocidura suaveolens group): an unexpected puzzle of Pleistocene survivors and prehistoric introductions.

An earlier study revealed the strong phylogeographical structure of the lesser white-toothed shrew (Crocidura suaveolens group) within the northern Palaearctic. Here, we aim to reconstruct the colonization history of Mediterranean islands and to clarify the biogeography and phylogeographical relationships of the poorly documented Middle East region with the northern Palaearctic. We performed analyses on 998-bp-long haplotypes of the mitochondrial cytochrome b gene of 143 samples collected around the Mediterranean basin, including islands and the Middle East. The analyses suggest that the Cypriot shrew belongs to the rare group of relict insular Pleistocene mammal taxa that have survived to the present day. In contrast, the Cretan, Corsican and Menorcan populations were independently introduced from the Middle East during the Holocene. The phylogeographical structure of this temperate Palaearctic species within the Middle East appears to be complex and rich in diversity, probably reflecting fragmentation of the area by numerous mountain chains. Four deeply divergent clades of the C. suaveolens group occur in the area, meaning that a hypothetical contact zone remains to be located in central western Iran.

Social chromosome variants differentially affect queen determination and the survival of workers in the fire ant Solenopsis invicta.

Intraspecific variation in social organization is common, yet the underlying causes are rarely known. An exception is the fire ant Solenopsis invicta in which the existence of two distinct forms of social colony organization is under the control of the two variants of a pair of social chromosomes, SB and Sb. Colonies containing exclusively SB/SB workers accept only one single queen and she must be SB/SB. By contrast, when colonies contain more than 10% of SB/Sb workers, they accept several queens but only SB/Sb queens. The variants of the social chromosome are associated with several additional important phenotypic differences, including the size, fecundity and dispersal strategies of queens, aggressiveness of workers, and sperm count in males. However, little is known about whether social chromosome variants affect fitness in other life stages. Here, we perform experiments to determine whether differential selection occurs during development and in adult workers. We find evidence that the Sb variant of the social chromosome increases the likelihood of female brood to develop into queens and that adult SB/Sb workers, the workers that cull SB/SB queens, are overrepresented in comparison to SB/SB workers. This demonstrates that supergenes such as the social chromosome can have complex effects on phenotypes at various stages of development.

A simple genetic basis for complex social behaviour mediates widespread gene expression differences.

A remarkable social polymorphism is controlled by a single Mendelian factor in the fire ant Solenopsis invicta. A genomic element marked by the gene Gp-9 determines whether workers tolerate one or many fertile queens in their colony. Gp-9 was recently shown to be part of a supergene with two nonrecombining variants, SB and Sb. SB/SB and SB/Sb queens differ in how they initiate new colonies, and in many physiological traits, for example odour and maturation rate. To understand how a single genetic element can affect all these traits, we used a microarray to compare gene expression patterns between SB/SB and SB/Sb queens of three different age classes: 1-day-old unmated queens, 11-day-old unmated queens and mated, fully reproductive queens collected from mature field colonies. The number of genes that were differentially expressed between SB/SB and SB/Sb queens of the same age class was smallest in 1-day-old queens, maximal in 11-day-old queens and intermediate in reproductive queens. Gene ontology analysis showed that SB/SB queens upregulate reproductive genes faster than SB/Sb queens. For all age classes, genes inside the supergene were overrepresented among the differentially expressed genes. Consistent with the hypothesized greater number of transposons in the Sb supergene, 13 transposon genes were upregulated in SB/Sb queens. Viral genes were also upregulated in SB/Sb mature queens, consistent with the known greater parasite load in colonies headed by SB/Sb queens compared with colonies headed by SB/SB queens. Eighteen differentially expressed genes between reproductive queens were involved in chemical signalling. Our results suggest that many genes in the supergene are involved in regulating social organization and queen phenotypes in fire ants.

Discordances between phylogenetic and morphological patterns in alpine leaf beetles attest to an intricate biogeographic history of lineages in postglacial Europe.

Pleistocene glacial and interglacial periods have moulded the evolutionary history of European cold-adapted organisms. The role of the different mountain massifs has, however, not been accurately investigated in the case of high-altitude insect species. Here, we focus on three closely related species of non-flying leaf beetles of the genus Oreina (Coleoptera, Chrysomelidae), which are often found in sympatry within the mountain ranges of Europe. After showing that the species concept as currently applied does not match barcoding results, we show, based on more than 700 sequences from one nuclear and three mitochondrial genes, the role of biogeography in shaping the phylogenetic hypothesis. Dating the phylogeny using an insect molecular clock, we show that the earliest lineages diverged more than 1 Mya and that the main shift in diversification rate occurred between 0.36 and 0.18 Mya. By using a probabilistic approach on the parsimony-based dispersal/vicariance framework (MP-DIVA) as well as a direct likelihood method of state change optimization, we show that the Alps acted as a cross-roads with multiple events of dispersal to and reinvasion from neighbouring mountains. However, the relative importance of vicariance vs. dispersal events on the process of rapid diversification remains difficult to evaluate because of a bias towards overestimation of vicariance in the DIVA algorithm. Parallels are drawn with recent studies of cold-adapted species, although our study reveals novel patterns in diversity and genetic links between European mountains, and highlights the importance of neglected regions, such as the Jura and the Balkanic range.

Genetic diversity and differentiation processes in the ploidy series of Olea europaea L.: a multiscale approach from subspecies to insular populations.

Geographical isolation and polyploidization are central concepts in plant evolution. The hierarchical organization of archipelagos in this study provides a framework for testing the evolutionary consequences for polyploid taxa and populations occurring in isolation. Using amplified fragment length polymorphism and simple sequence repeat markers, we determined the genetic diversity and differentiation patterns at three levels of geographical isolation in Olea europaea: mainland-archipelagos, islands within an archipelago, and populations within an island. At the subspecies scale, the hexaploid ssp. maroccana (southwest Morocco) exhibited higher genetic diversity than the insular counterparts. In contrast, the tetraploid ssp. cerasiformis (Madeira) displayed values similar to those obtained for the diploid ssp. guanchica (Canary Islands). Geographical isolation was associated with a high genetic differentiation at this scale. In the Canarian archipelago, the stepping-stone model of differentiation suggested in a previous study was partially supported. Within the western lineage, an east-to-west differentiation pattern was confirmed. Conversely, the easternmost populations were more related to the mainland ssp. europaea than to the western guanchica lineage. Genetic diversity across the Canarian archipelago was significantly correlated with the date of the last volcanic activity in the area/island where each population occurs. At the island scale, this pattern was not confirmed in older islands (Tenerife and Madeira), where populations were genetically homogeneous. In contrast, founder effects resulted in low genetic diversity and marked genetic differentiation among populations of the youngest island, La Palma.

Evolution of the mating system in colonizing plants.

Colonization is likely to be more successful for species with an ability to self-fertilize and thus to establish new populations as single individuals. As a result, self-compatibility should be common among colonizing species. This idea, labelled 'Baker's law', has been influential in discussions of sexual-system and mating-system evolution. However, its generality has been questioned, because models of the evolution of dispersal and the mating system predict an association between high dispersal rates and outcrossing rather than selfing, and because of many apparent counter examples to the law. The contrasting predictions made by models invoking Baker's law versus those for the evolution of the mating system and dispersal urges a reassessment of how we should view both these traits. Here, I review the literature on the evolution of mating and dispersal in colonizing species, with a focus on conceptual issues. I argue for the importance of distinguishing between the selfing or outcrossing rate and a simple ability to self-fertilize, as well as for the need for a more nuanced consideration of dispersal. Colonizing species will be characterized by different phases in their life pattern: dispersal to new habitat, implying an ecological sieve on dispersal traits; establishment and a phase of growth following colonization, implying a sieve on reproductive traits; and a phase of demographic stasis at high density, during which new trait associations can evolve through local adaptation. This dynamic means that the sorting of mating-system and dispersal traits should change over time, making simple predictions difficult.

Molecular estimation of dispersal for ecology and population genetics

The dispersal process, by which individuals or other dispersing agents such as gametes or seeds move from birthplace to a new settlement locality, has important consequences for the dynamics of genes, individuals, and species. Many of the questions addressed by ecology and evolutionary biology require a good understanding of species' dispersal patterns. Much effort has thus been devoted to overcoming the difficulties associated with dispersal measurement. In this context, genetic tools have long been the focus of intensive research, providing a great variety of potential solutions to measuring dispersal. This methodological diversity is reviewed here to help (molecular) ecologists find their way toward dispersal inference and interpretation and to stimulate further developments.