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.

Impact of sampling schemes on demographic inference: an empirical study in two species with different mating systems and demographic histories.

Most species have at least some level of genetic structure. Recent simulation studies have shown that it is important to consider population structure when sampling individuals to infer past population history. The relevance of the results of these computer simulations for empirical studies, however, remains unclear. In the present study, we use DNA sequence datasets collected from two closely related species with very different histories, the selfing species Capsella rubella and its outcrossing relative C. grandiflora, to assess the impact of different sampling strategies on summary statistics and the inference of historical demography. Sampling strategy did not strongly influence the mean values of Tajima's D in either species, but it had some impact on the variance. The general conclusions about demographic history were comparable across sampling schemes even when resampled data were analyzed with approximate Bayesian computation (ABC). We used simulations to explore the effects of sampling scheme under different demographic models. We conclude that when sequences from modest numbers of loci (<60) are analyzed, the sampling strategy is generally of limited importance. The same is true under intermediate or high levels of gene flow (4Nm > 2-10) in models in which global expansion is combined with either local expansion or hierarchical population structure. Although we observe a less severe effect of sampling than predicted under some earlier simulation models, our results should not be seen as an encouragement to neglect this issue. In general, a good coverage of the natural range, both within and between populations, will be needed to obtain a reliable reconstruction of a species's demographic history, and in fact, the effect of sampling scheme on polymorphism patterns may itself provide important information about demographic history.

Permanent Genetic Resources added to Molecular Ecology Resources Database 1 October 2009-30 November 2009.

This article documents the addition of 411 microsatellite marker loci and 15 pairs of Single Nucleotide Polymorphism (SNP) sequencing primers to the Molecular Ecology Resources Database. Loci were developed for the following species: Acanthopagrus schlegeli, Anopheles lesteri, Aspergillus clavatus, Aspergillus flavus, Aspergillus fumigatus, Aspergillus oryzae, Aspergillus terreus, Branchiostoma japonicum, Branchiostoma belcheri, Colias behrii, Coryphopterus personatus, Cynogolssus semilaevis, Cynoglossus semilaevis, Dendrobium officinale, Dendrobium officinale, Dysoxylum malabaricum, Metrioptera roeselii, Myrmeciza exsul, Ochotona thibetana, Neosartorya fischeri, Nothofagus pumilio, Onychodactylus fischeri, Phoenicopterus roseus, Salvia officinalis L., Scylla paramamosain, Silene latifo, Sula sula, and Vulpes vulpes. These loci were cross-tested on the following species: Aspergillus giganteus, Colias pelidne, Colias interior, Colias meadii, Colias eurytheme, Coryphopterus lipernes, Coryphopterus glaucofrenum, Coryphopterus eidolon, Gnatholepis thompsoni, Elacatinus evelynae, Dendrobium loddigesii Dendrobium devonianum, Dysoxylum binectariferum, Nothofagus antarctica, Nothofagus dombeyii, Nothofagus nervosa, Nothofagus obliqua, Sula nebouxii, and Sula variegata. This article also documents the addition of 39 sequencing primer pairs and 15 allele specific primers or probes for Paralithodes camtschaticus.

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.

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.

Genetic relationship among eurasian and american larix species based on allozymes

Genetic variation at 16 allozyme loci was studied in both American (Larix occidentalis Nutt., L. laricina (Du Roi) C. Koch, L. lyallii Parl.) and Eurasian (L. sibirica Ledeb., L. gmelinii Rupr., L. olgensis A. Henry, L. kaempferi (Lamb.) Carr. (=L. leptolepis (Sieb. et Zucc.) Endl.), L. kamtschatica (Rupr.) Carr. and L. decidua (Mill. )) larch species. Species with a limited range, such as L. olgensis and L. lyallii, had lower genetic variation than species with a wider range. Population differentiation within species was of the same order of magnitude among species. The resulting phylogeny indicates a clear separation between American and Eurasian species. This result is in agreement with recent palaeontological findings that suggest that gene flow between American and Eurasian species has been unlikely since the last glaciation.

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.

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.

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.