Sex-specific selection patterns in a dioecious insect-pollinated plant.

Competition for mate acquisition is the hallmark of any sexual organism. In insect-pollinated plants, competition to attract pollinators is expected to result in pollinator-mediated selection on attractive floral traits. This could overlap with sexual selection if the number of mating partners increases with pollinator attraction, resulting in an improved reproductive success. In this study, we measured a set of floral traits and estimated individual fitness in male and female Silene dioica in an experimental population. Results align with the predictions of Bateman's principles, in the absence of pollen limitation. In females, natural selection acted on traits that are typically linked with fertility (number of flowers and number of gametes), and selection strength was similar in open- and hand-pollinated females, suggesting a limited role of pollinator-mediated selection. In males, flowering duration and corolla width were positively associated with both reproductive success and number of mates, suggesting that sexual selection has played a role in the evolution of these traits. The use of Bateman's metrics further confirmed stronger sexual selection in males than in females. Taken together, our results shed light on the occurrence of sex-specific patterns of selection in an insect-pollinated plant population.

On the function of flower number: disentangling fertility from pollinator-mediated selection.

In animal-pollinated angiosperms, the 'male-function' hypothesis claims that male reproductive success (RS) should benefit from large floral displays, through pollinator attraction, while female RS is expected to be mainly limited by resource availability. As appealing as this theory might be, studies comparing selection strength on flower number in both sexes rarely document the expected asymmetry. This discrepancy could arise because flower number impacts both pollinator attraction and overall gamete number. In this study, we artificially manipulate floral displays to disentangle the fertility versus pollinator attraction components of selection, both in terms of mating and RS. In females, flower number was under strong fertility selection, as predicted in the absence of pollen limitation. By contrast, in males, flower number was mainly under sexual selection, which in turn increased male RS. However, these selection patterns were not different in males with artificially increased floral displays. This suggests that sexual selection acting on flower number in males does not occur because flower number increases pollinator attraction, but rather because more pollen is available to disperse on more mates. Our study illustrates the power of disentangling various components of selection with potentially sex-specific effects for understanding the evolution of sexual dimorphism.

Widespread coexistence of self-compatible and self-incompatible phenotypes in a diallelic self-incompatibility system in Ligustrum vulgare (Oleaceae).

The breakdown of self-incompatibility (SI) in angiosperms is one of the most commonly observed evolutionary transitions. While multiple examples of SI breakdown have been documented in natural populations, there is strikingly little evidence of stable within-population polymorphism with both inbreeding (self-compatible) and outcrossing (self-incompatible) individuals. This absence of breeding system polymorphism corroborates theoretical expectations that predict that in/outbreeding polymorphism is possible only under very restricted conditions. However, theory also predicts that a diallelic sporophytic SI system should facilitate the maintenance of such polymorphism. We tested this prediction by studying the breeding system of Ligustrum vulgare L., an insect-pollinated hermaphroditic species of the Oleaceae family. Using stigma tests with controlled pollination and paternity assignment of open-pollinated progenies, we confirmed the existence of two self-incompatibility groups in this species. We also demonstrated the occurrence of self-compatible individuals in different populations of Western Europe arising from a mutation affecting the functioning of the pollen component of SI. Our results show that the observed low frequency of self-compatible individuals in natural populations is compatible with theoretical predictions only if inbreeding depression is very high.

Diallelic self-incompatibility is the main determinant of fertilization patterns in olive orchards.

Self-incompatibility (SI) in flowering plants potentially represents a major obstacle for sexual reproduction, especially when the number of S-alleles is low. The situation is extreme in the commercially important olive tree, where in vitro pollination assays suggested the existence of a diallelic SI (DSI) system involving only two groups (G1 and G2). Varieties belonging to the same SI group cannot fertilize each other, such that successful fruit production is predicted to require pollination between varieties of different groups. To test this prediction, we explored the extent to which the DSI system determines fertilization patterns under field conditions. One hundred and seventeen olive cultivars were first genotyped using 10 highly polymorphic dinucleotide Simple Sequence Repeat (SSR) markers to ascertain varietal identity. Cultivars were then phenotyped through controlled pollination tests to assign each of them to one of the two SI groups. We then collected and genotyped 1440 open pollinated embryos from five different orchards constituted of seven local cultivars with known group of incompatibility groups. Embryos genotype information were used: (i) to assign embryos to the most likely pollen donor genotype in the neighbourhood using paternity analysis, and (ii) to compare the composition of the pollen cloud genetic among recipient trees in the five sites. The paternity analysis showed that the DSI system is the main determinant of fertilization success under field open pollination conditions: G1 cultivars sired seeds exclusively on G2 cultivars, and reciprocally. No self-fertilization events were observed. Our results demonstrate that DSI is a potent force determining pollination success among varieties within olive orchards used for production. They have the potential to improve management practices by guiding the selection of compatible varieties to avoid planting orchards containing sets of varieties with strongly unbalanced SI groups, as these would lead to suboptimal olive production.

CMS-G from Beta vulgaris ssp. maritima is maintained in natural populations despite containing an atypical cytochrome c oxidase.

While mitochondrial mutants of the respiratory machinery are rare and often lethal, cytoplasmic male sterility (CMS), a mitochondrially inherited trait that results in pollen abortion, is frequently encountered in wild populations. It generates a breeding system called gynodioecy. In Beta vulgaris ssp. maritima, a gynodioecious species, we found CMS-G to be widespread across the distribution range of the species. Despite the sequencing of the mitochondrial genome of CMS-G, the mitochondrial sterilizing factor causing CMS-G is still unknown. By characterizing biochemically CMS-G, we found that the expression of several mitochondrial proteins is altered in CMS-G plants. In particular, Cox1, a core subunit of the cytochrome c oxidase (complex IV), is larger but can still assemble into complex IV. However, the CMS-G-specific complex IV was only detected as a stabilized dimer. We did not observe any alteration of the affinity of complex IV for cytochrome c; however, in CMS-G, complex IV capacity is reduced. Our results show that CMS-G is maintained in many natural populations despite being associated with an atypical complex IV. We suggest that the modified complex IV could incur the associated cost predicted by theoretical models to maintain gynodioecy in wild populations.

Social and spatial effects on genetic variation between foraging flocks in a wild bird population.

Social interactions are rarely random. In some instances, animals exhibit homophily or heterophily, the tendency to interact with similar or dissimilar conspecifics, respectively. Genetic homophily and heterophily influence the evolutionary dynamics of populations, because they potentially affect sexual and social selection. Here, we investigate the link between social interactions and allele frequencies in foraging flocks of great tits (Parus major) over three consecutive years. We constructed co-occurrence networks which explicitly described the splitting and merging of 85,602 flocks through time (fission-fusion dynamics), at 60 feeding sites. Of the 1,711 birds in those flocks, we genotyped 962 individuals at 4,701 autosomal single nucleotide polymorphisms (SNPs). By combining genomewide genotyping with repeated field observations of the same individuals, we were able to investigate links between social structure and allele frequencies at a much finer scale than was previously possible. We explicitly accounted for potential spatial effects underlying genetic structure at the population level. We modelled social structure and spatial configuration of great tit fission-fusion dynamics with eigenvector maps. Variance partitioning revealed that allele frequencies were strongly affected by group fidelity (explaining 27%-45% of variance) as individuals tended to maintain associations with the same conspecifics. These conspecifics were genetically more dissimilar than expected, shown by genomewide heterophily for pure social (i.e., space-independent) grouping preferences. Genomewide homophily was linked to spatial configuration, indicating spatial segregation of genotypes. We did not find evidence for homophily or heterophily for putative socially relevant candidate genes or any other SNP markers. Together, these results demonstrate the importance of distinguishing social and spatial processes in determining population structure.

Replicated analysis of the genetic architecture of quantitative traits in two wild great tit populations.

Currently, there is much debate on the genetic architecture of quantitative traits in wild populations. Is trait variation influenced by many genes of small effect or by a few genes of major effect? Where is additive genetic variation located in the genome? Do the same loci cause similar phenotypic variation in different populations? Great tits (Parus major) have been studied extensively in long-term studies across Europe and consequently are considered an ecological 'model organism'. Recently, genomic resources have been developed for the great tit, including a custom SNP chip and genetic linkage map. In this study, we used a suite of approaches to investigate the genetic architecture of eight quantitative traits in two long-term study populations of great tits--one in the Netherlands and the other in the United Kingdom. Overall, we found little evidence for the presence of genes of large effects in either population. Instead, traits appeared to be influenced by many genes of small effect, with conservative estimates of the number of contributing loci ranging from 31 to 310. Despite concordance between population-specific heritabilities, we found no evidence for the presence of loci having similar effects in both populations. While population-specific genetic architectures are possible, an undetected shared architecture cannot be rejected because of limited power to map loci of small and moderate effects. This study is one of few examples of genetic architecture analysis in replicated wild populations and highlights some of the challenges and limitations researchers will face when attempting similar molecular quantitative genetic studies in free-living populations.

Fine-scale genetic structure in a wild bird population: the role of limited dispersal and environmentally based selection as causal factors.

Individuals are typically not randomly distributed in space; consequently ecological and evolutionary theory depends heavily on understanding the spatial structure of populations. The central challenge of landscape genetics is therefore to link spatial heterogeneity of environments to population genetic structure. Here, we employ multivariate spatial analyses to identify environmentally induced genetic structures in a single breeding population of 1174 great tits Parus major genotyped at 4701 single-nucleotide polymorphism (SNP) loci. Despite the small spatial scale of the study relative to natal dispersal, we found multiple axes of genetic structure. We built distance-based Moran's eigenvector maps to identify axes of pure spatial variation, which we used for spatial correction of regressions between SNPs and various external traits known to be related to fitness components (avian malaria infection risk, local density of conspecifics, oak tree density, and altitude). We found clear evidence of fine-scale genetic structure, with 21, seven, and nine significant SNPs, respectively, associated with infection risk by two species of avian malaria (Plasmodium circumflexum and P. relictum) and local conspecific density. Such fine-scale genetic structure relative to dispersal capabilities suggests ecological and evolutionary mechanisms maintain within-population genetic diversity in this population with the potential to drive microevolutionary change.

Genomic dissection of variation in clutch size and egg mass in a wild great tit (Parus major) population.

Clutch size and egg mass are life history traits that have been extensively studied in wild bird populations, as life history theory predicts a negative trade-off between them, either at the phenotypic or at the genetic level. Here, we analyse the genomic architecture of these heritable traits in a wild great tit (Parus major) population, using three marker-based approaches - chromosome partitioning, quantitative trait locus (QTL) mapping and a genome-wide association study (GWAS). The variance explained by each great tit chromosome scales with predicted chromosome size, no location in the genome contains genome-wide significant QTL, and no individual SNPs are associated with a large proportion of phenotypic variation, all of which may suggest that variation in both traits is due to many loci of small effect, located across the genome. There is no evidence that any regions of the genome contribute significantly to both traits, which combined with a small, nonsignificant negative genetic covariance between the traits, suggests the absence of genetic constraints on the independent evolution of these traits. Our findings support the hypothesis that variation in life history traits in natural populations is likely to be determined by many loci of small effect spread throughout the genome, which are subject to continued input of variation by mutation and migration, although we cannot exclude the possibility of an additional input of major effect genes influencing either trait.

The design and cross-population application of a genome-wide SNP chip for the great tit Parus major.

The vast amount of phenotypic information collected in some wild animal populations makes them extremely valuable for unravelling the genetics of ecologically important traits and understanding how populations adapt to changes in their environment. Next generation sequencing has revolutionized the development of large marker panels in species previously lacking genomic resources. In this study, a unique genomics toolkit was developed for the great tit (Parus major), a model species in ecology and behavioural biology. This toolkit consists of nearly 100,000 SNPs, over 250 million nucleotides of assembled genomic DNA and more than 80 million nucleotides of assembled expressed sequences. A SNP chip with 9193 SNP markers expected to be spaced evenly along the great tit genome was used to genotype 4702 birds from two of the most intensively studied natural vertebrate populations [Wytham Woods/Bagley Woods (United Kingdom) and de Hoge Veluwe/Westerheide (The Netherlands)]. We show that (i) SNPs identified in either of the two populations have a high genotyping success in the other population, (ii) the minor allele frequencies of the SNPs are highly correlated between the two populations and (iii) despite this high correlation, a large number of SNPs display significant differentiation (F(ST) ) between the populations, with an overrepresentation of genes involved in cardiovascular development close to these SNPs. The developed resources provide the basis for unravelling the genetics of important traits in many long-term studies of great tits. More generally, the protocols and pitfalls encountered will be of use for those developing similar resources.

Gynodioecy in structured populations: understanding fine-scale sex ratio variation in Beta vulgaris ssp. maritima.

Natural selection, random processes and gene flow are known to generate sex ratio variations among sexually polymorphic plant populations. In gynodioecious species, in which hermaphrodites and females coexist, the relative effect of these processes on the maintenance of sex polymorphism is still up for debate. The aim of this study was to document sex ratio and cytonuclear genetic variation at a very local scale in wind-pollinated gynodioecious Beta vulgaris ssp. maritima and attempt to elucidate which processes explained the observed variation. The study sites were characterized by geographically distinct patches of individuals and appeared to be dynamic entities, with recurrent establishment of distinct haplotypes through independent founder events. Along with substantial variation in sex ratio and unexpectedly low gene flow within study sites, our results showed a high genetic differentiation among a mosaic of genetically distinct demes, with isolation by distance or abrupt genetic discontinuities taking place within a few tens of metres. Overall, random founder events with restricted gene flow could be primary determinants of sex structure, by promoting the clumping of sex-determining genes. Such high levels of sex structure provide a landscape for differential selection acting on sex-determining genes, which could modify the conditions of maintenance of gynodioecy in structured populations.

Effects of fine-scale genetic structure on male mating success in gynodioecious Beta vulgaris ssp. maritima.

Plant mating systems are known to influence population genetic structure because pollen and seed dispersal are often spatially restricted. However, the reciprocal outcomes of population structure on the dynamics of polymorphic mating systems have received little attention. In gynodioecious sea beet (Beta vulgaris ssp. maritima), three sexual types co-occur: females carrying a cytoplasmic male sterility (CMS) gene, hermaphrodites carrying a non-CMS cytoplasm and restored hermaphrodites that carry CMS genes and nuclear restorer alleles. This study investigated the effects of fine-scale genetic structure on male reproductive success of the two hermaphroditic forms. Our study population was strongly structured and characterized by contrasting local sex-ratios. Pollen flow was constrained over short distances and depended on local plant density. Interestingly, restored hermaphrodites sired significantly more seedlings than non-CMS hermaphrodites, despite the previous observation that the former produce pollen of lower quality than the latter. This result was explained by the higher frequency of females in the local vicinity of restored (CMS) hermaphrodites as compared to non-CMS hermaphrodites. Population structure thus strongly influences individual fitness and may locally counteract the expected effects of selection, suggesting that understanding fine scale population processes is central to predicting the evolution of gender polymorphism in angiosperms.

Nuclear and cytoplasmic genetic diversity in weed beet and sugar beet accessions compared to wild relatives: new insights into the genetic relationships within the Beta vulgaris complex species.

Hybridization between cultivated species and their wild relatives is now widely considered to be common. In the Beta vulgaris complex, the sugar beet seed multiplication areas have been the scene of inadvertent pollination of sugar beet seed bearers by wild ruderal pollen donors, generating a weedy form of beet which infests sugar beet fields in European countries. Up to now, investigations of evolutionary dynamics of genetic diversity within the B. vulgaris complex were addressed using few genetical markers and few accessions. In this study, we tackled this issue using a panel of complementary markers: five nuclear microsatellite loci, four mitochondrial minisatellite loci and one chloroplastic PCR-RFLP marker. We sampled 1,640 individuals that illustrate the actual distribution of inland ruderal beets of South Western France, weed beets and wild sea beets of northern France as well as the diversity of 35 contemporary European diploid cultivars. Nuclear genetic diversity in weed beets appeared to be as high as those of ruderal beets and sea beets, whereas the narrowness of cultivar accessions was confirmed. This genetic bottleneck in cultivars is even more important in the cytoplasmic genome as only one haplotype was found among all sugar beet cultivars. The large majority of weed beet populations also presented this unique cytoplasmic haplotype, as expected owing to their maternal cultivated origin. Nonetheless, various cytoplasmic haplotypes were found within three populations of weed beets, implying wild-to-weed seed flows. Finally, our findings gave new insights into the genetical relationships between the components of the B. vulgaris complex: (1) we found a very strong genetic divergence between wild sea beet and other relatives, which was unexpected given the recent evolutionary history and the full cross-compatibility of all taxa and (2) we definitely confirmed that the classification into cultivated, wild, ruderal and weed forms according to their geographical location, phenotype or their domesticated status is clearly in accordance with genetic clustering despite the very recent domestication process of sugar beet.