Search

A transition to selfing can be beneficial when mating partners are scarce, for example, due to ploidy changes or at species range edges. Here we explain how self-compatibility evolved in diploid SiberianArabidopsis lyrata,and how it contributed to the establishment of allotetraploidA. kamchatica. First, we provide chromosome-level genome assemblies for two self-fertilizing diploidA. lyrataaccessions, one from North America and one from Siberia, including a fully assembled S-locus for the latter. We then propose a sequence of events leading to the loss of self-incompatibility in SiberianA. lyrata,date this independent transition to ∼90 Kya, and infer evolutionary relationships between Siberian and North AmericanA. lyrata,showing an independent transition to selfing in Siberia. Finally, we provide evidence that this selfing SiberianA. lyratalineage contributed to the formation of the allotetraploidA. kamchaticaand propose that the selfing of the latter is mediated by the loss-of-function mutation in a dominantS-allele inherited fromA. lyrata.

Balancing selection is a form of natural selection maintaining diversity at the sites it targets and at linked nucleotide sites. Due to selection favouring heterozygosity, it has the potential to facilitate the accumulation of a “sheltered” load of tightly linked recessive deleterious mutations. However, precisely evaluating the extent of these effects has remained challenging. Taking advantage of plant self-incompatibility as one of the best-understood examples of long-term balancing selection, we provide a highly resolved picture of the genomic extent of balancing selection on the sheltered genetic load. We used targeted genome resequencing to reveal polymorphism of the genomic region flanking the self-incompatibility locus in three sample sets in each of the two closely related plant speciesArabidopsis halleriandA. lyrata, and used 100 control regions from throughout the genome to factor out differences in demographic histories and/or sample structure. Nucleotide polymorphism increased strongly around theS-locus in all sample sets, but only over a limited genomic region, as it became indistinguishable from the genomic background beyond the first 25-30kb. Genes in this chromosomal interval exhibited no excess of mutations at 0-fold degenerated sites relative to putatively neutral sites, hence revealing no detectable weakening of the efficacy of purifying selection even for these most tightly linked genes. Overall, our results are consistent with the predictions of a narrow genomic influence of linkage to theS-locus, and clarify how natural selection in one genomic region affects the evolution of the adjacent genomic regions.

Despite the increasing accessibility of high-throughput sequencing, obtaining high-quality genomic data on non-model organisms without proximate well-assembled and annotated genomes remains challenging. Here we describe a workflow that takes advantage of distant genomic resources and ingroup transcriptomes to select and jointly enrich long open reading frames (ORFs) and ultraconserved elements (UCEs) from genomic samples for integrative studies of microevolutionary and macroevolutionary dynamics. This workflow is applied to samples of the African unionid bivalve tribe Coelaturini (Parreysiinae) at basin and continent-wide scales. Our results indicate that ORFs are efficiently captured without prior identification of intron-exon boundaries. The enrichment of UCEs was less successful, but nevertheless produced a substantial dataset. Exploratory continent-wide phylogenetic analyses with ORF supercontigs (>515,000 parsimony informative sites) resulted in a fully resolved phylogeny, the backbone of which was also retrieved with UCEs (>11,000 informative sites), although some branches lack support in the latter case. Variant calling on the exome of Coelaturini from the Malawi Basin produced ~2,000 SNPs per population pair. Nucleotide diversity and population differentiation was low compared to previous estimates in mollusks, but comparable to those in recently diversifying Malawi cichlids and other taxa at an early stage of speciation. Skimming non-specific sequence data obtained for Coelaturini of the Malawi Basin, we reconstructed the maternally-inherited mitogenome, which displays an identical gene order to that of the most recent common ancestor of Unionidae. Overall, our workflow and results provide exciting perspectives for the development of integrative genomic studies on micro- and macroevolutionary dynamics in non-model organisms.

International audience; Whole-genome sequencing of non-model organisms is now widely accessible and has allowed a range of questions in the field of molecular ecology to be investigated with greater power. However, some genomic regions that are of high biological interest remain problematic for assembly and data-handling. Three such regions are the major histocompatibility complex (MHC), sex-determining regions (SDRs) and the plant self-incompatibility locus (S-locus). Using these as examples, we illustrate the challenges of both assembling and resequencing these highly polymorphic regions and how bioinformatic and technological developments are enabling new approaches to their study. Mapping short-read sequences against multiple alternative references improves genotyping comprehensiveness at the S-locus thereby contributing to more accurate assessments of allelic frequencies. Long-read sequencing, producing reads of several tens to hundreds of kilobase pairs in length, facilitates the assembly of such regions as single sequences can span the multiple duplicated gene copies of the MHC region, and sequence through repetitive stretches and translocations in SDRs and S-locus haplotypes. These advances are adding value to short-read genome resequencing approaches by allowing, for example, more accurate haplotype phasing across longer regions. Finally, we assessed further technical improvements, such as nanopore adaptive sequencing and bioinformatic tools using pangenomes, which have the potential to further expand our knowledge of a number of genomic regions that remain challenging to study with classical resequencing approaches.

We present DILS, a deployable statistical analysis platform for conducting demographic inferences with linked selection from population genomic data using an Approximate Bayesian Computation framework. DILS takes as input single-population or two-population datasets (multilocus fasta sequences) and performs three types of analyses in a hierarchical manner, identifying: 1) the best demographic model to study the importance of gene flow and population size change on the genetic patterns of polymorphism and divergence, 2) the best genomic model to determine whether the effective sizeNeand migration rateN.mare heterogeneously distributed along the genome (implying linked selection) and 3) loci in genomic regions most associated with barriers to gene flow. Also availableviaa web interface, an objective of DILS is to facilitate collaborative research in speciation genomics. Here, we show the performance and limitations of DILS by using simulations, and finally apply the method to published data on a divergence continuum composed by 28 pairs ofMytilusmussel populations/species.

International audience; The Leavenworthia self-incompatibility locus (S locus) consists of paralogs (Lal2, SCRL) of the canonical Brassicaceae S locus genes (SRK, SCR), and is situated in a genomic position that differs from the ancestral one in the Brassicaceae. Unexpectedly, in a small number of Leavenworthia alabamica plants examined, sequences closely resembling exon 1 of SRK have been found, but the function of these has remained unclear. BAC cloning and expression analyses were employed to characterize these SRK-like sequences. An SRK-positive Bacterial Artificial Chromosome clone was found to contain complete SRK and SCR sequences located close by one another in the derived genomic position of the Leavenworthia S locus, and in place of the more typical Lal2 and SCRL sequences. These sequences are expressed in stigmas and anthers, respectively, and crossing data show that the SRK/SCR haplotype is functional in self-incompatibility. Population surveys indicate that < 5% of Leavenworthia S loci possess such alleles. An ancestral translocation or recombination event involving SRK/SCR and Lal2/SCRL likely occurred, together with neofunctionalization of Lal2/SCRL, and both haplotype groups now function as Leavenworthia S locus alleles. These findings suggest that S locus alleles can have distinctly different evolutionary origins.

UMR AGAP - équipe AFEF - Architecture et fonctionnement des espèces fruitières; International audience; The olive (Olea europaea L.) is a typical important perennial crop species for which the genetic determination and even functionality of self-incompatibility (SI) are still largely unresolved. It is still not known whether SI is under gametophytic or sporophytic genetic control, yet fruit production in orchards depends critically on successful ovule fertilization. We studied the genetic determination of SI in olive in light of recent discoveries in other genera of the Oleaceae family. Using intra- and inter-specific stigma tests on 89 genotypes representative of species-wide olive diversity and the compatibility/incompatibility reactions of progeny plants from controlled crosses, we confirmed that O. europaea shares the same homomorphic diallelic self-incompatibility (DSI) system as the one recently identified in Phillyrea angustifolia and Fraxinus ornus. SI is sporophytic in olive. The incompatibility response differs between the two SI groups in terms of how far pollen tubes grow before growth is arrested within stigma tissues. As a consequence of this DSI system, the chance of cross-incompatibility between pairs of varieties in an orchard is high (50%) and fruit production may be limited by the availability of compatible pollen. The discovery of the DSI system in O. europaea will undoubtedly offer opportunities to optimize fruit production.

Whole genome duplication events are notably widespread in plants and this poses particular challenges for phylogenetic inference in allopolyploid lineages, i.e. lineages that result from the merging of two or more diverged genomes after interspecific hybridization. The nuclear genomes resulting from allopolyploidization contain homologous gene copies from different evolutionary origins called homoeologs, whose orthologs must be sorted out in order to reconstruct the evolutionary history of polyploid clades. In this study, we propose a methodological approach to resolve the phylogeny of allopolyploid clades focusing on mesopolyploid genomes, which experienced some level of genome reshuffling and gene fractionation across their subgenomes. To illustrate our methodological framework we applied it to a clade belonging to the model Brassicaceae plant family, the Brassiceae tribe, that experienced a mesohexaploidy event. The dataset analysed consists of both publically available genomic sequences and new transcriptomic data according to taxa. The present methodology requires a well-annotated reference genome, for which the identification of the parental subgenome fragments has been performed (e.g. Brassica rapa and Brassica oleracea). Focusing on fully retained genes (i.e. genes for which all homoeologous gene copies inherited from the parental lineages are still present in the reference genome), the method constructs multi-labelled gene trees that allow subsequent assignment of each gene copy to its diploid parental lineage. Once the orthologous copies are identified, genes from the same parental origin are concatenated and tree-building methods are used to reconstruct the species tree. This method allows resolving the phylogenetic relationships (i) among extant species within a mesopolyploid clade, (ii) among the parental lineages of a mesopolyploid lineage, and (iii) between the parental lineages and closely related extant species. We report here the first well-resolved nuclear-based phylogeny of the Brassiceae tribe.

Small non-coding RNAs are central regulators of genome activity and stability. Their regulatory function typically involves sequence similarity with their target sites, but understanding the criteria by which they specifically recognize and regulate their targets across the genome remains a major challenge in the field, especially in the face of the diversity of silencing pathways involved. The dominance hierarchy among self-incompatibility alleles in Brassicaceae is controlled by interactions between a highly diversified set of small non-coding RNAs produced by dominant S-alleles and their corresponding target sites on recessive S-alleles. By controlled crosses, we created numerous heterozygous combinations of S-alleles inArabidopsis halleriand developed an RT-qPCR assay to compare allele-specific transcript levels for the pollen determinant of self-incompatibility (SCR). This provides the unique opportunity to evaluate the precise base-pairing requirements for effective transcriptional regulation of this target gene. We found strong transcriptional silencing of recessiveSCRalleles in all heterozygote combinations examined. A simple threshold model of base-pairing for the sRNA-target interaction captures most of the variation inSCRtranscript levels. For a subset of S-alleles, we also measured allele-specific transcript levels of the determinant of pistil specificity (SRK) and found sharply distinct expression dynamics throughout flower development betweenSCRandSRK. In contrast toSCR, bothSRKalleles were expressed at similar levels in the heterozygote genotypes examined, suggesting no transcriptional control of dominance for this gene. We discuss the implications for the evolutionary processes associated with the origin and maintenance of the dominance hierarchy among self-incompatibility alleles.

How two-components genetic systems accumulate evolutionary novelty and become diversified in the course of evolution is a fundamental problem in evolutionary systems biology. In the Brassicaceae, self-incompatibility (SI) is a spectacular example of a diversified allelic series in which numerous highly diverged receptor-ligand combinations are segregating in natural populations. However, the evolutionary mechanisms by which new SI specificities arise in the first place have remained elusive. Usingin plantaancestral protein resurrection, we demonstrate that two allelic variants currently segregating as distinct receptor-ligand combinations diverged through an asymmetrical process whereby one variant has retained the same recognition specificity as the (now extinct) ancestor, while the other has functionally diverged and now represents a novel specificity no longer recognized by the ancestor. Examination of the structural determinants of the shift in binding specificity suggests that allosteric changes may be an important source of evolutionary novelty in this highly diversified receptor-ligand system.

Plant self-incompatibility (SI) is a genetic system that prevents selfing and enforces outcrossing. Because of strong balancing selection, the genes encoding SI are predicted to maintain extraordinarily high levels of polymorphism, both in terms of the number of functionally distinct S-alleles that segregate in SI species and in terms of their nucleotide sequence divergence. However, because of these two combined features, documenting polymorphism of these genes also presents important methodological challenges that have so far largely prevented the comprehensive analysis of complete allelic series in natural populations, and also precluded the obtention of complete genic sequences for many S-alleles. Here, we develop a powerful methodological approach based on a computationally optimized comparison of short Illumina sequencing reads from genomic DNA to a database of known nucleotide sequences of the extracellular domain of SRK (eSRK). By examining mapping patterns along the reference sequences, we obtain highly reliable predictions of S-genotypes from individuals collected from natural populations of Arabidopsis halleri. Furthermore, using a de novo assembly approach of the filtered short reads, we obtain full-length sequences of eSRK even when the initial sequence in the database was only partial, and we discover putative new SRK alleles that were not initially present in the database. When including those new alleles in the reference database, we were able to resolve the complete diploid SI genotypes of all individuals. Beyond the specific case of Brassicaceae S-alleles, our approach can be readily applied to other polymorphic loci, given reference allelic sequences are available.

Invertebrates are exceptionally diverse, but many are in decline because of anthropogenic changes to their habitat. This situation is particularly problematic for taxa that are not well monitored or taxonomically poorly understood, because the lack of knowledge hampers conservation. Despite their important functional role in freshwater ecosystems, African bivalves of the family Unionidae remain poorly studied compared to their highly threatened relatives in Europe, the U.S.A. and Canada. To resolve relationships and to study diversification dynamics in space and time, we performed time-calibrated phylogenetic studies and biogeographical modeling on the unionids from the East African Rift System and surroundings, including representatives of all currently recognized Afrotropical genera except for Brazzaea (and Unio from southern Africa). Our analyses indicate that all sampled Afrotropical unionids belong to the tribe Coelaturini (subfamily Parreysiinae), as does the genus Moncetia from Lake Tanganyika, which is currently attributed to the family Iridinidae. Colonization of Africa from Eurasia by Parreysiinae occurred ~17 Ma ago, and the subsequent diversification of Coelaturini in Africa continued at a steady pace, although net diversification decreased over time as more niches and ecoregions became occupied. Clades in Coelaturini largely reflect drainage basins, with the oldest lineages and highest regional diversity occurring in Lake Tanganyika, followed by the Congo Basin watershed in general. The species assemblage of Lake Tanganyika reflects multiple independent events of colonization and intralacustrine diversification since the Late Miocene or Early Pliocene. The clades of other regions, including that containing the species from Lake Malawi, are comparatively young. Biogeographical analyses indicate that the colonization history was mainly driven by cladogenesis in sympatry, whereas few anagenetic events contributed to the modern distribution of Coelaturini. Ancestral range estimations demonstrate that Coelaturini originated in the Victoria and/or Tanganyika ecoregions, and that the Congo Basin played an essential role in the colonization of Africa by Coelaturini.

Plant self-incompatibility (SI) is a genetic system that prevents selfing and enforces outcrossing. Because of strong balancing selection, the genes encoding SI are predicted to maintain extraordinary high levels of polymorphism, both in terms of the number of S-alleles that segregate in SI species and in terms of nucleotide sequence divergence among distinct S-allelic lines. However, because of these two combined features, documenting polymorphism of these genes also presents important methodological challenges that have so far largely prevented the comprehensive analysis of complete allelic series in natural populations, and also precluded the obtention of complete genic sequences for many S-alleles. Here, we present a novel methodological approach based on a computationally optimized comparison of short Illumina sequencing reads from genomic DNA to a database of known nucleotide sequences of the extracellular domain ofSRK (eSRK). By examining mapping patterns along the reference sequences, we obtain highly reliable predictions of S-genotypes from individuals collected in natural populations ofArabidopsis halleri. Furthermore, using ade novoassembly approach of the filtered short reads, we obtain full length sequences of eSRK even when the initial sequence in the database was only partial, and we discover newSRKalleles that were not initially present in the database. When including those new alleles in the reference database, we were able to resolve the complete diploid SI genotypes of all individuals. Beyond the specific case of Brassicaceae S-alleles, our approach can be readily applied to other polymorphic loci, given reference allelic sequences are available.

Plant genomes are often characterized by a high level of repetitiveness and polyploid nature. Consequently, creating genome assemblies for plant genomes is challenging. The introduction of short-read technologies 10 years ago substantially increased the number of available plant genomes. Generally, these assemblies are incomplete and fragmented, and only a few are at the chromosome scale. Recently, Pacific Biosciences and Oxford Nanopore sequencing technologies were commercialized that can sequence long DNA fragments (kilobases to megabase) and, using efficient algorithms, provide high-quality assemblies in terms of contiguity and completeness of repetitive regions1–4. However, even though genome assemblies based on long reads exhibit high contig N50s (>1 Mb), these methods are still insufficient to decipher genome organization at the chromosome level. Here, we describe a strategy based on long reads (MinION or PromethION sequencers) and optical maps (Saphyr system) that can produce chromosome-level assemblies and demonstrate applicability by generating high-quality genome sequences for two new dicotyledon morphotypes, Brassica rapa Z1 (yellow sarson) and Brassica oleracea HDEM (broccoli), and one new monocotyledon, Musa schizocarpa (banana). All three assemblies show contig N50s of >5 Mb and contain scaffolds that represent entire chromosomes or chromosome arms. Assembling genomes to chromosome scale remains a challenge. Now, a study reports a strategy based on nanopore long reads and optical maps and uses it to produce high-quality chromosome-scale assemblies for the genomes of yellow sarson, broccoli and banana.

International audience; Human activities on the shoreline generate a growing pollution, creating deleterious habitats in coastal zones. Some species nevertheless succeed in such harsh milieus, raising the question of their tolerance to environmental stress. The annelid Hediste diversicolor lives buried in the sediments, directly exposed to contaminants trapped in the mud. After verifying the similarity of their genetic contexts, we compared reproductive output and individual immune resistance measures of populations living in polluted vs. ‘clean’ sediments, and related these assessments with measures of phthalates and metal pollution, and associated toxicity indices. Chemical analyses predicted no toxicity to the local infauna, and phenological studies evidenced no direct cost of living in noxious habitats. However, populations exposed to pollutants showed a significantly reduced survival upon infection with a local pathogen. Surprisingly, physiological studies evidenced a basal overinflammatory state in the most exposed populations. This over-activated baseline immune phenotype likely generates self-damage leading to enhanced immune cell death rate and immune failure. Monitoring the immune status of individual worms living in anthropic areas could thus be used as a reliable source of information regarding the actual health of wild populations.

Although the transition to selfing in the model plantArabidopsis thalianainvolved the loss of the self-incompatibility (SI) system, it clearly did not occur due to the fixation of a single inactivating mutation at the locus determining the specificities of SI (the S-locus). At least three groups of divergent haplotypes (haplogroups), corresponding to ancient functional S-alleles, have been maintained at this locus, and extensive functional studies have shown that all three carry distinct inactivating mutations. However, the historical process of loss of SI is not well understood, in particular its relation with the last glaciation. Here, we took advantage of recently published genomic re-sequencing data in 1,083Arabidopsis thalianaaccessions that we combined with BAC sequencing to obtain polymorphism information for the whole S-locus region at a species-wide scale. The accessions differed by several major rearrangements including large deletions and inter-haplogroup recombinations, forming a set of haplogroups that are widely distributed throughout the native range and largely overlap geographically. ‘Relict’A.thalianaaccessions that directly derive from glacial refugia are polymorphic at the S-locus, suggesting that the three haplogroups were already present when glacial refugia from the last Ice Age became isolated. Inter-haplogroup recombinant haplotypes were highly frequent, and detailed analysis of recombination breakpoints suggested multiple independent origins. These findings suggest that the complete loss of SI inA.thalianainvolved independent self-compatible mutants that arose prior to the last Ice Age, and experienced further rearrangements during post-glacial colonization.

International audience; The use of molecular data to reconstruct the history of divergence and gene flow between populations of closely related taxa represents a challenging problem. It has been proposed that the long-standing debate about the geography of speciation can be resolved by comparing the likelihoods of a model of isolation with migration and a model of secondary contact. However, data are commonly only fit to a model of isolation with migration and rarely tested against the secondary contact alternative. Furthermore, most demographic inference methods have neglected variation in introgression rates and assume that the gene flow parameter (Nm) is similar among loci. Here, we show that neglecting this source of variation can give misleading results. We analysed DNA sequences sampled from populations of the marine mussels, Mytilus edulis and M. galloprovincialis, across a well-studied mosaic hybrid zone in Europe and evaluated various scenarios of speciation, with or without variation in introgression rates, using an Approximate Bayesian Computation (ABC) approach. Models with heterogeneous gene flow across loci always outperformed models assuming equal migration rates irrespective of the history of gene flow being considered. By incorporating this heterogeneity, the best-supported scenario was a long period of allopatric isolation during the first three-quarters of the time since divergence followed by secondary contact and introgression during the last quarter. By contrast, constraining migration to be homogeneous failed to discriminate among any of the different models of gene flow tested. Our simulations thus provide statistical support for the secondary contact scenario in the European Mytilus hybrid zone that the standard coalescent approach failed to confirm. Our results demonstrate that genomic variation in introgression rates can have profound impacts on the biological conclusions drawn from inference methods and needs to be incorporated in future studies.

International audience; Evolutionary transitions between mating systems have occurred repetitively and independently in flowering plants. One of the most spectacular advances of the recent empirical literature in the field was the discovery of the underlying genetic machinery, which provides the opportunity to retrospectively document the scenario of the outcrossing to selfing transitions in a phylogenetic perspective. In this review, we explore the literature describing patterns of polymorphism and molecular evolution of the locus controlling self-incompatibility (S-locus) in selfing species of the Brassicaceae family in order to document the transition from outcrossing to selfing, a retrospective approach that we describe as the 'mating system genes approach'. The data point to strikingly contrasted scenarios of transition from outcrossing to selfing. We also perform original analyses of the fully sequenced genomes of four species showing self-compatibility, to compare the orthologous S-locus region with that of functional S-locus haplotypes. Phylogenetic analyses suggest that all species we investigated evolved independently towards loss of self-incompatibility, and in most cases almost intact sequences of either of the two S-locus genes suggest that these transitions occurred relatively recently. The S-locus region in Aethionema arabicum, representing the most basal lineage of Brassicaceae, showed unusual patterns so that our analysis could not determine whether self-incompatibility was lost secondarily, or evolved in the core Brassicaceae after the split with this basal lineage. Although the approach we detail can only be used when mating system genes have been identified in a clade, we suggest that its integration with phylogenetic and population genetic approaches should help determine the main routes of this predominant mating system shift in plants.

UMR AGAP - équipe AFEF - Architecture et fonctionnement des espèces fruitières; International audience; Bervillé et al. express concern about the existence of the diallelic self-incompatibility (DSI) system in Olea europaea, mainly because our model does not account for results from previous studies from their group that claimed to have documented asymmetry of the incompatibility response in reciprocal crosses. In this answer to their comment, we present original results based on reciprocal stigma tests that contradict conclusions from these studies. We show that, in our hands, not a single case of asymmetry was confirmed, endorsing that symmetry of incompatibility reactions seems to be the rule in Olive. We discuss three important aspects that were not taken into account in the studies cited in their comments and that can explain the discrepancy: (i) the vast uncertainty around the actual genetic identity of vernacular varieties, (ii) the risk of massive contamination associated with the pollination protocols that they used and (iii) the importance of checking for stigma receptivity in controlled crosses. These studies were thus poorly genetically controlled, and we stand by our original conclusion that Olive tree exhibits DSI.

Summary •Arabidopsis halleri is a pseudometallophyte involved in numerous molecular studies of the adaptation to anthropogenic metal stress. In order to test the representativeness of genetic accessions commonly used in these studies, we investigated the A. halleri population genetic structure in Europe. •Microsatellite and nucleotide polymorphisms from the nuclear and chloroplast genomes, respectively, were used to genotype 65 populations scattered over Europe. •The large-scale population structure was characterized by a significant phylogeographic signal between two major genetic units. The localization of the phylogeographic break was assumed to result from vicariance between large populations isolated in southern and central Europe, on either side of ice sheets covering the Alps during the Quaternary ice ages. Genetic isolation was shown to be maintained in western Europe by the high summits of the Alps, whereas admixture was detected in the Carpathians. •Considering the phylogeographic literature, our results suggest a distinct phylogeographic pattern for European species occurring in both mountain and lowland habitats. Considering the evolution of metal adaptation in A. halleri, it appears that recent adaptations to anthropogenic metal stress that have occurred within either phylogeographic unit should be regarded as independent events that potentially have involved the evolution of a variety of genetic mechanisms.

International audience

Homomorphic self-incompatibility (SI) evolved in many plant families to enforce selfing avoidance, and is controlled by a single multiallelic locus (the S-locus). In a fragmented landscape, strong variation in population size and in local density is expected to cause strong variation in allelic diversity at the S-locus, which could generate an Allee effect on female reproductive success by constraining compatible pollen availability. In this experimental study, we aimed at detecting this SI-specific Allee effect (or S-Allee effect) in the endangered species Biscutella neustriaca. We demonstrated the occurrence of a SI mating system in the species and determined compatibility relationships among genotypes through a large set of controlled pollinations. For the experiment, we chose three different pollen receptor genotypes, each compatible with respectively 100, 75 and 25% of four other genotypes, which constituted the pollen sources. We placed different ramets of each receptor at different distances from the pollen sources to control for pollen limitation due to low local density, and we measured the seed set on each receptor plant three times consecutively. Analyses performed with generalized linear mixed models showed that both the distance to the pollen sources and the mate availability due to SI had a significant effect on seed set, with a strong reduction observed when mate availability was limited to 25%. Our results suggest that pollen limitation due to a restriction in compatible mate availability could occur in small or scattered populations exhibiting low allelic diversity at the S-locus.

Recent theoretical advances have suggested that various forms of balancing selection may promote the evolution of dominance through an increase of the proportion of heterozygote genotypes. We test whether dominance can evolve in the sporophytic self-incompatibility (SSI) system in plants. SSI prevents mating between individuals expressing identical SI phenotypes by recognition of pollen by pistils, which avoids selfing and inbreeding depression. SI phenotypes depend on a complex network of dominance relationships between alleles at the self-incompatibility locus (S-locus). Empirical studies suggest that these relationships are not random, but the exact evolutionary processes shaping these relationships remain unclear. We investigate the expected patterns of dominance under the hypothesis that dominance is a direct target of natural selection. We follow the fate of a mutant allele at the S-locus whose dominance relationships are changed but whose specificity remains unaltered. We show that strict codominance is not evolutionarily stable in SSI, and that inbreeding depression due to deleterious mutations linked or unlinked to the S-locus exerts strong constraints on changes in relative levels of dominance in pollen and pistil. Our results provide a general adaptive explanation for most patterns of dominance relationships empirically observed in natural plant populations.

Our understanding of the impact of recombination, mutation, genetic drift, and selection on the evolution of a single gene is still limited. Here we investigate the impact of all these evolutionary forces at the complementary sex determiner (csd) gene that evolves under a balancing mode of selection. Females are heterozygous at the csd gene and males are hemizygous; diploid males are lethal and occur when csd is homozygous. Rare alleles thus have a selective advantage, are seldom lost by the effect of genetic drift, and are maintained over extended periods of time when compared with neutral polymorphisms. Here, we report on the analysis of 17, 19, and 15 csd alleles of Apis cerana, Apis dorsata, and Apis mellifera honeybees, respectively. We observed great heterogeneity of synonymous (piS) and nonsynonymous (piN) polymorphisms across the gene, with a consistent peak in exons 6 and 7. We propose that exons 6 and 7 encode the potential specifying domain (csd-PSD) that has accumulated elevated nucleotide polymorphisms over time by balancing selection. We observed no direct evidence that balancing selection favors the accumulation of nonsynonymous changes at csd-PSD (piN/piS ratios are all

International audience; A recent investigation found evidence that the transition of Arabidopsis thaliana from ancestral self-incompatibility (SI) to full self-compatibility occurred very recently and suggested that this occurred through a selective fixation of a nonfunctional allele (WSCR1) at the SCR gene, which determines pollen specificity in the incompatibility response. The main evidence is the lack of polymorphism at the SCR locus in A. thaliana. However, the nearby SRK gene, which determines stigma specificity in self-incompatible Brassicaceae species, has extremely high sequence diversity, with 3 very divergent SRK haplotypes, 2 of them present inmultiple strains. Such high diversity is extremely unusual in this species, and it suggests the possibility that multiple, different SRK haplotypes may have been preserved fromA. thaliana's self-incompatible ancestor. To study the evolution of S-haplotypes in the A. thaliana lineage, we searched the 2 most closely related Arabidopsis species Arabidopsis lyrata and Arabidopsis halleri, in which most populations have retained SI, and found SRK sequences corresponding to all 3 A. thaliana haplogroup sequences. Our molecular evolutionary analyses of these 3 S-haplotypes provide an independent estimate of the timing of the breakdown of SI and again exclude an ancient transition to selfing in A. thaliana. Comparing sequences of each of the 3 haplogroups between species, we find that 2 of the 3 SRK sequences (haplogroups A and B) are similar throughout their length, suggesting that little or no recombination with other SRK alleles has occurred since these species diverged. The diversity difference between the SCR and SRK loci in A. thaliana, however, suggests crossingover, either within SRK or between the SCR and SRK loci. If the loss of SI involved fixation of the WSCR1 sequence, the exchange must have occurred during its fixation. Divergence between the species is much lower at the S-locus, compared with reference loci, andwe discuss two contributory possibilities. Introgression may have occurred between A. lyrata and A. halleri and between their ancestral lineage and A. thaliana, at least for some period after their split. In addition, the coalescence times of sequences of individual S-haplogroups are expected to be less than those of alleles at non–S-loci.

We used chloroplast polymerase chain reaction-restriction-fragment length polymorphism (PCR-RFLP) and chloroplast microsatellites to assess the structure of genetic variation and postglacial history across the entire natural range of the common ash (Fraxinus excelsior L.), a broad-leaved wind-pollinated and wind-dispersed European forest tree. A low level of polymorphism was observed, with only 12 haplotypes at four polymorphic microsatellites in 201 populations, and two PCR-RFLP haplotypes in a subset of 62 populations. The clear geographical pattern displayed by the five most common haplotypes was in agreement with glacial refugia for ash being located in Iberia, Italy, the eastern Alps and the Balkan Peninsula, as had been suggested from fossil pollen data. A low chloroplast DNA mutation rate, a low effective population size in glacial refugia related to ash’s life history traits, as well as features of postglacial expansion were put forward to explain the low level of polymorphism. Differentiation among populations was high (G ST = 0.89), reflecting poor mixing among recolonizing lineages. Therefore, the responsible factor for the highly homogeneous genetic pattern previously identified at nuclear microsatellites throughout western and central Europe (Heuertz et al. 2004) must have been efficient postglacial pollen flow. Further comparison of variation patterns at both marker systems revealed that nuclear microsatellites identified complex differentiation patterns in south-eastern Europe which remained undetected with chloroplast microsatellites. The results suggest that data from different markers should be combined in order to capture the most important genetic patterns in a species.

Self-incompatibility systems in plants are genetic systems that prevent self-fertilization in hermaphrodites through recognition and rejection of pollen expressing the same allelic specificity as that expressed in the pistils. The evolutionary properties of these self-recognition systems have been revealed through a fascinating interplay between empirical advances and theoretical developments. In 1939, Wright suggested that the main evolutionary force driving the genetic and molecular properties of these systems was strong negative frequency-dependent selection acting on pollination success. The empirical observation of high allelic diversity at the self-incompatibility locus in several species, followed by the discovery of very high molecular divergence among alleles in all plant families where the locus has been identified, supported Wright's initial theoretical predictions as well as many of its later developments. In the last decade, however, advances in the molecular characterization of the incompatibility reaction and in the analysis of allelic frequencies and allelic divergence from natural populations have stimulated new theoretical investigations that challenged some important assumptions of Wright's model of gametophytic self-incompatibility. We here review some of these recent empirical and theoretical advances that investigated: (i) the hypothesis that S-alleles are selectively equivalent, and the evolutionary consequences of genetic interactions between alleles; (ii) the occurrence of frequency-dependent selection in female fertility; (iii) the evolutionary genetics of self-incompatibility systems in subdivided populations; (iv) the evolutionary implications of the self-incompatibility locus's genetic architecture; and (v) of its interactions with the genomic environment.

Intraspecific variation in nuclear DNA content in 19 populations of Armeria maritima was investigated. Significant differences in DNA amount among ecogeographic groups and among populations were observed. These differences are related to geographic origin, with Californian populations having 7 % more DNA than their European counterparts. Within Europe, the four ecogeographic groups tested are not distinct, whatever their geographic origin (maritime versus continental) and their ecological characteristics (heavy-metal polluted versus non-toxic soils). However, significant differences in DNA amount appeared among population clusters identified from biochemical classification of European populations based on flavonoid compounds and isozyme markers. Variation in genome size among populations was found to be significantly correlated to variation in pollen size. Weak interpopulational variation has been detected for the nucleotidic base ratio, but was not correlated to variation in DNA amount.

Many empirical studies have assessed fine-scale spatial genetic structure (SGS), i.e. the nonrandom spatial distribution of genotypes, within plant populations using genetic markers and spatial autocorrelation techniques. These studies mostly provided qualitative descriptions of SGS, rendering quantitative comparisons among studies difficult. The theory of isolation by distance can predict the pattern of SGS under limited gene dispersal, suggesting new approaches, based on the relationship between pairwise relatedness coefficients and the spatial distance between individuals, to quantify SGS and infer gene dispersal parameters. Here we review the theory underlying such methods and discuss issues about their application to plant populations, such as the choice of the relatedness statistics, the sampling scheme to adopt, the procedure to test SGS, and the interpretation of spatial autocorrelograms. We propose to quantify SGS by an ‘ Sp ’ statistic primarily dependent upon the rate of decrease of pairwise kinship coefficients between individuals with the logarithm of the distance in two dimensions. Under certain conditions, this statistic estimates the reciprocal of the neighbourhood size. Reanalysing data from, mostly, published studies, the Sp statistic was assessed for 47 plant species. It was found to be significantly related to the mating system (higher in selfing species) and to the life form (higher in herbs than trees), as well as to the population density (higher under low density). We discuss the necessity for comparing SGS with direct estimates of gene dispersal distances, and show how the approach presented can be extended to assess (i) the level of biparental inbreeding, and (ii) the kurtosis of the gene dispersal distribution.

In Belgium, at the western border of its continental distribution range, the perennial herb S. nutans has evolved two parapatric ecotypes (calcicolous or silicicolous), which show contrasted allozyme, morphological and adaptive patterns, and isolating reproductive barriers. We examined allozyme variation in 21 populations close to Belgium (from France and Luxemburg) in relation to their peripheral or central geographical position, and investigated their genetic relationships with the Belgian ecotypes. Both peripheral and central populations showed high genetic variation. They were differentiated from each other, but not in relation to the substrate (calcareous or siliceous). The peripheral and central populations were related to the Belgian silicicolous and calcicolous ecotype, respectively. This suggests an ancient differentiation and different past histories for the Belgian ecotypes, and parapatric speciation processes. This study exemplifies the high evolutionary potential of populations at range peripheries and the need for considering them in conservation strategies.

The cultivation area of industrial chicory, Cichorium intybus L. cv Sativum, coincides with the natural distribution area of its wild relative, C. intybus L., which could lead to gene flow between wild and cultivated types. The genetic diversity within and between the two types has therefore been studied using AFLP genotyping of samples from 12 wild populations collected in Belgium and ten commercial varieties. The genotyping of 233 individuals allowed the identification of 254 AFLP markers. Similar levels of genetic diversity were observed within wild populations and cultivated varieties, suggesting the absence of any strong bottleneck in the history of the cultivated types. The phylogenetic analysis pointed to a monophyletic origin of cultivated varieties as compared to the local wild populations studied, hence the two types of chicory form two separate gene pools. The genotyping of some individuals sampled in ruderal sites clearly showed that they belong to the cultivated gene pool, which suggests the existence of feral or weedy types. The low differentiation observed among wild populations indicates that gene flow might be important in this species.

An extensive allozyme survey was conducted within a natural ''meta'' population of the native North American annual legume, Chamaecrista fasciculata (Leguminosae) to quantify genetic structure at different spatial scales. Gene flow was then estimated by a recently developed indirect method based on a continuous population model, using pairwise kinship coefficients between individuals. The indirect estimates of gene flow, quantified in terms of neighborhood size, with an average value on the order of 150 individuals, were concordant among different spatial scales (subpopulation, population, metapopulation). This gene-flow value lies within the range of direct estimates previously documented from observations of pollen and seed dispersal for the same metapopulation. Monte Carlo simulations using the direct measures of gene flow as parameters further demonstrated that the observed spatial pattern of allozyme variation was congruent with a model of isolation by distance. Combining previously published estimates of pollen dispersal distances with kinship coefficients from this study, we quantified biparental inbreeding relative to either a single subpopulation or the whole metapopulation. At the level of a neighborhood, little biparental inbreeding was observed and most departure from Hardy-Weinberg genotypic proportions was explained by self-fertilization, whereas both selfing and biparental inbreeding contributed to nonrandom mating at the metapopulation level. Gene flow was also estimated from indirect methods based on a discontinuous population structure model. We discuss these results with respect to the effect of a patchy population structure on estimation of gene flow.

The Centaurea jacea aggregate is a polymorphic polyploid complex whose taxonomic treatment is still controversial. A numerical taxonomic approach was applied to 394 individuals of known ploidy level, from 19 populations, based on the main diagnostic characters proposed in earlier revisions. Populations from xeric grasslands were not considered. Principal Component Analysis shows that variation within the complex is continuous. UPGMA Cluster Analysis based on population means supports the recognition of three groups of populations. However, the limits between these groups are blurred to a considerable extent due to extensive within-population polymorphism. It is argued that the Belgian populations of Centaurea jacea occurring in mesic grasslands should be treated as a single species, with three subspecies. The two extremes of the morphological gradient can be referred to as C. jacea subsp. jacea and C. jacea subsp. nigra, with C. jacea subsp. pratensis occupying an intermediate position. Most populations from Belgium are tetraploid, a diploid chromosome number being found only in populations of C. jacea subsp. nigra from the Ardennes massif. On average, diploids grow at higher altitude and on more acidic soils than tetraploids. Finally, a key to the three subspecies is provided., SCOPUS: ar.j, FLWIN, info:eu-repo/semantics/published

We investigate the distribution of sizes of fragments obtained from the amplified fragment length polymorphism (AFLP) marker technique. We find that empirical distributions obtained in two plant species, Phaseolus lunatus and Lolium perenne, are consistent with the expected distributions obtained from analytical theory and from numerical simulations. Our results indicate that the size distribution is strongly asymmetrical, with a much higher proportion of small than large fragments, that it is not influenced by the number of selective nucleotides nor by genome size but that it may vary with genome-wide GC-content, with a higher proportion of small fragments in cases of lower GC-content when considering the standard AFLP protocol with the enzyme MseI. Results from population samples of the two plant species show that there is a negative relationship between AFLP fragment size and fragment population frequency. Monte Carlo simulations reveal that size homoplasy, arising from pulling together nonhomologous fragments of the same size, generates patterns similar to those observed in P. lunatus and L. perenne because of the asymmetry of the size distribution. We discuss the implications of these results in the context of estimating genetic diversity with AFLP markers.

We analysed genetic variation within and between populations of the common ash from Bulgaria in order to extract biological information useful in the context of conservation management of eastern European genetic resources of noble hardwood species. A total of 321 trees from three regions of Bulgaria were typed at six highly polymorphic microsatellite loci. Analysis of within-population inbreeding suggests an upper boundary value of 2.7% for the selfing rate. Significant spatial genetic structure consistent with models of isolation by distance was detected within four out of 10 populations as well as among populations. Estimates of neighbourhood size in the range 38-126 individual trees were obtained based on spatial genetic structure analyses at either the intrapopulation or interpopulation level. Differentiation among populations explained only about 8.7% of total genetic diversity. These results are discussed in comparison with data from social broad-leaved species such as oak and beech.

Simulations are used to investigate the expected pattern of variation at loci under different forms of multi-allelic balancing selection in a finite island model of a subdivided population. The objective is to evaluate the effect of restricted migration among demes on the distribution of polymorphism at the selected loci at equilibrium, and to compare the results with those expected for a neutral locus. The results show that the expected number of alleles maintained, and numbers of nucleotide differences between alleles, are relatively insensitive to the migration rate, and differentiation remains low even under very restricted migration. However, nucleotide divergence between copies of functionally identical alleles increases sharply when migration decreases. These results are discussed in relation to published surveys of allelic diversity in MHC and plant self-incompatibility systems, and to the possibility of inferring ancient population genetic events and processes. In addition, it is shown that, for sporophytic self-incompatibility systems, it is not necessarily true in a subdivided population that recessive alleles are more frequent than dominant ones. Udgivelsesdato: 2000-Aug

International audience; Microsatellite markers are increasingly being used in crop plants to discriminate among genotypes and as tools in marker-assisted selection. Here we evaluated the use of microsatellite markers to quantify the genetic diversity within as well as among accessions sampled from the world germplasm collection of sorghum. Considerable variation was found at the five microsatellite loci analysed, with an average number of alleles per locus equal to 2.4 within accessions and 19.2 in the overall sample of 25 accessions. The collection of sorghum appeared highly structured genetically with about 70% of the total genetic diversity occurring among accessions. However, differentiation among morphologically defined races of sorghum, or among geographic origins, accounted for less than 15% of the total genetic diversity. Our results are in global agreement with those obtained previously with allozyme markers. We were also able to show that microsatellite data are useful in identifying individual accessions with a high relative contribution to the overall allelic diversity of the collection.

Analysis of the spatial genetic structure within continuous populations in their natural habitat can reveal acting evolutionary processes. Spatial autocorrelation statistics are often used for this purpose, but their relationships with population genetics models have not been thoroughly established. Moreover, it has been argued that the dependency of these statistics on variation in mutation rates among loci strongly limits their interest for inferential purposes. In the context of an isolation by distance process, we describe relationships between a descriptor of the spatial genetic structure used in empirical studies, Moran's I statistic and population genetics parameters. In particular, we point out that, when Moran's I statistic is used to describe correlation in allele frequencies at the individual level, it provides an estimator of Wright's coefficient of relationship. We also show that the latter parameter, as a descriptor of genetic structure, is not influenced by selfing rate or ploidy level. Under specific finite population models, numerical simulations show that values of Moran's I statistic can be predicted from analytical theory. These simulations are also used to estimate the time taken to approach a structure at equilibrium. Finally, we discuss the conditions under which spatial autocorrelation statistics are little influenced by variation in mutation rates, so that they could be used to estimate gene dispersal parameters.