Hunting for beneficial mutations: conditioning on SIFT scores when estimating the distribution of fitness effect of new mutations

The Distribution of Fitness Effects (DFE) of new mutations is a key parameter of molecular evolution. The DFE can in principle be estimated by comparing the Site Frequency Spectra (SFS) of putatively neutral and functional polymorphisms. Unfortunately the DFE is intrinsically hard to estimate, especially for beneficial mutations since these tend to be exceedingly rare. There is therefore a strong incentive to find out whether conditioning on properties of mutations that are independent of the SFS could provide additional information. In the present study, we developed a new measure based on SIFT scores. SIFT scores are assigned to nucleotide sites based on their level of conservation across a multi species alignment: the more conserved a site, the more likely mutations occurring at this site are deleterious and the lower the SIFT score. If one knows the ancestral state at a given site, one can assign a value to new mutations occurring at the site based on the change of SIFT score associated with the mutation. We called this new measure δ. We show that properties of the DFE as well as the flux of beneficial mutations across classes covary with δ and, hence, that SIFT scores are informative when estimating the fitness effect of new mutations. In particular, conditioning on SIFT scores can help to characterize beneficial mutations.

Genomic signature of successful colonization of Eurasia by the allopolyploid shepherd's purse (Capsella bursa-pastoris)

Polyploidization is a dominant feature of flowering plant evolution. However, detailed genomic analyses of the interpopulation diversification of polyploids following genome duplication are still in their infancy, mainly because of methodological limits, both in terms of sequencing and computational analyses. The shepherd's purse (Capsella bursa-pastoris) is one of the most common weed species in the world. It is highly self-fertilizing, and recent genomic data indicate that it is an allopolyploid, resulting from hybridization between the ancestors of the diploid species Capsella grandiflora and Capsella orientalis. Here, we investigated the genomic diversity of C. bursa-pastoris, its population structure and demographic history, following allopolyploidization in Eurasia. To that end, we genotyped 261 C. bursa-pastoris accessions spread across Europe, the Middle East and Asia, using genotyping-by-sequencing, leading to a total of 4274 SNPs after quality control. Bayesian clustering analyses revealed three distinct genetic clusters in Eurasia: one cluster grouping samples from Western Europe and Southeastern Siberia, the second one centred on Eastern Asia and the third one in the Middle East. Approximate Bayesian computation (ABC) supported the hypothesis that C. bursa-pastoris underwent a typical colonization history involving low gene flow among colonizing populations, likely starting from the Middle East towards Europe and followed by successive human-mediated expansions into Eastern Asia. Altogether, these findings bring new insights into the recent multistage colonization history of the allotetraploid C. bursa-pastoris and highlight ABC and genotyping-by-sequencing data as promising but still challenging tools to infer demographic histories of selfing allopolyploids.

Speciation of two desert poplar species triggered by Pleistocene climatic oscillations

Despite the evidence that the Pleistocene climatic fluctuations have seriously affected the distribution of intraspecific diversity, less is known on its impact on interspecific divergence. In this study, we aimed to test the hypothesis that the divergence of two desert poplar species Populus euphratica Oliv. and P. pruinosa Schrenk. occurred during the Pleistocene. We sequenced 11 nuclear loci in 60 individuals from the two species to estimate the divergence time between them and to test whether gene flow occurred after species separation. Divergence time between the two species was estimated to be 0.66-1.37 million years ago (Ma), a time at which glaciation was at its maximum in China and deserts developed widely in central Asia. Isolation-with-Migration model also indicated that the two species had diverged in the presence of gene flow. We also detected evidence of selection at GO in P. euphratica and to a lesser extent at PhyB2. Together, these results underscore the importance of Pleistocene climate oscillations in triggering plant speciation as a result of habitats divergence.

New insights on the speciation history and nucleotide diversity of three boreal spruce species and a Tertiary relict

In all, 10 nuclear loci were re-sequenced in four spruce species. Three of the species are boreal species with very large natural ranges: Picea mariana and P. glauca are North American, and P. abies, is Eurasian. The fourth species, P. breweriana, is a Tertiary relict from Northern California, with a very small natural range. Although the boreal species population sizes have fluctuated through the Ice Ages, P. breweriana is believed to have had a rather stable population size through the Quaternary. Indeed, the average Tajima's D was close to zero in this species and negative in the three boreal ones. Reflecting differences in current population sizes, nucleotide diversity was an order of magnitude lower in P. breweriana than in the boreal species. This is in contrast to the similar and high levels of heterozygosity observed in previous studies at allozyme loci across species. As the species have very different histories and effective population sizes, selection at allozyme loci rather than demography appears to be a better explanation for this discrepancy. Parameters of Isolation-with-Migration (IM) models were also estimated for pairs of species. Shared polymorphisms were extensive and fixed polymorphisms few. Divergence times were much shorter than those previously reported. There was also evidence of historical gene flow between P. abies and P. glauca. The latter was more closely related to P. abies than to its sympatric relative P. mariana. This last result suggests that North American and Eurasian species might have been geographically much closer in the recent past than they are today.

Recent degeneration of an old duplicated flowering time gene in Brassica nigra

Gene and genome duplications play a major role in the evolution of plant species. The Brassica nigra genome is highly replicated as a result of ancient polyploidization events. Two copies of the flowering time gene CONSTANS (COa and COb) have been identified in B. nigra, and previous studies showed that COa is functional. In the present study, the polymorphism of 92 COb alleles sampled in seven populations was analyzed. Both polymorphism and recombination levels were elevated and varied strongly among populations and 8% of COb alleles exhibit apparently disabling mutations. Sequence data, however, do not provide unambiguous support for the presence of relaxed selective constraint on COb as compared to known functional CO genes. On the one hand, some of the disabling mutations reached high-frequency arguing for a loss of function but, on the other hand, the ratio of nonsynonymous to synonymous nucleotide polymorphism and diversity is low and similar to that observed in other B. nigra CO and CO-like genes, supporting the conservation of some function. We also showed that COb is still transcribed. Finally, the flowering time of Arabidopsis thaliana co mutant plants transformed with COb alleles with and without apparent disabling mutations was similar. We propose that COb was retained for a long period after duplication, but a recent fixation of a detrimental mutation, possibly as an effect of a bottleneck, resulted in its nonfunctionalization. We also speculate as to the presence of subsequent selection for rapid degeneration of the gene.

Intrageneric phylogeny of Capsella (Brassicaceae) and the origin of the tetraploid C. bursa-pastoris based on chloroplast and nuclear DNA sequences

Polyploidization, often accompanied by hybridization, has been of major importance in flowering plant evolution. Here we investigate the importance of these processes for the evolution of the tetraploid crucifer Capsella bursa-pastoris using DNA sequences from two chloroplast loci as well as from three nuclear low-copy genes. The near-absence of variation at the C. bursa-pastoris chloroplast markers suggests a single and recent origin of the tetraploid. However, despite supporting a single phylogeny, chloroplast data indicate that neither of the extant Capsella diploids is the maternal parent of the tetraploid. Combined with data from the three nuclear loci, our results do not lend support to previous hypotheses on the origin of C. bursa-pastoris as an allopolyploid between the diploids C. grandiflora and C. rubella or an autopolyploid of C. grandiflora. Nevertheless, for each locus, some of the C. bursa-pastoris accessions harbored C. rubella alleles, indicating that C. rubella contributed to the gene pool of C. bursa-pastoris, either through allopolyploid speciation or, more likely, through hybridization and introgression. To our knowledge, this study is the first of a wild, nonmodel plant genus that uses a combination of chloroplast and multiple low-copy nuclear loci for phylogenetic inference of polyploid evolution.

One step beyond lethal equivalents: characterization of deleterious loci in the rapid cycling Brassica rapa L. base population

The total number of lethal equivalents as defined by Morton, Crow and Muller (1956) is a function of three parameters: M, the number of loci at which deleterious mutations can occur, q, the frequency of the deleterious alleles at each locus, and s, their selective value. A new approach based on multi-generation inbreeding data is outlined and used to infer these three parameters as well as the dominance coefficient, h, in a self-incompatible species, Brassica rapa L. Germination and flowering data from thirty bud-selfed lines of fast-cycling B.rapa were assessed over three generations. Germination and flowering were significantly postponed by inbreeding but germination and flowering success were not so strongly decreased. Estimates of the average s values were obtained but it was not possible to get separate estimates of M and q. For both characters, the average dominance coefficient was particularly low. The number of lethal equivalents at the zygotic level was around two for germination and three for flowering, which, owing to the self-incompatibility of B.rapa, is an unexpectedly low value. These results may be explained by past biparental inbreeding which in turn may have increased self-compatibility thus allowing the purging of more deleterious alleles than under strict self-incompatibility.

Ancient DNA from pollen: a genetic record of population history in Scots pine

Assessments of plant population dynamics in space and time have depended on dated records of fossil pollen synthesized on a subcontinental scale. Genetic analyses of extant populations have revealed spatial relationships that are indicative of past spatial dynamics, but lack an explicit timescale. Synthesis of these data requires genetic analyses from abundant dated fossil material, and this has hitherto been lacking. Fossil pollen is the most abundant material with which to fill this data gap. Here we report genetic analyses of fossil pollen retrieved from Holtjärnen postglacial lake sediment in Sweden and show that plastid DNA is recoverable from Scots Pine and Norway spruce pollen grains that are 100 and 10 000 years old. By sequencing clones from two short plastid PCR products and by using multiple controls we show that the ancient sequences were endogenous to the fossil grains. Comparison of ancient sequences and those obtained from an extant population of Scots pine establishes the first genetic link between extant and fossil samples in this species, providing genetic continuity through time. The finding of one common haplotype present in modern, 100-year old and 10 000-year old samples suggests that it may have persisted near Holtjärnen throughout the postglacial period. This retrieval of ancient DNA from pollen has major implications for plant palaeoecology in conifer species by allowing direct estimates of population dynamics in space and time.

On the meaning and existence of an effective population size

We investigate conditions under which a model with stochastic demography or population structure converges to the coalescent with a linear change in timescale. We argue that this is a necessary condition for the existence of a meaningful effective population size. We find that such a linear timescale change is obtained when demographic fluctuations and coalescence events occur on different timescales. Simple models of population structure and randomly fluctuating population size are used to exemplify the ideas and provide an intuitive feel for the meaning of the conditions.

Spatial genetic structure in two congeneric epiphytes with different dispersal strategies analysed by three different methods

Three different approaches were used to assess the kinship structure of two epiphytic bryophytes, Orthotrichum speciosum and O. obtusifolium, that have different dispersal strategies. The two species were sampled in a 200 ha landscape where species occurrence and host trees had been mapped previously. Local environmental conditions at sampled trees were recorded and kinship between individuals was calculated based on amplified fragment length polymorphism (AFLP)-marker data. We did not detect any association between AFLP-markers and investigated environmental conditions. In both species, significant kinship coefficients were found between individuals up to 300-350 m apart which shows that both species have a restricted dispersal range. The spatial kinship structure was detected with both autocorrelation analysis and generalized additive models (GAMs), but linear regression failed to detect any structure in O. speciosum. Although the dioecious O. obtusifolium is currently the more common species it may, none the less, due to its restricted dispersal range and reproduction mode, become threatened in the future by current silvicultural practices which enhance the distance between host trees and decrease their life span. Finally, GAMs seem most appropriate for analysing spatial genetic structure because the effects of local environmental conditions and spatial structure can be analysed simultaneously, no assumption of a parametric form between kinship coefficient and distance is required, and spatial data resolution is not lost in the arbitrary choice of distance classes characterizing autocorrelation analysis.

Extensive sharing of chloroplast haplotypes among European birches indicates hybridization among Betula pendula, B. pubescens and B. nana

Extensive sharing of chloroplast haplotypes among the silver birch, Betula pendula Roth., the downy birch, B. pubescens Ehrh., and the dwarf birch, B. nana L., was discovered using polymerase chain reaction-restriction fragment length polymporphism markers. The geographical component of the genetic variation was stronger than the species component: the species were not significantly different while 11% of the variation could be attributed to differentiation between the two main regions studied, Scandinavia and western Russia. All haplotypes occurring in more than 2% of the individuals were shared among the species and the introgression ratios were quite large: 0.79 between B. pubescens and B. pendula and 0.67 between B. pubescens and B. nana. The data also indicate that B. pendula individuals are more similar to sympatric B. pubescens than to B. pendula individuals from nearby forests. However, this trend is not as pronounced when B. pubescens is considered, suggesting that introgression is not symmetrical. The haplotype sharing among the three Betula species is most likely caused by hybridization and subsequent cytoplasmic introgression.

Absence of geographical structure of chloroplast DNA variation in sallow, Salix caprea L

In the present study, we have used PCR-RFLP markers to investigate the chloroplast DNA variation in 24 European populations of Salix caprea L. A subset of these populations has also been analysed with chloroplast microsatellites. The main feature of both markers is the absence of a clear geographic structure (G(ST(PCR-RFLP))=0.090, G(ST(microsatellites))=-0.017) and high levels of variation within populations. This lack of phylogeographic structure in S. caprea is suggested to be the consequence of the joint action of several factors: (i) presence of intermediate latitude refugia with large population sizes during the last glacial maximum, (ii) a high speed of recolonisation and dispersal ability, (iii) a high mutation rate and (iv) extensive hybridisation with other willow species. In addition to the S. caprea samples, a limited number of individuals from several other Salix species were also analysed with PCR-RFLP: S. cinerea, S. aurita, S. purpurea, S. atrocinerea, S. appendiculata, S. elaeagnos, S. fragilis and S. alba. Many of the haplotypes found in Salix caprea were also detected in S. cinerea, S. aurita, S. purpurea, S. atrocinerea and/ or S. appendiculata but not in S. alba, S. elaeagnos or S. fragilis. Our data suggest that hybridisation and gene flow have occurred within these two groups but not between them.

Genetic mapping of sex-linked markers in Salix viminalis L

A total of 88 selective primer combinations were screened using bulked males and females sampled from four families of Salix viminalis. A total of more than 1000 polymorphic fragments was obtained, of which only four cosegregated with sex. These four sex-linked markers were subsequently scored in individuals that were used for bulked sample preparation in additional individuals of the same families, and in individuals in other families. A pair of primers that amplified the sex-linked fragments was constructed from one of the sex-linked amplified fragment length polymorphism (AFLP) fragments. In hybridization of Southern blot filters with the sex-linked DNA fragments, the band was present in females and absent in males, but the opposite pattern of band segregation (a band found in males and no band in females) was never observed in either the AFLP or RFLP experiments. Two of the sex-linked markers were placed on a linkage map. They both map at the same location in a linkage group comprising other markers not segregating with sex. Our data suggest that a single locus governs the sex determination and that nonrecombining sex chromosomes are absent in S. viminalis. A close association was found between skewed sex ratio and segregation distortion at this locus.

Postglacial recolonization and cpDNA variation of silver birch, Betula pendula

Chloroplast PCR-RFLP markers were used to reconstruct the history of the silver birch, Betula pendula Roth, in Europe since the last glacial maximum (LGM). In birch, fossil pollen maps do not reveal a clear chronological sequence of postglacial spread. If anything, the pollen record suggests that most of Europe was recolonized by birches as early as 10000 bp, probably from populations that remained close to the ice sheets during the LGM. The geographical distribution of haplotypes supports a scenario of early colonization. Two of the 13 haplotypes that were observed are common, representing 35% and 49% of the total sample, respectively. Although one of the common haplotypes is predominant in the NW and the other in the SE, both are present throughout most of the investigated geographical area. Rare haplotypes are geographically restricted. The distribution of the haplotypes reveals five genetic boundaries between groups of haplotypes and allows us to infer patterns of postglacial recolonization. Europe was re-occupied by two main waves of recolonization: one eastern and one western, with origins at intermediate latitudes. Populations in the Iberian Peninsula and in Italy did not take part in the postglacial recolonization of Europe.

The origin of Russian cultivars of red clover ( Trifolium pratense L.) and their genetic relationships to wild populations in the Urals

Propagation and breeding of red clover ( Trifolium pratense L.) in Russia was initiated about 200 years ago but the origin of present-day cultivars is disputed. Some authors argue that most Russian cultivars were derived from western European ones, whereas others support a Russian origin of the cultivars from local wild populations. In the present study we assessed the genetic variation at 17 allozyme loci in seven Russian cultivars, bearing the names of localities of the Urals, two American ones that have been used in Russia for scientific experiments and seven wild populations of the Urals and Western Siberia. Variation at the 17 protein loci supports the western European origin of the cultivars and also indicates that gene flow between cultivars and wild populations was limited or has not acted sufficiently long to affect the genetic composition of the red clover wild populations of the Urals.

Probability of identity by descent in metapopulations

Equilibrium probabilities of identity by descent (IBD), for pairs of genes within individuals, for genes between individuals within subpopulations, and for genes between subpopulations are calculated in metapopulation models with fixed or varying colony sizes. A continuous-time analog to the Moran model was used in either case. For fixed-colony size both propagule and migrant pool models were considered. The varying population size model is based on a birth-death-immigration (BDI) process, to which migration between colonies is added. Wright's F statistics are calculated and compared to previous results. Adding between-island migration to the BDI model can have an important effect on the equilibrium probabilities of IBD and on Wright's index.

A study of deleterious gene structure in plants using Markov chain Monte Carlo

The characteristics of deleterious genes have been of great interest in both theory and practice in genetics. Because of the complex genetic mechanism of these deleterious genes, most current studies try to estimate the overall magnitude of mortality effects on a population, which is characterized classically by the number of lethal equivalents. This number is a combination of several parameters, each of which has a distinct biological effect on genetic mortality. In conservation and breeding programs, it is important to be able to distinguish among different combinations of these parameters that lead to the same number of lethal equivalents, such as a large number of mildly deleterious genes or a few lethal genes, The ability to distinguish such parameter combinations requires more than one generation of mating. We propose a model for survival data from a two-generation mating experiment on the plant species Brassica rapa, and we enable inference with Markov chain Monte Carlo. This computational strategy is effective because a vast amount of missing genotype information must be accounted for. In addition to the lethal equivalents, the two-generation data provide separate information on the average intensity of mortality and the average number of deleterious genes carried by an individual. In our Markov chain Monte Carlo algorithm, we use a vector proposal distribution to overcome inefficiency of a single-site Gibbs sampler. Information about environmental effects is obtained from an outcrossing experiment conducted in parallel with the two-generation mating experiments.

Components of resistance of leaf rust (Melamspora laricii epitea Kleb./Melampsora ribesii-viminalis Kleb.) in Salix viminalis L

Selection on partial resistance components, namely latent period, development rate of uredia and number and size of uredia, has been suggested as a means to achieve durable resistance. Three experiments were carried out in growth chambers to assess the impact of environmental and genetic factors on these components in the Salix-Melampsora host-pathogen system. They confirmed the environmental ability of the Melampsora-Salicaceae relationship and provided no definite answer on the possibility of attaining durable resistance through selection on partial resistance components. On the one hand, there is a large amount of heritable variation among clones for most components; on the other hand, all components were extremely sensitive to environmental conditions and isolates, making the development of standard selection methods difficult.

Number of lethal equivalents in human populations: how good are the previous estimates?

Lee et al. (1996) recently developed a method for interval estimation of the number of lethal equivalents by using a hierarchical structure of likelihood functions. This hierarchical model consists of two multinomial trials: one of the sampling process of the parents from the population of interest, and the other for the survival of the offspring of the families obtained by mating the parents. The method, initially developed for selfing and full-sib mating, is extended here to include more general mating systems as well as mixtures of mating systems. We applied it to human data sets for which confidence intervals were previously not available. Our point estimates were close to previous ones, and the standard deviations were generally quite small. Thus, even if debate over the meaning of the concept of lethal equivalents has not been entirely resolved, our results showed that the previous estimates are at least statistically meaningful.

Impact of mating design on selection response in Brassica rapa L

The impact of four mating designs on selection response for leaf area was assessed at four different population sizes, using fast-cycling Brassica rapa L. Mating designs were either balanced (partial diallel or pair mating) or unbalanced (factorial mating designs with either one or two testers). When balanced, the mating designs required different numbers of crossings for the same number of parents: the partial diallel design, in the configuration retained here, required three times as many crossings as pair mating. Population sizes were 4, 8, 16, and 32. The percentage of selected individuals was kept constant at 25%. Despite an average estimated heritability around 0.4, the overall response to selection after five generations was fairly weak in all three replicates. For a given population size, selection response was larger under balanced mating designs than under unbalanced ones. There was no difference among balanced mating designs. Both results indicate that effective population size is more important than population size or the number of crossings in maintaining genetic gain.