1
|
Kunisch F, Campobasso C, Wagemans J, Yildirim S, Chan BK, Schaudinn C, Lavigne R, Turner PE, Raschke MJ, Trampuz A, Gonzalez Moreno M. Targeting Pseudomonas aeruginosa biofilm with an evolutionary trained bacteriophage cocktail exploiting phage resistance trade-offs. Nat Commun 2024; 15:8572. [PMID: 39362854 PMCID: PMC11450229 DOI: 10.1038/s41467-024-52595-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/29/2023] [Accepted: 09/12/2024] [Indexed: 10/05/2024] Open
Abstract
Spread of multidrug-resistant Pseudomonas aeruginosa strains threatens to render currently available antibiotics obsolete, with limited prospects for the development of new antibiotics. Lytic bacteriophages, the viruses of bacteria, represent a path to combat this threat. In vitro-directed evolution is traditionally applied to expand the bacteriophage host range or increase bacterial suppression in planktonic cultures. However, while up to 80% of human microbial infections are biofilm-associated, research towards targeted improvement of bacteriophages' ability to combat biofilms remains scarce. This study aims at an in vitro biofilm evolution assay to improve multiple bacteriophage parameters in parallel and the optimisation of bacteriophage cocktail design by exploiting a bacterial bacteriophage resistance trade-off. The evolved bacteriophages show an expanded host spectrum, improved antimicrobial efficacy and enhanced antibiofilm performance, as assessed by isothermal microcalorimetry and quantitative polymerase chain reaction, respectively. Our two-phage cocktail reveals further improved antimicrobial efficacy without incurring dual-bacteriophage-resistance in treated bacteria. We anticipate this assay will allow a better understanding of phenotypic-genomic relationships in bacteriophages and enable the training of bacteriophages against other desired pathogens. This, in turn, will strengthen bacteriophage therapy as a treatment adjunct to improve clinical outcomes of multidrug-resistant bacterial infections.
Collapse
Affiliation(s)
- Fabian Kunisch
- Faculty of Medicine, Universität Münster, Münster, Germany
- Center for Musculoskeletal Surgery, Charité - Universitätsmedizin Berlin, Corporate Member of Freie Universität Berlin and Humboldt-Universität zu Berlin, Berlin, Germany
- Department of Ecology and Evolutionary Biology, Yale University, New Haven, CT, USA
- Center for Phage Biology and Therapy, Yale University, New Haven, CT, USA
| | - Claudia Campobasso
- Department of Biosystems, KU Leuven, Leuven, Belgium
- Department of Biology, Università di Pisa, Pisa, Italy
| | | | - Selma Yildirim
- Berlin Institute of Health at Charité - Universitätsmedizin Berlin, BIH Center for Regenerative Therapies (BCRT), Berlin, Germany
| | - Benjamin K Chan
- Department of Ecology and Evolutionary Biology, Yale University, New Haven, CT, USA
- Center for Phage Biology and Therapy, Yale University, New Haven, CT, USA
| | - Christoph Schaudinn
- Advanced Light and Electron Microscopy (Zentrum für Biologische Gefahren und Spezielle Pathogene 4), Robert Koch Institute, Berlin, Germany
| | - Rob Lavigne
- Department of Biosystems, KU Leuven, Leuven, Belgium
| | - Paul E Turner
- Department of Ecology and Evolutionary Biology, Yale University, New Haven, CT, USA
- Center for Phage Biology and Therapy, Yale University, New Haven, CT, USA
- Program in Microbiology, Yale School of Medicine, New Haven, CT, USA
| | - Michael J Raschke
- Faculty of Medicine, Universität Münster, Münster, Germany
- Department of Trauma, Hand and Reconstructive Surgery, Universitätsklinikum Münster, Münster, Germany
| | - Andrej Trampuz
- Center for Musculoskeletal Surgery, Charité - Universitätsmedizin Berlin, Corporate Member of Freie Universität Berlin and Humboldt-Universität zu Berlin, Berlin, Germany.
- Berlin Institute of Health at Charité - Universitätsmedizin Berlin, BIH Center for Regenerative Therapies (BCRT), Berlin, Germany.
| | - Mercedes Gonzalez Moreno
- Center for Musculoskeletal Surgery, Charité - Universitätsmedizin Berlin, Corporate Member of Freie Universität Berlin and Humboldt-Universität zu Berlin, Berlin, Germany
- Berlin Institute of Health at Charité - Universitätsmedizin Berlin, BIH Center for Regenerative Therapies (BCRT), Berlin, Germany
| |
Collapse
|
2
|
Perez LJ, Cloherty GA, Berg MG. Parallel evolution of picobirnaviruses from distinct ancestral origins. Microbiol Spectr 2023; 11:e0269323. [PMID: 37888988 PMCID: PMC10714727 DOI: 10.1128/spectrum.02693-23] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2023] [Accepted: 09/19/2023] [Indexed: 10/28/2023] Open
Abstract
IMPORTANCE Picobirnaviruses (PBVs) are highly heterogeneous viruses encoding a capsid and RdRp. Detected in a wide variety of animals with and without disease, their association with gastrointestinal and respiratory infections, and consequently their public health importance, has rightly been questioned. Determining the "true" host of Picobirnavirus lies at the center of this debate, as evidence exists for them having both vertebrate and prokaryotic origins. Using integrated and time-stamped phylogenetic approaches, we show they are contemporaneous viruses descending from two different ancestors: avian Reovirus and fungal Partitivirus. The fungal PBV-R2 species emerged with a single segment (RdRp) until it acquired a capsid from vertebrate PBV-R1 and PBV-R3 species. Protein and RNA folding analyses revealed how the former came to resemble the latter over time. Thus, parallel evolution from disparate hosts has driven the adaptation and genetic diversification of the Picobirnaviridae family.
Collapse
Affiliation(s)
- Lester J. Perez
- Infectious Disease Core Research, Abbott Diagnostics Division, Abbott Laboratories, Abbott Park, Illinois, USA
- Abbott Pandemic Defense Coalition (APDC), Chicago, Illinois, USA
| | - Gavin A. Cloherty
- Infectious Disease Core Research, Abbott Diagnostics Division, Abbott Laboratories, Abbott Park, Illinois, USA
- Abbott Pandemic Defense Coalition (APDC), Chicago, Illinois, USA
| | - Michael G. Berg
- Infectious Disease Core Research, Abbott Diagnostics Division, Abbott Laboratories, Abbott Park, Illinois, USA
- Abbott Pandemic Defense Coalition (APDC), Chicago, Illinois, USA
| |
Collapse
|
3
|
Horton JS, Taylor TB. Mutation bias and adaptation in bacteria. MICROBIOLOGY (READING, ENGLAND) 2023; 169:001404. [PMID: 37943288 PMCID: PMC10710837 DOI: 10.1099/mic.0.001404] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/02/2023] [Accepted: 10/11/2023] [Indexed: 11/10/2023]
Abstract
Genetic mutation, which provides the raw material for evolutionary adaptation, is largely a stochastic force. However, there is ample evidence showing that mutations can also exhibit strong biases, with some mutation types and certain genomic positions mutating more often than others. It is becoming increasingly clear that mutational bias can play a role in determining adaptive outcomes in bacteria in both the laboratory and the clinic. As such, understanding the causes and consequences of mutation bias can help microbiologists to anticipate and predict adaptive outcomes. In this review, we provide an overview of the mechanisms and features of the bacterial genome that cause mutational biases to occur. We then describe the environmental triggers that drive these mechanisms to be more potent and outline the adaptive scenarios where mutation bias can synergize with natural selection to define evolutionary outcomes. We conclude by describing how understanding mutagenic genomic features can help microbiologists predict areas sensitive to mutational bias, and finish by outlining future work that will help us achieve more accurate evolutionary forecasts.
Collapse
Affiliation(s)
- James S. Horton
- Milner Centre for Evolution, Department of Life Sciences, University of Bath, BA2 7AY, UK
| | - Tiffany B. Taylor
- Milner Centre for Evolution, Department of Life Sciences, University of Bath, BA2 7AY, UK
| |
Collapse
|
4
|
Gitschlag BL, Cano AV, Payne JL, McCandlish DM, Stoltzfus A. Mutation and Selection Induce Correlations between Selection Coefficients and Mutation Rates. Am Nat 2023; 202:534-557. [PMID: 37792926 DOI: 10.1086/726014] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/06/2023]
Abstract
AbstractThe joint distribution of selection coefficients and mutation rates is a key determinant of the genetic architecture of molecular adaptation. Three different distributions are of immediate interest: (1) the "nominal" distribution of possible changes, prior to mutation or selection; (2) the "de novo" distribution of realized mutations; and (3) the "fixed" distribution of selectively established mutations. Here, we formally characterize the relationships between these joint distributions under the strong-selection/weak-mutation (SSWM) regime. The de novo distribution is enriched relative to the nominal distribution for the highest rate mutations, and the fixed distribution is further enriched for the most highly beneficial mutations. Whereas mutation rates and selection coefficients are often assumed to be uncorrelated, we show that even with no correlation in the nominal distribution, the resulting de novo and fixed distributions can have correlations with any combination of signs. Nonetheless, we suggest that natural systems with a finite number of beneficial mutations will frequently have the kind of nominal distribution that induces negative correlations in the fixed distribution. We apply our mathematical framework, along with population simulations, to explore joint distributions of selection coefficients and mutation rates from deep mutational scanning and cancer informatics. Finally, we consider the evolutionary implications of these joint distributions together with two additional joint distributions relevant to parallelism and the rate of adaptation.
Collapse
|
5
|
Arefinia N, Yaghobi R, Ramezani A, Sarvari J. Sequence Analysis of Hot Spot Regions of Spike and RNA-dependent-RNA polymerase (RdRp) Genes of SARS-CoV-2 in Kerman, Iran. Mediterr J Hematol Infect Dis 2023; 15:e2023042. [PMID: 37435034 PMCID: PMC10332355 DOI: 10.4084/mjhid.2023.042] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/26/2023] [Accepted: 06/18/2023] [Indexed: 07/13/2023] Open
Abstract
Background Mutations in the SARS-CoV-2 genome might influence pathogenicity, transmission rate, and evasion of the host immune system. Therefore, the purpose of the present study was to investigate the genetic alteration as well as assess their effects on the receptor binding domain (RBD) of the spike and the putative RNA binding site of the RdRp genes of SARS-CoV-2 using bioinformatics tools. Materials and Method In this cross-sectional study, 45 confirmed COVID-19 patients using qRT-PCR were included and divided into mild, severe, and critical groups based on the severity of the disease. RNA was extracted from nasopharyngeal swab samples using a commercial kit. RT-PCR was performed to amplify the target sequences of the spike and RdRp genes and sequence them by the Sanger method. Clustal OMEGA, MEGA 11 software, I-mutant tools, SWISS-MODEL, and HDOCK web servers were used for bioinformatics analyses. Results The mean age of the patients was 50.68±2.73. The results showed that four of six mutations (L452R, T478K, N501Y, and D614G) in RBD and three of eight in the putative RNA binding site (P314L, E1084D, V1883T) were missense. In the putative RNA binding site, another deletion was discovered. Among missense mutations, N501Y and V1883T were responsible for increasing structural stability, while others were responsible for decreasing it. The various homology models designed showed that these homologies were like the Wuhan model. The molecular docking analysis revealed that the T478K mutation in RBD had the highest binding affinity. In addition, 35 RBD samples (89.7%) and 33 putative RNA binding site samples (84.6%) were similar to the Delta variant. Conclusion Our results indicated that double mutations (T478K and N501Y) in the S protein might increase the binding affinity of SARS-CoV-2 to human ACE2 compared to the wild-type (WT) strain. Moreover, variations in the spike and RdRp genes might influence the stability of encoded proteins.
Collapse
Affiliation(s)
- Nasir Arefinia
- Department of Bacteriology and Virology, School of Medicine, Shiraz University of Medical Sciences, Shiraz, Iran
- School of Medicine, Jiroft University of Medical Sciences, Jiroft, Iran
| | - Ramin Yaghobi
- Shiraz Transplant Research Center, Shiraz University of Medical Sciences, Shiraz, Iran
| | - Amin Ramezani
- Department of Medical Biotechnology, School of Advanced Medical Sciences and Technologies, Shiraz University of Medical Sciences, Shiraz, Iran
- Shiraz Institute for Cancer Research, School of Medicine, Shiraz University of Medical Sciences, Shiraz, Iran
| | - Jamal Sarvari
- Department of Bacteriology and Virology, School of Medicine, Shiraz University of Medical Sciences, Shiraz, Iran
- Gastroenterohepatology Research Center, Shiraz University of Medical Sciences, Shiraz, Iran
| |
Collapse
|
6
|
Cano AV, Gitschlag BL, Rozhoňová H, Stoltzfus A, McCandlish DM, Payne JL. Mutation bias and the predictability of evolution. Philos Trans R Soc Lond B Biol Sci 2023; 378:20220055. [PMID: 37004719 PMCID: PMC10067271 DOI: 10.1098/rstb.2022.0055] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2022] [Accepted: 02/16/2023] [Indexed: 04/04/2023] Open
Abstract
Predicting evolutionary outcomes is an important research goal in a diversity of contexts. The focus of evolutionary forecasting is usually on adaptive processes, and efforts to improve prediction typically focus on selection. However, adaptive processes often rely on new mutations, which can be strongly influenced by predictable biases in mutation. Here, we provide an overview of existing theory and evidence for such mutation-biased adaptation and consider the implications of these results for the problem of prediction, in regard to topics such as the evolution of infectious diseases, resistance to biochemical agents, as well as cancer and other kinds of somatic evolution. We argue that empirical knowledge of mutational biases is likely to improve in the near future, and that this knowledge is readily applicable to the challenges of short-term prediction. This article is part of the theme issue 'Interdisciplinary approaches to predicting evolutionary biology'.
Collapse
Affiliation(s)
- Alejandro V. Cano
- Institute of Integrative Biology, ETH Zurich, 8092 Zurich, Switzerland
- Swiss Institute of Bioinformatics, 1015 Lausanne, Switzerland
| | - Bryan L. Gitschlag
- Simons Center for Quantitative Biology, Cold Spring Harbor Laboratory, Cold Spring Harbor, NY 11724, USA
| | - Hana Rozhoňová
- Institute of Integrative Biology, ETH Zurich, 8092 Zurich, Switzerland
- Swiss Institute of Bioinformatics, 1015 Lausanne, Switzerland
| | - Arlin Stoltzfus
- Office of Data and Informatics, Material Measurement Laboratory, National Institute of Standards and Technology, Rockville, MD 20899, USA
- Institute for Bioscience and Biotechnology Research, Rockville, MD 20850, USA
| | - David M. McCandlish
- Simons Center for Quantitative Biology, Cold Spring Harbor Laboratory, Cold Spring Harbor, NY 11724, USA
| | - Joshua L. Payne
- Institute of Integrative Biology, ETH Zurich, 8092 Zurich, Switzerland
- Swiss Institute of Bioinformatics, 1015 Lausanne, Switzerland
| |
Collapse
|
7
|
Venkataram S, Kryazhimskiy S. Evolutionary repeatability of emergent properties of ecological communities. Philos Trans R Soc Lond B Biol Sci 2023; 378:20220047. [PMID: 37004728 PMCID: PMC10067272 DOI: 10.1098/rstb.2022.0047] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/12/2022] [Accepted: 12/07/2022] [Indexed: 04/04/2023] Open
Abstract
Most species belong to ecological communities where their interactions give rise to emergent community-level properties, such as diversity and productivity. Understanding and predicting how these properties change over time has been a major goal in ecology, with important practical implications for sustainability and human health. Less attention has been paid to the fact that community-level properties can also change because member species evolve. Yet, our ability to predict long-term eco-evolutionary dynamics hinges on how repeatably community-level properties change as a result of species evolution. Here, we review studies of evolution of both natural and experimental communities and make the case that community-level properties at least sometimes evolve repeatably. We discuss challenges faced in investigations of evolutionary repeatability. In particular, only a handful of studies enable us to quantify repeatability. We argue that quantifying repeatability at the community level is critical for approaching what we see as three major open questions in the field: (i) Is the observed degree of repeatability surprising? (ii) How is evolutionary repeatability at the community level related to repeatability at the level of traits of member species? (iii) What factors affect repeatability? We outline some theoretical and empirical approaches to addressing these questions. Advances in these directions will not only enrich our basic understanding of evolution and ecology but will also help us predict eco-evolutionary dynamics. This article is part of the theme issue 'Interdisciplinary approaches to predicting evolutionary biology'.
Collapse
Affiliation(s)
- Sandeep Venkataram
- Department of Ecology, Behavior and Evolution, UC San Diego, La Jolla, CA 92093, USA
| | - Sergey Kryazhimskiy
- Department of Ecology, Behavior and Evolution, UC San Diego, La Jolla, CA 92093, USA
| |
Collapse
|
8
|
Wortel MT, Agashe D, Bailey SF, Bank C, Bisschop K, Blankers T, Cairns J, Colizzi ES, Cusseddu D, Desai MM, van Dijk B, Egas M, Ellers J, Groot AT, Heckel DG, Johnson ML, Kraaijeveld K, Krug J, Laan L, Lässig M, Lind PA, Meijer J, Noble LM, Okasha S, Rainey PB, Rozen DE, Shitut S, Tans SJ, Tenaillon O, Teotónio H, de Visser JAGM, Visser ME, Vroomans RMA, Werner GDA, Wertheim B, Pennings PS. Towards evolutionary predictions: Current promises and challenges. Evol Appl 2023; 16:3-21. [PMID: 36699126 PMCID: PMC9850016 DOI: 10.1111/eva.13513] [Citation(s) in RCA: 15] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2022] [Revised: 11/11/2022] [Accepted: 11/14/2022] [Indexed: 12/14/2022] Open
Abstract
Evolution has traditionally been a historical and descriptive science, and predicting future evolutionary processes has long been considered impossible. However, evolutionary predictions are increasingly being developed and used in medicine, agriculture, biotechnology and conservation biology. Evolutionary predictions may be used for different purposes, such as to prepare for the future, to try and change the course of evolution or to determine how well we understand evolutionary processes. Similarly, the exact aspect of the evolved population that we want to predict may also differ. For example, we could try to predict which genotype will dominate, the fitness of the population or the extinction probability of a population. In addition, there are many uses of evolutionary predictions that may not always be recognized as such. The main goal of this review is to increase awareness of methods and data in different research fields by showing the breadth of situations in which evolutionary predictions are made. We describe how diverse evolutionary predictions share a common structure described by the predictive scope, time scale and precision. Then, by using examples ranging from SARS-CoV2 and influenza to CRISPR-based gene drives and sustainable product formation in biotechnology, we discuss the methods for predicting evolution, the factors that affect predictability and how predictions can be used to prevent evolution in undesirable directions or to promote beneficial evolution (i.e. evolutionary control). We hope that this review will stimulate collaboration between fields by establishing a common language for evolutionary predictions.
Collapse
Affiliation(s)
- Meike T. Wortel
- Swammerdam Institute for Life SciencesUniversity of AmsterdamAmsterdamThe Netherlands
| | - Deepa Agashe
- National Centre for Biological SciencesBangaloreIndia
| | | | - Claudia Bank
- Institute of Ecology and EvolutionUniversity of BernBernSwitzerland
- Swiss Institute of BioinformaticsLausanneSwitzerland
- Gulbenkian Science InstituteOeirasPortugal
| | - Karen Bisschop
- Institute for Biodiversity and Ecosystem DynamicsUniversity of AmsterdamAmsterdamThe Netherlands
- Origins CenterGroningenThe Netherlands
- Laboratory of Aquatic Biology, KU Leuven KulakKortrijkBelgium
| | - Thomas Blankers
- Institute for Biodiversity and Ecosystem DynamicsUniversity of AmsterdamAmsterdamThe Netherlands
- Origins CenterGroningenThe Netherlands
| | | | - Enrico Sandro Colizzi
- Origins CenterGroningenThe Netherlands
- Mathematical InstituteLeiden UniversityLeidenThe Netherlands
| | | | | | - Bram van Dijk
- Max Planck Institute for Evolutionary BiologyPlönGermany
| | - Martijn Egas
- Institute for Biodiversity and Ecosystem DynamicsUniversity of AmsterdamAmsterdamThe Netherlands
| | - Jacintha Ellers
- Department of Ecological ScienceVrije Universiteit AmsterdamAmsterdamThe Netherlands
| | - Astrid T. Groot
- Institute for Biodiversity and Ecosystem DynamicsUniversity of AmsterdamAmsterdamThe Netherlands
| | | | | | - Ken Kraaijeveld
- Leiden Centre for Applied BioscienceUniversity of Applied Sciences LeidenLeidenThe Netherlands
| | - Joachim Krug
- Institute for Biological PhysicsUniversity of CologneCologneGermany
| | - Liedewij Laan
- Department of Bionanoscience, Kavli Institute of NanoscienceTU DelftDelftThe Netherlands
| | - Michael Lässig
- Institute for Biological PhysicsUniversity of CologneCologneGermany
| | - Peter A. Lind
- Department Molecular BiologyUmeå UniversityUmeåSweden
| | - Jeroen Meijer
- Theoretical Biology and Bioinformatics, Department of BiologyUtrecht UniversityUtrechtThe Netherlands
| | - Luke M. Noble
- Institute de Biologie, École Normale Supérieure, CNRS, InsermParisFrance
| | | | - Paul B. Rainey
- Department of Microbial Population BiologyMax Planck Institute for Evolutionary BiologyPlönGermany
- Laboratoire Biophysique et Évolution, CBI, ESPCI Paris, Université PSL, CNRSParisFrance
| | - Daniel E. Rozen
- Institute of Biology, Leiden UniversityLeidenThe Netherlands
| | - Shraddha Shitut
- Origins CenterGroningenThe Netherlands
- Institute of Biology, Leiden UniversityLeidenThe Netherlands
| | | | | | | | | | - Marcel E. Visser
- Department of Animal EcologyNetherlands Institute of Ecology (NIOO‐KNAW)WageningenThe Netherlands
| | - Renske M. A. Vroomans
- Origins CenterGroningenThe Netherlands
- Informatics InstituteUniversity of AmsterdamAmsterdamThe Netherlands
| | | | - Bregje Wertheim
- Groningen Institute for Evolutionary Life SciencesUniversity of GroningenGroningenThe Netherlands
| | | |
Collapse
|
9
|
Coradduzza E, Sanna D, Scarpa F, Azzena I, Fiori MS, Scivoli R, Rocchigiani AM, Bechere R, Dettori MA, Pintus D, Evangelista E, Casu M, Ligios C, Puggioni G. A Deeper Insight into Evolutionary Patterns and Phylogenetic History of ORF Virus through the Whole Genome Sequencing of the First Italian Strains. Viruses 2022; 14:v14071473. [PMID: 35891452 PMCID: PMC9318404 DOI: 10.3390/v14071473] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2022] [Revised: 06/21/2022] [Accepted: 06/28/2022] [Indexed: 02/01/2023] Open
Abstract
Orf virus (ORFV) is distributed worldwide and is the causative agent of contagious ecthyma that mainly occurs in sheep and goats. This disease was reported for the first time at the end of 18th century in Europe but very little is currently known about the temporal and geographic origins of this virus. In the present study, the use of new Italian whole genomes allowed for better inference on the evolutionary history of ORFV. In accordance with previous studies, two genome types (S and G) were described for infection of sheep and goats, respectively. These two well-differentiated groups of genomes originated for evolutive convergence in the late 1800s in two different areas of the world (Europe for S type and Asia for G type), but it was only in the early 1900s that the effective size of ORFV increased among hosts and the virus spread across the whole European continent. The Italian strains which were sequenced in the present study were isolated on the Mediterranean island of Sardinian and showed to be exclusive to this geographic area. One of them is likely representative of the early European forms of ORFV which infected sheep and became extinct about one century ago. Such an ancient Sardinian strain may have reached the island simple by chance, where it quickly adapted to the new habitat.
Collapse
Affiliation(s)
- Elisabetta Coradduzza
- Istituto Zooprofilattico Sperimentale della Sardegna, 07100 Sassari, Italy; (E.C.); (M.S.F.); (R.S.); (A.M.R.); (R.B.); (M.A.D.); (D.P.); (C.L.); (G.P.)
| | - Daria Sanna
- Dipartimento di Scienze Biomediche, Università di Sassari, 07100 Sassari, Italy; (F.S.); (I.A.)
- Correspondence:
| | - Fabio Scarpa
- Dipartimento di Scienze Biomediche, Università di Sassari, 07100 Sassari, Italy; (F.S.); (I.A.)
| | - Ilenia Azzena
- Dipartimento di Scienze Biomediche, Università di Sassari, 07100 Sassari, Italy; (F.S.); (I.A.)
- Dipartimento di Medicina Veterinaria, Università di Sassari, 07100 Sassari, Italy;
| | - Mariangela S. Fiori
- Istituto Zooprofilattico Sperimentale della Sardegna, 07100 Sassari, Italy; (E.C.); (M.S.F.); (R.S.); (A.M.R.); (R.B.); (M.A.D.); (D.P.); (C.L.); (G.P.)
| | - Rosario Scivoli
- Istituto Zooprofilattico Sperimentale della Sardegna, 07100 Sassari, Italy; (E.C.); (M.S.F.); (R.S.); (A.M.R.); (R.B.); (M.A.D.); (D.P.); (C.L.); (G.P.)
| | - Angela M. Rocchigiani
- Istituto Zooprofilattico Sperimentale della Sardegna, 07100 Sassari, Italy; (E.C.); (M.S.F.); (R.S.); (A.M.R.); (R.B.); (M.A.D.); (D.P.); (C.L.); (G.P.)
| | - Roberto Bechere
- Istituto Zooprofilattico Sperimentale della Sardegna, 07100 Sassari, Italy; (E.C.); (M.S.F.); (R.S.); (A.M.R.); (R.B.); (M.A.D.); (D.P.); (C.L.); (G.P.)
| | - Maria A. Dettori
- Istituto Zooprofilattico Sperimentale della Sardegna, 07100 Sassari, Italy; (E.C.); (M.S.F.); (R.S.); (A.M.R.); (R.B.); (M.A.D.); (D.P.); (C.L.); (G.P.)
| | - Davide Pintus
- Istituto Zooprofilattico Sperimentale della Sardegna, 07100 Sassari, Italy; (E.C.); (M.S.F.); (R.S.); (A.M.R.); (R.B.); (M.A.D.); (D.P.); (C.L.); (G.P.)
| | - Eloisa Evangelista
- Ames Polydiagnostic Group Center SRL, Casalnuovo di Napoli, 80013 Naples, Italy;
| | - Marco Casu
- Dipartimento di Medicina Veterinaria, Università di Sassari, 07100 Sassari, Italy;
| | - Ciriaco Ligios
- Istituto Zooprofilattico Sperimentale della Sardegna, 07100 Sassari, Italy; (E.C.); (M.S.F.); (R.S.); (A.M.R.); (R.B.); (M.A.D.); (D.P.); (C.L.); (G.P.)
| | - Giantonella Puggioni
- Istituto Zooprofilattico Sperimentale della Sardegna, 07100 Sassari, Italy; (E.C.); (M.S.F.); (R.S.); (A.M.R.); (R.B.); (M.A.D.); (D.P.); (C.L.); (G.P.)
| |
Collapse
|
10
|
Yashvardhini N, Jha DK, Kumar A, Gaurav M, Sayrav K. Genome sequence analysis of nsp15 from SARS-CoV-2. Bioinformation 2022; 18:432-437. [PMID: 36909703 PMCID: PMC9997503 DOI: 10.6026/97320630018432] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2022] [Revised: 04/30/2022] [Accepted: 04/30/2022] [Indexed: 11/23/2022] Open
Abstract
SARS-CoV-2 (Severe Acute Respiratory Syndrome), a causative agent of COVID-19 disease created a pandemic situation worldwide. Nsp15 is a uridine specific endoribonuclease encoded by the genome of SARS-CoV-2. It plays important role in processing viral RNA and, thus evades the host immune system. Therefore, it is of interest to identify mutants of nsp15 amongst Asian SARS-CoV-2 isolates, where a total of 1795 mutations, from 7793 sequences of Asia submitted till 31st January 2022, amongst which A231V, H234Y, K109N, K259R and S261A mutations were found frequent. Hence, we report data on the predicted secondary structure of wild type form followed by hydropathy plot, physiochemical properties, Ramachandran plot, B-cell epitopes prediction and protein modeling of wild type and mutant of nsp15 protein. Data shows that nsp15 of SARS-CoV-2 is a pontential candidate for the development of vaccine to control the infections of SARS-CoV-2.
Collapse
Affiliation(s)
- Niti Yashvardhini
- Department of Microbiology, Patna Women’s College, Patna, 800 001, Bihar, India
| | - Deepak Kumar Jha
- Department of Zoology, S.M.P. Girls Degree College, Ballia, 277401, Uttar Pradesh, India
| | - Amit Kumar
- Department of Botany, Patna University, Patna-800 005, Bihar, India
| | - Manjush Gaurav
- Department of Botany, Patna University, Patna-800 005, Bihar, India
| | - Kumar Sayrav
- Department of Chemistry, V.K.S. University, Ara-802301, Bihar India
| |
Collapse
|
11
|
Gaunt MW, Pettersson JHO, Kuno G, Gaunt B, de Lamballerie X, Gould EA. Widespread Interspecific Phylogenetic Tree Incongruence Between Mosquito-Borne and Insect-Specific Flaviviruses at Hotspots Originally Identified in Zika Virus. Virus Evol 2022; 8:veac027. [PMID: 35591877 PMCID: PMC9113262 DOI: 10.1093/ve/veac027] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/24/2020] [Revised: 10/22/2021] [Accepted: 04/17/2022] [Indexed: 11/16/2022] Open
Abstract
Intraspecies (homologous) phylogenetic incongruence, or ‘tree conflict’ between different loci within the same genome of mosquito-borne flaviviruses (MBFV), was first identified in dengue virus (DENV) and subsequently in Japanese encephalitis virus (JEV), St Louis encephalitis virus, and Zika virus (ZIKV). Recently, the first evidence of phylogenetic incongruence between interspecific members of the MBFV was reported in ZIKV and its close relative, Spondweni virus. Uniquely, these hybrid proteomes were derived from four incongruent trees involving an Aedes-associated DENV node (1 tree) and three different Culex-associated flavivirus nodes (3 trees). This analysis has now been extended across a wider spectrum of viruses within the MBFV lineage targeting the breakpoints between phylogenetic incongruent loci originally identified in ZIKV. Interspecies phylogenetic incongruence at these breakpoints was identified in 10 of 50 viruses within the MBFV lineage, representing emergent Aedes and Culex-associated viruses including JEV, West Nile virus, yellow fever virus, and insect-specific viruses. Thus, interspecies phylogenetic incongruence is widespread amongst the flaviviruses and is robustly associated with the specific breakpoints that coincide with the interspecific phylogenetic incongruence previously identified, inferring they are ‘hotspots’. The incongruence amongst the emergent MBFV group was restricted to viruses within their respective associated epidemiological boundaries. This MBFV group was RY-coded at the third codon position (‘wobble codon’) to remove transition saturation. The resulting ‘wobble codon’ trees presented a single topology for the entire genome that lacked any robust evidence of phylogenetic incongruence between loci. Phylogenetic interspecific incongruence was therefore observed for exactly the same loci between amino acid and the RY-coded ‘wobble codon’ alignments and this incongruence represented either a major part, or the entire genomes. Maximum likelihood codon analysis revealed positive selection for the incongruent lineages. Positive selection could result in the same locus producing two opposing trees. These analyses for the clinically important MBFV suggest that robust interspecific phylogenetic incongruence resulted from amino acid selection. Convergent or parallel evolutions are evolutionary processes that would explain the observation, whilst interspecific recombination is unlikely.
Collapse
Affiliation(s)
- Michael W Gaunt
- London School of Hygiene and Tropical Medicine, Keppel Street, London WC1E 7HT, United Kingdom
| | - John H-O Pettersson
- Department of Medical Biochemistry and Microbiology, Zoonosis Science Center, Uppsala University, Uppsala, Sweden
- Sydney Institute for Infectious Diseases, School of Life and Environmental Sciences and School of Medical Sciences, the University of Sydney, Sydney, New South Wales 2006, Australia
| | - Goro Kuno
- Formerly, Centers for Disease Control, Fort Collins, CO 80521, USA
| | - Bill Gaunt
- Aeon-sys, MBCS Kensington Road, Barnsley S75 2TU, UK
| | - Xavier de Lamballerie
- UMR “Unité des Virus Emergents”, Aix-Marseille Université-IRD 190-Inserm 1207-IHU Méditerranée Infection, Marseille, France
- APHM Public Hospitals of Marseille, Institut Hospitalo-Universitaire Méditerranée Infection, Marseille, France
| | - Ernest A Gould
- UMR “Unité des Virus Emergents”, Aix-Marseille Université-IRD 190-Inserm 1207-IHU Méditerranée Infection, Marseille, France
| |
Collapse
|
12
|
Population size mediates the contribution of high-rate and large-benefit mutations to parallel evolution. Nat Ecol Evol 2022; 6:439-447. [PMID: 35241808 DOI: 10.1038/s41559-022-01669-3] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2021] [Accepted: 01/11/2022] [Indexed: 12/15/2022]
Abstract
Mutations with large fitness benefits and mutations occurring at high rates may both cause parallel evolution, but their contribution is predicted to depend on population size. Moreover, high-rate and large-benefit mutations may have different long-term adaptive consequences. We show that small and 100-fold larger bacterial populations evolve resistance to a β-lactam antibiotic by using similar numbers, but different types of mutations. Small populations frequently substitute similar high-rate structural variants and loss-of-function point mutations, including the deletion of a low-activity β-lactamase, and evolve modest resistance levels. Large populations more often use low-rate, large-benefit point mutations affecting the same targets, including mutations activating the β-lactamase and other gain-of-function mutations, leading to much higher resistance levels. Our results demonstrate the separation by clonal interference of mutation classes with divergent adaptive consequences, causing a shift from high-rate to large-benefit mutations with increases in population size.
Collapse
|
13
|
Abstract
How do mutational biases influence the process of adaptation? A common assumption is that selection alone determines the course of adaptation from abundant preexisting variation. Yet, theoretical work shows broad conditions under which the mutation rate to a given type of variant strongly influences its probability of contributing to adaptation. Here we introduce a statistical approach to analyzing how mutation shapes protein sequence adaptation. Using large datasets from three different species, we show that the mutation spectrum has a proportional influence on the types of changes fixed in adaptation. We also show via computer simulations that a variety of factors can influence how closely the spectrum of adaptive substitutions reflects the spectrum of variants introduced by mutation. Evolutionary adaptation often occurs by the fixation of beneficial mutations. This mode of adaptation can be characterized quantitatively by a spectrum of adaptive substitutions, i.e., a distribution for types of changes fixed in adaptation. Recent work establishes that the changes involved in adaptation reflect common types of mutations, raising the question of how strongly the mutation spectrum shapes the spectrum of adaptive substitutions. We address this question with a codon-based model for the spectrum of adaptive amino acid substitutions, applied to three large datasets covering thousands of amino acid changes identified in natural and experimental adaptation in Saccharomyces cerevisiae, Escherichia coli, and Mycobacterium tuberculosis. Using species-specific mutation spectra based on prior knowledge, we find that the mutation spectrum has a proportional influence on the spectrum of adaptive substitutions in all three species. Indeed, we find that by inferring the mutation rates that best explain the spectrum of adaptive substitutions, we can accurately recover the species-specific mutation spectra. However, we also find that the predictive power of the model differs substantially between the three species. To better understand these differences, we use population simulations to explore the factors that influence how closely the spectrum of adaptive substitutions mirrors the mutation spectrum. The results show that the influence of the mutation spectrum decreases with increasing mutational supply (Nμ) and that predictive power is strongly affected by the number and diversity of beneficial mutations.
Collapse
|
14
|
Horton JS, Flanagan LM, Jackson RW, Priest NK, Taylor TB. A mutational hotspot that determines highly repeatable evolution can be built and broken by silent genetic changes. Nat Commun 2021; 12:6092. [PMID: 34667151 PMCID: PMC8526746 DOI: 10.1038/s41467-021-26286-9] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2021] [Accepted: 09/28/2021] [Indexed: 11/08/2022] Open
Abstract
Mutational hotspots can determine evolutionary outcomes and make evolution repeatable. Hotspots are products of multiple evolutionary forces including mutation rate heterogeneity, but this variable is often hard to identify. In this work, we reveal that a near-deterministic genetic hotspot can be built and broken by a handful of silent mutations. We observe this when studying homologous immotile variants of the bacteria Pseudomonas fluorescens, AR2 and Pf0-2x. AR2 resurrects motility through highly repeatable de novo mutation of the same nucleotide in >95% lines in minimal media (ntrB A289C). Pf0-2x, however, evolves via a number of mutations meaning the two strains diverge significantly during adaptation. We determine that this evolutionary disparity is owed to just 6 synonymous variations within the ntrB locus, which we demonstrate by swapping the sites and observing that we are able to both break (>95% to 0%) and build (0% to 80%) a deterministic mutational hotspot. Our work reveals a key role for silent genetic variation in determining adaptive outcomes.
Collapse
Affiliation(s)
- James S Horton
- Milner Centre for Evolution, Department of Biology & Biochemistry, University of Bath, Claverton Down, Bath, BA2 7AY, UK.
| | - Louise M Flanagan
- Milner Centre for Evolution, Department of Biology & Biochemistry, University of Bath, Claverton Down, Bath, BA2 7AY, UK
| | - Robert W Jackson
- School of Biosciences and Birmingham Institute of Forest Research (BIFoR), University of Birmingham, Edgbaston, Birmingham, B15 2TT, UK
| | - Nicholas K Priest
- Milner Centre for Evolution, Department of Biology & Biochemistry, University of Bath, Claverton Down, Bath, BA2 7AY, UK
| | - Tiffany B Taylor
- Milner Centre for Evolution, Department of Biology & Biochemistry, University of Bath, Claverton Down, Bath, BA2 7AY, UK.
| |
Collapse
|
15
|
Sarkar R, Mitra S, Chandra P, Saha P, Banerjee A, Dutta S, Chawla-Sarkar M. Comprehensive analysis of genomic diversity of SARS-CoV-2 in different geographic regions of India: an endeavour to classify Indian SARS-CoV-2 strains on the basis of co-existing mutations. Arch Virol 2021; 166:801-812. [PMID: 33464421 PMCID: PMC7814186 DOI: 10.1007/s00705-020-04911-0] [Citation(s) in RCA: 36] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/22/2020] [Accepted: 10/21/2020] [Indexed: 01/24/2023]
Abstract
Accumulation of mutations within the genome is the primary driving force in viral evolution within an endemic setting. This inherent feature often leads to altered virulence, infectivity and transmissibility, and antigenic shifts to escape host immunity, which might compromise the efficacy of vaccines and antiviral drugs. Therefore, we carried out a genome-wide analysis of circulating SARS-CoV-2 strains to detect the emergence of novel co-existing mutations and trace their geographical distribution within India. Comprehensive analysis of whole genome sequences of 837 Indian SARS-CoV-2 strains revealed the occurrence of 33 different mutations, 18 of which were unique to India. Novel mutations were observed in the S glycoprotein (6/33), NSP3 (5/33), RdRp/NSP12 (4/33), NSP2 (2/33), and N (1/33). Non-synonymous mutations were found to be 3.07 times more prevalent than synonymous mutations. We classified the Indian isolates into 22 groups based on their co-existing mutations. Phylogenetic analysis revealed that the representative strains of each group were divided into various sub-clades within their respective clades, based on the presence of unique co-existing mutations. The A2a clade was found to be dominant in India (71.34%), followed by A3 (23.29%) and B (5.36%), but a heterogeneous distribution was observed among various geographical regions. The A2a clade was highly predominant in East India, Western India, and Central India, whereas the A2a and A3 clades were nearly equal in prevalence in South and North India. This study highlights the divergent evolution of SARS-CoV-2 strains and co-circulation of multiple clades in India. Monitoring of the emerging mutations will pave the way for vaccine formulation and the design of antiviral drugs.
Collapse
Affiliation(s)
- Rakesh Sarkar
- Division of Virology, National Institute of Cholera and Enteric Diseases, P-33, C.I.T. Road, Scheme-XM, Beliaghata, Kolkata, West Bengal, 700010, India
| | - Suvrotoa Mitra
- Division of Virology, National Institute of Cholera and Enteric Diseases, P-33, C.I.T. Road, Scheme-XM, Beliaghata, Kolkata, West Bengal, 700010, India
| | - Pritam Chandra
- Division of Virology, National Institute of Cholera and Enteric Diseases, P-33, C.I.T. Road, Scheme-XM, Beliaghata, Kolkata, West Bengal, 700010, India
| | - Priyanka Saha
- Division of Virology, National Institute of Cholera and Enteric Diseases, P-33, C.I.T. Road, Scheme-XM, Beliaghata, Kolkata, West Bengal, 700010, India
| | - Anindita Banerjee
- Division of Virology, National Institute of Cholera and Enteric Diseases, P-33, C.I.T. Road, Scheme-XM, Beliaghata, Kolkata, West Bengal, 700010, India
| | - Shanta Dutta
- Division of Virology, National Institute of Cholera and Enteric Diseases, P-33, C.I.T. Road, Scheme-XM, Beliaghata, Kolkata, West Bengal, 700010, India
| | - Mamta Chawla-Sarkar
- Division of Virology, National Institute of Cholera and Enteric Diseases, P-33, C.I.T. Road, Scheme-XM, Beliaghata, Kolkata, West Bengal, 700010, India.
| |
Collapse
|
16
|
Yashvardhini N, Jha DK, Bhattacharya S. Identification and characterization of mutations in the SARS-CoV-2 RNA-dependent RNA polymerase as a promising antiviral therapeutic target. Arch Microbiol 2021; 203:5463-5473. [PMID: 34410443 PMCID: PMC8374121 DOI: 10.1007/s00203-021-02527-9] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/24/2020] [Revised: 08/11/2021] [Accepted: 08/12/2021] [Indexed: 12/29/2022]
Abstract
The causative agent of COVID-19 is a novel betacoronavirus or severe acute respiratory syndrome coronavirus (SARS-CoV-2), which has emerged as a pandemic of global concern. Considering its rapid transmission, WHO has declared public health emergency on 11th March 2020 worldwide. SARS-CoV-2 is a genetically diverse positive sense RNA virus that typically exhibit high rates of mutation than DNA viruses. Higher rates of mutation bring higher genomic variability which may lead to viral evolution and enabling viruses to evade the pre-existing immunity of host and quickly acquire drug resistance properties. The objective of our study was to compare the SARS-CoV-2 RdRp sequences of Indian SARS-CoV-2 isolates with those of Wuhan type virus. A total of 384 point mutations were detected from 488 sequence of the RdRp protein of Indian SARS-CoV-2 genome, out of which seven were used for subsequent study. Furthermore, prediction of secondary structure, protein modeling and its dynamics were performed which revealed that seven mutations (R118C, T148I, Y149C, E802A, Q822H, V880I and D893Y) significantly altered the stability and flexibility of RdRp protein. Present study was therefore, undertaken to analyze the variations occurring in RdRp due to multiple mutations leading to the alterations in the structure and function of RNA-dependent RNA polymerase which is essential for the replication /transcription of this virus and hence can be utilized as a promising therapeutic target to curb SARS-CoV-2 infections.
Collapse
Affiliation(s)
- Niti Yashvardhini
- grid.412457.10000 0001 1276 6626Department of Microbiology, Patna Women’s College, Patna, 800 001 India
| | - Deepak Kumar Jha
- Department of Zoology, P. C. Vigyan Mahavidyalaya, J. P. University, Chapra, 841 301 India
| | - Saurav Bhattacharya
- grid.449713.c0000 0004 5944 7827Department of Biotechnology, Techno India University, Kolkata, 700 091 India
| |
Collapse
|
17
|
Gomez K, Bertram J, Masel J. Mutation bias can shape adaptation in large asexual populations experiencing clonal interference. Proc Biol Sci 2020; 287:20201503. [PMID: 33081612 DOI: 10.1098/rspb.2020.1503] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
Abstract
The extended evolutionary synthesis invokes a role for development in shaping adaptive evolution, which in population genetics terms corresponds to mutation-biased adaptation. Critics have claimed that clonal interference makes mutation-biased adaptation rare. We consider the behaviour of two simultaneously adapting traits, one with larger mutation rate U, the other with larger selection coefficient s, using asexual travelling wave models. We find that adaptation is dominated by whichever trait has the faster rate of adaptation v in isolation, with the other trait subject to evolutionary stalling. Reviewing empirical claims for mutation-biased adaptation, we find that not all occur in the 'origin-fixation' regime of population genetics where v is only twice as sensitive to s as to U. In some cases, differences in U are at least ten to twelve times larger than differences in s, as needed to cause mutation-biased adaptation even in the 'multiple mutations' regime. Surprisingly, when U > s in the 'diffusive-mutation' regime, the required sensitivity ratio is also only two, despite pervasive clonal interference. Given two traits with identical v, the benefit of having higher s is surprisingly small, occurring largely when one trait is at the boundary between the origin-fixation and multiple mutations regimes.
Collapse
Affiliation(s)
- Kevin Gomez
- Graduate Interdisciplinary Program in Applied Mathematics, University of Arizona, Tucson, AZ, USA
| | - Jason Bertram
- Environmental Resilience Institute, Indiana University, Bloomington, IN, USA.,Department of Biology, Indiana University, Bloomington, IN, USA
| | - Joanna Masel
- Department of Ecology and Evolutionary Biology, University of Arizona, Tucson, AZ, USA
| |
Collapse
|
18
|
Chand GB, Banerjee A, Azad GK. Identification of novel mutations in RNA-dependent RNA polymerases of SARS-CoV-2 and their implications on its protein structure. PeerJ 2020; 8:e9492. [PMID: 32685291 PMCID: PMC7337032 DOI: 10.7717/peerj.9492] [Citation(s) in RCA: 44] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/10/2020] [Accepted: 06/16/2020] [Indexed: 01/18/2023] Open
Abstract
The rapid development of the SARS-CoV-2 mediated COVID-19 pandemic has been the cause of significant health concern, highlighting the immediate need for effective antivirals. SARS-CoV-2 is an RNA virus that has an inherently high mutation rate. These mutations drive viral evolution and genome variability, thereby facilitating viruses to have rapid antigenic shifting to evade host immunity and to develop drug resistance. Viral RNA-dependent RNA polymerases (RdRp) perform viral genome duplication and RNA synthesis. Therefore, we compared the available RdRp sequences of SARS-CoV-2 from Indian isolates and the ‘Wuhan wet sea food market virus’ sequence to identify, if any, variation between them. Our data revealed the occurrence of seven mutations in Indian isolates of SARS-CoV-2. The secondary structure prediction analysis of these seven mutations shows that three of them cause alteration in the structure of RdRp. Furthermore, we did protein modelling studies to show that these mutations can potentially alter the stability of the RdRp protein. Therefore, we propose that RdRp mutations in Indian SARS-CoV-2 isolates might have functional consequences that can interfere with RdRp targeting pharmacological agents.
Collapse
Affiliation(s)
| | - Atanu Banerjee
- Department of Zoology, Samastipur College, Samastipur, Bihar, India
| | | |
Collapse
|
19
|
Bertels F, Leemann C, Metzner KJ, Regoes R. Parallel evolution of HIV-1 in a long-term experiment. Mol Biol Evol 2019; 36:2400-2414. [PMID: 31251344 PMCID: PMC6805227 DOI: 10.1093/molbev/msz155] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/12/2018] [Revised: 05/06/2019] [Accepted: 06/22/2019] [Indexed: 12/15/2022] Open
Abstract
One of the most intriguing puzzles in biology is the degree to which evolution is repeatable. The repeatability of evolution, or parallel evolution, has been studied in a variety of model systems, but has rarely been investigated with clinically relevant viruses. To investigate parallel evolution of HIV-1, we passaged two replicate HIV-1 populations for almost 1 year in each of two human T-cell lines. For each of the four evolution lines, we determined the genetic composition of the viral population at nine time points by deep sequencing the entire genome. Mutations that were carried by the majority of the viral population accumulated continuously over 1 year in each evolution line. Many majority mutations appeared in more than one evolution line, that is, our experiments showed an extreme degree of parallel evolution. In one of the evolution lines, 62% of the majority mutations also occur in another line. The parallelism impairs our ability to reconstruct the evolutionary history by phylogenetic methods. We show that one can infer the correct phylogenetic topology by including minority mutations in our analysis. We also find that mutation diversity at the beginning of the experiment is predictive of the frequency of majority mutations at the end of the experiment.
Collapse
Affiliation(s)
- Frederic Bertels
- Department of Environmental Systems Sciences, ETH Zurich, Zurich.,Max-Planck-Institute for Evolutionary Biology, Department of Microbial Population Biology
| | - Christine Leemann
- Division of Infectious Diseases and Hospital Epidemiology, University Hospital Zurich, Zurich.,Insitute of Medical Virology, University of Zurich, Zurich
| | - Karin J Metzner
- Division of Infectious Diseases and Hospital Epidemiology, University Hospital Zurich, Zurich.,Insitute of Medical Virology, University of Zurich, Zurich
| | - Roland Regoes
- Department of Environmental Systems Sciences, ETH Zurich, Zurich
| |
Collapse
|
20
|
Payne JL, Menardo F, Trauner A, Borrell S, Gygli SM, Loiseau C, Gagneux S, Hall AR. Transition bias influences the evolution of antibiotic resistance in Mycobacterium tuberculosis. PLoS Biol 2019; 17:e3000265. [PMID: 31083647 PMCID: PMC6532934 DOI: 10.1371/journal.pbio.3000265] [Citation(s) in RCA: 33] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/26/2018] [Revised: 05/23/2019] [Accepted: 04/26/2019] [Indexed: 11/21/2022] Open
Abstract
Transition bias, an overabundance of transitions relative to transversions, has been widely reported among studies of the rates and spectra of spontaneous mutations. However, demonstrating the role of transition bias in adaptive evolution remains challenging. In particular, it is unclear whether such biases direct the evolution of bacterial pathogens adapting to treatment. We addressed this challenge by analyzing adaptive antibiotic-resistance mutations in the major human pathogen Mycobacterium tuberculosis (MTB). We found strong evidence for transition bias in two independently curated data sets comprising 152 and 208 antibiotic-resistance mutations. This was true at the level of mutational paths (distinct adaptive DNA sequence changes) and events (individual instances of the adaptive DNA sequence changes) and across different genes and gene promoters conferring resistance to a diversity of antibiotics. It was also true for mutations that do not code for amino acid changes (in gene promoters and the 16S ribosomal RNA gene rrs) and for mutations that are synonymous to each other and are therefore likely to have similar fitness effects, suggesting that transition bias can be caused by a bias in mutation supply. These results point to a central role for transition bias in determining which mutations drive adaptive antibiotic resistance evolution in a key pathogen. Some types of mutations occur more frequently than expected. This study shows that such bias —an excess of transitions over transversions—influences the evolution of antibiotic resistance in a key global pathogen, Mycobacterium tuberculosis.
Collapse
Affiliation(s)
- Joshua L. Payne
- Institute of Integrative Biology, ETH Zurich, Switzerland
- Swiss Institute of Bioinformatics, Lausanne, Switzerland
- * E-mail:
| | - Fabrizio Menardo
- Swiss Tropical and Public Health Institute, Basel, Switzerland
- University of Basel, Basel, Switzerland
| | - Andrej Trauner
- Swiss Tropical and Public Health Institute, Basel, Switzerland
- University of Basel, Basel, Switzerland
| | - Sonia Borrell
- Swiss Tropical and Public Health Institute, Basel, Switzerland
- University of Basel, Basel, Switzerland
| | - Sebastian M. Gygli
- Swiss Tropical and Public Health Institute, Basel, Switzerland
- University of Basel, Basel, Switzerland
| | - Chloe Loiseau
- Swiss Tropical and Public Health Institute, Basel, Switzerland
- University of Basel, Basel, Switzerland
| | - Sebastien Gagneux
- Swiss Tropical and Public Health Institute, Basel, Switzerland
- University of Basel, Basel, Switzerland
| | - Alex R. Hall
- Institute of Integrative Biology, ETH Zurich, Switzerland
| |
Collapse
|
21
|
Guterres A, de Oliveira RC, Fernandes J, de Lemos ERS. The mystery of the phylogeographic structural pattern in rodent-borne hantaviruses. Mol Phylogenet Evol 2019; 136:35-43. [PMID: 30914396 DOI: 10.1016/j.ympev.2019.03.020] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/25/2019] [Revised: 03/22/2019] [Accepted: 03/23/2019] [Indexed: 12/19/2022]
Abstract
Hantaviruses (order Bunyavirales, family Hantaviridae) are important zoonotic pathogens. Because of the great diversity of their reservoir hosts, hantaviruses are excellent models to evaluate the dynamics of virus-host co-evolution. To understand the mechanisms behind the evolutionary history of hantaviruses through virus-reservoir interactions, it is important to know how the radiation and diversity of hantaviruses occurred. In this paper, we evaluate the pattern of hantavirus diversification based on a complete S segment representing major groups of hantaviruses found in the Americas. Phylogenetic analyses revealed a high degree of phylogeographic structure and a surprising pattern of geographical distribution of New World hantaviruses. The available data suggest that hantaviruses related to the Arvicolinae rodent subfamily in North America probably emerged and initially adapted from a shared common ancestor of the Tula virus. The first clade of hantaviruses associated with Neotominae occupied a stem lineage, especially those that emerged in Central America or Mexico. Hantaviruses from Central America and Mexico found in Neotominae rodents spread northward and probably gave rise to the first phylogroup of hantaviruses associated with Sigmodontinae in North America. Two preferential host-switching transmissions in hantaviruses apparently gave rise to two different paraphyletic group in Neotominae and Sigmodontinae. Our study supports a probable epicenter of diversification in Central America and/or Mexico for hantaviruses related to both the Neotominae and Sigmodontinae subfamilies.
Collapse
Affiliation(s)
- Alexandro Guterres
- Laboratório de Hantaviroses e Rickettsioses, Instituto Oswaldo Cruz, Fundação Oswaldo Cruz, Rio de Janeiro, RJ, Brazil.
| | - Renata Carvalho de Oliveira
- Laboratório de Hantaviroses e Rickettsioses, Instituto Oswaldo Cruz, Fundação Oswaldo Cruz, Rio de Janeiro, RJ, Brazil
| | - Jorlan Fernandes
- Laboratório de Hantaviroses e Rickettsioses, Instituto Oswaldo Cruz, Fundação Oswaldo Cruz, Rio de Janeiro, RJ, Brazil
| | - Elba Regina Sampaio de Lemos
- Laboratório de Hantaviroses e Rickettsioses, Instituto Oswaldo Cruz, Fundação Oswaldo Cruz, Rio de Janeiro, RJ, Brazil
| |
Collapse
|
22
|
Rowe KL, Minaya MA, Belshe RB, Morrison LA. Molecular analyses and phylogeny of the herpes simplex virus 2 US9 and glycoproteins gE/gI obtained from infected subjects during the Herpevac Trial for Women. PLoS One 2019; 14:e0212877. [PMID: 30849089 PMCID: PMC6407778 DOI: 10.1371/journal.pone.0212877] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2018] [Accepted: 02/11/2019] [Indexed: 11/18/2022] Open
Abstract
Herpes simplex virus 2 (HSV-2) is a large double-stranded DNA virus that causes genital sores when spread by sexual contact and is a principal cause of viral encephalitis in newborns and infants. Viral glycoproteins enable virion entry into and spread between cells, making glycoproteins a prime target for vaccine development. A truncated glycoprotein D2 (gD2) vaccine candidate, recently tested in the phase 3 Herpevac Trial for Women, did not prevent HSV-2 infection in initially seronegative women. Some women who became infected experienced multiple recurrences during the trial. The HSV US7, US8, and US9 genes encode glycoprotein I (gI), glycoprotein E (gE), and the US9 type II membrane protein, respectively. These proteins participate in viral spread across cell junctions and facilitate anterograde transport of virion components in neurons, prompting us to investigate whether sequence variants in these genes could be associated with frequent recurrence. The nucleotide sequences and dN/dS ratios of the US7-US9 region from viral isolates of individuals who experienced multiple recurrences were compared with those who had had a single episode of disease. No consistent polymorphism(s) distinguished the recurrent isolates. In frequently recurring isolates, the dN/dS ratio of US7 was low while greater variation (higher dN/dS ratio) occurred in US8, suggesting conserved function of the former during reactivation. Phylogenetic reconstruction of the US7-US9 region revealed eight strongly supported clusters within the 55 U.S. HSV-2 strains sampled, which were preserved in a second global phylogeny. Thus, although we have demonstrated evolutionary diversity in the US7-US9 complex, we found no molecular evidence of sequence variation in US7-US9 that distinguishes isolates from subjects with frequently recurrent episodes of disease.
Collapse
Affiliation(s)
- Kelsey L. Rowe
- Department of Molecular Microbiology and Immunology, Saint Louis University School of Medicine, St. Louis, Missouri, United States of America
| | - Miguel A. Minaya
- Department of Molecular Microbiology and Immunology, Saint Louis University School of Medicine, St. Louis, Missouri, United States of America
| | - Robert B. Belshe
- Department of Internal Medicine, Division of Infectious Diseases, Saint Louis University School of Medicine, St. Louis, Missouri, United States of America
| | - Lynda A. Morrison
- Department of Molecular Microbiology and Immunology, Saint Louis University School of Medicine, St. Louis, Missouri, United States of America
- Department of Internal Medicine, Division of Infectious Diseases, Saint Louis University School of Medicine, St. Louis, Missouri, United States of America
| |
Collapse
|
23
|
Abstract
How virulence evolves after a virus jumps to a new host species is central to disease emergence. Our current understanding of virulence evolution is based on insights drawn from two perspectives that have developed largely independently: long-standing evolutionary theory based on limited real data examples that often lack a genomic basis, and experimental studies of virulence-determining mutations using cell culture or animal models. A more comprehensive understanding of virulence mutations and their evolution can be achieved by bridging the gap between these two research pathways through the phylogenomic analysis of virus genome sequence data as a guide to experimental study.
Collapse
Affiliation(s)
- Jemma L Geoghegan
- Department of Biological Sciences, Macquarie University, Sydney, New South Wales, Australia
| | - Edward C Holmes
- Marie Bashir Institute for Infectious Diseases and Biosecurity, Charles Perkins Centre, School of Life and Environmental Sciences and Sydney Medical School, The University of Sydney, Sydney, New South Wales, Australia.
| |
Collapse
|
24
|
Van den Bergh B, Swings T, Fauvart M, Michiels J. Experimental Design, Population Dynamics, and Diversity in Microbial Experimental Evolution. Microbiol Mol Biol Rev 2018; 82:e00008-18. [PMID: 30045954 PMCID: PMC6094045 DOI: 10.1128/mmbr.00008-18] [Citation(s) in RCA: 90] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
In experimental evolution, laboratory-controlled conditions select for the adaptation of species, which can be monitored in real time. Despite the current popularity of such experiments, nature's most pervasive biological force was long believed to be observable only on time scales that transcend a researcher's life-span, and studying evolution by natural selection was therefore carried out solely by comparative means. Eventually, microorganisms' propensity for fast evolutionary changes proved us wrong, displaying strong evolutionary adaptations over a limited time, nowadays massively exploited in laboratory evolution experiments. Here, we formulate a guide to experimental evolution with microorganisms, explaining experimental design and discussing evolutionary dynamics and outcomes and how it is used to assess ecoevolutionary theories, improve industrially important traits, and untangle complex phenotypes. Specifically, we give a comprehensive overview of the setups used in experimental evolution. Additionally, we address population dynamics and genetic or phenotypic diversity during evolution experiments and expand upon contributing factors, such as epistasis and the consequences of (a)sexual reproduction. Dynamics and outcomes of evolution are most profoundly affected by the spatiotemporal nature of the selective environment, where changing environments might lead to generalists and structured environments could foster diversity, aided by, for example, clonal interference and negative frequency-dependent selection. We conclude with future perspectives, with an emphasis on possibilities offered by fast-paced technological progress. This work is meant to serve as an introduction to those new to the field of experimental evolution, as a guide to the budding experimentalist, and as a reference work to the seasoned expert.
Collapse
Affiliation(s)
- Bram Van den Bergh
- Laboratory of Symbiotic and Pathogenic Interactions, Centre of Microbial and Plant Genetics, KU Leuven-University of Leuven, Leuven, Belgium
- Michiels Lab, Center for Microbiology, VIB, Leuven, Belgium
- Douglas Lab, Department of Entomology, Cornell University, Ithaca, New York, USA
| | - Toon Swings
- Laboratory of Symbiotic and Pathogenic Interactions, Centre of Microbial and Plant Genetics, KU Leuven-University of Leuven, Leuven, Belgium
- Michiels Lab, Center for Microbiology, VIB, Leuven, Belgium
| | - Maarten Fauvart
- Laboratory of Symbiotic and Pathogenic Interactions, Centre of Microbial and Plant Genetics, KU Leuven-University of Leuven, Leuven, Belgium
- Michiels Lab, Center for Microbiology, VIB, Leuven, Belgium
- imec, Leuven, Belgium
| | - Jan Michiels
- Laboratory of Symbiotic and Pathogenic Interactions, Centre of Microbial and Plant Genetics, KU Leuven-University of Leuven, Leuven, Belgium
- Michiels Lab, Center for Microbiology, VIB, Leuven, Belgium
| |
Collapse
|
25
|
Willis S, Masel J. Gene Birth Contributes to Structural Disorder Encoded by Overlapping Genes. Genetics 2018; 210:303-313. [PMID: 30026186 PMCID: PMC6116962 DOI: 10.1534/genetics.118.301249] [Citation(s) in RCA: 25] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2018] [Accepted: 07/18/2018] [Indexed: 11/18/2022] Open
Abstract
The same nucleotide sequence can encode two protein products in different reading frames. Overlapping gene regions encode higher levels of intrinsic structural disorder (ISD) than nonoverlapping genes (39% vs. 25% in our viral dataset). This might be because of the intrinsic properties of the genetic code, because one member per pair was recently born de novo in a process that favors high ISD, or because high ISD relieves increased evolutionary constraint imposed by dual-coding. Here, we quantify the relative contributions of these three alternative hypotheses. We estimate that the recency of de novo gene birth explains [Formula: see text] or more of the elevation in ISD in overlapping regions of viral genes. While the two reading frames within a same-strand overlapping gene pair have markedly different ISD tendencies that must be controlled for, their effects cancel out to make no net contribution to ISD. The remaining elevation of ISD in the older members of overlapping gene pairs, presumed due to the need to alleviate evolutionary constraint, was already present prior to the origin of the overlap. Same-strand overlapping gene birth events can occur in two different frames, favoring high ISD either in the ancestral gene or in the novel gene; surprisingly, most de novo gene birth events contained completely within the body of an ancestral gene favor high ISD in the ancestral gene (23 phylogenetically independent events vs. 1). This can be explained by mutation bias favoring the frame with more start codons and fewer stop codons.
Collapse
Affiliation(s)
- Sara Willis
- Department of Ecology and Evolutionary Biology, University of Arizona, Tucson, Arizona 85721
| | - Joanna Masel
- Department of Ecology and Evolutionary Biology, University of Arizona, Tucson, Arizona 85721
| |
Collapse
|
26
|
Parallel Evolution of Host-Attachment Proteins in Phage PP01 Populations Adapting to Escherichia coli O157:H7. Pharmaceuticals (Basel) 2018; 11:ph11020060. [PMID: 29925767 PMCID: PMC6027323 DOI: 10.3390/ph11020060] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/12/2018] [Revised: 06/07/2018] [Accepted: 06/11/2018] [Indexed: 12/14/2022] Open
Abstract
The emergence of antibiotic resistance has sparked interest in phage therapy, which uses virulent phages as antibacterial agents. Bacteriophage PP01 has been studied for potential bio-control of Escherichia coli O157:H7, its natural host, but in the laboratory, PP01 can be inefficient at killing this bacterium. Thus, the goal of this study was to improve the therapeutic potential of PP01 through short-term experimental evolution. Four replicate populations of PP01 were serially passaged 21 times on non-evolving E. coli O157:H7 with the prediction that the evolved phage populations would adsorb faster and more efficiently kill the host bacteria. Dead-cell adsorption assays and in vitro killing assays confirmed that evolved viruses improved their adsorption ability on E. coli O157:H7, and adapted to kill host bacteria faster than the wildtype ancestor. Sequencing of candidate tail-fiber genes revealed that the phage populations evolved in parallel; the lineages shared two point mutations in gp38 that encodes a host recognition protein, and surprisingly shared a ~600 bp deletion in gp37 that encodes the distal tail fibers. In contrast, no mutations were observed in the gp12 gene encoding PP01’s short tail fibers. We discuss the functional role of the observed mutations, including the possible adaptive role of the evolved deletions. This study demonstrates how experimental evolution can be used to select for viral traits that improve phage attack of an important bacterial pathogen, and that the molecular targets of selection include loci contributing to cell attachment and phage virulence.
Collapse
|
27
|
Abstract
The deterministic force of natural selection and stochastic influence of drift shape RNA virus evolution. New deep-sequencing and microfluidics technologies allow us to quantify the effect of mutations and trace the evolution of viral populations with single-genome and single-nucleotide resolution. Such experiments can reveal the topography of the genotype-fitness landscapes that shape the path of viral evolution. By combining historical analyses, like phylogenetic approaches, with high-throughput and high-resolution evolutionary experiments, we can observe parallel patterns of evolution that drive important phenotypic transitions. These developments provide a framework for quantifying and anticipating potential evolutionary events. Here, we examine emerging technologies that can map the selective landscapes of viruses, focusing on their application to pathogenic viruses. We identify areas where these technologies can bolster our ability to study the evolution of viruses and to anticipate and possibly intervene in evolutionary events and prevent viral disease.
Collapse
Affiliation(s)
- Patrick T Dolan
- Department of Biology, Stanford University, E200 Clark Center, 318 Campus Drive, Stanford, CA 94305, USA; Department of Microbiology and Immunology, University of California, San Francisco, 600 16th Street, GH-S572, UCSF Box 2280, San Francisco, CA 94143-2280, USA
| | - Zachary J Whitfield
- Department of Microbiology and Immunology, University of California, San Francisco, 600 16th Street, GH-S572, UCSF Box 2280, San Francisco, CA 94143-2280, USA
| | - Raul Andino
- Department of Microbiology and Immunology, University of California, San Francisco, 600 16th Street, GH-S572, UCSF Box 2280, San Francisco, CA 94143-2280, USA.
| |
Collapse
|
28
|
Divergent and parallel routes of biochemical adaptation in high-altitude passerine birds from the Qinghai-Tibet Plateau. Proc Natl Acad Sci U S A 2018; 115:1865-1870. [PMID: 29432191 DOI: 10.1073/pnas.1720487115] [Citation(s) in RCA: 57] [Impact Index Per Article: 9.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023] Open
Abstract
When different species experience similar selection pressures, the probability of evolving similar adaptive solutions may be influenced by legacies of evolutionary history, such as lineage-specific changes in genetic background. Here we test for adaptive convergence in hemoglobin (Hb) function among high-altitude passerine birds that are native to the Qinghai-Tibet Plateau, and we examine whether convergent increases in Hb-O2 affinity have a similar molecular basis in different species. We documented that high-altitude parid and aegithalid species from the Qinghai-Tibet Plateau have evolved derived increases in Hb-O2 affinity in comparison with their closest lowland relatives in East Asia. However, convergent increases in Hb-O2 affinity and convergence in underlying functional mechanisms were seldom attributable to the same amino acid substitutions in different species. Using ancestral protein resurrection and site-directed mutagenesis, we experimentally confirmed two cases in which parallel substitutions contributed to convergent increases in Hb-O2 affinity in codistributed high-altitude species. In one case involving the ground tit (Parus humilis) and gray-crested tit (Lophophanes dichrous), parallel amino acid replacements with affinity-enhancing effects were attributable to nonsynonymous substitutions at a CpG dinucleotide, suggesting a possible role for mutation bias in promoting recurrent changes at the same site. Overall, most altitude-related changes in Hb function were caused by divergent amino acid substitutions, and a select few were caused by parallel substitutions that produced similar phenotypic effects on the divergent genetic backgrounds of different species.
Collapse
|