Establishment process of a magic trait allele subject to both divergent selection and assortative mating.

Sexual selection and divergent selection are among the major driving forces of reproductive isolation, which could eventually result in speciation. A magic trait is defined such that a single trait is subject to both divergent selection and mate choice through phenotype-based assortative mating. We are here interested in the evolutionary behavior of alleles at a genetic locus responsible for a magic trait in a finite population. We assume that, in a pair of homogeneous subpopulations, a mutant allele arises at the magic trait locus, and theoretically obtain the probability that the new allele establishes in the population, or the establishment probability. We also show an analytical expression for the trajectory of allele frequency along the establishment, from which the time required for the establishment is obtained, or the establishment time. Under this model, divergent selection simply favors the new allele to fix where it is beneficial, whereas assortative mating works against rare alleles. It is theoretically demonstrated that the fate of the new allele is determined by the relative contributions of the two selective forces, divergent selection and assortative mating, when the allele is rare so that the two selective forces counteract. Our theoretical results for the establishment probability and time allow us to understand the relative role of random genetic drift in the establishment process of a magic trait allele in a finite population.

Divergent selection for opsin gene variation in guppy (Poecilia reticulata) populations of Trinidad and Tobago.

The guppy is known to exhibit remarkable interindividual variations in spectral sensitivity of middle to long wavelength-sensitive (M/LWS) cone photoreceptor cells. The guppy has four M/LWS-type opsin genes (LWS-1, LWS-2, LWS-3 and LWS-4) that are considered to be responsible for this sensory variation. However, the allelic variation of the opsin genes, particularly in terms of their absorption spectrum, has not been explored in wild populations. Thus, we examined nucleotide variations in the four M/LWS opsin genes as well as blue-sensitive SWS2-B and ultraviolet-sensitive SWS1 opsin genes for comparison and seven non-opsin nuclear loci as reference genes in 10 guppy populations from various light environments in Trinidad and Tobago. For the first time, we discovered a potential spectral variation (180 Ser/Ala) in LWS-1 that differed at an amino acid site known to affect the absorption spectra of opsins. Based on a coalescent simulation of the nucleotide variation of the reference genes, we showed that the interpopulation genetic differentiation of two opsin genes was significantly larger than the neutral expectation. Furthermore, this genetic differentiation was significantly related to differences in dissolved oxygen (DO) level, and it was not explained by the spatial distance between populations. The DO levels are correlated with eutrophication that possibly affects the color of aquatic environments. These results suggest that the population diversity of opsin genes is significantly driven by natural selection and that the guppy could adapt to various light environments through color vision changes.

Theoretical framework of population genetics with somatic mutations taken into account: application to copy number variations in humans.

Traditionally, population genetics focuses on the dynamics of frequencies of alleles acquired by mutations on germ-lines, because only such mutations are heritable. Typical genotyping experiments, however, use DNA from some somatic tissues such as blood, which harbors somatic mutations at the current generation in addition to germ-line mutations accumulated since the most recent common ancestor of the sample. This common practice may sometimes cause erroneous interpretations of polymorphism data, unless we properly understand the role of somatic mutations in population genetics. We here introduce a very basic theoretical framework of population genetics with somatic mutations taken into account. It is easy to imagine that somatic mutations at the current generation simply add individual-specific variations, as errors in mutation detection do. Our theory quantifies this increment under various conditions. We find that the major contribution of somatic mutations plus errors is to very rare variants, particularly to singletons. The relative contribution is markedly large when mutations are deleterious. Because negative selection also increases rare variants, it is important to distinguish the roles of these mutually confounding factors when we interpret the data, even after correcting for demography. We apply this theory to human copy number variations (CNVs), for which the composite effect of somatic mutations and errors may not be negligible. Using genome-wide CNV data, we demonstrate how the joint action of the two factors, selection and somatic mutations plus errors, shapes the observed pattern of polymorphism.

Efficient inference of bacterial strain trees from genome-scale multilocus data.

MOTIVATION: In bacterial evolution, inferring a strain tree, which is the evolutionary history of different strains of the same bacterium, plays a major role in analyzing and understanding the evolution of strongly isolated populations, population divergence and various evolutionary events, such as horizontal gene transfer and homologous recombination. Inferring a strain tree from multilocus data of these strains is exceptionally hard since, at this scale of evolution, processes such as homologous recombination result in a very high degree of gene tree incongruence. RESULTS: In this article we present a novel computational method for inferring the strain tree despite massive gene tree incongruence caused by homologous recombination. Our method operates in three phases, where in phase I a set of candidate strain-tree topologies is computed using the maximal cliques concept, in phase II divergence times for each of the topologies are estimated using mixed integer linear programming (MILP) and in phase III the optimal tree (or trees) is selected based on an optimality criterion. We have analyzed 1898 genes from nine strains of the Staphylococcus aureus bacteria, and identified a fully resolved (binary) strain tree with estimated divergence times, despite the high degrees of sequence identity at the nucleotide level and gene tree incongruence. Our method's efficiency makes it particularly suitable for analysis of genome-scale datasets, including those of strongly isolated populations which are usually very challenging to analyze. AVAILABILITY: We have implemented the algorithms in the PhyloNet software package, which is available publicly at http://bioinfo.cs.rice.edu/phylonet/.

The amounts of nucleotide variation within and between allelic classes and the reconstruction of the common ancestral sequence in a population.

The amounts of nucleotide variation within and between allelic classes were studied. The expectation and variance of the number of segregating sites and the expectation of the average number of pairwise differences among a sample of DNA sequences were obtained by using the theory of gene genealogy with no recombination. When the ancestral allelic class is unknown, it was found that the amount of variation within an allelic class increases with its frequency in the sample, while the amount of variation between two allelic classes is the largest when the two allelic classes exist equally. On the other hand, if we know the ancestral allelic class, as the frequency of the mutant allelic class increases, the amounts of variation within the mutant allelic class and between two allelic classes increase and the amount of variation within the ancestral allelic class decreases. As an example, we analyzed the polymorphism in the ND5 gene of Drosophila melanogaster and constructed the common ancestral sequence with high confidence, suggesting that the pattern of polymorphism within species gives useful information to know the ancestral sequence of the species.

Selection intensity against deleterious mutations in RNA secondary structures and rate of compensatory nucleotide substitutions.

A two-locus model of reversible mutations with compensatory fitness interactions is presented; single mutations are assumed to be deleterious but neutral in appropriate combinations. The expectation of the time of compensatory nucleotide substitutions is calculated analytically for the case of tight linkage between sites. It is shown that selection increases the substitution time dramatically when selection intensity Ns > 1, where N is the diploid population size and s the selection coefficient. Computer simulations demonstrate that recombination increases the substitution time, but the effect of recombination is small when selection is weak. The amount of linkage disequilibrium generated in the process of compensatory substitution is also investigated. It is shown that significant linkage disequilibrium is expected to be rare in natural populations. The model is applied to the mRNA secondary structure of the bicoid 3' untranslated region of Drosophila. It is concluded that average selection intensity Ns against single deleterious mutations is not likely to be much larger than 1.

Microsatellite polymorphism in natural populations of the wild plant Arabidopsis thaliana.

Variation in repeat number at 20 microsatellite loci of Arabidopsis thaliana was studied in a worldwide sample of 42 ecotypes to investigate the pattern and level of polymorphism in repetitive sequences in natural plant populations. There is a substantial amount of variation at microsatellite loci despite the selfing nature of this plant species. The average gene diversity was 0.794 and the average number of alleles per locus was 10.6. The distribution of alleles was centered around the mean of repeat number at most loci, but could not be regarded as normal. There was a significantly positive correlation between the number of repeats and the amount of variation. For most loci, the observed number of alleles was between the expected values of the infinite allele and stepwise mutation models. The two models were rejected by the sign test. Linkage disequilibrium was detected in 12.1% of the pairwise comparisons between loci. In phylogenetic tree, there was no association between ecotype and geographic origin. This result is consistent with the recent expansion of A. thaliana throughout the world.

DNA variation in the wild plant Arabidopsis thaliana revealed by amplified fragment length polymorphism analysis.

To investigate the level and pattern of DNA variation of Arabidopsis thaliana at the entire genome level, AFLP analysis was conducted for 38 ecotypes distributed throughout the world. Ten pairs of selective primers were used to detect a total of 472 bands, of which 374 (79. 2%) were polymorphic. The frequency distribution of polymorphic bands was skewed toward an excess of singleton variation. On the basis of AFLP variation, nucleotide diversity for the entire genome was estimated to be 0.0106, which was within the range reported previously for specific nuclear genes. The frequency distribution of pairwise distance was bimodal because of an ecotype (Fl-3) with a large number of unique bands. Linkage disequilibrium between polymorphic AFLPs was tested. The proportion of significant linkage disequilibria was close to random expectation after neglecting the ecotype Fl-3. This result indicates that the effect of recombination could not be ignored in this selfing species. A neighbor-joining tree was constructed on the basis of the AFLP variation. This tree has a star-like topology and shows no clear association between ecotype and geographic origin, suggesting a recent spread of this plant species and limited migration between its habitats.

A method for estimating nucleotide diversity from AFLP data.

A method for estimating the nucleotide diversity from AFLP data is developed by using the relationship between the number of nucleotide changes and the proportion of shared bands. The estimation equation is based on the assumption that GC-content is 0.5. Computer simulations, however, show that this method gives a reasonably accurate estimate even when GC-content deviates from 0.5, as long as the number of nucleotide changes per site (nucleotide diversity) is small. As an example, the nucleotide diversity of the wild yam, Dioscorea tokoro, was estimated. The estimated nucleotide diversity is 0.0055, which is larger than estimations from nucleotide sequence data for Adh and Pgi.

Species and recombination effects on DNA variability in the tomato genus.

Population genetics theory predicts that strong selection for rare, beneficial mutations or against frequent, deleterious mutations reduces polymorphism at linked neutral (or weakly selected) sites. The reduction of genetic variation is expected to be more severe when recombination rates are lower. In outbreeding species, low recombination rates are usually confined to certain chromosomal regions, such as centromeres and telomeres. In contrast, in predominantly selfing species, the rarity of double heterozygotes leads to a reduced effective recombination rate in the whole genome. We investigated the effects of restricted recombination on DNA polymorphism in these two cases, analyzing five Lycopersicon species with contrasting mating systems: L. chilense, L. hirsutum, L. peruvianum, L. chmielewskii, and L. pimpinellifolium, of which only the first three species have self-incompatibility alleles. In each species, we determined DNA sequence variation of five single-copy genes located in chromosomal regions with either high or low recombination rate. We found that the mating system has a highly significant effect on the level of polymorphism, whereas recombination has only a weak influence. The effect of recombination on levels of polymorphism in Lycopersicon is much weaker than in other well-studied species, including Drosophila. To explain these observations, we discuss a number of hypotheses, invoking selection, recombination, and demographic factors associated with the mating system. We also provide evidence that L. peruvianum, showing a level of polymorphism (almost 3%) that is comparable to the level of divergence in the whole genus, is the ancestral species from which the other species of the genus Lycopersicon have originated relatively recently.

Intragenic recombination in the Adh locus of the wild plant Arabidopsis thaliana.

Nucleotide variation in the Adh region of the wild plant Arabidopsis thaliana was analyzed in 17 ecotypes sampled worldwide to investigate DNA polymorphism in natural plant populations. The investigated 2.4-kb Adh region was divided into four blocks by intragenic recombinations between two parental sequence types that diverged 6.3 million years (Myr) ago, if the nucleotide mutation rate mu = 10(-9) is assumed. Within each block, dimorphism of segregating variations was observed with intermediate frequencies, which caused a substantial amount of nucleotide variation in A. thaliana at the species level. The first recombination introduced the divergent variation that resulted in dimorphism in this plant species approximately 3.3 Myr ago, and three subsequent intragenic combinations have occurred sporadically in approximately 1.1-Myr intervals. It was shown that there was only a limited number (six) of sequence types in this species and that no clear association was observed between sequence type and geographic origin. Taken together, these results suggest that A. thaliana has spread over the world only recently. It can be concluded that recombination played an important role in the evolutionary history of A. thaliana, especially through the generation of DNA polymorphism in the natural populations of this plant species.

The amount and pattern of DNA polymorphism under the neutral mutation hypothesis.

The amount and pattern of DNA polymorphism can give useful information on the maintenance mechanism of genetic variation at the DNA level. In this note we have shown the amount and pattern of DNA polymorphism expected under the neutral theory. The amount of DNA polymorphism can be estimated from the average number of nucleotide differences per site, the proportion of segregating sites, and so on. We have shown how to estimate theta from these quantities, where theta = 4Nv, N is the effective population size and v is the mutation rate per site per generation. We have also shown the expectations of the nucleotide variation within and between allelic classes.

DNA fingerprinting study on the intraspecific variation and the origin of Prunus yedoensis (Someiyoshino).

In order to investigate the intraspecific variation of Prunus yedoensis (Someiyoshino) and interspecific relationship among P. yedoensis, P. lannesiana (Oshimazakura) and P. pendula (Edohigan), DNA fingerprinting study was conducted by using two different kinds of probes, M13 repeat sequence and (GACA)4 synthetic oligonucleotide. In this study, 68 plants of P. yedoensis grown in 46 prefectures in Japan were investigated. All the P. yedoensis individuals investigated showed the completely same banding pattern, indicating their clonal origin from a single plant. On the other hand, each of P. lannesiana and P. pendula individuals investigated showed a unique banding pattern, suggesting a considerable amount of genetic variation in these two species. About 90% of bands in DNA fingerprints of P. yedoensis were detected in either P. lannesiana or P. pendula. This result supports the hypothesis that P. yedoensis is an interspecific hybrid between P. lannesiana and P. pendula. From those results, it is concluded that P. yedoensis was produced only once through hybridization between P. lannesiana and P. pendula, and that this particular hybrid plant has been spread vegetatively all over Japan.

Nucleotide polymorphism in the acidic chitinase locus (ChiA) region of the wild plant Arabidopsis thaliana.

To investigate DNA variation in natural plant populations, a 1.8-kb region of the acidic chitinase locus (ChiA)was analyzed for 17 ecotypes of Arabidopsis thaliana sampled worldwide and 3 Arabis species in Japan. As in the Adh region, dimorphism was detected throughout the investigated ChiA region, suggesting the possibility that dimorphic DNA variation exists in the entire nuclear genome of A. thaliana. The ChiA region was divided into two blocks by an intragenic recombination between two parental sequence types, which diverged 7.4 MYA under the assumption that nucleotide mutation rate per site per year is mu = 10(-9). Nucleotide diversity in the entire ChiA region was 0.0104. Tajima's test was significantly negative for both nucleotide and indel variations, which was manifested as an excess of unique polymorphisms. However, the level and pattern of polymorphism in the ChiA region were inconsistent with simple theoretical explanations. The HKA test detected no difference in the levels of intra- and interspecific variations between the ChiA and Adh regions. In the ChiA coding region, no difference in the patterns of synonymous and replacement variation was found in intra- and interspecific comparisons by the MK test. Although it was difficult to determine the exact genetic mechanism acting on the ChA locus, these results suggested that the ChA locus region was under the same genetic mechanism before and after the establishment of A. thaliana as a species.

Intra- and interspecific DNA variation and codon bias of the alcohol dehydrogenase (Adh) locus in Arabis and Arabidopsis species.

Sequence variation at the alcohol dehydrogenase (Adh) locus was analyzed for six species each of the genera Arabis and Arabidopsis. Phylogenetic analysis showed that investigated species were grouped into three clusters, and the generic classification did not correspond to the clusterings. The results indicated that the genera could not be distinguished on the basis of the Adh variation. A significant difference in the ratio of silent to replacement sites was detected by MK test in two comparisons, with Arabidopsis thaliana polymorphism due to excess silent divergence. Silent changes were predominant in the evolution of the Adh locus in Arabis and Arabidopsis. To infer evolutionary significance of silent substitutions, codon bias was studied. The degree of codon bias of the Adh region was relatively constant over Arabis and Arabidopsis species. "Preferred" codons of A. thaliana were determined. No evidence of natural selection on codon change was detected in the Adh regions of A. thaliana and Arabis gemmifera.