DOCSLIB.ORG

Analysis of Overlapping Reading Frames in Viral Genomes

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Analysis of Overlapping Reading Frames in Viral Genomes From the Institute of Biochemistry of the University of Lübeck Director: Prof. Dr. Rolf Hilgenfeld Analysis of overlapping reading frames in viral genomes DissertAtion in Fulfillment of the Requirements for the DoctorAl Degree of the University of Lübeck from the DepArtment of NAturAl Sciences Submitted by Aditi Shukla from RAnchi, IndiA Lübeck, 2015 1. First Referee: Prof. Dr. Rolf Hilgenfeld 2. Second Referee: Prof. Dr. Bernhard Haubold DAte of OrAl ExAminAtion: 20.10.2015 Approved for printing: Lübeck, on 20.10.2015 Acknowledgements I thank, Prof. Hilgenfeld with deep gratitude for providing me with this reseArch opportunity And making me A pArt of your esteemed Institute. I thAnk you for sharing your knowledge And ideAs, providing support And encourAgement throughout, and for An extensive review of this work. Your attention to detail towArds every project, every proposAl, every manuscript, Amazes me. I wish to imbibe this immaculAteness in my professionAl life. Prof. MArtinetz And Prof. MAmlouk for interesting discussions, useful critiques for this research work, And for helping to keep my progress on schedule. The friendly manAgement of the grAduAte school for your sustAined support throughout. Special thanks to Olga and Katja for organizing interesting workshops, seminArs And get-togethers. Achim for your technical AssistAnce; Mrs. Rosenfeld for making my initiAl ArrivAl in Germany very eAsy And pleAsAnt And Mrs. SchwAb for work behind the scenes! Raffaele, Monarin, Raspudin, Zhenggang, Linlin, Friedrich And Caroline for your kind fraternity; Jeroen and Ksenia for your suggestions And encourAgements whenever required; for being there as a support! Other stAff members at the institute for making the institute so interesting. UpAsAnA And SusmitA for your friendship. My fAmily for immense patience, undying love and constant support. Contents Abstract Declaration 1. Introduction ....................................................................................................................... 1 1.1 Overlapping reading frames ........................................................................................ 1 1.1.1 Types of overlap ....................................................................................................... 2 1.2 Direction and phase of overlapping reading frames ................................................. 3 1.2.1 Direction of overlap .................................................................................................. 3 1.2.2 Phase of overlap in overlapping reading frames ...................................................... 4 1.3 Selection pressure on overlapping reading frames .................................................. 5 1.3.1 Degeneracy of the genetic code .............................................................................. 5 1.3.2 Synonymous and non-synonymous substitutions .................................................... 7 1.3.3 Different kinds of evolutionary strategy adopted by overlapping genes ................... 7 1.4 Viruses with overlapping reading frames ................................................................... 8 1.5 Objective of this study ................................................................................................. 8 2. Biological Background .................................................................................................... 11 2.1 Accessory proteins of coronaviruses ....................................................................... 11 2.1.1 Alphacoronaviruses ................................................................................................ 13 2.1.2 Betacoronaviruses ................................................................................................. 16 2.1.3 Gammacoronaviruses ............................................................................................ 19 2.1.4 Deltacoronaviruses ................................................................................................ 19 2.1.5 SARS-coronavirus: a member of betacoronavirus lineage b; genome organization and overlapping proteins .............................................................. 21 2.2 Overlapping protein in Murine norovirus ................................................................. 27 2.3 Overlapping proteins in Hepatitis B virus ................................................................ 27 3. Special sequence properties of overlapping genes ...................................................... 29 3.1 Codon usage bias ....................................................................................................... 29 3.2 Methods ....................................................................................................................... 30 3.2.1 Relative Synonymous Codon Usage (RSCU) ........................................................ 30 3.2.2 Correlation analysis and codon usage ................................................................... 31 3.2.3 Datasets ................................................................................................................. 31 3.2.4 Comparison of codon usage of overlapping and non-overlapping gene sets ........ 33 3.3 Results ......................................................................................................................... 34 3.4 Discussion ................................................................................................................... 37 4. Effect of overlaps on the protein products .................................................................... 39 4.1 Intrinsically disordered proteins ............................................................................... 39 4.2 Disorder predictors ..................................................................................................... 40 4.3 Methods ....................................................................................................................... 41 4.3.1 Comparison of disorder content of overlapping proteins sets ................................ 41 4.4 Results ......................................................................................................................... 42 4.5 Discussion ................................................................................................................... 46 5. Prediction of RNA secondary structure at the site of initiation of alternative reading frames ...................................................................................................................... 49 5.1 Molecular mechanism for the initiation of the alternative reading frame in overlapping genes ............................................................................................................ 49 5.1.1 Leaky scanning ...................................................................................................... 50 5.1.2 Internal Ribosomal Entry ........................................................................................ 52 5.2 Methods ....................................................................................................................... 53 5.2.1 Vienna RNA package ............................................................................................. 53 5.2.2 RNAComposer ....................................................................................................... 54 5.2.3 Computer prediction of RNA secondary and tertiary structural elements .............. 54 5.3 Results ......................................................................................................................... 56 5.4 Discussion ................................................................................................................... 57 6. Mutational model for the evolution of SARS-CoV overlapping accessory protein 9b .............................................................................................................................. 59 6.1 Introduction ................................................................................................................. 59 6.2 Methods ....................................................................................................................... 59 6.2.1 Dataset ................................................................................................................... 59 6.2.2 Mutation rate analysis ............................................................................................ 60 6.2.3 Entropy-plotting ...................................................................................................... 60 6.3 Results ......................................................................................................................... 62 6.3.1 The effect of the overlap on the mutation rate in the N and orf9b genes ............... 62 6.3.2 Evolutionary strategy adopted by the overlapping N and orf9b gene set ............... 63 6.3.3 Effect of mutations on the three-dimensional structures of the overlapping proteins ........................................................................................................................... 67 6.4 Discussion ................................................................................................................... 68 7. Conclusions ...................................................................................................................... 70 Bibliography .........................................................................................................................
Load more

Recommended publications
  • Long Overlapping Genes in Pseudomonas Aeruginosa
    bioRxiv preprint doi: https://doi.org/10.1101/2021.02.09.430400; this version posted February 10, 2021. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-NC-ND 4.0 International license. 1 Shadow ORFs illuminated: long overlapping genes in Pseudomonas 2 aeruginosa are translated and under purifying selection 3 4 Michaela Kreitmeier1, Zachary Ardern1*, Miriam Abele2, Christina Ludwig2, Siegfried Scherer1, 5 Klaus Neuhaus3* 6 7 1Chair for Microbial Ecology, TUM School of Life Sciences, Technische Universität München 8 Weihenstephaner Berg 3, 85354 Freising, Germany. 9 2Bavarian Center for Biomolecular Mass Spectrometry (BayBioMS), TUM School of Life 10 Sciences, Technische Universität München, Gregor-Mendel-Strasse 4, 85354 Freising, 11 Germany. 12 3Core Facility Microbiome, ZIEL – Institute for Food & Health, Technische Universität München 13 Weihenstephaner Berg 3, 85354 Freising, Germany. 14 15 *email: [email protected]; [email protected] 16 17 Abstract 18 The existence of overlapping genes (OLGs) with significant coding overlaps revolutionises our 19 understanding of genomic complexity. We report two exceptionally long (957 nt and 1536 nt), 20 evolutionarily novel, translated antisense open reading frames (ORFs) embedded within 21 annotated genes in the medically important Gram-negative bacterium Pseudomonas 22 aeruginosa. Both OLG pairs show sequence features consistent with being genes and 23 transcriptional signals in RNA sequencing data. Translation of both OLGs was confirmed by 24 ribosome profiling and mass spectrometry. Quantitative proteomics of samples taken during 25 different phases of growth revealed regulation of protein abundances, implying biological 26 functionality.
    [Show full text]
  • A Novel Short L-Arginine Responsive Protein-Coding Gene (Laob)
    Erschienen in: BMC evolutionary biology ; 18 (2018), 1. - 21 https://dx.doi.org/10.1186/s12862-018-1134-0 Hücker et al. BMC Evolutionary Biology (2018) 18:21 https://doi.org/10.1186/s12862-018-1134-0 RESEARCH ARTICLE Open Access A novel short L-arginine responsive protein-coding gene (laoB) antiparallel overlapping to a CadC-like transcriptional regulator in Escherichia coli O157:H7 Sakai originated by overprinting Sarah M. Hücker1,5, Sonja Vanderhaeghen1, Isabel Abellan-Schneyder1,4, Romy Wecko1, Svenja Simon2, Siegfried Scherer1,3 and Klaus Neuhaus1,4* Abstract Background: Due to the DNA triplet code, it is possible that the sequences of two or more protein-coding genes overlap to a large degree. However, such non-trivial overlaps are usually excluded by genome annotation pipelines and, thus, only a few overlapping gene pairs have been described in bacteria. In contrast, transcriptome and translatome sequencing reveals many signals originated from the antisense strand of annotated genes, of which we analyzed an example gene pair in more detail. Results: A small open reading frame of Escherichia coli O157:H7 strain Sakai (EHEC), designated laoB (L-arginine responsive overlapping gene), is embedded in reading frame −2 in the antisense strand of ECs5115, encoding a CadC-like transcriptional regulator. This overlapping gene shows evidence of transcription and translation in Luria-Bertani (LB) and brain-heart infusion (BHI) medium based on RNA sequencing (RNAseq) and ribosomal- footprint sequencing (RIBOseq). The transcriptional start site is 289 base pairs (bp) upstream of the start codon and transcription termination is 155 bp downstream of the stop codon.
    [Show full text]
  • 1. a 6-Frame Translation Map of a Segment of DNA Is Shown, with Three Open Reading Frames (A, B, and C). Orfs a and B Are Known to Be in Separate Genes
    1. A 6-frame translation map of a segment of DNA is shown, with three open reading frames (A, B, and C). ORFs A and B are known to be in separate genes. 1a. Two transcription bubbles are shown, one in ORF A and one in ORF B. In the transcription bubble diagram, mark the following: • the location of RNA polymerase on the appropriate strand in each bubble • the RNA transcripts to show the relative lengths of RNA made by those two polymerases 1b. Are the promoters for ORFs A and/or B present in the DNA region shown in this diagram? Promoter for A: Present? Yes No (circle one) If present, mark its location and label it. Promoter for B: Present? Yes No (circle one) If present, mark its location and label it. 1c. Electron microscopy experiments failed to show RNA polymerases over the ORF "C" region of DNA. State whether each of the three explanations listed below is valid or not, explaining as necessary: Explanation If valid, just write “Valid.” If invalid, BRIEFLY explain why. _______________________________________________________________________ ORFs B and C are exons of the same gene and a splicing error causes ORF C to be left out. Splicing occurs after transcription. Incorrect splicing doesn't explain why transcription didn't happen. ORF C has a mutation in its start (ATG) codon, preventing transcription. The start codon is used for translation, not transcription... whether or not the start codon is intact, transcription could still happen. ORF C has a promoter mutation preventing transcription. VALID. (A promoter mutation is consistent with failure to transcribe the gene.) 2.
    [Show full text]
  • Insights Into the Codon Usage Bias of 13 Severe Acute
    bioRxiv preprint doi: https://doi.org/10.1101/2020.04.01.019463; this version posted April 4, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-NC-ND 4.0 International license. Insights into The Codon Usage Bias of 13 Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Isolates from Different Geo- locations Ali Mostafa Anwar *1, Saif M. Khodary 1 1 Department of Genetics, Faculty of Agriculture, Cairo University, Giza, 12613, Egypt *Correspondence: [email protected] Abstract Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of Coronavirus disease 2019 (COVID-19) which is an infectious disease that spread throughout the world and was declared as a pandemic by the World Health Organization (WHO). In the present study, we analyzed genome-wide codon usage patterns in 13 SARS-CoV-2 isolates from different geo-locations (countries) by utilizing different CUB measurements. Nucleotide and di-nucleotide compositions displayed bias toward A/U content in all codon positions and CpU-ended codons preference, respectively. Relative Synonymous Codon Usage (RSCU) analysis revealed 8 common putative preferred codons among all the 13 isolates. Interestingly, all of the latter codons are A/U-ended (U-ended: 7, A-ended: 1). Cluster analysis (based on RSCU values) was performed and showed comparable results to the phylogenetic analysis (based on their whole genome sequences) indicating that the CUB pattern may reflect the evolutionary relationship between the tested isolates.
    [Show full text]
  • Codon Usage and Adenovirus Fitness: Implications for Vaccine Development
    fmicb-12-633946 February 8, 2021 Time: 11:48 # 1 REVIEW published: 10 February 2021 doi: 10.3389/fmicb.2021.633946 Codon Usage and Adenovirus Fitness: Implications for Vaccine Development Judit Giménez-Roig1, Estela Núñez-Manchón1, Ramon Alemany2, Eneko Villanueva3 and Cristina Fillat1,4,5* 1 Institut d’Investigacions Biomèdiques August Pi i Sunyer (IDIBAPS), Barcelona, Spain, 2 Procure Program, Institut Català d’Oncologia- Oncobell Program, IDIBELL, L’Hospitalet de Llobregat, Barcelona, Spain, 3 Cambridge Centre for Proteomics, Department of Biochemistry, University of Cambridge, Cambridge, United Kingdom, 4 Centro de Investigación Biomédica en Red de Enfermedades Raras (CIBERER), Barcelona, Spain, 5 Facultat de Medicina i Ciències de la Salut, Universitat de Barcelona (UB), Barcelona, Spain Vaccination is the most effective method to date to prevent viral diseases. It intends to mimic a naturally occurring infection while avoiding the disease, exposing our bodies to viral antigens to trigger an immune response that will protect us from future infections. Among different strategies for vaccine development, recombinant vaccines are one of the most efficient ones. Recombinant vaccines use safe viral vectors as vehicles and incorporate a transgenic antigen of the pathogen against which we intend to Edited by: Rosa Maria Pintó, generate an immune response. These vaccines can be based on replication-deficient University of Barcelona, Spain viruses or replication-competent viruses. While the most effective strategy involves Reviewed by: replication-competent viruses, they must be attenuated to prevent any health hazard Kai Li, Harbin Veterinary Research Institute, while guaranteeing a strong humoral and cellular immune response. Several attenuation Chinese Academy of Agricultural strategies for adenoviral-based vaccine development have been contemplated over Sciences, China time.
    [Show full text]
  • Dissertation
    Dissertation submitted to the Combined Faculty of Natural Sciences and Mathematics of the Ruperto-Carola University of Heidelberg, Germany for the degree of Doctor of Natural Sciences Presented by M.Sc. Matilde Bertolini born in Parma, Italy Oral examination: 02/03/2021 Profiling interactions of proximal nascent chains reveals a general co-translational mechanism of protein complex assembly Referees: Prof. Dr. Bernd Bukau Prof. Dr. Caludio Joazeiro Preface Contributions The experiments and analyses presented in this Thesis were conducted by myself unless otherwise indicated, under the supervision of Dr. Günter Kramer and Prof. Dr. Bernd Bukau. The Disome Selective Profiling (DiSP) technology was developed in collaboration with my colleague Kai Fenzl, with whom I also generated initial datasets of human HEK293-T and U2OS cells. Dr. Ilia Kats developed important bioinformatics tools for the analysis of DiSP data, including RiboSeqTools and the sigmoid fitting algorithm. He also offered great input on statistical analyses. Dr. Frank Tippmann performed analysis of crystal structures. Table of Contents TABLE OF CONTENTS List of Figures ......................................................................................................V List of Tables ...................................................................................................... VII List of Equations .............................................................................................. VIII Abbreviations ....................................................................................................
    [Show full text]
  • The Design of Maximally Compressed Coding Sequences
    Two Proteins for the Price of One: The Design of Maximally Compressed Coding Sequences Bei Wang1, Dimitris Papamichail2, Steffen Mueller3, and Steven Skiena2 1 Dept. of Computer Science, Duke University, Durham, NC 27708 [email protected] 2 Dept. of Computer Science, State University of New York, Stony Brook, NY 11794 {dimitris, skiena}@cs.sunysb.edu 3 Dept. of Microbiology, State University of New York, Stony Brook, NY 11794 [email protected] Abstract. The emerging field of synthetic biology moves beyond con- ventional genetic manipulation to construct novel life forms which do not originate in nature. We explore the problem of designing the provably shortest genomic sequence to encode a given set of genes by exploit- ing alternate reading frames. We present an algorithm for designing the shortest DNA sequence simultaneously encoding two given amino acid sequences. We show that the coding sequence of naturally occurring pairs of overlapping genes approach maximum compression. We also investi- gate the impact of alternate coding matrices on overlapping sequence de- sign. Finally, we discuss an interesting application for overlapping gene design, namely the interleaving of an antibiotic resistance gene into a target gene inserted into a virus or plasmid for amplification. 1 Introduction The emerging field of synthetic biology moves beyond conventional genetic ma- nipulation to construct novel life forms which do not originate in nature. The synthesis of poliovirus from off-to-shelf components [1] attracted worldwide at- tention when announced in July 2002. Subsequently, the bacteriophage PhiX174 was synthesized using different techniques in only three weeks [2], and Kodu- mal, et al.
    [Show full text]
  • Phylogeny Recapitulates Learning: Self-Optimization of Genetic Code
    bioRxiv preprint doi: https://doi.org/10.1101/260877; this version posted February 6, 2018. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-NC 4.0 International license. Phylogeny Recapitulates Learning: Self-Optimization of Genetic Code Oliver Attie1, Brian Sulkow1, Chong Di1, and Wei-Gang Qiu1,2,3,* 1Department of Biological Sciences, Hunter College & 2Graduate Center, City University of New York; 3Department of Physiology and Biophysics & Institute for Computational Biomedicine, Weil Cornell Medical College, New York, New York, USA Emails: Oliver Attie [email protected]; Brian Sulkow [email protected]; Chong Di [email protected]; *Correspondence: [email protected] Abstract Learning algorithms have been proposed as a non-selective mechanism capable of creating complex adaptive systems in life. Evolutionary learning however has not been demonstrated to be a plausible cause for the origin of a specific molecular system. Here we show that genetic codes as optimal as the Standard Genetic Code (SGC) emerge readily by following a molecular analog of the Hebb’s rule (“neurons fire together, wire together”). Specifically, error-minimizing genetic codes are obtained by maximizing the number of physio-chemically similar amino acids assigned to evolutionarily similar codons. Formulating genetic code as a Traveling Salesman Problem (TSP) with amino acids as “cities” and codons as “tour positions” and implemented with a Hopfield neural network, the unsupervised learning algorithm efficiently finds an abundance of genetic codes that are more error- minimizing than SGC.
    [Show full text]
  • Codon Clusters with Biased Synonymous Codon Usage Represent
    bioRxiv preprint doi: https://doi.org/10.1101/530345; this version posted January 25, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission. 1 Codon clusters with biased synonymous codon usage represent 2 hidden functional domains in protein-coding DNA sequences 3 4 Authors 5 Zhen Peng1*, Yehuda Ben-Shahar1 6 7 Affiliations 8 1Department of Biology, Washington University in St. Louis, MO 63130, USA. 9 *Correspondence to: Zhen Peng ([email protected]). 10 11 Outline 12 1. Abstract 13 2. Introduction 14 3. Results 15 3.1. Identifying putatively functional codon clusters (PFCCs) 16 3.2. Codon usage patterns of PFCCs are diverse 17 3.3. PFCC distribution is not restricted to specific regions of protein-coding sequences 18 3.4. Specific protein functional classes are overrepresented in genes carrying PFCCs while most PFCCs 19 are not associated with known protein domains 20 3.5. Voltage-gated sodium channels include a conserved rare-codon cluster associated with the 21 inactivation gate 22 4. Discussion 23 5. Materials and Methods 24 5.1. Reference genome and data pre-procession Page 1 of 28 bioRxiv preprint doi: https://doi.org/10.1101/530345; this version posted January 25, 2019. The copyright holder for this preprint (which was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission. 25 5.2. Identifying PFCCs 26 5.3. Calculating TCAI 27 5.4. K-mean clustering of PFCCs 28 5.5.
    [Show full text]
  • Classification and Function of Small Open-Reading Frames Abstract
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by Sussex Research Online 1 Classification and function of small open-reading frames Juan-Pablo Couso1,2* and Pedro Patraquim2 1Centro Andaluz de Biologia del Desarrollo, CSIC-UPO, Sevilla, Spain and 2Brighton and Sussex Medical School, University of Sussex, Brighton, United Kingdom. *Author for correspondence: [email protected] Abstract Small open-reading frames (smORFs or sORFs) of 100 codons or less are usually - if arbitrarily - excluded from canonical proteome annotations. Despite this, the genomes of a wide range of metazoans, including humans, contain hundreds of smORFs, some of which fulfil key physiological functions. Recently, ribosomal profiling has been employed to show that the transcriptome of the model organism Drosophila melanogaster contains thousands of smORFs of different classes actively undergoing translation which produces peptides of mostly unknown function. Here we present a comprehensive analysis of the smORF repertoire in flies, mice and humans. We propose the existence of several classes of smORFs with different functions, from inert DNA sequences to transcribed and translated cis- regulators of translation, and finally to expression of functional peptides with a propensity to act as regulators of canonical membrane-associated proteins, or as components of ancestral protein complexes in the cytoplasm. We suggest that the different smORF classes could represent steps during the evolution of novel peptide and protein sequences. Our analysis introduces a distinction between different peptide-coding classes in animal genomes, and highlights the role of Drosophila melanogaster as a model organism for the study of small peptide biology in the context of development, physiology and human disease.
    [Show full text]
  • Effect of Correlated TRAN Abundances on Translation Errors and Evolution of Codon Usage Bias
    University of Tennessee, Knoxville TRACE: Tennessee Research and Creative Exchange Faculty Publications and Other Works -- Ecology and Evolutionary Biology Ecology and Evolutionary Biology 9-2010 Effect of Correlated TRAN Abundances on Translation Errors and Evolution of Codon Usage Bias Premal Shah University of Tennessee - Knoxville Michael A. Gilchrist Follow this and additional works at: https://trace.tennessee.edu/utk_ecolpubs Part of the Ecology and Evolutionary Biology Commons Recommended Citation Shah, Premal and Gilchrist, Michael A., "Effect of Correlated TRAN Abundances on Translation Errors and Evolution of Codon Usage Bias" (2010). Faculty Publications and Other Works -- Ecology and Evolutionary Biology. https://trace.tennessee.edu/utk_ecolpubs/37 This Article is brought to you for free and open access by the Ecology and Evolutionary Biology at TRACE: Tennessee Research and Creative Exchange. It has been accepted for inclusion in Faculty Publications and Other Works -- Ecology and Evolutionary Biology by an authorized administrator of TRACE: Tennessee Research and Creative Exchange. For more information, please contact [email protected]. Effect of Correlated tRNA Abundances on Translation Errors and Evolution of Codon Usage Bias Premal Shah1,2*, Michael A. Gilchrist1,2 1 Department of Ecology & Evolutionary Biology, University of Tennessee, Knoxville, Tennessee, United States of America, 2 National Institute for Mathematical and Biological Synthesis, University of Tennessee, Knoxville, Tennessee, United States of America Abstract Despite the fact that tRNA abundances are thought to play a major role in determining translation error rates, their distribution across the genetic code and the resulting implications have received little attention. In general, studies of codon usage bias (CUB) assume that codons with higher tRNA abundance have lower missense error rates.
    [Show full text]
  • Analysis of Codon Usage Patterns in Giardia Duodenalis Based on Transcriptome Data from Giardiadb
    G C A T T A C G G C A T genes Article Analysis of Codon Usage Patterns in Giardia duodenalis Based on Transcriptome Data from GiardiaDB Xin Li, Xiaocen Wang, Pengtao Gong, Nan Zhang, Xichen Zhang and Jianhua Li * Key Laboratory of Zoonosis Research, Ministry of Education, College of Veterinary Medicine, Jilin University, Changchun 130062, China; [email protected] (X.L.); [email protected] (X.W.); [email protected] (P.G.); [email protected] (N.Z.); [email protected] (X.Z.) * Correspondence: [email protected]; Tel.: +86-431-8783-6172; Fax: +86-431-8798-1351 Abstract: Giardia duodenalis, a flagellated parasitic protozoan, the most common cause of parasite- induced diarrheal diseases worldwide. Codon usage bias (CUB) is an important evolutionary character in most species. However, G. duodenalis CUB remains unclear. Thus, this study analyzes codon usage patterns to assess the restriction factors and obtain useful information in shaping G. duo- denalis CUB. The neutrality analysis result indicates that G. duodenalis has a wide GC3 distribution, which significantly correlates with GC12. ENC-plot result—suggesting that most genes were close to the expected curve with only a few strayed away points. This indicates that mutational pressure and natural selection played an important role in the development of CUB. The Parity Rule 2 plot (PR2) result demonstrates that the usage of GC and AT was out of proportion. Interestingly, we identified 26 optimal codons in the G. duodenalis genome, ending with G or C. In addition, GC content, gene expression, and protein size also influence G.
    [Show full text]