Vincent Castric

Does Speciation between Arabidopsis halleri and Arabidopsis lyrata Coincide with Major Changes in a Molecular Target of Adaptation

Ever since Darwin proposed natural selection as the driving force for the origin of species, the role of adaptive processes in speciation has remained controversial. In particular, a largely unsolved issue is whether key divergent ecological adaptations are associated with speciation events or evolve secondarily within sister species after the split. The plant Arabidopsis halleri is one of the few species able to colonize soils highly enriched in zinc and cadmium. Recent advances in the molecular genetics of adaptation show that the physiology of this derived ecological trait involves copy number expansions of the AhHMA4 gene, for which orthologs are found in single copy in the closely related A. lyrata and the outgroup A. thaliana. To gain insight into the speciation process, we ask whether adaptive molecular changes at this candidate gene were contemporary with important stages of the speciation process. We first inferred the scenario and timescale of speciation by comparing patterns of variation across the genomic backgrounds of A. halleri and A. lyrata. Then, we estimated the timing of the first duplication of AhHMA4 in A. halleri. Our analysis suggests that the historical split between the two species closely coincides with major changes in this molecular target of adaptation in the A. halleri lineage. These results clearly indicate that these changes evolved in A. halleri well before industrial activities fostered the spread of Zn- and Cd-polluted areas, and suggest that adaptive processes related to heavy-metal homeostasis played a major role in the speciation process.

Nuclear and chloroplast DNA phylogeography reveals vicariance among European populations of the model species for the study of metal tolerance, Arabidopsis halleri (Brassicaceae)

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.

DNA binding of the cell cycle transcriptional regulator GcrA depends on N6-adenosine methylation in Caulobacter crescentus and other Alphaproteobacteria

Several regulators are involved in the control of cell cycle progression in the bacterial model system Caulobacter crescentus, which divides asymmetrically into a vegetative G1-phase (swarmer) cell and a replicative S-phase (stalked) cell. Here we report a novel functional interaction between the enigmatic cell cycle regulator GcrA and the N6-adenosine methyltransferase CcrM, both highly conserved proteins among Alphaproteobacteria, that are activated early and at the end of S-phase, respectively. As no direct biochemical and regulatory relationship between GcrA and CcrM were known, we used a combination of ChIP (chromatin-immunoprecipitation), biochemical and biophysical experimentation, and genetics to show that GcrA is a dimeric DNA–binding protein that preferentially targets promoters harbouring CcrM methylation sites. After tracing CcrM-dependent N6-methyl-adenosine promoter marks at a genome-wide scale, we show that these marks recruit GcrA in vitro and in vivo. Moreover, we found that, in the presence of a methylated target, GcrA recruits the RNA polymerase to the promoter, consistent with its role in transcriptional activation. Since methylation-dependent DNA binding is also observed with GcrA orthologs from other Alphaproteobacteria, we conclude that GcrA is the founding member of a new and conserved class of transcriptional regulators that function as molecular effectors of a methylation-dependent (non-heritable) epigenetic switch that regulates gene expression during the cell cycle.

A general model to explore complex dominance patterns in plant sporophytic self-incompatibility systems

We developed a general model of sporophytic self-incompatibility under negative frequency-dependent selection allowing complex patterns of dominance among alleles. We used this model deterministically to investigate the effects on equilibrium allelic frequencies of the number of dominance classes, the number of alleles per dominance class, the asymmetry in dominance expression between pollen and pistil, and whether selection acts on male fitness only or both on male and on female fitnesses. We show that the so-called “recessive effect” occurs under a wide variety of situations. We found emerging properties of finite population models with several alleles per dominance class such as that higher numbers of alleles are maintained in more dominant classes and that the number of dominance classes can evolve. We also investigated the occurrence of homozygous genotypes and found that substantial proportions of those can occur for the most recessive alleles. We used the model for two species with complex dominance patterns to test whether allelic frequencies in natural populations are in agreement with the distribution predicted by our model. We suggest that the model can be used to test explicitly for additional, allele-specific, selective forces.

Recent and Ancient Signature of Balancing Selection around the S-Locus in Arabidopsis halleri and A. lyrata

Balancing selection can maintain different alleles over long evolutionary times. Beyond this direct effect on the molecular targets of selection, balancing selection is also expected to increase neutral polymorphism in linked genome regions, in inverse proportion to their genetic map distances from the selected sites. The genes controlling plant self-incompatibility are subject to one of the strongest forms of balancing selection, and they show clear signatures of balancing selection. The genome region containing those genes (the S-locus) is generally described as nonrecombining, and the physical size of the region with low recombination has recently been established in a few species. However, the size of the region showing the indirect footprints of selection due to linkage to the S-locus is only roughly known. Here, we improved estimates of this region by surveying synonymous polymorphism and estimating recombination rates at 12 flanking region loci at known physical distances from the S-locus region boundary, in two closely related self-incompatible plants Arabidopsis halleri and A. lyrata. In addition to studying more loci than previous studies and using known physical distances, we simulated an explicit demographic scenario for the divergence between the two species, to evaluate the extent of the genomic region whose diversity departs significantly from neutral expectations. At the closest flanking loci, we detected signatures of both recent and ancient indirect effects of selection on the S-locus flanking genes, finding ancestral polymorphisms shared by both species, as well as an excess of derived mutations private to either species. However, these effects are detected only in a physically small region, suggesting that recombination in the flanking regions is sufficient to quickly break up linkage disequilibrium with the S-locus. Our approach may be useful for distinguishing cases of ancient versus recently evolved balancing selection in other systems.

Contrasted patterns of molecular evolution in dominant and recessive self-incompatibility haplotypes in Arabidopsis.

Self-incompatibility has been considered by geneticists a model system for reproductive biology and balancing selection, but our understanding of the genetic basis and evolution of this molecular lock-and-key system has remained limited by the extreme level of sequence divergence among haplotypes, resulting in a lack of appropriate genomic sequences. In this study, we report and analyze the full sequence of eleven distinct haplotypes of the self-incompatibility locus (S-locus) in two closely related Arabidopsis species, obtained from individual BAC libraries. We use this extensive dataset to highlight sharply contrasted patterns of molecular evolution of each of the two genes controlling self-incompatibility themselves, as well as of the genomic region surrounding them. We find strong collinearity of the flanking regions among haplotypes on each side of the S-locus together with high levels of sequence similarity. In contrast, the S-locus region itself shows spectacularly deep gene genealogies, high variability in size and gene organization, as well as complete absence of sequence similarity in intergenic sequences and striking accumulation of transposable elements. Of particular interest, we demonstrate that dominant and recessive S-haplotypes experience sharply contrasted patterns of molecular evolution. Indeed, dominant haplotypes exhibit larger size and a much higher density of transposable elements, being matched only by that in the centromere. Overall, these properties highlight that the S-locus presents many striking similarities with other regions involved in the determination of mating-types, such as sex chromosomes in animals or in plants, or the mating-type locus in fungi and green algae.

Can we continue to neglect genomic variation in introgression rates when inferring the history of speciation? A case study in a Mytilus hybrid zone.

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.

Does Frequency-Dependent Selection with Complex Dominance Interactions Accurately Predict Allelic Frequencies at the Self-Incompatibility Locus in Arabidopsis halleri?

Abstract Frequency-dependent selection is a major force determining the evolutionary dynamics of alleles at the self-incompatibility locus (S-locus) in flowering plants. We introduce a general method using numerical simulations to test several alternative models of frequency-dependent selection on S-locus data from sporophytic systems, taking into account both genetic drift and observed patterns of dominance interactions among S-locus haplotypes (S-haplotypes). Using a molecular typing method, we estimated S-haplotype frequencies in a sample of 322 adult plants and of 245 offspring obtained from seeds sampled on 22 maternal plants, collected in a single population of Arabidopsis halleri (Brassicaceae). We found eight different S-haplotypes and characterized their dominance interactions by controlled pollinations. We then compared the likelihood of different models of frequency-dependent selection: we found that the observed haplotype frequencies and observed frequency changes in one generation best fitted a model with (1) the observed dominance interactions and (2) no pollen limitation. Overall, our population genetic models of frequency-dependent selection, including patterns of dominance interactions among S-haplotypes and genetic drift, can reliably predict polymorphism at the S-locus. We discuss how these approaches allow detecting additional processes influencing the evolutionary dynamics of the S-locus, such as purifying selection on linked loci.

Dominance hierarchy arising from the evolution of a complex small RNA regulatory network

The prevention of fertilization through self-pollination (or pollination by a close relative) in the Brassicaceae plant family is determined by the genotype of the plant at the self-incompatibility locus (S locus). The many alleles at this locus exhibit a dominance hierarchy that determines which of the two allelic specificities of a heterozygous genotype is expressed at the phenotypic level. Here, we uncover the evolution of how at least 17 small RNA (sRNA)–producing loci and their multiple target sites collectively control the dominance hierarchy among alleles within the gene controlling the pollen S-locus phenotype in a self-incompatible Arabidopsis species. Selection has created a dynamic repertoire of sRNA-target interactions by jointly acting on sRNA genes and their target sites, which has resulted in a complex system of regulation among alleles. Hierarchical interactions among alleles in a mustard plant explain how self-incompatibility evolved and is maintained. Dominance cascades in self-incompatibility Plants often cannot use their own pollen to set seed. This is known as self-incompatibility. Although some of the underlying genetics and mechanisms of self-incompatibility are understood, the evolution and maintenance of the system have remained mysterious. Durand et al. identified a collection of small RNAs and their respective matching targets within a self-incompatibility locus in Arabidopsis halleri. A subset of these alleles functioned in a dominant manner, which helps to explain how self-incompatibility is maintained. Science, this issue p. 1200

Origin and Diversification Dynamics of Self-Incompatibility Haplotypes

Self-incompatibility (SI) is a genetic system found in some hermaphrodite plants. Recognition of pollen by pistils expressing cognate specificities at two linked genes leads to rejection of self pollen and pollen from close relatives, i.e., to avoidance of self-fertilization and inbred matings, and thus increased outcrossing. These genes generally have many alleles, yet the conditions allowing the evolution of new alleles remain mysterious. Evolutionary changes are clearly necessary in both genes, since any mutation affecting only one of them would result in a nonfunctional self-compatible haplotype. Here, we study diversification at the S-locus (i.e., a stable increase in the total number of SI haplotypes in the population, through the incorporation of new SI haplotypes), both deterministically (by investigating analytically the fate of mutations in an infinite population) and by simulations of finite populations. We show that the conditions allowing diversification are far less stringent in finite populations with recurrent mutations of the pollen and pistil genes, suggesting that diversification is possible in a panmictic population. We find that new SI haplotypes emerge fastest in populations with few SI haplotypes, and we discuss some implications for empirical data on S-alleles. However, allele numbers in our simulations never reach values as high as observed in plants whose SI systems have been studied, and we suggest extensions of our models that may reconcile the theory and data.