Naoki Takebayashi

The frequency of polyploid speciation in vascular plants

Since its discovery in 1907, polyploidy has been recognized as an important phenomenon in vascular plants, and several lines of evidence indicate that most, if not all, plant species ultimately have a polyploid ancestry. However, previous estimates of the frequency of polyploid speciation suggest that the formation and establishment of neopolyploid species is rare. By combining information from the botanical communitys vast cytogenetic and phylogenetic databases, we establish that 15% of angiosperm and 31% of fern speciation events are accompanied by ploidy increase. These frequency estimates are higher by a factor of four than earlier estimates and lead to a standing incidence of polyploid species within genera of 35% (n = 1,506). Despite this high incidence, we find no direct evidence that polyploid lines, once established, enjoy greater net species diversification. Thus, the widespread occurrence of polyploid taxa appears to result from the substantial contribution of polyploidy to cladogenesis, but not from subsequent increases in diversification rates of polyploid lines.

MTML-msBayes: Approximate Bayesian comparative phylogeographic inference from multiple taxa and multiple loci with rate heterogeneity

BackgroundMTML-msBayes uses hierarchical approximate Bayesian computation (HABC) under a coalescent model to infer temporal patterns of divergence and gene flow across codistributed taxon-pairs. Under a model of multiple codistributed taxa that diverge into taxon-pairs with subsequent gene flow or isolation, one can estimate hyper-parameters that quantify the mean and variability in divergence times or test models of migration and isolation. The software uses multi-locus DNA sequence data collected from multiple taxon-pairs and allows variation across taxa in demographic parameters as well as heterogeneity in DNA mutation rates across loci. The method also allows a flexible sampling scheme: different numbers of loci of varying length can be sampled from different taxon-pairs.ResultsSimulation tests reveal increasing power with increasing numbers of loci when attempting to distinguish temporal congruence from incongruence in divergence times across taxon-pairs. These results are robust to DNA mutation rate heterogeneity. Estimating mean divergence times and testing simultaneous divergence was less accurate with migration, but improved if one specified the correct migration model. Simulation validation tests demonstrated that one can detect the correct migration or isolation model with high probability, and that this HABC model testing procedure was greatly improved by incorporating a summary statistic originally developed for this task (Wakeleys ΨW ). The method is applied to an empirical data set of three Australian avian taxon-pairs and a result of simultaneous divergence with some subsequent gene flow is inferred.ConclusionsTo retain flexibility and compatibility with existing bioinformatics tools, MTML-msBayes is a pipeline software package consisting of Perl, C and R programs that are executed via the command line. Source code and binaries are available for download at http://msbayes.sourceforge.net/ under an open source license (GNU Public License).

Pollen Limitation and the Evolution of Androdioecy from Dioecy

Androdioecy is an unusual breeding system in which populations consist of separate male and hermaphrodite individuals. The evolution of androdioecy is still poorly understood; however, there is evidence from several androdioecious species that the breeding system may have evolved from dioecy (males and females). This article presents a simple deterministic model showing that androdioecy can evolve from dioecy under a broad range of realistic conditions. For the evolution of androdioecy from dioecy, hermaphrodites must be able to invade the dioecious population. Then, males must be maintained, while females are eliminated. Hermaphrodite invasion is favored when females are pollen limited and hermaphrodites have high overall fertility and are self‐fertile. Male maintenance is favored when hermaphrodites resemble females, having high seed production and low pollen fitness, and when the selfing rate is not too high. These conditions were satisfied over a broad and realistic range of parameter values, suggesting that the evolution of androdioecy from dioecy is highly plausible.

Nucleotide diversity and linkage disequilibrium in balsam poplar (Populus balsamifera).

*Current perceptions that poplars have high levels of nucleotide variation, large effective population sizes, and rapid decay of linkage disequilibrium are based primarily on studies from one poplar species, Populus tremula. *We analysed 590 gene fragments (average length 565 bp) from each of 15 individuals from different populations from throughout the range of Populus balsamifera. *Nucleotide diversity (theta(total) = 0.0028, pi = 0.0027) was low compared with other trees and model agricultural systems. Patterns of nucleotide diversity and site frequency spectra were consistent with purifying selection on replacement and intron sites. When averaged across all loci we found no evidence for decay of linkage disequilibrium across 750 bp, consistent with the low estimates of the scaled recombination parameter, rho = 0.0092. *Compared with P. tremula, a well studied congener with a similar distribution, P. balsamifera has low diversity and low effective recombination, both of which indicate a lower effective population size in P. balsamifera. Patterns of diversity and linkage indicate that there is considerable variation in population genomic patterns among poplar species and unlike P. tremula, association mapping techniques in balsam poplar should consider sampling single nucleotide polymorphisms (SNPs) at well-spaced intervals.

Possible Diversifying Selection in the Imprinted Gene, MEDEA, in Arabidopsis

Coevolutionary conflict among imprinted genes that influence traits such as offspring growth may arise when maternal and paternal genomes have different evolutionary optima. This conflict is expected in outcrossing taxa with multiple paternity, but not self-fertilizing taxa. MEDEA (MEA) is an imprinted plant gene that influences seed growth. Disagreement exists regarding the type of selection acting on this gene. We present new data and analyses of sequence diversity of MEA in self-fertilizing and outcrossing Arabidopsis and its relatives, to help clarify the form of selection acting on this gene. Codon-based branch analysis among taxa (PAML) suggests that selection on the coding region is changing over time, and nonsynonymous substitution is elevated in at least one outcrossing branch. Codon-based analysis of diversity within outcrossing Arabidopsis lyrata ssp. petraea (OmegaMap) suggests that diversifying selection is acting on a portion of the gene, to cause elevated nonsynonymous polymorphism. Providing further support for balancing selection in A. lyrata, Hudson, Kreitman and Aguadé analysis indicates that diversity/divergence at silent sites in the MEA promoter and genic region is elevated relative to reference genes, and there are deviations from the neutral frequency spectrum. This combination of positive selection as well as balancing and diversifying selection in outcrossing lineages is consistent with other genes influence by evolutionary conflict, such as disease resistance genes. Consistent with predictions that conflict would be eliminated in self-fertilizing taxa, we found no evidence of positive, balancing, or diversifying selection in A. thaliana promoter or genic region.

CALIBRATING A MOLECULAR CLOCK FROM PHYLOGEOGRAPHIC DATA: MOMENTS AND LIKELIHOOD ESTIMATORS

Abstract We present moments and likelihood methods that estimate a DNA substitution rate from a group of closely related sister species pairs separated at an assumed time, and we test these methods with simulations. The methods also estimate ancestral population size and can test whether there is a significant difference among the ancestral population sizes of the sister species pairs. Estimates presented in the literature often ignore the ancestral coalescent prior to speciation and therefore should be biased upward. The simulations show that both methods yield accurate estimates given sample sizes of five or more species pairs and that better likelihood estimates are obtained if there is no significant difference among ancestral population sizes. The model presented here indicates that the larger than expected variation found in multitaxa datasets can be explained by variation in the ancestral coalescence and the Poisson mutation process. In this context, observed variation can often be accounted for by variation in ancestral population sizes rather than invoking variation in other parameters, such as divergence time or mutation rate. The methods are applied to data from two groups of species pairs (sea urchins and Alpheus snapping shrimp) that are thought to have separated by the rise of Panama three million years ago.

RECOMMENDATIONS FOR USING MSBAYES TO INCORPORATE UNCERTAINTY IN SELECTING AN ABC MODEL PRIOR: A RESPONSE TO OAKS ET AL.

Prior specification is an essential component of parameter estimation and model comparison in Approximate Bayesian computation (ABC). Oaks et al. present a simulation‐based power analysis of msBayes and conclude that msBayes has low power to detect genuinely random divergence times across taxa, and suggest the cause is Lindleys paradox. Although the predictions are similar, we show that their findings are more fundamentally explained by insufficient prior sampling that arises with poorly chosen wide priors that critically undersample nonsimultaneous divergence histories of high likelihood. In a reanalysis of their data on Philippine Island vertebrates, we show how this problem can be circumvented by expanding upon a previously developed procedure that accommodates uncertainty in prior selection using Bayesian model averaging. When these procedures are used, msBayes supports recent divergences without support for synchronous divergence in the Oaks et al. data and we further present a simulation analysis that demonstrates that msBayes can have high power to detect asynchronous divergence under narrower priors for divergence time. Our findings highlight the need for exploration of plausible parameter space and prior sampling efficiency for ABC samplers in high dimensions. We discus potential improvements to msBayes and conclude that when used appropriately with model averaging, msBayes remains an effective and powerful tool.

Cold tolerance in Arabidopsis kamchatica

UNLABELLED • PREMISE OF THE STUDY Cold tolerance is a critically important factor determining how plants will be influenced by climate change, including changes in snowcover and extreme weather events. Although a great deal is known about cold tolerance in Arabidopsis thaliana, it is not highly cold tolerant. This study examined cold tolerance and its genetic diversity in an herbaceous subarctic relative, Arabidopsis kamchatica, which generally occurs in much colder climates.• METHODS Thermal analysis and electrolyte leakage were used to estimate supercooling points and lethal temperatures (LT50) in cold-acclimated and nonacclimated families from three populations of A. kamchatica.• KEY RESULTS Arabidopsis kamchatica was highly cold tolerant, with a mean LT50 of -10.8°C when actively growing, and -21.8°C when cold acclimated. It also was able to supercool to very low temperatures. Surprisingly, actively growing plants supercooled more than acclimated plants (-14.7 vs. -12.7°C). There was significant genetic variation for cold tolerance both within and among populations. However, both cold tolerance and genetic diversity were highest in the midlatitude population rather than in the far north, indicating that adaptations to climate change are most likely to arise in the center of the species range rather than at the edges.• CONCLUSIONS Arabidopsis kamchatica is highly cold tolerant throughout its range. It is far more freeze tolerant than A. thaliana, and supercooled to lower temperatures, suggesting that A. kamchatica provides a valuable complement to A. thaliana for cold tolerance research.

Genome size variation and evolution in allotetraploid Arabidopsis kamchatica and its parents, Arabidopsis lyrata and Arabidopsis halleri

Genome doubling and changes in genome size are fundamental evolutionary processes. Arabidopsis kamchatica has been reported to contain both diploid and tetraploid individuals (2 or 4 copies of every chromosome). We did find genome size differences among populations, and among populations genome size varied 7%. However, all sampled A. kamchatica plants from a wide geographic range were tetraploids. This level of intraspecific genome size variation in A. kamchatica is lower than in other Arabidopsis taxa. Due to its close relationship to A. thaliana, A. kamchatica has the potential to be very useful in the study of polyploidy and genome evolution.

Evolutionary Genetics of an S-Like Polymorphism in Papaveraceae with Putative Function in Self-Incompatibility

Background Papaver rhoeas possesses a gametophytic self-incompatibility (SI) system not homologous to any other SI mechanism characterized at the molecular level. Four previously published full length stigmatic S-alleles from the genus Papaver exhibited remarkable sequence divergence, but these studies failed to amplify additional S-alleles despite crossing evidence for more than 60 S-alleles in Papaver rhoeas alone. Methodology/Principal Findings Using RT-PCR we identified 87 unique putative stigmatic S-allele sequences from the Papaveraceae Argemone munita, Papaver mcconnellii, P. nudicuale, Platystemon californicus and Romneya coulteri. Hand pollinations among two full-sib families of both A. munita and P. californicus indicate a strong correlation between the putative S-genotype and observed incompatibility phenotype. However, we also found more than two S-like sequences in some individuals of A. munita and P. californicus, with two products co-segregating in both full-sib families of P. californicus. Pairwise sequence divergence estimates within and among taxa show Papaver stigmatic S-alleles to be the most variable with lower divergence among putative S-alleles from other Papaveraceae. Genealogical analysis indicates little shared ancestral polymorphism among S-like sequences from different genera. Lack of shared ancestral polymorphism could be due to long divergence times among genera studied, reduced levels of balancing selection if some or all S-like sequences do not function in incompatibility, population bottlenecks, or different levels of recombination among taxa. Preliminary estimates of positive selection find many sites under selective constraint with a few undergoing positive selection, suggesting that self-recognition may depend on amino acid substitutions at only a few sites. Conclusions/Significance Because of the strong correlation between genotype and SI phenotype, sequences reported here represent either functional stylar S-alleles, tightly linked paralogs of the S-locus or a combination of both. The considerable complexity revealed in this study shows we have much to learn about the evolutionary dynamics of self-incompatibility systems.