The Nadeau-Taylor Chromosomal Breakage Model Revisited

Total Page:16

File Type:pdf, Size:1020Kb

The Nadeau-Taylor Chromosomal Breakage Model Revisited Transforming Men into Mice: the Nadeau-Taylor Chromosomal Breakage Model Revisited Pavel Pevzner Glenn Tesler Department of Computer Science Department of Computer Science University of California, San Diego University of California, San Diego La Jolla, CA 92093-0114 La Jolla, CA 92093-0114 [email protected] [email protected] ABSTRACT order) and published a milestone paper with a rearrange- Although analysis of genome rearrangements was pioneered ment scenario for the species D. pseudoobscura and D. mi- by Dobzhansky and Sturtevant 65 years ago, we still know randa. Every genome rearrangement study involves solving very little about the rearrangement events that produced the a combinatorial puzzle to find a series of genome rearrange- existing varieties of genomic architectures. The genomic se- ments to transform one genome into another. The problem quences of human and mouse provide evidence for a larger of finding the minimum number of reversals to transform one number of rearrangements than previously thought and shed unichromosomal genome into another is called the reversal some light on previously unknown features of mammalian distance problem. evolution. In particular, they reveal extensive re-use of Kececioglu and Sankoff, 1993 [17] were the first to recog- breakpoints from the same relatively short regions. Our nize the importance of combinatorial dependencies between analysis implies the existence of a large number of very short different breakpoints and to come up with an approximation “hidden” synteny blocks that were invisible in comparative algorithm for the reversal distance problem. Hannenhalli mapping data and were not taken into account in previous and Pevzner, 1995, 1999 [14] developed a polynomial algo- studies of chromosome evolution. These blocks are defined rithm for the reversal distance problem and further extended by closely located breakpoints and are often hard to detect. it to the genomic distance problem, i.e., finding a most par- Our result is in conflict with the widely accepted random simonious scenario for multichromosomal genomes under re- breakage model of chromosomal evolution. We suggest a versals, translocations, fusions, and fissions of chromosomes new “fragile breakage” model of chromosome evolution that [13, 43]. Later on, the Hannenhalli-Pevzner algorithm was postulates that breakpoints are chosen from relatively short further optimized and extended for other applications in [5, fragile regions that have much higher propensity for rear- 16, 11, 1, 4]. rangements than the rest of the genome. Even before Sankoff and colleagues introduced the com- binatorial approach to rearrangement studies, Nadeau and Taylor pioneered a statistical approach. In a landmark pa- Categories and Subject Descriptors per, Nadeau and Taylor, 1984 [29] introduced the notion J.3 [Life and Medical Sciences]: Biology and Genetics of conserved segments (segments with preserved gene or- ders) and estimated that there are roughly 180 conserved General Terms segments in human and mouse. In the same paper they pro- vided convincing arguments in favor of the random breakage Algorithms, Theory model of genomic evolution postulated by Ohno, 1973 [30]. The model assumes a random (i.e, uniform and independent) Keywords distribution of chromosome rearrangement breakpoints and is supported by the observation that the lengths of synteny breakpoint re-use, evolution, genome rearrangements blocks shared by human and mouse are well fitted by the pre- dicted distribution imposed by the random breakage model. 1. INTRODUCTION Since the model was first introduced in [30], it has been Analysis of genome rearrangements in molecular evolution analyzed by Nadeau and others [29, 27, 28, 39], and has was pioneered by Dobzhansky and Sturtevant, 1938 [10], become widely accepted. It was further supported by stud- who introduced the notion of a breakpoint (disruption of gene ies of significantly larger datasets that confirmed that newly discovered synteny blocks still fit the predicted exponential distribution very well [26, 9, 20, 24, 33]. These studies, with progressively increasing levels of resolution, transformed the Permission to make digital or hard copies of all or part of this work for random breakage model into the de facto theory of chromo- personal or classroom use is granted without fee provided that copies are some evolution. not made or distributed for profit or commercial advantage and that copies The arguments in favor of the random breakage model bear this notice and the full citation on the first page. To copy otherwise, to republish, to post on servers or to redistribute to lists, requires prior specific usually proceed as follows. One first constructs the distri- permission and/or a fee. bution of lengths of conserved segments and fits the resulting RECOMB’03, April 10–13, 2003, Berlin, Germany. histogram with the theoretical distribution predicted by the Copyright 2003 ACM 1-58113-635-8/03/0004 ...$5.00. 247 random breakage model. An important implication of this (a) (b) Lengths of synteny blocks (w/o hidden blocks) Lengths of synteny blocks (with hidden blocks) model is that the segment lengths approximate an expo- 80 300 1 −x/L nential distribution with density function f(x)= L e , 250 60 where L is the average length of all segments. Techni- 200 cally, Nadeau and Taylor, 1984 [29] did not have information 150 about all segments since most of them were still undiscov- 40 Frequency Frequency 100 ered in 1984. However, they were able to estimate L (and 20 therefore the number of still undiscovered segments) from 50 0 0 the small set of already discovered segments. The relatively 0 10 20 30 40 50 60 70 80 0 10 20 30 40 50 60 70 80 small departure from an exponential distribution was at- Block length (Mb) Block length (Mb) tributed to missing information about some conserved seg- ments. Of course, there was always a danger that newly dis- Figure 1: (a) Histogram of synteny block lengths covered segments would shift this estimate and even deviate in human for 281 synteny blocks of length at least 1 1 from the exponential distribution predicted by the model. Mb, fitted by an exponential distribution with mean However, this did not happen in the past, and the random block length L =9.6 Mb. (b) The same histogram breakage model was reinforced in a number of influential superimposed with the 190 “hidden” synteny blocks studies in the last decade. As a result, the Nadeau-Taylor revealed by genome rearrangement analysis, under predictions are viewed as among the most significant results the assumption that all “hidden” blocks are short, in “...the history and development of the mouse as a research i.e., less than 1 Mb in length. tool” (Pennisi, 2000 [31]). There is a conceptual difference between the Nadeau and Taylor, 1984 [29] statistical approach to studies of chro- closely located breakpoints that cannot be explained by the mosomal history (that is not concerned with the details random breakage model. of rearrangement history) and the combinatorial approach The surprisingly large number of breakpoint clumps is an that attempts to infer the rearrangement scenario. Sankoff, argument in favor of a different model of chromosome evo- 1999 [34] raised the problem of integrating these approaches, lution that we call the fragile breakage model.Thismodel which had not been done before. In this paper we attempt postulates that the breakpoints mainly occur within rela- to combine the statistical and combinatorial approaches and tively short fragile regions (hot spots of rearrangements). demonstrate that the combined analysis reveals evidence The existence of some fragile regions at the population level against the random breakage model. was supported by previous studies of cancer and infertility The draft human and mouse sequences reveal many pre- [8, 36], but the extent of this phenomenon in molecular evo- viously undiscovered synteny blocks and put the random lution became clear only after the human and mouse DNA breakage model to a new test. In particular, they reveal sequences became available. If one assumes the fragile re- 281 synteny blocks shared by human and mouse of size at gions are uniformly distributed through the genome then least 1 Mb (Pevzner and Tesler, 2003 [32]). Although the the fragile and random breakage models lead to identical number of synteny blocks is higher than the Nadeau-Taylor estimates for the number of long segments (e.g., segments predictions, the lengths of the blocks still fit the exponential longer than 1–2 Mb). In some sense, the random break- distribution (Fig. 1a), another argument in favor of the ran- age model can be viewed as an excellent null hypothesis dom breakage model. However, a different type of evidence for a certain level of resolution and genome heterogeneities. derived from genome rearrangement studies reveals an un- However, the random breakage and fragile breakage model expectedly large number of closely located breakpoints that generate very different predictions when it comes to short cannot be explained by the random breakage model. This segments that were below the granularity level of previous analysis implies that in addition to the segments shown in comparative mapping studies. Fig. 1a, there are another 190 “short” synteny blocks, typ- ically below 1 Mb in length. These blocks were never dis- 2. SYNTENY BLOCKS covered in the comparative mapping studies, and moreover, DNA sequences provide evidence that the human and most of them are hard to find even with available human 2 mouse genomes are significantly more rearranged than pre- and mouse sequences. The existence of these blocks im- viously thought. Moreover, they indicate that a large pro- mediately implies that an exponential distribution is not a portion of previously identified conserved segments are not good fit to reality (Fig. 1b). In other words, rearrangement really conserved since there is evidence of multiple micro- analysis of the human and mouse genomes reveals clumps of rearrangements (Mural et al., 2002 [24]).
Recommended publications
  • Grammar String: a Novel Ncrna Secondary Structure Representation
    Grammar string: a novel ncRNA secondary structure representation Rujira Achawanantakun, Seyedeh Shohreh Takyar, and Yanni Sun∗ Department of Computer Science and Engineering, Michigan State University, East Lansing, MI 48824 , USA ∗Email: [email protected] Multiple ncRNA alignment has important applications in homologous ncRNA consensus structure derivation, novel ncRNA identification, and known ncRNA classification. As many ncRNAs’ functions are determined by both their sequences and secondary structures, accurate ncRNA alignment algorithms must maximize both sequence and struc- tural similarity simultaneously, incurring high computational cost. Faster secondary structure modeling and alignment methods using trees, graphs, probability matrices have thus been developed. Despite promising results from existing ncRNA alignment tools, there is a need for more efficient and accurate ncRNA secondary structure modeling and alignment methods. In this work, we introduce grammar string, a novel ncRNA secondary structure representation that encodes an ncRNA’s sequence and secondary structure in the parameter space of a context-free grammar (CFG). Being a string defined on a special alphabet constructed from a CFG, it converts ncRNA alignment into sequence alignment with O(n2) complexity. We align hundreds of ncRNA families from BraliBase 2.1 using grammar strings and compare their consensus structure with Murlet using the structures extracted from Rfam as reference. Our experiments have shown that grammar string based multiple sequence alignment competes favorably in consensus structure quality with Murlet. Source codes and experimental data are available at http://www.cse.msu.edu/~yannisun/grammar-string. 1. INTRODUCTION both the sequence and structural conservations. A successful application of SCFG is ncRNA classifica- Annotating noncoding RNAs (ncRNAs), which are tion, which classifies query sequences into annotated not translated into protein but function directly as ncRNA families such as tRNA, rRNA, riboswitch RNA, is highly important to modern biology.
    [Show full text]
  • 120421-24Recombschedule FINAL.Xlsx
    Friday 20 April 18:00 20:00 REGISTRATION OPENS in Fira Palace 20:00 21:30 WELCOME RECEPTION in CaixaForum (access map) Saturday 21 April 8:00 8:50 REGISTRATION 8:50 9:00 Opening Remarks (Roderic GUIGÓ and Benny CHOR) Session 1. Chair: Roderic GUIGÓ (CRG, Barcelona ES) 9:00 10:00 Richard DURBIN The Wellcome Trust Sanger Institute, Hinxton UK "Computational analysis of population genome sequencing data" 10:00 10:20 44 Yaw-Ling Lin, Charles Ward and Steven Skiena Synthetic Sequence Design for Signal Location Search 10:20 10:40 62 Kai Song, Jie Ren, Zhiyuan Zhai, Xuemei Liu, Minghua Deng and Fengzhu Sun Alignment-Free Sequence Comparison Based on Next Generation Sequencing Reads 10:40 11:00 178 Yang Li, Hong-Mei Li, Paul Burns, Mark Borodovsky, Gene Robinson and Jian Ma TrueSight: Self-training Algorithm for Splice Junction Detection using RNA-seq 11:00 11:30 coffee break Session 2. Chair: Bonnie BERGER (MIT, Cambrige US) 11:30 11:50 139 Son Pham, Dmitry Antipov, Alexander Sirotkin, Glenn Tesler, Pavel Pevzner and Max Alekseyev PATH-SETS: A Novel Approach for Comprehensive Utilization of Mate-Pairs in Genome Assembly 11:50 12:10 171 Yan Huang, Yin Hu and Jinze Liu A Robust Method for Transcript Quantification with RNA-seq Data 12:10 12:30 120 Zhanyong Wang, Farhad Hormozdiari, Wen-Yun Yang, Eran Halperin and Eleazar Eskin CNVeM: Copy Number Variation detection Using Uncertainty of Read Mapping 12:30 12:50 205 Dmitri Pervouchine Evidence for widespread association of mammalian splicing and conserved long range RNA structures 12:50 13:10 169 Melissa Gymrek, David Golan, Saharon Rosset and Yaniv Erlich lobSTR: A Novel Pipeline for Short Tandem Repeats Profiling in Personal Genomes 13:10 13:30 217 Rory Stark Differential oestrogen receptor binding is associated with clinical outcome in breast cancer 13:30 15:00 lunch break Session 3.
    [Show full text]
  • Michael S. Waterman: Breathing Mathematics Into Genes >>>
    ISSUE 13 Newsletter of Institute for Mathematical Sciences, NUS 2008 Michael S. Waterman: Breathing Mathematics into Genes >>> setting up of the Center for Computational and Experimental Genomics in 2001, Waterman and his collaborators and students continue to provide a road map for the solution of post-genomic computational problems. For his scientific contributions he was elected fellow or member of prestigious learned bodies like the American Academy of Arts and Sciences, National Academy of Sciences, American Association for the Advancement of Science, Institute of Mathematical Statistics, Celera Genomics and French Acadèmie des Sciences. He was awarded a Gairdner Foundation International Award and the Senior Scientist Accomplishment Award of the International Society of Computational Biology. He currently holds an Endowed Chair at USC and has held numerous visiting positions in major universities. In addition to research, he is actively involved in the academic and social activities of students as faculty master Michael Waterman of USC’s International Residential College at Parkside. Interview of Michael S. Waterman by Y.K. Leong Waterman has served as advisor to NUS on genomic research and was a member of the organizational committee Michael Waterman is world acclaimed for pioneering and of the Institute’s thematic program Post-Genome Knowledge 16 fundamental work in probability and algorithms that has Discovery (Jan – June 2002). On one of his advisory tremendous impact on molecular biology, genomics and visits to NUS, Imprints took the opportunity to interview bioinformatics. He was a founding member of the Santa him on 7 February 2007. The following is an edited and Cruz group that launched the Human Genome Project in enhanced version of the interview in which he describes the 1990, and his work was instrumental in bringing the public excitement of participating in one of the greatest modern and private efforts of mapping the human genome to their scientific adventures and of unlocking the mystery behind completion in 2003, two years ahead of schedule.
    [Show full text]
  • Gene and Genome Duplication David Sankoff
    681 Gene and genome duplication David Sankoff Genomic sequencing projects have revealed the productivity of tetraploidization. I also summarize some mathematical processes duplicating genes or entire chromosome segments. modeling and algorithmics inspired by duplication phenomena. Substantial proportions of the yeast, Arabidopsis and human gene complements are made up of duplicates. This has prompted much Gene duplication interest in the processes of duplication, functional divergence and Li et al. [2] find that duplicated genes, as identified through loss of genes, has renewed the debate on whether an early fairly selective criteria, account for ~15% of the protein genes vertebrate genome was tetraploid, and has inspired mathematical in the human genome (counting both genes in each pair). In a models and algorithms in computational biology. survey of eukaryotic genome sequences, Lynch and Conery [3••], using a somewhat different filter, accounted for ~8%, Addresses 10% and 20% of the gene complement of the fly, yeast and Centre de recherches mathématiques, Université de Montréal, worm genomes, respectively. (Other estimates put the figure CP 6128 succursale Centre-Ville, Montreal, Québec H3C 3J7, Canada; at 16% for yeast and 25% for Arabidopsis [4•].) They estimated e-mail: [email protected] highly variable rates of gene duplication, averaging ~0.01 per Current Opinion in Genetics & Development 2001, 11:681–684 gene per Myr (million years). On the basis of ratios of silent and replacement rearrangements, they found that there is 0959-437X/01/$ — see front matter © 2001 Elsevier Science Ltd. All rights reserved. typically a period of neutral or (occasionally) even slightly accelerated evolution, lasting a few Myr at most, with one of Abbreviation the copies eventually being silenced in a large majority of Myr million years cases, and the remaining ones undergoing relatively stringent purifying selection.
    [Show full text]
  • Curriculum Vitae
    Curriculum Vitae Tandy Warnow Grainger Distinguished Chair in Engineering 1 Contact Information Department of Computer Science The University of Illinois at Urbana-Champaign Email: [email protected] Homepage: http://tandy.cs.illinois.edu 2 Research Interests Phylogenetic tree inference in biology and historical linguistics, multiple sequence alignment, metage- nomic analysis, big data, statistical inference, probabilistic analysis of algorithms, machine learning, combinatorial and graph-theoretic algorithms, and experimental performance studies of algorithms. 3 Professional Appointments • Co-chief scientist, C3.ai Digital Transformation Institute, 2020-present • Grainger Distinguished Chair in Engineering, 2020-present • Associate Head for Computer Science, 2019-present • Special advisor to the Head of the Department of Computer Science, 2016-present • Associate Head for the Department of Computer Science, 2017-2018. • Founder Professor of Computer Science, the University of Illinois at Urbana-Champaign, 2014- 2019 • Member, Carl R. Woese Institute for Genomic Biology. Affiliate of the National Center for Supercomputing Applications (NCSA), Coordinated Sciences Laboratory, and the Unit for Criticism and Interpretive Theory. Affiliate faculty member in the Departments of Mathe- matics, Electrical and Computer Engineering, Bioengineering, Statistics, Entomology, Plant Biology, and Evolution, Ecology, and Behavior, 2014-present. • National Science Foundation, Program Director for Big Data, July 2012-July 2013. • Member, Big Data Senior Steering Group of NITRD (The Networking and Information Tech- nology Research and Development Program), subcomittee of the National Technology Council (coordinating federal agencies), 2012-2013 • Departmental Scholar, Institute for Pure and Applied Mathematics, UCLA, Fall 2011 • Visiting Researcher, University of Maryland, Spring and Summer 2011. 1 • Visiting Researcher, Smithsonian Institute, Spring and Summer 2011. • Professeur Invit´e,Ecole Polytechnique F´ed´erale de Lausanne (EPFL), Summer 2010.
    [Show full text]
  • Structure-Based Realignment of Non-Coding Rnas in Multiple Whole Genome Alignments
    Structure-based Realignment of Non-coding RNAs in Multiple Whole Genome Alignments. by Michael Ku Yu Submitted to the Department of Electrical Engineering and Computer Science in partial fulfillment of the requirements for the degree of ARCHIVES Masters of Engineering in Computer Science and Engineering MASSACHUE N U TE at the OF TECH IOLOY MASSACHUSETTS INSTITUTE OF TECHNOLOGY JUN 2 1 2011 June 2011 LIBRARI ES @ Massachusetts Institute of Technology 2011. All rights reserved. '$7 A uthor ............ .. .. ... ............. Department of Electrical Wgineering and Computer Science May 20, 2011 Certified by..................................... ...... Bonnie Berger Professor of Applied Mathematics and Computer Science Thesis Supervisor Accepted by.... ....................................... Christopher J. Terman Chairman, Department Committee on Graduate Theses 2 Structure-based Realignment of Non-coding RNAs in Multiple Whole Genome Alignments by Michael Ku Yu Submitted to the Department of Electrical Engineering and Computer Science on May 20, 2011, in partial fulfillment of the requirements for the degree of Masters of Engineering in Computer Science and Engineering Abstract Whole genome alignments have become a central tool in biological sequence analy- sis. A major application is the de novo prediction of non-coding RNAs (ncRNAs) from structural conservation visible in the alignment. However, current methods for constructing genome alignments do so by explicitly optimizing for sequence simi- larity but not structural similarity. Therefore, de novo prediction of ncRNAs with high structural but low sequence conservation is intrinsically challenging in a genome alignment because the conservation signal is typically hidden. This study addresses this problem with a method for genome-wide realignment of potential ncRNAs ac- cording to structural similarity.
    [Show full text]
  • Algorithms for High-Throughput Sequencing Data
    Univerzita Komenského v Bratislave Fakulta matematiky, fyziky a informatiky Vladimír Boºa Autoreferát dizerta£nej práce Algorithms for high-throughput sequencing data na získanie akademického titulu philosophiae doctor v odbore doktorandského ²túdia: 9.2.1. - informatika Bratislava 2017 Dizerta£ná práca bola vypracovaná v dennej forme doktorandského ²túdia na Katedre ap- likovanej informatiky Fakulty matematiky, fyziky a informatiky Univerzity Komenského v Bratislave. Predkladate©: Vladimír Boºa Katedra aplikovanej informatiky Fakulta matematiky, fyziky a informatiky Univerzita Komenského Bratislava kolite©: Mgr. Tomá² Vina°, PhD. Katedra aplikovanej informatiky Fakulta matematiky, fyziky a informatiky Univerzita Komenského Bratislava Odbor: 9.2.1. - informatika Predseda odborovej komisie: prof. RNDr. Rastislav Krá©ovi£, PhD. Katedra informatiky Fakulta matematiky, fyziky a informatiky Univerzita Komenského Bratislava 1 Overview Knowledge of DNA sequences has become indispensable for basic biological research and in numerous applied elds like diagnostics, forensic biology, etc. It is also important for understanding cancer, ghting antibiotic resistant bacteria, etc. With decreasing cost of the DNA sequencing we are able to sequence more organisms, but consequently we have to handle much bigger amounts of data. Usual length of DNA sequence of an organism is between millions and billions of characters. For example, human DNA has approximatelly three billions bases. Unfortunately, we are still not able to read the whole DNA sequence at once, we can only read it in small pieces called reads. Usual length of reads vary between hundred to tens of thousands of bases depending on the sequencing technology. Thus, our only option is to try to computationaly reconstruct original sequence, using overlap information between reads. This process is called sequence assembly.
    [Show full text]
  • ISBRA 2012 Short Abstracts
    1 ISBRA 20 2 SHORT ABSTRACTS 8TH INTERNATIONAL SYMPOSIUM ON BIOINFORMATICS RESEARCH AND APPLICATIONS May 21-23, 2012 University of Texas at Dallas, Dallas, TX http://www.cs.gsu.edu/isbra12/ Symposium Organizers Steering Committee Dan Gusfield, University of California, Davis Ion Mandoiu, University of Connecticutt Yi Pan, Georgia State University Marie-France Sagot, INRIA Alex Zelikovsky, Georgia State University General Chairs Ovidiu Daesku, University of Texas at Dallas Raj Sunderraman, Georgia State University Program Chairs Leonidas Bleris, University of Texas at Dallas Ion Mandoiu, University of Connecticut Russell Schwartz, Carnegie Mellon University Jianxin Wang, Central South University Publicity Chair Sahar Al Seesi, University of Connecticut Finance Chairs Anu Bourgeois, Georgia State University Raj Sunderraman, Georgia State University Web Master, Web Design Piyaphol Phoungphol J. Steven Kirtzic Sponsors NATIONAL SCIENCE DEPARTMENT OF COMPUTER SCIENCE DEPARTMENT OF COMPUTER SCIENCE FOUNDATION GEORGIA STATE UNIVESITY UNIVERSITY OF TEXAS AT DALLAS i Program Committee Members Srinivas Aluru, Iowa State University Allen Holder, Rose-Hulman Istitute of S. Cenk Sahinalp, Simon Fraser Danny Barash, Ben-Gurion Technology University University Jinling Huang, Eastern Carolina David Sankoff, University of Ottawa Robert Beiko, Dalhousie University University Russell Schwartz, Carnegie Mellon Anne Bergeron, Universite du Lars Kaderali, University of University Quebec a Montreal Heidelberg Joao Setubal, Virginia Bioinformatics Iyad Kanj,
    [Show full text]
  • LANGUAGE CHANGE ACROSS the LIFESPAN: /R/ in MONTREAL FRENCH
    LANGUAGE CHANGE ACROSS THE LIFESPAN: /r/ IN MONTREAL FRENCH GILLIAN SANKOFF HE´ LE` NE BLONDEAU University of Pennsylvania University of Florida We address the articulation between language change in the historical sense and language change as experienced by individual speakers through a trend and panel study of the change from apical to dorsal /r/ in Montreal French. The community as a whole rapidly advanced its use of dorsal [R]. Most individual speakers followed across time were stable after the critical period, with phonological patterns set by the end of adolescence. A sizeable minority, however, made substantial changes. The window of opportunity for linguistic modification in later life may be expanded with rapid change in progress when linguistic variables take on social significance.* 1. INTRODUCTION. Reintegrating synchronic and diachronic linguistics by studying language change in progress, Weinreich, Labov, and Herzog (1968) demonstrated that change necessarily involves diffusion. They showed that language change implies varia- tion (though not the reverse), and that as an innovation diffuses through the speech community, the social structure of that community conditions its path. Subsequent research has not, however, directly addressed the question of the articulation between language change in the historical sense and language change as experienced by individ- ual speakers. We take this as our central concern in analyzing a phonological change in Montreal French—the community shift in the pronunciation of /r/—on the basis of data collected in 1971 and 1984. We address three major goals in our exploration of the problem. First, we look at the two major types of longitudinal research: trend and panel studies, and evaluate the relative contribution of each.
    [Show full text]
  • Copyright by Siavash Mir Arabbaygi 2015 the Dissertation Committee for Siavash Mir Arabbaygi Certifies That This Is the Approved Version of the Following Dissertation
    Copyright by Siavash Mir arabbaygi 2015 The Dissertation Committee for Siavash Mir arabbaygi certifies that this is the approved version of the following dissertation: Novel scalable approaches for multiple sequence alignment and phylogenomic reconstruction Committee: Keshav Pingali, Supervisor Tandy Warnow, Co-Supervisor David Hillis Bonnie Berger Joydeep Ghosh Ray Mooney Novel scalable approaches for multiple sequence alignment and phylogenomic reconstruction by Siavash Mir arabbaygi, B.S.; M. APPL S. DISSERTATION Presented to the Faculty of the Graduate School of The University of Texas at Austin in Partial Fulfillment of the Requirements for the Degree of DOCTOR OF PHILOSOPHY THE UNIVERSITY OF TEXAS AT AUSTIN August 2015 Dedicated to my mother, and the memory of my father. Acknowledgments I wish to thank the countless people who helped me throughout my PhD. I enjoyed working with many collaborators, most of whom I never met, and not all of whom I can name here. I had the good fortune of working with my labmates, especially Nam Nguyen, Shamsuzzoha Md. Bayzid, and Th´eo Zimmermann, and our interactions have always been smooth and fruitful. I benefited tremendously from my involvement in both the avian phylogenomics and the thousand transcripts projects, both of which included many top re- searchers from across the world. I want to especially thank Bastien Boussau, Edward Braun, Tom Gilbert, Erich Jarvis, Jim Leebens-Mack, Gane Wong, Norman Wickett, and Guojie Zhang. I also thank Kevin Liu, Mark Holder, and Jeet Sukumaran for allowing me to use their code throughout my work. I would also like to thank all members of my committee.
    [Show full text]
  • ISCB Honors Michael S. Waterman and Mathieu Blanchette Merry Maisel Ach Year, the International Fu¨ R Informatik and Chair of the ISCB Early 1970S
    Message From ISCB ISCB Honors Michael S. Waterman and Mathieu Blanchette Merry Maisel ach year, the International fu¨ r Informatik and chair of the ISCB early 1970s. In a memoir published on Society for Computational Awards Committee. ‘‘So much of our his Web site, Waterman has written: ‘‘I E Biology (ISCB) takes work is based on methods for finding was an innocent mathematician until nominations for its two major awards. sequence homology that, if we weren’t the summer of 1974. It was then that I met Temple Ferris Smith and for two An awards committee, composed of a constantly citing it, it would amaze us that the methodology was first devised months was cooped up with him in an group of current and past directors of only 25 years ago by Mike Waterman and office at Los Alamos. that the Society along with previous Temple Smith. Since then, Waterman experience transformed my research, recipients, evaluates the nominations developed the dynamic programming my life, and perhaps my sanity.’’ and selects the winners. In 2006, the approach to RNA structure prediction, Smith, now director of the awards committee honors two began the combinatorial study of RNA BioMolecular Engineering Research outstanding scientists. secondary structure, improved the Center at Boston University, was also statistical tests incorporated into BLAST visiting Los Alamos from a small Senior Scientist Accomplishment and related tools, and worked on the university in Michigan. Award: Michael S. Waterman assembly problem for genomes past and At Los Alamos, their fellow scientists present. He has also made important and friends included Stanslaw Ulam, contributions to phylogeny, tree Nick Metropolis, Marc Kac, and Gian- comparison, motif searching, cryptogene Carlo Rota, all towering names in analysis, parametric alignment, gapped computational science, at a time when alignment, optical mapping, haplotype the lab was a hotbed of intellectual estimation, gene family evolution, and a ferment.
    [Show full text]
  • Comparative Genomics Via Phylogenetic Invariants for Jukes-Cantor Semigroups
    Canadian Mathematical Society Conference Proceedings Comparative Genomics Via Phylogenetic Invariants For Jukes-Cantor Semigroups David Sankoff and Mathieu Blanchette Dedicated to Professor Donald A. Dawson Abstract. We review the theory of invariants as it has been developed for comparing the DNA sequences of homologous genes from phylogenetically re- lated species, with particular attention to the semigroups used to model se- quence evolution. We also outline the computational theory of genome rear- rangements, including the optimization problems in calculating edit distances between genomes and the simpler notion of breakpoint distance. The combi- natorics of rearrangements, involving non-local changes in the relative order of genes in the genome, are more complex than the base substitutions responsi- ble for gene sequence evolution. Nevertheless we can construct models of gene order evolution through symmetry assumptions about disruptions of gene ad- jacencies. Based on the extended Jukes-Cantor semigroup that emerges from this modeling, we derive a complete set of linear phylogenetic invariants. We use these invariants to relate mitochondrial genomes from a number of an- imal phyla and compare the results to parsimony trees also based on gene adjacencies and to minimal breakpoint trees. 1. Invariants for models of sequence evolution (1) (1) (N) (N) Consider the aligned DNA sequences: X1 ··· Xn , ··· ,X1 ··· Xn , all of length n, representing N species whose history of evolutionary divergence, or phylogeny, is represented by a tree T with vertex set V and edge set E, as in Figure 1. The terminal vertices represent observed, or present-day, species. The (J) non-terminal vertices represent hypothetical ancestral species.
    [Show full text]