Allele frequency spectra in structured populations: Novel-allele probabilities under the labelled coalescent.

We address the effect of population structure on key properties of the Ewens sampling formula. We use our previously-introduced inductive method for determining exact allele frequency spectrum (AFS) probabilities under the infinite-allele model of mutation and population structure for samples of arbitrary size. Fundamental to the sampling distribution is the novel-allele probability, the probability that given the pattern of variation in the present sample, the next gene sampled belongs to an as-yet-unobserved allelic class. Unlike the case for panmictic populations, the novel-allele probability depends on the AFS of the present sample. We derive a recursion that directly provides the marginal novel-allele probability across AFSs, obviating the need first to determine the probability of each AFS. Our explorations suggest that the marginal novel-allele probability tends to be greater for initial samples comprising fewer alleles and for sampling configurations in which the next-observed gene derives from a deme different from that of the majority of the present sample. Comparison to the efficient importance sampling proposals developed by De Iorio and Griffiths and colleagues indicates that their approximation for the novel-allele probability generally agrees with the true marginal, although it may tend to overestimate the marginal in cases in which the novel-allele probability is high and migration rates are low.

Cold tolerance in the genus Arabidopsis.

PREMISE: Cold tolerance is an important factor limiting the geographic distribution and growing season for many plant species, yet few studies have examined variation in cold tolerance extensively within and among closely related species and compared that to their geographic distribution. METHODS: This study examines cold tolerance within and among species in the genus Arabidopsis. We assessed cold tolerance by measuring electrolyte leakage from detached leaves in multiple populations of five Arabidopsis taxa. The temperature at which 50% of cells were lysed was considered the lethal temperature (LT50 ). RESULTS: We found variability within and among taxa in cold tolerance. There was no significant within-species relationship between latitude and cold tolerance. However, the northern taxa, A. kamchatica, A. lyrata subsp. petraea, and A. lyrata subsp. lyrata, were more cold tolerant than A. thaliana and A. halleri subsp. gemmifera both before and after cold acclimation. Cold tolerance increased after cold acclimation (exposure to low, but nonfreezing temperatures) for all taxa, although the difference was not significant for A. halleri subsp. gemmifera. For all taxa except A. lyrata subsp. lyrata, the LT50 values for cold-acclimated plants were higher than the January mean daily minimum temperature (Tmin ), indicating that if plants were not insulated by snow cover, they would not likely survive winter at the northern edge of their range. CONCLUSIONS: Arabidopsis lyrata and A. kamchatica were far more cold tolerant than A. thaliana. These extremely cold-tolerant taxa are excellent candidates for studying both the molecular and ecological aspects of cold tolerance.

Inductive determination of allele frequency spectrum probabilities in structured populations.

We present a method for inductively determining exact allele frequency spectrum (AFS) probabilities for samples derived from a population comprising two demes under the infinite-allele model of mutation. This method builds on a labeled coalescent argument to extend the Ewens sampling formula (ESF) to structured populations. A key departure from the panmictic case is that the AFS conditioned on the number of alleles in the sample is no longer independent of the scaled mutation rate (θ). In particular, biallelic site frequency spectra, widely-used in explorations of genome-wide patterns of variation, depend on the mutation rate in structured populations. Variation in the rate of substitution across loci and through time may contribute to apparent distortions of site frequency spectra exhibited by samples derived from structured populations.

Evolution of the sex ratio and effective number under gynodioecy and androdioecy.

We address the evolution of effective number of individuals under androdioecy and gynodioecy. We analyze dynamic models of autosomal modifiers of weak effect on sex expression. In our zygote control models, the sex expressed by a zygote depends on its own genotype, while in our maternal control models, it depends on the genotype of its maternal parent. Our analysis unifies full multi-dimensional local stability analysis with the Li-Price equation, which for all its heuristic appeal, describes evolutionary change over a single generation. We define a point in the neighborhood of a fixation state from which a single-generation step indicates the asymptotic behavior of the frequency of a modifier allele initiated at an arbitrary point near the fixation state. A concept of heritability appropriate for the evolutionary modification of sex emerges from the Li-Priceframework. We incorporate our theoretical analysis into our previously-developed Bayesian inference framework to develop a new method for inferring the viability of gonochores (males or females) relative to hermaphrodites. Applying this approach to microsatellite data derived from natural populations of the gynodioecious plant Schiedea salicaria and the androdioecious killifish Kryptolebias marmoratus, we find that while female and hermaphrodite S. salicaria appear to have similar viabilities, male K. marmoratus appear to survive to reproductive age at less than half the rate of hermaphrodites.

Beringian sub-refugia revealed in blackfish (Dallia): implications for understanding the effects of Pleistocene glaciations on Beringian taxa and other Arctic aquatic fauna.

BACKGROUND: Pleistocene climatic instability had profound and diverse effects on the distribution and abundance of Arctic organisms revealed by variation in phylogeographic patterns documented in extant Arctic populations. To better understand the effects of geography and paleoclimate on Beringian freshwater fishes, we examined genetic variability in the genus Dallia (blackfish: Esociformes: Esocidae). The genus Dallia groups between one and three nominal species of small, cold- and hypoxia-tolerant freshwater fishes restricted entirely in distribution to Beringia from the Yukon River basin near Fairbanks, Alaska westward including the Kuskokwim River basin and low-lying areas of Western Alaska to the Amguema River on the north side of the Chukotka Peninsula and Mechigmen Bay on the south side of the Chukotka Peninsula. The genus has a non-continuous distribution divided by the Bering Strait and the Brooks Range. We examined the distribution of genetic variation across this range to determine the number and location of potential sub-refugia within the greater Beringian refugium as well as the roles of the Bering land bridge, Brooks Range, and large rivers within Beringia in shaping the current distribution of populations of Dallia. Our analyses were based on DNA sequence data from two nuclear gene introns (S7 and RAG1) and two mitochondrial genome fragments from nineteen sampling locations. These data were examined under genetic clustering and coalescent frameworks to identify sub-refugia within the greater Beringia refugium and to infer the demographic history of different populations of Dallia. RESULTS: We identified up to five distinct genetic clusters of Dallia. Four of these genetic clusters are present in Alaska: (1) Arctic Coastal Plain genetic cluster found north of the Brooks Range, (2) interior Alaska genetic cluster placed in upstream locations in the Kuskokwim and Yukon river basins, (3) a genetic cluster found on the Seward Peninsula, and (4) a coastal Alaska genetic cluster encompassing downstream Kuskokwim River and Yukon River basin sample locations and samples from Southwest Alaska not in either of these drainages. The Chukotka samples are assigned to their own genetic cluster (5) similar to the coastal Alaska genetic cluster. The clustering and ordination analyses implemented in Discriminant Analysis of Principal Components (DAPC) and STRUCTURE showed mostly concordant groupings and a high degree of differentiation among groups. The groups of sampling locations identified as genetic clusters correspond to geographic areas divided by likely biogeographic barriers including the Brooks Range and the Bering Strait. Estimates of sequence diversity (θ) are highest in the Yukon River and Kuskokwim River drainages near the Bering Sea. We also infer asymmetric migration rates between genetic clusters. The isolation of Dallia on the Arctic Coastal Plain of Alaska is associated with very low estimated migration rates between the coastal Alaska genetic cluster and the Arctic Coastal Plain genetic cluster. CONCLUSIONS: Our results support a scenario with multiple aquatic sub-refugia in Beringia during the Pleistocene and the preservation of that structure in extant populations of Dallia. An inferred historical presence of Dallia across the Bering land bridge explains the similarities in the genetic composition of Dallia in West Beringia and western coastal Alaska. In contrast, historic and contemporary isolation across the Brooks Range shaped the distinctiveness of present day Arctic Coastal Plain Dallia. Overall this study uncovered a high degree of genetic structuring among populations of Dallia supporting the idea of multiple Beringian sub-refugia during the Pleistocene and which appears to be maintained to the present due to the strictly freshwater nature and low dispersal ability of this genus.

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.

Polyploidization and subsequent changes in genome size are fundamental processes in evolution and diversification. Little is currently known about the extent of genome size variation within taxa and the evolutionary forces acting on this variation. Arabidopsis kamchatica has been reported to contain both diploid and tetraploid individuals. The aim of this study was to determine the genome size of A. kamchatica, whether there is variation in ploidy and/or genome size in A. kamchatica and to study how genome size has evolved. We used propidium iodide flow cytometry to measure 2C DNA content of 73 plants from 25 geographically diverse populations of the putative allotetraploid A. kamchatica and its parents, Arabidopsis lyrata and Arabidopsis halleri. All A. kamchatica plants appear to be tetraploids. The mean 2C DNA content of A. kamchatica was 1.034 pg (1011 Mbp), which is slightly smaller than the sum of its diploid parents (A. lyrata: 0.502 pg; A. halleri: 0.571 pg). Arabidopsis kamchatica appears to have lost ∼37.594 Mbp (3.6 %) of DNA from its 2C genome. Tetraploid A. lyrata from Germany and Austria appears to have lost ∼70.366 Mbp (7.2 %) of DNA from the 2C genome, possibly due to hybridization with A. arenosa, which has a smaller genome than A. lyrata. We did find genome size differences among A. kamchatica populations, which varied up to 7 %. Arabidopsis kamchatica ssp. kawasakiana from Japan appears to have a slightly larger genome than A. kamchatica ssp. kamchatica from North America, perhaps due to multiple allopolyploid origins or hybridization with A. halleri. However, the among-population coefficient of variation in 2C DNA content is lower in A. kamchatica 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.

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 Lindley's 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.

Evaluation of the authenticity of a highly novel environmental sequence from boreal forest soil using ribosomal RNA secondary structure modeling.

The number of sequences from both formally described taxa and uncultured environmental DNA deposited in the International Nucleotide Sequence Databases has increased substantially over the last two decades. Although the majority of these sequences represent authentic gene copies, there is evidence of DNA artifacts in these databases as well. These include lab artifacts, such as PCR chimeras, and biological artifacts such as pseudogenes or other paralogous sequences. Sequences that fall in basal positions in phylogenetic trees and appear distant from known sequences are particularly suspect. Phylogenetic analyses suggest that a novel sequence type (NS1) found in two boreal forest soil clone libraries belongs to the fungal kingdom but does not fall unambiguously within any known phylum. We have evaluated this sequence type using an array of secondary-structure analyses. To our knowledge, such analyses have never been used on environmental ribosomal sequences. Ribosomal secondary structure was modeled for four rRNA loci (ITS1, 5.8S, ITS2, 5' LSU). These models were analyzed for the presence of conserved domains, conserved nucleotide motifs, and compensatory base changes. Minimal free energy (MFE) foldings and GC contents of sequences representing the major fungal clades, as well as NS1, were also compared. NS1 displays secondary rRNA structures consistent with other fungi and many, but not all, conserved nucleotide motifs found across eukaryotes. However, our analyses show that many other authentic sequences from basal fungi lack more of these conserved motifs than does NS1. Together our findings suggest that NS1 represents an authentic gene copy. The methods described here can be used on any rRNA-coding sequence, not just environmental fungal sequences. As new-generation sequencing methods that yield shorter sequences become more widely implemented, methods that evaluate sequence authenticity should also be more widely implemented. For fungi, the adjacent 5.8S and ITS2 loci should be prioritized. This region is not only suited to distinguishing between closely related species, but it is also more informative in terms of expected secondary structure.

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.

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

BACKGROUND: MTML-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. RESULTS: Simulation 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 (Wakeley's Ψ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. CONCLUSIONS: To 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).

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.

Heterogeneous selection on trichome production in Alaskan Arabidopsis kamchatica (Brassicaceae).

PREMISE OF THE STUDY: Environmental heterogeneity is thought to be one of the primary factors in the evolutionary maintenance of morphological variation. Here, we explore the role of environmental heterogeneity in the maintenance of variation in leaf hair (trichome) production in Arabidopsis kamchatica. • METHODS: We investigate abiotic correlates of trichome production in A. kamchatica via surveys of both herbarium specimens and wild populations. In addition, we examine patterns of phenotypic selection on trichome production among populations that differ in environmental characteristics. • KEY RESULTS: Trichome-producing herbarium specimens were more likely to occur at lower latitudes and in locations with lower mean annual precipitation and less annual variation in temperature than glabrous specimens. In surveys of wild populations, frequencies of trichome-producing plants were higher in drier habitats than in wetter environments. Using phenotypic selection analysis, we found divergent selection through female fitness (fruit production) on trichome number in populations that differ in environmental characteristics; there was selection for reduced trichome number in one population and selection for increased trichome number in another population. In a population containing both glabrous and trichome-producing plants, glabrous plants produced significantly more fruits than trichome-producing individuals, which indicates selection against the trichome morph. • CONCLUSIONS: Our results demonstrate that there is heterogeneity in selection among populations, which could be responsible for the maintenance of trichome variation in Alaskan populations of A. kamchatica.

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 community's 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.

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.

msBayes: pipeline for testing comparative phylogeographic histories using hierarchical approximate Bayesian computation.

BACKGROUND: Although testing for simultaneous divergence (vicariance) across different population-pairs that span the same barrier to gene flow is of central importance to evolutionary biology, researchers often equate the gene tree and population/species tree thereby ignoring stochastic coalescent variance in their conclusions of temporal incongruence. In contrast to other available phylogeographic software packages, msBayes is the only one that analyses data from multiple species/population pairs under a hierarchical model. RESULTS: msBayes employs approximate Bayesian computation (ABC) under a hierarchical coalescent model to test for simultaneous divergence (TSD) in multiple co-distributed population-pairs. Simultaneous isolation is tested by estimating three hyper-parameters that characterize the degree of variability in divergence times across co-distributed population pairs while allowing for variation in various within population-pair demographic parameters (sub-parameters) that can affect the coalescent. msBayes is a software package consisting of several C and R programs that are run with a Perl "front-end". CONCLUSION: The method reasonably distinguishes simultaneous isolation from temporal incongruence in the divergence of co-distributed population pairs, even with sparse sampling of individuals. Because the estimate step is decoupled from the simulation step, one can rapidly evaluate different ABC acceptance/rejection conditions and the choice of summary statistics. Given the complex and idiosyncratic nature of testing multi-species biogeographic hypotheses, we envision msBayes as a powerful and flexible tool for tackling a wide array of difficult research questions that use population genetic data from multiple co-distributed species. The msBayes pipeline is available for download at http://msbayes.sourceforge.net/ under an open source license (GNU Public License). The msBayes pipeline is comprised of several C and R programs that are run with a Perl "front-end" and runs on Linux, Mac OS-X, and most POSIX systems. Although the current implementation is for a single locus per species-pair, future implementations will allow analysis of multi-loci data per species pair.

Evaluation of psychophysiological asymmetry in patients with fibromyalgia syndrome.

Fibromyalgia syndrome (FMS) is characterized by systemic pain of unknown etiology, and is often accompanied by various psychological symptoms. In the present study, differences in surface electromyographic (SEMG) levels of the trapezius muscle, skin temperature (TEMP) and skin conductance level (SCL) were compared between the right and left side of the body in 31 FMS and 47 control subjects (Control Group). We observed significant asymmetries of SEMG level, TEMP and SCL in the FMS Group. These asymmetries might be related to central, peripheral and autonomic nervous system dysfunctions. Marked increase of SEMG levels, and a decrease of TEMP and SCL were observed at the dominant side in the FMS Group, and a negative correlation of SEMG levels with TEMP and SCL was found. These results suggest that continued antalgic postures in response to pain at the dominant side in FMS patients might lead to asymmetries of SEMG level, TEMP and SCL. Thus, a focus on pain related behaviors and muscle asymmetry might be a useful therapeutic approach.

Maximum-likelihood estimation of rates of recombination within mating-type regions.

Features common to many mating-type regions include recombination suppression over large genomic tracts and cosegregation of genes of various functions, not necessarily related to reproduction. Model systems for homomorphic self-incompatibility (SI) in flowering plants share these characteristics. We introduce a method for the exact computation of the joint probability of numbers of neutral mutations segregating at the determinant of mating type and at a linked marker locus. The underlying Markov model incorporates strong balancing selection into a two-locus coalescent. We apply the method to obtain a maximum-likelihood estimate of the rate of recombination between a marker locus, 48A, and S-RNase, the determinant of SI specificity in pistils of Nicotiana alata. Even though the sampled haplotypes show complete allelic linkage disequilibrium and recombinants have never been detected, a highly significant deficiency of synonymous substitutions at 48A compared to S-RNase suggests a history of recombination. Our maximum-likelihood estimate indicates a rate of recombination of perhaps 3 orders of magnitude greater than the rate of synonymous mutation. This approach may facilitate the construction of genetic maps of regions tightly linked to targets of strong balancing selection.

A simple method for computing exact probabilities of mutation numbers.

We describe a method for the recursive computation of exact probability distributions for the number of neutral mutations segregating in samples of arbitrary size and configuration. Construction of the recursions requires only characterization of evolutionary changes as a Markov process and determination of one-step transition matrices. We address the pattern of nucleotide diversity at a neutral marker locus linked to a determinant of mating type. Under a reformulation of parameters, the method also applies directly to metapopulation models with island migration among demes. Characterization of complete probability distributions facilitates parameter estimation and hypothesis testing by likelihood- as well as moment-based approaches.

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.