Satoshi Horai

Bio: Satoshi Horai is an academic researcher from Graduate University for Advanced Studies. The author has contributed to research in topics: Mitochondrial DNA & Mitochondrial myopathy. The author has an hindex of 41, co-authored 81 publications receiving 7994 citations. Previous affiliations of Satoshi Horai include University of Arizona & Nagoya Institute of Technology.

A mutation in the tRNA Leu(UUR) gene associated with the MELAS subgroup of mitochondrial encephalomyopathies

Abstract: MITOCHONDRIAL encephalomyopathies are usually divided into three distinct clinical subgroups: (1) mitochondrial myopathy, encephalopathy, lactic acidosis and stroke-like episodes (MELAS); (2) myoclonus epilepsy associated with ragged-red fibres (MERRF); and (3) chronic progressive external ophthalmoplegia (CPEO) including Kearns-Sayre syndrome1–5. Large deletions of human mitochondrial DNA and a transition mutation at the mitochondrial transfer RNAlys gene give rise to CPEO including Kearns–Sayre syndrome6–8 and MERRF9,10, respectively. Here we report an A-to-G transition mutation at nucleotide pair 3,243 in the dihydrouridine loop of mitochondrial tRNALeu(UUR) that is specific to patients with MELAS. Because this mutation creates an Apal restriction site, we could perform a simple molecular diagnostic test for the disease. The mutation was present in 26 out of 31 independent MELAS patients and 1 out of 29 CPEO patients, but absent in the 5 MERRF and 50 controls tested. Southern blot analysis confirmed that the mutant DNA always coexists with the wild-type DNA (heteroplasmy).

Recent African origin of modern humans revealed by complete sequences of hominoid mitochondrial DNAs.

Abstract: We analyzed the complete mitochondrial DNA (mtDNA) sequences of three humans (African, European, and Japanese), three African apes (common and pygmy chimpanzees, and gorilla), and one orangutan in an attempt to estimate most accurately the substitution rates and divergence times of hominoid mtDNAs. Nonsynonymous substitutions and substitutions in RNA genes have accumulated with an approximately clock-like regularity. From these substitutions and under the assumption that the orangutan and African apes diverged 13 million years ago, we obtained a divergence time for humans and chimpanzees of 4.9 million years. This divergence time permitted calibration of the synonymous substitution rate (3.89 x 10(-8)/site per year). To obtain the substitution rate in the displacement (D)-loop region, we compared the three human mtDNAs and measured the relative abundance of substitutions in the D-loop region and at synonymous sites. The estimated substitution rate in the D-loop region was 7.00 x 10(-8)/site per year. Using both synonymous and D-loop substitutions, we inferred the age of the last common ancestor of the human mtDNAs as 143,000 +/- 18,000 years. The shallow ancestry of human mtDNAs, together with the observation that the African sequence is the most diverged among humans, strongly supports the recent African origin of modern humans, Homo sapiens sapiens.

Mitochondrial myopathy, encephalopathy, lactic acidosis, and stroke‐like episodes (MELAS) A correlative study of the clinical features and mitochondrial DNA mutation

Abstract: We studied 40 MELAS patients (21 male and 19 female) to characterize the clinical features and biochemical and muscle biopsy findings related to the mtDNA mutation at the nucleotide position of 3,243, the most common genetic defect in MELAS. The most frequent symptom was episodic sudden headache with vomiting and convulsions, which commonly affected patients aged 5 to 15 years (80%). Biochemical defects in the muscle were variable; 13 patients had complex I, seven complex IV, and four complexes I + IV deficiencies. In four muscle biopsies without ragged-red fibers or any enzyme defect, we based the diagnosis on the identification of strongly SDH-reactive blood vessels, which occurred in 87.5% of the biopsies. The mtDNA mutation was present in 32 of 40 patients (80%). We conclude that there are no clinical and pathologic differences between the patients with and without this mtDNA mutation.

Y chromosomal DNA variation and the peopling of Japan.

Abstract: Four loci mapping to the nonrecombining portion of the Y chromosome were genotyped in Japanese populations from Okinawa, the southernmost island of Japan; Shizuoka and Aomori on the main island of Honshu; and a small sample of Taiwanese. The Y Alu polymorphic (YAP) element is present in 42% of the Japanese and absent in the Taiwanese, confirming the irregular distribution of this polymorphism in Asia. Data from the four loci were used to determine genetic distances among populations, construct Y chromosome haplotypes, and estimate the degree of genetic diversity in each population and on different Y chromosome haplotypes. Evolutionary analysis of Y haplotypes suggests that polymorphisms at the YAP (DYS287) and DXYS5Y loci originated a single time, whereas restriction patterns at the DYS1 locus and microsatellite alleles at the DYS19 locus arose more than once. Genetic distance analysis indicated that the Okinawans are differentiated from Japanese living on Honshu. The data support the hypotheses that modern Japanese populations have resulted from distinctive genetic contributions involving the ancient Jomon people and Yayoi immigrants from Korea or mainland China, with Okinawans experiencing the least amount of admixture with the Yayoi. It is suggested that YAP+ chromosomes migrated to Japan with the Jomon people > 10,000 years ago and that a large infusion of YAP- chromosomes entered Japan with the Yayoi migration starting 2,300 years ago. Different degrees of genetic diversity carried by these two ancient chromosomal lineages may be explained by the different life-styles (hunter-gatherer versus agriculturalist). of the migrant groups, the size of the founding populations, and the antiquities of the founding events.

PAML 4: Phylogenetic Analysis by Maximum Likelihood

Abstract: PAML, currently in version 4, is a package of programs for phylogenetic analyses of DNA and protein sequences using maximum likelihood (ML). The programs may be used to compare and test phylogenetic trees, but their main strengths lie in the rich repertoire of evolutionary models implemented, which can be used to estimate parameters in models of sequence evolution and to test interesting biological hypotheses. Uses of the programs include estimation of synonymous and nonsynonymous rates (d(N) and d(S)) between two protein-coding DNA sequences, inference of positive Darwinian selection through phylogenetic comparison of protein-coding genes, reconstruction of ancestral genes and proteins for molecular restoration studies of extinct life forms, combined analysis of heterogeneous data sets from multiple gene loci, and estimation of species divergence times incorporating uncertainties in fossil calibrations. This note discusses some of the major applications of the package, which includes example data sets to demonstrate their use. The package is written in ANSI C, and runs under Windows, Mac OSX, and UNIX systems. It is available at -- (http://abacus.gene.ucl.ac.uk/software/paml.html).

Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees.

Abstract: Examining the pattern of nucleotide substitution for the control region of mitochondrial DNA ( mtDNA ) in humans and chimpanzees, we developed a new mathematical method for estimating the number of transitional and transversional substitutions per site, as well as the total number of nucleotide substitutions. In this method, excess transitions, unequal nucleotide frequencies, and variation of substitution rate among different sites are all taken into account. Application of this method to human and chimpanzee data suggested that the transition / transversion ratio for the entire control region was - 15 and nearly the same for the two species. The 95% confidence interval of the age of the common ancestral mtDNA was estimated to be 80,000-480,000 years in humans and 0.57-2.72 Myr in common chimpanzees.

BEAST 2: A Software Platform for Bayesian Evolutionary Analysis

Abstract: We present a new open source, extensible and flexible software platform for Bayesian evolutionary analysis called BEAST 2. This software platform is a re-design of the popular BEAST 1 platform to correct structural deficiencies that became evident as the BEAST 1 software evolved. Key among those deficiencies was the lack of post-deployment extensibility. BEAST 2 now has a fully developed package management system that allows third party developers to write additional functionality that can be directly installed to the BEAST 2 analysis platform via a package manager without requiring a new software release of the platform. This package architecture is showcased with a number of recently published new models encompassing birth-death-sampling tree priors, phylodynamics and model averaging for substitution models and site partitioning. A second major improvement is the ability to read/write the entire state of the MCMC chain to/from disk allowing it to be easily shared between multiple instances of the BEAST software. This facilitates checkpointing and better support for multi-processor and high-end computing extensions. Finally, the functionality in new packages can be easily added to the user interface (BEAUti 2) by a simple XML template-based mechanism because BEAST 2 has been re-designed to provide greater integration between the analysis engine and the user interface so that, for example BEAST and BEAUti use exactly the same XML file format.

Population growth makes waves in the distribution of pairwise genetic differences.

Abstract: Episodes of population growth and decline leave characteristic signatures in the distribution of nucleotide (or restriction) site differences between pairs of individuals. These signatures appear in histograms showing the relative frequencies of pairs of individuals who differ by i sites, where i = 0, 1, .... In this distribution an episode of growth generates a wave that travels to the right, traversing 1 unit of the horizontal axis in each 1/2u generations, where u is the mutation rate. The smaller the initial population, the steeper will be the leading face of the wave. The larger the increase in population size, the smaller will be the distribution's vertical intercept. The implications of continued exponential growth are indistinguishable from those of a sudden burst of population growth Bottlenecks in population size also generate waves similar to those produced by a sudden expansion, but with elevated uppertail probabilities. Reductions in population size initially generate L-shaped distributions with high probability of identity, but these converge rapidly to a new equilibrium. In equilibrium populations the theoretical curves are free of waves. However, computer simulations of such populations generate empirical distributions with many peaks and little resemblance to the theory. On the other hand, agreement is better in the transient (nonequilibrium) case, where simulated empirical distributions typically exhibit waves very similar to those predicted by theory. Thus, waves in empirical distributions may be rich in information about the history of population dynamics.

Model Selection and Model Averaging in Phylogenetics: Advantages of Akaike Information Criterion and Bayesian Approaches Over Likelihood Ratio Tests

Abstract: Model selection is a topic of special relevance in molecular phylogenetics that affects many, if not all, stages of phylogenetic inference. Here we discuss some fundamental concepts and techniques of model selection in the context of phylogenetics. We start by reviewing different aspects of the selection of substitution models in phylogenetics from a theoretical, philosophical and practical point of view, and summarize this comparison in table format. We argue that the most commonly implemented model selection approach, the hierarchical likelihood ratio test, is not the optimal strategy for model selection in phylogenetics, and that approaches like the Akaike Information Criterion (AIC) and Bayesian methods offer important advantages. In particular, the latter two methods are able to simultaneously compare multiple nested or nonnested models, assess model selection uncertainty, and allow for the estimation of phylogenies and model parameters using all available models (model-averaged inference or multimodel inference). We also describe how the relative importance of the different parameters included in substitution models can be depicted. To illustrate some of these points, we have applied AIC-based model averaging to 37 mitochondrial DNA sequences from the subgenus Ohomopterus (genus Carabus) ground beetles described by Sota and Vogler (2001). (AIC; Bayes factors; BIC; likelihood ratio tests; model averaging; model uncertainty; model selection; multimodel inference.) It is clear that models of nucleotide substitution (henceforth models of evolution) play a significant role in molecular phylogenetics, particularly in the context of distance, maximum likelihood (ML), and Bayesian es- timation. We know that the use of one or other model affects many, if not all, stages of phylogenetic inference. For example, estimates of phylogeny, substitution rates, bootstrap values, posterior probabilities, or tests of the molecular clock are clearly influenced by the model of evolution used in the analysis (Buckley, 2002; Buckley