A multiple sequence alignment program, MAFFT, has been developed. The CPU time is drastically reduced as compared with existing methods. MAFFT includes two novel techniques. (i) Homologous regions are rapidly identified by the fast Fourier transform (FFT), in which an amino acid sequence is converted to a sequence composed of volume and polarity values of each amino acid residue. (ii) We propose a simplified scoring system that performs well for reducing CPU time and increasing the accuracy of alignments even for sequences having large insertions or extensions as well as distantly related sequences of similar length. Two different heuristics, the progressive method (FFT-NS-2) and the iterative refinement method (FFT-NS-i), are implemented in MAFFT. The performances of FFT-NS-2 and FFT-NS-i were compared with other methods by computer simulations and benchmark tests; the CPU time of FFT-NS-2 is drastically reduced as compared with CLUSTALW with comparable accuracy. FFT-NS-i is over 100 times faster than T-COFFEE, when the number of input sequences exceeds 60, without sacrificing the accuracy.

/pdf/mafft-a-novel-method-for-rapid-multiple-sequence-alignment-485azc27cn.pdf

MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform

Statistical method for testing the neutral mutation hypothesis by DNA polymorphism.

I reexamine the use of isolation by distance models as a basis for the estimation of demographic parameters from measures of population subdivision. To that aim, I first provide results for values of F-statistics in one-dimensional models and coalescence times in two-dimensional models, and make more precise earlier results for F-statistics in two-dimensional models and coalescence times in one-dimensional models. Based on these results, I propose a method of data analysis involving the regression of F(ST)/(1 - F(ST)) estimates for pairs of subpopulations on geographic distance for populations along linear habitats or logarithm of distance for populations in two-dimensional habitats. This regression provides in principle an estimate of the product of population density and second moment of parental axial distance. In two cases where comparison to direct estimates is possible, the method proposed here is more satisfactory than previous indirect methods.

https://www.genetics.org/content/genetics/145/4/1219.full.pdf

Genetic Differentiation and Estimation of Gene Flow from F -Statistics Under Isolation by Distance

We previously developed a cladistic approach to identify subsets of haplotypes defined by restriction endonuclease mapping or DNA sequencing that are associated with significant phenotypic deviations. Our approach was limited to segments of DNA in which little recombination occurs. In such cases, a cladogram can be constructed from the restriction site or sequence data that represents the evolutionary steps that interrelate the observed haplotypes. The cladogram is used to define a nested statistical design to identify mutational steps associated with significant phenotypic deviations. The central assumption behind this strategy is that any undetected mutation causing a phenotypic effect is embedded within the same evolutionary history that is represented by the cladogram. The power of this approach depends upon the confidence one has in the particular cladogram used to draw inferences. In this paper, we present a strategy for estimating the set of cladograms that are consistent with a particular sample of either restriction site or nucleotide sequence data and that includes the possibility of recombination. We first evaluate the limits of parsimony in constructing cladograms. Once these limits have been determined, we construct the set of parsimonious and nonparsimonious cladograms that is consistent with these limits. Our estimation procedure also identifies haplotypes that are candidates for being products of recombination. If recombination is extensive, our algorithm subdivides the DNA region into two or more subsections, each having little or no internal recombination. We apply this estimation procedure to three data sets to illustrate varying degrees of cladogram ambiguity and recombination.

https://www.genetics.org/content/genetics/132/2/619.full.pdf

A cladistic analysis of phenotypic associations with haplotypes inferred from restriction endonuclease mapping and DNA sequence data. III. Cladogram estimation.

A multilocus sequence typing (MLST) scheme has been developed for Staphylococcus aureus. The sequences of internal fragments of seven housekeeping genes were obtained for 155 S. aureus isolates from patients with community-acquired and hospital-acquired invasive disease in the Oxford, United Kingdom, area. Fifty-three different allelic profiles were identified, and 17 of these were represented by at least two isolates. The MLST scheme was highly discriminatory and was validated by showing that pairs of isolates with the same allelic profile produced very similar SmaI restriction fragment patterns by pulsed-field gel electrophoresis. All 22 isolates with the most prevalent allelic profile were methicillin-resistant S. aureus (MRSA) isolates and had allelic profiles identical to that of a reference strain of the epidemic MRSA clone 16 (EMRSA-16). Four MRSA isolates that were identical in allelic profile to the other major epidemic MRSA clone prevalent in British hospitals (clone EMRSA-15) were also identified. The majority of isolates (81%) were methicillin-susceptible S. aureus (MSSA) isolates, and seven MSSA clones included five or more isolates. Three of the MSSA clones included at least five isolates from patients with community-acquired invasive disease and may represent virulent clones with an increased ability to cause disease in otherwise healthy individuals. The most prevalent MSSA clone (17 isolates) was very closely related to EMRSA-16, and the success of the latter clone at causing disease in hospitals may be due to its emergence from a virulent MSSA clone that was already a major cause of invasive disease in both the community and hospital settings. MLST provides an unambiguous method for assigning MRSA and MSSA isolates to known clones or assigning them as novel clones via the Internet.

Multilocus sequence typing for characterization of methicillin-resistant and methicillin-susceptible clones of Staphylococcus aureus.

Possible emergence of new geminiviruses by frequent recombination.

Statistical tests for detecting gene conversion are described for a sample of homologous DNA sequences. The tests are based on imbalances in the distribution of segments on which some pair of sequences agrees. The methods automatically control for variable mutation rates along the genome and do not depend on a priori choices of potentially monophyletic subsets of the sample. The tests show strong evidence for multiple intragenic conversion events at two loci in Escherichia coli. The gnd locus in E. coli shows a highly significant excess of maximal segments of length 70-200 bp, which suggests conversion events of that size. The data also indicate that the rate of these short conversion events might be of the order of neutral mutation rate. There is also evidence for correlated mutation in adjacent codon positions. The same tests applied to a locus in an RNA virus were negative.

https://www.math.wustl.edu/~sawyer/resint2/recostat.pdf

Statistical tests for detecting gene conversion.

Frequencies of mutant sites are modeled as a Poisson random field in two species that share a sufficiently recent common ancestor. The selective effect of the new alleles can be favorable, neutral, or detrimental. The model is applied to the sample configurations of nucleotides in the alcohol dehydrogenase gene (Adh) in Drosophila simulans and Drosophila yakuba. Assuming a synonymous mutation rate of 1.5 x 10(-8) per site per year and 10 generations per year, we obtain estimates for the effective population size (N(e) = 6.5 x 10(6)), the species divergence time (tdiv = 3.74 million years), and an average selection coefficient (sigma = 1.53 x 10(-6) per generation for advantageous or mildly detrimental replacements), although it is conceivable that only two of the amino acid replacements were selected and the rest neutral. The analysis, which includes a sampling theory for the independent infinite sites model with selection, also suggests the estimate that the number of amino acids in the enzyme that are susceptible to favorable mutation is in the range 2-23 at any one time. The approach provides a theoretical basis for the use of a 2 x 2 contingency table to compare fixed differences and polymorphic sites with silent sites and amino acid replacements.

https://www.genetics.org/content/genetics/132/4/1161.full.pdf

Population genetics of polymorphism and divergence.

We present a worldwide analysis of humid tropical forest dynamics and tree species richness. New tree mortality, recruitment, and species richness data include the most dynamic and diverse mature tropical forests known. Twenty-five sites show a strong tendency for the most species-rich forests to be dynamic and aseasonal. Mean annual tree mortality and recruitment-turnover-is the most predictive factor of species richness, implying that small-scale disturbance helps regulate tropical forest diversity. Turnover rates are also closely related to the amount of basal area turnover in mature tropical forests. Therefore the contribution of small-scale disturbance to maintaining tropical forest diversity may ultimately be driven by ecosystem productivity.

https://www.pnas.org/content/pnas/91/7/2805.full.pdf

Dynamics and species richness of tropical rain forests.

Population geneticists have long sought to estimate the distribution of selection intensities among genes of diverse function across the genome. Only recently have DNA sequencing and analytical techniques converged to make this possible. Important advances have come from comparing genetic variation within species (polymorphism) with fixed differences between species (divergence). These approaches have been used to examine individual genes for evidence of selection. Here we use the fact that the time since species divergence allows combination of data across genes. In a comparison of amino-acid replacements among species of the mustard weed Arabidopsis with those among species of the fruitfly Drosophila, we find evidence for predominantly beneficial gene substitutions in Drosophila but predominantly detrimental substitutions in Arabidopsis. We attribute this difference to the Arabidopsis mating system of partial self-fertilization, which corroborates a prediction of population genetics theory that species with a high frequency of inbreeding are less efficient in eliminating deleterious mutations owing to their reduced effective population size.

Stanley Sawyer

Papers

Possible emergence of new geminiviruses by frequent recombination.

Statistical tests for detecting gene conversion.

Population genetics of polymorphism and divergence.

Dynamics and species richness of tropical rain forests.

The cost of inbreeding in Arabidopsis